How Winds and Ocean Currents Influence the Drift of Floating Objects

Till J. W. Wagner aUniversity of Wisconsin–Madison, Madison, Wisconsin
bUniversity of North Carolina Wilmington, Wilmington, North Carolina

Search for other papers by Till J. W. Wagner in
Current site
Google Scholar
PubMed
Close
https://orcid.org/0000-0003-4572-1285
,
Ian Eisenman cScripps Institution of Oceanography, University of California, San Diego, La Jolla, California

Search for other papers by Ian Eisenman in
Current site
Google Scholar
PubMed
Close
,
Amanda M. Ceroli bUniversity of North Carolina Wilmington, Wilmington, North Carolina
dOcean Sciences Division, U.S. Naval Research Laboratory, Stennis Space Center, Mississippi

Search for other papers by Amanda M. Ceroli in
Current site
Google Scholar
PubMed
Close
, and
Navid C. Constantinou eAustralian National University, Canberra, Australian Capital Territory, Australia

Search for other papers by Navid C. Constantinou in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

Arctic icebergs, unconstrained sea ice floes, oil slicks, mangrove drifters, lost cargo containers, and other flotsam are known to move at 2%–4% of the prevailing wind velocity relative to the water, despite vast differences in the material properties, shapes, and sizes of objects. Here, we revisit the roles of density, aspect ratio, and skin and form drag in determining how an object is driven by winds and water currents. Idealized theoretical considerations show that although substantial differences exist for end members of the parameter space (either very thin or thick and very light or dense objects), most realistic cases of floating objects drift at approximately 3% of the free-stream wind velocity (measured outside an object’s surface boundary layer) relative to the water. This relationship, known as a long-standing rule of thumb for the drift of various types of floating objects, arises from the square root of the ratio of the density of air to that of water. We support our theoretical findings with flume experiments using floating objects with a range of densities and shapes.

© 2022 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Till Wagner, till.wagner@wisc.edu

Abstract

Arctic icebergs, unconstrained sea ice floes, oil slicks, mangrove drifters, lost cargo containers, and other flotsam are known to move at 2%–4% of the prevailing wind velocity relative to the water, despite vast differences in the material properties, shapes, and sizes of objects. Here, we revisit the roles of density, aspect ratio, and skin and form drag in determining how an object is driven by winds and water currents. Idealized theoretical considerations show that although substantial differences exist for end members of the parameter space (either very thin or thick and very light or dense objects), most realistic cases of floating objects drift at approximately 3% of the free-stream wind velocity (measured outside an object’s surface boundary layer) relative to the water. This relationship, known as a long-standing rule of thumb for the drift of various types of floating objects, arises from the square root of the ratio of the density of air to that of water. We support our theoretical findings with flume experiments using floating objects with a range of densities and shapes.

© 2022 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Till Wagner, till.wagner@wisc.edu

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.

Following Martin and Adcroft (2010), we write the drag force on the floating object due to water as
Fw=12(CwFAwV+CwSAwH)ρw|Δvw|Δvw,
where Δvwvwv denotes the velocity of the object v relative to the free-stream water velocity vw, CwF and CwS are the form and skin drag coefficients for water, and AwV and AwH are the cross-sectional vertical and tangential horizontal surface areas of the object facing the relative flow of the water. Similarly, the drag force due to air is
Fa=12(CaFAaV+CaSAaH)ρa|Δva|Δva,
where the definitions of all terms are as given above for (1) except that the subscripts have been switched from w for water to a for air.
In geophysical settings, other forces including the Coriolis force may also play important roles in determining the drift of floating objects (as discussed in section 1). The momentum equation for the drifting object can be written in the form
mdvdt=Fw+Fa+Fc+F,
where m is the mass of the object, Fc=mfk^×v is the Coriolis force (with f being the Coriolis parameter and k^ being the vertical unit vector), and the term F is a placeholder to represent any other forces acting on the object. The Maxey–Riley set mentioned above considers the force balance (3) for spherical particles that are immersed in unsteady and nonuniform flow (Beron-Vera et al. 2019).
In this study, we consider the limit where inertial effects, the Coriolis force, and other forces are small relative to air and water drag. In this case, the drift velocity is determined by a balance of drag forces due to near-surface water currents and winds. To provide a broad-strokes assessment of when this limit is applicable, we consider the ratio of the water drag force to the Coriolis force:
|Fw||Fc|ρwρ|Δvw|lf,
where ρ is the density of the object and l is its length. Here we have approximated that CwF1, CwS=0, m=ρAwVl, and |v| ∼ |Δvw|, and we have neglected a factor of 2. For objects floating in the ocean, commonly ρw/ρ ∼ 1 and |Δvw|0.1ms1. With a typical midlatitude value for the Coriolis parameter of f ∼ 10−4 s−1, we find that the drag force dominates the Coriolis effect when l1km. Similarly, considering the ratio of air drag to Coriolis, we have
|Fa||Fc|ρaρ|Δva|2lf|v|2.
In this case ρa/ρ ∼ 10−3 and |Δva| ∼ 10 m s−1, resulting in dominant drag when l10km. In summary, the effects due to the rotation of Earth are typically small when the size of the floating object is much less than O(1 km). For larger objects, the Coriolis force becomes important. The force balance for drift of giant tabular icebergs, which is dominated by the Coriolis and pressure gradient forces, was previously considered in a similar framework by Wagner et al. (2017).
Based on the discussion above, to determine the steady drift of small floating objects—with lO(1 km)—the air and water drag forces need to balance:
mdvdt=Fw+Fa=0.
Here, we consider idealized rectangular shapes such that
AaV=wb,AwV=wd,andAaH=AwH=wl,
where w is the object’s width, b is its freeboard (height above the waterline), d is its draft (depth below the waterline), and l is its along-flow length. The across-flow width w is the same for all surfaces and therefore cancels out in all calculations below. Henceforth, we only consider a two-dimensional framework, with vertical and along-flow dimensions. The freeboard and draft of the object can be calculated in terms of the object’s height h, its density ρ, and the density of water ρw as
b=h(1ρ/ρw)andd=hρ/ρw.

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.

Since the air velocity υa is typically much larger than the velocity of the object υ, we can approximate Δυaυaυυa (e.g., Stolzenbach et al. 1977). Substituting the drag expressions (1) and (2) into the force balance (4) and solving for υ, we find
υ=υw+γυa,
where γ is the leeway factor or wind factor and can be interpreted as the object’s sensitivity to wind forces. Here, the wind factor is defined as
γρaρwbCaF+lCaSdCwF+lCwS.
The four terms under the second square root represent the relative contributions of form and skin air drag (numerator) and form and skin water drag (denominator), respectively.

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 lh), and that of dominant form drag (where lh).

a. The limit of dominant skin drag

The limit where skin drag dominates over form drag, which we refer to as the “membrane limit,” has been investigated previously (e.g., Stolzenbach et al. 1977). In this limit, lh and thus the wind factor γ simplifies to
γρaρwCaSCwS.
If the skin drag coefficients are replaced by effective drag coefficients Ca and Cw that include form and skin drag, this expression for γ becomes equivalent to the Nansen number, Naρa/ρwCa/Cw (Leppäranta 2011).

For many typical materials the air and water skin drag coefficients are approximately equal, CaSCwS (Stolzenbach et al. 1977), in which case the wind factor (8) reduces to γ ≈ (ρa/ρw)1/2 ≈ 3% (since ρw ≈ 103 kg m−3 and ρa ≈ 1 kg m−3). This explains why sufficiently thin objects typically drift at approximately 3% of the wind speed relative to the water, regardless of the object’s density and length-to-height ratio, which do not appear in (8). This membrane limit is indicated by the solid blue horizontal line in Fig. 1. For many floating objects the membrane limit is not appropriate, but we show below that a wide range of objects nonetheless have a wind factor that is similar to 3%.

Fig. 1.
Fig. 1.

Dependence of the wind factor γ on the density ratio ρ/ρw, for different aspect ratios l/h. Shown are the membrane limit (l/h → ∞; purple horizontal line) and the case of a square block (l/h = 1; red solid curve). The general sliver limit is shown (l/h → 0; green dashed), as well as the light sliver case where ρ/ρw → 0 (dotted black line). In this limit γ scales as (ρ/ρw)−1/2, as indicated by the black-outlined triangle. The red dashed line indicates the limit of a light block where the force balance is between water skin drag and air form drag. Also shown is the limit of a heavy square block (gray star). The insets indicate the dominant force balance in each case: red and blue arrows represent water and air velocity, and red and blue vertical and horizontal lines denote the surfaces on which the dominant drag components act. Note that for ρ < ρa (gray shaded region), physical objects would float into the air.

Citation: Journal of Physical Oceanography 52, 5; 10.1175/JPO-D-20-0275.1

In a recent study, Samelson (2020) showed that the flow velocity υint at the interface of two fluids is in general given by
|υintυ2|=ρ1ρ1+ρ2|υ1υ2|,
where ρ1 is the density of the lighter fluid, ρ2 is the density of the heavier fluid, and υ1,2 are their corresponding free-stream velocities. When ρ2ρ1 (such as is the case for air and water), the interface velocity reduces to υintυ2 + (ρ1/ρ2)1/2υ1 [(20) in Samelson (2020)], which is identical to the expression for the membrane limit drift velocity above. We note that the approach of Samelson (2020) is different from the drag force balance we use here. Samelson’s result arises from a condition of stress continuity at the fluid interface, while our result derives from balancing the drag forces acting on an object located at the fluid interface. However, the stress continuity condition at the fluid interface (alongside a dimensional analysis argument that relates the surface stress to a velocity scale) implies that ρaυ˜a2=ρwυ˜w2, where υ˜a and υ˜w are the friction velocities for air and water (Samelson 2020). This condition is remarkably similar to the drag force balance in the membrane limit. Note that a similar dimensional analysis argument is also used in deriving the drag force expressions (1) and (2). This explains the similarity between our wind factor expression in the membrane limit (and ignoring any lateral boundary layer effects associated with the finite horizontal extent of the floating membrane), given by (8) and expression (20) in Samelson (2020).

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).

In the limit l → 0, the wind factor γ reduces to
γ=ρaρwbCaFdCwF=ρaρwCaFCwF1ρ/ρwρ/ρw.

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 CaF=CwF in Fig. 1. Given the assumption of equal drag coefficients, (9) is equivalent to expression (16) in Daniel et al. (2002).

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 γ = ρa/ρw(1ρ/ρw) (not shown).

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 γ=1ρa/ρw1, since ρa/ρw ≈ 10−3 ≪ 1. This is in agreement with basic physical intuition: when an object’s density approaches that of air, the object becomes no longer submerged in the water and its drift is fully determined by wind forcing (neglecting capillary effects and surface tension, which are not included in this analysis). Considering instead the ratio of the density of the object and the density of water, in the light limit ρ/ρw ≪ 1 we can approximate the last term in (9) as (ρ/ρw)−1/2. This scaling is illustrated in Fig. 1 (dotted black line). The wind factor then reduces to γ = (ρa/ρ)1/2.

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 ρ/ρw102 does the underwater skin drag become important and γ diverges notably from the sliver approximation. In the limit ρ/ρw → 0, the force balance is purely between the above-water form drag and the underwater skin drag, which gives γ=ρa/ρwCaF/CwS1 using CaF/CwS103 (red dashed line). For the physical lower bound ρ = ρa, we find γ = 0.71. That is to say, an idealized square block that is as light as air but still feels skin drag from the water will drift at 71% of the wind velocity. Therefore, even in this extreme scenario the water skin drag still plays a substantial role in determining the object’s drift.

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 γ=ρa/ρwCaS/CwF 103 (star symbol in Fig. 1).

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 CaF=CwF, we find that this corresponds to an upper bound for the wind factor for sea ice γ = 2.5%, and a lower bound γ = 1.3% for pure ice (when l/h = 1). As l/h increases, these upper and lower bounds converge to the membrane limit γ = (ρa/ρw)1/2 (Fig. 2). An estimate of typical iceberg density, accounting for the presence of a firn layer, is ρ = 850 kg m−3 (Bigg et al. 1997), which gives γ = 1.8% (Fig. 2). These theoretical values agree with observed ranges of γ = 1.6%–1.8% for icebergs (Smith and Banke 1983; Garrett et al. 1985; Bigg et al. 1997) and γ = 2.0%–2.5% for sea ice (Nansen 1902; Zubov 1945; Browne and Crary 1958; Thorndike and Colony 1982).

Fig. 2.
Fig. 2.

The wind factor γ as a function of the aspect ratio l/h, for a range of object densities ρ. In the membrane limit (l/h → ∞), the wind factor reduces to γ = (ρa/ρw)1/2 = 3.4%, independent of ρ. The lower bound for aspect ratios shown here is the square block, with l/h = 1. In this case γ is sensitive to density ρ. Values of γ for ice of different densities (see the text) are indicated, as well as the example of a material with the density of cork (250 kg m−3). The gray-shaded areas show the limits for which objects would start floating up into the air (ρ < ρa) or sinking down into the water (ρ > ρw).

Citation: Journal of Physical Oceanography 52, 5; 10.1175/JPO-D-20-0275.1

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.

Fig. 3.
Fig. 3.

(a) Experimental flume setup with fan and Styrofoam block. The opposing directions of the water flow (dark-blue arrow) and airflow (light-blue arrow) are indicated. Inset: examples of two Styrofoam blocks used, one of which has a few weights resting next to the cutout, which has been temporarily removed for the weights to be inserted into the hole. (b) Measured mean air velocity at 10 cm above the water surface as a function of distance along the flume (blue circles). The black curve represents a visual best-fit sine function, added purely for illustrative purposes. Air and water drag forces are balanced when the airspeed is equal (and opposite) to υw/γ (red line; see text). In the gray-shaded regions (where the black curve crosses below the red line) water drag is greater than air drag (Fw > Fa) and the object is advected to the right. When the black curve is above the red line, drift occurs to the left, as indicated by the black arrows. The filled black circle therefore represents a stable equilibrium location, whereas the open black circle represents an unstable equilibrium location.

Citation: Journal of Physical Oceanography 52, 5; 10.1175/JPO-D-20-0275.1

The Styrofoam blocks were then placed in the flume and would assume equilibrium resting positions at different distances from the fan where the drags due to the flume current and the fan wind balanced. From (6), we find that in this setting (for v =0) the wind factor is simply given by
γ=|υw|/|υa|.

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.

Fig. 4.
Fig. 4.

Dependence of wind factor on density for flume experiments. Red and blue symbols represent experiments with the flume set to different water speeds as indicated in the figure legend. The symbols correspond to different object geometries: l/h is 1 (diamonds and squares) or 2 (circles and triangles). The across-drag dimension of the blocks w/l is 1 (squares and circles), 2 (diamonds), or 4 (triangles). The theoretical curves are as in Fig. 1, using the best-fit ratio of form drag coefficients, CaF/CwF=1.4.

Citation: Journal of Physical Oceanography 52, 5; 10.1175/JPO-D-20-0275.1

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 CaF/CwF=1. Previous work on drag coefficients for icebergs (Kubat et al. 2005; Keghouche et al. 2009) and other objects have given values for this ratio, which tends to be similar to 1. Here, we can use the observed wind factor dependence, γ(ρ/ρw), to estimate the ratio CaF/CwF for the blocks in the flume experiments. This is done by varying CaF/CwF in the theoretical estimate for γ and computing the root-mean-square error between the resulting theoretical curve and the observed values of γ. We find that for our experimental configuration, the root-mean-square error is minimized when we take the ratio to be CaF/CwF=1.4. For this drag ratio, we observe strikingly close correspondence between the observed and theoretical estimates for the full range of block densities in the experiments (Fig. 4). Notably, this ratio for the Styrofoam blocks is in close agreement with the findings for icebergs of Bigg et al. (1996), who derive this ratio from iceberg trajectories and report CaF/CwF=1.44.

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 l/h100 (see Fig. 1). The experimental setup was not suitable for sheets of such small thickness: surface ripples and capillary waves did not allow for top surfaces that were consistently splash free or bottom surfaces that did not feature bubbles, thus distorting the drag balance considerably. Hence, we were not able to approach the membrane limit in the flume experiments.

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.

Fig. 5.
Fig. 5.

Previous flume experiment results derived from data digitized from Van der Stocken et al. (2015). Each black marker corresponds to one species of mangrove. Error bars are ±1σ for all measurements for a given species. The different colors of the small symbols correspond to different wind strengths in the flume experiments: low 2.7 (blue), medium 4.5 (green), and high 6.0 (blue) m s−1. The different symbols correspond to different water velocities: low 0 (circle), medium 0.15 (square), and high 0.3 (triangle) m s−1. Gray-filled symbols denote experiments in which the wind direction was opposite the water flow; white-filled symbols denote experiments in which it was aligned with the water flow. For most species, the theoretical solution with CaF/CwF=1.7 (red curve) yields a relatively close fit.

Citation: Journal of Physical Oceanography 52, 5; 10.1175/JPO-D-20-0275.1

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 CaF/CwF1.7 (Fig. 5). Outliers are the lightest propagules and the densest ones. We note that for the densest species, the surface area and freeboard exposed to the wind becomes very small, on the order of the capillary length of water (∼3 mm), which means that effects such as surface tension and wind-induced ripples in the water surface are expected to play a nonnegligible role. The length-to-height aspect ratios range from vertically drifting, sliver-like C. tagal with l/h ≪ 1 to thin, flat R. mucronata with l/h ≈ 30 (see Van der Stocken et al. 2015, their Fig. S1A). All of these fall nevertheless near the limit of dominant form drag in the theory presented above.

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 CaF/CwF, under the assumptions of these theoretical calculations. The theoretical considerations presented here provide a physically intuitive framework to parse out the relative importance of water currents and winds in determining the drift of floating objects.

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, 34473454, https://doi.org/10.1029/JC080i024p03447.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 35873590, https://doi.org/10.1029/96GL03369.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 113135, https://doi.org/10.1016/S0165-232X(97)00012-8.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 100109, https://doi.org/10.1016/j.apor.2011.01.005.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 741752, https://doi.org/10.1007/s10236-012-0522-z.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Breivik, Ø., A. A. Allenb, C. Maisondieuc, and M. Olagnonc, 2013: Advances in search and rescue at sea. Ocean Dyn., 63, 8388, https://doi.org/10.1007/s10236-012-0581-1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 279288, https://doi.org/10.1016/S1353-2561(02)00075-0.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Garrett, C., J. Middleton, M. Hazen, and F. Majaess, 1985: Tidal currents and eddy statistics from iceberg trajectories off Labrador. Science, 227, 13331335, https://doi.org/10.1126/science.227.4692.1333.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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 90319 915, https://doi.org/10.1029/2000JC000347.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 22162227, https://doi.org/10.1175/2009JTECHO678.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kubat, I., M. Sayed, S. B. Savage, and T. Carrieres, 2005: An operational model of iceberg drift. Int. J. Offshore Polar Eng., 15, 125131.

    • Search Google Scholar
    • Export Citation
  • Leppäranta, M., 2011: The Drift of Sea Ice. Springer, 350 pp.

    • Crossref
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 15471562, https://doi.org/10.5194/gmd-8-1547-2015.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 111124, https://doi.org/10.1016/j.ocemod.2010.05.001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McPhee, M., 2008: Air-Ice-Ocean Interaction: Turbulent Ocean Boundary Layer Exchange Processes. Springer, 215 pp.

    • Crossref
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 387402, https://doi.org/10.5194/os-14-387-2018.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Export Citation
  • 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, 15191533, https://doi.org/10.1007/s10236-012-0576-y.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 241255, https://doi.org/10.1016/0165-232X(83)90045-9.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 19431954, https://doi.org/10.1175/JTECH-D-20-0013.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thorndike, A., and R. Colony, 1982: Sea ice motion in response to geostrophic winds. J. Geophys. Res. Oceans, 87, 58455852, https://doi.org/10.1029/JC087iC08p05845.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Timco, G., and R. Frederking, 1996: A review of sea ice density. Cold Reg. Sci. Technol., 24, 16, https://doi.org/10.1016/0165-232X(95)00007-X.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wagner, T. J., R. W. Dell, and I. Eisenman, 2017: An analytical model of iceberg drift. J. Phys. Oceanogr., 47, 16051616, https://doi.org/10.1175/JPO-D-16-0262.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zubov, N., 1945: L’dy Arktiki [Ice of the Arctic]. Izd-vo Glavsevmorputi, 491 pp.

Save
  • 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, 34473454, https://doi.org/10.1029/JC080i024p03447.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 35873590, https://doi.org/10.1029/96GL03369.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 113135, https://doi.org/10.1016/S0165-232X(97)00012-8.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 100109, https://doi.org/10.1016/j.apor.2011.01.005.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 741752, https://doi.org/10.1007/s10236-012-0522-z.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Breivik, Ø., A. A. Allenb, C. Maisondieuc, and M. Olagnonc, 2013: Advances in search and rescue at sea. Ocean Dyn., 63, 8388, https://doi.org/10.1007/s10236-012-0581-1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 279288, https://doi.org/10.1016/S1353-2561(02)00075-0.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Garrett, C., J. Middleton, M. Hazen, and F. Majaess, 1985: Tidal currents and eddy statistics from iceberg trajectories off Labrador. Science, 227, 13331335, https://doi.org/10.1126/science.227.4692.1333.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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 90319 915, https://doi.org/10.1029/2000JC000347.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 22162227, https://doi.org/10.1175/2009JTECHO678.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Kubat, I., M. Sayed, S. B. Savage, and T. Carrieres, 2005: An operational model of iceberg drift. Int. J. Offshore Polar Eng., 15, 125131.

    • Search Google Scholar
    • Export Citation
  • Leppäranta, M., 2011: The Drift of Sea Ice. Springer, 350 pp.

    • Crossref
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 15471562, https://doi.org/10.5194/gmd-8-1547-2015.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 111124, https://doi.org/10.1016/j.ocemod.2010.05.001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • McPhee, M., 2008: Air-Ice-Ocean Interaction: Turbulent Ocean Boundary Layer Exchange Processes. Springer, 215 pp.

    • Crossref
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 387402, https://doi.org/10.5194/os-14-387-2018.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Export Citation
  • 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, 15191533, https://doi.org/10.1007/s10236-012-0576-y.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 241255, https://doi.org/10.1016/0165-232X(83)90045-9.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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, 19431954, https://doi.org/10.1175/JTECH-D-20-0013.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Thorndike, A., and R. Colony, 1982: Sea ice motion in response to geostrophic winds. J. Geophys. Res. Oceans, 87, 58455852, https://doi.org/10.1029/JC087iC08p05845.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Timco, G., and R. Frederking, 1996: A review of sea ice density. Cold Reg. Sci. Technol., 24, 16, https://doi.org/10.1016/0165-232X(95)00007-X.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • 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.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wagner, T. J., R. W. Dell, and I. Eisenman, 2017: An analytical model of iceberg drift. J. Phys. Oceanogr., 47, 16051616, https://doi.org/10.1175/JPO-D-16-0262.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zubov, N., 1945: L’dy Arktiki [Ice of the Arctic]. Izd-vo Glavsevmorputi, 491 pp.

  • Fig. 1.

    Dependence of the wind factor γ on the density ratio ρ/ρw, for different aspect ratios l/h. Shown are the membrane limit (l/h → ∞; purple horizontal line) and the case of a square block (l/h = 1; red solid curve). The general sliver limit is shown (l/h → 0; green dashed), as well as the light sliver case where ρ/ρw → 0 (dotted black line). In this limit γ scales as (ρ/ρw)−1/2, as indicated by the black-outlined triangle. The red dashed line indicates the limit of a light block where the force balance is between water skin drag and air form drag. Also shown is the limit of a heavy square block (gray star). The insets indicate the dominant force balance in each case: red and blue arrows represent water and air velocity, and red and blue vertical and horizontal lines denote the surfaces on which the dominant drag components act. Note that for ρ < ρa (gray shaded region), physical objects would float into the air.

  • Fig. 2.

    The wind factor γ as a function of the aspect ratio l/h, for a range of object densities ρ. In the membrane limit (l/h → ∞), the wind factor reduces to γ = (ρa/ρw)1/2 = 3.4%, independent of ρ. The lower bound for aspect ratios shown here is the square block, with l/h = 1. In this case γ is sensitive to density ρ. Values of γ for ice of different densities (see the text) are indicated, as well as the example of a material with the density of cork (250 kg m−3). The gray-shaded areas show the limits for which objects would start floating up into the air (ρ < ρa) or sinking down into the water (ρ > ρw).

  • Fig. 3.

    (a) Experimental flume setup with fan and Styrofoam block. The opposing directions of the water flow (dark-blue arrow) and airflow (light-blue arrow) are indicated. Inset: examples of two Styrofoam blocks used, one of which has a few weights resting next to the cutout, which has been temporarily removed for the weights to be inserted into the hole. (b) Measured mean air velocity at 10 cm above the water surface as a function of distance along the flume (blue circles). The black curve represents a visual best-fit sine function, added purely for illustrative purposes. Air and water drag forces are balanced when the airspeed is equal (and opposite) to υw/γ (red line; see text). In the gray-shaded regions (where the black curve crosses below the red line) water drag is greater than air drag (Fw > Fa) and the object is advected to the right. When the black curve is above the red line, drift occurs to the left, as indicated by the black arrows. The filled black circle therefore represents a stable equilibrium location, whereas the open black circle represents an unstable equilibrium location.

  • Fig. 4.

    Dependence of wind factor on density for flume experiments. Red and blue symbols represent experiments with the flume set to different water speeds as indicated in the figure legend. The symbols correspond to different object geometries: l/h is 1 (diamonds and squares) or 2 (circles and triangles). The across-drag dimension of the blocks w/l is 1 (squares and circles), 2 (diamonds), or 4 (triangles). The theoretical curves are as in Fig. 1, using the best-fit ratio of form drag coefficients, CaF/CwF=1.4.

  • Fig. 5.

    Previous flume experiment results derived from data digitized from Van der Stocken et al. (2015). Each black marker corresponds to one species of mangrove. Error bars are ±1σ for all measurements for a given species. The different colors of the small symbols correspond to different wind strengths in the flume experiments: low 2.7 (blue), medium 4.5 (green), and high 6.0 (blue) m s−1. The different symbols correspond to different water velocities: low 0 (circle), medium 0.15 (square), and high 0.3 (triangle) m s−1. Gray-filled symbols denote experiments in which the wind direction was opposite the water flow; white-filled symbols denote experiments in which it was aligned with the water flow. For most species, the theoretical solution with CaF/CwF=1.7 (red curve) yields a relatively close fit.

All Time Past Year Past 30 Days
Abstract Views 734 0 0
Full Text Views 1927 1222 129
PDF Downloads 1419 896 90