This study investigates the heavy metals (HM) and sulfurs concentration, enrichment level and potential sources in the bottom sediments of Nador lagoon, North-East Morocco. 50 samples of sediments from Nador lagoon heavy metals were analyzed by ICP-MS, while their spatial distribution and source identification were evaluated by Enrichment Factor (EF), QGis Inverse Distances Weighted (IDW) interpolation models, multivariate statistical techniques (Principal Component Analysis - Agglomerative Hierarchical Clustering - Pearson correlation matrix), Positive matrix factorization (PMF model) and grain size/organic matter analysis. The average heavy metal concentrations in the sediments are in the decreasing order: Sr > Ba > V > Zr > Zn > Cr > Rb > La > Cu > Pb > Ni > Ce > Nd > Co > Sc > Nb > Ga > Th > Y > Hf, and the Sulfurs. The results showed that the average concentrations of the S, Sr, Pb, and Nd in bottom sediments were higher than their background contents. The EF values demonstrated that HM contents were enriched by anthropogenic activities, especially for Sr and S, and secondly for Th, V, Ce, Nd, and Pb. The potential HM and sulfurs identified sources showed that: 1) Ga, Y, Hf, Nb, Rb, Sc, Zr, Nb, Co, Zn, Th, Ni, V, Ba, Cr, and Cu were derived industrial activities related to stealing and mining; 2) Wastewater treatment station and lagoon watershed were the main sources of Sr in the lagoon; 3) Zn, Cu, and Pb also linked to sulfurs, mostly origin from agriculture inputs and mining activities; 4) Ce, Cr, and Ni linked to sulfurs, mainly origin form metallurgical activities and solid wastes. The results suggest that heavy metal concentrations were significantly affected anthropogenic activities. Severe enrichment of Sr was found in the lagoon ecosystem with no correlation with carbonates, with four sources of heavy metals apportioned by Geostatistics and Positive matrix factorization. Solid waste, water waste and industrial activities were the main anthropogenic source. Lutite fraction and organic matter show a high accord with heavy metals concentration distribution.