Complex networks are characterized by latent geometries induced by their topology or by the dynamics on the top of them. In the latter case, different network-driven processes induce distinct geometric features that can be captured by adequate metrics. Random walks, a proxy for a broad spectrum of processes, from simple contagion to metastable synchronization and consensus, have been recently used in [Phys. Rev. Lett. 118, 168301 (2017)] to define the class of diffusion geometry and pinpoint the functional mesoscale organization of complex networks from a genuine geometric perspective. Here, we firstly extend this class to families of distinct random walk dynamics -- including local and non-local information -- on the top of multilayer networks -- a paradigm for biological, neural, social, transportation, biological and financial systems -- overcoming limitations such as the presence of isolated nodes and disconnected components, typical of real-world networks. Secondly, we characterize the multilayer diffusion geometry of synthetic and empirical systems, highlighting the role played by different random search dynamics in shaping the geometric features of the corresponding diffusion manifolds.