Whole-brain computational modelling of emergent brain dynamics has become an increasingly important tool for understanding the systems-level mechanisms that govern large-scale brain activity in health and disease. Neuronal synchronization and brain criticality are key candidates of such mechanisms both for healthy cognitive functions and in constituting mechanistic biomarkers for several brain disorders. Despite significant advances, there is not an abundance of modeling approaches that yield both in-vivo-like critical synchronization dynamics and observables that directly match those obtainable with multi-modal neuroimaging. Here, we advance a framework for hierarchical Kuramoto models where the simplest two-layer model comprises a network of pairwise coupled nodes which each contain a large number of oscillators. Already at two levels of hierarchy, this enables explicit representation of local and inter-areal coupling and observations of emergent multi-scale synchronization dynamics. We show here that the model produces observables of meso- and macro-scale synchronization dynamics that are physiologically plausible when compared to, e.g., magneto- and electroencephalographic (M/EEG) data. We also present here a novel approach for fitting this model with individual experimental data that uses not only functional connectivity but also critical dynamics to converge model parameters that are personalized with individual M/EEG and structural-connectivity data. We posit that brain-dynamics matched personalized models yield a basis for 'digital twins' that may support both clinical applications and basic research on the systems-level mechanisms of cognitive functions. The model and model fitting approach advanced here support both new mechanistic understanding of experimental observations and derivation of predictions for future research.