We present a quantum Monte Carlo study of the "quantum glass" phase of the 2D Bose-Hubbard model with random potentials at filling ρ = 1. In the narrow region between the Mott and superfluid phases the compressibility has the form κ ∼ exp(−b/T α ) + c with α < 1 and c vanishing or very small. Thus, at T = 0 the system is either incompressible (a Mott glass) or nearly incompressible (a Mott-glass-like anomalous Bose glass). At stronger disorder, where a glass reappears from the superfluid, we find a conventional highly compressible Bose glass. On a path connecting these states, away from the superfluid at larger Hubbard repulsion, a change of the disorder strength by only 10% changes the low-temperature compressibility by more than four orders of magnitude, lending support to two types of glass states separated by a phase transition or a sharp cross-over. PACS numbers: 64.70.Tg, 67.85.Hj, 67.10.Fj There are two types of ground states of interacting lattice bosons in the absence of disorder; the superfluid (SF) and the Mott-insulator (MI). In the Bose-Hubbard model (BHM) with repulsive on-site interactions [1,2] an MI state has an integer number of particles per site and there is a gap to states with added or removed particles. The gapless SF can have any filling fraction. These phases and the quantum phase transitions between them are well understood [1][2][3][4][5][6] and have been realized experimentally with ultracold atoms in optical lattices [7,8].If disorder in the form of random site potentials is introduced in the BHM (which can also be accomplished in optical lattices [9,10]) a third state appears-an insulating but gapless quantum glass. This state has been the subject of numerous studies [1][2][3][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27] but many of its properties are still not well understood. Two types of glass states are known; the compressible Bose glass (BG) and the incompressible Mott glass (MG), with the latter commonly believed to appear only at commensurate filling fractions in systems with particle-hole symmetry [18][19][20][26][27][28][29]. The currently prevailing notion is that the glass state in the 2D BHM with random potentials is always of the compressible BG type [20][21][22]25].We here present quantum Monte Carlo (QMC) results for the two-dimensional (2D) site-disordered BHM, showing that there is actually an extended parameter region in which the BG is either replaced by an MG or has an anomalously small (in practice undetectable) compressibility. The system is described by the Hamiltonianwhere ij are nearest neighbors on the square lattice, bsite occupation numbers, and ǫ i random potentials uniformly distributed in the range [−Λ − µ, Λ − µ] about the average chemical potential µ. We study the model using the stochastic series expansion (SSE) QMC method with directed loop updates [30]. We adjust the chemical potential so that the mean filling-fraction ρ = n i = 1 (to within < 10 −5 ) when averaged over sites i, disorder realizations, quantum and thermal ...