Six “lab-scale” model simulations (L1–L6) are completed to facilitate an examination of velocity fields near the ice boundary. The simulations feature two different monochromatic wave types (col. 2 and 3) and three different ice drag settings (col. 6). For all of these tests, modeled ice thickness (col. 4) is fixed at 1 cm, and tank depth (col. 5) is set to 1.5 m. These six scenarios include a combination of high resolution (for resolving waves, boundary layers and ice thickness) and moderate to large horizontal extent (to avoid numerical effects near the wavemaker and ice edges). The horizontal and vertical grid-spacing of these smaller-scale simulations are cm and mm, respectively, the tank length, , is 26.2 m, and output is sampled at s.
Two additional “field-scale” simulations (F1–F2) feature different wave types in a much larger domain ( km) with significantly lower resolution ( m, m, and s). Modeled ice thickness is fixed at 50 cm, and water depth is set to 50 m (F1) and 100 m (F2) as required to maintain a deep-water wave environment. These lower-resolution tests are less computationally expensive than the lab-scale simulations, but they do not provide detailed information about the flow field near the ice boundary.
In
Section 3.1, examples of modeled flow fields are presented for two of the lab-scale cases to illustrate the variations in the ice layer’s effect on fluid velocities. In
Section 3.2 and
Section 3.3, model-estimated velocity and Reynolds stress profiles for lab-scale cases L1–L6 are directly compared to theoretical predictions and laboratory measurements. The lower resolution cases F1–F2 are not sufficiently resolved to capture the details of the flow field near the ice, so they are excluded from these results. In
Section 3.4 and
Section 3.5, both the lab-scale and the field-scale model estimates of wave attenuation and wavenumber in ice are compared with theory and measurements from the lab and field. As noted earlier, this modeling study is part of a series of investigations of wave–ice interactions and the water–ice boundary layer. Some results from the theoretical analysis of [
26] and laboratory experiments by [
6] are reproduced below for comparison with model output.
3.1. Flow Field Effects
Sample velocity fields extracted from a near-surface, ice-covered location in model cases L1 and L6 are presented in
Figure 2. Vorticity is elevated in the vicinity of the ice layer in both panels, as the velocity damping generally has reduced flow in the direction of wave propagation. However, vorticity is considerably lower for case L1 (top panel), in which shorter waves (
T = 1 s,
H = 2 cm) propagate through a 1 cm ice layer with a relatively low drag coefficient,
. For case L6 (lower panel), the ice layer has a significantly larger drag coefficient,
, and imposes a visible negative correction on horizontal velocities near the surface, affecting both positive and negative components for the larger waves (
T = 1.34 s,
H = 6 cm). This correction is too small to be seen in the flow field for case L1. There is little or no damping of the vertical velocity components in either case.
It is worth noting that, because the surface ice boundary in these simulations is moving vertically with the waves, it produces a different effect on the fluid than would be experienced along a relatively fixed bottom boundary at the seabed (for example, under shoaling waves near the beach). This surface boundary moves higher under wave crests and lower under wave troughs. Over several wave periods, the net shear and drag effects of the ice on the water are also spread out vertically over the water column, to an extent determined by the amplitude of the propagating waves. This spreading effect plays a role in how we define the “thickness” of the ice layer and is addressed further in the following section.
3.2. RMS Velocity Profiles
Vertical profiles of RMS velocity from the adapted wave–ice model are computed for cases L1–L6 along a transect at a location corresponding to
in
Figure 1. Profiles are obtained after each simulation has reached a steady state (i.e., wave height at location
remains approximately constant). In each case, along-tank RMS velocities are computed as
where
is horizontal velocity at depth
z and time
t,
is the sample rate, and
is equivalent to the time for 10 wave periods
1. Near-surface sections of the profiles for the six lab-scale wave configurations are presented in
Figure 3 (solid lines), along with corresponding profiles that would be obtained in open water for each wave type (dashed lines). A cyan-shaded region is included near the top of each panel to indicate the approximate depth range that was primarily or wholly occupied by ice during each wave period (arbitrarily set to wave amplitude plus half of ice thickness in each case).
As expected, the magnitudes of the RMS velocities are significantly greater for the larger waves in cases L4–L6. The divergence of the modeled profiles (solid lines) from the open-water profile (dashed line) in the upper several centimeters of the water column provides evidence of a wave-generated boundary layer immediately below the ice interface. Due to the eddy generation and drag effects of the semi-permeable ice cover, this type of boundary layer is expected to be significantly larger in magnitude and vertical extent than that seen in open water when ice is not present [
26], and such an under-ice boundary layer can often be turbulent [
23]. The increased velocity variance is expected to be greatest near the water–ice interface and taper with increasing depth. In these model results, the boundary layer region extends to roughly three times greater depth for the L4–L6 waves (
H = 6 cm) than it does for the L1–L3 waves (
H = 2 cm), implying that the vertical depth to which the boundary layer affects RMS velocities is dependent on the wave amplitude. The magnitudes of the boundary layer velocities also increase with increasing values of the ice drag coefficient,
, although this effect appears somewhat more muted for the larger coefficient values (i.e., between L2–L3 or L5–L6).
As detailed field data are not presently available, results from Yu [
26] and Orzech et al. [
6] are instead used here to establish “ground truth” values for the magnitude and vertical extent of RMS velocity profiles and boundary layers that would be expected under these conditions. A comprehensive, detailed picture of the velocity field properties expected in this near-surface region is provided by the authors of [
26], who completed a theoretical wave–ice analysis in which the surface ice was represented as a continuous viscous layer. Direct measurements of these under-ice velocity fields were also obtained in lab tests conducted by the authors of [
6], who employed a particle imaging velocimetry (PIV) system to record subsurface water velocities as waves propagated through ice floes in a refrigerated laboratory wave flume.
Sample velocity profiles from both the theoretical analysis and the laboratory experiment are replotted in
Figure 4 below for comparison with the model results for cases L1–L6. As the theoretical profiles were originally calculated for a non-dimensional domain, both theoretical and laboratory results are presented in a non-dimensional form to allow them to be directly compared. Following the conventions of [
26], depths in each panel are normalized by the nominal ice thickness,
, and velocities are normalized by the linear open-water velocity at the water–ice interface. Depths are also adjusted so that
corresponds to the normalized depth of the interface. Both the theoretical and lab results of
Figure 4 appear qualitatively similar to the dimensional profiles presented for the model in
Figure 3. In all of the profiles (solid lines), the boundary layer is visible as a divergence of RMS velocity from the (dashed) open-water profile in the region near the water–ice interface. In the left panel of
Figure 4, the three different theoretical curves demonstrate how the vertical extent of each boundary layer is affected by both the relative viscosity of the water and the relative thickness of the ice layer.
To more directly compare modeled RMS velocity profiles with the theory and lab results, the model profiles of
Figure 3 are replotted in a normalized form in
Figure 5. As in
Figure 4, the plotted velocities for cases L1–L6 are now normalized by the theoretical open-water velocity at the interface. Depths are normalized by a representative ice region thickness,
which is set equal to the wave amplitude plus half the modeled ice thickness, as in
Figure 3. Normalized depths are also adjusted so that zero depth now corresponds to the bottom of the representative ice region, as with the theoretical and experimental profiles.
Examining the first two theoretical cases in
Figure 4 (with viscosity ratios
1 × 10
4 and 1 × 10
2, respectively), it is apparent that when effective water viscosity
increases relative to ice viscosity
, the vertical extent of the boundary layer also increases. Comparing the second and third theoretical cases (both with
1 × 10
2, but with the latter case featuring a 3× thicker ice layer), it can be seen that an increase in the ice layer thickness
(relative to wavelength) can produce a further increase in the boundary layer’s relative vertical extent.
The RMS velocity profile obtained from lab measurements by [
6] (
Figure 4, right panel) contrasts significantly with the theoretical profiles. This profile was measured during a simulation in which small, monochromatic waves (
T = 1 s,
H = 4 cm) propagated through a 3.8 cm ice layer. The measured boundary layer exhibits a considerably greater relative vertical extent than those predicted by viscous layer theory. Based only on the results described for the theoretical profiles, it is logical to conclude that the velocity variance of the water immediately beneath the ice was significantly enhanced by turbulence in the experiment, and that the experimental ice thickness was unusually large compared to the waves (both of which were likely true). However, it seems doubtful that these characteristics of the lab simulation would be sufficient in themselves to extend the measured velocity profile so much deeper than the theoretical profiles.
In normalized form, the results for the six model cases (
Figure 5) are surprisingly consistent in both magnitude and the vertical extent of the boundary layer region. The normalized velocity magnitudes in the boundary layer are still affected by the size of the ice drag coefficient,
, but they do not appear to be strongly affected by the specific wave period. Employing the wave amplitude in the normalization of the water depth appears to have largely eliminated the variation in the vertical extent of the boundary layer that was seen for the two different wave configurations in
Figure 3. The model profiles are qualitatively comparable to both the theoretical and the measured RMS profiles of
Figure 4. However, though the vertical extent of the theoretical boundary region is only a fraction of the ice thickness, for both the measured and modeled boundary layers, this normalized extent is roughly of order one, i.e., the boundary layer depth is roughly the same as the nominal ice thickness value.
To explain this difference, we must consider that, unlike the theoretical analysis, the water–ice interfaces in the lab and model environments are both vertically oscillating, so that the boundary layers themselves are moving up and down several centimeters as each wave passes. The region occupied by ice in the lab and the model is determined by the amplitude of the waves, which carry the ice higher in their crests and deeper in their troughs. Within this expanded ice region, locations near mean sea level have ice in them at all times, whereas other neighboring depths are only occupied by ice for part of a wave cycle. As described in
Section 3.1, this oscillation effectively spreads out the measured boundary layer region over a larger range, roughly proportional to the wave amplitude. When averaged over multiple wave periods, this vertical displacement of the existing boundary layer affects RMS velocities deeper in the water column in a manner that is not considered in the theoretical analysis.
In addition, both the modeled and lab ice include some degree of permeability, whereas the theoretical approach of Yu [
26] enforces an impermeable, fixed boundary at the mean water–ice interface. The flux of water into and out of the ice layer in the lab and the model augments the ice-generated turbulence, an additional effect that is not considered by a viscous layer representation. As wave height and/or ice thickness increase, this turbulent contribution grows in magnitude and is advected deeper into the water column below the ice.
3.3. Reynolds Stresses
The wave-induced Reynolds stress,
is computed from the time-averaged product of the wave-induced velocities,
, and represents the mean momentum flux that is generated by the wave fluctuations. For rotational wave motion,
, but when this motion becomes irrotational (as in the vicinity of surface ice), the Reynolds stress becomes non-zero. Several examples of Reynolds stress profiles computed from viscous layer theory by Yu [
26] and a single profile measured in the lab by Orzech et al. [
6] are reproduced in
Figure 6 below to again establish a baseline for what might be expected from the adapted wave–ice model. As in the preceding section, these profiles are displayed in a non-dimensional form, with the water depth
z normalized by the ice thickness
h, and the Reynolds stress computed with
u and
v velocities that have been normalized by the open-water velocity amplitude at the level of the water–ice interface. The lab-measured stress profile is computed for the same test case as was used in
Section 3.2, with monochromatic waves of period
T = 1 s and height
H = 4 cm and a surface ice thickness of 3.8 cm.
Within the surface ice layer
, the values of both theoretical and experimental Reynolds stresses are positive, which is consistent with the expectation that the phase of the
v oscillation is leading that of the
u oscillation by less than 90° in the viscous fluid [
26]. As we move downward and away from the wave–ice interface inside the wave boundary layer
, the stress first decreases rapidly in all profiles, then reaches a negative peak
, before finally returning to zero at greater depths. As seen earlier with the RMS velocity profiles, the experimental profile again appears vertically elongated relative to the theoretical profiles, likely due to both different conditions and differing treatments of the ice layer. The experimental profile also has a much larger negative peak and almost no positive Reynolds stress values.
Profiles of normalized Reynolds stresses are presented in
Figure 7 for datasets from cases L1–L6 of the adapted wave–ice model, using the same normalization methods as above. These model profiles, each computed over 10 wave periods, again exhibit behavior and magnitudes that are qualitatively similar to the theoretical results. In all but one of the six model profiles, stress values consistently decrease in the region immediately below the water–ice interface
, usually reaching a negative minimum value before returning to zero at greater depths. The vertical extent of the boundary layer and negative peak region is somewhat greater in the model results than in Reynolds stress profiles obtained from viscous layer theory, and in this respect, the modeled profiles are once again more like the experimental profile. In magnitude, however, the range of the modeled stress profiles is this time much closer to that of the theoretical profiles, with consistently positive values in the shallowest region close to the ice interface and only a slightly negative minimum stress value
2.
3.4. Wave Attenuation
Wave attenuation is estimated through an analysis of surface elevation time series at two selected ice-covered locations in the model domain. The model is allowed to reach a steady state at both locations (i.e., height remaining approximately constant) before the time series data are extracted. The wave amplitude attenuation rate,
, is calculated over a central portion of the ice region via the formulation
where
is the distance between the two measured locations (four times the open-water wavelength), and
is the RMS wave height at location
, with
j =
a or
b (see
Figure 1). In cases L1–L3 from
Table 1,
= 6.24 m, and in cases L4–L6,
= 11.2 m. In the cases F1 and F2, the values of
are 224.5 m and 400 m, respectively. The RMS wave heights are determined from the standard deviation,
of the surface time series at the two locations as
where
is free surface elevations,
is their mean, and
M is number of timesteps averaged.
Attenuation coefficients obtained for the eight ice and wave configurations are summarized in
Table 2. These results are consistent with expectations for changing values of either the ice drag coefficient or the wave properties. As
is increased, the wave attenuation rate
also increases for a given wave type. When ice drag is held constant, the shorter period waves are consistently more attenuated than the longer waves.
As with the velocity and stress profiles of the preceding sections, it is again helpful to express wave attenuation rates in terms of non-dimensional dynamic parameters, so that data from different scales may be directly compared. When the effective ice thickness,
, is used as a normalizing quantity for these data (again following the conventions of [
26]), the non-dimensional wave frequency,
, can be computed as
where
is the open-water wave frequency, and
m/s is taken as the gravitational acceleration in polar regions. Correspondingly, the non-dimensional wave attenuation rate,
, is obtained from
This representation allows for small-scale lab results to be examined together with dynamically similar field data on the same curve. Ultimately, this capability can facilitate the extraction of generalized attenuation information from lab and/or model simulations, which can then be used for the configuration of field-scale forecasts in preparation for future field expeditions.
Normalized attenuation results for the eight model configurations of
Table 2 are presented graphically (in log–log format) in
Figure 8, along with additional results from both lab and field experiments. The model results line up reasonably well with the linear best fit to the other lab and field data that was determined by [
26], with modeled attenuation rates slightly lower than the trend line in all eight cases.
3.5. Wavenumber in Ice
The wavenumber in ice,
, is also expected to vary with different wave and ice conditions. Experimental results from [
6] indicate that the value of
for a given open-water wavenumber generally increases with increasing ice thickness. These findings are supported by theoretical predictions of wave dispersion in ice based on a mass loading ice representation [
34]. However, an opposite trend (i.e., decreasing
with increasing ice thickness) is predicted when an elastic plate ice representation is employed [
35,
36].
Using the surface elevation data at the two ends of the test range shown in
Figure 1, the wavenumber in ice is estimated by following the procedure outlined in Parra et al. [
35]. For each trial, the time series from the two locations are compared, and the time required
for a selected wave crest to travel the length of the test range (
) is determined. The wave phase speed is then computed as
The wavenumber is then determined by making use of the relationship
in which
f = 1/
T is the wave frequency (1 Hz or 0.75 Hz, depending on the specific trial).
Values of
are summarized for the eight model test cases in
Table 3, along with their corresponding open-water values,
, where
is the open-water wavelength. For each case, the value of
is fixed at
, as indicated in
Figure 1. The values measured for
increase consistently with increasing ice drag
. In cases L1–L3 as well as L4–L6, this results in a phase speed that decreases slightly as
increases. This decreasing relative phase speed is reflected in the increasing values of
in
Table 3. For cases F1 and F2, both
and
are significantly greater for the 10 s waves, and the drag
remains constant, producing the expected smaller value of
for the larger F2 waves. The wavenumber magnitudes for the cases where drag is weak (L1, L4, F1, and F2) are all close to their open-water values, whereas those for ice with greater drag (L2, L3, L5, L6) are somewhat larger.
When normalized by ice thickness in the same manner as the attenuation coefficients, the wavenumber values for all eight cases are reasonably close to those predicted by theory for waves of their frequency passing through a viscous surface layer of comparable thickness, and they align well with several sets of lab measurements (
Figure 9). Note, however, that for this lower range of normalized frequency, the wavenumber curve based on viscous-layer theory is essentially indistinguishable from an open-water wavenumber curve (see, for example,
Figure 2a in Yu [
26]).