  This paper is concerned with the processes of spatial propagation and
penetration of turbulence from the regions where it is locally excited into
initially laminar regions. The phenomenon has come to be known as "turbulence
spreading" and witnessed a renewed attention in the literature recently. Here,
we propose a comprehensive theory of turbulence spreading based on fractional
kinetics. We argue that the use of fractional-derivative equations permits a
general approach focussing on fundamentals of the spreading process regardless
of a specific turbulence model and/or specific instability type. The starting
point is the Hamiltonian of the resonant wave-wave interactions, from which a
family of scaling laws for the asymptotic spreading is derived. Both three- and
four-wave interactions are considered. The results span from a subdiffusive
spreading in the parameter range of weak chaos to avalanche propagation in
regimes with population inversion. Attention is paid to how non-ergodicity
introduces weak mixing, memory and intermittency into spreading dynamics, and
how the properties of non-Markovianity and nonlocality emerge from the presence
of islands of regular dynamics in phase space. Also we resolve an existing
question concerning turbulence spillover into gap regions, where the
instability growth is locally suppressed, and show that the spillover occurs
through exponential (Anderson like) localization in case of four-wave
interactions and through an algebraic (weak) localization in case of triad
interactions. In the latter case an inverse-cubic behavior of the spillover
function is found. Wherever relevant, we contrast our findings against the
available observational and numerical evidence, and we also commit ourselves to
establish connections with the models of turbulence spreading proposed
previously.
