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.

Introducing thresholds to analyze time series of emission from the Sun enables a new and simple definition of solar flare events, and their interoccurrence times. Rescaling time by the rate of events, the waiting and quiet time distributions both conform to scaling functions that are independent of the intensity threshold over a wide range. The scaling functions are well described by a two parameter function, with parameters that depend on the phase of the solar cycle. For flares identified according to the current, standard definition, similar behavior is found.

A sequence of bursts observed in an intermittent time series may be caused by a single avalanche, even though these bursts appear as distinct events when noise and/or instrument resolution impose a detection threshold. In the Bak-Tang-Wiesenfeld sandpile, the statistics of quiet times between bursts switches from Poissonian to scale invariant on raising the threshold for detecting instantaneous activity, since each zero-threshold avalanche breaks into a hierarchy of correlated bursts. Calibrating the model with the time resolution of GOES data, qualitative agreement with the inter-occurrence time statistics of solar flares at different intensity thresholds is found.

Networks are paradigms for describing complex biological, social and technological systems. Here I argue that networks provide a coherent framework to construct coarse-grained models for many different physical systems. To elucidate these ideas, I discuss two long-standing problems. The first concerns the structure and dynamics of magnetic fields in the solar corona, as exemplified by sunspots that startled Galileo almost 400 years ago. We discovered that the magnetic structure of the corona embodies a scale free network, with spots at all scales. A network model representing the three-dimensional geometry of magnetic fields, where links rewire and nodes merge when they collide in space, gives quantitative agreement with available data, and suggests new measurements. Seismicity is addressed in terms of relations between events without imposing space-time windows. A metric estimates the correlation between any two earthquakes. Linking strongly correlated pairs, and ignoring pairs with weak correlation organizes the spatio-temporal process into a sparse, directed, weighted network. New scaling laws for seismicity are found. For instance, the aftershock decay rate decreases as 1/t in time up to a correlation time, t[omori]. An estimate from the data gives t[omori] to be about one year for small magnitude 3 earthquakes, about 1400 years for the Landers event, and roughly 26,000 years for the earthquake causing the 2004 Asian tsunami. Our results confirm Kagan's conjecture that aftershocks can rumble on for centuries.

The sun provides an explosive, heavenly example of self-organized criticality. Sudden bursts of intense radiation emanate from rapid rearrangements of the magnetic field network in the corona. Avalanches are triggered by loops of flux that reconnect or snap into lower energy configurations when they are overly stressed. Our recent analysis of observational data reveals that the loops (links) and footpoints (nodes), where they attach on the photosphere, embody a scale free network. The statistics of the avalanches and of the network structure are unified through a simple dynamical model where the avalanches and network co-generate each other into a complex, critical state. This particular example points toward a general dynamical mechanism for self-generation of complex networks.

We propose that the coronal magnetic field, linking concentrations on the photosphere through an interwoven web of flux, embodies a scale-free network. It arises from a self-organized critical dynamics including flux emergence, the diffusion and merging of magnetic concentrations, as well as avalanches of reconnecting flux tubes. Magnetic concentrations such as fragments, pores and sunspots, are `nodes' joined by flux tubes or `links'. The number of links emanating from a node is scale-free. We reanalyze the quiet-Sun data of Close et al and show that the distribution of magnetic concentration strengths is a power law with an index $\gamma= 1.7 \pm 0.3$, over the entire range of the measurement, about $(2-500)\times 10^{17}$ Mx. This distribution is compatible with that for the sizes of active regions reported by Harvey and Schwaan. Thus magnetic concentrations may be scale-free from the smallest measurable fragments to the large active regions. Numerical simulations of a self-organized critical model give the same index $\gamma$, within statistical uncertainty. The exponential distribution of flux tube lengths also agrees quantitatively with results from the model. Calibration with the measured diffusion constant of magnetic concentrations allows us to calculate a flux turnover time in the model to be of order 10 hours and the total solar flux to be of order $10^{23}$Mx, agreeing with observations. We introduce two other statistical quantities to characterize scale-free networks. The probability distribution for the amount of flux connecting a pair of concentrations, and the number of distinct concentrations linked to a given one are predicted to be scale-free, with different indices. Our approach unifies the observation of scale free flare energies with the coronal magnetic field structure.

We show that large, slowly driven systems can evolve to a self-organized critical state where long range temporal correlations between bursts or avalanches produce low frequency $1/f^{\alpha}$ noise. The avalanches can occur instantaneously in the external time scale of the slow drive, and their event statistics are described by power law distributions. A specific example of this behavior is provided by numerical simulations of a deterministic ``sandpile'' model.

The almost perfect uniformity of the cosmic microwave background (CMB) radiation, discovered by Penzias and Wilson in 1965 appears to present clearcut evidence that the universe was uniform and in equilibrium at the decoupling transition when a plasma of protons and electrons condensed into a gas of Hydrogen. COBE indicates that only very small ripples of order 10^{-5} existed at decoupling. Gravity then caused hydrogen to cluster and possibly reheat parts of the universe to form the luminous matter that we observe today. We suggest an alternative scenario, where a spatially intermittent structure of extremely hot matter already existed in an otherwise uniform plasma state at the decoupling transition. The plasma was not in equilibrium but in a very high Reynolds number turbulent state. The sparse bursts would not affect the uniformity of the CMB radiation. Luminous matter originates from localized hot bursts already present in the plasma state prior to decoupling. No reheating, and no exotic matter is needed to get luminous matter.