Near-Field Analysis of Turbidity Flows Generated by Polymetallic Nodule Mining Tools

: The interest in polymetallic nodule mining has considerably increased in the last few decades. This has been largely driven by population growth and the need to move towards a green future, which requires strategic raw materials. Deep-Sea Mining (DSM) is a potential source of such key materials. While harvesting the ore from the deep sea by a Polymetallic Nodule Mining Tool (PNMT), some bed sediment is unavoidably collected. Within the PNMT, the ore is separated from the sediment, and the remaining sediment–water mixture is discharged behind the PNMT, forming an environmental concern. This paper begins with surveying the state-of-the-art knowledge of the evolution of the discharge from a PNMT, in which the discharge characteristics and generation of turbidity currents are discussed. Moreover, the existing water entrainment theories and coefﬁcients are analyzed. It is shown how plumes and jets can be classiﬁed using the ﬂux balance approach. Following that, the models of Lee et al. (2013) and Parker et al. (1986) are combined and utilized to study the evolution of both the generated sediment plume and the subsequent turbidity current. The results showed that a smaller sediment ﬂux at the impingement point, where the plume is transformed into a turbidity current, results in a shorter run-out distance of the turbidity current, consequently being more favorable from an environmental point of view.


Introduction
The industrial revolution and the rising living standards around the world require an unconventional search for resources for many key materials (e.g., Cobalt (Co), Nickel (Ni), Tellurium (Te), Titanium (Ti), Platinum (Pt), and rare Earth metals). These raw materials are vital to many modern industries (e.g., wind turbines, telecommunications, and computers) [1,2]. However, the continuous decrease of the grade of these minerals on land makes it worthwhile for mining companies to go into the deep ocean, which has the advantage of having better ore grades [3], for instance researchers estimated that the largest reserves of cobalt, nickel, and manganese are found on the sea bed [4]. However, Deep-Sea Mining (DSM) requires the consideration of several technical aspects, and these should be compatible with the targeted mineral deposits, seabed conditions, water depth, collection method, and deep-sea fauna [5] (see Figure 1).
The DSM system for polymetallic nodules consists of three parts, namely the Hydraulic Polymetallic Nodule Mining Tool (PNMT), the Vertical Transport System (VTS), and the Production Support Vessel (PSV). Deep-ocean polymetallic nodules form on or just below the abyssal plains of the ocean. An enormous amount of nodules is located on the ocean bed at a depth of approximately 4300-5500 [2]. Regarding polymetallic nodule mining operation, the PNMT collects the nodules from the sea floor and primarily separates the nodules from the excess water and fine sediments, which are discharged directly at the seafloor. Next, the VTS transports the collected nodules from the sea bed to the PSV at the sea surface, where the collected ore is dewatered. Finally, the VTS is used again to return the transported water, containing Sediments, Waste, and Other Effluents (SWOEs) to the deep sea. Figure 1 shows a typical scenario for a DSM system. There are two main sources of sediment plumes generated by DSM activities: (i) The discharge of the sediment-water mixture from a PNMT, which remains after the separation process [6]; (ii) The SWOE discharge from the VTS after the end of the dewatering and ore handling process on the PSV [7,8].
The released suspended sediment, whether it is from the VTS or from the PNMT, forms a sediment cloud termed a "sediment plume". This plume increases the turbidity level of deep-sea water [9,10], which negatively impacts the marine ecosystem [11]. In other words, fish behavior is affected, and the mortality rate of zooplankton species increases [11]. Additionally, the fresh deposited sediments from the discharged plume on the sea bed clog the feeding paths to sea bed benthic organisms. Evidently, DSM requires an accurate evaluation of its environmental consequences [2,12].  [12]).
The effect of sediment disturbance resulting from DSM applications on the local physicochemical conditions, i.e., soil composition and organic carbon content, of the disturbed area was investigated experimentally by [13]. The physical parameters (e.g., velocity and concentration profiles) of sediment plumes resulting from dredging operations were studied by [14,15]. The work of [7] documented numerical and experimental results of plume discharge from a VTS. A CFD analysis was performed by [6] for a horizontal discharge from a PNMT to investigate the effect of different discharge conditions (see the controlling parameters of a discharge process in Section 2.2) on the plume dispersion. A preprototype PNMT was developed and tested in the Belgian and German exploration areas; based on the results, a comprehensive monitoring program for the local environment was developed [16]. The work of [17,18] showed that the Particle Size Distribution (PSD) governs the traveling distance of fine suspended sediment. They found that the particle size changes as a result of the aggregation-breakup process between sediments (see Section 6 for more details about the flocculation process).
In the last few decades, valuable insight has been gained into DSM processes through various projects, where the environmental impacts of DSM were investigated, e.g., Refs. [5,9,10,13,17,[19][20][21][22][23]. Some of the DSM research projects are mentioned as follows: • The Deep Ocean Mining Environmental Study (DOMES) project was one of the leading projects, which aimed to obtain and investigate the necessary data for an independent impact assessment of DSM activities. DOMES was divided into two main phases as follows [19]: (i) Gaining quantitative data on the biological communities prior to mining and developing a framework to predict the impact of manganese nodule mining on the marine environment; (ii) Determining the accuracy of the environmental impact predictions obtained in the first phase through monitoring of pilot mining tests.
As a result of this project, a quantitative baseline of environmental parameters and a predictive framework were developed to characterize DSM environmental impacts. Furthermore, preliminary environmental guidelines for DSM were defined; • The Marine E-tech project took place in the Tropic Seamount in the north east Atlantic near the Canary Islands and aimed to study the environment sensitivity to Fe-Mn crust DSM operations. The Fe-Mn crusts are differentiated from polymetallic nodules as a result of different locations, depths, and DSM operation techniques. The local influencing variables (e.g., temperature, pressure) that govern the composition and formation of the Fe-Mn nodules were investigated. A Remotely Operated underwater Vehicle (ROV) was used to generate sediment plumes, and the plume dispersion was studied. The measurements showed that the plumes were significantly smaller than predicted because the effect of flocculation was not taken into account [17]; • The Towards Responsible Extraction Of Submarine Mineral Resources (TREASURE) project studied SWOE discharge from a VTS [22]. Within the scope of TREASURE, numerical simulations were performed using the drift-flux modeling approach (for details on drift-flux modeling, see [24]) to predict the VTS discharged plume characteristics (e.g., velocity and concentration) [25]. In addition, Reference [7] carried out detailed laboratory experiments to test different discharge parameters such as initial concentrations, momentum, and distance from the bed on the plume dispersion. The numerical results (e.g., velocity and concentration profiles) were compared against the experimental results, and a good agreement was found between them; • The PLUMEX2018 field experiments were conducted in the Southern California Bight at the beginning of 2018 [26]. The Multiresolution primitive equation regional ocean modeling system (MIT-MSEAS) model was used to predict the plume dispersion, i.e., direction, velocity, and concentration. Good agreement was found between the MIT-MSEAS model and the PLUMEX field experiments [23]; • The JPI-Oceans Mining Impact II research project aimed to develop a new framework for environmental monitoring and predictions for the environmental impacts of mining operations [27]. Global Sea Mineral Resources (GSR), a subsidiary of DEME group, designed and tested a PNMT in the Clarion-Clipperton Fracture Zone (CCFZ) within the Belgian license area to assess its environmental impact [16]. GSR performed detailed measurements on the current environment (climate, geomorphological, physical oceanographic, seabed substrate characteristics, natural hazards, noise, and light). Moreover, a biological baseline was assessed based on habitat heterogeneity. The potential environmental impact was divided into six categories: Biogeochemical changes of the sediment particles; (iv) Potential release of toxic sediment into the lower water column; (v) Emissions to the air; (vi) Natural hazards (weather condition, storms).
The objective of this paper is to present an analysis of the characterizations of plumes generated by a PNMT in terms of the transport, spreading, and deposition of sediments. This will help determine the optimal discharge conditions to minimize the environmental impact caused by sediment plumes. In this regard, the emphasis will be on the nearfield effects (see Section 2.1 for the nearfield/farfield definitions), where engineering and the design of equipment can significantly influence the spread and deposition of the generated plume.
This paper proceeds as follows. Section 2 provides an overview of the discharge process and the key parameters controlling it. In Section 3, we present the relevant length and time scales associated with polymetallic nodule mining process and an overview of sediment-laden jets and plumes. Section 4 presents the most important physical aspects of the jets, plumes, and turbidity currents analyzed herein. In Section 5, the model of [28] and the four-equation model of [29] are employed to study the mining-generated plume and subsequent turbidity currents. Section 6 discusses the flocculation effect on the discharged turbidity current. Finally, Section 7 provides key directions for future research.

Discharge Process
In this section, we first present an overview of the PNMT discharge process in order to study the turbidity flows generated by DSM. Secondly, a classification of the discharge properties and the physical parameters of the PNMT discharge process are described.

Overview
The main differences between the two discharge sources in a polymetallic nodule mining operation are the orientation (vertical from the VTS and mostly horizontal from the PNMT), distance to the seabed, flow rate, and concentration of the suspended sediments. Some studies showed that the SWOE discharge location ranges from near the seabed to just below the oxygen minimum layer (800-1000 m water depth) [5,7,8,21], while the distance to the seabed of a PNMT depends on the discharge position on the PNMT, ranging from 1-3 m. The SWOE discharge flow rate and concentration are estimated to be about 0.56 m 3 /s and 8 g/L, respectively [21,23]. Table 1 shows the estimated solid flux rate and sediment concentration range of a PNMT discharge. Bed disturbances mainly result from the movement of a PNMT and the pick-up process [13], which could be hydraulic, mechanical, or hybrid. A benthic disturbance experiment mimicking PNMT disturbances was conducted to investigate particles' resettlement on marine ecosystems during DSM operation [13]. They found that the particles migrate to the adjacent areas outside the mining zone, causing a change in the physicochemical bed properties. These migrated particles potentially clog the feeding paths of the benthic organisms.
Discharging these sediment-water mixtures without carefully optimizing the discharge parameters might unnecessarily enlarge the area affected by plume dispersion. A few studies are found in the current body of research about the discharge process from the PNMT, e.g., [6]. In this respect, we present in Section 5 our investigation of the sedimentwater discharge from a PNMT. We divided the horizontal discharge of a sediment-water mixture from a PNMT into four main parts of interest as follows (see Figure 2):

1.
Discharge source: This contains the initial conditions such as the momentum, concentration of suspended sediments, and distance from the sea bed z. The physical parameters depend on the design of the PNMT (e.g., methods of collection and separation); 2.
Jet or plume regime: In this region, depending on the flow discharge parameters, the flow can be a jet or plume. Later, when the buoyancy force is dominant, the flow becomes a plume (see Sections 2.2 and 4.2); 3.
Impingement region: This region is located on the sea bed. Here, the negative buoyant plume changes its direction due to the direct interaction with the seabed. Sediment deposition and possible sea bed erosion are expected to take place within this region; 4.
Turbidity current: This current is formed beyond the impingement region.  In the literature, classifying of sediment plumes is often related to length and time scales. Accordingly, terms such as nearfield or farfield regions are defined as follows: (see Nearfield region: This defined as the region close to the discharge apparatus, and it is mostly controlled by the discharge conditions. The flows in this region have a typical length scale up to few hundreds of meters and a time scale in the range of seconds to minutes; • Farfield region: This is defined as the region where the plume trajectory is dominated by the environmental parameters, such as the currents and seabed topology. The flows in this region have large time and length scales, which are typically in the range of days and kilometers, respectively.

Physical Parameters
The discharge properties can be divided into three main groups as follows [34] (see Concerning the mixture properties, the sediment-water discharge parameters (momentum, buoyancy, and volume) control the flow regime at the discharge source. The sediment plume dispersion is highly dependent on the ambient conditions, such as the turbulence level and background stratification, which are of high importance to quantify the plume dispersion within the nearfield/farfield regions. The discharge geometry and its orientation are two of the key aspects that need to be well considered during the design process of a PNMT.
The discharge source is defined as the release point of the sediment-water mixture at the back of a PNMT. The source condition is characterized by the momentum flux M o , buoyancy flux B o , volume flux Q o , and reduced gravity g o of a mixture. These parameters are determined as follows: where A is the cross-sectional area of the discharge geometry; ∆ρ = ρ − ρ o , where ρ o is the ambient density, ρ is the discharge density, and u o is the discharge velocity. The Boussinesq approach neglects the density difference, while the non-Boussinesq approach assumes a difference in density between the discharged and ambient flows. The turbulence level is indicated by the Reynolds number (Re). Furthermore, the ratio between the inertial force and the gravitational force represents the densimetric Froude number (Fr), which classifies the discharged flow into different flow regimes (Fr = 1, critical flow, Fr > 1, supercritical flow, Fr < 1, subcritical flow). The mixing process of a buoyant jet can be classified based on the Richardson number (Ri); if Ri > 1, the flow is dominated by buoyancy, whereas if Ri < 1, the flow is dominated by momentum [35].  [34].
where µ is the dynamic viscosity and L is the characteristic length, which is calculated based on the discharge geometry in the discharge origin (see Figure 2). A design framework characterizing the discharge process of the sediment-water mixture generated by PNMT is needed to analyze the sediment plume dispersion and the characteristics of the formed deposition layer. Besides the discharge defining, the nodule pick-up and separation processes have a great importance in controlling the initial conditions of the sedimentwater discharge parameters (see Equations (1)- (5)). These processes highly affect the mixture discharge, i.e., increased erosion of the sea bed leads to a high flux of sediments at the discharge source.

Flow Specification
Different flow regimes (e.g., jet, plume, and turbidity current) are expected in the nearfield area, and these depend on the discharge parameters described in the previous section. These flow regimes occur within certain length and time scales. We here present the expected length and time scales for the underlying physical processes associated with a polymetallic nodule mining process. Moreover, we provide an overview of the physics governing sediment particle suspensions and a brief review of sediment-laden jets and plumes.

Length and Time Scales
The characteristic length and time scales in DSM applications vary over a wide range of magnitudes. Figures 3 and 4 represent the expected length and time scales associated with a polymetallic nodule mining process. For a dimensional analysis, the most important length scales of a discharge process are the momentum length scale L m and the buoyancy length scale L b [34,36], where L m is defined as the distance where the momentum of the discharge is dominant, i.e., jet-like, and L b is defined as the distance where buoyancy is dominant, i.e., plume-like [37]. In other words, L b represents the distance at which the flow velocity decreases to the ambient velocity [34]. Regarding jet flow, there is an established jet core near the discharge source, which is not affected by the ambient entrainment. Therefore, Reference [36] divided the jet region into two zones: Zone of Flow Establishment (ZFE) and Zone of Established Flow (ZEF) (see Figure 6). A Lagrangian model was presented by [36] for a buoyant jet, which can predict the ZEF and ZFE. Far from the discharge source, buoyancy forces dominate the flow. Mining operations affect the mined areas; however, the turbidity current resulting from the discharge process might also affect the region beside the mined area, as it could potentially travel for a long distance. Recently, using an industrial-scale particle transport model, Reference [18] estimated this distance to be in the range of 4-9 km.

Particle Physics
Different ranges of particle size are encountered in DSM projects. This section outlines the main physics of a moving particle in a flow and the hindered settling concept. A single spherical particle settling in a quiescent and an infinite domain is governed by five forces: the drag F d , gravitational F g , buoyancy F b , added mass F a , and history F h forces [24,38] (see Figure 7).
The added mass force and history force can be omitted in the scope of this work, because these forces are considered unsteady terms and less significant compared to the other forces. Hence, the particle motion equation can be expressed as follows [24]: where A p and V p are the projected area of the particle and the particle size, respectively, C d is the drag coefficient, v p is the velocity of the particle, u is the carrier fluid velocity in the vertical direction, ρ 0 is the carrier fluid density, m p is the particle mass, and g is the gravitational acceleration. The drag coefficient has many formulas in the literature depending on the particle shape and roughness. These formulas were collected and classified according to the particle Reynolds number Re p [39]. The particle Reynolds number is the ratio between the particle inertial and viscous forces; moreover, it is a key parameter in determining particle motion. It is worth mentioning that these formulas were derived for spherical particles. The shape of the particle has a strong influence on the motion of the particle. A general expression for the settling velocity was derived by [40] as follows: where µ f is the kinematic viscosity, R = (ρ p − ρ)/ρ p , ρ p is the particle density, ρ is the carrier fluid density, C 1 is reported according to particle shape and lies between the value of (18 < C 1 < 24), and C 2 is the drag coefficient, C 2 = 0.4 for spherical particles and 1 < C 2 < 1.2 for nonspherical particles. Figure 7. Forces acting on a nonrotating settling particle in a quiescent flow, where fg is the gravitational force, fb is the buoyancy force, Fd is the drag force, and ws is the settling velocity (Goeree, 2016) [24].
The Stokes number (St) represents the ratio between the particle response time τ p (the time taken by a particle to adapt to the fluid motion) and the hydrodynamic time scale τ h (τ h = L/u, where L is the typical length scale and u is the flow velocity at this length scale).
In the case of small particles, St << 1, the particle is guided by the flow motion, whilst for big particles, St >> 1, the flow does not affect the particle motion [24].
Hindered settling represents the reduction of the settling velocity due to the interaction of the particle with the neighboring particles (i.e., collisions, displaced water, group settling). Based on sedimentation and fluidization experiments, Reference [41] calculated the hindered settling velocity for individual particle suspended within a mixture as follows: where V m (φ) is the hindered settling function depending on Re p and w s (φ) is the hindered settling velocity, which is related to the volumetric mixture concentration φ. According to [41], the hindered settling function considers the particle Reynolds number and the concentration value as follows: where n is called the Richardson and Zaki index, and it depends on the Reynolds number of the particles.
Polymetallic nodules are abundantly available in the Clarion-Clipperton Fracture Zone (CCFZ) region, which is located in the Pacific Ocean [21]. Sediment in the CCFZ consists of several fractions such as biogenic ooze, terrigenous or pelagic clay, volcano debris, hydrogenous material, and metalliferous sediment.
The Blue Nodules project is a Horizon 2020 project (EU Research and Innovation Program), which ran between February 2016 and July 2020. The project aimed to develop and test a novel PNMT in the CCFZ region. Within the scope of the Blue Nodules project, NTNU (Norwegian University of Science and Technology), a partner in the project, analyzed the obtained CCFZ sediment samples. The sediment compositions for different sites in the CCFZ are compared in Tables 2 and 3, where BC064 and BC062 are two box core samples obtained from the GSR contract area. Interoceanmetal Joint Organization (IOM) data were collected from IOM license area, which is located in the eastern part of the CCFZ [42,43]. For further details on CCFZ sediment composition, the reader is referred to the Blue Nodules project public reports and GSR public reports [16,[42][43][44].

Sediment-Laden Jets and Plumes
A fundamental understanding of the physical parameters of discharge, e.g., velocity, concentration, and turbulence, is needed to investigate the effect of the PNMT discharge on the sea environment. An increasing particle deposition rate and reducing plume dispersion rate would minimize the environmental impact. Bleninger et al. (2002) [46] and Neves et al. (2002) [47] carried out small-scale experiments on horizontal jet discharges to measure the particles' deposition. A large scatter in the deposition rate measurements was noted due to difficulties in ensuring a constant concentration discharge and steady-state flow rate during the discharge process. Lab experiments were conducted to study the sedimentation from horizontal and inclined buoyant jets in a stationary environment [48]. They introduced an integral model to determine the deposition behavior from an inclined, turbulent, buoyant jet. A dimensionless fall speed parameter, i.e., defined as the ratio between the settling velocity and the entrainment velocity, was incorporated into the model to measure the dependency of the source velocity on the deposition near the discharge source. The model was validated with laboratory experiments where the deposition behavior, plume shape, and sedimentation near the source were assessed. The earlier experiments on sediment-laden jets and plumes had the potential pitfall of not sustaining the steadiness of concentration and uniform shape at the jet exit [28].
A theoretical Lagrangian model of horizontal jets discharged in a stationary environment was presented by [28]. They investigated the deposition mechanisms of a horizontal particle-laden jet in terms of longitudinal distance and spreading angle. Their Lagrangian model was used by [49] to validate the experimental work of [28,48,50]. This model is adopted in this work in Section 5.
The impingement region ( Figure 2) entails a complex behavior of the flow due to the interaction between the discharged mixture and the ocean bed. At this region, the flow makes a turn and converts from a plume to a wall-bounded flow (i.e., turbidity current [51]), due to the presence of the ocean bed. This current propagates downstream of the ambient flow and interacts with the bed.

Turbidity Current
Turbidity currents can be described as particle-laden underflows, which are driven by the excess hydrostatic pressure resulting from the density difference between the ambient water and the sediment-water mixture [52,53]. Within this current, turbulence is developed due to the continuous motion of the current over the bed and the generated shear stresses as a result of the mixing process with the ambient water at the upper boundary of the current.
Turbidity currents can be triggered by several phenomena, such as internal waves or tides [54], river plumes [55], and breaching flows slides [56]. The mining process by a PNMT is also expected to trigger a turbidity current at the impingement region by the sediment-water discharge. Here, we refer to this current as a mining-generated turbidity current. The interaction between the discharge and the seabed determines the behavior of the generated turbidity current; the hydraulic properties prior to the impingement region are the main information to characterize the generated turbidity current.
Turbidity currents can occur in turbulent and laminar regimes. The flow regime is classified with the Reynolds number (flow thickness is the characteristic length scale), where Re < 1000 is laminar; above that, the flow is turbulent. The deposition and resuspension of sediment is highly affected by the turbulence structure [53]. Turbidity currents are divided into three parts: head, body, and tail (see Figure 8). From the perspective of hydraulics, the head has different properties than the region behind it (body and tail) [57,58]. The mass and momentum at the head differ significantly from the body and tail, as the head displaces the ambient fluid. Hence, the head is the densest part, as it experiences a friction resistance. Furthermore, the head is considered a "locus of erosion" [57,59], impacting the deposition and bed morphology. The head forms a nose starting at the lower region due to the no-slip boundary condition at the bed and the friction resistance at the upper region [60]. At the back of the head, vortices start to take place due to the effect of the velocity shear and turbulence in the ambient fluid. These vortices define the dynamics of the head and can be identified as Kelvin-Helmholtz instabilities [61]. As a result of these instabilities, the back of the head forms a sharp discontinuity in the thickness of the current. However, the average velocity of the body region has to be larger than the forward velocity of the head to achieve a constant rate of advance. [57]. The entrainment of the ambient fluid into the turbidity current and the change in the amount of suspended solids, due to net erosion or net deposition at the bed, depend on the characteristics of the current, and these also lead to changes in velocity, concentration, and particle size over time. Based on the velocity and concentration profiles, it is possible to derive characterizing layer-averaged parameters of a turbidity current. Parker et al. (1987) [62] derived a set of equations to estimate the layer-averaged characteristics for a turbidity current as follows: Many mathematical models are available to describe the phenomenon of turbidity current. These models can be divided into two main categories [57]. The first category is based on the vertical averaging technique, which describes the velocity and concentration with a single point (layer-averaged magnitude) along the traveling distance of the current.
The characterization of the structure of each model depends on the assumptions made for the conservation of these flow parameters. The conservation of momentum as one equation model was used in the work of [63], while [64] involved four equations, i.e., conservation of mass, momentum, sediment mass, and energy equations, to derive his model.
The second category drives the whole velocity and concentration profiles from a turbulence model. The equations of the total kinetic energy k and its dissipation have been used as a turbulence closure. Spatial averaging (large eddy simulation) can also be used as a closure for the turbulence [65][66][67][68].

Flow Physics
The literature reports many physical approaches that characterize the flow physics in the nearfield region. Here, we describe the most important physics governing the discharged flow in the nearfield region. We also summarize the differences between the Gaussian and top-hat approaches, which were used in previous studies of plume integral models. Moreover, we outline the flux balance approach, the entrainment theory, and the turbidity current.

Gaussian and Top-Hat Profiles
Many researchers [69][70][71][72][73][74] used integral models to describe the physics of plumes. Gaussian and top-hat profiles are the main assumptions for these models. The study of [75] extensively compared the governing equations of plumes using Gaussian and top-hat assumptions. The top-hat profile uses a constant velocity and concentration distribution over the entire plume domain, where the velocity and concentration are the average values of the dominant eddies (see Figure 9). Additionally, the large eddies within the plume are responsible for the overall transport of momentum and mass, while the small eddies have no effect. Top-hat profiles are used mostly in the Lagrangian approaches [36]. On the other hand, the Gaussian profile uses a Gaussian function to describe the velocity and concentration profiles, which are more suitable for Eulerian approaches [36]. Using the tophat assumption would lead to a larger entrainment coefficient compared to the Gaussian profile [76].
Describing the plumes' spreading behavior using the top-hat profiles offers an analytical simplicity because of the simple equations that can be used directly at the source and do not need any prior treatment for the initial condition (momentum, buoyancy, volume fluxes). Gaussian and top-hat approaches are crude assumptions to describe the selfsimilarity of the plumes [75]. Furthermore, the difference between these two approaches is small in terms of the predictions of the physical parameters. Moreover, the additional physical parameters that come within the Gaussian analysis have a minor role in describing the mean plume characteristics.
Various cases of turbulent jet geometries, e.g., slot and round shape, were studied in the work of [28]. They illustrated the governing equations that describe the jet propagation, i.e., velocity and concentration profiles along the jet width, for each of these cases. They provided an analysis for a turbulent jet from a slot in which the velocity and the concentration profiles are described b a Gaussian function (2D, x-y, Cartesian coordinates).

Flux Balance Approach
Flux balance parameter Γ is an analytical parameter that introduces a theoretical framework of plumes [78]. The flux balance parameter Γ is defined as the ratio of buoyancy force to inertia force [69]. The flux balance parameter can be expressed as follows [78]: where Γ o , B o , G o , and Q o are the initial flux balance, buoyancy, mass, and volume fluxes, respectively, α is the entrainment coefficient (see the entrainment theory in Section 4.3), and M B,o is the initial momentum flux under the Boussinesq approximation. Using the flux balance parameter, based on the discharge conditions, it is possible to classify the discharged flow into different regimes as follows: The plume is considered pure when its momentum and buoyancy fluxes are the same at the discharge source [79]. The plume is regarded as lazy when its buoyancy forces are dominant at the discharge source [80], while the plume is regarded as a forced plume (also known as a buoyant jet) when its momentum forces are dominant at the discharge source [70]. The pure jet remains a jet until it dissipates due to viscous diffusion [79].
Additionally, Γ can be introduced as a square ratio between the source length L s and jet length L m , i.e., the momentum-dominated region [81].
where, for Boussinesq plumes: and for non-Boussinesq plumes:

Water Entrainment Theory
Jets and plumes are shear flows generated by buoyant and inertial sources. The interaction of jets and plumes with the surroundings is known as the entrainment process. The development of a jet or a plume causes a velocity shearing that creates turbulent eddies that entrain the surrounding fluid [79]. The integral models of turbulent jets and plumes use the entrainment hypothesis as a closure relation for the turbulence. The entrainment hypothesis links the entrainment velocity u e , the rate at which ambient water is entrained across the edge of turbulent plume, to the characteristic velocity of the plume by a single coefficient of proportionality, α, which is called the entrainment coefficient [69,82,83]. The entrainment velocity is calculated as follows: where η is the ratio of the plume density to the ambient density. The classical approach of the entrainment hypothesis is based on the macroscopic conservation of momentum, mass, and volume fluxes. It has also a self-similar behavior of the turbulence in a dimensionless form over the downstream distance from the source [69]. A range of entrainment coefficients values (see Table 4) for different authors have been reported in the literature [84]. The entrainment coefficient ranges are 0.10 < α < 0.16 for plumes and 0.065 < α < 0.080 for jets (these values were calculated assuming a top-hat profile and a self-similar behavior of the plume) [84]. The variations in α occur due to the different setups of the experiments and the resulting experimental error. The systematic differences between the reported values of α for plumes and jets suggest that α depends on the ratio between buoyancy and inertia (Γ parameter; see Section 4.2). Table 4. Different entrainment coefficients values reported in the literature; "J" is the Jet-type flow. "P" is the Plume-type flow. "L" is Liquid. "P" is Pipe. "c" is Constriction. For the method, "w" is the hot wire probe and "l" is the laser Doppler anemometer. z * is the ratio z/d in which z is the distance from the source. d is the discharge diameter [84]. As mentioned above, the entrainment coefficient is the only parameter that represents the turbulence effect on the mean flow parameters. In this regard, Reference [85] proposed imposing restrictions on the entrainment coefficient by the mean kinetic energy equation. These restrictions are referred to as the entrainment relation. This relation couples the entrainment coefficient to the physical process, such as the buoyancy effect and turbulence production, while the entrainment models are closure relations, which are obtained after all the coefficients of the entrainment relation are parameterized [69,86].

Source
The mathematical models of [69,86] were developed by the integration of timeaveraged Navier-Stokes equations across a plane perpendicular to the mean flow direction (axial cross-section). The result of this integration is a system of coupled ordinary differential equations. The mathematical model of [69] depends on the volume flux, momentum flux, and buoyancy flux, while [86] based their model on the momentum, buoyancy, and mean kinetic energy, where the equations of [86] rely on the turbulence kinetic energy and Reynolds stresses. A few years later, Reference [87] presented a theory on an isolated plume in still air, and a new hypothesis was formulated in which the turbulence intensity directly affects the entrainment coefficient. Hence, a new equation was added to solve the turbulence kinetic energy. Utilizing the works of [69,86], the momentum, buoyancy, volume, and kinetic energy conservation equations were combined [88]. This was the first attempt to impose the constraints on α based on the conservation of kinetic energy.
Plumes and jets are conical shapes that develop in a self-similar fashion [78]. Rescaling dependent variables on the radial coordinate using characteristic scales such as velocity w m , buoyancy b m , and local width r m is a way to view the self-similarity [89]. The effective β entrainment radius and effective density parameter ∆ for Boussinesq and non-Boussinesq flows were introduced as follows [90]: A self-similar solution for the plume characteristics from the governing equations in the z coordinate of steady-state plumes with the top-hat profile was derived as follows [90]: A method to calculate the characteristic scales, which does not rely directly on the Gaussian shape assumption, but on the flow integral quantities, was proposed as follows [89]: M and Q are the momentum and volume fluxes, respectively, and B is buoyancy in the integral form [91].
Integrating the continuity equation from the work of [69] over the radial direction leads to the dilution in jets and plumes: Here, [ru] ∞ is the entrainment volume flux into a plume or a jet per unit height and ζ is a dimensionless vertical coordinate equal to z 0 r −1 m dz . From the entrainment assumption (u = αu e ), the next equation can be derived: Substitution of (31) into (30) leads to: Equation (32) shows that the entrainment coefficient is defined as half of the increase of the volume flux over a plume or jet radius r m . This relation calculates the amount of fluid entraining jets or plumes per unit r m , where more entrainment leads to higher values of α [89].

Assessment of the Nearfield Generated Plume and Turbidity Currents
In this section, we investigate the effect of all properties of the sediment-water mixture discharged from a PNMT on the hydrodynamics of the generated plume and turbidity current. This eventually aims at selecting the optimal discharge scenarios from an environmental point of view.
In our modeling approach, we considered that we have two separate, but connected regions. The first region entails the trajectory of the plume starting from the discharge source and ending at the impingement point, while the second region includes the generation and evolution of a turbidity current downstream the impingement point (see Figure 2). The governing equations of these two regions and the numerical results are presented in the following subsections.

Lagrangian Plume Model
The two-layer Lagrangian jet model, which was developed by [28], was used here to simulate the evolution of the jet/plume from the discharge source until the impingement point at the seabed. This model was validated by [49] against the experimental results of [28,48,50], justifying our choice of this model. The Lagrangian formulation approaches a changing cross-section and a coordinate system that both march in time with discrete time steps. The model follows a slicewise approach, meaning that there is a plume slice at each time step (see Figure 10), which has its unique characteristics, assuming a uniform top-hat profile.
The discharge trajectory depends on water entrainment, plume bending, and plume growth. The entrainment flux was calculated using Equation (34), and the entrainment coefficient α was calculated based on the densimetric Froude number Fr using Equation (33) [92]. The initial angle of the plume α p,0 at the discharge source is 0 • (with the horizontal axis), and the initial mass of the plume can be calculated using M p,0 = ρπb 2 0 h 0 , where a is the slice width and b 0 = a 0 = d 0 /2 is the slice length. The time step ∆t = 0.1a 0 /V 0 was calculated based on the initial velocity magnitude V 0 and initial diameter, and it was consistent throughout the marching procedure. The subscripts 0, i, and i + 1 denote the quantities of the initial condition, previous time step, and current time step, respectively. Following [49], the governing equations are expressed and solved in the following sequence: where u and w are the velocities in the x and z coordinates, respectively, and φ is the volumetric concentration.

Four-Equation Model for Turbidity Currents
The layer-averaged, four-equation model of [29] was developed by integrating the four conservation equations of momentum, fluid mass, sediment mass, and turbulent kinetic energy over the height of the turbidity current. The water entrainment, shear stress, and sediment entrainment were included as source terms. Following [93], the four equations can be expressed as follows: where h is the current height, α s is the bed slope angle, Ri is the Richardson number, U is the layer-averaged velocity, C is the layer-averaged concentration, K is the layer-averaged turbulent kinetic energy, 0 is the layer-averaged mean rate of turbulent energy dissipation, n is the bed porosity, e w is the water entrainment coefficient for turbidity currents, v er,n is the net erosion velocity perpendicular to the bed surface resulting from the combined effects of deposition and erosion, and w s is the hindered settling velocity. The bed shear velocity is calculated by the relation u * = √ α k K in which α k is a constant (0.05 to 0.5). The water entrainment coefficient e w is calculated as follows [29]: and 0 is calculated using the relation: where β 0 is a dimensionless parameter and C * f is a pseudo friction coefficient. The erosion velocity for clay is calculated using the relation presented in [94] as follows.
where c v is the consolidation/swelling coefficient, α m is a factor equal to 10, c u is the undrained shear rate, τ b is the bed shear stress, and τ cr = γ cr PI 0.2 is the critical shear stress in which PI is the plasticity index and γ cr = 0.7 Pa for 0.35 < τ cr < 1.4 Pa.
To account for sediment deposition, the sedimentation velocity is calculated v s = (c b w s cos α s )/(1 − n), in which c b is the near-bed volumetric sediment concentration and is estimated by the relation c b = r 0 C in which r 0 = 1 + 31.5(u * /w s ) −1.46 is a dimensionless parameter. The settling velocity for a single particle w o is calculated based on Equation (10), and hindered velocity w s is calculated based Equation (12) as follows: Note that more advanced models can be used for the numerical computations of turbidity currents. Nonetheless, we used a simple model in our analysis, as our objective was to study the relationship between the initial conditions at the discharge source and the run-out distance of the mining-generated turbidity current, rather than obtaining precise results. Figure 10 shows the case considered in the numerical simulations. The origin of the z coordinate was at the center of the discharge source. The streamwise coordinate in the downstream region is denoted as s.

Model Application
A typical forward velocity of a PNMT ranges between 0.25 m/s and 0.5 m/s. To reduce the environmental effect of the discharge, engineers aim to have a discharge velocity equal to the forward velocity of the PNMT. In our analysis, therefore, we considered a stationary PNMT with two discharge velocities of 0.25 m/s and 0.5 m/s and three different scenarios for the volumetric discharge concentration: 1.5%, 2%, and 2.5%. The corresponding diameters of the discharge, which result in identical volumetric suspended sediment transport rates per unit width UCH, were 0.45 m, 0.34 m, and 0.27 m, respectively.
At the impingement point, the flow makes a sharp turn, and a turbidity current is formed, which travels downstream and interacts with the bed surface. Consequently, depending on the conditions, net erosion or deposition upon the bed could take place. For simplification, we assumed that we have a complete transfer of the plume characteristics to the downstream region, which has a constant bed slope angle of 3 • , in agreement with the average slope found in the eastern part of the Clarion-Clipperton Zone (CCZ) [16]. Table 5 summarizes the sediment properties used in the numerical computations. Most of these properties were obtained from the technical report of Global Sea Mineral Resources NV on the field experiments conducted in the eastern part of the CCZ [16]. Note that the non-Newtonian aspect was not considered in our numerical analysis, but it might be relevant for real cases.

Comparison of the Results
For simplicity, we assumed that the run-out distance of the turbidity current was the main criterion to assess its environmental harmfulness. Additionally, a turbidity current of a lower volumetric suspended sediment at the end of the considered numerical domain tends to die out earlier and thus has a shorter run-out distance. As this paper deals with the nearfield flows, we analyzed the results up to 100 m downstream the impingement point.
Various discharge scenarios were explored and compared, as illustrated in Table 6. The values of U, C, h, and UCh (volumetric suspended sediment transport rate per unit width) at s = 100 m are documented. The evolution of the plume and the turbidity current along the downstream region is shown in Figure 11 for a discharge source located at an elevation of z = 1.5 m. The cases involving a discharge velocity of 0.25 m/s and 0.5 m/s had an identical initial volumetric suspended sediment transport rate of 1.69 × 10 −3 m 2 /s and 3.38 × 10 −3 m 2 /s, respectively. Although the discharge had the same volumetric suspended sediment at the discharge source, the plume evolution resulted in different upstream boundary conditions for the turbidity current. In these cases, a discharge of a lower concentration resulted in a shorter run-out distance. The results also showed that a higher discharge velocity resulted in a longer plume trajectory and a higher initial velocity for the turbidity current, which in turn led to a longer run-out distance.
To study the effect of the elevation of the discharge source on the hydrodynamics of the turbidity current, two additional runs were conducted: a run with a discharge duct located directly on the sea bed and a run with a discharge source located at an elevation of z = 0.5 m (see Table 6). It was found that a smaller elevation resulted in a less diluted, less energetic turbidity current, and thus a shorter run-out distance.
The numerical results showed that the overall behavior of the turbidity current did not vary between the considered cases: the layer thickness increased downstream; the layeraveraged concentration decreased downstream; the turbidity current initially decelerated until it reached a constant layer-averaged velocity.
The overall comparison of the numerical results showed that the volumetric suspended sediment of the plume at the impingement point is crucial for the environmental impact of the turbidity current. In other words, a lower volumetric suspended sediment at the impingement point led to a shorter run-out distance of the turbidity current, being less harmful. Table 6. Initial conditions of the numerical runs and the numerical results; V e , φ e , d e , and V e φ e d e are the results at the impingement point, while U e , C e , h e , and U e C e h e are the results 100 m downstream the impingement point.

Initial Condition
Plume Model Turbidity Current Model   Figure 11. Comparison of the numerical results using different initial momentum fluxes; spatial evolution of a plume (left) and spatial evolution of a turbidity current (right).

Background
The ocean bed is mainly composed of cohesive sediments (or "mud"), which consist of granular, organic, and mineral solids in a liquid phase. The solids are clay, sand, and silt, and the liquid phase is water [95]. Due to the sticking properties of the cohesive sediment, which result from the presence of clay and organic materials, the particles may undergo break-up and aggregation processes; the latter is known as flocculation. This process occurs on the microscale level between sediment particles and causes them to aggregate into larger flocs. Flocculation occurs through three mechanisms: (i) Brownian motion, i.e., the random movement of the particles; (ii) Differential settling, i.e., the particles with high settling velocities collide with the particles with low settling velocities and aggregate together; (iii) Turbulent mixing.
The turbulent mixing and concentration values within the turbidity flows are the main factors governing the flocculation process, which are responsible for the cohesive particles coming into contact, resulting in large flocs [95]. The shape factor of the flocs is expected to play a major role in determining the settling velocities, as flocs tend to form irregular shapes [18].
To understand the flocculation process in dredging operations, Reference [96] conducted a field experiment to quantify the flocs' density and their volumetric fractions. Special attention was paid to the suspended flocs in order to study bed consolidation during the dredging process. It was found that the flocs of densities ρ bed f locs < 1200 kg/m 3 represent 0.5 to 0.8 of the suspended mass, while the flocs (1200 kg/m 3 > ρ bed f locs > 1800 kg/m 3 ) represent 0.2 to 0.5 of the total suspended mass. Based on these results, they concluded that the settling velocity and floc sizes increase gradually with time, and a high concentration value of the suspended sediments was found to be favorable for floc formation.
For polymetallic nodule mining applications, Reference [18] conducted a series of experiments with different flow conditions, e.g., shear rate and concentration, to study the flocculation process between sediment particles. They performed the experiments with CCFZ sediment, with different concentrations (105, 175, 500 mg L −1 ) and different shear rates (2.4, 5.7, 10.4 s −1 ). They found that a concentration of 500 mg L −1 and a shear rate of 2.4 s −1 led to a high flocculation efficiency, inducing a high settling flux. Furthermore, they calculated the settling velocities of the particle sizes of 70-1357 µm (floc), which were in the range of 7-355 md −1 . Based on these results, they concluded that the particles would deposit rapidly with the typical deep-sea flow conditions. In addition, in the nearfield region, a portion of the sediments settles, creating the so-called "blanketing", which is defined as a thin layer of deposited sediments over the ocean bed. Using the experimental observations, Reference [18] estimated that the blanketing effect would reach 9 km far away from the mining site.
Laboratory experiments were conducted to test the flocculation effect on the settling velocity of the particles and crust debris using in situ samples from sea water obtained near a seamount surface with concentrations of 20 mg/L and 100 mg/L [17]. It was found that flocculation takes place in both crust debris and sediment particles. Nonetheless, the measurements showed that the debris (1700 kg/m 3 ) has a lower density than the sediment (2600 kg/m 3 ). Moreover, the mean settling velocities were reported as 11 mm/s and 9 mm/s for the sediment and the debris, respectively. It was noticed that the mean settling velocities obtained by [18] were lower than [17] due to the existence of extracellular polymers and bacteria in the real sea water and the electrostatic properties of the crust particles.
The above-mentioned observations indicated that flocculation could be a key phenomenon in enhancing the particles' settling and reducing the plume concentration. The sediment concentration and turbulence levels are much higher in the nearfield region compared to the farfield region. As a result, optimizing the discharge concentration and shear rate is highly needed to enhance the flocculation potential. Furthermore, an assessment of the deposited layer parameters (e.g., porosity, consolidation, structure) is required for the environmental impact assessment. It is expected that the discharged flow would experience non-Newtonian behavior in the case of high sediment concentrations. The threshold of the transition between the Newtonian and non-Newtonian behavior of the discharged flow remains an open question.

Numerical Assessment of the Flocculation Effect
To explore the flocculation effect on the mixture discharged from a PNMT, we utilized the numerical model of turbidity currents described in Section 5.2. Note that the flocculation effect on the first region of the discharge (jet/plume) was neglected here due to the short residence time of particles. We took into account the increase of the particle size and settling velocity along the run-out distance of the turbidity current as a result of flocculation. Moreover, since the residence time plays a major role in forming flocs, we extended the numerical domain to 350 m in the s-direction to capture the flocculation effect.
For the sake of comparison, we carried out two additional runs using the initial conditions of Run 8 (Table 6): a run excluding the flocculation effect, Run 8 N.F , and another run including the flocculation effect, Run 8 F . To the best of our knowledge, Reference [18] was the only experimental study that investigated the flocculation effect using CCFZ sediments under shear rates comparable to those occurring in polymetallic mining operations. Therefore, we adopted the measurements of [18] in our run to account for flocculation; their observations showed that a particle size of 12 µm within 10-50 min and under shear rate of 2.4-10.4 s −1 could reach from 250-550 µm. Based on the results of Run 8 N.F , we calculated the traveling time needed for the turbidity current to reach s = 350 m, and it was nearly 25 min. Additionally, we calculated the shear rate, and it was comparable to the values of [18] mentioned above. Within 25 min, following the measurements of [18], the suspended particle size would gradually increase from 12 µm to 550 µm. Run 8 F was performed using the same initial condition of Run 8 N.F with the difference that the particle size was increased every step and the settling velocity was updated accordingly following the formula given by [18] for a 10.4 s −1 shear rate and a 500 mg/L particle concentration as follows: The evolution of the turbidity current along the downstream for Run 8 F and Run 8 N.F is shown in Figure 12. The results manifestly show the flocculation effect on the hydrodynamics of the turbidity current; flocculation results in a slower and less concentrated turbidity current. This suggests that the flocculation effect could be observed in the nearfield region. In practice, this implies that optimizing the flow conditions (e.g., concentration and shear rate) in the nearfield region to enhance the flocculation process would be effective in increasing the sedimentation rate of particles in the farfield region.

Synthesis and Outlook
This study was conducted to survey the state-of-the-art knowledge and physical processes of the sediment-water mixture discharged from a PNMT. This resulted in the identification of some relevant knowledge gaps. The numerical assessment conducted within this paper revealed that the volumetric suspended sediment of a plume at the impingement point primarily indicates the extent of the environmental hazard posed by the generated turbidity current. A smaller volumetric suspended sediment would produce a shorter run-out turbidity current. This finding can be taken into account in designing the PNMT to minimize the environmental impact of the mining process.
The dispersion of the mining plumes could be reduced by the effect of flocculation. The shear rates and residence time of particles within a turbidity flow play a major role in triggering the flocculation process. Therefore, it is recommended to conduct experimental research to point out the possibilities to enhance the flocculation phenomenon in the nearfield area.
It is expected that a wake would form behind the PNMT, as it moves forward. Further research is required to address the impact of this wake, as it might be a trigger mechanism to increase the turbulent shear rates within the turbidity flow, thereby increasing the flocculation probability. On the other hand, the wake is turbulent and promotes water entrainment and the dilution of the discharged mixture. Properly scaled lab experiments and validated numerical models will be beneficial to quantify this effect.
It is of great importance to study the interaction between the discharged mixture and the sediment bed, in particular at the impingement region, in order to define the upstream boundary conditions of the mining-generated turbidity current. An integrated research approach combining soil and fluid mechanics is required to develop an in-depth understanding of this interaction.

Conflicts of Interest:
The authors declare no conflict of interest.