Journal of Hazardous Materials 480 (2024) 135732 Contents lists available at ScienceDirect Journal of Hazardous Materials journal homepage: www.elsevier.com/locate/jhazmat Unveiling the organic nature of phosphogypsum foam: Insights into formation dynamics, pollution load, and contribution to marine pollution in the Southern Mediterranean Sea Radhouan Belgacem El Zrelli a,* , Sébastien Fabre b, Sylvie Castet c, Michel Grégoire c , Oussema Fersi d, Claudie Josse e, Anne-Marie Cousin c , Pierre Courjault-Radé c a International Environmental Consultant, Toulouse, France Institut de Recherche en Astrophysique et Planétologie (IRAP), Université de Toulouse, 14 Avenue Edouard Belin, 31400 Toulouse, France c Géosciences Environnement Toulouse (GET), Université de Toulouse, UMR 5563 CNRS/UPS/IRD/CNES, 14 Avenue Edouard Belin, 31400 Toulouse, France d National Institute of Oceanography and Applied Geophysics (OGS), Borgo Grotta Gigante 42/C, 34010 Sgonico, TS, Italy e UAR Raimond CASTAING, Université de Toulouse, 31400 Toulouse, France b H I G H L I G H T S G R A P H I C A L A B S T R A C T • This is the first detailed organic characterization of PGF. • PGF is artificially induced, surfactantstabilized, and ephemeral aqueous foam. • Short-chain amphiphilic OM and gypsum crystals are pivotal in PGF formation. • PGF’s amphiphilic nature and specific OM groups cause its high contaminant load. • PGF accumulates, carries, and disperses PG radiochemical contaminants in the sea. A R T I C L E I N F O A B S T R A C T Keywords: Phosphogypsum foam Organic characterization Amphiphilic organic matter Surfactant Southern Mediterranean Sea The foamability of dissolved phosphogypsum from the phosphate fertilizer factories of Gabes (SE Tunisia) is a spectacular phenomenon that has not yet been thoroughly studied. The main objective of this research was to investigate the organic properties of phosphogypsum foam (PGF) to understand its formation process, determine the origin of its enhanced radiochemical contaminants load, and identify its role in pollutants dispersion in marine environment of the Southern Mediterranean Sea. This study identified PGF as an unnatural, surfactantstabilized, and ephemeral aqueous foam. PGF-forming process comprises three main steps: (i) formation (through phosphogypsum dissolution), (ii) stabilization (facilitated by organic surfactants and gypsum crystals), and (iii) destabilization (geochemical (involving the dissolution of the PGF skeleton gypsum) and/or mechanical (influenced by wind and wave action)). The amphiphilic nature of PGF organic matter and the presence of specific organic groups are responsible for its high toxic contaminants load. PGF contributes, through its elevated * Corresponding author. E-mail address:
[email protected](R.B. El Zrelli). https://doi.org/10.1016/j.jhazmat.2024.135732 Available online 5 September 2024 0304-3894/© 2024 Elsevier B.V. All rights are reserved, including those for text and data mining, AI training, and similar technologies. Journal of Hazardous Materials 480 (2024) 135732 R.B. El Zrelli et al. pollutants content and its ability to migrate far from its source, to the marine dispersion of industrial toxic radiochemical contaminants. It is therefore recommended to mitigate the environmental and health risks associated with PGF, including banning the discharge of untreated phosphogypsum and other industrial wastes into the coastal environment of Gabes. 1. Introduction small fraction recycled (̴ 15 %; [43]) in specific sectors (agricultural fertilizers, soil amendment for saline soils, road construction, building materials…; [46-50,43,51]). However, the majority of the produced PG (̴ 255 Mt.y−1) is stored in ‘wilderness’ stacks (e.g., Florida (USA; [52, 53]), Cubatão (Brazil; [54]), Sfax and Gafsa (Tunisia; [55,56]), Tébessa (Algeria; [57]), Nea Karvali (Greece; [58]), Qattinah (Syria; [59]), Huelva (Spain; [60,61]), Vasiliko (Cyprus; [62])) or directly dumped into aquatic environments without any treatment (sea or watercourses; e.g., Gabes (Tunisia; [63]), Jorf Lasfar and Safi (Morocco; [64]; Belahbib et al., [65]; [66]), Annaba (Algeria; [67]), Batroun (Lebanon; [68,69]), and São Paulo State (Brazil; [70])). Ecological, health, economic, and even social disasters/catastrophes are observed worldwide wherever the phosphate industry operates (e.g., Tunisia, Spain, Algeria, Lebanon, Greece, USA, Morocco, China…; [71,57,72,38,73–86,42,50,54,43,87]). Consequently, PG is not solely a problem of developing countries (Tunisia, Morocco, Algeria, Lebanon, Syria…) but also a global issue affecting even well-developed countries such as Spain, Canada, USA, Finland, Greece and China. Ultimately, relocating this polluting industry to third world ‘lawless’ countries or to countries with lax regulations will not resolve the environmental and health challenges associated with this industry in a world where pollution is a global concern. Since the 1970′s, the Central Gulf of Gabes acts as a repository for nearly all untreated solid and liquid industrial wastes of the Tunisian Chemical Group (Groupe Chimique Tunisien; GCT) factories (as well as other Gabesian chemical factories). This marine environment is characterized by a minimal agitation [88] with a moderate hydrodynamic regime [89]. These natural conditions significantly reduce the dissolution of gypsum, which represents the main component of phosphogypsum (more than 95 wt%; [90,41-43]). Knowing that PG constitutes the primary industrial waste released into the Gabesian littoral environment accumulating to over 500 million tons of untreated wet PG since 1972 (with daily quantities of ~ 30 ×103 tons; [39,91]). To facilitate the disposal of these substantial quantities into the Gabesian coastal environment, GCT dissolves the PG wastes in seawater, leading to the formation of the phosphogypsum foam ([92]; Fig. 1). The foamability of the dissolved PG, known as gypseous water, from the phosphate fertilizer factories in Gabes, is remarkable and visually striking, as depicted in Fig. 1. El Zrelli et al. [93] conducted comprehensive chemical and mineralogical characterizations of GCT’s PGF, nevertheless numerous aspects of this man-made foam remain unexplored and require further investigation. In this context, we employ chemical organic proxies for the first time to elucidate the PGF formation process, its elevated pollutants content, and its contribution to marine environmental pollution in the Central Gulf of Gabes (Southern Mediterranean Sea). Additionally, we discuss the health and the environmental impacts of PGF. For a long time, natural foams (e.g., foam nests of frogs [1,2], fish [3], insects [4,5], sea foams [6,7]), and man-made foams (e.g., soap foams [8,9], fire-fighting [10,11], beer [12,13] have intrigued the international scientific community. Consequently, numerous studies have focused on foams (e.g., [14-17]); each providing its own definition (Table 1). Generally, foam is defined as a colloidal dispersion [18] in the form of a gas cells enclosed in a liquid [19]. It forms at the surface of a liquid (water, oil, milk…) due to the dispersion of gas bubbles in a continuous liquid phase [18] or solid matrices [20]. Today, these foams find several important industrial applications (food processing, petrochemicals, fire protection, mineral flotation…; [21-24]) and constitute an entire science discipline known as ‘Aqueous Foam Science’, dating back to the 19th century, with its founding father being the blind Belgian physicist Joseph Plateau [25]. In spite of its great economic importance (agriculture, food industry, pharmaceutical industry…), the phosphate fertilizer industry is widely regarded as one of the most polluting industries globally [35-37]. This industrial activity is characterized by its diverse, hazardous, and substantial atmospheric emissions (e.g., Sulfur oxides (SO2), fluorinated gases, nitrogen oxides (NOx), ammonia gas (NH3), malodorous gases, dust…; [38]), liquid discharges (seawater used for the acid cooling system and gas washing…; EC, 2017; [39]), and solid waste outputs (phosphogypsum, sulfur dross, spent catalysts…; [40]). Phosphogypsum (PG) remains the primary mismanaged radiochemical waste by-product of the phosphate fertilizer industry globally [41-44]. Annually, approximately 300 Mt of PG are produced worldwide [45], with only a Table 1 Some definitions of foam. Foam definitions References Liquid foam is an example of soft matter (or a complex fluid) with a very well-defined structure. Agglomeration of gas bubbles separated from each other by thin liquid films. An impermanent form of matter in which a gas, often air, is dispersed in an agglomeration of bubbles that are separated from one another by films of liquid that is almost but not entirely water. Two-phase colloidal system in which the gas bubbles are dispersed in the continuous liquid phases. Temporary dilute dispersion of bubbles in the liquid, but on ageing, the structure gradually changes, and the bubbles transform into polyhedral gas cells with thin flat walls. Highly concentrated dispersions of gas (dispersed phase) in a liquid (continuous phase). Dispersion of gas bubbles in liquid or solid matrices. Dispersion of gas in a liquid or a solid, whereas the volume fraction of gas in the foam is mostly between 0.5 and 0.9. Structures of bubbles in contact. Liquid with gas cells enclosed. Dispersion in which a large proportion of gas by volume in the form of bubbles is dispersed in a liquid, solid or gel, hence forming closed cell structures. Dispersion of a high-volume fraction of gas bubbles is a small volume of liquid with some additives. Dispersion of gas in water. Bingham pseudoplastic fluids, which have yield stress, and the viscosity of foams decreases with increasing shear rate. Fluid mixtures of gas and liquid. Foam is a colloidal dispersion in the form of gas cells enclosed on a liquid. It’s formed on the surface of a liquid (water, oil, milk…) by a gas bubbles dispersion in a continuous liquid phase or solid matrices. [25] [14] [26] [27] [28] [17] [20] [29] 2. Materials and methods [30] [19] [15] Five samples of phosphogypsum foam, each weighing approximately 200 g, were collected from the littoral industrial discharge channel adjacent to the GCT factories in Gabes (Fig. 1). These samples were combined to create a composite sample of approximately 1 kg. The latter was then dried, homogenized, and stored according to the sampling protocol outlined by El Zrelli et al. [93]. In the laboratory, the PGF composite sample underwent successive extractions of its organic compounds. To ensure maximum extraction efficiency, we employed three organic solvents of varying polarities simultaneously: cyclohexane (Cyc, a nonpolar solvent), tetrahydrofuran (THF, a low-polarity solvent), and methanol (MeOH, a high polar [31] [32] [33] [34] This study 2 R.B. El Zrelli et al. Journal of Hazardous Materials 480 (2024) 135732 Fig. 1. Phosphogypsum foam layers of Tunisian Chemical Group in Chatt Sidi Abd Essalam Beach. solvent). For this purpose, 10 g of PGF was treated with 200 ml of each solvent (Cyc, THF, and MeOH). Then, the organic analyses of the different extraction products from the PGF sample (PGF-Cyc, PGF-THF, and PGF-MeOH) underwent separate analyses using five complementary techniques: Proton Nuclear Magnetic Resonance spectroscopy (1H NMR), High Performance Liquid Chromatography with Ultraviolet detection (HPLC/UV), Gas Chromatography-Mass Spectrometry (GCMS), Fourier Transform Infrared spectroscopy (FTIR), and Size Exclusion Chromatography (SEC) equipped with UltraViolet (UV) and Differential Refractometric (DR) detectors. (pH=2.65) with pure acetic acid, alongside a water/acetonitrile mixture (20:80 v/v) at the same pH (2.65). Ensuring the detection of the majority of compounds, the three PGF extracts were standardized to a concentration of 20 mg.ml−1 in acetonitrile. This technique is generally used to detect compounds with a molar mass of less than 3000 g.mol−1. 2.3. Gas chromatography-mass spectrometry The Gas chromatography-mass spectrometry was used to identify volatile and thermostable organic compounds. GC-MS analysis was performed using a Varian Saturn 2000 ion trap GC-MS system, coupled with a CP-3800 GC system equipped with a fused silica capillary DB-5MS column (30 ×0.25 mm, 5 % phenyl methylpolysiloxane, film thickness 0.25 µm). The analytical chromatographic conditions involved a temperature increase from 60 to 260 ◦ C at a gradient of 5 ◦ C.min−1, followed by 15 min of isothermal conditions at 260 ◦ C. A secondary gradient was applied, raising the temperature to 340 ◦ C at a rate of 40 ◦ C.min−1. The trap temperature was maintained at 250 ◦ C, while the transfer line temperature was set at 270 ◦ C. Mass scanning ranged from 40 to 650 m.z−1. Extracts were dissolved in their respective extraction solvents (Cyc, THF, and MeOH) at a concentration of 3 mg.ml−1, with 2 µl injected. We used two methods to identify the compounds in each extract. The first identification method involves comparing their retention index (RI) relative to C5-C24 n-alkanes obtained on a nonpolar DB-5MS column with literature values, while the second method entails comparing their mass spectra with those recorded in NIST/EPA/NIH Mass Spectral Library (NIST 08). 2.1. Proton nuclear magnetic resonance spectroscopy The Proton nuclear magnetic resonance spectroscopy technique involves the application of nuclear magnetic resonance principles, specifically targeting the hydrogen-1 nuclei within a substance’s molecules [94]. Its primary objective is to elucidate the molecular structure of the substance under investigation [94]. The 1H NMR analyses were carried out on a Bruker NMR-300 spectrometer operating at 300 MHz. This spectrometer uses the commonly used deuterated NMR solvents (CDCl3, THF-d8, and CD3OD) which allow to identify the most important molecules (polar or nonpolar) present in each PGF extract. A 5 mm diameter tube was filled with each sample, comprising 20 mg, along with 500 µl of deuterated solvent. 2.2. High performance liquid chromatography with ultraviolet detection The High performance liquid chromatography is a method using Ultraviolet detection (UV) to determine the chromatographic profile of UV–visible absorbing compounds. The HPLC analysis involved the utilization of a Thermo Scientific Dionex UltiMate 3000 Pump (Model HPG-3400RS), and an ultraviolet detector (UV-150 model detector) with a 280 nm detection wavelength. An RP-C18 column (Phenomenex) measuring 25 cm × 4.6 mm with 5 µm particle size, was used for separation at ambient temperature. Elution proceeded at a flow rate of 1.2 ml.min−1, utilizing a mobile phase composed of acidified water 2.4. Fourier transform infrared The Fourier transform infrared technique is used to determine the characteristic function of organic and inorganic compounds. The FT-IR analyses were obtained by using an Alpha-P spectrometer (Bruker). The measurements were performed on a Diamond Attenuated total Reflectance (ATR) Crystal, which was covered with a flow-through cell to facilitate sample analysis. The instrument was equipped with OPUS 3 R.B. El Zrelli et al. Journal of Hazardous Materials 480 (2024) 135732 Table 2 Extraction yields (% (w/w)) of cyclohexane (Cyc), tetrahydrofuran (THF), and methanol (MeOH) extracts from the phosphogypsum foam (PGF). PGF (%) Cyc THF MeOH 0.19 5.18 2.47 Table 3 Identification of 1H NMR chemical shifts functional group in cyclohexane (Cyc), tetrahydrofuran (THF), and methanol (MeOH) extracts from phosphogypsum foam (PGF). (Sm: small peak; Md: Medium peak; Gr: Great peak). Type of protons Chemical shift δ (ppm) Proton Formula Methyl CH3 0.8 – 2 Methyl, Methylene, Methine CH3, CH2, CH 2–5 Hydroxy, amine, thiol Alkene Aromatic Aldehydes Carboxylic acid OH, NH, SH 1.5 – 5 –C H-C– HCHO COOH 5.5 – 7 7–9 9 – 10 10 – 12 PGF Cyc THF MeOH Sm (1.7, 1.5, 1.1) Md (0.8) Gr (1.28) Sm (4.05, 2.3, 2, Md (0.9) Sm (1.1, 1.6, 1,8, 1,9, 1.95) Gr (1.32) Sm (1.82) Md (1.63) Sm (2.1, 2.3, 2.4, 3.32, 3.40, 3.42, 3.7, 3.8, 3.9, 4.3,5.1) Sm (2.18, 2.26,) Md (2.02) Gr (4.23, 3.96) Md (5.1) Fig. 2. (A) HPLC chromatograms of cyclohexane (Cyc), tetrahydrofuran (THF), and methanol (MeOH) extracts from the phosphogypsum foam (PGF). (B) Infrared spectrums of cyclohexane (Cyc), tetrahydrofuran (THF), and methanol (MeOH) extracts from the phosphogypsum foam (PGF). 3. Results Sm (5.4, 5.3) Sm (8.1) Sm (9.8) 3.1. Extraction yields software (version 7.5, Bruker). All PGF extracts spectra were recorded from 400 to 4000 cm−1 with a resolution exceeding 4 cm−1, co-adding 32 scans, achieving a frequency data accuracy better than 0.01 cm−1, and at room temperature. The extraction yields of organic solvents (Cyc, THF, and MeOH) from PGF samples are presented in Table 2. THF exhibited a significantly higher extraction yield (5.18 %) compared to both polar (MeOH) and nonpolar solvents (Cyc). The order of decreasing performance in extraction yield is as follows: THF˃MeOH˃Cyc (Table 2). Furthermore, PGF contains a higher proportion of polar compounds (2.47 %) compared to nonpolar compounds (0.19 %). 2.5. Size exclusion chromatography 3.2. Proton nuclear magnetic resonance spectroscopy analysis The Size exclusion chromatography was used for separating compounds based on their size, specifically their hydrodynamic volume, which was determined by how effectively they pass through the pores of the stationary phase [95]. The SEC system that was used included an HXL PROD01N pre-column mounted in series with three Waters Styragel columns (HR1 (separation range 100–5000), HR3 (separation range 500–30000), and HR4 (separation range 5000–600000)) with tetrahydrofuran as an eluent, at a flow rate of 1 ml.mn−1, and an injection volume of 100 µl. It was equipped with a differential refractometric detector (Optilab Rex Wyatt, T = 35 ◦ C) and a Varian UV detector (λ = 254 nm). To determine the average molar weights of sample compounds, mass calibration was done with increasing mass polymethyl methacrylate (PMMA) standards from Agilent (960, 1780, 4640, 10640, 26080, 66650, 265300, 675500 and 1568000 g.mol−1). All samples were dissolved in THF, filtered through 0.2 µm PTFE filters, and then analyzed in the SEC system under the same conditions as the PMMA standards. The molecular sample masses were calculated using the retention volume of compounds by ASTRA V Software (Wyatt Technology). For each extract, the weight average molecular weight (Mw), the number average Molecular weight (Mn), and the poly-dispersity index (Đ=Mw/Mn) were determined. Table 3 illustrates the identification outcomes of 1H NMR functional group chemical shifts in the three PGF extracts. The table reveals that traces of aromatics (8.1 ppm) and alkenes (5.3 and 5.4 ppm) are only detectable in the PGF-THF extract. Additionally, a broad peak at 5.1 ppm is solely present in the PGF-MeOH extract, indicating the likely presence of a hydroxy, amine, or thiol group. Nevertheless, for the methyl, methylene and methine functions (2–5 ppm), traces are observed in all three PGF extracts, with higher concentrations characterizing the PGF-MeOH extract (2.02, 3.96 and 4.23 ppm). While all three PGF extracts contain traces and medium quantities of methyl groups (0.8–2 ppm), major quantities are found only in the PGF-THF extract. 3.3. High performance liquid chromatography analysis The HPLC analysis of PGF-Cyc, THF, and MeOH extracts reveals detection by UV at 280 nm, as illustrated in Fig. 2.A. The HPLC spectra of these three extracts indicate the presence of few polar compounds (retention times: 1–10 min), and a larger quantity of nonpolar compounds (retention time: 40–60 min; Fig. 2). Furthermore, the chromatograms of the three PGF extracts exhibit distinct appearances between 40 and 60 min, with MeOH showing more absorbable compounds compared to the Cyc and THF. This discrepancy may stem from 4 Journal of Hazardous Materials 480 (2024) 135732 R.B. El Zrelli et al. Table 4 List of compounds identified in Gas Chromatography-Mass Spectrometry (GC-MS) analysis of cyclohexane (Cyc), tetrahydrofuran (THF), and methanol (MeOH) extracts from phosphogypsum foam (PGF). N◦ TR (min) Name Molecular formula Class/ Family 1 8.77 m-Cymene C10H14 Aromatic 2 9.14 3.3.5Trimethylcyclohexanone C9H16O Cyclic ketone Structure PGF Cyc THF MeOH £ £ 3 12.3 1.1′-Bicyclohexyl C12H22 Bicyclic saturated 4 13.68 3-Ethyl−5-methyl−1.2.4trithiolane C5H10S3 Sulfide £ 5 14.60 Hexathiane S6 Sulfide £ 6 15.72 2.6.10-Trimethyltetradecane C17H36 Alkane £ 7 15.95 Tetradecanoic acid C14H28O2 Fatty acid £ 8 16.84 Phthalic acid, butyl tetradecyl ester C26H42O4 Phtalate £ £ (continued on next page) 5 R.B. El Zrelli et al. Journal of Hazardous Materials 480 (2024) 135732 Table 4 (continued ) N◦ TR (min) Name Molecular formula Class/ Family 9 17.1 10 Structure PGF Heptadecanoic acid, 16methyl, methyl ester C19H38O2 Fatty acid ester 17.39 Hexadecanoic acid C16H30O4 Fatty acid £ 11 18.34 Cyclic-octaatomic sulfur S8 Sulfide £ 12 18.39 Methyl stearate C19H38O2 Fatty acid ester 13 22.96 Ethanediamide, N-(2ethoxyphenyl)-N′-(2ethylphenyl)- C18H20N2O3 Alcaloids Cyc THF MeOH £ £ £ £ differences in the chemical structures of their nonpolar compounds. albeit with different chemical functionalities: PGF-Cyc (1705 cm−1), PGF-THF (1702 cm−1), and PGF-MeOH (1698 cm−1). 3.4. Gas chromatography-mass spectrometry analysis 3.6. Size exclusion chromatography analysis Table 4 provides a summary of the compounds identified by GC-MS analysis of the three PGF extracts. Generally, the majority of identified compounds are nonpolar, with the presence of a few polar compounds. In total, 13 compounds were identified: 7 compounds for Cyc (1.1′bicyclohexyl; 2.6.10-trimethyltetradecane; tetradecanoic acid; phthalic acid, butyl tetradecyl ester; hexadecanoic acid; cyclic-octaatomic sulfur; and ethanediamide, N-(2-ethoxyphenyl)-N′-(2-ethylphenyl)-), 4 for THF (m-cymene; 3.3.5-trimethylcyclohexanone; 3-ethyl-5-methyl-1.2.4-trithiolane; and hexathiane), and 3 for MeOH (Heptadecanoic acid, 16methyl, methyl ester; hexadecanoic acid; and methyl stearate). However, the detection of hexadecanoic acid in both Cyc and MeOH GC-MS lists (Table 4) is attributed to its primary nonpolar function (majority), and its polar function (minority). Table 5 presents the calculated weight average molecular mass (Mw), number average molecular mass (Mn), and polydispersity index or dispersity (Đ= Mw/Mn) of the three studied PGF extracts. The PGF-THF extract contains compounds with higher Mw (3498 g.mol−1) and Mn (4694 g.mol−1) compared to those extracted by Cyc (777 and 230 g. mol−1, respectively) and MeOH (834 and 250 g.mol−1, respectively). The THF extract molecule size is around 20 times (Mn ratio) larger than the two other PGF extracts (Cyc and MeOH). Despite the high ratio of Mn, the PGF-THF molecules remain relatively small compared to polymers. The results obtained from all these five complementary techniques (1H NMR, HPLC/UV, GC-MS, FTIR, and SEC) provide an initial understanding of the nature of the PGF organic matter. In fact, it appears to be an amphiphilic (polar and nonpolar) organic matter composed of various organic families including aromatic, cyclic ketone, bicyclic saturated, sulfide, alkane, fatty acid, phthalate… Moreover, it is primarily formed by oligomers with molecular weights generally less than 5000 g.mol−1. 3.5. Fourier transform infrared spectroscopy analysis Fig. 2. B depicts the infrared spectra of Cyc, THF and MeOH extracts from the PGF. In the 3000 to 3500 cm−1 part of these spectra, the PGFMeOH extract exhibits the highest concentration of OH, NH, or SH, followed by PGF-THF, and finally the least concentrated, PGF-Cyc. Between 2800 and 3000 cm−1, the spectra of the three extracts are comparable, showing the presence of CH, CH2 and CH3. Lastly, in the 1650 to 1800 cm−1, carbonyl functions (C– –O) are observed in all samples, 6 R.B. El Zrelli et al. Journal of Hazardous Materials 480 (2024) 135732 Fig. 3. (A and B) Organic matter deposited on the surface of gypsum crystals in phosphogypsum (PG) wastes. (C) The beginning of phosphogypsum foam (PGF) formation in the form of the blackish-brown color organic suspension in gypseous water in the inside part of industrial channel. (D) The end of formation of PGF on the surface of gypseous water at the exit of the outside part of the industrial channel. 4. Discussion 4.1. Role of organic matter in the PGF formation process To prevent gypsum deposition at the entrance of the industrial littoral discharge channel and subsequent mouth clogging, PG is initially dissolved in seawater (utilized in sulfuric acid cooling process) forming gypseous water after its vacuum separation from phosphoric acid [93]. Then, the gypseous water is dumped into an artificial channel ⁓820 m long (within the GCT’s factories), where it is mixed with seawater used for sulfuric acid cooling process and gas washing (40 ×103 m3.h−1; [38]). Following this, the mixture is discharged into the sea via the extension of the channel outside the GCT’s factories (⁓665 m). In general, PG solubility increases with temperature and stirring time in the solvent medium [96]. Thereby, the dissolution of PG is facilitated by the strong artificial current within the industrial channel and the elevated temperature of the discharged seawater originated from sulfuric acid cooling process. The dissolution of gypsum releases the organic matter deposited on the surface of PG crystals (Fig. 3. A and B), and thus leads to the formation of an organic suspension in gypseous water which gives it this blackish-brown color (Fig. 3. C). Furthermore, the artificial agitation within the industrial channel generates air bubbles on the surface of the gypseous water solution (Fig. 3. C). As the solution moves through the channel towards its mouth, the PGF’s inter-bubbles liquid drains rapidly due the gravity difference between air and foam liquid [16], and thus a second less wet foam (Fig. 3. D) is progressively obtained from the first more wet foam (Fig. 3. C). Typically, during this foam formation stage (after the PGF drainage), and in the absence of stabilizing agents (e.g., proteins, organic matter, salts, solid particles…; [97-100]), the formed surface bubbles will burst when the thin liquid film separating them ruptures, causing to bubbles coalescence [16,101]. However, in our case, the amphiphilic nature of PGF organic matter allows it to adhere to the bubbles’ surface, rigidify them, accumulate them against each other, and thus forming a first wet foam on the surface of gypseous water (Fig. 3. C). Therefore, the PGF amphiphilic organic matter (surfactant) decreases the surface tension at Fig. 4. (A) Gypsum crystals skeleton at the bubble interfaces and fine gypsum crystals at the surface of phosphogypsum foam (PGF), in both backscattered (B) and secondary (C) electron image. bubbles’ interface [102], resulting in what we refer to as surfactant-stabilized foams [15,97]. The fine gypsum crystals within the PGF (Fig. 4. A) also synergistically contribute with surfactant organic molecules to the PGF stabilization, resembling other instances of foam stabilization by solid particles [103,98,104,100]. Fig. 4. B and C clearly illustrate that, in addition to the surfactant organic matter, the fine gypsum crystals on the PGF thin film surfaces contributes to consolidating the bubble interfaces, reduce the surface tension to avoid their rupture, and therefore the destruction of the PGF. At the end of this process, the foam gradually evolves from a less stable wet form to a much more stable less wet one (Fig. 3. D). However, as the PGF transitions from a gypsum saturated environment (industrial channel) to an undersaturated one (the open sea), the gypsum PGF skeleton crystals (Fig. 4. A) begin to dissolve, marking the onset of PGF geochemical destabilization. Moreover, waves and winds actively contribute to this destruction, hence we refer to it as PGF mechanical destabilization (Fig. 5). 7 R.B. El Zrelli et al. Journal of Hazardous Materials 480 (2024) 135732 Fig. 5. The various stages of geochemical and mechanical phosphogypsum foam (PGF) destruction during its marine migration: (A) Well- formed PGF on the seawater surface, (B) PGF degradation initiation, (C) Advanced PGF degradation, and (D) Complete PGF degradation. (Fig. 1). Kruglyakov et al. [104] elucidate the formation of such spectacular foams through the utilization of a mixture of nano/micro solid particles (small gypsum crystals in the PGF case; Fig. 4) and short-chain amphiphiles (oligomers in the PGF case; Table 5). The low density of the PGF’s organic matter (0.96 ± 0.01 g.cm−3) is less than the seawater average density (1.02–1.07 g.cm−3; [109]). The lightness allows PGF to float on the surface of seawater, to be carried long distances, and thus to disperse its elevated toxic pollutant load (trace elements, radionuclides, rare earth elements…) in the marine environment, after its degradation (geochemical and/or mechanical). Upon exiting the industrial waste discharge channel, the PGF can reach a thickness of ~8.5 to < 12 cm. Once in the marine environment, the PGF evolution depends primarily on the direction of the winds and, consequently, the amplitude of the waves. Specifically, if the winds blow (slowly) from the western sector, the PGF gradually migrates across the sea surface (owing to weak waves) over distances ranging from a few to several tens of kilometers, occasionally, reaching the Island of Djerba, as reported by numerous fishermen in Gabes and Zarzis (following a field survey carried out in 2013). During this marine migration, the PGF lightness allows its aerial part to function as a ‘sail’, facilitating forward movement aided by wind action. At this step, the two main antifoam agents (winds and waves) start to act, by breaking the thin films occurring between PGF bubbles, initiating mechanical foam destruction (destabilization). Occasionally, waves generated by fishermen’s boats can also contribute to this PGF mechanical destruction as they traverse the floating foam layers (Fig. 5. A). Throughout the PGF marine migration, the synergy between mechanical and geochemical foam destruction (Fig. 5. A to D) is the driving force behind the release of PGF’s radiochemical pollutant load into the marine environment of the Southern Mediterranean Sea. Thus, PGF plays the role of the main accumulating, transporting, and dispersing agent of PG radiochemical pollutants in the Central Gulf of Gabes [93]. Table 5 Calculated weight average molecular mass (Mw), number average molecular mass (Mn), and poly-dispersity index (Đ) of cyclohexane (Cyc), tetrahydrofuran (THF), and methanol (MeOH) extracts from the phosphogypsum foam (PGF). Samples PGF Cyc THF MeOH Mw (g.mol−1) Mn (g.mol−1) Đ 777 3482 834 230 4694 250 3.4 0.7 3.3 4.2. Role of organic matter in the PGF pollution load PGF stands out from other natural and man-made foams by its richness in organic matter (24.7 %; [92]). In addition to its main role in the PGF formation process, this organic matter also plays another role in relation with its strong pollutants load. In fact, previous studies have revealed that PGF contains high concentrations of trace elements (Zn, Cr, Cu, Cd, U, As…; [92]), radioelements (226Ra, 238U, 40K, 232Th…; [93]), and rare earth elements (Ce, Y, La, Nd, Pr, Sm…; [105]). The amphiphilic organic matter found in PGF (Table 2), is known for its strong affinity for toxic environmental pollutants [106-108]. Indeed, Chuang et al. [106] link this high affinity to the strong chelating binding sites associated to some of the hydrophilic (or charged) groups of amphipathic organic matter (e.g., R-N-R, R-S-R, -OH, -COOH, -O=POOO-, and -CONHO- groups). For instance, Pb prefers S (or N) over O binding sites, unlike Th, which prefers O over S and N binding sites [106]. In our study, we identified several of these sites in the three PGF extracts (e.g., -OH (tetradecanoic acid (in PGF-Cyc), hexadecanoic acid (in PGF-Cyc and -MeOH), and methyl stearate (in PGF-MeOH)), R-S-R (3-Ethyl-5-methyl-1.2.4-trithiolane (in PGF-THF), hexathiane (in PGF-THF), and cyclic-octaatomic sulfur (in PGF-Cyc), and R-N-R (Ethanediamide, N-(2-ethoxyphenyl)-N′-(2-ethylphenyl)- (in PGF-Cyc)); Table 4). Therefore, the amphiphilic nature of PGF organic matter and the presence of these organic groups most likely explain the potential PGF toxic contaminant load. 4.4. Health and environmental impacts of PGF Firstly, it’s important to state that the health and environmental impacts of PGF are closely related to its high contaminant load and the potential release of this toxic load, particularly after PGF geochemical and/or mechanical destruction. In Gabes, the coastal harbor developments play a key role in the spatial limitation of the environmental and health impacts of PGF, and other untreated industrial wastes discharged into the littoral 4.3. Role of organic matter, wind and wave in the diffusion of pollutants in the marine environment We have already noted the excessive foamability of the GCT’s PGF 8 R.B. El Zrelli et al. Journal of Hazardous Materials 480 (2024) 135732 Fig. 6. (A) Phosphogypsum foam in the beach of Casemate (July 2019). (B) Coastal artisanal fishing on the beach of Chatt Sidi Abd Essalam and in the middle of phosphogypsum foam waves (October 2022), with a crate of fish caught before (C) and after seawater washing (D). environment ([63,84] and 2019b). Confronted with the escalating concerns regarding industrial pollution in Gabes following the post-revolutionary events [110], authorities have opted to contain its impact geographically, instead of permanently halting coastal discharges of PG and other industrial wastes. Indeed, the piers of commercial harbor already ensure the protection of the northern beaches (Ghannouche, Ouedhref, Metouia…) against GCT marine pollution impacts. However, following the completion of the fishing harbor developments in 2018, the PGF seldom reaches the southern beaches of Gabes (Corniche, Casmate, M’tourrech, Tebelbou…; Fig. 6. A). Really, these works have transformed the inter-harbor zone (Chatt Sidi Abd Essalam Beach; Fig. 1) into a marine dumping ground for various industrial effluents (gypseous water, seawater used for the acid cooling process and gas washing (hydrogen fluoride (HF)) …; [39]), and solid wastes (phosphate grains, sulfur dross, spent catalysts…; [40]). As a result, PGF is often trapped in the inter-port area (Fig. 1). Consequently, at Chatt Sidi Abd Essalam Beach, local fishermen catch PGF-soaked fish normally inedible by humans, as depicted in Fig. 6. B, C and D. Through various and different pathways, people living in the vicinity of this coastal zone are chronically exposed to the toxic radiochemical pollutants of PGF (trace elements, radionuclides, rare earth elements…). These exposure routes include the dermal contact with seawater through fishing and bathing in the sea (Fig. 6. A and B), ingestion of contaminated fish (Fig. 6. C and D), inhalation of dried volatile PGF microfragments stranded on the beach, and direct exposure to radiation emitted by PGF ([93] and c). The occurrence of several diseases in the area, such as congenital malformations, cancers (e.g., lung, kidney, prostate, liver…), respiratory and cardiovascular diseases, infertility, bone fluorosis, and premature death ([40,41] and 2021a), may be closely related to such exposure. Nevertheless, it is worth to note that the PGF toxic pollutants act synergistically with the other untreated industrial toxic pollutants from GCT (e.g., sulfide oxides (SO2; 21.02 ×103 t. y−1), fluorinated gases (131.4 t.y−1), nitrogen oxides (NOx; 876 t.y−1), industrial dusts (1.6 ×106 t.y−1), sulfide dross, used catalysts…; [38, 40];). This cocktail of potentially toxic pollutants is likely exacerbating the health situation of Gabesian people [71,73,111,112]. However, it should be noted that only an epidemiological analysis could confirm with certainty the origin of the various illnesses suffered mainly by the people living near Chatt Sidi Abd Essalam beach. Similarly, marine animals inhabiting on this polluted coastal zone are also constantly exposed to the various toxic pollutants contained in the PGF via contact with polluted seawater, ingestion of contaminated organisms, inhalation of PGF volatile micro-fragments, and direct exposition to radiation emitted by PGF [92,93]. This is why, Chatt Sidi Abd Essalam Beach is where the majority of paralyzed (e.g., Larus michahellis, Sterna hirundo…), malformed (e.g., Liza aurata, Portunus segnis…), and dead marine animals (e.g., Caretta caretta (mass mortality), Phalacrocorax carbo (mass mortality), Liza saliens…; Fig. 7) are discovered. Locals have aptly dubbed this beach ‘the cemetery’ of the Mediterranean Sea. 9 R.B. El Zrelli et al. Journal of Hazardous Materials 480 (2024) 135732 Fig. 7. Photographs of some of paralyzed (A: Chroicocephalus ridibundus; B: Charadrius alexandrinus; C: Larus michahellis), malformed (D: Liza aurata; E: Portunus segnis; F: Belone belone), and dead animals (E, F, and G: Cartetta caretta); H: Sterna hirundo and its chick; I: Mola mola; J; Liza saliens) found in Chatt Sidi Abd Essalam Beach. 5. Conclusions may contribute to the dispersion of radiochemical pollutants of phosphate fertilizer factories in the Southern Mediterranean Sea. The dispersion of PGF toxic pollutants has significant impacts on human health and wildlife, especially in Chatt Sidi Abd Essalam littoral communities. One may note that other chemical and physical aspects pertaining to phosphogypsum foam formation warrant further investigation. These include its structure (bubbles shape, and arrangement…), rheology, gravity and capillary driven drainage, bubbles coarsening and coalescence… Additionally, the dynamic evolution of phosphogypsum foam from wet to dry forms involves myriad unknown processes occurring across several scales of length (molecular, mesoscopic, and In summary, phosphogypsum foam represents an artificial, surfactant-stabilized, and ephemeral aqueous foam. The formation process of this synthetic foam can be outlined in three steps: (i) formation (by dissolution of phosphogypsum wastes), (ii) stabilization (by organic surfactants and small gypsum crystals), and (iii) destabilization (geochemical (skeleton gypsum dissolution by seawater) and/or mechanical (resulting from the action of waves and winds)). The amphiphilic nature of PGF organic matter and its richness in specific organic groups (e.g., -OH, R-S-R, and R-N-R) are responsible of its high toxic pollutant load. Furthermore, the PGF low density and excessive foaming 10 Journal of Hazardous Materials 480 (2024) 135732 R.B. El Zrelli et al. macroscopic), and time (microseconds to hours). Further research in these areas is imperative for a comprehensive understanding of the behavior and implications of phosphogypsum foam in aquatic environments. [4] Szterk, A., Flis, S., Ofiara, K., Strus, B., 2024. Chemical composition of the foam enfolding juveniles of Aphrophora alni (Hemiptera: Aphrophoridae). J Asia-Pac Entomol 27 (1), 102185. https://doi.org/10.1016/j.aspen.2023.102185. [5] Tonelli, M., Gomes, G., Silva, W.D., Magri, N.T.C., Vieira, D.M., Aguiar, C.L., et al., 2018. Spittlebugs produce foam as a thermoregulatory adaptation. Sci Rep 8, 4729. https://doi.org/10.1038/s41598-018-23031-z. [6] Rahlff, J., Stolle, C., Giebel, H.-A., Mustaffa, N.I.H., Wurl, O., Herlemann, D.P.R., 2021. Sea foams are ephemeral hotspots for distinctive bacterial communities contrasting sea-surface microlayer and underlying surface water. FEMS Microbiol Ecol 97 (4), fiab035. https://doi.org/10.1093/femsec/fiab035. [7] Stogryn, A., 1971. The emissivity of sea foam at microwave frequencies. Antennas Propag Soc Int Symp. https://doi.org/10.1109/APS.1971.1150923. [8] Dixon, N., Morgan, M., Equils, O., 2017. Foam soap is not as effective as liquid soap in eliminating hand microbial flora. Am J Infect Control 45 (7), 813–814. https://doi.org/10.1016/j.ajic.2017.01.020. [9] Kraynik, A.M., Reinelt, D.A., 1996. Linear elastic behavior of dry soap foams. J Colloid Interface Sci 181 (2), 511–520. https://doi.org/10.1006/ jcis.1996.0408. [10] Gardiner, B.S., Dlugogorski, B.Z., Jameson, G.J., 1998. Rheology of fire-fighting foams. Fire Saf J 31 (1), 61–75. https://doi.org/10.1016/S0379-7112(97)000490. [11] Peshoria, S., Nandini, D., Tanwar, R.K., Narang, N., 2020. Short-chain and longchain fluorosurfactants in firefighting foam: a review. Environ Chem Lett 18, 1277–1300. https://doi.org/10.1007/s10311-020-01015-8. [12] Bamforth, C.W., 2004. The relative significance of physics and chemistry for beer foam excellence: theory and practice. J Inst Brew 110, 259–266. https://doi.org/ 10.1002/j.2050-0416.2004.tb00620.x. [13] Lyu, W., Bauer, T., Jahn, A., Gatternig, B., Delgado, A., Schellin, T.E., 2023. Experimental and numerical investigation of beer foam. Phys Fluids 35, 023318. https://doi.org/10.1063/5.0132657. [14] Bickerman, J.J., 1973. Foam films, in. Foams. Springer-Verlag, NY,, pp. 1–31. https://doi.org/10.1007/978-3-642-86734-7. [15] Hill, C., Eastoe, J., 2017. Foams: from nature to industry. Adv Colloid Interface Sci 247, 496–513. https://doi.org/10.1016/j.cis.2017.05.013. [16] Langevin, D., 2017. Aqueous foams and foam films stabilised by surfactants. Gravity-free studies. Comptes Rendus Mécanique 345 (1), 47–55. https://doi.org/ 10.1016/j.crme.2016.10.009. [17] Bhakta, A., Eli Ruckenstein, E., 1997. Drainage and Coalescence in Standing Foams. Journal of Colloid and Interface Science 191 (1), 184–201. https://doi. org/10.1006/jcis.1997.4953. [18] Schramm, L.L., 2005. Emulsions, Foams, and Suspensions: Fundamentals and Applications. Wiley-VCH Verlag GmbH & Co. KGaA,, Weinheim, Germany. 〈https://doi.org/10.1002/3527606750〉. [19] Fink, J.K., 2015. Dispersions, Emulsions, and Foams. Petroleum Engineer’s Guide to Oil Field Chemicals and Fluids, Chapter 21, 741–774. https://doi.org/ 10.1016/B978–0-12–383844-5.00021–0. [20] Weaire, D., Hutzler, S., 1999. The Physics of Foams. Clarendon Press,, Oxford. [21] Hargrave, J.M., Miles, N.J., Hall, S.T., 1996. The use of grey level measurement in predicting coal flotation performance. Miner Eng 9 (6), 667–674. https://doi.org/ 10.1016/0892-6875(96)00054-4. [22] Jäsberg, A., Selenius, P., Koponen, A., 2015. Experimental results on the flow rheology of fiber-laden aqueous foams. Colloids Surf A Physicochem Eng Asp 473, 147–155. https://doi.org/10.1016/j.colsurfa.2014.11.041. [23] Narchi, I., Vial, C.H., Djelveh, G., 2009. Effect of protein-polysaccharide mixtures on the continuous manufacturing of foamed food products. Food Hydrocoll 23 (1), 188–201. https://doi.org/10.1016/j.foodhyd.2007.12.010. [24] Rohani, M.R., Ghotbi, C., Badakhshan, A., 2014. Foam stability and foam-oil interactions. Pet Sci Technol 32 (15), 1843–1850. https://doi.org/10.1080/ 10916466.2012.683920. [25] Plateau, J., 1873. Experimental and theoretical statics of liquids subject to molecular forces only. Gauthier-Villars Paris,. [26] Aubert, J.H., Kraynik, A.M., Rand, P.B., 1986. Aqueous Foams. Scientific American 254 (5), 74–83. http://www.jstor.org/stable/24975954. [27] Walstra, P. ,1989. Principles of Foam Formation and Stability. In: Wilson, A.J., Ed., Foams: Physics, Chemistry and Structure, Springer, London, 1-51. https:// doi.org/10.1007/978-1-4471-3807-5_1. [28] Pugh, R.J., 1996. Foaming, foam films, antifoaming and defoaming. Advances in Colloid and Interface Science. 64, 67–142. https://doi.org/10.1016/0001-8686 (95)00280-4. [29] Arzhavitina, A., Steckel, H., 2010. Foams for pharmaceutical and cosmetic application. International Journal of Pharmaceutics 394 (1-2), 1–17. https://doi. org/10.1016/j.ijpharm.2010.04.028. [30] Karakashev, S.I., Grozdanova, M.V., 2012. Foams and antifoams. Advances in Colloid and Interface Science 176-177, 1–17. https://doi.org/10.1016/j. cis.2012.04.001. [31] Anazadehsayed, A., Rezaee, N., Naser, J., Nguyen, A.V., 2018. A review of aqueous foam in microscale. Adv Colloid Interface Sci 256, 203–229. https://doi. org/10.1016/j.cis.2018.04.004. [32] Yu, K., Li, B., Zhang, H., Wang, Z., Zhang, W., Wang, D., Xu, H., Harbottle, D., Wang, J., Pan, J., 2021. Critical role of nanocomposites at air–water interface: From aqueous foams to foam-based lightweight functional materials. Chemical Engineering Journal 416, 129121. https://doi.org/10.1016/j.cej.2021.129121. [33] Kouko, J., Prakash, B., Luukkainen, V.-M., Jäsberg, A., Koponen, A.I., 2021. Generation of aqueous foams and fiber foams in a stirred tank. Chemical Engineering Research and Design 167, 15–24. https://doi.org/10.1016/j. cherd.2020.12.013. Environmental implication The formation of PGF, which contains high levels of toxic trace, radioactive and rare earth elements, has not been investigated to date. Our study showed that the PGF formation process comprises three main stages: (i) formation, (ii) stabilization, and (iii) destabilization (geochemical and/or mechanical). The potential PGF toxic contaminants load threatens the marine environment and the coastal populations health. Indeed, it plays the role of the main accumulating, transporting, and dispersing agent of phosphogypsum radiochemical pollutants in the Gabes Gulf. This research is an important part of the response to the global pollution problem caused by the phosphate fertilizer industry. CRediT authorship contribution statement Pierre Courjault-Radé: Writing – review & editing, Supervision, Methodology, Conceptualization. Claudie Josse: Writing – review & editing, Formal analysis, Data curation. Anne-Marie Cousin: Writing – review & editing, Resources. Michel Grégoire: Writing – review & editing, Supervision, Funding acquisition, Conceptualization. Oussema Fersi: Writing – review & editing, Formal analysis, Data curation. Sébastien Fabre: Writing – review & editing, Formal analysis, Data curation. Sylvie Castet: Writing – review & editing, Formal analysis, Data curation. Radhouan Belgacem El Zrelli: Writing – original draft, Validation, Supervision, Resources, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization. Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Data availability Data will be made available on request. Acknowledgments In loving memory of my beloved father, Belgacem Ben Mohamed Ben Amor Zrelli (22.11.1946/05.05.2023), this work is dedicated. The authors extend heartfelt appreciation to all individuals whose support was invaluable in conducting fieldwork sampling and laboratory analyses. Special gratitude is expressed to the editor, Dr. Joao Pinto da Costa, and the three anonymous reviewers whose insightful comments and feedback greatly contributed to enhancing the quality of this manuscript. We are particularly grateful to Dr. Jessica K. Klar from the Interdisciplinary Laboratory of Numerical Sciences (LISN, Paris-Saclay University) for her editorial expertise in revising the manuscript. References [1] Cooper, A., Vance, S.J., Smith, B.O., Malcolm, W., Kennedy, M.W., 2017. Frog foams and natural protein surfactants. Colloids Surf A: Physicochem Eng Asp 534, 120–129. https://doi.org/10.1016/j.colsurfa.2017.01.049. [2] Shigeri, Y., Nakata, M., Kubota, H.Y., Tomari, N., Yamamoto, Y., Uegaki, K., et al., 2020. Identification of Novel Proteins in Foam Nests of the Japanese Forest Green Tree Frog, Rhacophorus arboreus. Zool Sci 38 (1), 8–19. https://doi.org/ 10.2108/zs200113. [3] Andrade, D.V., Abe, A.S., 1997. Foam nest production in the armoured catfish. J Fish Biol 50, 665–667. https://doi.org/10.1111/j.1095-8649.1997.tb01957.x. 11 R.B. El Zrelli et al. Journal of Hazardous Materials 480 (2024) 135732 [34] Yu, X., Li, F., Wang, J., Lin, Y., Zong, R., Lu, S., 2022. Effects of Fe (II) on stability of aqueous foam prepared by hydrolyzed rice protein in the presence of oil. Journal of Molecular Liquids 345, 117666. https://doi.org/10.1016/j. molliq.2021.117666. [35] Ahmad, N., Usman, M., Ahmad, H.R., Sabir, M., Farooqi, Z.U.R., Shehzad, M.T., 2023. Environmental implications of phosphate-based fertilizer industrial waste and its management practices. Environ Monit Assess 195, 1326. https://doi.org/ 10.1007/s10661-023-11958-4. [36] El Zrelli, R., 2017. Metallic trace element transfer modalities in the central part of Gabes Gulf, Tunisia: a geochemical, mineralogical, sedimentological and biological approach. (PhD dissertation). Paul Sabatier University, Toulouse, France. Available from: 〈https://theses.fr/2017TOU30232〉 (Accessed date: 07 August 2024). [37] Mishra, C.S.K., Nayak, S., Guru, B.C., Rath, B.C., 2010, Environmental impact and management of wastes from phosphate fertilizer plants. Journal of Industrial Pollution Control, 26(1), 57–60. Available from: 〈https://www.icontrolpollution. com/articles/environmental-impact-and-management-of-wastes-from-phosphat e-fertilizer-plants.pdf〉 (Accessed date: 07 August 2024). [38] EC, The European Commission, 2017a. Impact study of industrial pollution on the economy of the Gabès region, final report. Available from: 〈http://www.ods.nat. tn/upload/Rapport_Final.pdf〉 (Accessed date: 07 August 2024). [39] El Zrelli, R., Rabaoui, L., Ben Alaya, M., Daghbouj, N., Castet, S., Besson, P., et al., 2018. Seawater quality assessment and identification of pollution sources along the central coastal area of Gabes Gulf (SE Tunisia): evidence of industrial impact and implications for marine environment protection. Mar Pollut Bull 127, 445–452. https://doi.org/10.1016/j.marpolbul.2017.12.012. [40] El Zrelli, R., Rabaoui, L., Ben Alaya, M., Castet, S., Zouiten, C., Bejaoui, N., et al., 2019. Decadal effects of solid industrial wastes on the coast: Gulf of Gabes (Tunisia, Southern Mediterranean Sea) as an example. Estuarine Coastal Shelf Sci 224, 281–288. https://doi.org/10.1016/j.ecss.2019.04.021. [41] El Zrelli, R., Rabaoui, L., Daghbouj, N., Abda, H., Castet, S., Josse, C., et al., 2018. Characterization of phosphate rock and phosphogypsum from Gabes phosphate fertilizer factories (SE Tunisia): high mining potential and implications for environmental protection. Environ Sci Pollut Res 25, 14690–14702. https://doi. org/10.1007/s11356-018-1648-4. [42] Rutherford, P.M., Dudas, M.J., Samek, R.A., 1994. Environmental impacts of phosphogypsum. Sci Total Environ 149, 1–38. https://doi.org/10.1016/00489697(94)90002-7. [43] Silva, L.F.O., Oliveira, M.L.S., Crissien, T.J., Santosh, M., Bolivar, J., Shao, L., et al., 2022. A review on the environmental impact of phosphogypsum and potential health impacts through the release of nanoparticles. Chemosphere 286 (1), 131513. https://doi.org/10.1016/j.chemosphere.2021.131513. [44] Vásconez-Maza, M.D., Martínez-Segura, M.A., Bueso, M.C., Faz, Á., GarcíaNieto, M.C., Gabarrón, M., et al., 2019. Predicting spatial distribution of heavy metals in an abandoned phosphogypsum pond combining geochemistry, electrical resistivity tomography and statistical methods. J Hazard Mater 374, 392–400. https://doi.org/10.1016/j.jhazmat.2019.04.045. [45] Yang, L., Zhang, Y., Yan, Y., 2016. Utilization of Original Phosphogypsum as Raw Material for the Preparation of Self-Leveling Mortar. J Clean Prod 127, 204–213. https://doi.org/10.1016/j.jclepro.2016.04.054. [46] Cuadri, A.A., Navarro, F.J., García-Morales, M., Bolívar, J.P., 2014. Valorization of phosphogypsum waste as asphaltic bitumen modifier. J Hazard Mater 279, 11–16. https://doi.org/10.1016/j.jhazmat.2014.06.058. [47] Hentati, O., Abrantes, N., Caetano, A.L., Bouguerra, S., Gonçalves, F., Römbke, J., et al., 2015. Phosphogypsum as a soil fertilizer: Ecotoxicity of amended soil and elutriates to bacteria, invertebrates, algae and plants. J Hazard Mater 294, 80–89. https://doi.org/10.1016/j.jhazmat.2015.03.034. [48] Kuryatnyk, T., Angulski da Luz, C., Ambroise, J., Pera, J., 2008. Valorization of phosphogypsum as hydraulic binder. J Hazard Mater 160 (2-3), 681–687. https:// doi.org/10.1016/j.jhazmat.2008.03.014. [49] Kuzmanović, P., Todorović, N., Mrđa, D., Forkapić, S., Petrović, L.F., Miljević, B., et al., 2021. The possibility of the phosphogypsum use in the production of brick: radiological and structural characterization. J Hazard Mater 413, 125343. https://doi.org/10.1016/j.jhazmat.2021.125343. [50] Saadaoui, E., Ghazel, N., Ben Romdhane, C., Massoudi, N., 2017. Phosphogypsum: potential uses and problems – a review. Int J Environ Stud 74 (4), 558–567. https://doi.org/10.1080/00207233.2017.1330582. [51] Zrelli, A., Metoui, E., Doucouré, A., 2023. Studying the effect of phosphogypsum addition on ceramic membrane properties. Eng Technol J 41 (9), 1. https://doi. org/10.30684/etj.2023.139309.1426. [52] Burnett, W.C., Elzerman, A.W., 2001. Nuclide migration and the environmental radiochemistry of Florida phosphogypsum. J Environ Radioact 54 (1), 27–51. https://doi.org/10.1016/S0265-931X(00)00164-8. [53] Rutherford, P.M., Dudas, M.J., Arocena, J.M., 1995. Radioactivity and elemental composition of phosphogypsum produced from three phosphate rock sources. Waste Manag Res 13, 407–423. https://doi.org/10.1016/S0734-242X(05)800217. [54] Santos, A.J.G., Silva, P.S.C., Mazzilli, B.P., Fávaro, D.I.T., 2006. Radiological characterization of disposed phosphogypsum in Brazil: evaluation of the occupational exposure and environmental impact, Radiation. Prot Dosim 121 (2), 179–185. https://doi.org/10.1093/rpd/ncl011. [55] Jalali, J., Gaudin, P., Capiaux, H., Ammar, E., Lebeau, T., 2019. Fate and transport of metal trace elements from phosphogypsum piles in Tunisia and their impact on soil bacteria and wild plants. Ecotoxicol Environ Saf 174, 12–25. https://doi.org/10.1016/j.ecoenv.2019.02.051. [56] Zairi, M., Rouis, M.J., 1999. Impacts environnementaux du stockage du phosphogypse à Sfax (Tunisie). Bulletin des laboratoires des ponts et chaussées, 219 (4145), 29–40. Available from: 〈https://www.ifsttar.fr/collections/BLPCpdf s/blpc__219_29–40.pdf〉 (Accessed date: 07 August 2024). [57] Boumaza, B., Kechiched, R., Chekushina, T.V., Benabdeslam, N., Senouci, K., Hamitouche, A., et al., 2024. Geochemical distribution and environmental assessment of potentially toxic elements in farmland soils, sediments, and tailings from phosphate industrial area (NE Algeria). J Hazard Mater 465, 133110. https://doi.org/10.1016/j.jhazmat.2023.133110. [58] Noli, F., Sidirelli, M., Tsamos, P., 2024. The impact of phosphate fertilizer factory on the chemical and radiological pollution of the surrounding marine area (seawater and sediments) in northwestern Greece. Reg Stud Mar Sci, 103458. https://doi.org/10.1016/j.rsma.2024.103458. [59] Al Attar, L., Al-Oudat, M., Kanakri, S., Budeir, Y., Khalily, H., Al Hamwi, A., 2011. Radiological impacts of phosphogypsum. J Environ Manag 92 (9), 2151–2158. https://doi.org/10.1016/j.jenvman.2011.03.041. [60] Pérez-López, R., Castillo, J., Sarmiento, A.M., Nieto, J.M., 2011. Assessment of phosphogypsum impact on the salt-marshes of the Tinto River (SW Spain): role of natural attenuation processes. Mar Pollut Bull 62, 2787–2796. https://doi.org/ 10.1016/j.marpolbul.2011.09.008. [61] Rentería-Villalobos, M., Vioque, I., Mantero, J., Manjón, G., 2010. Radiological, chemical and morphological characterizations of phosphate rock and phosphogypsum from phosphoric acid factories in SW Spain. J Hazard Mater 181 (1-3), 193–203. https://doi.org/10.1016/j.jhazmat.2010.04.116. [62] Liatsou, I., Pashalidis, P., 2016. Radio-environmental impacts and uranium radiochemistry of phosphogypsum disposed at a coastal area in Cyprus, 4th International Conference on Sustainable Solid Waste Management, Limassol, Cyprus. Available from: 〈https://api.semanticscholar.org/CorpusID:23463978〉 (Accessed date: 07 August 2024). [63] El Zrelli, R., Courjault-Radé, P., Rabaoui, L., Castet, S., Michel, S., Bejaoui, N., 2015. Heavy metal contamination and ecological risk assessment in the surface sediments of the coastal area surrounding the industrial complex of Gabes city Gulf of Gabes, SE Tunisia. Mar Pollut Bull 101, 922–929. https://doi.org/ 10.1016/j.marpolbul.2015.10.047. [64] Cheggour, M., Langston, W.J., Chafik, A., Texier, H., Idrissi, H., Boumezzough, A., 1999. Phosphate industry discharges and their impact on metal contamination and intertidal macrobenthos: Jorf Lasfar and Safi coastline (Morocco). Toxicol Environ Chem 70, 159–179. https://doi.org/10.1080/02772249909358747. [65] Belahbib, L., Arhouni, F.E., Boukhair, A., Essadaoui, A., Ouakkas, S., Hakkar, M., et al., 2021. Impact of phosphate industry on natural radioactivity in sediment, seawater, and coastal marine fauna of El Jadida Province, Morocco. J Hazard Toxic Radioact Waste 25 (1), 04020010. https://doi.org/10.1061/(ASCE) HZ.2153-5515.0000563. [66] Akfas, F., Elghali, A., Aboulaich, A., Munoz, M., Benzaazoua, M., Bodinier, J.-L., 2024. Exploring the potential reuse of phosphogypsum: a waste or a resource? Sci Total Environ 908, 168196. https://doi.org/10.1016/j.scitotenv.2023.168196. [67] Chaalal, O., Madhuranthakam, C.M.R., Moussa, B., Hossain, M.M., 2020. Sustainable approach for recovery of sulfur from phophogypsum. ACS Omega 5 (14), 8151–8157. https://doi.org/10.1021/acsomega.0c00420. [68] El Samad, O., Aoun, M., Nsouli, B., Khalaf, G., Hamze, M., 2014. Investigation of the radiological impact on the coastal environment surrounding a fertilizer plant. J Environ Radioact 133, 69–74. https://doi.org/10.1016/j.jenvrad.2013.05.009. [69] Fakhri, M., Abboud - Abi Saab, M., Romano, J-C., 2008. The use of sediments to assess the impact of Selaata phosphate plant on Batroun coastal area (Lebanon, Levantine Basin). Lebanese Science Journal 9 (1), 29–42. Available from: 〈https ://lsj.cnrs.edu.lb/wp-content/uploads/2015/12/fakhri1.pdf〉 (Accessed date: 09 August 2024). [70] Sanders, L.M., Luiz-Silva, W., Machado, W., Sanders, C.J., Marotta, H., EnrichPrast, A., et al., 2013. Rare earth element and radionuclide distribution in surface sediments along an estuarine system affected by fertilizer industry contamination. Water Air Soil Pollut 224, 1742. https://doi.org/10.1007/s11270-013-1742-7. [71] Abdenneji, S., Fehri, S.M., Kwass, H., 2024. Pollution dans la région de Gabès: corrélation avec l′évolution des pneumonies communautaires aiguës. Rev Des Mal Respir Actual 16 (1), 211–212. https://doi.org/10.1016/j.rmra.2023.11.435. [72] Darmoul, B., 1988. Pollution dans le Golfe de Gabès (Tunisie): bilan des six années de surveillance (1976–1981). Bull. Inst. Natl. Sci. Technol. Mer Salammbô 15, 61–85. Available from: 〈https://n2t.net/ark:/68747/INSTM.Bulletin.v15. 1010〉 (Accessed date: 07 August 2024). [73] EC, The European Commission, 2017b. Impact study of industrial pollution on human health in Gabès, final report. [74] El Kateb, A., Stalder, C., Rüggeberg, A., Neururer, C., Spangenberg, J.E., Spezzaferri, S., 2018. Impact of industrial phosphate waste discharge on the marine environment in the Gulf of Gabes (Tunisia). PLOS ONE 13 (5), e0197731. https://doi.org/10.1371/journal.pone.0197731. [75] El Zrelli, R.B., Yacoubi, L., Castet, S., Grégoire, M., Lin, Y.-J., Attia, F., et al., 2023. Compartmentation of trace metals in Cymodocea nodosa from a heavily polluted area (Central Gulf of Gabes; Southern Mediterranean Sea): potential use of the seagrass as environmental monitoring and bioremediation tool. Reg Stud Mar Sci 65, 103056. https://doi.org/10.1016/j.rsma.2023.103056. [76] El Zrelli, R., Yacoubi, L., Castet, S., Grégoire, M., Josse, C., Olive, J.-F., et al., 2023. PET plastics as a Trojan horse for radionuclides. J Hazard Mater 441, 129886. https://doi.org/10.1016/j.jhazmat.2022.129886. [77] El Zrelli, R., Rabaoui, L., Roa-Ureta, R.H., Gallai, N., Castet, S., Grégoire, M., et al., 2020. Economic impact of human-induced shrinkage of Posidonia oceanica meadows on coastal fisheries in the Gabes Gulf (Tunisia, Southern Mediterranean 12 R.B. El Zrelli et al. [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] Journal of Hazardous Materials 480 (2024) 135732 Sea). Mar Pollut Bull 155, 111124. https://doi.org/10.1016/j. marpolbul.2020.111124. El Zrelli, R., Courjault-Radé, P., Rabaoui, L., Daghbouj, N., Mansour, L., Balti, R., et al., 2017. Biomonitoring of coastal pollution in the Gulf of Gabes (SE. Tunisia): use of Posidonia oceanica seagrass as a bioindicator and its mat as an archive of coastal metallic contamination. Environ Sci Pollut Res 24, 22214–22225. https:// doi.org/10.1007/s11356-017-9856-x. Hammouda, A., Ayadi, T., Selmi, S., 2024. Long-term exposure to industrial chemical contamination affects the magnitude of predator-induced immunosuppression in a free-living passerine. Bull Environ Contam Toxicol 112, 42. https://doi.org/10.1007/s00128-024-03857-2. Hattab, S., Boughattas, I., Cappello, T., Zitouni, N., Touil, G., Romdhani, I., et al., 2023. Heavy metal accumulation, biochemical and transcriptomic biomarkers in earthworms Eisenia andrei exposed to industrially contaminated soils from southeastern Tunisia (Gabes Governorate). Sci Total Environ 887, 163950. https://doi. org/10.1016/j.scitotenv.2023.163950. May, A., Sweeney, J.W., 1982. Assessment of Environmental Impacts Associated with Phosphogypsum in Florida. United States Department of the Interior. The Fertilizer Institute, Bureau of Mines Report of Investigations, RI: 8639. 〈https:// doi.org/10.1520/STP30276S〉. Mouawad, R., Khalaf, G., Salameh, Y., 2009. Impact of phosphogypsum and other factory effluents on meiofauna communities of Batroun coastal region. Lebanese Science Journal, 10(1), 23–34. Available from: 〈https://lsj.cnrs.edu.lb/wp-conten t/uploads/2015/12/mouawad1.pdf〉 (Accessed date: 07 August 2024). Pérez-López, R., Álvarez-Valero, A.M., Nieto, J.M., 2007. Changes in mobility of toxic elements during the production of phosphoric acid in the fertilizer industry of Huelva (SW Spain) and environmental impact of phosphogypsum wastes. J Hazard Mater 148, 745–750. https://doi.org/10.1016/j.jhazmat.2007.06.068. Rabaoui, L., Balti, R., El Zrelli, R., Tlig-Zouari, S., 2014. Assessment of heavy metals pollution in the gulf of Gabes (Tunisia) using four mollusk species. Mediterr Mar Sci 15, 45–58. https://doi.org/10.12681/mms.504. Rabaoui, L., El Zrelli, R., Ben Mansour, M., Balti, R., Mansour, L., Tlig-Zouari, S., et al., 2015. On the relationship between the diversity and structure of benthic macroinvertebrate communities and sediment enrichment with heavy metals in Gulf of Gabes Tunisia. J Mar Biol Assoc UK 95, 233–245. https://doi.org/ 10.1017/S0025315414001489. Rabaoui, L., El Zrelli, R., Balti, R., Mansour, L., Courjault-Radé, P., Daghbouj, N., et al., 2017. Metal bioaccumulation in two edible cephalopods in the Gulf of Gabes, South-eastern Tunisia: environmental and human health risk assessment. Environ Sci Pollut Res 24, 1686–1699. https://doi.org/10.1007/s11356-0167945-x. Wang, Z., Ma, X., Pan, H., Yang, X., Zhang, X., Lyu, Y., et al., 2023. Investigating effects of phosphogypsum disposal practices on the environmental performance of phosphate fertilizer production using emergy analysis and carbon emission amounting: a case study from China. J Clean Prod 409, 137248. https://doi.org/ 10.1016/j.jclepro.2023.137248. Darmoul, B., Hadj Ali Salem, M., Vitiello, P., 1980. Effets des rejets industriels de la région de Gabès (Tunisie) sur le milieu récepteur. Bulletin de l′Institut National des Sciences et Technologie de la Mer de Salammbô, 5–61. Available from: 〈htt ps://n2t.net/ark:/68747/INSTM.Bulletin.v7.1138〉 (Accessed date: 07 August 2024). Poizat, C., 1970. Hydrodynamisme et sédimentation dans le Golfe de Gabès (Tunisie). Téthys 2 (1), 267–296. Available from: 〈http://paleopolis.rediris.es/ benthos/REF/som/T-pdf/1970_2–1-267.pdf〉 (Accessed date: 07 August 2024). Bilal, E., Bellefqih, H., Bourgier, V., Mazouz, H., Dumitraş, D.-G., Bard, F., et al., 2023. Phosphogypsum circular economy considerations: a critical review from more than 65 storage sites worldwide. J Clean Prod 414, 137561. https://doi. org/10.1016/j.jclepro.2023.137561. El Zrelli, R., Hcine, A., Yacoubi, L., Roa-Ureta, R.H., Gallai, N., Castet, S., et al., 2023. Economic losses related to the reduction of Posidonia ecosystem services in the Gulf of Gabes (Southern Mediterranean Sea). Mar Pollut Bull 186, 114418. 〈https://doi.org/10.1016/j.marpolbul.2022.114418〉. El Zrelli, R., Rabaoui, L., Van Beek, P., Castet, S., Souhaut, M., Grégoire, M., et al., 2019. Natural radioactivity and radiation hazard assessment of industrial wastes from the coastal phosphate treatment plants of Gabes (Tunisia, Southern Mediterranean Sea). Mar Pollut Bull 146, 454–461. https://doi.org/10.1016/j. marpolbul.2019.06.075. El Zrelli, R., Rabaoui, L., Abda, H., Daghbouj, N., Pérez-López, R., Castet, S., et al., 2019. Characterization of the role of phosphogypsum foam in the transport of [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] 13 metals and radionuclides in the Southern Mediterranean Sea. J Hazard Mater 363, 258–267. https://doi.org/10.1016/j.jhazmat.2018.09.083. Nikolić, V., Ilić-Stojanović, S., Petrović, S., Tačić, A., Nikolić, L., 2019. Chapter 21 - Administration Routes for Nano Drugs and Characterization of Nano Drug Loading, In Micro and Nano Technologies. Charact Biol Nanomater Drug Deliv 587–625. https://doi.org/10.1016/B978-0-12-814031-4.00021-0. Nagy, K., Vékey, K., 2008. Chapter 5 - Separation methods. Med Appl Mass Spectrom 61–92. https://doi.org/10.1016/B978-044451980-1.50007-0. Bejaoui, B., Rais, S., Koutitonsky, V., 2004. Modélisation de la dispersion du phosphogypse dans le golfe de Gabès. Bulletin de l′Institut National des Sciences et Technologie de la Mer de Salammbô 31,103–109. Available from: 〈https://n2t. net/ark:/68747/INSTM.Bulletin.v31.773〉 (Accessed date: 07 August 2024). Amani, P., Karakashev, S.I., Grozev, N.A., Simeonova, S.S., Miller, R., Rudolph, V., et al., 2021. Effect of selected monovalent salts on surfactant stabilized foams. Adv Colloid Interface Sci 295, 102490. https://doi.org/ 10.1016/j.cis.2021.102490. Kaptay, G., 2003. Interfacial criteria for stabilization of liquid foams by solid particles. Colloids Surf A: Physicochem Eng Asp 230 (1-3), 67–80. https://doi. org/10.1016/j.colsurfa.2003.09.016. Zhan, F., Youssef, M., Shah, B.R., Li, J., Li, B., 2022. Overview of foam system: natural material-based foam, stabilization, characterization, and applications. Food Hydrocoll 125, 107435. https://doi.org/10.1016/j.foodhyd.2021.107435. Zhang, Y., Liu, Q., Ye, H., Yang, L., Luo, D., Peng, B., 2021. Nanoparticles as foam stabilizer: mechanism, control parameters and application in foam flooding for enhanced oil recovery. J Pet Sci Eng 202, 108561. https://doi.org/10.1016/j. petrol.2021.108561. Rio, E., Drenckhan, W., Salonen, A., Langevin, D., 2014. Unusually stable liquid foams. Adv Colloid Interface Sci 205, 74–86. https://doi.org/10.1016/j. cis.2013.10.023. Burdette, T.C., Bramblett, R.L., Deegan, A.M., Coffey, N.R., Wozniak, A.S., Frossard, A.A., 2022. Organic signatures of surfactants and organic molecules in surface microlayer and subsurface water of delaware bay. ACS Earth Space Chem. https://doi.org/10.1021/acsearthspacechem.2c00220. Dickinson, E., 2010. Food emulsions and foams: stabilization by particles. Curr Opin Colloid Interface Sci 15 (1-2), 40–49. https://doi.org/10.1016/j. cocis.2009.11.001. Kruglyakov, P.M., Elaneva, S.I., Vilkova, N.G., 2011. About mechanism of foam stabilization by solid particles. Adv Colloid Interface Sci 165 (2), 108–116. https://doi.org/10.1016/j.cis.2011.02.003. El Zrelli, R., Baliteau, J.Y., Yacoubi, L., Castet, S., Grégoire, M., Fabre, S., et al., 2021. Rare earth elements characterization associated to the phosphate fertilizer plants of Gabes (Tunisia, Central Mediterranean Sea): geochemical properties and behavior, related economic losses, and potential hazards. Sci Total Environ 791, 148268. https://doi.org/10.1016/j.scitotenv.2021.148268. Chuang, C.Y., Santschi, P.H., Wen, L.S., Guo, L., Xu, C., Zhang, S., et al., 2015. Binding of Th, Pa, Pb, Po and Be radionuclides to marine colloidal macromolecular organic matter. Mar Chem 173, 320–329. https://doi.org/ 10.1016/j.marchem.2014.10.014. Santschi, P.H., Xu, C., Zhang, S., Schwehr, K.A., Lin, P., Yeager, C.M., et al., 2017. Recent advances in the detection of specific natural organic compounds as carriers for radionuclides in soil and water environments, with examples of radioiodine and plutonium. J Environ Radioact 171, 226–233. https://doi.org/ 10.1016/j.jenvrad.2017.02.023. Zhao, L.Y.L., Schulin, R., Weng, L., Nowack, B., 2007. Coupled mobilization of dissolved organic matter and metals (Cu and Zn) in soil columns. Geochim Cosmochim Acta 71, 3407–3418. https://doi.org/10.1016/j.gca.2007.04.020. Jing, X., Luo, Q., Cui, X., Wang, Q., Liu, Y., Fu, Z., 2022. Molecular dynamics simulation of CO2 hydrate growth in salt water. J Mol Liq 366, 120237. https:// doi.org/10.1016/j.molliq.2022.120237. Robert, D., Crisp, J., 2021. Protest movements against industry-related environmental burdens and territorial injustice in Gabès and Kerkennah (Tunisia). Spatial justice, 16. Available from: 〈https://www.jssj.org/wp-content/ uploads/2021/07/JSSJ_16_Robert_EN.pdf〉 (Accessed date: 07 August 2024). Kwas, H., Rangareddy, H., Rajhi, H.H., 2024. Impact of outdoor air pollutants exposure on the severity and outcomes of community-acquired pneumonia in Gabes Region, Tunisia. Cureus 16 (8), e66578. https://doi.org/10.7759/ cureus.66578. Nacef, T., Nasfi, F., Bouacha, H., Hessairi, M., 1989. La pollution aérienne à Gabès, ses effets sur la santé des écoliers. Tunis Med 67 (12), 785–794.