We describe a classification system for daily phytoplankton primary productivity models based on four implicit levels of mathematical integration. Depth-integrated productivity models have appeared in the literature on average once every 2 years over the past four decades. All of these models can be related to a single formulation equating depth-integrated primary production (Sigma PP) to surface phytoplankton biomass (C(surf)), a photoadaptive variable (P(opt)(b)), euphotic depth (Z(eu)), an irradiance-dependent function (F), and daylength (DL), The primary difference between models is the description of F, yet we found that irradiance has a relatively minor effect on variability in Sigma PP. We also found that only a small fraction of variability in Sigma PP can be attributed to vertical variability in phytoplankton biomass or variability in the light-limited slope for photosynthesis. Our results indicate that (1) differences between or within any model category have the potential to improve estimates of Sigma PP by <10%, so long as equivalent parameterizations are used for C(surf) and P(opt)(b), and (2) differences in estimates of global annual primary production are due almost entirely to differences in input biomass fields and estimates of the photoadaptive variable, P(opt)(b), not to fundamental differences between model constructs.