En este documento, proponemos un modelo estocástico de fitoplancton-zooplancton impulsado por ruidos coloreados correlacionados, que contiene tanto toxinas antropogénicas como naturales. Utilizando la transformación de Khasminskii y el método de promediado estocástico, primero transformamos el sistema original en un sistema de difusión de Itô. Posteriormente, derivamos la densidad de probabilidad estacionaria de la ecuación de amplitud de promediado utilizando la correspondiente ecuación de Fokker-Planck-Kolmogorov. Luego, se estudia la estabilidad de la amplitud de promediado y se proporciona la densidad de probabilidad conjunta del sistema original bidimensional. Finalmente, los resultados teóricos se verifican mediante simulaciones numéricas, y se ilustran aún más los efectos de las características del ruido y las toxinas en la dinámica del sistema.