Systematic CFD-based evaluation of physical factors influencing the spatiotemporal distribution patterns of microplastic particles in lakes

Spatiotemporal distribution patterns of microplastic (MP) particles in lakes hinge on both the physical conditions in the lake and particle properties. Using numerical simulations, we systematically investigated the influence of lake depth and bathymetry, wind and temperature conditions, MP particle release location and timing, as well as particle diameter (10, 20, and 50 μm). Our results indicate that maximum lake depth had the greatest effect on the residence time in the water column, as it determines the settling timescale and occurrence of hydrodynamic complexity such as density-driven flows in the lake. Increasing particle size from 10 to 20 and 50 μm also significantly reduced the residence time making particle size the factor with the second strongest effect on the residence time and, in turn, on the availability of MP particles for uptake by organisms. Changing bathymetry from a uniform to a non-uniform had a less pronounced effect on particle residence time compared to maximum depth and particle size. Release location, wind conditions, and release time had comparably little effect on particle behavior but became more important as MP particle size decreased. The release of the 10 μm MP particles in the deeper lakes with uniform bathymetry during summer with stable thermal stratification, resulted in a nearly month-long turnover phase in the fall in which both settling and rising of particles occurred simultaneously. This was caused by convective heat and water transport during this period. In these scenarios about 2.6 to 5.4 % of the released MP particles were held in or returned to the water layers near the lake surface. While acknowledging the dominant role of lake depth and MP particle size on the particle residence time, this study further emphasizes that it is ultimately a particular combination of different factors and their interactions that shape MP distribution patterns in lakes.

Storm-induced coastward expansion of Margalefidinium polykrikoides bloom in Chesapeake Bay

An unusual coastward expansion of the toxic dinoflagellate species Margalefidinium polykrikoides was observed in 2020 summer after a tropical storm passing Chesapeake Bay. Such coastward expansion was only recorded in 2007. A newly developed coupled Lagrangian particle tracking and harmful algal bloom model driven by environmental variables was used to investigate the underlying mechanisms and successfully reproduced the expansion patterns. Persistent pre-storm southerly winds favored the delivery of bloom source water originated inside the bay to the coast. Storm-induced strong upwelling of denser subsurface water interacted with the after-storm outflow plume (steered southward as the storm's impacts waned), forming a transport barrier to accumulate algae and delineate the coastwide bloom extent. Algal diel vertical migrations and transport barrier enable algae to stay in the nearshore regions. The storm-induced coastward expansion of M. polykrikoides might increase future bloom possibility in the coastal area.

CFD Simulation of Particle-Laden Flow in a 3D Differentially Heated Cavity Using Coarse Large Eddy Simulation

Particulate flow in closed space is involved in many engineering applications. In this paper, the prediction of particle removal is investigated in a thermally driven 3D cavity at turbulent Rayleigh number Ra = 10⁹ using Coarse Large Eddy Simulation (CLES). The depletion dynamics of SiO₂ aerosol with aerodynamic diameters between 1.4 and 14 µm is reported in an Euler/Lagrange framework. The main focus of this work is therefore to assess the effect of the subgrid-scale motions on the prediction of the particulate flow in a buoyancy driven 3D cavity flow when the mesh resolution is coarse and below optimal LES standards. The research is motivated by the feasibility of modeling more complex particulate flows with reduced CPU cost. The cubical cavity of 0.7 m side-length is set to have a temperature difference of 39 K between the two facing cold and hot vertical walls. As a first step, the carrier fluid flow was validated by comparing the first and second-moment statistics against both previous well-resolved LES and experimental databases [Kalilainen (J. Aero Sci. 100:73-87, 2016); Dehbi (J. Aero. Sci. 103:67-82, 2017)]. First moment Eulerian statistics show a very good match with the reference data both qualitatively and quantitatively, whereas higher moments show underprediction due to the lesser spatial resolution. In a second step, six particle swarms spanning a wide range of particle Stokes numbers were computed to predict particle depletion. In particular, predictions of 1.4 and 3.5 µm particles were compared to LES and available experimental data. Particles of low inertia i.e. dp < 3.5 µm are more affected by the SGS effects, while bigger ones i.e. dp = 3.5-14 µm exhibit much less grid-dependency. Lagrangian statistics reported in both qualitative and quantitative fashions show globally a very good agreement with reference LES and experimental databases at a fraction of the CPU power needed for optimal LES.

A cell-resolved, Lagrangian solver for modeling red blood cell dynamics in macroscale flows

When red blood cells (RBCs) experience non-physiologically high stresses, e.g., in medical devices, they can rupture in a process called hemolysis. Directly simulating this process is computationally unaffordable given that the length scales of a medical device are several orders of magnitude larger than that of a RBC. To overcome this separation of scales, the present work introduces an affordable computational framework that accurately resolves the stress and deformation of a RBC in a spatially and temporally varying macroscale flow field such as those found in a typical medical device. The underlying idea of the present framework is to treat RBCs as one-way coupled tracers in the macroscale flow by capturing the effect of the flow on their dynamics but neglecting their effect on the flow at the macroscale. As a result, the RBC dynamics are simulated after those of the flow in a postprocessing step by receiving the fluid velocity gradient tensor measured along the RBC trajectory as the input. To resolve the fluid velocity in the immediate vicinity of the RBC as well as the motion of the membrane, we employ the boundary integral method coupled to a structural solver. The governing equations are discretized in space using spherical harmonics, yielding spectral integration accuracy. The predictions produced by this formulation are in good agreement with those obtained from simulations of spherical capsules in shear flows and optical tweezers experiments. The accuracy of the present method is evaluated using unbounded shear flow as a benchmark. Its computational cost grows proportional to p ⁵, where p is the degree of the spherical harmonic. It also exhibits a fast convergence rate that is approximately O ( p 6 ) for p ⪅ 20.

Comparison between fully resolved and time-averaged simulations of particle cloud dispersion produced by a violent expiratory event

In this work we compare the DNS results (Fabregat et al. 2021, Fabregat et al. 2021) for a mild cough already reported in the literarure with those obtained with a compressible URANS equations with a k-ϵ turbulence model. In both cases, the dispersed phase has been modelled as spherical Lagrangian particles using the one-way coupling assumption. Overall, the URANS model is capable of reproducing the observed tendency of light particles under 64 µm in diameter to rise due to the action of the drag exerted by the buoyant puff generated by the cough. Both DNS and URANS found that particles above 64 µm will tend to describe parabolic trajectories under the action of gravitational forces. Grid independence analysis allows to qualify the impact of increasing mesh resolution on the particle cloud statistics as flow evolves. Results suggest that the k-ϵ model overpredicts the horizontal displacement of the particles smaller than 64 µm while the opposite occurs for the particles larger than 64 µm.

The role of wind in controlling the connectivity of blue mussels (Mytilus edulis L.) populations

BACKGROUND

Larval connectivity between distinct benthic populations is essential for their persistence. Although connectivity is difficult to measure in situ, it can be predicted via models that simulate biophysical interactions between larval behaviour and ocean currents. The blue mussel (Mytilus Edulis L.) is widespread throughout the northern hemisphere and extensively commercialised worldwide. In the Irish Sea, this industry represents ~ 50% of Welsh shellfisheries, where cultivation is mainly based on wild spat. However, the main sources and amount of spat varied interannually (1100 tonnes harvest in 2014 against zero in 2018). The aim of this study is to characterise the structure and dynamics of the blue mussel metapopulation within the northern part of the Irish Sea.

METHODS

We develop a Lagrangian particle tracking model, driven by a high-resolution (from 30 to 5000 m) validated unstructured coastal hydrodynamic model of the Irish Sea, to simulate spatial and temporal variability of larval dispersal and connectivity between distinct mussel populations and potential settlement areas.

RESULTS

Our results showed that: (1) larvae positioned near the surface were strongly influenced by wind-driven currents suggesting that connectivity networks had the potential to span hundreds of kilometres; (2) in contrast, larvae positioned deeper in the water column were driven by tidal currents, producing intricate spatial patterns of connectivity between mussel beds over tens of kilometres that were consistent over time.

CONCLUSIONS

Dispersal of mussel larvae in the tidally energetic Irish Sea during the April-May spawning season is potentially driven by wind-driven surface currents, as confirmed by fisherman observations of inter-annual variability in wild spat collection. These results have important implications for metapopulation dynamics within the context of climate change and sustainable shellfisheries management (i.e. gain and loss of populations and harvest areas according to wind conditions).

In-Silico Conceptualisation of Continuous Millifluidic Separators for Magnetic Nanoparticles

Magnetic nanoparticles are researched intensively not only for biomedical applications, but also for industrial applications including wastewater treatment and catalytic processes. Although these particles have been shown to have interesting surface properties in their bare form, their magnetisation remains a key feature, as it allows for magnetic separation. This makes them a promising carrier for precious materials and enables recovery via magnetic fields that can be turned on and off on demand, rather than using complex (nano)filtration strategies. However, designing a magnetic separator is by no means trivial, as the magnetic field and its gradient, the separator dimensions, the particle properties (such as size and susceptibility), and the throughput must be coordinated. This is showcased here for a simple continuous electromagnetic separator design requiring no expensive materials or equipment and facilitating continuous operation. The continuous electromagnetic separator chosen was based on a current-carrying wire in the centre of a capillary, which generated a radially symmetric magnetic field that could be described using cylindrical coordinates. The electromagnetic separator design was tested in-silico using a Lagrangian particle-tracking model accounting for hydrodynamics, magnetophoresis, as well as particle diffusion. This computational approach enabled the determination of separation efficiencies for varying particle sizes, magnetic field strengths, separator geometries, and flow rates, which provided insights into the complex interplay between these design parameters. In addition, the model identified the separator design allowing for the highest separation efficiency and determined the retention potential in both single and multiple separators in series. The work demonstrated that throughputs of ~1/4 L/h could be achieved for 250–500 nm iron oxide nanoparticle solutions, using less than 10 separator units in series.

Investigating Pseudo-nitzschia australis introduction to the Gulf of Maine with observations and models

In 2016, an unprecedented Pseudo-nitzschia australis bloom in the Gulf of Maine led to the first shellfishery closures due to domoic acid in the region's history. In this paper, potential introduction routes of P. australis are explored through observations, a hydrodynamic model, and a Lagrangian particle tracking model. Based on particle tracking experiments, the most likely source of P. australis to the Gulf of Maine was the Scotian Shelf. However, in 2016, connectivity between the Scotian Shelf and the bloom region was not significantly different from the other years between 2012 and 2019, nor were temperature conditions more favorable for P. australis growth. Observations indicated changes on the Scotian Shelf in 2016 preceded the introduction of P. australis: increased bottom salinity and decreased surface salinity. The increased bottom salinity on the shelf may be linked to anomalously saline water observed near the coast of Maine in 2016 via transport through Northeast Channel. The changes in upstream water mass properties may be related to the introduction of P. australis, and could be the result of either increased influence of the Labrador Current or increased outflow from the Gulf of St. Lawrence. The ultimate source of P. australis remains unknown, although the species has previously been observed in the eastern North Atlantic, and connectivity across the ocean is possible via a subpolar route. Continued and increased monitoring is warranted to track interannual Pseudo-nitzschia persistence in the Gulf of Maine, and sampling on the Scotian Shelf should be conducted to map upstream P. australis populations.

Three-dimensional numerical modelling of transport, fate and distribution of microplastics in the northwestern Arabian/Persian Gulf

Marine plastic litter has been a major concern over the past decade particularly in semi-enclosed seas such as the Arabian/Persian Gulf, which are likely to impose a relatively higher threat to ecosystem and human health. In this work, we have focused our efforts on the transport features of marine surface microplastics (MPs) in the Gulf. The assessment utilizes a 3D hydrodynamic model of the northern Gulf which was coupled with a particle tracking model. We have considered five release locations and investigated two dominant wind conditions by applying different numerical scenarios. The results revealed that the northerly winds result in high dispersion and seaward transport of MPs in the open coastal zones, while in semi-enclosed regions they result in high trapping and beaching verified by visual investigation. The study shows that further detailed field investigations are warranted to enable the models to better parameterize the fate and distributions of MPs.

An optimal 𝒪(N) scheme for simulations of colliding, particle-laden flows on unstructured grids

The cost of tracking Lagrangian particles in a domain discretized on an unstructured grid can become prohibitively expensive as the number of particles or elements grows. A major part of the cost in these calculations is spent on locating the element that hosts a particle and detecting binary collisions, with the latter traditionally requiring 𝒪 ( N 2 ) operations, N being the number of particles. This paper introduces an optimal search box strategy to significantly reduce the cost of these two operations, ensuring a nearly 𝒪 ( N ) scaling of the cost of collision detection for large-scale simulations. The particle localization strategy is constructed by obtaining an a priori estimate for the optimal number of search boxes as a function of the number of elements, particles, and time steps. The introduced method is generic, as it must be tuned only once for a given implementation and element type. The optimal number of search boxes for collision detection, although complex in form, can be reasonably approximated as the number of particles. The optimality of our method is shown using three drastically varying geometries.

Impact of tidally induced residual circulations on chemical oxygen demand (COD) distribution in Laizhou Bay, China

To understand the impact of hydrodynamics on pollutant transport in Laizhou Bay, China, we conducted numerical simulations using Mike 21. The model was calibrated with good agreements to field monitoring data at various monitoring stations. The simulation results show a clockwise and an anti-clockwise tidally-induced residual circulation in the western and eastern bay, respectively. Historical COD monitoring data also indicate two rings of high COD concentration in the same regions of the bay. This suggests that the hydrodynamics of tidal and residual currents is the main cause of the ring-shaped high COD concentration field in the bay. Pollutant inputs from inland rivers are also important for the COD distribution, making the near-shore side of the COD ring higher than the offshore side. Regions with higher retention time in the bay are usually associated with higher COD concentrations. This study is useful in understanding the mechanism of pollutant spatial distribution and subsequent pollution control in a sea bay.

The influence of cascade hydropower development on the hydrodynamic conditions impacting the reproductive process of fish with semi-buoyant eggs

Hydropower projects have changed the physical habitat of rivers, which has a serious impact on the survival of local fish. The reproduction of fish producing semi-buoyant eggs requires a specific hydrodynamic condition. To predict the influence of cascade hydropower development on the reproductive process of fish with semi-buoyant eggs, a MIKE 21 Flow Model was applied. Lagrangian particle tracking was used to simulate the movement of semi-buoyant eggs and larvae using the Agent-based Modeling (ABM) lab module. The calibrated model showed good agreement between the simulated and observed data for the hydrodynamic process in the reservoir. Twelve scenarios were defined to fully understand whether fish with semi-buoyant eggs can reproduce naturally. The results showed the following: (1) It is difficult to form a fish migration passage in the reservoir in all potential scenarios. (2) Semi-buoyant eggs and larvae sank to the bottom and perished before they hatched and were old enough to survive, since the hydrodynamic conditions could not meet the minimum flow velocity required to keep them in the drift. (3) Even if the hydrodynamic conditions can keep them in the drift in impossible high-discharge conditions, there was not enough drifting time and distance in the reservoir. The results implied that fishes with semi-buoyant eggs cannot reproduce naturally in the main stream, but it is possible that they can reproduce successfully after the protection of the tributary. The method is transferrable to other locations via establishment of models with relevant data to a particular area.

A novel mechanism of mixing by pulsing corals

The dynamic pulsation of xeniid corals is one of the most fascinating phenomena observed in coral reefs. We quantify for the first time the flow near the tentacles of these soft corals, the active pulsations of which are thought to enhance their symbionts' photosynthetic rates by up to an order of magnitude. These polyps are approximately 1 cm in diameter and pulse at frequencies between approximately 0.5 and 1 Hz. As a result, the frequency-based Reynolds number calculated using the tentacle length and pulse frequency is on the order of 10 and rapidly decays as with distance from the polyp. This introduces the question of how these corals minimize the reversibility of the flow and bring in new volumes of fluid during each pulse. We estimate the Péclet number of the bulk flow generated by the coral as being on the order of 100-1000 whereas the flow between the bristles of the tentacles is on the order of 10. This illustrates the importance of advective transport in removing oxygen waste. Flow measurements using particle image velocimetry reveal that the individual polyps generate a jet of water with positive vertical velocities that do not go below 0.1 cm s^-1 and with average volumetric flow rates of approximately 0.71 cm³ s^-1 Our results show that there is nearly continual flow in the radial direction towards the polyp with only approximately 3.3% back flow. 3D numerical simulations uncover a region of slow mixing between the tentacles during expansion. We estimate that the average flow that moves through the bristles of the tentacles is approximately 0.03 cm s^-1 The combination of nearly continual flow towards the polyp, slow mixing between the bristles, and the subsequent ejection of this fluid volume into an upward jet ensures the polyp continually samples new water with sufficient time for exchange to occur.

In Vitro Study of Particle Transport in Successively Bifurcating Vessels

To reach a predictive understanding of how particles travel through bifurcating vessels is of paramount importance in many biomedical settings, including embolization, thromboembolism, and drug delivery. Here we utilize an in vitro model in which solid particles are injected through a rigid vessel that symmetrically bifurcates in successive branching generations. The geometric proportion and fluid dynamics parameters are relevant to the liver embolization. The volumetric flow field is reconstructed via phase-contrast magnetic resonance imaging, from which the particle trajectories are calculated for a range of size and density using the particle equation of motion. The method is validated by directly tracking the injected particles via optical imaging. The results indicate that, opposite to the common assumption, the particles distribution is fundamentally different from the volumetric flow partition. In fact, the amount of delivered particles vary substantially between adjacent branches even when the flow is uniformly distributed. This is not due to the inertia of the particles, nor to gravity. The particle distribution is rather rooted in their different pathways, which in turn are linked to their release origin along the main vessel cross-section. Therefore, the tree geometry and the associated flow streamlines are the prime determinant of the particle fate, while local changes of volumetric flow rate to selected branches do not generally produce proportional changes of particle delivery.

Global-scale dispersal and connectivity in mangroves

Mangroves are of considerable ecological and socioeconomical importance; however, substantial areal losses have been recorded in many regions, driven primarily by anthropogenic disturbances and sea level rise. Oceanic dispersal of mangrove propagules provides a key mechanism for shifting distributions in response to environmental change. Here we provide a model framework for describing global dispersal and connectivity in mangroves. We identify important dispersal routes, barriers, and stepping-stones and quantify the influence of minimum and maximum floating periods on simulated connectivity patterns. Our study provides a baseline to improve our understanding of observed mangrove species distributions and, in combination with climate data, the expected range shifts under climate change.

Dispersal provides a key mechanism for geographical range shifts in response to changing environmental conditions. For mangroves, which are highly susceptible to climate change, the spatial scale of dispersal remains largely unknown. Here we use a high-resolution, eddy- and tide-resolving numerical ocean model to simulate mangrove propagule dispersal across the global ocean and generate connectivity matrices between mangrove habitats using a range of floating periods. We find high rates of along-coast transport and transoceanic dispersal across the Atlantic, Pacific, and Indian Oceans. No connectivity is observed between populations on either side of the American and African continents. Archipelagos, such as the Galapagos and those found in Polynesia, Micronesia, and Melanesia, act as critical stepping-stones for dispersal across the Pacific Ocean. Direct and reciprocal dispersal routes across the Indian Ocean via the South Equatorial Current and seasonally reversing monsoon currents, respectively, allow connectivity between western Indian Ocean and Indo-West Pacific sites. We demonstrate the isolation of the Hawaii Islands and help explain the presence of mangroves on the latitudinal outlier Bermuda. Finally, we find that dispersal distance and connectivity are highly sensitive to the minimum and maximum floating periods. We anticipate that our findings will guide future research agendas to quantify biophysical factors that determine mangrove dispersal and connectivity, including the influence of ocean surface water properties on metabolic processes and buoyancy behavior, which may determine the potential of viably reaching a suitable habitat. Ultimately, this will lead to a better understanding of global mangrove species distributions and their response to changing climate conditions.

Accounting for residence-time in blood rheology models: do we really need non-Newtonian blood flow modelling in large arteries?

Patient-specific computational fluid dynamics (CFD) is a promising tool that provides highly resolved haemodynamics information. The choice of blood rheology is an assumption in CFD models that has been subject to extensive debate. Blood is known to exhibit shear-thinning behaviour, and non-Newtonian modelling has been recommended for aneurysmal flows. Current non-Newtonian models ignore rouleaux formation, which is the key player in blood's shear-thinning behaviour. Experimental data suggest that red blood cell aggregation and rouleaux formation require notable red blood cell residence-time (RT) in a low shear rate regime. This study proposes a novel hybrid Newtonian and non-Newtonian rheology model where the shear-thinning behaviour is activated in high RT regions based on experimental data. Image-based abdominal aortic and cerebral aneurysm models are considered and highly resolved CFD simulations are performed using a minimally dissipative solver. Lagrangian particle tracking is used to define a backward particle RT measure and detect stagnant regions with increased rouleaux formation likelihood. Our novel RT-based non-Newtonian model shows a significant reduction in shear-thinning effects and provides haemodynamic results qualitatively identical and quantitatively close to the Newtonian model. Our results have important implications in patient-specific CFD modelling and suggest that non-Newtonian models should be revisited in large artery flows.

Modeling oil spill trajectory in Bosphorus for contingency planning

Bosphorus, is a strongly driven international maritime route between the Black Sea and the Sea of Marmara and is a high risk area for oil spill due to the heavy tanker traffic. In this study an oil spill trajectory model was developed for investigating the potential risks of accidental oil spills in Bosphorus. The proposed oil spill trajectory model combines the surface current velocity data obtained from a calibrated hydrodynamic model with the advection, spreading, and evaporation processes that are effective only on the sea surface and dominant for a couple of hours after the oil spill. Model simulations revealed that spilled oil reaches the shoreline on both sides of Bosphorus in <4h following the spill. We proposed locations for emergency intervention stations in Bosphorus which can be used to devise a suitable oil spill contingency plan to keep the adverse impacts of oil spills at minimum.

Combining hydrodynamic modelling with genetics: can passive larval drift shape the genetic structure of Baltic Mytilus populations?

While secondary contact between Mytilus edulis and Mytilus trossulus in North America results in mosaic hybrid zone formation, both species form a hybrid swarm in the Baltic. Despite pervasive gene flow, Baltic Mytilus species maintain substantial genetic and phenotypic differentiation. Exploring mechanisms underlying the contrasting genetic composition in Baltic Mytilus species will allow insights into processes such as speciation or adaptation to extremely low salinity. Previous studies in the Baltic indicated that only weak interspecific reproductive barriers exist and discussed the putative role of adaptation to environmental conditions. Using a combination of hydrodynamic modelling and multilocus genotyping, we investigate how oceanographic conditions influence passive larval dispersal and hybrid swarm formation in the Baltic. By combining our analyses with previous knowledge, we show a genetic transition of Baltic Mytilus species along longitude 12°-13°E, that is a virtual line between Malmö (Sweden) and Stralsund (Germany). Although larval transport only occurs over short distances (10-30 km), limited larval dispersal could not explain the position of this genetic transition zone. Instead, the genetic transition zone is located at the area of maximum salinity change (15-10 psu). Thus, we argue that selection results in weak reproductive barriers and local adaptation. This scenario could maintain genetic and phenotypic differences between Baltic Mytilus species despite pervasive introgressive hybridization.

Dredging for dilution: A simulation based case study in a Tidal River

A 2-D hydrodynamic finite element model with a Lagrangian particle module is used to investigate the effects of dredging on the hydrodynamics and the horizontal dilution of pollutant particles originating from a wastewater treatment facility (WWTF) in tidal Oyster River in New Hampshire, USA. The model is driven by the semi-diurnal (M2) tidal component and includes the effect of flooding and drying of riverine mud flats. The particle tracking method consists of tidal advection plus a horizontal random walk model of sub-grid scale turbulent processes. Our approach is to perform continuous pollutant particle releases from the outfall, simulating three different scenarios: a base-case representing the present conditions and two different dredged channel/outfall location configurations. Hydrodynamics are investigated in an Eulerian framework and Lagrangian particle dilution improvement ratios are calculated for all cases. Results show that the simulated hydrodynamics are consistent with observed conditions. Eulerian and Lagrangian residuals predict an outward path suggesting flushing of pollutants on longer (>M2) time scales. Simulated dilution maps show that, in addition to dredging, the relocation of the WWTF outfall into the dredged main channel is required for increased dilution performance. The methodology presented here can be applied to similar managerial problems in all similar systems worldwide with relatively little effort, with the combination of Lagrangian and Eulerian methods working together towards a better solution. The statistical significance brought into methodology, by using a large number of particles (16000 in this case), is to be emphasized, especially with the growing number of networked parallel computer clusters worldwide. This paper improves on the study presented in Bilgili et al., 2006b, by adding an Eulerian analysis.

Intercomparison of oil spill prediction models for accidental blowout scenarios with and without subsea chemical dispersant injection

We compare oil spill model predictions for a prototype subsea blowout with and without subsea injection of chemical dispersants in deep and shallow water, for high and low gas-oil ratio, and in weak to strong crossflows. Model results are compared for initial oil droplet size distribution, the nearfield plume, and the farfield Lagrangian particle tracking stage of hydrocarbon transport. For the conditions tested (a blowout with oil flow rate of 20,000 bbl/d, about 1/3 of the Deepwater Horizon), the models predict the volume median droplet diameter at the source to range from 0.3 to 6mm without dispersant and 0.01 to 0.8 mm with dispersant. This reduced droplet size owing to reduced interfacial tension results in a one to two order of magnitude increase in the downstream displacement of the initial oil surfacing zone and may lead to a significant fraction of the spilled oil not reaching the sea surface.