Understanding drivers of phytoplankton biomass in dynamic coastal regions is key to predicting present and future ecosystem functioning. Using a clustering-based method, we objectively determined biophysical provinces in a complex estuarine sea. The Salish Sea contains three major distinct provinces where phytoplankton dynamics are controlled by diverse stratification regimes. Our method is simple to implement and broadly applicable for identifying structure in large model-derived datasets.