Primary Productivity (PP) has been recently investigated using remote sensing based models over quite limited geographical areas of the Red Sea. This work sheds light on how phytoplankton and primary production would react to the effects of global warming in the extreme environment of the Red Sea and, hence, illuminates how similar regions may behave in the context of climate variability. Our study focuses on using satellite observations to conduct an intercomparison of three net primary production (NPP) models-the VGPM (Vertically Generalized Production Model), the Eppley-VGPM and the CbPM (Carbonbased Production Model)-produced over the Red Sea domain for the 1998-2018 time period. A detailed investigation is conducted using multilinear regression analysis, multivariate visualization and moving averages correlative analysis to uncover the models' responses to various climate factors. Here we use the models' 8day composite and monthly averages compared with satellitebased variables including chlorophyll-a (Chla), mixed layer depth (MLD) and sea surface temperature (SST). Seasonal anomalies of NPP are analyzed against different climate indices, namely, the North Pacific Gyre Oscillation (NPGO), the Multivariate ENSO Index (MEI), the Pacific Decadal Oscillation (PDO), the North Atlantic Oscillation (NAO) and the Dipole Mode Index (DMI). In our study, only the CbPM showed significant correlations with NPGO, MEI and PDO, with disagreements relative to the other two NPP models. This can be attributed to the models' connection to oceanographic and atmospheric parameters, as well as the Manuscript received.