Next Article in Journal
From Plates to Baits: Using a Remote Video Foraging System to Study the Impact of Foraging on Fouling Non-Indigenous Species
Next Article in Special Issue
Environmental Impacts and Challenges Associated with Oil Spills on Shorelines
Previous Article in Journal
Maximum Power Control Algorithm for Power Take-Off System Based on Hydraulic System for Floating Wave Energy Converters
Previous Article in Special Issue
Recent Advances in Functional Materials for Wastewater Treatment: From Materials to Technological Innovations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

An Overview of Oil-Mineral-Aggregate Formation, Settling, and Transport Processes in Marine Oil Spill Models

1
Department of Civil and Resource Engineering, Faculty of Engineering, Dalhousie University, Halifax, NS B3H 4R2, Canada
2
Department of Engineering, Faculty of Agriculture, Dalhousie University, Truro, NS B2N 5E3, Canada
3
School of Marine Science, Sun Yat-Sen University (Zhuhai Campus), Zhuhai 519082, China
4
Pearl River Estuary Marine Ecosystem Research Station, Ministry of Education, Zhuhai 519082, China
5
Ocean and Ecosystem Science Division, Fisheries and Ocean Canada, Bedford Institute of Oceanography, Dartmouth, NS B2Y 4A2, Canada
*
Authors to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(5), 610; https://doi.org/10.3390/jmse10050610
Submission received: 31 January 2022 / Revised: 19 April 2022 / Accepted: 21 April 2022 / Published: 30 April 2022
(This article belongs to the Special Issue Reviews in Marine Environmental Science and Engineering)

Abstract

:
An oil spill is considered one of the most serious polluting disasters for a marine environment. When oil is spilled into a marine environment, it is dispersed into the water column as oil droplets which often interact with suspended particles to form oil-mineral-aggregate (OMA). Knowing how OMA form, settle, and are transported is critical to oil spill modelling which can determine the fate and mass balance of the spilled volumes. This review introduces oil weathering and movement, and the commonly used numerical models that oil spill specialists use to determine how a spill will evolve. We conduct in-depth reviews of the environmental factors that influence how OMA form and their settling velocity, and we review how OMA formation and transport are modelled. We point out the existing gaps in current knowledge and the challenges of studying OMA. Such challenges include having to systematically conduct laboratory experiments to investigate how the environment affects OMA formation and settling velocities, and the need for a comprehensive algorithm that can estimate an OMA settling velocity.

1. Introduction

Although it is critically important to reduce the consumption of petroleum hydrocarbons, oil products are still important to society due to the lack of alternative energy sources such as fuel for marine transportation. Producing and transporting oil products means that oil spills can potentially occur due to accidents and/or operational errors [1]. Spills of crude oils and other petroleum products in a marine environment is a major environmental problem. A report by the International Tanker Owners Pollution Federation (ITOPF) stated that there were about 10,000 large oil spills (>7 tonnes) from tankers and oil carriers from 1970 to 2018 [2]. Although the frequency of large spills has decreased over the past five decades [2], the number of small spills (for example, the Marathassa oil spill in English Bay, British Columbia in 2008) is still high [3]. These small oil spills usually occur in the open ocean and are caused by inappropriate onboard operations. Furthermore, because of their relatively small environmental impact, the small spills receive less public concern than large oil spills [4]. However, regardless of their size, oil spills are some of the most critical forms of marine pollution, leading to significant long-term impacts on the environment and socio-economy of affected areas. Therefore, local authorities and oil producers develop contingency plans to prevent, control, and reduce the negative impacts of oil spills. Accurately predicting oil spill behavior is critical for an immediate, quick, and efficient spill response.
Oil spill modelling involves simulating oil spill trajectories and behavior based on mathematical algorithms controlled by hydrodynamics and wind forcing. Once oil is spilled into water, two common processes occur simultaneously: weathering, which alters the physico-chemical properties of the spilled oil; and transport, which increases the size of the contamination area. Three factors affect these two processes: the volume of oil spilled, where the oil is spilled, and the type of oil [5]. Numerical modelling is commonly employed to predict the fate and behavior of spilled oil, and modelling is an important aspect of oil spill contingency planning and coastal management. For example, French-McCay et al. [6] used the Spill Impact Model Application Package (SIMAP) to simulate the trajectories and fate of the Deepwater Horizon oil spill. Over the past few decades, other numerical models have been developed to simulate oil spills in marine environments including: OSCAR (Oil Spill Contingency and Response) model by SINTEF (Stiftelsen for INdustriell of TEknisk Forskning ved NTH—Foundation for Industrial and Technical Research) [7,8], OILMAP and SIMAP by RP-ASA (Applied Science Associates, Inc., Narragansett, RI, USA) [9,10], the MOHID water model by MARETEC (Marine and Environmental Technology Research Center) [11], GNOME (General NOAA Operational Modeling Environment) by NOAA (National Oceanic and Atmospheric Administration) [12], and OpenDrift model and its subclass, OpenOil, by the Norwegian Meteorological Institute [13,14].
When oil is spilled into a marine environment, a significant portion of the oil is dispersed into the water column as oil droplets [15,16]. These oil droplets interact with suspended sediments to form oil-mineral-aggregate (OMA) as shown in Figure 1. OMA have been observed to form in many oil spills. For example, in 2010, more than 3.2 million liters of diluted bitumen spilled from the Enbridge Line 6B pipeline into the Kalamazoo River in Michigan, and a large amount of spilled oil (heavy component) interacted with the sediments to form OMA along channel margins, backwaters, and oxbows [17,18]. The formation, settling, and transport of the OMA are key to the transport and fate of the spilled oil [19,20,21].
Although several reviews on oil spill modelling and studies about OMA have been published in recent years, an updated and comprehensive review for OMA formation, settling, and transport is lacking. Our paper, therefore, aims to provide a dedicated review of oil spill trajectory modelling, OMA formation and settling processes, and modelling of OMA processes. This review is organized as follows: (1) Section 2 describes oil weathering and movement, and commonly used oil spill models to give an overview of current oil spill modelling; (2) Section 3 details how OMA form and settle, and we describe how OMA formation and transport have been and are currently modelled; (3) Section 4 discusses existing challenges and future work on modelling OMA formation, settling, and transport; (4) in Section 5, we summarize our review with key findings. Figure 2 illustrates the structure and themes of this review.

2. Oil Spill Weathering and Movement and How They Are Modelled

2.1. Oil Weathering

Oil weathering refers to oil spreading, evaporation, emulsification, and natural dispersion, followed by dissolution, photo-oxidation, biodegradation, and sedimentation [23]. Each process happens at different times at varying rates. For example, oil evaporation begins immediately when the oil slick first presents on the water surface. Oil evaporation is a relatively short process. In contrast, oil emulsification is a long process that begins relatively slowly [24]. Furthermore, different weathering processes have different impacts on the physico-chemical properties of the spilled oil. Among these processes, photo-oxidation and biodegradation are the most critical for modifying the chemical compositions of spilled oil [23,25]. Furthermore, as evaporation decreases the volume of the spilled oil, it also increases the oil’s density [26]. After evaporation, the remaining heavy components of a lighter oil can combine with sediments in the water column and subsequently sink to the bottom. The oil remaining at the surface naturally disperses into the water column [27]. Lastly, oil emulsification depends on the types of oil and conditions in the surrounding environment. Emulsification changes the oil’s physical properties such as water content and viscosity [28].

2.2. The Movement of Spilled Oil

Knowing how spilled oil moves on the water surface is important to oil spill modelling. Spreading is one of the most critical movement processes because spreading oil can increase the size of the contaminated area [24]. Furthermore, gravity and oil-water interfacial tension can help the spilled oil spread into a slick over the water surface [23]. The effects of gravity could gradually diminish over time, but the oil would continue to spread under the effects of interfacial tension [29]. The spreading rate of spilled oil is closely related to environmental conditions, such as tidal streams, water currents, coastlines, and wind speeds, as well as being related to oil properties [24,30,31,32]. For instance, in a confined waterway, it is difficult for spilled oil to spread quickly to reach the shoreline [33], but, in the open ocean, spilled oil can spread rapidly and widely to cover a large area in a short period [24]. Moreover, spilled oil that has lower molecular weight components presumably spreads more quickly [24]. Additionally, spreading depends on oil weathering. For example, evaporation and dissolution can reduce the area covered by the spreading oil [32].
Oil advection at the sea surface and into the water column is another crucial movement process and mainly depends on ocean currents, wind, and waves. The effects of the wind and currents on spilled oil depend on the spill location and hydrographic conditions.

2.3. Commonly Used Models

Various oil spill models have been developed to track the fate and behavior of spilled oil. These oil spill models are classified as oil weathering models or trajectory/deterministic models. Oil weathering models (OWM) do not consider oil transport and only simulate oil weathering processes.
The ADIOS2 (Automated Data Inquiry for Oil Spills, Silver Spring, MD 20910, USA) oil weathering model, which developed by NOAA, mainly includes the weathering processes of evaporation, dispersion, emulsification, spreading, and sedimentation. ADIOS can also assess how effective clean-up techniques would be [34,35].
The IKU Oil Weathering Model, which developed by SINTEF, is also commonly used. The IKU OWM is also known as SINTEF OWM. SINTEF OWM is based on small and mesoscale laboratory tests and full-scale field experiments [7,8], and SINTEF OWM considers emulsification, natural dispersion, and evaporation as its main weathering processes. The model applies a pseudo-component method to calculate how much oil is lost to evaporation, and Mackay’s equation evaluates the viscosity of the emulsified oil [36,37]. It is also worth mentioning that nearly 200 different types of oil have been characterized by the SINTEF laboratory in Norway, all of which are included in SINTEF OWM [38].
More specialized weathering models, such as the Diluted Bitumen Weathering Model (DBWM), have been developed and recently validated against meso-scale experimental data to simulate the weathering of unconventional crude products such as diluted bitumen [39].
Trajectory/deterministic models are grouped as two-dimensional (2D) and three-dimensional (3D) models [40]. A 2D model only simulates the trajectory of spilled oil at the water’s surface, whereas a 3D model tracks the spilled oil at the surface and in the water column. Various algorithms have been developed for trajectory/deterministic models, including Eulerian methods, Lagrangian methods, and Eulerian–Lagrangian methods (ELMs). According to state-to-the-art reviews [40,41,42], Lagrangian methods are the most popular algorithms for 3D oil spill models. For instance, the COSMoS model developed by Environment and Climate Change Canada [43] and SPILLCALC developed by Tetra Tech EBA [44] use Lagrangian algorithms. Table 1 lists examples of 3D models that use Lagrangian algorithms.
Oil spill models can also be classified based on their purpose. After the Deepwater Horizon oil spill (2010 to 2020), Ainsworth et al. [63] categorized and summarized 330 published modelling applications developed by the Gulf of Mexico Research Initiative (GOMRI), the Natural Resource Damage Assessment (NRDA), and other researchers. Ainsworth et al. classified the models by (1) circulation/mixing, (2) abiotic transport (far field), (3) oil fate, (4) biotic transport, (5) biological impacts, (6) other plume dynamics, (7) turbulence/local mixing, (8) water chemistry, (9) atmosphere, (10) oil spill response support, and (11) other [63]. Furthermore, Keramea et al. [64] summarized commonly used models based on model capabilities and purpose. Keramea et al. classified the oil spill models into: (1) surface oil spill models and blowout/buoyant plume models, (2) operational response models, (3) deep sea blowout/buoyant plume models, and (4) models with spill response and environmental impact assessment support [64]. Zhao et al. [65] evaluated algorithms for oil spreading, evaporation, emulsification, dispersion, dissolution, biodegradation, and photo-oxidation [65]. We note that sedimentation (OMA formation and transport) is absent in all the 2D and most 3D oil spill models. The 3D models that include sedimentation usually describe it with simple coefficients and do not have detailed sedimentation algorithms.

3. Oil-Mineral-Aggregate (OMA)

3.1. Factors Influencing OMA Formation

OMA formation is affected by many factors, including oil properties (type, droplet size, and concentration), particle properties (size and shape, organic matter content, density, and concentration), and ambient conditions such as mixing energy, salinity, and temperature [66]. It has been reported that the average aggregate size is positively correlated to the oil-particle ratio (g oil/g sediment) [67]. However, when only considering sediment concentration, the relationship between particle concentration and the average aggregate size fits a Gaussian distribution with the highest average aggregate size of ~1400 um at a particle concentration of 0.8 g/L [67]. Kehlifa et al. [68] showed that normalized cumulative size distributions of OMA are like the size distributions of mineral-free oil droplets formed in seawater.
In terms of the oil viscosity’s influence on how OMA form, Wood et al. [69]) observed that the aggregation rate of OMA is greater as the oil viscosity and density decrease. However, Khelifa et al. [70] indicated that oil viscosity has a negligible effect on the maximum and mean sizes of OMA.
The types of clay in the suspended sediment impacts how OMA form. Zhang et al. [71] and Wang et al. [72] found that the size of an OMA increases as the minerals change from hydrophilic to hydrophobic. The OMA’s size increases from a few micrometers for a natural kaolin (hydrophilic) to tens of micrometers for a modified kaolin (hydrophobic). The change in size occurs because the hydrophobicity of minerals promotes the affinity of minerals to oil, which enhances the formation of OMA. Khelifa et al. [68] also indicated that the clay type is crucial to the OMA’s size, which has even more influence on aggregate size than the type of oil does.
The effect of salinity on how OMA form is controversial. Several studies illustrate that low salinity (0–5 ppt) strongly influences OMA formation, and no significant effect exists for greater salinity increases [68,73,74,75]. For instance, Le-Floch et al. [74] reported that the proportion of oil in an OMA increases as the salinity increases, but the oil content stabilizes when salinity is greater than ~2 ppt. Khelifa et al. [68] illustrated that the median and maximum size of OMA increases significantly as the salinity increases from 0 ppt to 1.2 ppt, but a further increase in salinity to 3.5 ppt decreases the OMA’s size [68]. However, Payne [76] found the opposite result: an OMA forms more quickly with salinities of 15 and 30 ppt than in freshwater. Furthermore, Guyomarch et al. [77] reported that the OMA formation rate slows as salinity increases from 10 to 50 ppt. Theoretically, the effects of salinity on OMA formation occur because of the change in the electrical double layer thickness of the oil and mineral particles. As salinity increases, the layer thickness compresses, which reduces surface repulsion and increases the inter-particle attraction between the oil and mineral. However, the influence of salinity on OMA formation is minimal once the layer thickness is below a threshold as the level of salinity in the water rises. Therefore, OMA formation eventually peaks and stabilizes as the salinity increases [74].
Mixing energy is an important environmental factor that influences the size of OMA [58,71,78,79,80]. An adequate mixing energy promotes the interaction between dispersed oil droplets and suspended particles so that OMA can form [79]. On the other hand, an OMA could break apart if the mixing energy is too aggressive [71]. Ji et al. [58] conducted a series of laboratory experiments to investigate the impact of mixing energy on the time it takes for OMA to form. They suggested that the collision efficiency between oil and sediments under a higher mixing energy rate (200 rpm) is three times higher than under a lower mixing energy rate (150 rpm). High mixing energy results in rapidly forming OMA during the first 3 h, followed by the OMA breaking apart. Ji et al. found that, after 24 h of mixing, only in low energy mixing cases would the dispersed OMA reform because there would be more oil droplets on the water’s surface.
As a human initiated influence, chemical dispersants are key to OMA formation. Page et al. [15] carried out a series of wave tank experiments to study the effectiveness of using dispersants as an oil spill chemical countermeasure in the surf zone. Their experiment results indicated that chemically dispersed oil associates less intensively with mineral matter than physically dispersed oil [15]. Zhang et al. [71] found that applying chemical dispersants is a dominant factor in how OMA form and behave. Several wave tank experiments conducted by Lee et al. [81] indicated that chemical dispersants reduce surface tension between the oil and water and stimulate interaction between oil and fine mineral, eventually increasing the concentration of OMA. Lee et al. also indicated that the chemical dispersants reduce the size of OMA to a mean diameter of 15–25 µm [81]. However, Khelifa et al. [82] observed that chemical dispersants enhance the stickiness of oil to lead to larger OMA, probably due to the surfactant coating the oil droplets. Fu et al. [83] carried out roller table experiments to explore how marine oil snow forms when a stereotype oil dispersant (Corexit EC9500A) is present. Fu et al. found that adding the chemical dispersant enhances particle aggregation and formation of marine snow. O’Laughlin et al. [84] claimed that dispersants delay flocculation of natural sediments and create a surplus of available sediment to interact with oil droplets and potentially form OMA.

3.2. OMA Settling

The settling velocity of OMA mainly depends on their composition densities. The buoyancy of OMA should cause OMA to rise to the surface if the composition is mainly oil; otherwise, the OMA composed primarily of sediment should sink to the benthic layer [85]. A better understanding of OMA settling is critical to predicting the fate and transport of OMA.
Table 2 summarizes laboratory experiments of OMA settling. Muschenheim and Lee [86] investigated the amount and rates of OMA formation and settling velocities by using a focused flow reactor. They reported that large (100 to 200 µm) OMA have settling velocities ranging from 2.2 to 10.4 mm/s. Waterman and Garcia [87] tested settling velocities by using a 1.6 m height Plexiglass settling column and observed that the OMA had settling velocities between 1.0 and 11.2 mm/s, with the majority being between 1.0 and 3.0 mm/s. O’Laughlin et al. [84] conducted experiments to measure the variability in particle size and settling velocity of OMA in response to sediment concentrations and presence/absence of chemical dispersants. O’Laughlin et al. found that OMA size, settling velocity, and effective particle density increase under a higher concentration of suspended sediment (comparing 10 mg/L with 50 mg/L), which indicates that high concentration suspension produces large, inorganic-derived flocs that settle rapidly. Their results also suggested that dispersants inhibit natural sediment flocculation. Ye et al. [88] conducted laboratory experiments to investigate the influence of clay types on OMA structures and settling velocities by using the LabSFLOC-2 system and digital microscopy. Ye et al. found that for low stickiness Kaolinite clay, droplet OMA form with much smaller settling velocities than the settling velocities of pure Kaolinite flocs. Furthermore, Ye et al. found that high stickiness Bentonite clay, generates flack/solid OMA with settling velocities greater (by 25%) than settling velocities of pure Bentonite flocs.

3.3. Modelling OMA Size Distribution

The size of OMA varies from a few micrometers to hundreds, even thousands, of micrometers [67,71,78,79]. The range of sizes of OMA illustrates the complexity of OMA formation. This complexity encouraged scientists to develop sophisticated numerical models to predict how OMA form and subsequently transported. However, modelling OMA formation and transport is still at an early stage and the ability of the models to accurately predict either process is limited. In most oil spill models, only the ratio of mineral-stabilized oil is computed from standard equilibrium partitioning theory (by using a dimensionless coefficient), and the concentration of suspended particle matter in the water column is just estimated [42]. Therefore, there is still a need to develop even more advanced OMA models that can be incorporated into existing oil spill models.
With regards to modelling OMA formation, the population balance equation, based on collision theory between oil droplets and suspended sediment materials, is currently widely used [68,89,90,91,92,93]. The population balance equation for particle collision due to Brownian motion was first proposed by Smoluchowski [94] and is listed as Equation (1) in Table 3. After Smoluchowski developed their equation, other physical mechanisms, such as shear turbulence and differential settling, were included in the modelling (Equation (2) in Table 3) [95]. Sterling et al. [92] extended the population balance equation to simulate the changes in particle size distribution and density due to aggregation (Equations (3) and (4) in Table 3). Sterling et al. conducted batch flocculation experiments and introduced an algorithm to estimate homogeneous collision efficiency values ( α H O M O ) to account for the effects of particle types on OMA formation. Furthermore, Sterling et al. found that α H O M O is greater for clay (0.7) and crude oil (0.3) compared to silica (0.01); thus, they classified clay and crude oil as cohesive particles and classified silica as a non-cohesive particle. In addition, they found that α H O M O is similar for oil-clay (0.4) and oil-silica (0.3), suggesting that crude oil increases the aggregation of non-cohesive silica particles [92]. However, the Sterling et al. model utilized a looped term to represent oil–oil, oil–particle, and particle–particle interactions, mainly because the laboratory operations probably could not distinguish between the interactions. Moreover, Sterling et al. did not include an OMA breakup process in their model.
Just like Sterling et al.’s model, breakup is neglected in several other OMA formation models [89,93,102,103]. For instance, Bandara et al. [89] developed a numerical model by using a three-dimensional (3D) advection-diffusion equation to simulate oil-sediment interaction (Equation (5) in Table 3). They used the Lagrangian Parcel (LP) method to reduce their program size and to operate more efficiently. Bandara et al.’s simulations showed that up to 65% of released oil can be removed from the water column as OMA. When oil droplets and sediment particles are smaller than 0.1 mm, more OMA form. Bandara et al. also stated that their lacking knowledge of oil sediment aggregation collision efficiencies, sediment aggregation efficiencies, and oil partitioning led to some uncertainty in their results, and more laboratory and field experimental work/data would further test and improve their model’s adequacy [89]. Zhao et al. [93] developed the A-DROP model, based on the population balance equation, by introducing a new formulation of oil-mineral coagulation efficiency (Equation (6) in Table 3) to account for the coasted areas on the surface of oil droplets, the effects of hydrophobicity, and the ratio of particle to droplet size. The Zhao et al. formulation satisfactorily estimates oil trapping efficiency by incorporating the effects of shape and packing on OMA coagulation. They suggested that increasing particle concentration in the swash zone quickens oil–particle interaction, but the amount of oil trapped in the OMA does not correspond to the increasing particle concentration [93]. Recently, Cui et al. [102] developed an OMA formation model, OPAMOD, within the Coupled Ocean-Atmosphere-Wave-and-Sediment Transport (COAWST) modelling system, by modifying the existing population balance equation (Figure 3). The authors performed sensitivity tests on fractal dimension and collision efficiency by using the OPAMOD model. They stated that fractal dimension is important to OMA size distribution because fractal dimension influences the effective particle density; however, collision efficiency has less impact on the size distribution [102].
Although the size distribution of OMA is successfully simulated by the models that do not include disaggregation, the OMA breakup process is still valuable when modelling OMA formation. Khelifa et al. [90,91,97] developed a Monte Carlo simulation model (involving disaggregation) to predict OMA size distribution based on a population balance equation. The Monte Carlo method (Equation (7) in Table 3) is applied as a probabilistic tool to solve the model. The specific event, aggregation or disaggregation, is selected randomly at each step during the simulation. The simulation stops once equilibrium is reached. Khelifa et al. suggested that it is appropriate to integrate their new breakage model to describe oil droplet formation. The maximum permissible size of the oil droplets in their simulation is the key physical input of the model and can be upgraded by integrating an empirical-theoretical model [97]. We note that this model predicts the particle size distribution well, but information about oil (mass) in OMA, needed for further risk assessment, is unavailable. In addition, the conceptual time to run the Khelifa et al. model using the Monte Carlo method is unknown.
In addition to the size distribution of the OMA, other parameters important to OMA formation, such as the time scale of formation [98,104], oil trapping efficiency [74,99,100], and oil sinking rate [101], have been studied and predicted. Hill et al. [98] presented Equation (8) (Table 3) to define the time required to form stable OMA [98]. They stated that OMA formation can occur over a short time scale. Among all their simulations, about 7% of the required times were shorter than 30 s, and 10% were longer than a day. In their modelling, when sediment concentration was high and droplets were large, the OMA formation time was the shortest. Hill et al. also suggested that their equation is suitable for the coastal zone [98]. Ajijolaiya et al. [99] developed Equation (9) (Table 3) to describe oil trapping efficiency. They found that trapping efficiency increases with improving sediment concentration and decreasing sediment size. Ajijolaiya et al. also stated that maximum trapping efficiency occurs when the range of sediment concentrations is near unity. Wang et al. [100] conducted several experiments and developed an empirical equation to calculate oil trapping efficiency as a function of sediment concentrations (Equation (10) in Table 3). Their simulation results indicated that the formation rate of OMA depends on sediment concentration, mixing time, salinity, and the use of a chemical dispersant in the water column. Wu et al. [101] developed a simple approach to estimate the sinking rate of spilled diluted bitumen (Equation (11) in Table 3). They suggested that the sinking rate is impacted by the density of oil, sediment, water, and oil size distribution.

3.4. OMA Transport Modelling

The transport of OMA is critically important to oil spill modelling. Some models have already been developed to simulate OMA transport [85,89,92,96,103,105,106,107,108]. The fate and transport of OMA have been typically modeled by implementing the advection-diffusion equation and a random walk model. For example, as we noted before, Bandara et al. [89] used the Lagrangian Parcel (LP) method to develop a model based on the 3D advection-diffusion equation to predict the fate and transport of OMA. The model simulates six processes: advection and diffusion of oil and sediments, dissolution of oil, aggregation of sediments, aggregation of oil-sediment particles, oil partitioning, and deposition of sediment and OMA. Unfortunately, Bandara et al. made some assumptions due to their lacking knowledge of oil sediment aggregation efficiencies, sediment aggregation efficiencies, and oil partitioning.
Niu and Lee [106] used a 3D random walk model to simulate the transport of OMA under hydrodynamic conditions involving wave-induced velocities, random velocities due to turbulence, and a settling velocity due to gravity (Figure 4). In their model, they use a fixed number of particles to represent the OMA at the spill site at the beginning of a simulation. Then, the particles move on each subsequent time-step, according to Lagrangian motion, while the OMA size distribution is artificially defined. Jones and Garcia [103] avoided having to artificially define the size distribution by combining the random walk model with the A-DROP model to simulate the fate and transport of OMA. The Jones and Garcia’s model simulates the interactions between oil droplets and sediments during their residence time in a river, which informs the user under what conditions the oil is more likely to coagulate with suspended sediment. However, Jones and Garcia neglected the re-entrainment of settled OMA.
The settling velocity of OMA is one of the most important parameters when modelling OMA transport. There are equations that calculate settling velocity [89,93,96,106,109]. Winterwerp [109] developed the most widely used settling velocity equation basing it on Stokes’ formula (Equation (12) in Table 4). Winterwerp’s equation is suitable for spherical, massive particles in Stokes’s regime ( R e 1 ). This settling velocity equation integrates a particle tracking model or an advection-diffusion equation. Integrated model’s simulation results compare favorably with laboratory observations [89,110]. Later, Zhao et al. [102] improved the Winterwerp equation by expanding the applicable flows ( 1 < R e < 10 4 and R e > 10 4 ) as shown in Equation (13) (Table 4). Zhao et al.’s predicted settling velocity is also consistent with experimental data [103,111]. Niu et al. [103] developed a simple regression equation (Equation (14) in Table 4) to simulate the settling velocity relative to OMA diameter. Niu et al. also provide equations to estimate the concentration of the settled OMA and the oil content in each particle as shown in Equation (15) in Table 4 [106,112]. However, the Niu et al. equations are only suitable for specific conditions.

4. Challenges and Recommendations

We described different factors that influence OMA formation and settling velocity, and the varied laboratory operations that different research groups conducted. Unfortunately, no one has conducted replicated experiments, and, thus, the experimental errors are unknown, making the reported data less comparable and reliable. Therefore, in the future, the experimental runs need to be replicated. Although many factors have been observed that affect OMA formation and settling, the significance of these factors remains unknown. We highly recommended designing a comprehensive statistical analysis that systematically investigates the significance of each factor on OMA formation and settling. In addition, effects of factors interacting with each other, such as dependencies, could exist. For example, the influence of temperature on OMA formation and settling could depend on the salinity. Identifying such dependencies is important for more accurate modelling of OMA formation and settling.
As for OMA formation modelling itself, oil-oil, sediment-sediment, OMA-OMA, oil-OMA, and sediment-OMA interactions are commonly neglected in current modelling. However, some of these processes could be important in predicting the OMA size distribution. In future modelling, we suggest including these interactive processes. Collision theory is widely applied to simulate OMA formation, and the collision efficiency/frequency is crucial to adequately simulating the process. However, available collision frequency/efficiency data are currently limited. More comprehensive empirical datasets are required. Moreover, although oil trapping efficiency is simulated by some of current models, the simulation results are barely validated due to insufficient experimental data. Current laboratory/modelling studies only focus on OMA or marine oil snow, but both OMA and marine oil snow likely cause oil sedimentation. Simultaneously studying OMA and marine oil snow deserves more research attention.
We also note that most OMA formation and transport models assume the oil droplets, sediment particles, and OMA are spherical. However, these particles’ shapes vary significantly, and the shapes of the particles might be key to how OMA form and settle. Attempting to develop a shape factor for different particles is, therefore, meaningful. Lastly, an empirical equation for estimating an OMA’s settling velocity is not always available, hindering the accuracy of OMA transport modelling. It is therefore valuable to combine experimental data with statistical regression to develop an adequate OMA settling velocity empirical equation.

5. Conclusions

Accurately tracking the transport and behavior of spilled oil is important to local authorities and oil companies because the knowledge makes the oil spill response more efficient and timelier, which reduces negative impacts. Numerous efforts have been made in oil spill modelling, in which oil-mineral-aggregate (OMA) formation, settling, and transport are modelled as crucial processes. Our review summarized oil spill weathering processes, movement, and commonly used oil spill models. We further focused on current modelling of OMA processes. Briefly, we find:
  • Many environmental factors have been reported to influence OMA size distribution and settling velocity, such as temperature, salinity, sediment concentration, the presence of dispersants, and so on.
  • Statistical design and analysis should be used to determine the significance of each factor and their inter-dependencies.
  • Attempts have been made to measure settling velocities in laboratory experiments, and the reported settling velocities ranged from 1 to 10.4 mm/s depending on experimental conditions. However, the lack of an adequate empirical equation for estimating an OMA settling velocity hinders OMA transport and oil spill modelling.
  • Efforts at modelling the OMA size distribution have been made based upon collision theory. The Monte Carlo method has also been applied to model the size distribution of OMA; however, including the OMA breakup process or disregarding it could influence modelling results.

Author Contributions

Conceptualization, X.Z., Y.W. and H.N.; resources, L.L.; writing—original draft preparation, X.Z.; writing—review and editing, P.L. and H.N.; supervision, H.N.; project administration, H.N.; funding acquisition, H.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Marine Environment Observation Prediction and Response Work (MEOPAR) [Project: 2-02-03-037.3] and Multi-Partner Research Initiative (MPRI) [Project: 6.02].

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge the Marine Environment Observation Prediction and Response Work (MEOPAR) and Multi-Partner Research Initiative (MPRI).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mishra, A.K.; Kumar, G.S. Weathering of Oil Spill: Modeling and Analysis. Aquat. Procedia 2015, 4, 435–442. [Google Scholar] [CrossRef]
  2. ITOPF. ITOPF Handbook. 2019. Available online: https://www.itopf.org/knowledge-resources/documents-guides/document/itopf-handbook/ (accessed on 27 November 2019).
  3. Fingas, M. Oil Spill Science and Technology; Gulf Professional Publishing: Houston, TX, USA, 2016. [Google Scholar]
  4. Ivanov, A.Y.; Zatyagalova, V.V. A GIS Approach to Mapping Oil Spills in a Marine Environment. Int. J. Remote Sens. 2008, 29, 6297–6313. [Google Scholar] [CrossRef]
  5. Verma, P.; Wate, S.R.; Devotta, S. Simulation of Impact of Oil Spill in the Ocean—A Case Study of Arabian Gulf. Environ. Monit. Assess. 2008, 146, 191–201. [Google Scholar] [CrossRef] [PubMed]
  6. French-McCay, D.P.; Spaulding, M.L.; Crowley, D.; Mendelsohn, D.; Fontenault, J.; Horn, M. Validation of Oil Trajectory and Fate Modeling of the Deepwater Horizon Oil Spill. Front. Mar. Sci. 2021, 8, 136. [Google Scholar] [CrossRef]
  7. Daling, P.S.; Aamo, O.M.; Lewis, A.; Strøm-Kristiansen, T. Sintef/Iku Oil-Weathering Model: Predicting Oils’ Properties at Sea. Int. Oil Spill Conf. Proc. 1997, 1997, 297–307. [Google Scholar] [CrossRef]
  8. Reed, M.; Daling, P.; Lewis, A.; Ditlevsen, M.K.; Brørs, B.; Clark, J.; Aurand, D. Modelling of Dispersant Application to Oil Spills in Shallow Coastal Waters. Environ. Model. Softw. 2004, 19, 681–690. [Google Scholar] [CrossRef]
  9. RPS-ASA. Software—OILMAP. Available online: http://www.asascience.com/software/oilmap/ (accessed on 10 September 2018).
  10. RPS-ASA. Software—SIMAPTM. Available online: http://asascience.com/software/simap/ (accessed on 10 September 2018).
  11. Fernandes, R. Risk Management of Coastal Pollution from Oil Spills Supported by Operational Numerical Modelling. Ph.D. Thesis, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal, 2018. [Google Scholar]
  12. National Oceanic and Atmospheric Administration (NOAA). GNOME User’s Manual. Available online: http://www.webcitation.org/72P8XJ1SZ (accessed on 10 September 2018).
  13. Dagestad, K.-F.; Röhrs, J.; Breivik, Ø.; Ådlandsvik, B. OpenDrift v1.0: A Generic Framework for Trajectory Modelling. Geosci. Model Dev. 2018, 11, 1405–1420. [Google Scholar] [CrossRef] [Green Version]
  14. Röhrs, J.; Dagestad, K.-F.; Asbjørnsen, H.; Nordam, T.; Skancke, J.; Jones, C.E.; Brekke, C. The Effect of Vertical Mixing on the Horizontal Drift of Oil Spills. Ocean Sci. 2018, 14, 1581–1601. [Google Scholar] [CrossRef] [Green Version]
  15. Page, C.A.; Bonner, J.S.; Sumner, P.L.; McDonald, T.J.; Autenrieth, R.L.; Fuller, C.B. Behavior of a Chemically-Dispersed Oil and a Whole Oil on a near-Shore Environment. Water Res. 2000, 34, 2507–2516. [Google Scholar] [CrossRef]
  16. Sterling, M.C.; Bonner, J.S.; Page, C.A.; Fuller, C.B.; Ernest, A.N.S.; Autenrieth, R.L. Partitioning of Crude Oil Polycyclic Aromatic Hydrocarbons in Aquatic Systems. Environ. Sci. Technol. 2003, 37, 4429–4434. [Google Scholar] [CrossRef]
  17. Dollhopf, R.H.; Fitzpatrick, F.A.; Kimble1, J.W.; Capone, D.M.; Graan, T.P.; Zelt, R.B.; Johnson, R. Response to Heavy, Non-Floating Oil Spilled in a Great Lakes River Environment: A Multiple-Lines-Of-Evidence Approach for Submerged Oil Assessment and Recovery. Int. Oil Spill Conf. Proc. 2014, 2014, 434–448. [Google Scholar] [CrossRef]
  18. Lee, K.; Bugden, J.; Cobanli, S.; King, T.; McIntyre, C.; Robinson, B.; Ryan, S.; Wohlgeschaffen, G. UV-Epifluorescence Microscopy Analysis of Sediments Recovered from the Kalamazoo River; Centre for Offshore Oil, Gas and Energy Research (COOGER), Fisheries and Oceans Canada: Dartmouth, NS, Canada, 2012. [Google Scholar]
  19. Gordon, D.C., Jr.; Keizer, P.D.; Prouse, N.J. Laboratory Studies of the Accommodation of Some Crude and Residual Fuel Oils in Sea Water. J. Fish. Res. Board Can. 1973, 30, 1611–1618. [Google Scholar] [CrossRef]
  20. Johansson, S.; Larsson, U.; Boehm, P. The Tsesis Oil Spill Impact on the Pelagic Ecosystem. Mar. Pollut. Bull. 1980, 11, 284–293. [Google Scholar] [CrossRef]
  21. Muschenheim, D.K.; Lee, K. Removal of Oil from the Sea Surface through Particulate Interactions: Review and Prospectus. Spill Sci. Technol. Bull. 2002, 8, 9–18. [Google Scholar] [CrossRef]
  22. Gustitus, S.A.; Clement, T.P. Formation, Fate, and Impacts of Microscopic and Macroscopic Oil-Sediment Residues in Nearshore Marine Environments: A Critical Review. Rev. Geophys. 2017, 55, 1130–1157. [Google Scholar] [CrossRef]
  23. Fingas, M. Chapter 8—Introduction to Spill Modeling. In Oil Spill Science and Technology; Fingas, M., Ed.; Gulf Professional Publishing: Boston, MA, USA, 2011; pp. 187–200. [Google Scholar]
  24. Lee, K.; Boufadel, M.; Chen, B.; Foght, J.; Hodson, P.; Swanson, S.; Venosa, A. Expert Panel Report on the Behaviour and Environmental Impacts of Crude Oil Released into Aqueous Environments; Royal Society of Canada: Ottawa, ON, Canada, 2015. [Google Scholar]
  25. Dutta, T.K.; Harayama, S. Fate of Crude Oil by the Combination of Photooxidation and Biodegradation. Environ. Sci. Technol. 2000, 34, 1500–1505. [Google Scholar] [CrossRef]
  26. Fingas, M. The Evaporation of Oil Spills. In Arctic and Marine Oilspill Program Technical Seminar; Ministry of Supply and Services: Ottawa, ON, Canada, 1995; pp. 43–60. [Google Scholar]
  27. Walker, A.H.; Stern, C.; Scholz, D.; Nielsen, E.; Csulak, F.; Gaudiosi, R. Consensus Ecological Risk Assessment of Potential Transportation-Related Bakken and Dilbit Crude Oil Spills in the Delaware Bay Watershed, USA. J. Mar. Sci. Eng. 2016, 4, 23. [Google Scholar] [CrossRef] [Green Version]
  28. Fingas, M. Chapter 10—Models for Water-in-Oil Emulsion Formation. In Oil Spill Science and Technology; Fingas, M., Ed.; Gulf Professional Publishing: Boston, MA, USA, 2011; pp. 243–273. ISBN 978-1-85617-943-0. [Google Scholar]
  29. Fay, J.A. The Spread of Oil Slicks on a Calm Sea. In Oil on the Sea, Proceedings of a Symposium on the Scientific and Engineering Aspects of Oil Pollution of the Sea, Sponsored by Massachusetts Institute of Technology and Woods Hole Oceanographic Institution and Held at Cambridge, MA, USA, 16 May 1969; Hoult, D.P., Ed.; Ocean Technology; Springer: Boston, MA, USA, 1969; pp. 53–63. ISBN 978-1-4684-9019-0. [Google Scholar]
  30. Dew, W.A.; Hontela, A.; Rood, S.B.; Pyle, G.G. Biological Effects and Toxicity of Diluted Bitumen and Its Constituents in Freshwater Systems. J. Appl. Toxicol. 2015, 35, 1219–1227. [Google Scholar] [CrossRef]
  31. Lehr, W.J. Review of Modeling Procedures for Oil Spill Weathering Behavior. Adv. Ecol. Sci. 2001, 9, 51–90. [Google Scholar]
  32. Wang, S.-D.; Shen, Y.-M.; Guo, Y.-K.; Tang, J. Three-Dimensional Numerical Simulation for Transport of Oil Spills in Seas. Ocean Eng. 2008, 35, 503–510. [Google Scholar] [CrossRef]
  33. Michel, J.; Rutherford, N. Impacts, Recovery Rates, and Treatment Options for Spilled Oil in Marshes. Mar. Pollut. Bull. 2014, 82, 19–25. [Google Scholar] [CrossRef] [PubMed]
  34. Lehr, W.; Jones, R.; Evans, M.; Simecek-Beatty, D.; Overstreet, R. Revisions of the ADIOS Oil Spill Model. Environ. Model. Softw. 2002, 17, 189–197. [Google Scholar] [CrossRef]
  35. Lehr, W.J.; Overstreet, R.; Jones, R.; Watabayashi, G. ADIOS-Automated Data Inquiry for Oil Spills. In Proceedings of the Arctic and Marine Oilspill Program Technical Seminar, Edmonton, AB, Canada, 10–12 June 1992. [Google Scholar]
  36. Mackay, D.; Buist, I.; Mascarenhas, R.; Paterson, S. Oil Spill Processes and Models; Technical Report EE-8; Environment Canada: Ottawa, ON, Canada, 1980; p. 17. [Google Scholar]
  37. Mackay, D.; Paterson, S.; Trudel, K. A Mathematical Model of Oil Spill Behaviour; Department of Chemical and Applied Chemistry, University of Toronto: Toronto, ON, Canada, 1980; 39p. [Google Scholar]
  38. Aamo, O.M.; Reed, M.; Daling, P.S. A Laboratory-Based Weathering Model: PC Version for Coupling to Transport Models. In Proceedings of the Arctic and Marine Oilspill Program Technical Seminar, Calgary, AB, Canada, 7–9 June 1993. [Google Scholar]
  39. Li, P.; Niu, H.; Li, S.; King, T.L.; Zou, S.; Chen, X.; Lu, Z. DBWM: A Diluted Bitumen Weathering Model. Mar. Pollut. Bull. 2022, 175, 113372. [Google Scholar] [CrossRef] [PubMed]
  40. Reed, M.; Johansen, Ø.; Brandvik, P.J.; Daling, P.; Lewis, A.; Fiocco, R.; Mackay, D.; Prentki, R. Oil Spill Modeling towards the Close of the 20th Century: Overview of the State of the Art. Spill Sci. Technol. Bull. 1999, 5, 3–16. [Google Scholar] [CrossRef]
  41. American Society of Civil Engineers (ASCE). State-of-the-Art Review of Modeling Transport and Fate of Oil Spills. J. Hydraul. Eng. 1996, 122, 594–609. [Google Scholar] [CrossRef]
  42. Spaulding, M.L. State of the Art Review and Future Directions in Oil Spill Modeling. Mar. Pollut. Bull. 2017, 115, 7–19. [Google Scholar] [CrossRef]
  43. Marcotte, G.; Bourgouin, P.; Mercier, G.; Gauthier, J.P.; Pellerin, P.; Smith, G.; Brown, C.W. Canadian Oil Spill Modelling Suite: An Overview; Gulf Publishing Company: Halifax, NS, Canada, 2016; pp. 1026–1034. [Google Scholar]
  44. Tetra Tech. Spillcalc Oil and Contaminant Spill Model. Available online: http://www.tetratech.com/en/projects/spillcalc-oil-and-contaminant-spill-model (accessed on 10 September 2018).
  45. Reed, M.; Gundlach, E.; Kana, T. A Coastal Zone Oil Spill Model: Development and Sensitivity Studies. Oil Chem. Pollut. 1989, 5, 411–449. [Google Scholar] [CrossRef]
  46. NOAA. GNOME Suite for Oil Spill Modeling. Available online: https://response.restoration.noaa.gov/oil-and-chemical-spills/oil-spills/response-tools/gnome-suite-oil-spill-modeling.html (accessed on 29 January 2022).
  47. Al-Rabeh, A.H.; Lardner, R.W.; Gunay, N. Gulfspill Version 2.0: A Software Package for Oil Spills in the Arabian Gulf. Environ. Model. Softw. 2000, 15, 425–442. [Google Scholar] [CrossRef]
  48. Lardner, R.; Zodiatis, G.; Hayes, D.; Pinardi, N. Application of the MEDSLIK Oil Spill Model to the Lebanese Spill of July 2006; European Group of Experts on satellite monitoring of sea based oil pollution; European Communities: Maastricht, The Netherlands, 2006; pp. 1018–5593. [Google Scholar]
  49. Lardner, R.W.; Zodiatis, G.; Loizides, L.; Demetropoulos, A. An Operational Oil Spill Model for the Levantine Basin (Eastern Mediterranean Sea). In Proceedings of the International Symposium on Marine Pollution, Monaco, Monaco, 5–9 October 1999. [Google Scholar]
  50. De Dominicis, M.; Pinardi, N.; Zodiatis, G.; Archetti, R. MEDSLIK-II, a Lagrangian Marine Surface Oil Spill Model for Short-Term Forecasting—Part 2: Numerical Simulations and Validations. Geosci. Model Dev. 2013, 6, 1871–1888. [Google Scholar] [CrossRef] [Green Version]
  51. Danish Hydraulic Institute (DHI). DHI Oil Spill Model—Scientific Description. Available online: http://manuals.mikepoweredbydhi.help/2017/General/DHI_OilSpill_Model.pdf (accessed on 31 March 2022).
  52. Carracedo, P.; Torres-López, S.; Barreiro, M.; Montero, P.; Balseiro, C.F.; Penabad, E.; Leitao, P.C.; Pérez-Muñuzuri, V. Improvement of Pollutant Drift Forecast System Applied to the Prestige Oil Spills in Galicia Coast (NW of Spain): Development of an Operational System. Mar. Pollut. Bull. 2006, 53, 350–360. [Google Scholar] [CrossRef]
  53. Daniel, P.; Marty, F.; Josse, P.; Skandrani, C.; Benshila, R. Improvement of Drift Calculation in Mothy Operational Oil Spill Prediction System. Int. Oil Spill Conf. Proc. 2003, 2003, 1067–1072. [Google Scholar] [CrossRef] [Green Version]
  54. Spaulding, M.L.; Kolluru, V.S.; Anderson, E.; Howlett, E. Application of Three-Dimensional Oil Spill Model (WOSM/OILMAP) to Hindcast the Braer Spill. Spill Sci. Technol. Bull. 1994, 1, 23–35. [Google Scholar] [CrossRef]
  55. Garcia, M.C.; Sanz-Bobi, M.A.; del Pico, J. SIMAP: Intelligent System for Predictive Maintenance: Application to the Health Condition Monitoring of a Windturbine Gearbox. Comput. Ind. 2006, 57, 552–568. [Google Scholar] [CrossRef]
  56. Berry, A.; Dabrowski, T.; Lyons, K. The Oil Spill Model Oiltrans and Its Application to the Celtic Sea. Mar. Pollut. Bull. 2012, 64, 2489–2501. [Google Scholar] [CrossRef] [Green Version]
  57. Reed, M.; Aamo, O.M.; Daling, P.S. Quantitative Analysis of Alternate Oil Spill Response Strategies Using OSCAR. Spill Sci. Technol. Bull. 1995, 2, 67–74. [Google Scholar] [CrossRef]
  58. Ji, W.; Boufadel, M.; Zhao, L.; Robinson, B.; King, T.; An, C.; Zhang, B.; Lee, K. Formation of Oil-Particle Aggregates: Impacts of Mixing Energy and Duration. Sci. Total Environ. 2021, 795, 148781. [Google Scholar] [CrossRef]
  59. Smith, R.A.; Slack, J.R.; Wyant, T.; Lanfear, K.J. The Oilspill Risk Analysis Model of the US Geological Survey. Geological Survey Professional Paper 1227; US Geological Survey: Reston, VA, USA, 1982; Volume 36.
  60. Ji, Z.-G. Hydrodynamics and Water Quality: Modeling Rivers, Lakes, and Estuaries; John Wiley & Sons: Hoboken, NJ, USA, 2017. [Google Scholar]
  61. Annika, P.; George, T.; George, P.; Konstantinos, N.; Costas, D.; Koutitas, C. The Poseidon Operational Tool for the Prediction of Floating Pollutant Transport. Mar. Pollut. Bull. 2001, 43, 270–278. [Google Scholar] [CrossRef]
  62. Nittis, K.; Perivoliotis, L.; Korres, G.; Tziavos, C.; Thanos, I. Operational Monitoring and Forecasting for Marine Environmental Applications in the Aegean Sea. Environ. Model. Softw. 2006, 21, 243–257. [Google Scholar] [CrossRef]
  63. Ainsworth, C.H.; Chassignet, E.P.; French-McCay, D.; Beegle-Krause, C.J.; Berenshtein, I.; Englehardt, J.; Fiddaman, T.; Huang, H.; Huettel, M.; Justic, D.; et al. Ten Years of Modeling the Deepwater Horizon Oil Spill. Environ. Model. Softw. 2021, 142, 105070. [Google Scholar] [CrossRef]
  64. Keramea, P.; Spanoudaki, K.; Zodiatis, G.; Gikas, G.; Sylaios, G. Oil Spill Modeling: A Critical Review on Current Trends, Perspectives, and Challenges. J. Mar. Sci. Eng. 2021, 9, 181. [Google Scholar] [CrossRef]
  65. Zhao, L.; Nedwed, T.; Mitchell, D. Review of the Science behind Oil Spill Fate Models: Are Updates Needed? Int. Oil Spill Conf. Proc. 2021, 2021, 687874. [Google Scholar] [CrossRef]
  66. Gong, Y.; Zhao, X.; Cai, Z.; O’Reilly, S.E.; Hao, X.; Zhao, D. A Review of Oil, Dispersed Oil and Sediment Interactions in the Aquatic Environment: Influence on the Fate, Transport and Remediation of Oil Spills. Mar. Pollut. Bull. 2014, 79, 16–33. [Google Scholar] [CrossRef] [PubMed]
  67. Guyomarch, J.; Le Floch, S.; Merlin, F.-X. Effect of Suspended Mineral Load, Water Salinity and Oil Type on the Size of Oil–Mineral Aggregates in the Presence of Chemical Dispersant. Spill Sci. Technol. Bull. 2002, 8, 95–100. [Google Scholar] [CrossRef]
  68. Khelifa, A.; Stoffyn-Egli, P.; Hill, P.S.; Lee, K. Effects of Salinity and Clay Type on Oil–Mineral Aggregation. Mar. Environ. Res. 2005, 59, 235–254. [Google Scholar] [CrossRef]
  69. Wood, P.A.; Lunel, T.; Daniel, F.; Swannell, R.; Lee, K.; Stoffyn-Egli, P. Influence of Oil and Mineral Characteristics on Oil-Mineral Interaction. In Proceedings of the Arctic and Marine Oilspill Program Technical Seminar, Edmonton, AB, Canada, 10–12 June 1998. [Google Scholar]
  70. Khelifa, A.; Stoffyn-Egli, P.; Hill, P.S.; Lee, K. Characteristics of Oil Droplets Stabilized by Mineral Particles: Effects of Oil Type and Temperature. Spill Sci. Technol. Bull. 2002, 8, 19–30. [Google Scholar] [CrossRef]
  71. Zhang, H.; Khatibi, M.; Zheng, Y.; Lee, K.; Li, Z.; Mullin, J.V. Investigation of OMA Formation and the Effect of Minerals. Mar. Pollut. Bull. 2010, 60, 1433–1441. [Google Scholar] [CrossRef]
  72. Wang, W.; Zheng, Y.; Li, Z.; Lee, K. PIV Investigation of Oil–Mineral Interaction for an Oil Spill Application. Chem. Eng. J. 2011, 170, 241–249. [Google Scholar] [CrossRef]
  73. Khelifa, A.; Stoffyn-Egli, P.; Hill, P.S.; Lee, K. Characteristics of Oil Droplets Stabilized by Mineral Particles: The Effect of Salinity. Int. Oil Spill Conf. Proc. 2003, 2003, 963–970. [Google Scholar] [CrossRef]
  74. Le-Floch, S.; Guyomarch, J.; Merlin, F.-X.; Stoffyn-Egli, P.; Dixon, J.; Lee, K. The Influence of Salinity on Oil–Mineral Aggregate Formation. Spill Sci. Technol. Bull. 2002, 8, 65–71. [Google Scholar] [CrossRef]
  75. Kerebel, D. Study of the Influence of Salinity on the Flocculation Oil-Clay; Centre de Documentation de Recherche et D’Expérimentations sur les Pollutions Accidentelles des Eaux (Cèdre): Brest, France, 1997. [Google Scholar]
  76. Payne, J.R. Oil-Ice-Sediment Interactions during Freezeup and Breakup; US Department of Commerce, National Oceanic and Atmospheric Administration, Ocean Assessments Division, Alaska Office: Anchorage, AK, USA, 1989; Volume 64.
  77. Guyomarch, J.; Merlin, F.-X.; Bernanose, P. Oil Interaction with Mineral Fines and Chemical Dispersion: Behaviour of the Dispersed Oil in Coastal or Estuarine Conditions. In Proceedings of the Arctic and Marine Oilspill Program Technical Seminar, Ministry of Supply and Services, Ottawa, ON, Canada, 1 January 1999; Volume 1, pp. 137–150. [Google Scholar]
  78. Ma, X.; Cogswell, A.; Li, Z.; Lee, K. Particle Size Analysis of Dispersed Oil and Oil-mineral Aggregates with an Automated Ultraviolet Epi-fluorescence Microscopy System. Environ. Technol. 2008, 29, 739–748. [Google Scholar] [CrossRef]
  79. Sun, J.; Khelifa, A.; Zhao, C.; Zhao, D.; Wang, Z. Laboratory Investigation of Oil–Suspended Particulate Matter Aggregation under Different Mixing Conditions. Sci. Total Environ. 2014, 473–474, 742–749. [Google Scholar] [CrossRef] [PubMed]
  80. Sun, J.; Khelifa, A.; Zheng, X.; Wang, Z.; So, L.L.; Wong, S.; Yang, C.; Fieldhouse, B. A Laboratory Study on the Kinetics of the Formation of Oil-Suspended Particulate Matter Aggregates Using the NIST-1941b Sediment. Mar. Pollut. Bull. 2010, 60, 1701–1707. [Google Scholar] [CrossRef] [PubMed]
  81. Lee, K.; Li, Z.; King, T.; Kepkay, P.; Boufadel, M.C.; Venosa, A.D. Wave Tank Studies on Formation and Transport of OMA from the Chemically Dispersed Oil. In Oil Spill Response: A Global Perspective; Davidson, W.F., Lee, K., Cogswell, A., Eds.; Springer: Dordrecht, The Netherlands, 2008; pp. 159–177. [Google Scholar]
  82. Khelifa, A.; Fingas, M.; Brown, C. Effects of Dispersants on Oil-SPM Aggregation and Fate in US Coastal Waters; Final Report Grant Number: NA04NOS4190063; Environment Canada: Ottawa, ON, Canada, 2008. [Google Scholar]
  83. Fu, J.; Gong, Y.; Zhao, X.; O’Reilly, S.E.; Zhao, D. Effects of Oil and Dispersant on Formation of Marine Oil Snow and Transport of Oil Hydrocarbons. Environ. Sci. Technol. 2014, 48, 14392–14399. [Google Scholar] [CrossRef] [PubMed]
  84. O’Laughlin, C.M.; Law, B.A.; Zions, V.S.; King, T.L.; Robinson, B.; Wu, Y. Settling of Dilbit-Derived Oil-Mineral Aggregates (OMAs) & Transport Parameters for Oil Spill Modelling. Mar. Pollut. Bull. 2017, 124, 292–302. [Google Scholar] [CrossRef]
  85. Sterling, M.C.; Bonner, J.S.; Ernest, A.N.S.; Page, C.A.; Autenrieth, R.L. Application of Fractal Flocculation and Vertical Transport Model to Aquatic Sol–Sediment Systems. Water Res. 2005, 39, 1818–1830. [Google Scholar] [CrossRef]
  86. Lee, K. Oil–Particle Interactions in Aquatic Environments: Influence on the Transport, Fate, Effect and Remediation of Oil Spills. Spill Sci. Technol. Bull. 2002, 8, 3–8. [Google Scholar] [CrossRef]
  87. Waterman, D.M.; Garcia, M.H. Laboratory Tests of Oil-Particle Interactions in a Freshwater Riverine Environment with Cold Lake Blend Weathered Bitumen; Ven Te Chow Hydrosystems Laboratory, University of Illinois: Champaign, IL, USA, 2015. [Google Scholar]
  88. Ye, L.; Manning, A.J.; Hsu, T.-J. Oil-Mineral Flocculation and Settling Velocity in Saline Water. Water Res. 2020, 173, 115569. [Google Scholar] [CrossRef]
  89. Bandara, U.C.; Yapa, P.D.; Xie, H. Fate and Transport of Oil in Sediment Laden Marine Waters. J. Hydro-Environ. Res. 2011, 5, 145–156. [Google Scholar] [CrossRef]
  90. Khelifa, A.; Lee, K.; Hill, P.S. Prediction of Oil Droplet Size Distribution in Agitated Aquatic Environments. In Proceedings of the Arctic and Marine Oilspill Program Technical Seminar, Edmonton, AB, Canada, 8–10 June 2004. [Google Scholar]
  91. Khelifa, A.; Lee, K.; Hill, P.S.; Ajijolaiya, L.O. Modelling the Effect of Sediment Size on OMA Formation. In Proceedings of the Arctic and Marine Oilspill Program Technical Seminar, Edmonton, AB, Canada, 8–10 June 2004. [Google Scholar]
  92. Sterling, M.C.; Bonner, J.S.; Page, C.A.; Fuller, C.B.; Ernest, A.N.S.; Autenrieth, R.L. Modeling Crude Oil Droplet—Sediment Aggregation in Nearshore Waters. Environ. Sci. Technol. 2004, 38, 4627–4634. [Google Scholar] [CrossRef]
  93. Zhao, L.; Boufadel, M.C.; Geng, X.; Lee, K.; King, T.; Robinson, B.; Fitzpatrick, F. A-DROP: A Predictive Model for the Formation of Oil Particle Aggregates (OPAs). Mar. Pollut. Bull. 2016, 106, 245–259. [Google Scholar] [CrossRef]
  94. Smoluchowski, M.V. Versuch Einer Mathematischen Theorie Der Koagulationskinetik Kolloider Lösungen. Z. Phys. Chem. 1918, 92, 129–168. [Google Scholar] [CrossRef] [Green Version]
  95. Pruppacher, H.R.; Klett, J.D. Microphysics of Clouds and Precipitation; Atmospheric and Oceanographic Sciences Library; Springer: Dordrecht, The Netherlands, 1980; Volume 284. [Google Scholar]
  96. Ernest, A.N.; Bonner, J.S.; Autenrieth, R.L. Determination of Particle Collision Efficiencies for Flocculent Transport Models. J. Environ. Eng. 1995, 121, 320–329. [Google Scholar] [CrossRef]
  97. Khelifa, A.; Hill, P.S.; Lee, K. Assessment of Minimum Sediment Concentration for OMA Formation Using a Monte Carlo Model. In Developments in Earth and Environmental Sciences; Al-Azab, M., El-Shorbagy, W., Eds.; Oil Pollution and its Environmental Impact in the Arabian Gulf Region; Elsevier: Amsterdam, The Netherlands, 2005; Volume 3, pp. 93–104. [Google Scholar]
  98. Hill, P.S.; Khelifa, A.; Lee, K. Time Scale for Oil Droplet Stabilization by Mineral Particles in Turbulent Suspensions. Spill Sci. Technol. Bull. 2002, 8, 73–81. [Google Scholar] [CrossRef]
  99. Ajijolaiya, L.O.; Hill, P.S.; Khelifa, A.; Islam, R.M.; Lee, K. Laboratory Investigation of the Effects of Mineral Size and Concentration on the Formation of Oil–Mineral Aggregates. Mar. Pollut. Bull. 2006, 52, 920–927. [Google Scholar] [CrossRef] [PubMed]
  100. Wang, Z.; Yanqiu, Z.; Zhiyu, Y.; Bing, S.; Hui, L. Adsorption Model of Oil Droplets Interacting with Suspended Particulate Materials in Marine Coastal Environments. Cont. Shelf Res. 2019, 173, 87–92. [Google Scholar] [CrossRef]
  101. Wu, Y.; Hannah, C.; Law, B.; King, T.; Robinson, B. An Estimate of the Sinking Rate of Spilled Diluted Bitumen in Sediment-Laden Coastal Waters. In Proceedings of the 39th AMOP Technical Seminar, Environment and Climate Change Canada, Ottawa, ON, Canada, 7–9 June 2016; pp. 331–347. [Google Scholar]
  102. Cui, L.; Harris, C.K.; Tarpley, D.R.N. Formation of Oil-Particle-Aggregates: Numerical Model Formulation and Calibration. Front. Mar. Sci. 2021, 8, 510. [Google Scholar] [CrossRef]
  103. Jones, L.; Garcia, M.H. Development of a Rapid Response Riverine Oil—Particle Aggregate Formation, Transport, and Fate Model. J. Environ. Eng. 2018, 144, 04018125. [Google Scholar] [CrossRef]
  104. Wang, Z.; Yang, W.; Zhang, Y.; Yan, Z.; Liu, H.; Sun, B. A Practical Adsorption Model for the Formation of Submerged Oils under the Effect of Suspended Sediments. RSC Adv. 2019, 9, 15785–15790. [Google Scholar] [CrossRef] [Green Version]
  105. Niu, H.; Li, Z.; Lee, K.; Kepkay, P.; Mullin, J.V. Modelling the Transport of Oil–Mineral-Aggregates (OMAs) in the Marine Environment and Assessment of Their Potential Risks. Environ. Model. Assess. 2011, 16, 61–75. [Google Scholar] [CrossRef]
  106. Niu, H.; Lee, K. Study of the Dispersion/Settling of Oil-Mineral-Aggregates Using Particle Tracking Model and Assessment of Their Potential Risks. Int. J. Environ. Pollut. 2013, 52, 32–51. [Google Scholar] [CrossRef]
  107. Ouartassi, B.; Doyon, B.; Heniche, M. Numerical Prediction of Oil Mineral Aggregates Dispersion in the Estuary Of ST-Lawrence River. J. Phys. Conf. Ser. 2021, 1743, 012033. [Google Scholar] [CrossRef]
  108. Zhu, Z.; Waterman, D.M.; Garcia, M.H. Modeling the Transport of Oil–Particle Aggregates Resulting from an Oil Spill in a Freshwater Environment. Environ. Fluid Mech. 2018, 18, 967–984. [Google Scholar] [CrossRef]
  109. Winterwerp, J.C. A Simple Model for Turbulence Induced Flocculation of Cohesive Sediment. J. Hydraul. Res. 1998, 36, 309–326. [Google Scholar] [CrossRef]
  110. Li, Y.; Cao, R.; Chen, H.; Mu, L.; Lv, X. Impact of Oil−sediment Interaction on Transport of Underwater Spilled Oil in the Bohai Sea. Ocean Eng. 2022, 247, 110687. [Google Scholar] [CrossRef]
  111. Zhao, L.; Boufadel, M.C.; Adams, E.; Socolofsky, S.A.; King, T.; Lee, K.; Nedwed, T. Simulation of Scenarios of Oil Droplet Formation from the Deepwater Horizon Blowout. Mar. Pollut. Bull. 2015, 101, 304–319. [Google Scholar] [CrossRef]
  112. Niu, H.; Li, Z.; Lee, K.; Kepkay, P.; Mullin, J. A Method for Assessing Environmental Risks of Oil-Mineral-Aggregate to Benthic Organisms. Hum. Ecol. Risk Assess. Int. J. 2010, 16, 762–782. [Google Scholar] [CrossRef]
Figure 1. How OMA form in the nearshore environment (reproduced from Gustitus et al. [22], with kind permission Copyright 2020 from John Wiley & Sons, Inc.).
Figure 1. How OMA form in the nearshore environment (reproduced from Gustitus et al. [22], with kind permission Copyright 2020 from John Wiley & Sons, Inc.).
Jmse 10 00610 g001
Figure 2. General schematic of this review.
Figure 2. General schematic of this review.
Jmse 10 00610 g002
Figure 3. Schematic of OMA (OPA) formation model reproduced from Cui et al. (2021) [102], with kind permission Copyright from Frontiers Media S.A.
Figure 3. Schematic of OMA (OPA) formation model reproduced from Cui et al. (2021) [102], with kind permission Copyright from Frontiers Media S.A.
Jmse 10 00610 g003
Figure 4. Flowchart of the modelling system adapted from Niu et al. [105] with kind permission from the authors.
Figure 4. Flowchart of the modelling system adapted from Niu et al. [105] with kind permission from the authors.
Jmse 10 00610 g004
Table 1. Examples of 3D oil spill models that use a Lagrangian algorithm to simulate oil transport.
Table 1. Examples of 3D oil spill models that use a Lagrangian algorithm to simulate oil transport.
ModelDeveloperReferences
COSMoSEnvironment and Climate Change Canada (ECCC)[43]
COZOILDepartment of the Interior Minerals Management Service[45]
GNOME SuiteNational Academies of Sciences, Engineering, and Medicine (NOAA)[46]
GULFSPILLKFUPM/RI[47]
MEDSLIK/MEDSLIK-IIOceanography Centre of the University of Cyprus (OC-UCY)[48,49,50]
MIKE 21/3Danish Hydraulic Institute (DHI)[51]
MOHIDMARETEC (Marine and Environmental Technology Research Center)[52]
MOTHYMétéo-France[53]
OILMAP/SIMAPASA[54,55]
OILTRANSThe Atlantic Regions’ Coastal Pollution Response (ARCOPOL)[56]
OpenDrift/OpenOilNorwegian Meteorological Institute[13,14]
OSCARSINTEF[57]
OSRABureau of Ocean Energy Management (BOEM)[58,59,60]
POSEIDON-OSMHellenic Centre for Marine Research (HCMR)[61,62]
SPILLCALCTetra Tech EBA[44]
Table 2. Summary of laboratory experiments on the OMA settling velocity.
Table 2. Summary of laboratory experiments on the OMA settling velocity.
ReferencesMain Objectives/MethodsResults
[21]To measure OMA settling velocities using a focused flow reactorSettling velocities ranging from 2.2 to 10.4 mm/s for 100–200 μm OMA
[87]Settling velocity tests conducted in a 1.6 m height Plexiglass settling column.Settling velocities between 1.0–11.2 mm/s, with the most in the between 1.0 and 3.0 mm/s.
[84]Sediment concentration influences settling velocityHigher concentration of suspended sediment (10 vs. 50 mg/L) means greater settling velocity and effective particle density
[88]To explore effects of clay type on OMA structure and settling velocity using the LabSFLOC-2 system and digital microscopyFor low stickiness Kaolinite clay, OMA settling velocity was about twice as small as for pure Kaolinite flocs; for high stickiness Bentonite clay, the OMA settling velocity was greater than for pure Bentonite flocs
Table 3. A summary of equations relevant to modelling OMA formation.
Table 3. A summary of equations relevant to modelling OMA formation.
EquationEquationsDenoteReferences
(1) θ = d n k d t = 1 2 i + j = k α β ( i , j ) n i n j n k i = 1 α β ( i , k ) n i α   is   the   collision   efficiency .   β   is   the   collision   frequency ,   and   n i   and   n j are the particle concentrations for the particles of size i and j, respectively.[94]
(2) β = β B r + β s h + β d s
β B r = 2 k T 3 μ ( 1 d i + 1 d j ) ( d i + d j )
β s h = G 6 ( d i + d j ) 3
G = ( ε μ ) 0.5
β d s = π 4 ( d i + d j ) 2 | ( U i U j ) |
β B r ,   β s h ,   and   β d s   are   the   collision   efficiency   due   to   Brownian   motion ,   fluid   shear ,   differential   sedimentation ,   respectively .   k   is   the   Boltzman s   constant .   T   is   the   absolute   temperature ,   and   μ   is   the   dynamic   viscosity   of   the   media .   d i   and   d j   are   the   effective   particle   diameters .   G   is   the   root   mean   square   of   the   velocity   gradient .   ε   is   the   dissipation   rate .   U i   and   U j   are   the   settling   velocities   of   the   two   collided   entities   of   size   d i   and   d j .[92,93,96]
(3) d n k , q d t = D z 2 n k , q Z 2 w k , q n k , q Z + θ k , q n k , q   is   the   particle   size   distribution .   k   is   particle   volume .   q   is   effective   density .   z   is   the   vertical   distance .   D z   is   the   vertical   dispersion   coefficient .   w k , q   is   the   settling   velocity .   θ k , q is the interaction term due to coagulation.[96]
(4) α o b s ( o , p ) = α H O M O , 1 ( s f 1 , o ) ( s f 1 , p ) α H O M O , 2 ( s f 2 , o ) ( s f 2 , p ) α H E T , 1 2 ( s f 1 , o s f 2 , p s f 1 , p s f 2 , o ) s f 1 , o ,   s f 2 , o ,   s f 1 , p ,   and   s f 2 , p   are   surface   fractions   of   constituent   particles   1   and   2   in   aggregate   with   density   o   and   constituent   particles   1   and   2   in   aggregate   with   density   p .   α H O M O , 1 ,   α H O M O , 2 ,   and   α H E T are probabilities of successful aggregation through contacting floc constituent types 1-1, 2-2, and 1-2.[92]
(5) C i t + ( u k C i ) x k = ( E k C i x k ) x k + w s i C i x 3 + S i , A g g + S i , D e + S i , A b s For   i = 1 5 ,   corresponding   to   5   species   mentioned   in   assumptions .   k = 1 ,   2 ,   3   correspond   to   directions   of   x ,   y ,   z .   C i = volumetric   concentration   of   the   i th   species .   u k   = component   of   current   velocity   in   x ,   y ,   and   z   directions .   w s i   is   the   buoyant / settling   velocity   of   the   i th   species .   E k   is   the   diffusion   coefficient   in   k th   direction .   S i , A g g = C i t | A g g   = source / sin k   terms   of   the   i th   species   due   to   aggregation .   S i , D e = C i t | D e   = sin k   term   due   to   deposition   of   the   i th   species ,   and   S i , A b s = source/sink term for the ith species due to partitioning. [89]
(6) α   ( t ) = α s t a ( 1   A p p r o j   i n   O M A s F S P   A o ) α s t a   is   the   stability   ratio .   A o   is   the   surface   area   of   an   oil   droplet .   A p p r o j   is   the   projection   area   of   particles   on   the   droplet   surface   when   an   OMA   forms ,   and   F S P is a factor to account for particle shape and packing effects on the coagulation process.[93]
(7) P b r k = { 0 ,       n b = 0 0.5 ,   n b = 1     1 ,       n b > 1   n b   is   the   number   of   droplets   larger   than   the   maximum   allowed   D m a x .   A   breakage   process   is   chosen   when   two   or   more   droplets   are   larger   than   D m a x .   P a g g = 1 P b r k .   A   random   number   r l   is   selected   from   a   uniform   distribution   between   1   and   0 .   A   breakage   event   is   selected   if   P b r k   rl, otherwise an aggregation event is selected.[90,91,97]
(8) d N s d t = 0.16 α o s ( D s + D o ) 3 ( ε / v ) 1 2 N s N o D s   is   sediment   diameter .   D o   is   oil   droplet   diameter .   N s   is   OMA   number   concentration .   α o s is   the   collision   efficiency ;   ε is the dissipation rate. v is the kinematic viscosity of water.[98]
(9) E = E m a x ( C s C s 50 ) n 1 + ( C s C s 50 ) n E m a x   is   the   maximum   possible   trapping   efficiency .   C s   is   the   mass   concentration   of   the   sediment ,   C s 50   is   the   sediment   concentration   at   50 %   trapping   efficiency .   n   is   the   shape   of   the   trapping   efficiency   versus   sediment   concentration   curve .   Least - squares   fitting   of   the   equation   to   the   experiment   data   yields   n = 3 ,   and   E m a x = 85%.[99]
(10) E = ( K d / 10 3 ) × S P M 1 + ( K d / 10 3 ) × S P M E   is   the   oil   trapping   efficiency .   K d is a distribution coefficient. SPM is the sediment concentration.[100]
(11) S R = V s V = r c 4 s r m i n 4 s r m a x 4 s r m i n 4 s × 100 % SR   is   the   oil   sin king   ratio .   r c   is   50 %   of   the   critical   oil   droplets   size .   r m a x   and   r m i n are the maximum and minimum radii of the oil droplets. s = 2.3 based on laboratory data.[101]
Table 4. Equations used to estimate the settling velocity of oil-mineral-aggregate.
Table 4. Equations used to estimate the settling velocity of oil-mineral-aggregate.
EquationEquationsDenoteReferences
(12) w s = ( ρ s ρ w ) g D 2 18 μ w s   is   the   settling   velocity   of   OMA .   D   is   the   diameter   of   OMA .   ρ w   and   μ   is   the   density   and   viscosity   of   the   liquid   respectively .   ρ s   is   the   OMA   density .   and   g is the gravity. [109]
(13) w s = 4 g | ρ s ρ w | D 3 C D ρ w
If   R e < 1 ,   C D = 24 R e
If   1 < R e < 10 4 ,   C D = 24 R e + 3 R e + 0.34
If   R e > 10 4 ,   C D = 0.34   t o   0.4
C D   is   the   drag   coefficient .   R e is the Reynold number.[93,111]
(14) w s = e D f e   and   f are regression coefficients.[112]
(15) C s e t t l e d = i = 1 k N i × P M i A c e l l
P M i = M s p i l l e d × P s e t t l e d × p i n i
V o i l = ρ s e d i m e n t ρ O M A ρ s e d i m e n t ρ o i l V O M A
k   is   the   number   of   particle   size   classes .   N i   is   the   number   of   the   i t h - size   class   particles   in   the   user   specified   cell .   A c e l l   is   the   area   of   the   cell ,   and   P M i   is   the   amount   of   oil   per   particle   for   the   i t h - class   particles .   M s p i l l e d   is   the   total   mass   of   spilled   oil .   P s s t t l e d   is   the   percentage   of   spilled   oil   that   may   be   transferred   to   sediment .   n i   is   the   number   of   particles   used   in   the   simulation   for   class   i ,   and   p i   is   the   fraction   of   settled   oil   ( in   percentage )   carried   by   the   particle   class .   V o i l   is   the   volume   of   oil   in   an   OMA   of   class   i .   V O M A   is   the   volume   of   an   OMA .   ρ o i l ,   ρ s e d i m e n t ,   and   ρ O M A are the densities of oil, sediment, and the OMA, respectively.[105]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhong, X.; Niu, H.; Li, P.; Wu, Y.; Liu, L. An Overview of Oil-Mineral-Aggregate Formation, Settling, and Transport Processes in Marine Oil Spill Models. J. Mar. Sci. Eng. 2022, 10, 610. https://doi.org/10.3390/jmse10050610

AMA Style

Zhong X, Niu H, Li P, Wu Y, Liu L. An Overview of Oil-Mineral-Aggregate Formation, Settling, and Transport Processes in Marine Oil Spill Models. Journal of Marine Science and Engineering. 2022; 10(5):610. https://doi.org/10.3390/jmse10050610

Chicago/Turabian Style

Zhong, Xiaomei, Haibo Niu, Pu Li, Yongsheng Wu, and Lei Liu. 2022. "An Overview of Oil-Mineral-Aggregate Formation, Settling, and Transport Processes in Marine Oil Spill Models" Journal of Marine Science and Engineering 10, no. 5: 610. https://doi.org/10.3390/jmse10050610

APA Style

Zhong, X., Niu, H., Li, P., Wu, Y., & Liu, L. (2022). An Overview of Oil-Mineral-Aggregate Formation, Settling, and Transport Processes in Marine Oil Spill Models. Journal of Marine Science and Engineering, 10(5), 610. https://doi.org/10.3390/jmse10050610

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop