Assessment of Multiple Trace Metal Fluxes in a Semi-Arid Watershed Containing Mine Tailing, Using a Multiple Tool Approach (Zaida Mine, Upper Moulouya Watershed, Morocco)

: Few studies have quantified the complex flux of trace metals from mine tailings to rivers through water erosion, especially in the semi-arid region of North Morocco (Zaida mine) where soil erosion is a severe issue. This study applies (i) methods to understand and estimate the complex flux of trace metals from mine tailings to rivers, using the RUSLE model combined with the concentration of trace metals in the soil and additionally (ii) pollution indices and statistical analyses to assess the sediment contamination by Cd, Cu, Pb, and Zn. Our study revealed that the basin has a low erosion rate, with an average of 9.1 t/ha/yr. Moreover, the soil contamination is particularly high at the north of the mine tailings, as prevailing winds disperse particles across the basin. The assessment of the sediments indicated that Pb is the main contaminant, with concentrations exceeding 200 mg/kg specifically downstream of the tailings. This study also identified high a concentration of trace elements 14 km away from the tailings alongside the Moulouya river, due to the specific hydrological transport patterns in the area. This research contributes to a better understanding of the transport and fate of the trace metals in mining areas. It proposes a replicable method that can be applied in other regions to assess the contamination flows and thereby assist water resource management.


Introduction
Soil is one of the most important resources for human development (including food production while taking into account climate changes) [1].Topsoil is currently considered as a non-renewable resource that is subject to intense degradation worldwide [2].
Mediterranean areas are strongly impacted by soil erosion due to their hilly topography and prolonged periods of dryness followed by irregular, stormy rainfall [3].Several researchers have studied the water erosion in Mediterranean regions.They used different methodologies, including, among others, radioactive isotope dating and bathymetry [4,5].Other researchers have chosen to estimate soil erosion using modelling approaches such as USLE, SWAT, STREAM, etc. [6][7][8].The Revised Universal Soil Loss Equation (RUSLE) is one of the most widely used empirical models for predicting the average annual soil loss worldwide [9].RUSLE is a combination of five key factors for erosion, such as rainfall erosivity R, slope length and steepness LS, soil erodibility K, soil management C, and conservation practices P. Its strength lies in its ability to provide good results despite its limited use of input data compared to other mathematical models [10].
Soil erosion has a negative impact on water resources.Water erosion carries a large number of fertilizers, pesticides, and trace metals (TMs) to the surface and groundwater.This leads in turn to the pollution of water bodies [11].These processes affect agricultural soils, disrupt aquatic ecosystems, and present challenges for the sustainable management of water resources.
Human activities such as agriculture, deforestation, and mining activities make corresponding soils more vulnerable to excessive soil erosion [12].Specifically, opencast mining is a significant contributor to erosion [13].The soil removed during this process is often stored into hills, forming mining tailings [14].In general, mining tailings are a source of environmental contamination due to their high concentrations in trace metals.These sites are subject to erosion, particularly during periods of high water.The trace metals present in the tailings are then gradually released into the environment, mainly via the fine particles bonded to these metals [15,16].This is the case of the Zeïda mine, located in the Upper Moulouya, which started operating in 1972 and closed in 1986 [17].Mining was mainly conducted in quarries, and the cerussite (PbCO 3 ) ore was processed on site via gravimetry and flotation.This process generating around 12 Mt of tailings abandoned without any rehabilitation [18].These tailings are a significant source of pollution for the river and surrounding soils [19,20].
Table S1 contains a summary of the studies conducted on the transport of trace metals from the basin to the river.These studies have focused either on measuring the trace metal flows from mining areas or on assessing the sediment contamination by mine tailings (assessment of sediment, water, soil, etc.).However, none of them (Cf. the table in Table S1) have simultaneously integrated these two perspectives.
Our study, on the other hand, aims at (i) estimating the trace metal load linked to soil erosion.This involves quantifying the amount of trace metals that are mobilized from the soil surface due to erosive forces, particularly in the areas affected by anthropogenic activities like mining, and (ii) assess the resulting contamination of sediments in a semiarid climatic context.This assessment is important for determining how mobilized trace metals accumulate in riverbeds and downstream of the basin.To achieve these objectives, we used a spatial superposition method based on the RUSLE model and data on trace metal concentrations in the soil.We also evaluated the level of sediment contamination in the river using contamination indices.The river acts as a sink for soil particles eroded from contaminated lands, where trace metal concentrations (in sediments) can exceed the threshold limits.These high levels of trace metals' concentration are potentially hazardous to the river ecosystem, mainly its sediments.
In conclusion, our study makes a significant contribution by filling this gap in the literature, by integrating the assessments of metal fluxes through erosion with those of environmental impacts, thus providing a better understanding of the metal transport processes in semi-arid watersheds.

Study area
The abandoned Zeïda mine is located in the Upper Moulouya watershed in northeastern Morocco, as presented in Figure 1.The basin covers an area of 1218 km 2 , with altitudes ranging from 1313 to 2401 m.The Upper Moulouya forms the southwestern part of the eastern Meseta.It is bordered to the southwest by the High Atlas and to the northwest by the Middle Atlas [21].Geologically, the basin consists of a Paleozoic basement composed mainly of granite.It is covered in discordance by Triassic formations (arkoses, sandstones, conglomerates, and argillite's) and Jurassic formations (limestone and dolomite).The mineralization found in Zeïda is a stratiform deposit, located in Triassic arkoses and microconglomerates [22].The main mineral species exploited at Zaida mine include well-crystallized cerussite (PbCO 3 ) and galena (PbS) [23].
dolomite).The mineralization found in Zeïda is a stratiform deposit, located in Triassic arkoses and microconglomerates [22].The main mineral species exploited at Zaida mine include well-crystallized cerussite (PbCO3) and galena (PbS) [23].The climatic conditions in the basin are semi-arid.The annual precipitation ranges from 100 to 400 mm, and the average annual temperatures vary between 12 and 14 °C.The rainfall in the basin is usually short and violent, which emphasizes the erosion and leads to severe inundation [24].The Moulouya, which extends for around 520 km, is the main river in the basin.The flow of the Moulouya in the Zaida station is characterized by The climatic conditions in the basin are semi-arid.The annual precipitation ranges from 100 to 400 mm, and the average annual temperatures vary between 12 and 14 • C. The rainfall in the basin is usually short and violent, which emphasizes the erosion and leads to severe inundation [24].The Moulouya, which extends for around 520 km, is the main river in the basin.The flow of the Moulouya in the Zaida station is characterized by high floods during the rainy periods and low water levels in the summer [25].In addition to the Ansegmir river, the Moulouya is the main tributary of the Hassan II dam.The dam has a capacity of 400 Mm 3 dedicated to drinking water supply and irrigation [26].

Data Source and Pre-Processing
The data used in this study are derived from various origins.To ensure a clear representation, they have been categorized into two types: primary and secondary data.Table 1 provides a summary of the different datasets and their respective sources.

Primary data
To measure the concentrations of trace metals (Cd, Cu, Pb, and Zn) and major elements (Na, K, Ca, Mg, and Al) in the river sediments, a sampling campaign was carried out in March 2022.A total of 10 sediment samples (1-10 cm depth from the surface) were collected (Figure 1).This depth reflects the recently eroded material deposited on the riverbed.These deposits illustrate the risk of contamination of the sediment by loads coming from the basin.Samples were taken from the main Moulouya River (MS1, MS2, MS3, MS4, MS5, MS6, and MS7), its tributary Ansegmir (AS1 and AS2), and a reference sample RF located 35 km from the mining area.The stations were located to cover the entire catchment area (Cf Figure 1), positioned upstream and downstream of the mining area and alongside the major tributaries.This placement allows a precise assessment of the contaminant inputs from the catchment.
For each sampling location, 1 kg of sediment was collected using a mini shovel, then preserved in polyethylene bags at a temperature of 4 • C. In the laboratory, samples were dried in ambient air, ground into powder, and sieved through a 0.25 mm mesh.
The pH of the soil was measured in a soil-solution mixture with a 1:1 ratio, using a pH meter (Hach HQ4300, Hach Company, Berlin, Germany) [28].The soil conductivity was measured using an electronic conductivity meter (Hach HQ4300, Hach Company, Berlin, Germany)), following the method of [29].Meanwhile, the CaCO 3 was measured by using the method of Bernard [30].The samples for the trace metal and major element analysis were digested according to the three acid protocol [31].For this, 1 g of dry matter was placed in a pre-cleaned flask.A mixture of nitric acid (68%), hydrochloric acid (35%), and hydrogen peroxide was added to the flask.The blend was then placed on a sand-bath hot plate to boil lightly for digestion.In the final step, the mixture was diluted with 100 mL of ultra-pure water and filtered through an ashless cellulose membrane.The filtrate was then sent to the Regional University Interface Centre (CURI) for analysis via Inductively Coupled Plasma-Atomic Emission Spectrometry HORIBA Jobin Yvon (ICP-AES).Blank samples were also prepared with ultra-pure water, following the same protocol and under the same conditions.
ICP-AES analyses were carried out for 11 trace metals (As, Ca, Cd, Co, Cr, Cu, Ni, Pb, Si, Ti, and Zn).However, in this study, only the trace metals, Pb, Zn, Cu, and Cd, were retained.The concentrations of the other metals were below the detection threshold of the ICP-AES.
Secondary data In addition to the sampling campaign already conducted, supplementary data were carefully collected from secondary sources.These included trace metal concentrations in the soil around the mine tailings, rainfall, digital elevation models, satellite imagery, and soil texture.
The concentrations of trace elements (Pb, Cu, Zn, Cd) in the soil were obtained from the study of [27].A total of 51 soil samples located around the mine tailings were collected and analyzed via ICP-AES (for more details, see [27]).The various trace metal concentrations, together with the coordinates of the sampling points from [27], were collected and used in this study.
Several preliminary data processing operations were carried out in the study.The data on the metal concentrations in the soils and sediments were normalized.The satellite images and maps at different resolutions were matched using techniques such as resampling, geometric correction, and alignment.

Soil Erosion Estimation via the RUSLE Model
The average annual rate of soil erosion in the watershed was estimated via the RUSLE model [32] using Equation (1).
where A is the average annual soil loss (t•ha −1 •yr −1 ); R is the rainfall erosivity factor (MJ•mm•ha −1 •h −1 •yr −1 ); K is the soil erodibility factor (t•ha•h•ha −1 •MJ −1 •mm −1 ); L is the slope length factor; S is the slope steepness factor; C is the cover-management practice factor; and P is the support practice factor.L, S, C, and P are all dimensionless.

The Rainfall Erosivity Factor (R)
The rainfall erosion factor estimates the precipitation erosion capacity, considering the features like raindrop characteristics and rainfall properties [33].In this study, to estimate the R factor, we used the Modified Fournier Index (MFI) as calculated by Equation (2) and developed by [34].The studies conducted in the Mediterranean regions confirm the efficacy of MFI in estimating the R factor compared to alternative indices such as the Precipitation Concentration Index (PCI), and the Seasonality Index (SI) [35].
where R represents the annual rainfall erosivity factor, pi is the cumulated rainfall of the month i, and P is the cumulated annual rainfall.The Kriging interpolation method was implemented to generate a map (Figure 2a) representing the spatial distribution of the rainfall erosivity in the Zaïda watershed.
where R represents the annual rainfall erosivity factor, pi is the cumulated rainfall of the month i, and P is the cumulated annual rainfall.
The Kriging interpolation method was implemented to generate a map (Figure 2a) representing the spatial distribution of the rainfall erosivity in the Zaïda watershed.

The Topographical Factor (LS)
The topographical factor (LS) is a ratio of the soil loss under specific conditions that occur at a "standard" slope steepness of 9% and a slope length of 22.6 m [36].In this study, the topographical factor was calculated according to Equations ( 3) and (4) developed by [37].The LS factor was calculated through QGIS 3.28 software.

The Topographical Factor (LS)
The topographical factor (LS) is a ratio of the soil loss under specific conditions that occur at a "standard" slope steepness of 9% and a slope length of 22.6 m [36].In this study, the topographical factor was calculated according to Equations ( 3) and (4) developed by [37].The LS factor was calculated through QGIS 3.28 software.
where L is the slope length factor, S is the slope steepness factor, A i,j is the contribution area at the inlet mesh in (m2), D is the size of the mesh in (m), x i, j is the sum of the sine and cosine orientations, m represents the slope length exponent, b is the ratio between rill and inter-rill erosion, and β is the slope angle in degrees.

The Soil Erodibility Factor (K)
The erodibility factor (k) represents the soil loss rate affected by physical soil properties such as soil structure, soil texture, organic content, and hydraulic conductivity [38].The erodibility factor was determined using equations (5, 6, 7, 8, and 9).These equations were developed by [39].

The Cover Management Factor (C)
The (C) factor is an estimate of the ratio between the soil loss from land cultivated under specific conditions and the corresponding loss from continuous clean-tilled fallow [40].The (C) factor was determined using Equations (10) and (11) developed by [41], using an approach based on the Normalized Difference Vegetation Index (NDVI).
where N IR is near infra-red band (band 5), and R is red band (band 4).The values chosen for α and β are 2 and 1, respectively, as they provide good results [41].

The Support Practice Factor (P)
Due to the lack of detailed data on the soil conservation measures applied in the Zaïda watershed, the factor (P) was kept constant at a value of 1.
To calculate the erosion in the 1218 km 2 Zaïda catchment, the RUSLE model was applied using a spatial discretization per 30 m 2 grid cell.The calculations were carried out at the annual temporal resolution using ArcMap 10.8.1 and RStudio 2023.06.1 + 524 software.

Trace Elements' Load Due to the Soil Erosion Model
To calculate the amount of trace metals transported with the eroded soil, Equation (12) proposed by [42] was used.
where M TM is the particulate trace metal load due to erosion in (g•ha −1 •yr −1 ), C TM is the trace metal concentration in soil (mg•kg −1 ), A is the average annual soil loss in (t•ha −1 •yr −1 ) calculated by RUSLE, SDR is the sediment delivery ratio, and E f is the trace metal enrichment factor.
To obtain a complete spatial representation of trace metals' concentration in soils (C TM ), we applied the inverse distance spatial interpolation method to the available data [43].
The SDR is the ratio between the sediment delivery at the watershed outlet and the total erosion in the entire watershed.In this study, the SDR was determined using Equation ( 13) developed by [44].
where A represents the drainage area in mi 2 .
During erosion episodes, small particles are the easiest to be carried away, leading to a phenomenon of particle selection via erosion, especially for clays.Trace metals mainly attached to the fine fraction of the soil are more concentrated in the sediments than in the original soil [45].Therefore, an enrichment factor is calculated using Equation ( 14) developed by [46].
where Sed is the amount of eroded soil in Kg•ha −1 .

Sediment Assessment 2.4.1. Geo-Accumulation Index
To determine the contamination levels in the sediment, we used the geo-accumulation index (Igeo) developed by [47].It is calculated by Equation ( 15) as follows: where C n represents the concentration of the n metal in the sediment, and B n represents the geochemical background value of the n metal.Given the absence of a geochemical background for the study area, reference sediment concentrations were adopted.Table S2 shows the classification of metal contamination levels according to the geo-accumulation index (Igeo) [48].

Pollutant Load Index (PLI)
The PLI is a useful tool to assess the multi-metallic contamination of a single sample [49].It is calculated using the contamination factor (CF) according to Equations ( 16) and ( 17) developed by Tomlinson [49]: where Ci and C i0 are the concentration of the ith trace element in the collected sample and in the geochemical background, respectively.Based on the obtained results, the soil can be classified as unpolluted when PLI ≈ 1 or polluted when PLI > 1.

Sediment Quality Guidelines
To determine the potential risks posed by trace elements on the benthic ecosystem, a comparison is proposed.This involves comparing the Probable Effect Concentration (PEC) and Threshold Effect Concentration (TEC) with the trace element concentrations found in samples [50].The TEC refers to levels below which adverse effects on benthic organisms are not expected to occur.In addition, the PEC indicates the concentrations above which adverse effects are expected.

RUSLE Factors
The spatial distribution maps of the four factors influencing erosion are illustrated in Figure 2.These factors were calculated using ArcGIS and R Studio.
The R factor results value ranges from 31 to 58, with an average of 38 MJ 2a).The data analysis revealed that the highest R values were calculated during the months of November to January.This period of the year is characterized by intense thunderstorm events leading to significant erosion.The lowest values are located in the centre and east of the basin, while the highest values are calculated in the north of the basin, where the Middle Atlas Mountains dominate.
The LS factor in the watershed ranges from 0.03 to 25, as shown in Figure 2b.The highest values of this factor are found mainly in the northern part of the watershed, characterized by hilly slopes and rugged topography.In contrast, most LS factor values of the watershed fall in the range of 0.03 to 4. This distribution of LS factor values is due to the geomorphological nature of the watershed, which is essentially composed of the Oriental Meseta.This region is characterized by uplands surrounded by medium-altitude mountain ranges [21].In these uplands, the slopes are less marked, leading to relatively low values of the LS factor.
The spatial distributions of the K factor are shown in Figure 2c.The erodibility of the soils in the watershed is classified into four categories according to soil type.In the northwestern part of the watershed, the K factor values are the lowest, fluctuating between 0.01 and 0.18.This is due to the presence of an important amount of organic matter attributable to the forest cover [51].Furthermore, the mining areas exhibit the highest K factor values, which range from 0.22 to 0.27.
These high values for the erodibility factor are due to the alteration of the soil resulting from anthropogenic mining activities [52].In fact, the crushing and grinding of ores increase the amount of silt and clay in the tailings, contrasting with the surrounding soils, which are mainly composed of sand (Figure S3).
The average C factor values range from 0.02 to 0.95 as shown in Figure 2d.The spatial distribution of the values reveals that the highest values are observed in bare lands and grasslands.These areas are characterized by a lack of dense vegetation coverage, which contributes to the accentuation of the erosive processes.In contrast, the lowest values are calculated in agricultural areas.This is due to the region's reputation for arboriculture and interplant cropping, practices that preserve the soil and maintain the vegetation cover throughout the year [53].

Distribution and Classification of Soil Erosion
The examination of the soil erosion in the watershed of the Zaida mine reveals a significant variability.The erosion rate ranges from 0.01 to 158 (t•ha −1 •yr −1 ), with an average annual soil loss of 9.1 (t•ha −1 •yr −1 ).The statistical data on soil loss in the watershed are given in Figure 3.The distribution of soil erosion was categorized into five distinct classes: <7, 7-15, 15-30, 30-70, 70-160, as illustrated in Figure 4.The statistical results indicate that the low erosion class covers 64% of the total basin area.However, its contribution to the total erosion remains low, representing only 24% of the total erosion.In contrast, the classes associated with high and very high erosion rates occupy only 12% and 5% of the total surface area, respectively, but despite this, they contribute to 55% of the total erosion in the basin.
The assessment of soil erosion according to the land use indicates that bare soil, urban areas, and agricultural zones are mainly associated with erosion levels classified as low.Meanwhile, high and very high erosion classes are more frequent in steppe and forest areas, respectively (further information on this observation is provided in the following section).Figure 3 shows that the steppe is the major contributor to soil erosion, accounting for 52% of the total erosion.Moreover, the forest, contributes to 22% of the total soil erosion, although it only covers 10% of the total surface area of the basin.
The analysis of the mining area was conducted within the boundaries of the tailings zone, which is a significant source of metallic pollution in the watershed.The erosion rates in this area vary from 0.16 to 46 (t•ha -1 • yr -1 ), with an average of around 10.6 (t•ha -1 • yr -1 ).This value exceeds the average erosion rate of 4.2 (t•ha -1 • yr -1 ) observed in bare soils, including the mine tailings.Additionally, it surpasses the watershed's average erosion rate of 9.1 (t•ha -1 • yr -1 ).These observations highlight the considerable erosive potential associated with the mining area, which raises a serious environmental issue.
The observation of erosion rates presents some unusual results.The erosion averages estimated in the forest and steppe environments exceed those estimated in the bare soil areas.This occurs even though the bare soils are devoid of vegetation that usually protects the soil.Indeed, research has revealed that increasing the percentage of vegetation cover helps dissipate the kinetic energy of water drops, thereby reducing the speed of runoff [54].This situation can be explained by the degradation of vegetation cover being in forest and steppe areas, leaving a large part of the soil unprotected.Indeed, forests in the basin are characterized by sparse vegetation that does not completely cover the soil, contrary to forests in temperate or equatorial areas, which offer dense and continuous coverage [55].Furthermore, an analysis of the factors influencing soil erosion, using Pearson's correlation, is shown in Table 2.The result indicates that the soil erosion in the region is strongly linked to the topographical LS factor, with a correlation of 0.94.The correlation between the LS factor and the soil erosion clarifies the results observed in the watershed.The topographical factor is highly dominant in the forest and steppe areas, whereas it is less marked in the bare land areas.This difference explains the average erosion in forests and steppes even in the presence of vegetation cover.The statistical results indicate that the low erosion class covers 64% of the total basin area.However, its contribution to the total erosion remains low, representing only 24% of the total erosion.In contrast, the classes associated with high and very high erosion rates occupy only 12% and 5% of the total surface area, respectively, but despite this, they contribute to 55% of the total erosion in the basin.
The assessment of soil erosion according to the land use indicates that bare soil, urban areas, and agricultural zones are mainly associated with erosion levels classified as low.Meanwhile, high and very high erosion classes are more frequent in steppe and forest areas, respectively (further information on this observation is provided in the following section).Figure 3 shows that the steppe is the major contributor to soil erosion, accounting for 52% of the total erosion.Moreover, the forest, contributes to 22% of the total soil erosion, although it only covers 10% of the total surface area of the basin.
The analysis of the mining area was conducted within the boundaries of the tailings zone, which is a significant source of metallic pollution in the watershed.The erosion rates in this area vary from 0.16 to 46 (t•ha −1 •yr −1 ), with an average of around 10.6 (t•ha −1 •yr −1 ).This value exceeds the average erosion rate of 4.2 (t•ha −1 •yr −1 ) observed in bare soils, including the mine tailings.Additionally, it surpasses the watershed's average erosion rate of 9.1 (t•ha −1 •yr −1 ).These observations highlight the considerable erosive potential associated with the mining area, which raises a serious environmental issue.
The observation of erosion rates presents some unusual results.The erosion averages estimated in the forest and steppe environments exceed those estimated in the bare soil areas.This occurs even though the bare soils are devoid of vegetation that usually protects the soil.Indeed, research has revealed that increasing the percentage of vegetation cover helps dissipate the kinetic energy of water drops, thereby reducing the speed of runoff [54].This situation can be explained by the degradation of vegetation cover being in forest and steppe areas, leaving a large part of the soil unprotected.Indeed, forests in the basin are characterized by sparse vegetation that does not completely cover the soil, contrary to forests in temperate or equatorial areas, which offer dense and continuous coverage [55].Furthermore, an analysis of the factors influencing soil erosion, using Pearson's correlation, is shown in Table 2.The result indicates that the soil erosion in the region is strongly linked to the topographical LS factor, with a correlation of 0.94.The correlation between the LS factor and the soil erosion clarifies the results observed in the watershed.The topographical factor is highly dominant in the forest and steppe areas, whereas it is less marked in the bare land areas.This difference explains the average erosion in forests and steppes even in the presence of vegetation cover.The statistical data from the trace metal analyses in 51 samples are summarized in Table 3.The average metal concentrations follow the decreasing order, Pb > Zn > Cu > Cd.In fact, cerussite was the main mineral species extracted in the mine, composed mainly of lead carbonate.This composition explains the high concentrations of lead contamination in the soils.a Cited by [56]; b World Health Organization (cited by [57]).
A comparison between the average trace metal concentrations and the thresholds established by [56] highlights the contamination of the soil surrounding the tailing (Table 3).Furthermore, among the 51 soil samples analyzed, over 23% exceeded the limits set by the World Health Organization (WHO) for Pb, 47% for Zn, 15% for Cu, and 37% for Cd.These findings suggest that the inappropriate management of mine tailings after extraction activities have stopped has significantly polluted the soil with trace metals.The spatial distribution of trace metals in the soils of the Upper Moulouya catchment was analyzed using digital mapping techniques, based on the ordinary kriging interpolation.The results indicate that the presence of trace metals is mainly localized around the mine tailings, with a significant decrease in concentrations occurring 4 to 6 km away.
In the context of arid climates, research has shown that wind transport is the main factor influencing the transport of fine particles from mine tailings to the surrounding soils [49].The results obtained in the neighbourhood of the Zaida mine indicate a high dispersion of contamination towards the eastern, northeastern, and southwestern sectors of the tailings.This dispersion is correlated with the prevailing wind directions observed in the Zaida watershed [31].These results highlight the importance of wind transport in the propagation of metal pollution around the Zaida mine watershed.

Spatial Distribution of the Trace Metal Loads Caused by Soil Erosion
The spatial distribution of trace metal loads reveals similar trends among the various metals studied (Pb, Zn, Cu, Cd).However, we paid particular attention to the distribution of Pb.Pb shows the highest concentrations in soils, thus representing a major source of concern in terms of contamination in the catchment.
The Pb load distribution model presents a different pattern compared to that of soil erosion (Figure 5).In fact, soil erosion in the basin is mainly concentrated in areas with a high topography.However, the distribution of the Pb load is mainly located around the mine tailings.The rate of water erosion of Pb in the watershed shows a significant variation, ranging from 0.4 to over 5000 g/ha/yr, with an annual average flux of 166 g/ha.The highest Pb transport rates were observed in the eastern region of the mine tailings, over a distance of 2 km.The results also show that soils with an average annual Pb transport of 200 g/ha extended up to 6 km from the mine tailings, mainly along the Moulouya riverbed.This wide distribution along the river promotes the migration of Pb into the aquatic system during storm events.
The soil erosion and Pb load according to different types of land occupation are shown in Table S4.The land occupation in the watershed is classified into five categories according to their increasing surface area: steppe > bare land > forest > cropland > urban.This hierarchy based on surface area has no influence on the Pb transport rates.In fact, the distribution of the Pb transport rate to the river as a function of land use depends mainly on the nearness of the tailings to the land use.Bare land, which contains mine tailings, is the most contaminated land via metal pollution.These areas are also the main sources of Pb release, with a total migration rate of 9803 kg/yr.They are followed by the steppes, which account for a total transport of 2941 kg/yr of Pb.As regards the urban class, the city of Zaida is located 1.5 km from the mine tailings and is directly affected by the pollution that they generate.However, because of their small surface area, they only contribute 49 kg/yr to the Pb migration rate.The soil erosion and Pb load according to different types of land occupation are shown in Table S4.The land occupation in the watershed is classified into five categories according to their increasing surface area: steppe > bare land > forest > cropland > urban.This hierarchy based on surface area has no influence on the Pb transport rates.In fact, the distribution of the Pb transport rate to the river as a function of land use depends mainly on the nearness of the tailings to the land use.Bare land, which contains mine tailings, is the most contaminated land via metal pollution.These areas are also the main sources of Pb release, with a total migration rate of 9803 kg/yr.They are followed by the steppes, which account for a total transport of 2941 kg/yr of Pb.As regards the urban class, the city of Zaida is located 1.5 km from the mine tailings and is directly affected by the pollution that they generate.However, because of their small surface area, they only contribute 49 kg/yr to the Pb migration rate.

Physiochemical Properties of Sediment
The physiochemical properties of the sediments (pH, electrical conductivity, and carbonate content) are shown in Table S5.The characteristics of the sediments analyzed in this study (pH, electrical conductivity, and carbonate content), together with the percentage of organic matter and particle size, represent the important factors controlling the speciation and adsorption of trace metals in the environment [58,59].
Naturally (without any consideration of mining activities), the pH ranged from 8.12 to 9.3, with a mean of 8.61.The spatial variation in pH is limited, with a standard deviation of 0.42, indicating that the sediment is alkaline in nature.The alkaline pH of the sediments is mainly the result of the abundance of limestone bedrock and dolomite in the watershed.This is reflected in the high percentage of calcite (on average 32%), which acts as a counterbalancing factor against the acidification of the environment [60].The electrical conductivity (EC) depends on the dissolved ions that conduct the electrical current in the solution.In this study, the high conductivity values in the Moulouya River may be

Physiochemical Properties of Sediment
The physiochemical properties of the sediments (pH, electrical conductivity, and carbonate content) are shown in Table S5.The characteristics of the sediments analyzed in this study (pH, electrical conductivity, and carbonate content), together with the percentage of organic matter and particle size, represent the important factors controlling the speciation and adsorption of trace metals in the environment [58,59].
Naturally (without any consideration of mining activities), the pH ranged from 8.12 to 9.3, with a mean of 8.61.The spatial variation in pH is limited, with a standard deviation of 0.42, indicating that the sediment is alkaline in nature.The alkaline pH of the sediments is mainly the result of the abundance of limestone bedrock and dolomite in the watershed.This is reflected in the high percentage of calcite (on average 32%), which acts as a counterbalancing factor against the acidification of the environment [60].The electrical conductivity (EC) depends on the dissolved ions that conduct the electrical current in the solution.In this study, the high conductivity values in the Moulouya River may be attributed to the abundant carbonate substratum in the region.Alternatively, they could locally result from the input of untreated wastewater from urban agglomerations.

Trace Metal Spatial Distribution in Sediments
The trace element concentrations in the collected sediment samples varied significantly between the different locations.As shown in Figure 6, the ranges of variation in concentrations are as follows: 1.2-2.75mg/kg for Cd, 12-32.5 mg/kg for Cu, 25.38-215.21mg/kg for Pb, and 50.93-110.32mg/kg for Zn.In addition, the concentrations of As, Co, Cr and Ni remained below the ICP-AES detection limit set at 0.01 mg/L.The average concentrations of trace elements show a decreasing order as follows: Zn (83.30 mg/kg) > Pb (81.9 mg/kg) > Cu (20.1 mg/kg) > Cd (1.64 mg/kg), indicating that Pb and Zn have the highest concentrations.This could be due to the dominant extraction of cerussite, containing significant amounts of Pb and Zn.Consequently, these two trace metals were found in very high concentrations in the mine tailings, with values of 4822.90 and 112.7 mg/kg for Pb and Zn, respectively.concentrations are as follows: 1.2-2.75mg/kg for Cd, 12-32.5 mg/kg for Cu, 25.38-215.21mg/kg for Pb, and 50.93-110.32mg/kg for Zn.In addition, the concentrations of As, Co, Cr and Ni remained below the ICP-AES detection limit set at 0.01 mg/l.The average concentrations of trace elements show a decreasing order as follows: Zn (83.30 mg/kg) > Pb (81.9 mg/kg) > Cu (20.1 mg/kg) > Cd (1.64 mg/kg), indicating that Pb and Zn have the highest concentrations.This could be due to the dominant extraction of cerussite, containing significant amounts of Pb and Zn.Consequently, these two trace metals were found in very high concentrations in the mine tailings, with values of 4822.90 and 112.7 mg/kg for Pb and Zn, respectively.Furthermore, the analysis of the distribution of trace elements clearly reveals the contamination of the sediment quality due to the mine tailings.The highest concentrations were observed at sampling point MS3, immediately downstream of the tailings zone.More specifically, the concentrations of Pb and Zn reached more than 200 and 100 mg/kg, respectively.Additionally, a comparison of trace element concentrations reveals a significant increase in all elements (except for Cu) between the samples taken at point MS3 and the reference sample.For example, the concentrations determined for Cd, Pb, and Zn were 2.3, 7.7, and 1.4 times higher, respectively, than the reference sample.Considering that the region is devoid of any industrial activities likely to generate metal pollution (particularly lead pollution), these results indicate that the sediments in the vicinity of the mine tailings were significantly polluted by trace metals, transported from the tailings via erosion.
The river sediments were locally contaminated by the mine tailings.The concentrations of trace metals significantly decreased after point MS3 (except for MS7, which will be explained in Section 3.2.3).This observation is consistent with the sediment properties previously discussed, such as pH and CaCO3.These environmental characteristics increase the adsorption and limit the mobility of trace metals.In fact, the abundance of carbonates in the sediment acts as a buffer, enhancing the complexation of trace metals.This leads to the formation of secondary minerals such as carbonates and hydroxides, which adsorb the metals dissolved in water [61].In addition, the typical low flow of the Furthermore, the analysis of the distribution of trace elements clearly reveals the contamination of the sediment quality due to the mine tailings.The highest concentrations were observed at sampling point MS3, immediately downstream of the tailings zone.More specifically, the concentrations of Pb and Zn reached more than 200 and 100 mg/kg, respectively.Additionally, a comparison of trace element concentrations reveals a significant increase in all elements (except for Cu) between the samples taken at point MS3 and the reference sample.For example, the concentrations determined for Cd, Pb, and Zn were 2.3, 7.7, and 1.4 times higher, respectively, than the reference sample.Considering that the region is devoid of any industrial activities likely to generate metal pollution (particularly lead pollution), these results indicate that the sediments in the vicinity of the mine tailings were significantly polluted by trace metals, transported from the tailings via erosion.
The river sediments were locally contaminated by the mine tailings.The concentrations of trace metals significantly decreased after point MS3 (except for MS7, which will be explained in Section 3.2.3).This observation is consistent with the sediment properties previously discussed, such as pH and CaCO 3 .These environmental characteristics increase the adsorption and limit the mobility of trace metals.In fact, the abundance of carbonates in the sediment acts as a buffer, enhancing the complexation of trace metals.This leads to the formation of secondary minerals such as carbonates and hydroxides, which adsorb the metals dissolved in water [61].In addition, the typical low flow of the Moulouya River upstream of the watershed enhances the sedimentation of particles contaminated with trace metals and limits their transport downstream.

Assessment of the Sediment Pollution via Indexes
To assess the level of the metallic contamination of sediments, the Igeo index was applied (Figure 7).In the absence of the geochemical background, reference site values were used in this study.According to the Igeo index classification, the values range from −0.8 for Zn (uncontaminated) to 2.3 for Pb (moderately to heavily contaminated).At most locations, the metal indexes were either negative or nearly zero, indicating an absence of sediment contamination.However, for Pb, the contamination levels were higher, with a Igeo exceeding zero for six samples along the Moulouya.In addition, sampling points MS3 and MS7 showed the highest contamination indexes which was expected for MS3 given its proximity to the mine tailings.These results are similar to those obtained in different parts of the world, where the highest concentrations are determined in the vicinity of mine tailings [62,63].sediment contamination.However, for Pb, the contamination levels were higher, with a Igeo exceeding zero for six samples along the Moulouya.In addition, sampling points MS3 and MS7 showed the highest contamination indexes which was expected for MS3 given its proximity to the mine tailings.These results are similar to those obtained in different parts of the world, where the highest concentrations are determined in the vicinity of mine tailings [62,63].Furthermore, the watershed is devoid of industrial activities that could lead to the discharge of metal pollutants into the river, and the main activities in the area are agriculture (located upstream) and grazing.In addition, the mine tailings contain very high concentrations of trace metals (Table S6) and are subject to significant erosion, as shown in the previous results.Overall, these observations highlight the spread of contaminated particles from the tailings to the corresponding sediments.
In addition, we used the PLI, another contamination assessment index, which considers the multi-metallic contamination for each sediment sample.The PLI is calculated based on the contamination factors (CFs), providing a better understanding of the enrichment intensity.The results obtained using the CFs and PLI are shown in Table 4.The CF results reveal that most of the sediment samples show a low contamination of Cu, Zn, and Cd, except for two samples moderately contaminated with Cd.On the other hand, for Pb, the range of contamination factors extends from "low" to "very high", highlighting the contamination caused by the anthropogenic activities for this metal.The PLI ranges from 0.90 to 2.14, with an average of 1.24.The sediment samples located near mining tailings show the highest levels (MS3 = 2.14), revealing a significant trace metal contamination.This finding is consistent with our previous results, which show a high level of pollution around mining areas.Furthermore, the watershed is devoid of industrial activities that could lead to the discharge of metal pollutants into the river, and the main activities in the area are agriculture (located upstream) and grazing.In addition, the mine tailings contain very high concentrations of trace metals (Table S6) and are subject to significant erosion, as shown in the previous results.Overall, these observations highlight the spread of contaminated particles from the tailings to the corresponding sediments.
In addition, we used the PLI, another contamination assessment index, which considers the multi-metallic contamination for each sediment sample.The PLI is calculated based on the contamination factors (CFs), providing a better understanding of the enrichment intensity.The results obtained using the CFs and PLI are shown in Table 4.The CF results reveal that most of the sediment samples show a low contamination of Cu, Zn, and Cd, except for two samples moderately contaminated with Cd.On the other hand, for Pb, the range of contamination factors extends from "low" to "very high", highlighting the contamination caused by the anthropogenic activities for this metal.The PLI ranges from 0.90 to 2.14, with an average of 1.24.The sediment samples located near mining tailings show the highest levels (MS3 = 2.14), revealing a significant trace metal contamination.This finding is consistent with our previous results, which show a high level of pollution around mining areas.
However, despite the decrease in pollution as we move downstream from the mining areas, sampling point MS7, located 14 km from these wastes, shows significant pollution.The MS7 reveals a high concentration of metals, classifying it as the second most contaminated point in the river.It has a geo-accumulation index of 2 for Pb (moderately contaminated) and a metal load pollution index of 1.41.The abnormal pollution index values observed at this point can be attributed to two combined factors.Firstly, it is possible that this point is receiving pollution from contaminated soils located to the east of the tailings.As mentioned earlier, wind action carries particles from the unprotected tailings, contaminating the soil in that area.Additionally, the network of creeks drains this part of the basin and transports contaminated particles to this point in the river.Secondly, this area acts as a deposition and sedimentation zone for particles transported by the river during rainy periods.Indeed, the elevation profile (Figure 8) shows a very slight slope (0.16%) between points MS6 and MS7.This slope is significantly less than the gradient found between points MS3 and MS6.This difference in slope reduces the velocity of the flow, thereby increasing the deposition of particles released by the tailings in this location (MS7).Although the research in the Upper Moulouya suggests that contaminated particles will accumulate in the sediments around the tailings [18], it also highlighted the fact that flash floods can lead to the transport and deposition of sediments downstream of mining areas.Indeed, several studies [64,65] highlight the partial correlation between the migration of trace metals in sediments and the flow rate.This mechanism can lead to the resuspension of these particles, thereby affecting the spatial distribution of trace metals downstream of the river system.The impact of trace metal accumulation in sediments on the benthic fauna was assessed using SQGs.The results in Table S7 show that Cu and Zn levels in most samples did not exceed the Threshold Effect Concentration (TEC).This suggests that these two metals have a negligible impact on the benthic fauna.Furthermore, all the points show Cd concentrations ranging between the PEC and the TEC values.This indicates a potential risk of toxicity for bottom-dwelling organisms.On the other hand, Pb exceeded the PEC at 20% of the sampling points, surpassing the threshold above which negative effects are likely to occur.Furthermore, 40% of the sampling sites fall between the PEC and TEC values, indicating that Pb remains the most hazardous metal for the ecological system surrounding the tailings.A high level of toxicity indicated by the PEC can disrupt the ecological food chain by decreasing the abundance of primary organisms.These organisms are crucial suppliers of energy for the entire ecosystem [66].

Statistical Data Analysis
The results of the Shapiro-Wilk test [67], presented in Table S8 The impact of trace metal accumulation in sediments on the benthic fauna was assessed using SQGs.The results in Table S7 show that Cu and Zn levels in most samples did not exceed the Threshold Effect Concentration (TEC).This suggests that these two metals have a negligible impact on the benthic fauna.Furthermore, all the points show Cd concentrations ranging between the PEC and the TEC values.This indicates a potential risk of toxicity for bottom-dwelling organisms.On the other hand, Pb exceeded the PEC at 20% of the sampling points, surpassing the threshold above which negative effects are likely to occur.Furthermore, 40% of the sampling sites fall between the PEC and TEC values, indicating that Pb remains the most hazardous metal for the ecological system surrounding the tailings.A high level of toxicity indicated by the PEC can disrupt the ecological food chain by decreasing the abundance of primary organisms.These organisms are crucial suppliers of energy for the entire ecosystem [66].

Statistical Data Analysis
The results of the Shapiro-Wilk test [67], presented in Table S8, indicate that most of the data follow a normal distribution.These findings indicate that multivariate statistical methods could be used to analyze the relationships between trace metals and examine the similarities between their sources.In order to include the various factors that could influence the distribution of trace metals in sediments, a Pearson correlation analysis was carried out on trace metals, major elements, and physicochemical parameters (Table 5).The correlations (p < 0.05) between variables and factors are classified as highly significant for 0.7 < r < 1 and significant for 0.5 < r < 0.7 [43].The result of the study shows a positive correlation between the physicochemical parameters, pH and CaCO 3 (r = 0.660), while no significant correlation is observed between the pH/CaCO 3 % and the EC.Regarding the trace elements, significant correlations were observed between them: a very strong correlation between Cd and Cu (r = 0.901), as well as a strong correlation between Cd, Cu, and CaCO 3 .There was also a strong correlation between Pb and EC (r = 0.731).Zn showed strong correlations with Cd and Pb (r = 0.601 and r = 0.622, respectively), as well as a moderate correlation with Cu and CaCO 3 (r = 0.558 and r = 0.518, respectively).However, no correlation was found between Cu/Pb and Cd/Pb.Several studies have shown that a positive correlation between the variables may suggest that they follow the same transport mechanism, or come from the same source [68,69].In addition, the Pearson matrix reveals a strong positive correlation between the major elements (Fe, Mn, Na, K, Mg, Al, Ca).Finally, negative correlations are mainly observed between Cd, Cu, and most of the major elements.
To better understand the mechanisms underlying the trace metal enrichment in sediments and to identify the common sources of contamination, a Principal Component Analysis (PCA) was carried out (Figure 9).The PCA was performed on a dataset consisting of 14 variables and nine sampling points.The Kaiser-Meyer-Olkin (KMO) result of 0.612 indicates that the data were suitable for principal component analysis.According to the eigenvalue criterion, only factors with an eigenvalue greater than one have been retained in this study.The F1 and F2 factors were selected, with eigenvalues of 7.4 and 2.5, respectively.Together, they account for 71.2% of the variability in the data.
Analysis (PCA) was carried out (Figure 9).The PCA was performed on a dataset consisting of 14 variables and nine sampling points.The Kaiser-Meyer-Olkin (KMO) result of 0.612 indicates that the data were suitable for principal component analysis.According to the eigenvalue criterion, only factors with an eigenvalue greater than one have been retained in this study.The F1 and F2 factors were selected, with eigenvalues of 7.4 and 2.5, respectively.Together, they account for 71.2% of the variability in the data.The first factorial plane is mainly correlated with Cu, Cd, and CaCO 3 (correlation > 0.613) and shows a negative correlation with the major elements present in the sediments.The F2 is mainly influenced by Pb and electrical conductivity (EC), with values greater than 0.867.Zn, on the other hand, shows a moderate correlation with F1 (0.496) and a much higher correlation with F2 (0.643).
The correlation circle reveals that Cd and Cu involve a similar positive orientation, suggesting a potential common source for these metals.However, Pb does not show a positive correlation with these metals.This implies that mine tailings might not be the main source of metallic contamination for Cu and Cd.Indeed, it is possible that cadmium and copper originate from different sources unrelated to mining activities.Previous studies have indicated that cadmium and copper contamination may be linked to agricultural practices and the discharge of untreated wastewater [70,71].Moreover, the primary ore extracted from the Upper Moulouya is cerussite, which naturally contains significant amounts of Cu, Cd, and Zn [24].The contamination by Cu and Cd seems to be attributed to the rock alteration within the basin and anthropogenic sources of contamination.This suggests that Factor 1 is associated with non-mining-related metallic contamination.The F2 factor can be defined as the principal mining component, mainly due to the significant correlation with Pb.The high concentrations of lead in mine tailings, along with the contamination indices of sediments heavily affected by lead (Igeo, CF, and PLI), confirm that F2 is mainly correlated with the contamination from mine tailings.In addition, zinc showed a positive and moderate correlation with F1 and F2, which means that it is coming from multiple sources that are affected by both lithogenic and anthropogenic factors.
The positive correlation between CaCO 3 and trace metals highlights the major role played by CaCO 3 in the adsorption and the increase in trace metals' concentration in sediments.The positive correlation between EC and Pb on F2 suggests that their sources are common.This could be explained by the discharge of mine tailings loaded with minerals, which also leads to an increase in electrical conductivity.The major elements were negatively correlated with trace metals, indicating a competition pattern between them in the sediment adsorption sites.
The results of the PCA performed on the stations are shown in Figure 10.A zonation based on the physicochemical properties is observed, with four groups of stations identified.The first group, consisting of stations MS1, MS2, and MS4, is characterized by low metal pollution and represents non-contaminated sediments.The second group, composed of stations MS7, MS6, and MS5, has high concentrations of major elements but low concentrations of Cd and Cu.As a result, they are negatively positioned on the F1 axis, which reflects the Cd and Cu contamination.The third group is positively correlated with the F1 axis and includes two stations located on the Ansegmir (tributary of the Moulouya river).These stations are characterized by high Cd and Cu concentrations, which can be attributed to the intense agricultural activity prevalent in this region of the basin.Further studies carried out in the watershed have revealed that Cu and Cd may come from different sources than Pb and Zn [18,72].The results are consistent with the Midelt province report, which identifies the Ansegmir sub-catchment as one of the most extensive irrigated croplands in the watershed with an irrigated area of 1030 ha and a volume of 6 Mm/yr of water abstracted [53].These explanations support the suggestion that part of the Cd and Cu contamination in the watershed comes from agricultural activities or naturel sources.Finally, the last group consists of station MS3, located between both F1 and F2, indicating high contamination by all trace metals.Station MS3 is located downstream of the mine tailings, receiving large quantities of metals transported via erosion to the river.

Conclusions
At the first stage, the main objective of this research was to examine the transport flows of trace metals from a watershed containing mine tailings to the river, to determine their effect on the sediment quality.Our study also revealed that the basin has a low average erosion rate.This rate is mainly influenced by the topographical LS factor, which correlates strongly with soil erosion, as shown by the correlation coefficient of 0.94.In addition, the evaluation of the soil samples indicated that contamination is particularly severe in the vicinity of the mine tailings.Additionally, this contamination pattern aligns with the northeast and southwest directions, which are influenced by the prevailing winds.

Conclusions
At the first stage, the main objective of this research was to examine the transport flows of trace metals from a watershed containing mine tailings to the river, to determine their effect on the sediment quality.Our study also revealed that the basin has a low average erosion rate.This rate is mainly influenced by the topographical LS factor, which correlates strongly with soil erosion, as shown by the correlation coefficient of 0.94.In addition, the evaluation of the soil samples indicated that contamination is particularly severe in the vicinity of the mine tailings.Additionally, this contamination pattern aligns with the northeast and southwest directions, which are influenced by the prevailing winds.
Regarding the transport of metals to the river, it was observed that Pb exhibited an average transport rate of 166 g/ha/yr, with the highest concentrations detected 2 km from the mine tailings.Moreover, bare soil, identified as the primary land use impacting this flow, accounted for a significant contribution of 9803 kg/yr.The average concentrations of trace metals in the sediments were in an increasing order: Zn, Pb, Cu, then Cd.The levels of Pb and Zn exceeded concentrations of 200 mg/kg and 100 mg/kg, respectively.This contamination extended over more than 14 km from the mine tailings site.Contamination indices such as Igeo, PLI, and CF show that Pb is the most hazardous metal, particularly downstream of the tailings.Furthermore, statistical analysis shows that the sources of Cu and Cd pollution were different from those of Pb.This suggests that these metals are influenced by a variety of sources, both natural and anthropogenic, such as agriculture.This methodology demonstrates its ability to describe the complex transport of trace metals, providing useful explanations and a replicable approach for water resources management while addressing water quality surveys (which have potability purposes) and monitoring.Further analyses of sediments, water, and suspended solids are also planned, to identify the contamination caused by tailings in different parts of the ecosystem.These additional analyses will also be used to validate our model.

Figure 1 .
Figure 1.Position, topography, and land use distribution within the Zaïda watershed.

Figure 1 .
Figure 1.Position, topography, and land use distribution within the Zaïda watershed.

Figure 2 .
Figure 2. Maps illustrating the factors used to calculate soil erosion in the Zaïda catchment area: (a) the R factor, (b) the LS factor, (c) the K factor, and (d) the C factor.

Figure 2 .
Figure 2. Maps illustrating the factors used to calculate soil erosion in the Zaïda catchment area: (a) the R factor, (b) the LS factor, (c) the K factor, and (d) the C factor.

25 Figure 3 .
Figure 3. Statistical data on erosion based on land use in the Upper Moulouya Basin.

Figure 3 . 25 Figure 4 .
Figure 3. Statistical data on erosion based on land use in the Upper Moulouya Basin.Hydrology 2024, 11, x FOR PEER REVIEW 12 of 25

Figure 4 .
Figure 4. Spatial distribution of soil erosion in Zaida watershed.

Hydrology 2024 , 25 Figure 5 .
Figure 5. Spatial distribution of the average annual Pb load in the Zaida watershed.

Figure 5 .
Figure 5. Spatial distribution of the average annual Pb load in the Zaida watershed.

Figure 6 .
Figure 6.Spatial variation of Cd, Cu, Pb, and Zn contents in sediments of the Zaida watershed.

Figure 6 .
Figure 6.Spatial variation of Cd, Cu, Pb, and Zn contents in sediments of the Zaida watershed.

Figure 7 .
Figure 7. Geo-accumulation indexes for Cd, Cu, Pb, and Zn from various sediment samplings.

Figure 7 .
Figure 7. Geo-accumulation indexes for Cd, Cu, Pb, and Zn from various sediment samplings.

Figure 8 .
Figure 8. Altimetric profile of the Upper Moulouya River.

Figure 10 .
Figure 10.The biplot correlation of variables and stations on the F1-F2 plane.

Figure 10 .
Figure 10.The biplot correlation of variables and stations on the F1-F2 plane.

Table 1 .
Summary of the primary and secondary data used in this study.
4: International Soil Reference and Information Centre website.

Table 2 .
Pearson correlation matrix of the RUSLE factors and annual soil erosion.

Table 2 .
Pearson correlation matrix of the RUSLE factors and annual soil erosion.

Table 3 .
Statistical data on the concentration of four trace metals in the soil at the Zaida watershed.

Table 4 .
Contamination factor (CF) and pollution load index (PLI) of different trace metals in sediment samples from the Upper Moulouya.

Figure 8 .
, indicate that most of the data follow a normal distribution.These findings indicate that multivariate statistical Altimetric profile of the Upper Moulouya River.

Table 5 .
Pearson correlation matrix for physicochemical parameters, chemical elements, and trace metals.