Cyanotoxins, a kind of poisonous substances produced by cyanobacteria, are responsible for health risks in drinking and recreational water uses. The aim of this study is to improve our previous and successful work about cyanotoxins prediction from some experimental cyanobacteria concentrations in the Trasona reservoir (Asturias, Northern Spain) using the multivariate adaptive regression splines (MARS) technique at a local scale. In fact, this new improvement consists of using not only biological variables, but also the physical–chemical ones. As a result, the coefficient of determination has improved from 0.84 to 0.94, that is to say, more accurate predictive calculations and a better approximation to the real problem were obtained. Finally the agreement of the MARS model with experimental data confirmed the good performance. ► An improved MARS model is built as a predictive model of cyanotoxins presence. ► Cyanobacterial HABs are dangerous for environment and people in fresh waters. ► Biological and physical–chemical variables in this process are studied in depth.