Abstract. Diversity plays a key role in the adaptive capacities of marine ecosystems to environmental changes. However, modeling phytoplankton trait diversity remains challenging due to the strength of the competitive exclusion of sub-optimal phenotypes. Trait diffusion (TD) is a recently developed approach to sustain diversity in plankton models by allowing the evolution of functional traits at ecological timescales. In this study, we present a model for Simulating Plankton Evolution with Adaptive Dynamics (SPEAD), where phytoplankton phenotypes characterized by two traits, nitrogen half-saturation constant and optimal temperature, can mutate at each generation using the TD mechanism. SPEAD does not resolve the different phenotypes as discrete entities, computing instead six aggregate properties: total phytoplankton biomass, mean value of each trait, trait variances, and inter-trait covariance of a single population in a continuous trait space. Therefore SPEAD resolves the dynamics of the population's continuous trait distribution by solving its statistical moments, where the variances of trait values represent the diversity of ecotypes. The ecological model is coupled to a vertically-resolved (1D) physical environment, and therefore the adaptive dynamics of the simulated phytoplankton population are driven by seasonal variations in vertical mixing, nutrient concentration, water temperature, and solar irradiance. The simulated bulk properties are validated by observations from BATS in the Sargasso Sea. We find that moderate mutation rates sustain trait diversity at decadal timescales and soften the almost total inter-trait correlation induced by the environment alone, without reducing the annual primary production or promoting permanently maladapted phenotypes, as occur with high mutation rates. As a way to evaluate the performance of the continuous-trait approximation, we also compare the solutions of SPEAD to the solutions of a classical discrete entities approach, both approaches including TD as a mechanism to sustain trait variance. We only find minor discrepancies between the continuous model SPEAD and the discrete model, the computational cost of SPEAD being lower by two orders of magnitude. Therefore SPEAD should be an ideal eco-evolutionary plankton model to be coupled to a general circulation model (GCM) at the global ocean.