  Motivated by recent optical lattice experiments [J.-y. Choi et al., Science
352, 1547 (2016)], we study the dynamics of strongly interacting bosons in the
presence of disorder in two dimensions. We show that Gutzwiller mean-field
theory (GMFT) captures the main experimental observations, which are a result
of the competition between disorder and interactions. Our findings highlight
the difficulty in distinguishing glassy dynamics, which can be captured by
GMFT, and many-body localization, which cannot be captured by GMFT, and
indicate the need for further experimental studies of this system.
