Se caracteriza, cuantifica y confirma la afectación contaminante causada al ecosistema Jose mediante las técnicas estadísticas multivariantes de componentes principales y de conglomerados, en base al comportamiento y distribución de las fracciones granulométricas: arena-limo-arcilla, en relación al contenido de la materia orgánica (C-org, N-org y razón C/N) e inorgánico (metales trazas: Cd, Cu, Cr, Pb y Zn) en sedimentos superficiales de 17 estaciones, bajo la influencia biogenética in situ y operacional de la industria petrolera y petroquímica, aunado al influjo antropogénico de afluentes residuales de origen exógeno divergente, que determina las fuentes que ocasionan el paulatino deterioro del ecosistema, principalmente por metales pesados. Los tres primeros componentes principales exhiben el 74,63% de la varianza total explicada, Se encontraron grupos bien definidos entre la granulometría y las especies químicas así como en las estaciones mediante el análisis de conglomerados. Existe elevada correlación estadística positiva y significativa entre las variables razón C/N-C-org(r=0,89); N-orL,Cu(r=0,75); Zn-Cr(r=0,71); Zn-arcilla (r=0.68); con similar comportamiento entre Cu y C-org (r=0.62); Cu arcilla (r=0.61); Limo-razón C/N(r=0.57); N-org-C-org (r=0.56) y limo- C-org (r=0.53). Un comportamiento opuesto de significación negativa en la distribución de arena-arcilla (r=-0,80); Cu-arena (r=- 0,72); Zn-arena (r=-0,70); Zn-Cd (r=-0,69); limo-arena (r= -0,66) y con tendencia similar C-org- arena (r=-0,56) y razón C/N-arena (r=-0,53). Infiriendo que el Pb sigue un comportamiento independiente, indicando esto el grado de contaminación industrial y doméstica asociada a las fracciones granulométricas y al conferido orgánico.