Plankton influences biogeochemical and ecosystem processes, such as sequestration of atmospheric CO2, carbon export to the ocean floor, and the productivity of higher trophic levels. Body size is a proxy for many plankton functional traits, and one means of analyzing its community structure is through the distribution of biovolume across size classes (the size spectrum). To understand how climate forcing can affect plankton communities, we assessed the size spectra in the historical simulations of seven Earth System Models (ESMs) included in the 6th Coupled Model Intercomparison Project (CMIP6) and analyzed projected changes under a high emissions scenario (SSP5-8.5). We compared the historical estimates with the Pelagic Size Structure database (PSSdb), a novel size structure dataset from imaging systems. The median slope from models ranged from -1.66 to -1.07, with shallower slopes from this range falling near both the theoretical expectation and PSSdb observations (-1.05), with variations around the median representing differences in the total biovolume distribution across plankton functional groups. Consistent with the observations, most ESMs show steeper slopes and lower biovolume in oligotrophic subtropical gyres compared to productive ocean regions. There was a lack of agreement between models and observations in the size spectra seasonal cycle, possibly stemming from missing model processes and incomplete sampling. Despite these caveats, the size spectra from ESMs presented here, and their evaluation with PSSdb, provides insights on how climate change will affect ecological processes in the plankton, and highlights areas of improvement in model development and imaging data coverage.