There can be significant intraspecific individual-level heterogeneity in the intracellular P of phytoplankton, which can affect the population-level growth rate. Several mechanisms can create this heterogeneity, including phenotypic variability in various physiological functions (e.g., nutrient uptake rate). Here, we use modeling to explore the contribution of various mechanisms to the heterogeneity in phytoplankton grown in a laboratory culture. An agent-based model simulates individual cells and their intracellular P. Heterogeneity is introduced by randomizing parameters (e.g., maximum uptake rate) of daughter cells at division. The model was calibrated to observations of the P quota of individual cells of the centric diatom Thalassiosira pseudonana, which were obtained using synchrotron X-ray fluorescence (SXRF). A number of simulations, with individual mechanisms of heterogeneity turned off, then were performed. Comparison of the coefficient of variation (CV) of these and the baseline simulation (i.e., all mechanisms turned on) provides an estimate of the relative contribution of these mechanisms. The results show that the mechanism with the largest contribution to variability is the parameter characterizing the maximum intracellular P, which, when removed, results in a CV of 0.21 compared to a CV of 0.37 with all mechanisms turned on. This suggests that nutrient/element storage capabilities/mechanisms are important determinants of intrapopulation heterogeneity.