The generalized fractal-diffusive SOC model predicts the probability distribution functions for each parameter as a function of the dimensionality, diffusive spreading exponent, fractal dimension, and type of (coherent/incoherent) radiation process. The waiting time distributions are predicted by the FD-SOC model to follow a power law with a slope of during active and contiguously flaring episodes, while an exponential cutoff is predicted for the time intervals of quiescent periods. This dual regime of the waiting time distribution predict both persistence and memory during the active periods, and stochasticity during the quiescent periods. These predictions provide useful constraints of the physical parameters and underlying scaling laws. Significant deviations from the size distributions predicted by the FD-SOC model could indicate problems with the measurements or data analysis. The generic FD-SOC model is considered to have universal validity and explains the statistics and scaling between SOC parameters but does not reveal the detailed physical mechanism that governs the instabilities and energy dissipation in a particular SOC process.

The fractal nature in avalanching systems with SOC is investigated here for phenomena in the solar photosphere and transition region. In the standard SOC model, the fractal Hausdorff dimension is expected to cover the range of [1, 2], with a mean of for 2-D observations projected in the plane-of-sky, and the range of [2, 3], with a mean of for real-world 3-D structures. Observations of magnetograms and with IRIS reveal four groups: (i) photospheric granulation with a low fractal dimension of ; (ii) transition region plages with a low fractal dimension of ; (iii) sunspots at transition region heights with an average fractal dimension of ; and (iv) active regions at photospheric heights with an average fractal dimension of . Phenomena with a low fractal dimension indicate sparse curvilinear flows, while high fractal dimensions indicate near space-filling flows. Investigating the SOC parameters, we find a good agreement for the event areas and mean radiated fluxes in events in transition region plages.

The size distribution of waiting times are found to have an exponential distribution in the case of a stationary Poissonian process. In reality, however, the waiting time distributions reveal power law-like distribution functions, which can be modeled in terms of non-stationary Poisson processes by a superposition of Poissonian distribution functions with time-varying event rates. We model the time evolution of such waiting time distributions by polynomial, sinusoidal, and Gaussian functions, which have exact analytical solutions in terms of the incomplete Gamma function, as well as in terms of the Pareto type-II approximation, which has a power law slope of , where represents the linear time evolution, or with representing nonlinear growth rates, which have a power law slope of . Our mathematical modeling confirms the existence of significant deviations from ideal power law size distributions (of waiting times), but no correlation or significant interval–size relationship exists, as would be expected for a simple (linear) energy storage-dissipation model.

The occurrence frequency distributions (size distributions) are the most important diagnostics for self-organized criticality systems. There are at least three formats for size distributions: (i) the differential size distribution function, (ii) the cumulative size distribution function, and (iii) the rank-order plot. Each of the three formats (or methods) has at least three ranges of event sizes: (i) a range with statistically incomplete sampling; (ii) an inertial range or power law fitting range with statistically complete sampling; and (iii) a range bordering finite system sizes. Only the intermediate range with power law behavior should be used to determine the power law slope from fitting the observed size distributions. The establishment of power law functions in a given observed size distribution depends crucially on the choice of the fitting range, which should have a logarithmic range of at least 2–3 decades. Often the fitted distribution functions exhibit significant deviations from an ideal power law and can be fitted better with alternative functions, such as log-normal distributions, Pareto type-II distributions, and Weibull distributions.

Among stellar systems, we find many with applications of SOC, such as stellar flares or pulsar glitches. Stellar flares occur mostly in the wavelength ranges of ultraviolet, soft X-rays and UV, and in visible light. A breakthrough in new stellar data was accomplished with the Kepler spacecraft, which allowed unprecedented detections of exoplanets, while the same light curves could be searched for large stellar flares. Exploiting these promising new datasets, one finds that most stellar flare datasets exhibit dominant size distributions that converges to a power law slope of , regardless of the star type. The size distributions of pulsar glitches are mostly found outside of the valid range of the Standard FD-SOC model and thus require a different model. Power law fits are not always superior to fits with the log-normal function or Weibull function. This discrepancy between observed and modeled power law slopes in stellar SOC systems is mostly due to small-number statistics of the samples, incomplete sampling near the lower threshold, and due to ill-defined power law fitting ranges, which can cause significant deviations from ideal power laws.

From the statistics of solar radio bursts, we learn that we can discriminate between three diagnostic regimes: (i) the incoherent regime where the radio burst flux is essentially proportional to the flare volume (with a power law slope of ), as it occurs for gyroemission, gyroresonance emission, gyrosynchrotron emission; (ii) the coherent regime that implies a nonlinear scaling between the radio flux and the flare volume ; as it occurs for the electron beam instability, the loss-cone instability, or maser emission; and (iii) the exponential regime that does not display a power law function, but rather an exponential cutoff as expected for random noise distributions. Thus, the power law slopes offer a useful diagnostic to verify the flux–volume scaling law and to discriminate between coherent and incoherent radio emission processes, as well as to distinguish between SOC processes and non-SOC processes. An additional diagnostic comes from the inertial range of power law fits: SOC-related power law size distributions should extend over multiple decades, while power law ranges of less than one decade are most likely not related to SOC processes.

Can we claim that the dynamics of the solar wind is consistent with a SOC system? Observationally we find that magnetic field and kinetic energy fluctuations measured in the solar wind exhibit power law distributions, which is consistent with a SOC system. What about the driver, instability, and avalanches expected in a SOC system? The driver mechanism is the acceleration of the solar wind in the solar corona itself, a process that basically follows the hydrodynamic model of Parker (1958), and may be additionally complicated by the presence of nonlinear wave–particle interactions, such as ion-cyclotron resonance. Then, the instability threshold, triggering extreme bursts of magnetic field fluctuations, the avalanches of solar wind SOC events, can be caused by dissipation of Alfven waves, onset of turbulence, or by the ion-cyclotron instability. Thus, in principle the generalized SOC concept can be applied to the solar wind, if there is a system-wide threshold for an instability that causes extreme magnetic field fluctuations.

The size distribution of solar energetic particle (SEP) events, which represent a more energetic subset than flare events, is mostly found to follow power law distribution functions, rather than Poissonian random distribution functions. However, the numerical value of the power law slope is generally flatter than the slopes of the flare size distributions in hard X-rays, soft X-rays, and EUV (Hudson 1978), which can be explained in at least four different ways: (i) normal flares and proton flares are produced by two fundamentally different acceleration mechanisms; (ii) proton flares behave differently than normal flares; (iii) the fractal dimensionality of SEP events is different from normal flares; (iv) proton flares are subject to a selection bias toward the most energetic events and thus are not a representative sample of large flares. Nevertheless, the standard fractal-diffusive SOC model can explain the observed slopes of SEP size distributions, but observations reveal deviations from straight power law functions, or broken power law slopes, and thus are not unique and need to be modeled in more detail.

We focus on the statistics of SOC-related solar flare parameters in soft X-ray wavelengths, including their size and waiting time distributions. An early SOC model assumed a linear increase of the energy storage, but this pioneering model is not consistent with the expected correlation between the waiting time interval and the subsequently dissipated energy. The Neupert effect in solar flares implies a correlation between the hard X-ray fluence and the soft X-ray flux, which predicts identical size distributions for these two parameters. Quantifying of thermal flare energies in soft X-ray emitting plasma needs also to include radiative and conductive losses. The intermittency and bursty variability of the solar dynamo implies a nonstationary SOC driver, which yields a universal value for the power law slope of fluxes, but the power law slopes of waiting times vary with the flare rate. While our focus encompasses primarily SOC models, alternative models in terms of MHD turbulence can explain some characteristics of SOC features also, such as size distribution functions, Fourier spectra, and structure functions.

The total 2pN net shifts per orbit and the orbital precessions are calculated as the sum of two contributions: the direct ones due to the 2pN acceleration and the mixed, or indirect, ones caused by the 1pN instantaneous shifts during the orbital revolution. A comparison with other approaches existing in the literature is made.

This chapter explores a core question in astrobiology: what is the future of life on Earth and beyond? The first part describes the cessation of habitable conditions in Earth’s distant future (about a billion years hereafter), and the myriad risks that apparently confront humanity on shorter timescales, ranging from wars and artificial intelligence to asteroid impacts and massive volcanoes. The second segment outlines the possibility of humans migrating to other worlds in the solar system, and the numerous technological and logistical challenges expected to arise during this endeavour. The even more daunting notion of interstellar travel is also touched upon, and the propulsion systems and spacecraft advanced in this regard are sketched. The textbook comes to a close by taking stock of the fates that might await humankind.

Mars has always been one of the most promising targets in the search for current or extinct extraterrestrial life. The chapter commences with a brief summary of Mars’ basic characteristics, before describing its potential for instantaneous habitability (e.g., energy sources, bioessential elements), with an emphasis on the availability of water. This is followed by an exposition of how several aspects of Martian habitability have diminished over time, ranging from extensive atmospheric loss to the shutdown of its dynamo, both of which might have contributed to the emergence of its cold and arid climate today. Nevertheless, some specialised abodes where life may have persisted are touched upon (e.g., deep subsurface). In the last part of the chapter, the contentious history of life detection on Mars – the Viking mission experiments in the 1970s and the meteorite ALH84001 – is reviewed, and forthcoming missions to Mars are surveyed.

Icy worlds with subsurface oceans are potentially among the most common repositories of liquid water in the Universe. Moreover, the solar system is confirmed to host a number of such worlds, notably: Europa, Enceladus, and Titan. Motivated by these considerations, this chapter examines the habitability of icy worlds from a general standpoint. The oceanic properties of Europa, Enceladus, and Titan are reviewed, followed by a simple analysis of the physical conditions in which subsurface oceans may be supported. The pathways for the formation of the building blocks of life, their assembly into polymers, and subsequent delivery to the subsurface ocean are elucidated. The possible constraints on the availability of energy sources and bioessential elements are delineated, as well as the types of organisms and ecosystems that could exist. The chapter concludes by briefly speculating about the trajectories of biological evolution conceivable on icy worlds.