Cancer spheroids offer a valuable experimental model that mimics the complexity and heterogeneity of solid tumors. Characterizing their mechanical response is crucial for understanding tumor development, progression, and drug response. Currently, whole live spheroids are analyzed primarily using image analysis, which is challenging, requires extended incubation times, and has limited imaging depth. Here, we present a new label-free approach for characterizing sub-superficial structures of bladder cancer spheroids and measuring their mechanical response at three distinct stages of cancer progression. We study the microrheological changes induced by aging at the cellular and cluster levels by conducting a multi-physics characterization and modeling approach. We find that spheroids exhibit viscoelastic behavior that can be described by fractional models. We show that spheroids are mechanically heterogeneous, with strong depth and time-dependent variations associated with evolving structural features. Our approach opens new possibilities to study 3Din vitromodels, paving the way for the discovery of novel and more precise procedure in cancer diagnosis based on the use of mechanomarkers.