Algal blooms are global phenomena in marine and fresh waters. An understanding of the impacts of climate change and anthropogenic activities on marine algal blooms is essential for protecting marine ecosystems. We have studied the long-term changes in the spatial patterns of algal blooms over the northern Indian Ocean and assessed the influence of the ocean, atmosphere, and land-based processes on them. We have analyzed multi-source satellite data (2003 – 2020) of various parameters such as algal concentration, temperature, salinity, sea level, particulate organic carbon, precipitation, wind speed, and river discharge on an open-source cloud computing platform. The seasonality of variables is assessed using a harmonic model, and their spatial trend is estimated using the non-parametric Mann-Kendall test. The spatial relationship between various geophysical parameters and algal bloom concentration is assessed by the cross-correlation test, and the prediction is done using a multivariate autoregression model. Results suggest that each variable exhibits distinguishable seasonal patterns. Algal bloom duration showed a significant increase over time. However, algal bloom coverage showed a significant decline in the coastal waters in all seasons except post-monsoon. The cross-correlation test unveils that lower temperatures, salinity, and sea level promote algal bloom occurrences in coastal waters rather than open waters. Furthermore, the statistical model with a time lag of 2 months shows better performance for prediction of bloom intensity (r2 = 0.86).