Solar coronal jets are small, transient, collimated ejections most easily observed in coronal holes (CHs). The upcoming Parker Solar Probe (PSP) mission provides the first opportunity to encounter CH jets in situ near the Sun and examine their internal structure and dynamics. Using projected mission orbital parameters, we have simulated PSP encounters with a fully three-dimensional magnetohydrodynamic (MHD) model of a CH jet. We find that three internal jet regions, featuring different wave modes and levels of compressibility, have distinct identifying signatures detectable by PSP. The leading Alfv\'{e}n wave front and its immediate wake are characterized by trans-Alfv\'{e}nic plasma flows with mild density enhancements. This front exhibits characteristics of a fast switch-on MHD shock, whose arrival is signaled by the sudden onset of large-amplitude transverse velocity and magnetic-field oscillations highly correlated in space and time. The trailing portion is characterized by supersonic but sub-Alfv\'{e}nic outflows of dense plasma with uncorrelated velocity and magnetic-field oscillations. This compressible region contains most of the jet's mass. The volume between the immediate wake and dense jet, the remote wake, mixes and transitions the characteristics of the two other regions. In addition to probing each region separately, we also simulate a co-rotational PSP-jet encounter. In this scenario, the simulated spacecraft hovers over the jet-producing CH, as may occur during the mission's co-rotational phases, sampling each jet region in turn. We estimate that PSP will encounter numerous CH jets over the lifetime of the mission.

Numerous observations have revealed that power-law distributions are ubiquitous in energetic solar processes. Hard X-rays, soft X-rays, extreme ultraviolet radiation, and radio waves all display power-law frequency distributions. Since magnetic reconnection is the driving mechanism for many energetic solar phenomena, it is likely that reconnection events themselves display such power-law distributions. In this work, we perform numerical simulations of the solar corona driven by simple convective motions at the photospheric level. Using temperature changes, current distributions, and Poynting fluxes as proxies for heating, we demonstrate that energetic events occurring in our simulation display power-law frequency distributions, with slopes in good agreement with observations. We suggest that the braiding-associated reconnection in the corona can be understood in terms of a self-organized criticality model driven by convective rotational motions similar to those observed at the photosphere.

We have reported previously on a new method we are developing for using image-based information to improve global coronal magnetic field models. In that work we presented early tests of the method which proved its capability to improve global models based on flawed synoptic magnetograms, given excellent constraints on the field in the model volume. In this follow-up paper we present the results of similar tests given field constraints of a nature that could realistically be obtained from quality white-light coronagraph images of the lower corona. We pay particular attention to difficulties associated with the line-of-sight projection of features outside of the assumed coronagraph image plane, and the effect on the outcome of the optimization of errors in localization of constraints. We find that substantial improvement in the model field can be achieved with this type of constraints, even when magnetic features in the images are located outside of the image plane.

Extreme-ultraviolet and X-ray jets occur frequently in magnetically open coronal holes on the Sun, especially at high solar latitudes. Some of these jets are observed by white-light coronagraphs as they propagate through the outer corona toward the inner heliosphere, and it has been proposed that they give rise to microstreams and torsional Alfv\'{e}n waves detected in situ in the solar wind. To predict and understand the signatures of coronal-hole jets, we have performed a detailed statistical analysis of such a jet simulated with an adaptively refined magnetohydrodynamics model. The results confirm the generation and persistence of three-dimensional, reconnectiondriven magnetic turbulence in the simulation. We calculate the spatial correlations of magnetic fluctuations within the jet and find that they agree best with the M\"{u}ller - Biskamp scaling model including intermittent current sheets of various sizes coupled via hydrodynamic turbulent cascade. The anisotropy of the magnetic fluctuations and the spatial orientation of the current sheets are consistent with an ensemble of nonlinear Alfv\'{e}n waves. These properties also reflect the overall collimated jet structure imposed by the geometry of the reconnecting magnetic field. A comparison with Ulysses observations shows that turbulence in the jet wake is in quantitative agreement with that in the fast solar wind.

The coronal magnetic field directly or indirectly affects a majority of the phenomena studied in space physics. It provides energy for coronal heating, controls the release of coronal mass ejections (CMEs), and drives heliospheric and magnetospheric activity, yet the coronal magnetic field itself has proven difficult to measure. This difficulty has prompted a decades-long effort to develop accurate, timely, models of the field - an effort that continues today. We have developed a method for improving global coronal magnetic field models by incorporating the type of morphological constraints which could be derived from coronal images. Here we report promising initial tests of this approach on two theoretical problems, and discuss opportunities for application.

We present the results from the first ensemble prediction model for major solar flares (M and X classes). The primary aim of this investigation is to explore the construction of an ensemble for an initial prototyping of this new concept. Using the probabilistic forecasts from three models hosted at the Community Coordinated Modeling Center (NASA-GSFC) and the NOAA forecasts, we developed an ensemble forecast by linearly combining the flaring probabilities from all four methods. Performance-based combination weights were calculated using a Monte-Carlo-type algorithm that applies a decision threshold $P_{th}$ to the combined probabilities and maximizing the Heidke Skill Score (HSS). Using the data for 13 recent solar active regions between years 2012 - 2014, we found that linear combination methods can improve the overall probabilistic prediction and improve the categorical prediction for certain values of decision thresholds. Combination weights vary with the applied threshold and none of the tested individual forecasting models seem to provide more accurate predictions than the others for all values of $P_{th}$. According to the maximum values of HSS, a performance-based weights calculated by averaging over the sample, performed similarly to a equally weighted model. The values $P_{th}$ for which the ensemble forecast performs the best are 25 \% for M-class flares and 15 \% for X-class flares. When the human-adjusted probabilities from NOAA are excluded from the ensemble, the ensemble performance in terms of the Heidke score, is reduced.

We study structure and dynamics of the turbulent photospheric magnetic field in the active region NOAA 11158 by characterizing the spatial and temporal scaling properties of the line-of-sight (LOS) component. Using high-resolution high-cadence LOS magnetograms from SDO/HMI, we measured power-law exponents $\alpha$ and $\beta$ describing wavenumber- ($k$) and frequency-domain ($f$) Fourier power spectra, respectively, and investigated their evolution during the passage of the active region through the field of view of HMI. Flaring active region NOAA 11158 produces an average one-dimensional spatial power spectral density that follows approximately a $k^{-2}$ power law -- a spectrum that suggests parallel MHD fluctuations in an anisotropic turbulent medium. In addition, we found that values of $\alpha$ capture systematic changes in the configuration of LOS photospheric magnetic field during flaring activity in the corona. Position-dependent values of the temporal scaling exponent $\beta$ showed that, on average, the core of the active region scales with $\beta >$ 3 surrounded by a diffusive region with an approximately $f^{-2}$-type spectrum. Our results indicate that only about 1 - 3 \% of the studied LOS photospheric magnetic flux displays $\beta\approx\alpha$, implying that Taylor's hypothesis of frozen-in-flow turbulence is typically invalid for this scalar in presence of turbulent photospheric flows. In consequence, both spatial and temporal variations of the plasma and magnetic field must be included in a complete description of the turbulent evolution of active regions.

Using data from STEREO and SOHO spacecraft, we show that temporal organization of energy release events in the quiet solar corona is close to random, in contrast to the clustered behavior of flaring times in solar active regions. The locations of the quiet-Sun events follow the meso- and supergranulation pattern of the underling photosphere. Together with earlier reports of the scale-free event size statistics, our findings suggest that quiet solar regions responsible for bulk coronal heating operate in a driven self-organized critical state, possibly involving long-range Alfv\'{e}nic interactions.

Shortly after the seminal paper {\sl "Self-Organized Criticality: An explanation of 1/f noise"} by Bak, Tang, and Wiesenfeld (1987), the idea has been applied to solar physics, in {\sl "Avalanches and the Distribution of Solar Flares"} by Lu and Hamilton (1991). In the following years, an inspiring cross-fertilization from complexity theory to solar and astrophysics took place, where the SOC concept was initially applied to solar flares, stellar flares, and magnetospheric substorms, and later extended to the radiation belt, the heliosphere, lunar craters, the asteroid belt, the Saturn ring, pulsar glitches, soft X-ray repeaters, blazars, black-hole objects, cosmic rays, and boson clouds. The application of SOC concepts has been performed by numerical cellular automaton simulations, by analytical calculations of statistical (powerlaw-like) distributions based on physical scaling laws, and by observational tests of theoretically predicted size distributions and waiting time distributions. Attempts have been undertaken to import physical models into the numerical SOC toy models, such as the discretization of magneto-hydrodynamics (MHD) processes. The novel applications stimulated also vigorous debates about the discrimination between SOC models, SOC-like, and non-SOC processes, such as phase transitions, turbulence, random-walk diffusion, percolation, branching processes, network theory, chaos theory, fractality, multi-scale, and other complexity phenomena. We review SOC studies from the last 25 years and highlight new trends, open questions, and future challenges, as discussed during two recent ISSI workshops on this theme.

On 20 August 2010 an energetic disturbance triggered large-amplitude longitudinal oscillations in a nearby filament. The triggering mechanism appears to be episodic jets connecting the energetic event with the filament threads. In the present work we analyze this periodic motion in a large fraction of the filament to characterize the underlying physics of the oscillation as well as the filament properties. The results support our previous theoretical conclusions that the restoring force of large-amplitude longitudinal oscillations is solar gravity, and the damping mechanism is the ongoing accumulation of mass onto the oscillating threads. Based on our previous work, we used the fitted parameters to determine the magnitude and radius of curvature of the dipped magnetic field along the filament, as well as the mass accretion rate onto the filament threads. These derived properties are nearly uniform along the filament, indicating a remarkable degree of cohesiveness throughout the filament channel. Moreover, the estimated mass accretion rate implies that the footpoint heating responsible for the thread formation, according to the thermal nonequilibrium model, agrees with previous coronal heating estimates. We estimate the magnitude of the energy released in the nearby event by studying the dynamic response of the filament threads, and discuss the implications of our study for filament structure and heating.

A set of co-aligned high resolution images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) is used to investigate propagating disturbances (PDs) in warm fan loops at the periphery of a non-flaring active region NOAA AR 11082. To measure PD speeds at multiple coronal temperatures, a new data analysis methodology is proposed enabling quantitative description of subvisual coronal motions with low signal-to-noise ratios of the order of 0.1 %. The technique operates with a set of one-dimensional "surfing" signals extracted from position-time plots of several AIA channels through a modified version of Radon transform. The signals are used to evaluate a two-dimensional power spectral density distribution in the frequency - velocity space which exhibits a resonance in the presence of quasi-periodic PDs. By applying this analysis to the same fan loop structures observed in several AIA channels, we found that the traveling velocity of PDs increases with the temperature of the coronal plasma following the square root dependence predicted for the slow mode magneto-acoustic wave which seems to be the dominating wave mode in the studied loop structures. This result extends recent observations by Kiddie et al. (Solar Phys., 2012) to a more general class of fan loop systems not associated with sunspots and demonstrating consistent slow mode activity in up to four AIA channels.

Hot flow anomalies (HFAs) represent a subset of solar wind discontinuities interacting with collisionless bow shocks. They are typically formed when the normal component of motional (convective) electric field points toward the embedded current sheet on at least one of its sides. The core region of an HFA contains hot and highly deflected ion flows and rather low and turbulent magnetic field. In this paper, we report first observations of HFA-like events at Mercury identified over a course of two planetary years. Using data from the orbital phase of the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) mission, we identify a representative ensemble of active current sheets magnetically connected to Mercury's bow shock. We show that some of these events exhibit unambiguous magnetic and particle signatures of HFAs similar to those observed earlier at other planets, and present their key physical characteristics. Our analysis suggests that Mercury's bow shock does not only mediate the flow of supersonic solar wind plasma but also provides conditions for local particle acceleration and heating as predicted by previous numerical simulations. Together with earlier observations of HFA activity at Earth, Venus and Saturn, our results confirm that hot flow anomalies are a common property of planetary bow shocks, and show that the characteristic size of these events is of the order of one planetary radius.

The observed solar activity is believed to be driven by the dissipation of nonpotential magnetic energy injected into the corona by dynamic processes in the photosphere. The enormous range of scales involved in the interaction makes it difficult to track down the photospheric origin of each coronal dissipation event, especially in the presence of complex magnetic topologies. In this paper, we propose an ensemble-based approach for testing the photosphere - corona coupling in a quiet solar region as represented by intermittent activity in SOHO MDI and STEREO EUVI image sets. For properly adjusted detection thresholds corresponding to the same degree of intermittency in the photosphere and corona, the dynamics of the two solar regions is described by the same occurrence probability distributions of energy release events but significantly different geometric properties. We derive a set of scaling relations reconciling the two groups of results and enabling statistical description of coronal dynamics based on photospheric observations. Our analysis suggests that multiscale intermittent dissipation in the corona at spatial scales > 3 Mm is controlled by turbulent photospheric convection. Complex topology of the photospheric network makes this coupling essentially nonlocal and non-deterministic. Our results are in an agreement with the Parker's coupling scenario in which random photospheric shuffling generates marginally stable magnetic discontinuities at the coronal level, but they are also consistent with an impulsive wave heating involving multiscale Alfvenic wave packets and/or MHD turbulent cascade. A back reaction on the photosphere due to coronal magnetic reconfiguration can be a contributing factor.

Multiscale topological complexity of solar magnetic field is among the primary factors controlling energy release in the corona, including associated processes in the photospheric and chromospheric boundaries. We present a new approach for analyzing multiscale behavior of the photospheric magnetic flux underlying this dynamics as depicted by a sequence of high-resolution solar magnetograms. The approach involves two basic processing steps: (1) identification of timing and location of magnetic flux origin and demise events (as defined by DeForest et al., 2007) by tracking spatiotemporal evolution of unipolar and bipolar photospheric regions, and (2) analysis of collective behavior of the detected magnetic events using a generalized version of Grassberger - Procaccia correlation integral algorithm. The scale-free nature of the developed algorithms makes it possible to characterize the dynamics of the photospheric network across a wide range of distances and relaxation times. Three types of photospheric conditions are considered to test the method: a quiet photosphere, a solar active region (NOAA 10365) in a quiescent non-flaring state, and the same active region during a period of M-class flares. The results obtained show (1) the presence of a topologically complex asymmetrically fragmented magnetic network in the quiet photosphere driven by meso- and supergranulation, (2) the formation of non-potential magnetic structures with complex polarity separation lines inside the active region, and (3) statistical signatures of canceling bipolar magnetic structures coinciding with flaring activity in the active region. Each of these effects can represent an unstable magnetic configuration acting as an energy source for coronal dissipation and heating.

We report numerical evidence of a self-organized criticality (SOC) and intermittent turbulence (IT) symbiosis in a resistive magnetohydrodynamic (MHD) current sheet model that includes a local hysteretic switch to capture plasma physical processes outside of MHD that are described in the model as current-dependent resistivity. Results from numerical simulations show scale-free avalanches of magnetic energy dissipation characteristic of SOC, as well as multi-scaling in the velocity field numerically indistinguishable from certain hierarchical turbulence theories. We suggest that SOC and IT may be complementary descriptions of dynamical states realized by driven current sheets -- which occur ubiquitously in astrophysical and space plasmas.

An extended data set of extreme ultraviolet images of the solar corona provided by the SOHO spacecraft are analyzed using statistical methods common to studies of self-organized criticality (SOC) and intermittent turbulence (IT). The data exhibits simultaneous hallmarks of both regimes, namely power law avalanche statistics as well as multiscaling of structure functions for spatial activity. This implies that both SOC and IT may be manifestations of a single complex dynamical process entangling avalanches of magnetic energy dissipation with turbulent particle flows.