1. Introduction
The drift of a wide range of floating objects in geophysical settings follows separate empirical rules of thumb that all predict motion at 2%–4% of the wind velocity, relative to the velocity of the water. For example, observations consistently show that freely drifting Arctic icebergs travel at the water velocity plus 1.6%–1.8% of the near-surface wind velocity (Smith and Banke 1983; Garrett et al. 1985; Bigg et al. 1997). Unconstrained sea ice floes, similarly, have long been observed to typically drift at 2%–2.5% of the wind velocity relative to the ocean current (Nansen 1902; Zubov 1945; Browne and Crary 1958; Thorndike and Colony 1982). Comparable drift behavior has separately been found for mangrove propagules in laboratory flume experiments (Van der Stocken et al. 2015). The drift of human survivors of ship wreckages as well as that of man-made objects such as life rafts, cargo containers, and other types of flotsam have been the subject of an extensive body of research that is often referred to as “search and rescue” literature and that reports similar drift behavior (e.g., Allen and Plourde 1999; Daniel et al. 2002; Allen 2005; Breivik et al. 2011, 2012; Röhrs et al. 2012; Breivik et al. 2013; Nesterov 2018; Sutherland et al. 2020). Oil slicks have also been observed to drift at approximately 3% of the wind velocity relative to the ocean surface (Stolzenbach et al. 1977; Rasmussen 1985). Despite the wide range of shapes, materials, and aspect ratios, all of these types of objects demonstrate similar drift behavior.
Here, we consider the origin of these approximately equivalent wind sensitivities. We begin with an idealized theoretical analysis of the momentum balance, focusing on the water and air drag forces on a partially submerged rectangular object. This general approach, which results in the drift velocity of a floating object being given as that of the near-surface ocean current velocity plus a small percentage of the ambient wind velocity (the leeway factor), has been referred to as the “leeway modeling approach” (Olascoaga et al. 2020). Our idealized theoretical considerations (sections 2 and 3) thus share a number of underlying assumptions with some previous work on leeway drift modeling (e.g., Daniel et al. 2002; Nesterov 2018), as discussed below. The results provide a physically intuitive explanation for the 2%–4% drift laws. Our framework also allows us to assess under what conditions the sensitivity to wind forcing deviates substantially from this value. In section 4 we present the results of flume tank experiments, which support the theoretical results.
Note that the idealized considerations presented here are not intended for use in operational drift forecasting, but rather to gain fundamental physical understanding of how wind and ocean current drags influence the drift of floating objects. To accurately predict the actual trajectory of drifters, a framework is needed that accounts for a number of other pertinent factors, including inertial effects, lift and added mass forces, and the Coriolis force. Such a comprehensive account is given by the Maxey–Riley set for surface ocean inertial particles, and a series of recent studies have made substantial progress toward accurately representing the drift of floating objects in real-world settings (e.g., Beron-Vera et al. 2019; Olascoaga et al. 2020; Miron et al. 2020).
2. Balance between air and water drag
In contrast to previous leeway modeling studies, here we explicitly distinguish between two types of drag force, namely, skin friction drag and form drag, both of which can be relevant for floating objects in typical geophysical settings. This is in part motivated by the use of separate drag terms for skin and form drag in iceberg models (e.g., Gladstone et al. 2001; Martin and Adcroft 2010; Marsh et al. 2015). Skin drag arises from the object’s surface (its “skin”) being subject to shear stresses as the object moves through a fluid. This friction effect is due to laminar or turbulent flow in the boundary layer close to the object’s surface. Form drag, on the other hand, is determined by an object’s size and shape. It is proportional to the cross-sectional area of the object normal to the direction of flow of the fluid. Skin and form drag follow equivalent drag equations, scaling with the square of the relative velocity between fluid and object, but skin drag coefficients are typically three orders of magnitude smaller than form drag coefficients.
Membrane-like objects, such as oil slicks, are predominantly driven by skin drag. Form drag becomes dominant for objects with smaller length-to-height aspect ratios. Deeper keels will naturally give more importance to the role of water currents, whereas larger sails increase the sensitivity of an object’s drift to wind forcing. Skin and form drag coefficients, density, and aspect ratio are thus central characteristics that will determine an object’s sensitivity to winds and water currents.
We note that in real-world settings the drag balance is complicated by the existence of two distinct types of turbulent boundary layers, both above and below the water surface: boundary layers associated with the air–water interface and boundary layers associated with the object–fluid interfaces. Identifying the appropriate free-stream fluid velocities for the drag equations in a setting with such complex boundary layer structures is nontrivial. For simplicity, here we assume that the dominant turbulent boundary layer above and below the floating object is that associated with the fluid–object interfaces. We thus ignore the role of the boundary layers associated with the air–water interface. The free-stream fluid velocities are referenced at a height that is greater than the combined height of the sail of the object plus the thickness of the object–air boundary layer, or below the depth of the keel plus the thickness of the object–water boundary layer.
In general, it may not be straightforward to relate the free-stream wind velocity considered here to the commonly used 10-m surface wind velocity. For large objects such as icebergs with freeboard greater than 10 m, the free-stream wind velocity should naturally be measured at some height above the freeboard. For objects with small freeboards, on the other hand, the relevant free-stream velocity should be considered at a height just outside the surface boundary layer—potentially much lower than 10 m. While this free-stream velocity may in certain settings differ from the 10-m surface wind velocity, in many geophysical scenarios the two values are likely comparable, particularly when both values are measured outside the turbulent surface boundary layer.
We finally note that the exact structures of these turbulent boundary layers—particularly the structure of the near-surface ocean layer—remain a subject of ongoing research.
In what follows, we limit our analysis to one dimension. In some settings, in particular for sea ice, it has been observed that the floating object typically drifts at a turning angle relative to the direction of the surface wind, due to the Coriolis force (Nansen 1902; McPhee 2008). For sea ice, this angle is between 0° and 40° to the right of the wind in the Northern Hemisphere (Leppäranta 2011). For oil slicks, Stolzenbach et al. (1977) report turning angles smaller than 10°. In the cases of icebergs and mangrove drifters, turning angles are not often discussed; Garrett et al. (1985) argue that this is due to the draft of icebergs being deep enough to not be significantly affected by the Ekman spiral. In the search-and-rescue literature, the deviation of an object’s drift from the downwind direction is referred to as “leeway divergence,” which is not regarded to be primarily due to the Coriolis force, and ranges from −30° to beyond +30° (see, e.g., Fig. 1 in Breivik et al. 2011). Here, for simplicity we focus on the drift speed in a one-dimensional framework, and we do not address the issue of turning angles.
Note that, using the expressions (5) for freeboard and draft, we could express γ in (7) in terms of the four drag coefficients and three dimensionless parameters: the density ratio of air and water ρa/ρw, the length-to-height aspect ratio l/h, and the density of the object relative to that of water ρ/ρw.
In the following we consider two limits, that of dominant skin drag (where l ≫ h), and that of dominant form drag (where l ≲ h).
a. The limit of dominant skin drag
For many typical materials the air and water skin drag coefficients are approximately equal,
b. The limit of dominant form drag
Form drag coefficients are typically order 1, whereas skin drag coefficients are typically order 10−3. This implies that horizontal skin drag is the dominant force only when l > 103 h. In cases where l ≪ 103 h, which applies for example to most icebergs, mangrove drifters, and flotsam, the force balance on an object is predominantly a balance of the water and air form drags. Sea ice presents an intermediate case: while young, very thin, or uniformly grown sea ice (e.g., nilas) is predominantly influenced by skin drag, broken-up or deformed ice floes (e.g., featuring rafting or pressure ridges) are subject to substantial form drag (e.g., Arya 1975; Lu et al. 2011).
We refer here to an object in this limit as a vertical “sliver,” describing an (unphysical) object of finite height h and zero along-flow length l. In Fig. 1, the sliver limit is indicated by a green dashed curve, showing that this limit of (9) is essentially indistinguishable from the full solution of (7) (red solid curve) for l/h ≤ 1 when ρ/ρw > 10−2. The full solution shown in Fig. 1 is that for a square block with l/h = 1. For aspect ratios smaller than 1 this full solution will approximate the dashed sliver solution even more closely. Even for l/h = 10, the sliver is still a good approximation as long as ρ/ρw > 10−1 (not shown). This encompasses most floating solid objects in geophysical settings. Since the air and water form drag coefficients are both order 1 and show up as a ratio under a square root (which brings the value closer to unity) in the expression for γ, we make the approximation
Next, we consider two extremes for the object density: a light object where ρ ≪ ρw, and a heavy object where ρ ≥ ρw. If ρ ≥ ρw, the object is fully submerged or sinks, and the wind factor is γ = 0. More generally, as the object density increases toward the limit ρ/ρw → 1, the wind factor scales as
Turning to the light object limit, if ρ < ρa then the object would become airborne. As the object density decreases toward the limit ρ → ρa the wind factor reduces to
Next, we consider the low-density limit for the particular case of the square block with l/h = 1 (solid red curve). Since the form drag coefficients are three orders of magnitude larger than skin drag coefficients, aspect ratios of order 1 lead to a force balance that is dominated by form drag for a large range of densities. Only when
Last, for a heavy block with ρ = ρw traveling at the water surface, the force balance is between air skin drag and water form drag, giving
3. Geophysical parameter range
The higher the aspect ratio l/h for an object, the less important the object’s density becomes for determining the wind factor γ. In the membrane limit, γ is independent of ρ/ρw (Fig. 1), since the dominant force balance is between the above-water and below-water skin drags. The surfaces on which the skin stress acts do not change with density ρ, since the density only impacts how high in the water the object floats. On the other hand, for objects with small aspect ratios where form drags dominate the force balance, ρ determines the relative above-water and underwater surface areas that the form stresses act on.
Oil slicks and other thin membrane-like floating objects are therefore subject to a wind factor γ ≈ 3%, irrespective of their density. On the other hand, the wind factor of blocky icebergs and broken-up or deformed sea ice floes vary with the density of the ice. Measured sea ice densities are typically in the range from 740 to 917 kg m−3 (Timco and Frederking 1996). Using (9) and assuming
Naturally, objects that have small aspect ratios but are much lighter than sea ice will have substantially higher wind factors. For example, Fig. 2 illustrates the wind-factor dependence on l/h for an object with the density of cork, ρ = 250 kg m−3 (blue line). Here, the wind factor is γ = 7.4% for l/h = 1. In this case, the freeboard is b = 0.25h if the cork is floating in freshwater. Nevertheless, wind factors for objects investigated in search-and-rescue studies—ranging from persons in water to unballasted life rafts and drifting fishing vessels—are typically below 4% (Fig. 1 in Breivik et al. 2011).
4. Flume experiments
We performed a set of experiments to assess the validity of the theoretical considerations above in a laboratory setting. The experimental setup (see Fig. 3a) consists of a custom-built 4-m horizontal flume and a 20-in. household floor fan with adjustable speed located at the University of North Carolina Wilmington. The flume is filled with freshwater and has an adjustable water flow speed of up to 20 cm s−1. Beyond this upper limit significant ripples occur at the water surface. The flume speed is close to uniform for the upper 20 cm of the water column, and it can be set with an accuracy of about 1 cm s−1. The fan is positioned to blow in the opposite direction to the flume current. Styrofoam blocks of different geometries were sealed with impermeable tape, and we cut out openings for laboratory weights (Fig. 3a, inset). Weights of different masses were added to vary the density of the blocks.
While the flume was set to produce a known flow speed υw, the wind speed υa at the equilibrium location of the Styrofoam block was measured by hand, using a Kestrel 1000 wind anemometer at a height z = 10 cm above the water surface and taking 10-s average values. These readings were not sensitive to wind measurement height in the range z ≈ 1–20 cm above the water surface. This method was made possible by the airspeed varying considerably with distance from the fan: it would increase over a span of approximately 80 cm from 0 to almost 2.5 m s−1 (for the fan at its lowest setting), and then decrease gradually with increasing distance from the fan (see Fig. 3b). While the figure shows that the profile of increasing air velocities presented good agreement with a sinusoidal fit, the decreasing profile was less regular, and it was notably impacted by the wall near the end of the flume. The more reliable readings were therefore obtained by determining the location of the unstable equilibrium point closer to the fan (unfilled black circle in Fig. 3b), rather than the stable equilibrium point farther from the fan (solid black circle). This was done by placing the Styrofoam blocks repeatedly in the flume and observing whether they would eventually be advected by the water current to the right (when Fa < Fw), or by the wind to the left (when Fa > Fw). An approximate equilibrium location was readily found using this scheme.
We note that more sophisticated experimental setups, such as that of Miron et al. (2020), use a water flume–wind tunnel setup that features constant air and water speeds throughout. In such cases, the force balance is determined from the (approximately) constant velocity of the floating object, rather than considering a cancelling of the drag forces as we have done here.
The Styrofoam blocks had a thickness of h = 5 cm, lengths l = 5 or 10 cm, and widths w = 5, 10, or 20 cm. Here l is the along-flow dimension of the block, and w is the across-flow dimension. The different values of w were used to test the assumption that the across-flow dimension does not significantly impact the wind factor in this setting. Adding laboratory weights to the blocks allowed us to vary the densities in the range 0.03 < ρ/ρw < 0.91.
Since the shapes of the Styrofoam blocks were restricted to aspect ratios l/h of 1 and 2, these experiments are only suitable to test the validity of the theory in the sliver limit (where form drag dominates). The experimental results show that the wind factor was approximately insensitive to which aspect ratio was used (Fig. 4), consistent with the theoretical results for the sliver limit. Similarly, no discernible dependence of γ on the across-flow dimension w was observed. Two flume settings were used: a faster flow of 11 cm s−1 mostly for lighter blocks (since these were more sensitive to the fan’s wind forcing), and a slower flow of 6 cm s−1 for heavier blocks (otherwise the fan’s peak wind forcing was not sufficient to balance the water drag). The observed wind speeds at the location of stationary Styrofoam blocks fell in the range 0.5 < υa < 4.3 m s−1. This gives a span of wind factors from γ = 1.4% for the densest block to γ = 22% for the lightest block.
For ρ/ρw > 0.4, the observed wind factors deviate notably from the theoretical “light sliver” limit and fall closely on the curve of the general sliver limit. This highlights that in this limit the force balance truly is between the below-water and above-water form drags.
In section 2, we approximated that
The experimental setup imposes several limitations on the parameter range that can feasibly be explored: (i) The wind speed has an upper limit of ∼4.5 m s−1 and the flume speed has a lower limit of ∼5 cm s−1. According to the theory, this gives a lower limit for the experimental wind factor of γ ∼ 1.1%. This also implies an upper limit for the object density of ρ/ρw = 0.91 (which is roughly that of pure ice). For densities greater than that, the maximum fan strength would not be sufficient to balance the water drag, even at the lowest flume setting. (ii) To explore the low-density regime, Styrofoam blocks were hollowed out such that the effective density was as low as ρ/ρw = 0.03. It was found not to be experimentally feasible to find steady state positions for blocks with densities lower than that. (iii) To obtain observable deviations from the sliver limit (i.e., where skin drag becomes important), we estimate that the aspect ratio has to be of the order
Previous flume experiments with mangrove propagules
Van der Stocken et al. (2015) report on flume experiments with a setup similar to the one presented above, investigating the drift of different species of mangrove propagules. The authors carried out a set of 16 different experiments for each species with either no wind or one of three different wind velocities (low, 2.7 m s−1; medium, 4.5 m s−1; high, 6.0 m s−1), three different water velocities (low, 0 m s−1; medium 0.15 m s−1; high, 0.3 m s−1), and two wind directions (along the water flow and against the water flow). The densities of the propagules fall in the range ρ/ρw = 0.5–1. The shapes for the different species vary widely, from the relatively light H. littoralis whose propagules resemble “small sailboats” to the dense spherical X. granatum, which is known as the “cannonball” (Van der Stocken et al. 2015).
We compute the wind factor from the reported velocities as γ(ρ) = [υ(ρ) − υw]/υa, where υ(ρ) was digitized from Fig. 2 in Van der Stocken et al. (2015). One may have expected that γ would extend over a large range for this widely varied set of experiments. However, even for the lightest propagules we find γ < 5% (Fig. 5). For most of the propagules lighter than ρ/ρw = 0.85 we find γ = 1.5%–3.5%, whereas for ρ/ρw > 0.85 we typically find γ = 0.5%–2.5%. The mangrove propagules are thus subject to wind factors that are broadly similar to those for ice and oil.
One might further expect a large spread in the ratio of drag coefficients due to the very different shapes of the propagules. Yet, for most species the best fit to the theoretically computed γ is obtained using
5. Conclusions
We found that in many geophysical settings, the drift of floating objects is predominantly determined by a balance of near-surface winds and water currents, and the sensitivity to wind forcing falls in the relatively narrow range of γ ≃ 2%–4%. We explained the behavior using an analytical solution of an approximate force balance and showed that this solution holds for sea ice, icebergs, oil slicks, mangrove propagules, and a wide range of man-made floating objects. The behavior is similar despite the widely varying densities, shapes, and sizes of these floating objects. We showed that this dependence essentially occurs due to the wind sensitivity being approximately equal to (ρa/ρw)1/2 ≈ 3%. The different characteristics of these different objects typically only result in a small correction to this number. The drift of most objects is set by a balance of form drags, with skin drag playing a negligible role, except for the membrane-like oil slicks and to a certain degree sea ice. The flume experiments presented here, as well as those with mangrove drifters by Van der Stocken et al. (2015), can serve two practical purposes: (i) a straightforward way to establish the wind sensitivity of a given type of floating object and (ii) a low-tech approach to estimate the form drag ratio
Acknowledgments.
This study was supported by NSF OPP Grants 1643445 and 1744835. Author Constantinou was supported by the Australian Research Council DECRA Fellowship DE210100749. Without implying their endorsement, we thank Emma Beer, Roger Samelson, Juan Restrepo, and Tapio Schneider for helpful conversations about this project.
Data availability statement.
The experimental data are available online (https://github.com/tillwagner/drift).
REFERENCES
Allen, A. A., 2005: Leeway divergence. U.S. Coast Guard Research and Development Center Tech. Rep. CG-D-05-05, 128 pp.
Allen, A. A., and J. V. Plourde, 1999: Review of leeway: field experiments and implementation. U.S. Coast Guard Research and Development Center Tech. Rep. CG-D-08-99, 351 pp.
Arya, S., 1975: A drag partition theory for determining the large-scale roughness parameter and wind stress on the Arctic pack ice. J. Geophys. Res., 80, 3447–3454, https://doi.org/10.1029/JC080i024p03447.
Beron-Vera, F. J., M. J. Olascoaga, and P. Miron, 2019: Building a Maxey–Riley framework for surface ocean inertial particle dynamics. Phys. Fluids, 31, 096602, https://doi.org/10.1063/1.5110731.
Bigg, G. R., M. R. Wadley, D. P. Stevens, and J. A. Johnson, 1996: Prediction of iceberg trajectories for the north Atlantic and Arctic Oceans. Geophys. Res. Lett., 23, 3587–3590, https://doi.org/10.1029/96GL03369.
Bigg, G. R., M. R. Wadley, D. P. Stevens, and J. A. Johnson, 1997: Modelling the dynamics and thermodynamics of icebergs. Cold Reg. Sci. Technol., 26, 113–135, https://doi.org/10.1016/S0165-232X(97)00012-8.
Breivik, Ø., A. A. Allen, C. Maisondieu, and J. C. Roth, 2011: Wind-induced drift of objects at sea: The leeway field method. Appl. Ocean Res., 33, 100–109, https://doi.org/10.1016/j.apor.2011.01.005.
Breivik, Ø., A. A. Allen, C. Maisondieu, J.-C. Roth, and B. Forest, 2012: The leeway of shipping containers at different immersion levels. Ocean Dyn., 62, 741–752, https://doi.org/10.1007/s10236-012-0522-z.
Breivik, Ø., A. A. Allenb, C. Maisondieuc, and M. Olagnonc, 2013: Advances in search and rescue at sea. Ocean Dyn., 63, 83–88, https://doi.org/10.1007/s10236-012-0581-1.
Browne, I., and A. Crary, 1958: The movement of ice in the Arctic Ocean. Arctic Sea Ice, Publ. 598, National Academy of Sciences, 191–207.
Daniel, P., G. Jan, F. Cabiocâh, Y. Landau, and E. Loiseau, 2002: Drift modeling of cargo containers. Spill Sci. Technol. Bull., 7, 279–288, https://doi.org/10.1016/S1353-2561(02)00075-0.
Garrett, C., J. Middleton, M. Hazen, and F. Majaess, 1985: Tidal currents and eddy statistics from iceberg trajectories off Labrador. Science, 227, 1333–1335, https://doi.org/10.1126/science.227.4692.1333.
Gladstone, R. M., G. R. Bigg, and K. W. Nicholls, 2001: Iceberg trajectory modeling and meltwater injection in the southern ocean. J. Geophys. Res. Oceans, 106, 19 903–19 915, https://doi.org/10.1029/2000JC000347.
Keghouche, I., L. Bertino, and K. A. Lisæter, 2009: Parameterization of an iceberg drift model in the Barents sea. J. Atmos. Oceanic Technol., 26, 2216–2227, https://doi.org/10.1175/2009JTECHO678.1.
Kubat, I., M. Sayed, S. B. Savage, and T. Carrieres, 2005: An operational model of iceberg drift. Int. J. Offshore Polar Eng., 15, 125–131.
Leppäranta, M., 2011: The Drift of Sea Ice. Springer, 350 pp.
Lu, P., Z. Li, B. Cheng, and M. Leppäranta, 2011: A parameterization of the ice-ocean drag coefficient. J. Geophys. Res. Oceans, 116, C07019, https://doi.org/10.1029/2010JC006878.
Marsh, R., and Coauthors, 2015: NEMO–ICB (v1. 0): Interactive icebergs in the NEMO ocean model globally configured at eddy-permitting resolution. Geosci. Model Dev., 8, 1547–1562, https://doi.org/10.5194/gmd-8-1547-2015.
Martin, T., and A. Adcroft, 2010: Parameterizing the fresh-water flux from land ice to ocean with interactive icebergs in a coupled climate model. Ocean Modell., 34, 111–124, https://doi.org/10.1016/j.ocemod.2010.05.001.
McPhee, M., 2008: Air-Ice-Ocean Interaction: Turbulent Ocean Boundary Layer Exchange Processes. Springer, 215 pp.
Miron, P., S. Medina, M. Olascoaga, and F. Beron-Vera, 2020: Laboratory verification of the buoyancy dependence of the carrying flow in a Maxey–Riley theory for inertial ocean dynamics. Phys. Fluids, 32, 071703, https://doi.org/10.1063/5.0018272.
Nansen, F., 1902: The Oceanography of the North Polar Basin. Vol. 3, The Norwegian North Polar Expedition 1893–1896: Scientific Results, Longmans, Green and Co., 427 pp.
Nesterov, O., 2018: Consideration of various aspects in a drift study of MH370 debris. Ocean Sci., 14, 387–402, https://doi.org/10.5194/os-14-387-2018.
Olascoaga, M. J., F. J. Beron-Vera, P. Miron, J. Triñanes, N. Putman, R. Lumpkin, and G. Goni, 2020: Observation and quantification of inertial effects on the drift of floating objects at the ocean surface. Phys. Fluids, 32, 026601, https://doi.org/10.1063/1.5139045.
Rasmussen, D., 1985: Oil spill modeling—A tool for cleanup operations. International Oil Spill Conference Proceedings (1985), American Petroleum Institute, 243–249, https://doi.org/10.7901/2169-3358-1985-1-243.
Röhrs, J., K. H. Christensen, L. R. Hole, G. Broström, M. Drivdal, and S. Sundby, 2012: Observation-based evaluation of surface wave effects on currents and trajectory forecasts. Ocean Dyn., 62, 1519–1533, https://doi.org/10.1007/s10236-012-0576-y.
Samelson, R., 2020: Turbulent universality and the drift velocity at the interface between two homogeneous fluids. Phys. Fluids, 32, 082107, https://doi.org/10.1063/5.0019733.
Smith, S. D., and E. G. Banke, 1983: The influence of winds, currents and towing forces on the drift of icebergs. Cold Reg. Sci. Technol., 6, 241–255, https://doi.org/10.1016/0165-232X(83)90045-9.
Stolzenbach, K. D., O. S. Madsen, E. E. Adams, A. M. Pollack, and C. Cooper, 1977: A review and evaluation of basic techniques for predicting the behavior of surface oil slicks. NOAA Tech. Rep. MITSG 77-8, 324 pp., https://repository.library.noaa.gov/view/noaa/9623.
Sutherland, G., and Coauthors, 2020: Evaluating the leeway coefficient of ocean drifters using operational marine environmental prediction systems. J. Atmos. Oceanic Technol., 37, 1943–1954, https://doi.org/10.1175/JTECH-D-20-0013.1.
Thorndike, A., and R. Colony, 1982: Sea ice motion in response to geostrophic winds. J. Geophys. Res. Oceans, 87, 5845–5852, https://doi.org/10.1029/JC087iC08p05845.
Timco, G., and R. Frederking, 1996: A review of sea ice density. Cold Reg. Sci. Technol., 24, 1–6, https://doi.org/10.1016/0165-232X(95)00007-X.
Van der Stocken, T., B. Vanschoenwinkel, D. J. De Ryck, T. J. Bouma, F. Dahdouh-Guebas, and N. Koedam, 2015: Interaction between water and wind as a driver of passive dispersal in mangroves. PLOS ONE, 10, e0121593, https://doi.org/10.1371/journal.pone.0121593.
Wagner, T. J., R. W. Dell, and I. Eisenman, 2017: An analytical model of iceberg drift. J. Phys. Oceanogr., 47, 1605–1616, https://doi.org/10.1175/JPO-D-16-0262.1.
Zubov, N., 1945: L’dy Arktiki [Ice of the Arctic]. Izd-vo Glavsevmorputi, 491 pp.