Despite the elegance of Alan Turing’s reaction-diffusion (RD) model, its relevance for the precise description of morphogenesis in real organisms is largely disputed. We will investigate if subtle color sub-clustering predicted by the theoretical model are detectable in real lizards using hyperspectral imaging and extensive histological analyses. If these predictions are confirmed, these data will show that subtle mesoscopic properties of biological dynamical systems, as well as some of the underlying microscopic features, are captured by simple RD models without integrating the unmanageable profusion of variables at lower scales. In addition, we will investigate if the various stochastic (CA and Lenz-Ising) and RD models of patterning can be generalized to five species, belonging to five divergent lineages, exhibiting largely different adult patterns, and all exhibiting dynamics of post-hatching scale-by-scale color change. Taking advantage of their constitutive or effective spatial discretization, we will compare the respective efficiencies of Lenz-Ising, sCA, 2D-dRD, 2D-cRD, and 3D-cRD mathematical models to predict both actual patterns and their statistical attributes. Furthermore, we will develop a mechanical model to predict folding patterns in the developing skin of crocodiles and combine it with RD to produce an extensive model integrating signaling and mechanics. This biophysical model will constitute a powerful tool for predicting how evolution produces functional phenotypes.