The impact parameter dependent dipole model is an ideal tool for investigating the spatial structure of the proton. We investigate the incoherent ep cross section in exclusive J/ψ photoproduction as measured by HERA, and find that as |t| increases, the model needs several levels of substructure of gluon density fluctuations in order to describe the measured data well. This substructure is modelled as hotspots within hotspots. We find that three levels of substructure are adequate for a good description of all available data, up to |t| = 30 GeV 2 . We note that the gluonic density fluctuation structure follows a scaling behaviour, such that the logarithms of the number of hotspots and their size fall on a line. We further conclude that a more advanced model containing gluon-gluon correlations needs also will need increased complexity as |t| increases. Our findings in this paper systematically constrains which features such a model must have.