We present a theoretical and experimental study of light transport in disordered media with strongly heterogeneous distribution of scatterers formed via non-scattering regions.Step correlations induced by quenched disorder are found to prevent diffusivity from diverging with increasing heterogeneity scale, contrary to expectations from annealed models. Spectral diffusivity is measured for a porous ceramic where nanopores act as scatterers and macropores render their distribution heterogeneous. Results agree well with Monte Carlo simulations and a proposed analytical model.Series of incremental random changes govern the evolution of countless systems around us, from the movement of particles and molecules [1][2][3][4] and wave propagation in disordered media [5,6], to the foraging of animals [7] and spread of disease [8]. By virtue of the central limit theorem, macroscopic evolution of such random walk processes can often be explained in terms of classical diffusion. Although therefore ubiquitous, diffusion is in each particular case determined by unique microscopic mechanisms. These mechanisms are often complex and understanding the onset and speed of diffusion is generally a challenge [9][10][11][12][13][14]. The matter is particularly relevant to research in optics of disordered media, including the study of radiative transfer through planetary atmospheres [15][16][17], optical imaging and spectroscopy in biomedical [18] and material science [19] and, more recently, anomalous diffusion in engineered disordered materials [20][21][22][23][24].Multiple scattering of light is typically viewed as a Poissonian random walk of independent and exponentially-distributed steps. This viewpoint inherently assumes a uniform random distribution of scatterers throughout the medium and results in a well-known expression for the diffusion constant [25], D = v t /3, where v is the average transport velocity for light in the medium and t the transport mean free path. This diffusivity relation, however, breaks down in systems with an heterogeneous distribution of scatterers, such as clouds, biological tissues, porous materials and foams. The reason is two-fold. First, the presence of non-scattering regions in the scattering medium leads to a broader (nonexponential) distribution of step lengths, which, in turn, induces an increase of the diffusivity [26]. Second, the quenched (i.e. spatially frozen) heterogeneity induces step correlations that tend to counteract the increase in diffusivity caused by long steps [21,26]. Due to the complexity of these aspects, understanding of light transport in systems with heterogeneous distribution of scatterers remains rather limited. Important insight has, nonetheless, been reached with the development of generalized transport equations and homogenization theory [27][28][29][30][31][32][33][34][35] as well as probabilistic analysis of random walks [26,36,37]. When it comes to quenched disorder, most works are theoretical and fall within the context of anomalous diffusion [21,23,24,[38][39][40]....