The dynamics of physico-chemical, nutrient, and chlorophyll-a variables were studied in the bloom and non-bloom locations along the off-Gujarat coastal waters to understand the variability in biogeochemistry using multivariate analytical tests. The dissolved oxygen was significantly lower in the bloom stations (3.89 ± 0.44 mgL-1) than in the non-bloom stations (5.50 ± 0.70mg L-1), due to the biological degradation of organic matter in addition to anaerobic microbial respiration. Nutrients (PO4 and NO3) and Chl-a concentrations were recorded higher in the bloom locations at 0.83 ± 0.21µmol L-1, 4.47 ± 0.69µmol L-1, 4.14 ± 1.49mgm-3, respectively. PO4 and NO3 have shown a significantly higher positive correlation of r = 0.73 and r = 0.69 with Chl-a for bloom data than the non-bloom data. The percentage variance contributed by PC1 and PC2 for both bloom and non-bloom locations were estimated at 52.33%. The variable PO4 explains the highest 24.19% variability in PC1, followed by Chl-a (19.89%). The PO4 triggers the bloom formation and also correlates to the higher concentrations of Chl-a in the bloom locations. The bloom concentration ranges from 9553 to 12,235 trichomes L-1. The bloom intensity has shown a significant positive correlation with Chl-a (r = 0.77), NO3 (r = 0.56), and PO4 (r = 0.30), but a negative correlation was noticed with DO (r = - 0.63) and pH (r = - 0.49). The study also initiates a way forward research investigation on ocean-color technologies to identify and monitor blooms and climate change-driven factors for bloom formation. The occurrence of bloom and its influence on fishery resources and other marine biotas will open many research windows in marine fisheries, oceanography, remote sensing, marine biology, and trophodynamics.