Seasonal patterns of marine net primary production (NPP) are crucial for understanding the marine ecosystem and biogeochemical cycles. Uncoupling of seasonal variations between NPP and phytoplankton chlorophyll-a over different areas has attracted much attention. In this study, following a review of previous studies, monthly climatological NPP data from 2003 to 2020, estimated using the Size-fractioned Phytoplankton Pigment Absorption (aph)-based NPP Model (SABPM), were selected to study the seasonal variability of NPP in the South China Sea (SCS). Results showed the spatial differences of NPP seasonality and its departures from climatology in extreme El Niño-Southern Oscillation (ENSO) and Indian Ocean Dipole (IOD) years. Cluster analysis for climatological monthly data identified significant differences of NPP seasonality in five typical regions. In coastal regions along the northern SCS and off eastern Vietnam, NPP exhibited the most obvious seasonal cycle with maximum (minimum) values in summer (winter), attributable mainly to river discharge and summer upwelling. In regions off northwestern Luzon and coast of southern SCS, NPP showed peaks in winter, which were related to strong mixing and upwelling. In northwestern SCS, NPP was high during May–September in phase with sea surface temperature and the primary controlling factors were found to be shallow nutricline depth and wind-driven mixing. Owing to the deep nutricline depth in the central basin, NPP exhibited little seasonal variability; only a weak signal was observed in spring in phase with photosynthetically active radiation. Local dynamics on regulating the nutrient supply and light availability contribute to these regional differences in NPP seasonality, which could also be affected by extreme climate events. The largest anomalies of the NPP seasonal cycle coincide with 2015/2016 ENSO and super IOD in 2020. During these events, enhanced (weakened) westerly winds caused fall (rise) of SLA and increase (decrease) of NPP in coastal regions along the northern SCS and that off eastern Vietnam. Overall, the aph-based model shows a new perspective to study the spatiotemporal variations of NPP in the SCS.