In recent years, harmful algal blooms (HABs) have increased in their severity and extent in many parts of the world and pose serious threats to local aquaculture, fisheries, and public health. In many cases, the mechanisms triggering and regulating HAB events remain poorly understood. Using underwater microscopy and Residual Neural Network (ResNet‐18) to taxonomically classify imaged organisms, we developed a daily abundance record of four potentially harmful algae (Akashiwo sanguinea, Chattonella spp., Dinophysis spp., and Lingulodinium polyedra) and major grazer groups (ciliates, copepod nauplii, and copepods) from August 2017 to November 2020 at Scripps Institution of Oceanography pier, a coastal location in the Southern California Bight. Random Forest algorithms were used to identify the optimal combination of environmental and ecological variables that produced the most accurate abundance predictions for each taxon. We developed models with high prediction accuracy for A. sanguinea (R2=0.79±0.06), Chattonella spp. (R2=0.63±0.06), and L. polyedra (R2=0.72±0.08), whereas models for Dinophysis spp. showed lower prediction accuracy (R2=0.24±0.07). Offshore nutricline depth and indices describing climate variability, including El Niño Southern Oscillation, Pacific Decadal Oscillation, and North Pacific Gyre Oscillation, that influence regional‐scale ocean circulation patterns and environmental conditions, were key predictor variables for these HAB taxa. These metrics of regional‐scale processes were generally better predictors of HAB taxa abundances at this coastal location than the in situ environmental measurements. Ciliate abundance was an important predictor of Chattonella and Dinophysis spp., but not of A. sanguinea and L. polyedra. Our findings indicate that combining regional and local environmental factors with microzooplankton populations dynamics can improve real‐time HAB abundance forecasts.