Effect of Phosphogypsum on Soil Physical Properties in Moroccan Salt-Affected Soils

: Salinity is one of the most critical challenges for crop production and soil and water management in arid and semi-arid regions, such as a large area of Morocco. These regions are characterized by low rainfall with an erratic distribution, long drought periods, and high evaporation, resulting in salt accumulation in the superﬁcial layer of the soil and soil and water degradation. Therefore, phosphogypsum (PG) could be a promising amendment to reduce the salinity effect and improve soil quality in salt-affected soils. In this context, the present study aimed to evaluate the effect of PG on the physical properties of Luvisols and Cambisols collected from salt-affected soils in four regions in Morocco: Chichaoua, Ras El Ain, Sidi Zouine, and Sed El Masjoune. The treatments consisted of different rates of PG (15, 30, and 45 t/ha), natural Gypsum (G) (15 t/ha), and control. Our ﬁndings revealed that PG application improved soil structure by promoting ﬂocculant action provided by calcium. Linear regression indicated that Water Aggregate Stability (WAS) and PG doses were strongly correlated with a high coefﬁcient of determination (R2 = 93.41%, p value < 0.05). Compared to the control, the overall efﬁciency of 45 t/ha of PG amendment reached 53%, 95%, and 36%, respectively, in Chichaoua, Ras El Ain, and Sed El Masjoune soils. PG application presented a positive effect on other soil physical properties (soil hydraulic properties, total porosity, and bulk density), especially for the soils of Chichaoua and Ras El Ain regions. The total porosity was increased by 8% with 45 t PG/ha in Ras El Ain soil, and in Chichaoua soil, the bulk density was 5% lower in the pot treated with 45 t PG/ha compared to the control. This study supports the use of PG as an amendment for reclaiming salt-affected soils through monitoring agronomic and environmental impacts.


Introduction
Soil is one of the most important natural resources that may be prone to different forms of degradation. Soil degradation results in its quality and health deterioration and the reduction of its capacity to provide ecosystem services, and to ensure its main functions as well [1,2]. By 2050, population and income growth are expected to require a 70% increase in food production worldwide and up to 100% in developing countries [3], while the agricultural land surface is increasingly degraded. In fact, 33% of the land surface is subject to at least one soil degradation type [4]. The soil is susceptible to different types of degradation: physical, chemical, biological, and ecological. Degradation processes can be natural, but it has limited effects when compared to anthropogenic degradation by different factors, called accelerated degradation with large scale-impacts [5].
Salinity is a type of soil chemical degradation. It is one of the most limiting factors for crop production in arid and semi-arid regions, which are characterized by low and irregular annual rainfall, long drought periods, and high levels of evaporation, resulting in (1) Annual worldwide production of PG is about 300 Mt [28]. PG is recycled and reused in several sectors and applications [29], such as cement manufacturing, concrete, panels production [30,31], material for coastal protection [32], and, also, as daily landfill cover [32]. Diouri et al. reported that the use of phosphogypsum in the pavement structure showed promising results, good performance, and lower cost compared to other road construction materials [33]. PG has also been used in agriculture. It has shown positive effects on crop nutrition because it contains several nutrients (Ca, S, and P). Moreover, PG is used as an amendment for the reclamation of degraded soils: saline, sodic, acidic, and alkaline soils [26,34].
PG is an alternative to natural gypsum generally used to reduce salinity. Calcium sulfate can readily furnish soluble calcium to substitute exchangeable sodium [35]: CaSO 4 + 2Na-soil → Ca-Soil + Na 2 SO 4 (2) The PG showed its desalinization and desodification efficiencies by reducing electrical conductivity (EC) and exchangeable sodium percentage (ESP) [36]. Under saline conditions [37] found that P, Ca, and Fe soil contents were boosted by PG application when compared to control soil and gypsum-amended soil. In addition, PG improved rice and wheat yields. Miller reported that PG application improves soil's physical properties by increasing water infiltration and pore stabilization, and consequently, decreases soil erosion [38]. PG has also been approved for its ability to reduce soil dispersion [39] and enhance soil stabilization and aggregation [40]. The application of PG increases ion concentration, which leads to enhanced soil permeability [41]. Nayak et al. concluded that PG application improved hydraulic conductivity better than gypsum amendment [37]. On the other hand, PG contains many heavy metals. This aspect must be deeply investigated in future research.
Morocco is dominated by large arid and semi-arid conditions, and the salt-affected soil area is estimated to be 1.148 Mha [42]. This country is among the major PG producers in the world. However, it is essential to promote the use of PG as an amendment for the farmers by testing and evaluating the effect of PG in different soil classes and different saline-sodic soils [43]. Although several pieces of research were devoted to the agricultural uses of PG [44], few studies are interested in the evaluation of the use of PG as an amendment for reclaiming saline-sodic soils in Morocco [25]. In addition, the evaluations of soil physical properties were seldom considered in this type of study. The evaluation of PG as an amendment could contribute to understanding its mechanism of action and to define the optimal rate that could be used in different situations, to reduce the PG stocks and increase its valorization.
The objective of this study was to evaluate the impact of different phosphogypsum rates on soil physical properties of different salt-affected soils. We hypothesized that the application of phosphogypsum, as a source of calcium and sulfur besides its acidity that can contribute to soil lime solubilization could improve soil physical properties by reducing the dispersive effect of sodium and increasing cementing agents in the soil, which could improve soil aggregation, soil porosity, and soil water retention capacity.  Figure 1.

Materials and Methods
Sed el Masjoune, an area affected by salinity [24], has a semiarid climate with an annual average rainfall of 200 mm, a minimum temperature of 3.6 • C in winter, and a maximum temperature of 48 • C in summer [45]. The Chichaoua region is a part of a large watershed basin called Tensift. It is a semi-arid area with an average annual rainfall of 180 mm and an average temperature of 15 to 20 • C [46]. Sidi Zouine area is part of the Marrakech region; the minimum and maximum temperatures are, respectively, around 4.9 • C and 38 • C, with an average of 18 • C. The annual rainfall is about 190 mm [47]. Ras El Ain is located in the Rehamna region, which is characterized by a semi-arid climate with high precipitation interannual variability and an annual average rainfall of 189 mm. The temperature ranges between −3.4 and 35 • C [48] The soils were described and classified according to IUSS Working Group (2022) [49]. The soils were classified as Haplic Luvisol in Chichaoua and Ras El Ain regions, and as Haplic Cambisol in Sidi Zouine and Sed El Masjoune regions. According to different studies, the dominant clay mineral in these study regions is Montmorillonite [50,51].

Trials Setup
The experiment was conducted using pot trials under greenhouse conditions at Mohammed 6 Polytechnic University (UM6P) in Ben Guerir (Morocco). Each pot was filled with 10 kg of soil. Phosphogypsum (PG) was provided by the industrial complex of Jorf Lasfar, situated near El Jadida city. The second amendment consists of natural gypsum (G) for agricultural use. The experimental design consists of a completely randomized design (CRD) with three replicates. The treatments were: Control (0 t/ha), 15 t/ha of G, 15, 30, and 45 t/ha of PG. PG and G were mixed with the top 9 cm of the soil in the pot to ensure good homogenization.

Soil and Amendment Characterization
The soils were air-dried, crushed, and passed through a 2 mm sieve. pH was determined in 1:5 soil: water extract using a pH meter (InoLab pH 7310). The saturated paste was manually prepared and mixed [52], and electrical conductivity (ECe) was determined by an electrical conductivity meter (Mettler Toledo. SevenCompact). Determination of soil organic carbon was performed using sulfochromic oxidation of carbon in a mixture of potassium dichromate (K2Cr2O7) and sulfuric acid (H2SO4) at 135 °C according to Walkley and Black [53]. Calcium carbonate was analyzed using chlorohydric acid [54]. The available phosphorus was carried out by spectroscopy using the Olsen method [55]. Soil texture was determined using the hydrometer method proposed by Bouyoucos [56]. Spectrophotometry (Agilent Technologies. Cary 60 UV-Vis) was used to determine Sulfate content. The sodium, potassium, calcium, and magnesium contents were extracted with ammonium acetate at pH = 7 and determined by atomic absorption spectroscopy (Agilent Technologies. 200 Series AA). Soil total nitrogen was determined by distillation after mineralization using the Kjeldahl method [57]. The exchangeable sodium percentage and sodium adsorption ratio were calculated as indicated in equations 3 and 4: = Exchangeable {(Na)/(Ca + Mg + K + Na)} × 100 Na, Ca, Mg, and K expressed as mEq/100 g.

Trials Setup
The experiment was conducted using pot trials under greenhouse conditions at Mohammed 6 Polytechnic University (UM6P) in Ben Guerir (Morocco). Each pot was filled with 10 kg of soil. Phosphogypsum (PG) was provided by the industrial complex of Jorf Lasfar, situated near El Jadida city. The second amendment consists of natural gypsum (G) for agricultural use. The experimental design consists of a completely randomized design (CRD) with three replicates. The treatments were: Control (0 t/ha), 15 t/ha of G, 15, 30, and 45 t/ha of PG. PG and G were mixed with the top 9 cm of the soil in the pot to ensure good homogenization.

Soil and Amendment Characterization
The soils were air-dried, crushed, and passed through a 2 mm sieve. pH was determined in 1:5 soil: water extract using a pH meter (InoLab pH 7310). The saturated paste was manually prepared and mixed [52], and electrical conductivity (ECe) was determined by an electrical conductivity meter (Mettler Toledo. SevenCompact). Determination of soil organic carbon was performed using sulfochromic oxidation of carbon in a mixture of potassium dichromate (K 2 Cr 2 O 7 ) and sulfuric acid (H 2 SO 4 ) at 135 • C according to Walkley and Black [53]. Calcium carbonate was analyzed using chlorohydric acid [54]. The available phosphorus was carried out by spectroscopy using the Olsen method [55]. Soil texture was determined using the hydrometer method proposed by Bouyoucos [56]. Spectrophotometry (Agilent Technologies. Cary 60 UV-Vis) was used to determine Sulfate content. The sodium, potassium, calcium, and magnesium contents were extracted with ammonium acetate at pH = 7 and determined by atomic absorption spectroscopy (Agilent Technologies. 200 Series AA). Soil total nitrogen was determined by distillation after mineralization using the Kjeldahl method [57]. The exchangeable sodium percentage and sodium adsorption ratio were calculated as indicated in Equations (3) and (4): Na, Ca, Mg, and K expressed as mEq/100 g.
Na, Ca, and Mg expressed as mEq/l. Soils, PG, and G characterization was carried out at the Soil-Plant-Water Laboratory of the Agriculture Innovation Transfer Technology Center at the University Mohammed 6 Polytechnic at Ben Guerir (Morocco). Collected soils were classified for their salinity status following the classification proposed by Rengasamy [15]. The soil of Chichaoua is alkaline-saline soil (ECe > 4; SAR < 6; pH 8-9), while Ras Al Ain, Sidi Zouine, and Sed El Masjoune soils are very alkaline-saline-sodic soils (ECe > 4; SAR > 6; pH 8-9) ( Table 1). pH and EC were measured in the 1:5 PG or G: water extracts. The amendments were analyzed using aqua regia extraction followed by two hours of digestion at 95 • C. Ca, S, P, K, Mg, and Al were determined by ICP-OES (Agilent Technologies. 5110 ICP-OES) ( Table 2)

Soil Physical Evaluation
Soil physical attributes were evaluated after the Faba bean (Vicia faba L.) harvest. The sowing was on 22 December 2020, and the harvest was on 20 May 2021. The crop received Undisturbed soil samples were collected to determine the bulk density using rings (21.62 cm 3 ) according to the method proposed by Blake and Hartge [58]. These undisturbed samples were saturated by capillarity with water up to two-thirds of the ring height and submitted after saturation to the matric potentials of suction using 33 kPa and 1500 kPa to determine field capacity (FC) and permanent wilting point (PWP) [59], respectively. Available water capacity for plants was calculated using the Equation (5): Total porosity (TP) was estimated using saturated volumetric soil water content (θs). Aggregate stability was performed using the wet sieving apparatus (Eijkelkamp Agrisearch Equipment) [60]. Four grams of sieved aggregates from 1 to 2 mm air-dried soil was transferred to a sieve with 0.25-mm. Soil samples were pre-moistened for ten minutes before being submerged in distilled water. The apparatus sieved mechanically the samples in cans by an up and down movement for three minutes (42 cycles/min) while placing sufficient distilled water into the cans to cover the soil. Unstable aggregates, passed through the 0.25-mm sieve, were separated from the stable aggregates, and sand and organic matter were retained on the 0.25-mm sieve. The cans having unstable aggregates were oven-dried at 105 • C for 24 h and weighed. Retained materials in a 0.25-mm sieve were immersed and sieved again for ten minutes in cans filled with a dispersing solution (containing 2 g sodium hexametaphosphate/L). The water stable aggregate (WSA) was calculated as the weight of soil obtained in the dispersing solution cans divided by the sum of the weights obtained in the dispersing solution cans and distilled water cans: where: A − 0.2: is the oven-dried mass (g) of the water-stable aggregates (soil obtained in the dispersing solution cans) and B is the oven-dried mass (g) of water unstable aggregates. Overall efficiency (%) of the soil amendment was used to evaluate the effect of the rate of the amendment on the water stable aggregate, using Equation (7) [23].
Overall efficiency (%) = (WASt − WASc) WASc × 100 (7) WASt is the water aggregate stability in soil treated by amendment, and WASc is water aggregate stability in the control treatment.

Statistical Analysis
The different soil physical properties were used to evaluate the effects of phosphogypsum rates using analysis of variance (ANOVA). When ANOVA indicates a significant difference between treatments, post-hoc pair-wise comparisons were performed to investigate differences between treatment means using the Tukey test. All statistical analyses were performed using the R 3.4.1 software [61].

Water Aggregate Stability
The region and the amendment have a significant effect on the water aggregate stability (WAS) (Figure 2). Region effect differs from one soil to another, and the highest value of WAS was observed in the soil of Chichaoua region. Overall, PG application increases WAS, except for soil collected from the region of Sidi Zouine. The impact of PG on the improvement of soil structure was more remarkable in the soil of Chichaoua, and the overall efficiency of 45 t/ha of PG amendment reached 53% compared to the control. In soil from the Ras El Ain region, PG improved soil aggregation by 12% using 30 t PG/ha, and 95% using 45 t PG/ha. Concerning Sed El Masjoune soil, the PG improved the soil structure using a rate greater than 30t/ha, and the overall efficiency of PG was 36% and 64% for 45 t PG/ha and 30 t PG/ha, respectively (Table 3).
the soil structure using a rate greater than 30t/ha, and the overall efficiency of PG was 36% and 64% for 45 t PG/ha and 30 t PG/ha, respectively (Table 3).
There is a linear relationship between the PG application rate and WAS (Table 4 and Figure 3). This relationship in different soils showed two types of correlation. The first group presented a positive correlation in soil from the regions of Chichaoua, Ras El Ain, and Sad El Masjoune. The second group presented a negative correlation in soil from the region of Sidi Zouine (Table 4). Figure 3 represents the linear regression between the PG application rate and WAS for all regions. PG rates effect on Water Aggregate stability showed a high coefficient of determination (R 2 = 93.41%, p value <0.05).    There is a linear relationship between the PG application rate and WAS (Table 4 and Figure 3). This relationship in different soils showed two types of correlation. The first group presented a positive correlation in soil from the regions of Chichaoua, Ras El Ain, and Sad El Masjoune. The second group presented a negative correlation in soil from the region of Sidi Zouine. Figure 3 represents the linear regression between the PG application     Table 5 summarizes the phosphogypsum effect on soil hydraulic properties. Soils had significantly different field-capacity moisture (FC) (p value< 0.05). Values of FC in the soils from Chichaoua and Ras El Ain were significantly different from the others. For each soil, field capacity (FC), permanent wilting point (PWP), and available water capacity for plants (AWC) were generally lower in the control treatment compared to treated pots. However, the FC was not significantly affected by the applied PG rate. Regarding the AWC in the soil of Chichaoua, the effect of PG was significantly different from the treatment that received G and the amendment rate of 15 t/ha improved AWC by 16% compared to the control. In the same soil, treatment with 45 t/ha of PG presented unexcepted results with AWC lower than the control. Concerning Ras El Ain soil, the soil amendments (PG and G) used did not present a significant effect on the hydrological properties. However, the AWC was increased by 69% with 45 t/ha of PG compared to the control. Regarding the soils from the Sed El Masjoune and Sidi Zouine regions, there was no significant difference among the amendment rates.   Table 5 summarizes the phosphogypsum effect on soil hydraulic properties. Soils had significantly different field-capacity moisture (FC) (p value< 0.05). Values of FC in the soils from Chichaoua and Ras El Ain were significantly different from the others. For each soil, field capacity (FC), permanent wilting point (PWP), and available water capacity for plants (AWC) were generally lower in the control treatment compared to treated pots. However, the FC was not significantly affected by the applied PG rate. Regarding the AWC in the soil of Chichaoua, the effect of PG was significantly different from the treatment that received G and the amendment rate of 15 t/ha improved AWC by 16% compared to the control. In the same soil, treatment with 45 t/ha of PG presented unexcepted results with AWC lower than the control. Concerning Ras El Ain soil, the soil amendments (PG and G) used did not present a significant effect on the hydrological properties. However, the AWC was increased by 69% with 45 t/ha of PG compared to the control. Regarding the soils from the Sed El Masjoune and Sidi Zouine regions, there was no significant difference among the amendment rates. The values represent the mean (n = 3) ± the standard error. Lower-case letters compare the different treatments in same region and upper-case letters compare the same treatment in different region using the Tukey test (p < 0.05). Figure 4 shows the effects of phosphogypsum rates on the Total Porosity (TP). Soil region and amendment rate of PG presented a significant effect on TP (p value < 0.05). In general, the phosphogypsum increased TP. In Chichaoua Soil, TP was 6 and 8% greater in the 15 t PG/ha and 30 t PG/ha treatments than in the control treatment. In Ras El Ain soil, the increase in TP was proportional to the increase in phosphogypsum rates, and 45 t PG/ha resulted in an 8% increase in TP when compared to the control treatment. Regarding the Bulk density (BD), soil region and amendment rate of PG presented a significant effect on this soil property (p value < 0.05) ( Figure 5). A slight effect was observed of BD in Chichaoua Soil, and BD was 5% lower in the pot treated with 45 t PG/ha compared to the control. However, the soil of Sidi Zouine and Sed El Masjoune showed similar values of BD for the different phosphogypsum rates.

Porosity and Bulk Density
soil, the increase in TP was proportional to the increase in phosphogypsum rates, and 45 t PG/ha resulted in an 8% increase in TP when compared to the control treatment. Regarding the Bulk density (BD), soil region and amendment rate of PG presented a significant effect on this soil property (p value < 0.05) ( Figure 5). A slight effect was observed of BD in Chichaoua Soil, and BD was 5% lower in the pot treated with 45 t PG/ha compared to the control. However, the soil of Sidi Zouine and Sed El Masjoune showed similar values of BD for the different phosphogypsum rates.

Effect of Amendment by Phosphogypsum/Gypsum on Soil-Aggregate Stability
In general, soil structure is affected by soil salinization and sodification [62], which was supported by the results of this study, when comparing soil aggregation in different soils with different levels of salinization and sodification. The action of phosphogypsum on soil aggregation can occur by ensuring cation and electrolyte concentration balance   Figure 5. Effects of phosphogypsum rates on the bulk density. The values represent the mean (n = 3) and the vertical bars represent the standard error. Lower-case letters compare the different treatment in same region and upper-case letters compare the same treatment in different region using the Tukey test (p < 0.05).

Effect of Amendment by Phosphogypsum/Gypsum on Soil-Aggregate Stability
In general, soil structure is affected by soil salinization and sodification [62], which was supported by the results of this study, when comparing soil aggregation in different soils with different levels of salinization and sodification. The action of phosphogypsum on soil aggregation can occur by ensuring cation and electrolyte concentration balance and altering the soil susceptibility to dispersion, by affecting the balance between attractive and repulsive forces [38]. In addition, phosphogypsum application prevents soil structure deterioration during leaching and reduces sodium concentration and its dispersant effect on soil aggregation [63]. Phosphogypsum acted by bringing divalent cations (Calcium and Magnesium) in the soil solution to replace the monovalent cations (Sodium and Potassium) with Ca 2+ and Mg 2+ and contributes to the formation of aggregates and macro-pores. Therefore, improving soil structure and soil water retention.
According to Farahani et al. [64], the monovalent cation caused clay dispersion and reduced soil aggregation. In this study, we evaluated the effect of phosphogypsum rates on water aggregate stability (WAS). The results indicate that the application of PG explained 93% of WAS variation (coefficient of determination, R 2 = 0.93), as indicated in Figure 3 and Tables 3 and 4. This relationship was more pronounced in saline soils with a low level of salinity (SAR < 6), principally in the soil of Chichaoua and Ras El Ain (SAR = 9.9). Similar results were obtained by Farahani et al. and Levy et al. [64,65], who observed an effect of sodium on reducing aggregate stability and the amount of macro-aggregates. The response to amendment rate (gypsum or phosphogypsum) in different soils was different, as well as the reaction between sodium concentration and aggregate stability. This response can be explained by salinity and by the effect of other cement agents (soil organic matter and divalent cations) on soil aggregations. However, Armstrong and Tanton found that an application of 10 T/ha of gypsum did not have an effect on aggregate stability in saline-sodic clay soils [66].
Concerning the comparison between PG and G, chemical analysis revealed that PG has a lower pH ( Table 2). Ennaciri et al. attributed the acidity of the PG to the presence of residual acids like H 3 PO 4 , H 2 SO 4 , and HF [67]. PG has a higher calcium and sulfur content than natural gypsum. However, the latter is richer in Mg, K, and Al. Calcium and sulfur are the major constituents of PG and G. However, other elements present in low concentrations include magnesium, phosphorus, potassium, and aluminum [68,69]. Our finding confirmed that PG is more soluble than G (Table 2); as reported by Ennaciri et al. [67], due to the low pH of PG [70].
The ratio Ca: Mg and Ca: K are very low in the studied amendments. Indeed, K and Mg contents could not affect the amendment application. On the other hand, several authors reported that G and PG applications increased K and Mg leaching from the soil [71,72]. He et al. reported that Ca: Mg ratios did not influence clay dispersion of montmorillonite and kaolinite [73].

Effect of Amendment by Phosphogypsum/Gypsum on Soil Hydraulic Properties
The response to the application of phosphogypsum regarding hydraulic soil properties was observed in soil with a low-salinity level, principally in Chichaoua and Ras El Ain regions. According to Farahani et al., the presence of sodium in the soil increases clay dispersion, decreases soil aggregation, and reduces the macro-porosity in the soil, hence, reducing the soil water holding capacity [64]. Chawla et al. reported that soil water holding capacity is less for samples with higher salinity than those with lower soil salinity [74]. Overall, PG improved soil hydraulic properties compared to G and control, which could be explained by the effect of calcium and sulfate. PG promoted a greater concentration of calcium and sulfate compared to G ( Table 2). PG has a favorable effect in increasing Ca 2+ and reducing the pH and increasing the efficiency of reducing Na concentration in the soil. Furthermore, acidifying effect of PG contributed to the solubilization of soil calcium carbonates and enhanced Ca 2+ in the soil [71,75].
Melo et al. attributed the increase in soil water dynamic to the substitution of the sodium for the calcium contained in the gypsum amendment [76]. Soil hydraulic properties were considered by different studies to be invariant to variations of salt concentration in the soil [77]. Previous research by Aǧar showed that PG improved water infiltration over the 90 cm depth of a sodic soil and to leach 180 cm of water, and it took 1546 h in control without PG application, but only 368 and 225 h, respectively, for 7.5 and 15 t/ha of applied PG [78].
Cambisol and Luvisol, assessed in the presented study, are constituted by Montmorillonite. Montmorillonite is a 2:1 clay mineral characterized by a large specific surface area, high-cation-exchange capacity, and high ability to hold cations (Na + , Mg 2+ , K + ) [79]. In saline-sodic soils, the expandable interlayer space of Montmorillonite is occupied by sodium, causing soil swelling and clay dispersion. Sodium in the soil constituted a double layer due to low charge density [71], which reduces water flow in the soil, principally in soils with Montmorillonite [80]. However, the application of PG or G brings high content of calcium. Sodium is easily replaced by calcium due to valence, great hydrated ionic radius, and low adsorption selectivity (lyotropic series). Lyotropic series presented the relative strength of ion adsorption onto clay fraction. Lyotropic series for the principal cations in the soil is: Ca 2+ > Mg 2+ > K + ≈ NH4 + > Na + . Sodium ranked in the last position in this sequence, which means low electrical charges with clay fraction, which promotes its replacement [81], reducing the sodic effect in the soil ( Figure 6).

Effect of Amendment by Phosphogypsum/Gypsum on Total Porosity and Bulk Density
In the soil of Chichaoua, 45 t/ha of PG and 15 t/ha of G decreased Bulk density by 5% and 3%, respectively. A similar finding was reported by Müller et al., who indicated that the use of 12 t/ha of gypsum resulted in a 4% decrease in soil density compared to the control in Oxisols under no-till system [82]. The effect was clearer with the long-term application of gypsum (17 years) combined with green manure, limestone, and pasture. This strategy resulted in a 13.95% decrease in soil bulk density [83]. In this study, there is no effect of the amendment on the bulk density of the soils of Ras El Ain, Sidi Zouine, and Sed El Masjoune. This result can be explained by the short time of contact and reaction with the soil, the high-salinity level of those soils which probably require higher quantities of amendments. In fact, according to Gharaibeh et al., the positive effect of PG on salinesodic soils was correlated with rising PG rates and time of contact and reaction with the soil [84].
PG application improved the total porosity. This can be explained by the improvement of soil aggregation promoted by PG due to the flocculant action provided by calcium [76]. In line with our results, studies conducted by Nayak et al., comparing different soil

Effect of Amendment by Phosphogypsum/Gypsum on Total Porosity and Bulk Density
In the soil of Chichaoua, 45 t/ha of PG and 15 t/ha of G decreased Bulk density by 5% and 3%, respectively. A similar finding was reported by Müller et al., who indicated that the use of 12 t/ha of gypsum resulted in a 4% decrease in soil density compared to the control in Oxisols under no-till system [82]. The effect was clearer with the long-term application of gypsum (17 years) combined with green manure, limestone, and pasture. This strategy resulted in a 13.95% decrease in soil bulk density [83]. In this study, there is no effect of the amendment on the bulk density of the soils of Ras El Ain, Sidi Zouine, and Sed El Masjoune. This result can be explained by the short time of contact and reaction with the soil, the high-salinity level of those soils which probably require higher quantities of amendments. In fact, according to Gharaibeh et al., the positive effect of PG on saline-sodic soils was correlated with rising PG rates and time of contact and reaction with the soil [84].
PG application improved the total porosity. This can be explained by the improvement of soil aggregation promoted by PG due to the flocculant action provided by calcium [76]. In line with our results, studies conducted by Nayak et al., comparing different soil treated with gypsum, showed that gypsum-treated soils had higher porosity and soil aggregate stability [37]. Müller et al. demonstrated an increase in total porosity of 2.71% for 12 t/ha of gypsum and 2.38% for 4 t/ha [82].. This increase was mostly due to the increase in the total macropores.
The present study evaluated the effect of PG/G on the physical properties of Luvisols and Cambisols collected from salt-affected soils in four regions in Morocco. Overall, there is a clear difference between the soil of Chichaoua and the other soils, which could be explained by the low content of sodium in this soil (Soil no sodic: ESP = 7%). Overall efficiency is higher in Chichaoua soil by the fact that the calcium provided by the amendments (PG/G) was sufficient to replace the sodium in the clay fraction. Therefore, the reduction of the dispersive effect of sodium on clay and improvement of soil physical properties (aggregation, water retention, porosity, and bulk density).

Conclusions
This work presents the effect of phosphogypsum application as an amendment on some physical properties of Moroccan salt-affected soils. The overall efficiencies of phosphogypsum in different soils open the possibility of using phosphogypsum as an amendment in the reclamation process of marginated soils to improve soil structure, mitigate water stress, and reduce soil degradation processes and effects. Water stable aggregates improved with increasing phosphogypsum rate with a positive relationship and a high coefficient of determination (R 2 = 93.41%, p value < 0.05). Phosphogypsum application generally improved soil physical properties, principally for soils of Chichaoua and Ras El Ain regions. Phosphogypsum, a coproduct from the processing of phosphate rock, represents a promising solution for improving soil quality. The improvement of the soil's physical attributes enhances the root, vegetative and fructification development of the plants, which will positively affect agricultural yields, especially in arid and semiarid areas.