The resistance of bacterial biofilms to antibiotic treatment has been attributed to the emergence of structurally heterogeneous microenvironments containing metabolically inactive cell populations. In this study, we use a three-dimensional individual-based cellular automata model to investigate the influence of nutrient availability and quorum sensing on microbial heterogeneity in growing biofilms. Mature biofilms exhibited at least three structurally distinct strata: a high-volume, homogeneous region sandwiched between two compact sections of high heterogeneity. Cell death occurred preferentially in layers in close proximity to the substratum, resulting in increased heterogeneity in this section of the biofilm; the thickness and heterogeneity of this lowermost layer increased with time, ultimately leading to sloughing. The model predicted the formation of metabolically dormant cellular microniches embedded within faster-growing cell clusters. Biofilms utilizing quorum sensing were more heterogeneous compared to their non-quorum sensing counterparts, and resisted sloughing, featuring a cell-devoid layer of EPS atop the substratum upon which the remainder of the biofilm developed. Overall, our study provides a computational framework to analyze metabolic diversity and heterogeneity of biofilm-associated microorganisms and may pave the way toward gaining further insights into the biophysical mechanisms of antibiotic resistance. © 2017, Society for Mathematical Biology.