Copper Tailing Flocculation in Seawater: Relating the Yield Stress with Fractal Aggregates at Varied Mixing Conditions

: The implications of physical conditions of the feedwell on the rheological properties of synthetic copper tailings, ﬂocculated in seawater, were analysed. The mixing intensity of ﬂocculation was related to the structural characteristics of the aggregates, and the outcomes were linked to the yield stress of the pulp sediments. Tailings settling assays were conducted by using a 30 mm turbine type stirrer with an in-situ aggregate size characterisation. The structural characteristics of the aggregates were determined by using the focused beam reﬂectance measurement (FBRM). After a mixing time between the pulp and the ﬂocculant, the sample was allowed to settle for 2.5 h, where the variation of the sediment height was minimal. The sediment was gently removed and subjected to rheological characterisation. The yield stress was measured on an Anton Paar MCR 102 rheometer (ANAMIN Group, Santiago, Chile), with a vane-in-cup conﬁguration. The mixing intensity was related to the characteristics of the aggregates, and the outcomes were linked to the yield stress of the ﬂocculated pulp sediments. More aggressive hydrodynamics deteriorated the structure of the aggregates, promoting the reduction of both its size and the fractal dimension. This brought direct consequences on the rheological properties of the sediments: at higher mixing level, the yield stress was lower. The explanation lies in the structural changes of the aggregates, where at a ﬁxed mixing rate, the yield stress presented a seemingly exponential increase over the fractal dimension. Additionally, correlations were found between the rheological properties with settling rate and aggregate size. The composition of Cl and was / The of determined by acid–base volumetry, was showed of / at ◦ The pH of slurries was ﬁxed at pH with


Introduction
The recovery and succeeding reuse of processed water is essential to guarantee the sustainability of the copper mining industry, and it is directly associated with operational and environmental costs It is interesting to comprehend the interplay of concentrated slurries in saline environments since the use of seawater is an option that is being increasingly used by the copper mining industry. Salinity alone is not a significant problem for concentration operations; when working at natural pH (pH 8 approx.), impressive recoveries of copper and molybdenum can be achieved in flotation stages [21][22][23]. It is noteworthy considering that a desalination plant by reverse osmosis leads to diverse financial and environmental concerns. For example, reverse osmosis (RO) plants commonly return hypersaline water to the sea, where the consequences to the marine life can be significant considering the high concentrations of ions and precipitating complexes (CaCO 3 , CaSO 4 , BaSO 4 ). Additionally, the high temperature of the stream (30-40 • C) increases the toxicity of metals and chemical [24,25].
Recently, Reyes et al. [26] proved the notable implication of seawater on the yield stress of magnetite tailings. For this, they used mixtures of seawater with process water, and the results were associated with measurements of zeta potential. Jeldres et al. [27] studied the yield stress and viscoelasticity of kaolin sediments of slurries flocculated with an anionic polyelectrolyte of high molecular weight. The authors found that pulps fixed with low salinity led to a rise in sediment yield stress, while increasing the NaCl concentration the yield stress diminished. This behaviour was interpreted in terms of the DLVO theory, where compressions of the ionic cloud that surround the particles surfaces and functional groups of the polymer produced aggregates with a weaker structure at high salinity. This phenomenon was shown later by molecular dynamics simulations, considering the interaction between a polyacrylamide with a quartz surface, varying the salt concentrations [28].
To date, several studies show the influence of mixing degree during flocculation on suspension dewatering [29][30][31] and, consequently, the relationship between the structure of the aggregates with the settling rate [32,33]. However, few systematic studies deepen in the relationship between the microscopic properties over the rheological behaviour of the underflow in thickening stages. The most relevant work was recently developed by Benn et al. [34], who performed turbulent flocculation in a large-scale pipe reactor (slurry flows 25 L min −1 ). The authors analysed the yield stress profile from the bottom to the top of the sediment and related the outcomes to the solids concentration. Additional information was presented on the influence of flocculant management. In this context, our research aims to analyse the implications of physical conditions of the feedwell on the rheological properties of the underflow of flocculated clay-based tailings in seawater. The mixing intensity was related to the structural characteristics of the aggregates, and the outcomes were linked to the yield stress of the flocculated pulp sediment. The results provided in this study are of interest to the copper mining industry, especially those that use seawater in their operations.

Materials
Quartz particles purchased from DondeCapo (Santiago, Chile) were used, where the SiO 2 content estimated by X-Ray diffraction (XRD) analysis was greater than 99 wt.%. A D5000 X-ray diffractometer (Siemens S.A., Lac Condes, Chile) was used and the data were processed with Total Pattern Analysis Software (TOPAS) (Siemens S.A., Lac Condes, Chile). The mean square chord length, determined by the focused beam reflectance measurement (FBRM) probe, was 24 µm. The kaolin particles were acquired from Ward's Science (Rochester NY, USA), and the XRD indicated a composition of 84% Kaolinite and 16% Halloysite. The mean square chord length was 12 µm.
Seawater was collected from San Jorge Bay in Antofagasta (Chile) and was filtered at 1 µm using a UV filter system (BIO&NAM, Santiago, Chile), to defeat its bacterial activity. The ionic composition was determined with varied chemical methods depending on the ion type. By atomic absorption spectrometry: Na + 10.9 g/L, Mg 2+ 1.38 g/L, Ca 2+ 0.4 g/L, K + 0.39 g/L. The composition of Cl − was determined by argentometric method and it was 19.6 g/L. The concentration of HCO 3− was determined by acid-base volumetry, whose value was 0.15 mg/L. Seawater showed a conductivity of 50.4 mS/cm at 25 • C. The pH of the slurries was fixed at pH 8 with sodium hydroxide.

Flocculation Test: Aggregate Size and Settling
Flocculation experiments were conducted using a PTFE, 30 mm-diameter turbine type impeller at the end of vertical shaft (diameter 4 mm) within a 1 L capacity, 100 mm diameter vessel (see Figure 1). The base of the impeller was positioned 20 mm above the vessel's bottom surface. A volume of liquid (232.2 mL) was transferred to the vessel, with known masses of the solid phases then added as dry powders to give mixtures containing 92.3 wt.% quartz and 7.7 wt.% kaolin at the required total solid concentration. The suspension was preconditioned by mixing at 450 rpm with a laboratory stirring IKA EUROSTAR 60 Control (WReichmann, Santiago, Chile) for 30 min. After this period, the mixing rate was reduced according to the test (from 100 to 300 rpm), and diluted flocculant was added to give the required dosage, with the mixing maintained for a further required flocculation time (from 10 to 80 s). The flocculated suspension was then carefully transferred into a stoppered cylinder (300 cm 3 volume, 35 mm internal diameter), which was inverted twice by hand (the entire cylinder rotation process taking~4 s in all cases). The fall of the mudline, defined as the suspension-supernatant interface, was monitored for one hour and the settling rate taken as the initial linear slope of a plot of the mudline height against time. Flocculation experiments were conducted using a PTFE, 30 mm-diameter turbine type impeller at the end of vertical shaft (diameter 4 mm) within a 1 L capacity, 100 mm diameter vessel (see Figure  1). The base of the impeller was positioned 20 mm above the vessel's bottom surface. A volume of liquid (232.2 mL) was transferred to the vessel, with known masses of the solid phases then added as dry powders to give mixtures containing 92.3 wt.% quartz and 7.7 wt.% kaolin at the required total solid concentration. The suspension was preconditioned by mixing at 450 rpm with a laboratory stirring IKA EUROSTAR 60 Control (WReichmann, Santiago, Chile) for 30 min. After this period, the mixing rate was reduced according to the test (from 100 to 300 rpm), and diluted flocculant was added to give the required dosage, with the mixing maintained for a further required flocculation time (from 10 to 80 s). The flocculated suspension was then carefully transferred into a stoppered cylinder (300 cm 3 volume, 35 mm internal diameter), which was inverted twice by hand (the entire cylinder rotation process taking ~4 s in all cases). The fall of the mudline, defined as the suspensionsupernatant interface, was monitored for one hour and the settling rate taken as the initial linear slope of a plot of the mudline height against time. The size distribution of the primary and aggregate particles was recorded using the Particle Track G400 equipment that uses Mettler Toledo's focused beam reflectance measurement (FBRM) technology. The measurement began with the tracking of the primary particles of the copper tailings and culminated three minutes after the polymer addition.

Yield Stress of Sediment
After 2.5 h of sedimentation, the supernatant liquid was gently removed and the lower cylinder separated (containing the sediment of the column). The yield stress was determined by employing the vane-in-cup configuration, with the stress-strain method. The measurements were carried out on an Anton Paar MCR 102 rheometer (ANAMIN Group, Santiago, Chile) and the RheocompassTM Light version software processed the data (ANAMIN Group, Santiago, Chile). The diameter of the vane was 2.2 cm and that of the cup was 4.2 cm.

Fractal Dimension
The aggregates usually have scale-invariant properties, allowing them to be characterised by a simple parameter: the fractal dimension. This parameter is directly linked to the porosity of the aggregates, and its value may range from Df = 1 (for a one-dimensional line) up to 3 (solid body of three dimensions). As for flocculation, an aggregate with fractal dimension Df = 3 would match to a solid sphere. The calculation of the fractal dimension is made with Equation (1) [33] and requires entering the values of the hindered settling rate and aggregate size. Both 'input' variables were obtained through the experiments described in Section 2.2. The size distribution of the primary and aggregate particles was recorded using the Particle Track G400 equipment that uses Mettler Toledo's focused beam reflectance measurement (FBRM) technology. The measurement began with the tracking of the primary particles of the copper tailings and culminated three minutes after the polymer addition.

Yield Stress of Sediment
After 2.5 h of sedimentation, the supernatant liquid was gently removed and the lower cylinder separated (containing the sediment of the column). The yield stress was determined by employing the vane-in-cup configuration, with the stress-strain method. The measurements were carried out on an Anton Paar MCR 102 rheometer (ANAMIN Group, Santiago, Chile) and the RheocompassTM Light version software processed the data (ANAMIN Group, Santiago, Chile). The diameter of the vane was 2.2 cm and that of the cup was 4.2 cm.

Fractal Dimension
The aggregates usually have scale-invariant properties, allowing them to be characterised by a simple parameter: the fractal dimension. This parameter is directly linked to the porosity of the aggregates, and its value may range from D f = 1 (for a one-dimensional line) up to 3 (solid body of three dimensions). As for flocculation, an aggregate with fractal dimension D f = 3 would match to a solid sphere. The calculation of the fractal dimension is made with Equation (1) [33] and requires entering the values of the hindered settling rate and aggregate size. Both 'input' variables were obtained through the experiments described in Section 2.2.
(1)  Figure 2 shows the chord length distributions (CLD) of aggregates composed of synthetic tailings in seawater at pH 8. The square-weighted distributions reflect a volume-weighting and will be more influenced by the contribution of any aggregates formed. FBRM provides a distribution of particle chords, unlike a conventional distribution based on diameters, however Heath et al. [35] proved that the mean or mode averages of the square-weighted chord lengths are comparable to conventional sizing over the range from 50 to 400 µm. The onset of fine particle aggregation from flocculation would typically be expected to induce a shift to the right in CLDs; however, Figure 2 displayed a variation to the left while for increased mixing time. This is explained because of the irreversible fragmentation of the aggregates while they are subjected to the hydrodynamics shearing, and the finite duration that flocculants remain adsorbed on the particle surfaces. It is expected that the breakage will be more intense in larger structures; therefore, the differences are attenuated after a certain period when few large-aggregates remain in the system. For this reason, the shift to the left of the size distribution is attenuated as time progresses.  Figure 2 shows the chord length distributions (CLD) of aggregates composed of synthetic tailings in seawater at pH 8. The square-weighted distributions reflect a volume-weighting and will be more influenced by the contribution of any aggregates formed. FBRM provides a distribution of particle chords, unlike a conventional distribution based on diameters, however Heath et al. [35] proved that the mean or mode averages of the square-weighted chord lengths are comparable to conventional sizing over the range from 50 to 400 µ m. The onset of fine particle aggregation from flocculation would typically be expected to induce a shift to the right in CLDs; however, Figure 2 displayed a variation to the left while for increased mixing time. This is explained because of the irreversible fragmentation of the aggregates while they are subjected to the hydrodynamics shearing, and the finite duration that flocculants remain adsorbed on the particle surfaces. It is expected that the breakage will be more intense in larger structures; therefore, the differences are attenuated after a certain period when few large-aggregates remain in the system. For this reason, the shift to the left of the size distribution is attenuated as time progresses.

Analysis of Yield Stress
The vigour at which the slurry is mixed with the polymers is determined by the agitation rate, quantified by the impeller revolutions per minute (rpm) and the mixing time. These variables are crucial for the aggregation since they allow to ensure particle-particle collisions, and they form the basis for the proper design of thickener feedwells [36]. Figure 3 presents the effect of the mixing intensity through flocculation on the yield stress of the sediment separated after 2.5 h that the tailings sediment. The trend given is clear: higher mixing intensity leads to lower yield stress values. For example, the yield stress taken at 100 rpm for 10 s was higher than 300 Pa, while after an intensive stirring of 300 rpm for 80 s, the value was close to 150 Pa. Interestingly, no relevant differences in the solid concentration of sediment were distinguished concerning the various disturbances imposed to the system and mostly oscillates between 42-44 wt. % (see Table 1), without showing clear trends regarding the mixing variables. Accordingly, the performance of the yield stress observed in Figure 3 is explained by the structural changes of the particle aggregates, as a product of changing the mixing intensity. Stronger agitations cause a more significant breaking of bonds between the particle/aggregate networks that make up the slurry, before the sedimentation process. This causes the aggregate structures to be more fragile and open that is consistent with the characterisation of the fractal dimension, which is described later. The vigour at which the slurry is mixed with the polymers is determined by the agitation rate, quantified by the impeller revolutions per minute (rpm) and the mixing time. These variables are crucial for the aggregation since they allow to ensure particle-particle collisions, and they form the basis for the proper design of thickener feedwells [36]. Figure 3 presents the effect of the mixing intensity through flocculation on the yield stress of the sediment separated after 2.5 h that the tailings sediment. The trend given is clear: higher mixing intensity leads to lower yield stress values. For example, the yield stress taken at 100 rpm for 10 s was higher than 300 Pa, while after an intensive stirring of 300 rpm for 80 s, the value was close to 150 Pa. Interestingly, no relevant differences in the solid concentration of sediment were distinguished concerning the various disturbances imposed to the system and mostly oscillates between 42-44 wt.% (see Table 1), without showing clear trends regarding the mixing variables. Accordingly, the performance of the yield stress observed in Figure  3 is explained by the structural changes of the particle aggregates, as a product of changing the mixing intensity. Stronger agitations cause a more significant breaking of bonds between the particle/aggregate networks that make up the slurry, before the sedimentation process. This causes the aggregate structures to be more fragile and open that is consistent with the characterisation of the fractal dimension, which is described later.   Table 1. Figure 4 displays the relationship between the yield stress of the sediment and the average size of the aggregates that are formed after the flocculant addition. For a fixed mixing rate, a larger aggregate size leads to higher yield stress. The trend is similar for the all mixing rates examined in this research (from 100 to 300 rpm). However, care should be taken when relating only the aggregate size, since it is essential to consider the physical conditions in which they were formed. For example, aggregate with a mean chord length of 100 microns causes higher yield stress if it is built at a mixing rate of 300 rpm instead of 150 rpm. The reagent dosage is also relevant: at a higher amount of polymer, the yield stress is larger. This is illustrated in Figures 7 and 8.
The relationship between aggregate size and yield stress coincides with that proposed by Neelakantal et al. [20], who studied the effect of shearing kaolinite sediments, flocculated with polyacrylamide copolymers. The authors based their explanation on the elimination of retained water within the porous aggregates as a result of the fragmentation caused by shearing. This would create a dramatic reduction in the apparent volume of solids in the suspension. However, unlike the assays executed by Neelakantal et al., our work based on studying the mixing intensity on flocculation, that is, before the sedimentation begins; therefore, the structural behaviour between both studies differs.  Figure 4 displays the relationship between the yield stress of the sediment and the average size of the aggregates that are formed after the flocculant addition. For a fixed mixing rate, a larger aggregate size leads to higher yield stress. The trend is similar for the all mixing rates examined in this research (from 100 to 300 rpm). However, care should be taken when relating only the aggregate size, since it is essential to consider the physical conditions in which they were formed. For example, aggregate with a mean chord length of 100 microns causes higher yield stress if it is built at a mixing rate of 300 rpm instead of 150 rpm. The reagent dosage is also relevant: at a higher amount of polymer, the yield stress is larger. This is illustrated in Figures 7 and 8.
The relationship between aggregate size and yield stress coincides with that proposed by Neelakantal et al. [20], who studied the effect of shearing kaolinite sediments, flocculated with polyacrylamide copolymers. The authors based their explanation on the elimination of retained water within the porous aggregates as a result of the fragmentation caused by shearing. This would create a dramatic reduction in the apparent volume of solids in the suspension. However, unlike the assays executed by Neelakantal et al., our work based on studying the mixing intensity on flocculation, that is, before the sedimentation begins; therefore, the structural behaviour between both studies differs.   Table 1. Figure 5 gives the connection between fractal dimension and yield stress when the flocculation is carried out at different mixing rates. An apparent decrease in the fractal dimension is generated with the increase of the mixing, which is agreeable with that declared by Deng and Davè [37], who state that when aggregates are produced with a higher shear rate, they grow irregularly, resulting in lower fractal dimensions. On the other hand, more aggressive hydrodynamics deteriorates the structure of the aggregates, promoting the reduction of the fractal dimension.
At a fixed mixing rate, the yield stress presented a seemingly exponential increase over the fractal dimension. This agreement with the proposed by Deng and Davè [37], who established that a lower fractal dimension leads to a lower mechanical strength of the agglomerate, which was determined by the model proposed by Kendall and Stainton [38]. The authors considered the porosity, adhesion work and particle diameter. Zhou et al. presented two interesting works that directly links the rheological properties with the interaction forces [18,39]. The authors used an atomic force microscope to record the interaction forces between silica surfaces in solutions with cationic polyelectrolytes of different electric charge. The magnitude of adhesion force was related to the strength of the particle networks that make up a silica slurry. Interestingly, it was found that the strength of aggregates has a direct relationship with the rheological properties of the pulps.   Figure 5 gives the connection between fractal dimension and yield stress when the flocculation is carried out at different mixing rates. An apparent decrease in the fractal dimension is generated with the increase of the mixing, which is agreeable with that declared by Deng and Davè [37], who state that when aggregates are produced with a higher shear rate, they grow irregularly, resulting in lower fractal dimensions. On the other hand, more aggressive hydrodynamics deteriorates the structure of the aggregates, promoting the reduction of the fractal dimension.
At a fixed mixing rate, the yield stress presented a seemingly exponential increase over the fractal dimension. This agreement with the proposed by Deng and Davè [37], who established that a lower fractal dimension leads to a lower mechanical strength of the agglomerate, which was determined by the model proposed by Kendall and Stainton [38]. The authors considered the porosity, adhesion work and particle diameter. Zhou et al. presented two interesting works that directly links the rheological properties with the interaction forces [18,39]. The authors used an atomic force microscope to record the interaction forces between silica surfaces in solutions with cationic polyelectrolytes of different electric charge. The magnitude of adhesion force was related to the strength of the particle networks that make up a silica slurry. Interestingly, it was found that the strength of aggregates has a direct relationship with the rheological properties of the pulps.   Table 1. Figure 5 gives the connection between fractal dimension and yield stress when the flocculation is carried out at different mixing rates. An apparent decrease in the fractal dimension is generated with the increase of the mixing, which is agreeable with that declared by Deng and Davè [37], who state that when aggregates are produced with a higher shear rate, they grow irregularly, resulting in lower fractal dimensions. On the other hand, more aggressive hydrodynamics deteriorates the structure of the aggregates, promoting the reduction of the fractal dimension.
At a fixed mixing rate, the yield stress presented a seemingly exponential increase over the fractal dimension. This agreement with the proposed by Deng and Davè [37], who established that a lower fractal dimension leads to a lower mechanical strength of the agglomerate, which was determined by the model proposed by Kendall and Stainton [38]. The authors considered the porosity, adhesion work and particle diameter. Zhou et al. presented two interesting works that directly links the rheological properties with the interaction forces [18,39]. The authors used an atomic force microscope to record the interaction forces between silica surfaces in solutions with cationic polyelectrolytes of different electric charge. The magnitude of adhesion force was related to the strength of the particle networks that make up a silica slurry. Interestingly, it was found that the strength of aggregates has a direct relationship with the rheological properties of the pulps.  The sedimentation rate of a suspension is strongly determined by the size and structure of the aggregates that compose it. For the same size, more compact structures have higher settling rates than open and porous have [40]. In this sense, it is expected that aggregates with larger fractal dimensions settle faster. In parallel, the structural aggregate features define the rheological behaviour of a suspension, so it is assumed that there is a close correlation among the settling rate and the rheology of the sediment, as shown in Figure 6: An increase in the settling rate was accompanied by higher yield stress. This result has a significant practical utility since there are continuous attempts in the mining industry to optimise the recovery of water from thickening stages, whose search generally lies in raising the sedimentation rate while maintaining low yield stress values in the thickener underflow. Then, to sustain proper thickener control, it is crucial to understand that both parameters are closely linked and the decisions that operators make inevitably impact both. The sedimentation rate of a suspension is strongly determined by the size and structure of the aggregates that compose it. For the same size, more compact structures have higher settling rates than open and porous have [40]. In this sense, it is expected that aggregates with larger fractal dimensions settle faster. In parallel, the structural aggregate features define the rheological behaviour of a suspension, so it is assumed that there is a close correlation among the settling rate and the rheology of the sediment, as shown in Figure 6: An increase in the settling rate was accompanied by higher yield stress. This result has a significant practical utility since there are continuous attempts in the mining industry to optimise the recovery of water from thickening stages, whose search generally lies in raising the sedimentation rate while maintaining low yield stress values in the thickener underflow. Then, to sustain proper thickener control, it is crucial to understand that both parameters are closely linked and the decisions that operators make inevitably impact both.  Table 1.
The flocculant dose was a key part of all the cases examined. An extract of the experimental data developed in the present work is shown in Figure 7, which shows a clear tendency of the rheology regarding the dosage of the polymer. A higher amount of polymer molecules leads to a greater amount of bonds that bind the particles, generating higher yield stress. On the other hand, the higher polymer dosage produces aggregates with a more robust structure, causing an increase in its fractal dimension (Figure 8), as previously reported by Ofori et al. [41].  The flocculant dose was a key part of all the cases examined. An extract of the experimental data developed in the present work is shown in Figure 7, which shows a clear tendency of the rheology regarding the dosage of the polymer. A higher amount of polymer molecules leads to a greater amount of bonds that bind the particles, generating higher yield stress. On the other hand, the higher polymer dosage produces aggregates with a more robust structure, causing an increase in its fractal dimension (Figure 8), as previously reported by Ofori et al. [41]. The sedimentation rate of a suspension is strongly determined by the size and structure of the aggregates that compose it. For the same size, more compact structures have higher settling rates than open and porous have [40]. In this sense, it is expected that aggregates with larger fractal dimensions settle faster. In parallel, the structural aggregate features define the rheological behaviour of a suspension, so it is assumed that there is a close correlation among the settling rate and the rheology of the sediment, as shown in Figure 6: An increase in the settling rate was accompanied by higher yield stress. This result has a significant practical utility since there are continuous attempts in the mining industry to optimise the recovery of water from thickening stages, whose search generally lies in raising the sedimentation rate while maintaining low yield stress values in the thickener underflow. Then, to sustain proper thickener control, it is crucial to understand that both parameters are closely linked and the decisions that operators make inevitably impact both.  Table 1.
The flocculant dose was a key part of all the cases examined. An extract of the experimental data developed in the present work is shown in Figure 7, which shows a clear tendency of the rheology regarding the dosage of the polymer. A higher amount of polymer molecules leads to a greater amount of bonds that bind the particles, generating higher yield stress. On the other hand, the higher polymer dosage produces aggregates with a more robust structure, causing an increase in its fractal dimension (Figure 8), as previously reported by Ofori et al. [41].

Conclusions
This research aims to analyse the implications of physical conditions of the feedwell on the rheological properties of the underflow of flocculated clay-based tailings in seawater. The mixing intensity was related to the structural characteristics of the aggregates, and the outcomes were linked to the yield stress of the flocculated pulp sediment. Stronger agitations cause a more significant breaking of bonds between the particle/aggregate networks that make up the slurry, before the sedimentation process. This induces the aggregate structures to be more fragile and open; then, higher mixing intensity leads to lower yield stress values. In this context, at a fixed mixing rate, the yield stress presented a seemingly exponential increase over the fractal dimension. Additionally, a direct relationship between the yield stress and the settling rate was found since both parameters depend on the structural characteristic of the particle aggregates. The outcomes provided in this study are of interest to the copper mining industry, especially those that use seawater in their operations.