From one pattern into another: analysis of Turing patterns in heterogeneous domains via WKBJ

Pattern formation from homogeneity is well studied, but less is known concerning symmetry-breaking instabilities in heterogeneous media. It is non-trivial to separate observed spatial patterning due to inherent spatial heterogeneity from emergent patterning due to nonlinear instability. We employ WKBJ asymptotics to investigate Turing instabilities for a spatially heterogeneous reaction-diffusion system, and derive conditions for instability which are local versions of the classical Turing conditions. We find that the structure of unstable modes differs substantially from the typical trigonometric functions seen in the spatially homogeneous setting. Modes of different growth rates are localized to different spatial regions. This localization helps explain common amplitude modulations observed in simulations of Turing systems in heterogeneous settings. We numerically demonstrate this theory, giving an illustrative example of the emergent instabilities and the striking complexity arising from spatially heterogeneous reaction-diffusion systems. Our results give insight both into systems driven by exogenous heterogeneity, as well as successive pattern forming processes, noting that most scenarios in biology do not involve symmetry breaking from homogeneity, but instead consist of sequential evolutions of heterogeneous states. The instability mechanism reported here precisely captures such evolution, and extends Turing's original thesis to a far wider and more realistic class of systems.