Use of long‐term data and multivariate ordination techniques to identify environmental factors governing estuarine phytoplankton species dynamics