The cosmic background radiation has been observed to deviate from the Planck law expected from a blackbody at $\sim$ 2.7 K at frequencies below $\sim 3$ GHz. We discuss the abundance of the low-energy photons from the perspective of nonequilibrium statistical mechanics. We propose a mechanism of stochastic frequency-diffusion, the counterpart to stochastic acceleration for charged particles in a turbulent plasma, to modify the standard Kompaneets equation. The resulting violation of the Einstein relation allows to take advantage of low-frequency localization and finally leads to photon cooling. The new equation predicts a frequency distribution in agreement with the absolute temperature measurements of the cosmic background radiation down to about 20 MHz, for which we offer here an updated compilation. In that sense, the so called 'space roar' we observe today is interpreted as a nonequilibrium echo of the early universe, and of nonequilibrium conditions in the primordial plasma more specifically.

A model of the Lu-Hamilton kind is applied to the study of critical behavior of the magnetized solar atmosphere. The main novelty is that its driving is done via sources undergoing a diffusion. This mimics the effect of a virtual turbulent substrate forcing the system. The system exhibits power-law statistics not only in the size of the flares, but also in the distribution of the waiting times.

A ``sandpile'' cellular automaton achieves complex temporal correlations, like a $1/f$ spectrum, if the position where it is perturbed diffuses slowly rather than changing completely at random, showing that the spatial correlations of the driving are deeply related to the intermittent activity. Hence, recent arguments excluding the relevance of self-organized criticality in seismicity and in other contexts are inaccurate. As a toy model of single fault evolution, and despite of its simplicity, our automaton uniquely reproduces the scaling form of the broad distributions of waiting times between earthquakes.

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.