Macro-Turbulent Flow and Its Impacts on Sediment Transport Potential of a Subarctic River during Ice-Covered and Open-Channel Conditions

Macro-turbulent flows (i.e., coherent flow structures reaching through the whole water column), have not been studied widely in northern seasonally frozen rivers during both open-channel and ice-covered flow conditions. Thus, we aim: (1) to detect and compare the macro-turbulent flow, both at open-channel and ice-covered flow conditions; (2) to explore spatial variation of macro-turbulent flow characteristics within a meander bend; and (3) to detect the effects of near-bed layer velocity fluctuation on bedload transport during differing overall flow conditions. The analyses are based on 5–10 min-long acoustic Doppler current profiler (ADCP) measurements from a subarctic river. The ice-covered low flow, and open-channel higher and lower flow conditions were measured over the period of 2016 to 2020. This study found that macro-turbulent flow existed at all measurement locations under both open-channel and ice-covered flow conditions. Macro-turbulent flow was most consistent and obvious in the streamwise velocity component, and in particular at the inlet and outlet of the investigated meander bend. During all seasons, the near-bed velocities consistently exceeded the sufficient amount for sediment transport. At inlet and outlet areas, the greatest near-bed velocity fluctuation across the critical threshold for sediment transport coincided with the measurement times having frequent macro-turbulent flow.


Introduction
River ice can have a significant impact on the distribution and magnitude of flows in the river channels [1][2][3]. Under ice-covered conditions, increases are seen in the spatial variability and characteristics of depositional and erosional locations compared to the open-channel conditions [3]. However, these observations are based on studies conducted during only one season and on measurements that are too short (less than 1 min) to capture larger pulsating structures of high and low flows. The impact of this flow pulsation on the sediment transport potential and morphology of rivers during open-channel conditions has been studied over several decades (e.g., [4]), and it has been identified that this characteristic likely plays a major role in bedload sediment transport [5,6]. This has The underwater roughness of the ice-water interface was defined by lifting up sections of ice from drill holes and visually inspecting the underside. In addition, submerged camera photographs and videos were taken under water in 2017 and 2018 to see the ice-cover ( Figure 2). These photographs and footage also showed that the riverbed material was moving in ice-covered  Table 1. The measurement times, applied sensors, discharge, and notes. In the acronyms of the measurement times W = winter (low flow period, February), S = spring (snow-melt flood period, May), A = autumn (low flow period, September). These raw ADCP data sets (Matlab format) are shared with the paper (DOI: 10 Table 3. The basic flow characteristics of each of the analyzed ice-covered measurements from their measurement location (see Figures 1 and 2) within cross-sections (CS). The average, minimum, and maximum values of depth-averaged (D-A) velocities were defined. In some winters the CS1, CSA, and CS2 had mid-channel bars, where the ice was bottom-fast (i.e., reached the riverbed). During those winters, the channel width was measured from both sides (left and right bank) of the mid-channel bars. NMB = no mid-channel bar.  Section 3 and Table  1), and their locations were within the cluster of other CS1 measurements. The GoPro camera photograph shows the riverbed sediment and the low undulation of the ice-water interface at the CS1 on 10 February 2018.

Measurements of Flow Characteristics
According to Demers et al. [9], 10 min-long ADCP measurements with 1 s sampling frequency, (i.e., which have 600 samples), were sufficient to capture these macro-turbulent structures during icecovered flow conditions. It was found that a recording time of 160 s (i.e., samples) is enough for detecting the macro-turbulent coherent flow structures [11]. Based on the preliminary assessment of the local flow conditions at the Pulmanki River, the flow characteristics were measured with an ADCP (Sontek) for around a 5-10 min-long period at each measurement location. The sampling frequency was 1 s. The measurements were conducted in winter (February), spring (May), and autumn (September) during the years 2016 to 2020 (Table 1). Either the M9 or S5 ADCP sensor was used, depending on the measurement time and availability of the sensors ( Table 1). The S5 sensor has four 3.0 MHz beams, but the M9 sensor has both four 3.0 MHz beams but also four 1.5 MHz beams. Both sensors have accuracies up to +/-0.25% of measured velocity (+/− 0.2 cm/s) [16].
During ice-covered measurements, the sensor was lowered through drill-holes, with a pole attached to it. The pole was mounted to a tripod, and the Global Positioning System (GPS) antenna was attached on top of the pole or tripod, or as high as the wiring allowed. This enabled the sensor and its GPS to remain steady during the measurement period. The sensor's transducer was lowered at the level of the ice-water interface (i.e., the bottom of the ice cover), so that it would measure the turbulence of the entire water column, including the influence of the ice cover. The minimum water depth, where the sensors are able to measure, is circa 0.2 m [16], due to the riverbed (side lobe interference) or the water surface (blanking distance) interference (see [16,17]). Due to the M9 sensor's blanking distance (i.e., the distance between the transducer and the first measurement cell within vertical water column), the first cells were measured 0.06 m below the transducer. When the S5 sensor was used (note: only in open-channel conditions, Table 1), the first measurement cell was measured at 0.17 m depth below the transducer. As the minimum total depth value of all measurement times and locations (Tables 2 and 3) was 0.34 m (on 6 February 2020, at CS1), there was then more than 20  Table 1), and their locations were within the cluster of other CS1 measurements. The GoPro camera photograph shows the riverbed sediment and the low undulation of the ice-water interface at the CS1 on 10 February 2018.
The study is based on data from one meander bend, where the sub-water area is approximately 20 m wide during low flow open-channel conditions [13,14]. The slope of the channel has been defined at the study bend to be c. 0.00034 m/m during the autumn low flow period (September 2012) and 0.0003 m/m during the spring discharge peak (May 2011) [13]. During spring flow conditions in May 2011, the stream power values of 0.31-0.7 N/m 2 occurred at the high velocity core locations along the meander bend. Lotsari et al. [3] also defined the bedload's D50 grain size as 0.5 mm based on the measurements done in ice-covered discharge conditions of 1.1 m 3 /s.
The locations for macro-turbulent flow measurements were selected, based on initial cross-sectional and temporally short measurements, to be the high velocity core locations of each of the seven cross-sections (Figures 1 and 2). This location occurred in the middle of the channel always in the first cross-section (CS1) at the inlet area, and during some measurement times at the second cross-section (CSA), but at the outer bank (right bank) of the other five cross-sections. Due to the mid-channel bars at the inlet area (CS1 and CSA) during many years, the water depth was shallower there than in many other cross-sections ( Figure 2, Tables 1 and 2). The CS1 measurement location of the Pulmanki River represents flow coming to the meander bend from the straight river section.
The underwater roughness of the ice-water interface was defined by lifting up sections of ice from drill holes and visually inspecting the underside. In addition, submerged camera photographs and videos were taken under water in 2017 and 2018 to see the ice-cover ( Figure 2). These photographs and footage also showed that the riverbed material was moving in ice-covered conditions. The ice cover roughness was visually interpreted "as smooth-rough" (Figure 2). This visual classification followed the approach presented by Demers et al. [1]. The structure of the ice was also defined visually. The bottom layer of ice was translucent indicating black ice, similarly as Kämäri et al. [15] had observed in 2014 at the same study area.

Measurements of Flow Characteristics
According to Demers et al. [9], 10 min-long ADCP measurements with 1 s sampling frequency, (i.e., which have 600 samples), were sufficient to capture these macro-turbulent structures during ice-covered flow conditions. It was found that a recording time of 160 s (i.e., samples) is enough for detecting the macro-turbulent coherent flow structures [11]. Based on the preliminary assessment of the local flow conditions at the Pulmanki River, the flow characteristics were measured with an ADCP (Sontek) for around a 5-10 min-long period at each measurement location. The sampling frequency was 1 s. The measurements were conducted in winter (February), spring (May), and autumn (September) during the years 2016 to 2020 (Table 1). Either the M9 or S5 ADCP sensor was used, depending on the measurement time and availability of the sensors ( Table 1). The S5 sensor has four 3.0 MHz beams, but the M9 sensor has both four 3.0 MHz beams but also four 1.5 MHz beams. Both sensors have accuracies up to ±0.25% of measured velocity (±0.2 cm/s) [16].
During ice-covered measurements, the sensor was lowered through drill-holes, with a pole attached to it. The pole was mounted to a tripod, and the Global Positioning System (GPS) antenna was attached on top of the pole or tripod, or as high as the wiring allowed. This enabled the sensor and its GPS to remain steady during the measurement period. The sensor's transducer was lowered at the level of the ice-water interface (i.e., the bottom of the ice cover), so that it would measure the turbulence of the entire water column, including the influence of the ice cover. The minimum water depth, where the sensors are able to measure, is circa 0.2 m [16], due to the riverbed (side lobe interference) or the water surface (blanking distance) interference (see [16,17]). Due to the M9 sensor's blanking distance (i.e., the distance between the transducer and the first measurement cell within vertical water column), the first cells were measured 0.06 m below the transducer. When the S5 sensor was used (note: only in open-channel conditions, Table 1), the first measurement cell was measured at 0.17 m depth below the transducer. As the minimum total depth value of all measurement times and locations (Tables 2 and 3) was 0.34 m (on 6 February 2020, at CS1), there was then more than 20 cm of the water column to be measured. In all other measurement times and locations, the measured flow depth was deeper. Screening distance was defined as 0 m.
The ice thickness was also measured at each measurement location. Drill-holes were also made near to the banks of the cross-section of each ADCP measurement location to identify the edges of the flowing water area (i.e., where the surface ice reached the riverbed and confined the channel width). Measurements were used in Equation (1) [18,19] to determine whether or not the ice-covered flow was pressurized during the ADCP measurement times (Table 3).
where t i is the ice-cover thickness and W is the width of ice-cover (or channel) ( Table 3). According to Turcotte et al. [19] and their theoretical consideration, the pressurized conditions occur only when the equation holds true (i.e., t i is larger). Turcotte et al. [19] stated that in sub-arctic conditions (ice thickness up to 1.5 m) the confined channel conditions can occur in intermediate channels up to 60 m wide (i.e., W value). This equation is considered appropriate as the maximum river channel width of the study site measured during winter was approximately 24 m (Table 3), within the range where the equation has been previously applied [19]. During open-channel conditions, the ADCP measurements were conducted with the sensor installed on a floating platform ("HydroBoard" by Sontek). Due to the platform of the M9 sensor, the transducer depth was then 0.06 m below the water surface. Screening distance was 0 m. For the September 2017 measurements, a different M9 sensor platform was applied which floated slightly deeper. The platform for the S5 sensor floated deeper, and the transducer depth was 0.11 and screening distance was 0 m. During each open-channel flow measurement, the sensor was kept still either Water 2020, 12, 1874 8 of 30 manually or by attaching the platform to shore with poles and ropes. After initial comparison between the M9 and S5 data sets, ambiguities and malfunction of the S5 sensor were observed in most of the resulting measurements. CS1 was the only measurement location where data was good each measurement time. Thus, the S5 measurements of CS1 were included in the final analyses from 25 May 2016, 23 May 2018, and 8 September 2018 ( Table 1). Note that their detailed locations are not marked in Figure 1 as they were within the cluster of CS1 measurements already presented in that figure.

ADCP Data Analyses
For analyzing the differences between the measurement times, the ranges of depth-averaged velocities of each measurement (received directly from the measurement files) were reported, in addition to the flow depth and the ice-thickness above the water column (Tables 2 and 3). Each measurement file consisted of "samples" (i.e., one single vertical profile) measured at 1 s resolution. Each sample consists of depth "cells". As the sensor adjusts the cells based on the total depth of the water column, the vertical "thickness" of each cell varies between 2 and 6 cm in each sample. Thus, the number of cells can differ between measurement locations depending on the available depth. Based on all analyzed raw measurement files, the average, median, and mode values for the cell amount in vertical direction were 24, 20, and 18, respectively. Theoretically, each sample should have the same number of cells in vertical direction, when the sensor is kept static, as was done in our measurements. In practice, however, the transport of sand and ripples on the riverbed can cause variations of the riverbed depth by a couple of centimeters, while flow magnitude variations can similarly lift or lower the sensor by a couple of centimeters during the measurement. Therefore, in some measurement files there were one or two cells difference in the start depth and end depth of the cells between consecutive samples. This was taken into account in our further analyses. Each of the measurements contained approximately 300 to 600 samples to be analyzed (Tables 2 and 3). Only samples that showed a signal-to-noise ratio (SNR) of at least 5 dB were analyzed [9].
The data sets were further processed in R v. 0.3.6 [20], using the package 'adcp' v. 0.1.0. The Supplementary Material (see R package folders S1 and S2) contains the R code used to process the data. In addition, the package can be found online (https://github.com/coffeemuggler/adcp). The package allows importing the ADCP data exported as Matlab files (selecting the ENU components for export) and arranging the data in an appropriate object structure. Upon data import, the package calculates the overall speed value (i.e., the vector sum of all three spatial velocity components), by using equation √ E 2 + N 2 + U 2 as well as the rotated horizontal components, as along-stream (i.e., streamwise) and across stream data, using the geographic orientation data contained in the Matlab file. The package also allows linear interpolation of the measurement data to equal depth intervals throughout all samples, a prerequisite for subsequent analysis steps that depend on regularly spaced data.
Using the interpolated data, mean depth-wise velocities (for all components) were calculated (i.e., the arithmetic means of cells from all samples at a common depth, thus "depth-wise" mean). This interpolation converted the measurement data into a raster form, with time in horizontal dimension and depth in vertical direction. The number of cells was set to 20 for each measurement file, so that the original vertical variation of the flow was retained. The cell amount of the raster files was defined so that it best corresponded to all winter, autumn, and spring measurement data files (the average, median, and mode values for the cell amount in vertical direction were 24, 20, and 18, respectively). The deviation of velocity values was calculated for each cell from the average depth-wise velocity and used as a first-order proxy for cells with higher and lower than average values. To identify clusters of consistently high and low flow, we applied LISA (Local Indicators of Spatial Association [21]), using the R package ncf v. 1.2-9 [22]. LISA calculates the spatial autocorrelation in a focal window using Moran's I method and returns, among other output, the correlation value and two-sided p-value for each cell. We applied LISA with a 3 by 3 cell window [9] to the along-stream, across-stream, vertical and total velocities. In contrast to Demers et al. [9], we used the p-value (0.05 threshold) as criterion for the validity of the output. Finally, we classified the resulting data sets based on the following criteria (see [9]):

1.
Statistically significant high flow cluster: significant (p < 0.05) autocorrelation and flow velocity greater than the depth-averaged mean; 2.
Statistically significant low flow cluster: significant (p < 0.05) autocorrelation and flow velocity less than the depth-averaged mean; 3.
The occurrence of statistically significant macro-turbulent high and low flow wedges, which reached from surface layers to the near-bed layers, were analyzed from the data sets. The wedges were considered "strong" if they reached throughout the whole water column, and "weak" if there was a 5-20% gap between cells within the water column. This depended also on the temporal length of the cluster (i.e., the temporally longer the cluster, the greater the vertical gap allowed). If the cluster was slightly tilted, this was still considered as a macro-turbulent cluster. Thus, despite possible tilting, if the cells were in contact with each other and reached from the surface to the riverbed, this was considered a strong wedge. If there was a significant low flow cluster next to a high flow area, and vice versa, the macro-turbulent clusters were counted as separate clusters.
Finally, it was determined if near-bed streamwise velocities reached the critical velocity threshold needed for the incipient motion of the sediment particles (see Section 3.3). Note that near-bed velocities mean the values of the deepest recorded cells. These were located approximately 6-37 cm above the channel bed, depending on the ADCP measurement location and time. According to the Sontek M9 sensor's manual [16], potential data contamination in the cell that would be partially or fully touching the riverbed, or potential side-lobe interference at the end of the profile, leaves the section of water at the bottom of the water column unmeasured. As a result of the topography of the measurement location and the cell sizes adjusted by the sensor itself, the data of different measurement locations can have different gaps between the riverbed and the deepest recorded cell.

Sediment Transport Measurements
The bed-load samples were analyzed to determine if there is overall sediment transport under ice cover and to define the transported grain sizes. The bedload transport was measured on 6 February 2020 using a Helley-Smith (HS) [23] sampler with an intake aperture of 152 mm. The measurements were 6 min long. On 6 February 2020, three samples were collected from CS1 and CS4 ADCP measurement location, and two samples from CS2 and CS3 measurement locations ( Table 4). The samples were dried in 105 • C degrees. After drying, the dry mass of the sediment catch was used for calculating the bedload transport (mg/(m·s)). The largest sample of CS1 (HS3) and CS4 (HS1) were dry sieved, and their D10, D50, and D90 grain sizes were defined ( Table 4). The samples of CS2 and CS3 were small, and their sieving would have included uncertainties, as the scale with an accuracy of 0.1 mg would not have recorded small amounts of material on sieves.
The critical velocities needed for the material to be moved were analyzed based on the Hjulström [24] diagram and its velocity threshold for erosion. Those were defined for the D10, D50, and D90 sediment grain sizes. The critical velocity needed to move the largest D50 (0.318 m/s) and D90 (0.432 m/s) grain sizes of the bedload sample "CS4 HS1" (Table 4) were compared against the measured near-bed layer streamwise velocities and their fluctuation ranges. Note that based on [24], the critical velocity for the D90 sized particle (of "CS4 HS1" sample) to deposit is required to be only 0.07 m/s if the particle is already in motion.
Despite these possible sieving errors due to a small sample size, the corresponding author [3] had previously sieved one HS sample, which had been measured on 15 February 2017, and which had the dry weight (1.25 g) close to the ones of CS2 and CS3 measured on 6 February 2020. Lotsari et al. [3] measured the 2017 HS sample from the straight river reached upstream of the CS1 ADCP measurement location. The distance between CS1 and the measurement location of 15 February 2017 HS sample was 320 m. Bedload transport was 22.84 mg/(s·m) on 15 February 2017 [3]. The D50 grain size of that bedload sample was 0.5 mm, based on dry sieving [3]. The D10 and D90 grain sizes were 0.2192 and 0.9024 mm, respectively (Table 4). These D10 and D90 have not been published in Lotsari et al. [3]. Based on [24], the critical velocity for this smallest analyzed D10 value of 2017 sample to deposit is required to be only 0.016 m/s, if the particle is already in motion.
To test the amount of suspended load, we measured 500 mL sized water samples (Table 5)  In addition, the color (mg/L Pt) and turbidity (FTU: Formazine Turbidity Unit) were analyzed with YSI 9500 photometer (a Xylem brand) from the rest of the sample. Note, one sample from 21 May 2019 was 145 mL in total, and all of it was used for TSS analyses. Thus, no color or turbidity values were analyzed from that sample.
Lotsari et al. [3] also measured and analyzed TSS transport in February 2017 (same date as the HS bedload sample) via similar methods used in this study. The three samples of Lotsari et al. [3] showed TSS transport of 0.22, 0.66, and 0.67 mg/L. In addition, their fourth sample showed 0 mg/L TSS transport, which could have been an erroneous measurement (note that this fourth sample had not been included in previous work by Lotsari et al. [3]).

Spatial and Temporal Variation of the Depth-Averaged Velocities
The stark contrast between average and maximum depth-averaged velocities of each of the measurements indicates that there was pulsating flow (Tables 2 and 3). The average values of the depth-averaged velocities of each profile show that except for W2019 (see the acronyms from Table 1), the winter season velocity was highest at the inlet area (measurement location CS1): 0.57 m/s (W2016), 0.378 m/s (W2017), 0.514 (W2018), 0.369 m/s (W2020). In W2019, the highest average value of these depth-averaged velocities was at the CSA (0.25 m/s), which is close to the upstream end of the meander bend. The lowest average value of these depth-averaged velocities during winter measurements varied between years: they were at CSC in W2016 and W2018, at CS4 in W2017, at CS1 in W2019, and at CS2 in W2020. Overall, out of winter measurements, the highest average depth-averaged velocities were on W2016, and lowest were on W2017. The greatest maximum depth-averaged winter velocities located at CSA (1.787 m/s) in W2016, at CSC (0.555 m/s) in 2017, at CS3 (0.942 m/s) in W2018, at CSA (2.166 m/s) in W2019, and at CSA (2.124 m/s) in W2020. Therefore, the measurement location of the maximum depth-averaged velocity was mostly at CSA ( Table 3). The largest range between the average and maximum depth-averaged velocities was in W2019. Then the maximum depth-averaged velocity was the highest of all measurements, even greater than during open-channel conditions.
The open-channel conditions showed similarities to winter measurements regarding the highest and lowest average depth-averaged velocity locations (Table 2). In spring 2017, the highest flow (0.647 m/s) was in CS1, and slowest (0.101 m/s) was in CSC. In spring 2018, the highest was in CSB (0.422 m/s) and lowest in CS4 (0.064 m/s). Again, in autumn 2017 the highest flow was in CSB (0.667 m/s), but lowest was in CS3 (0.122 m/s). Regarding the maximum depth-averaged velocities, the highest were at CSB (0.827 m/s), CS1 (1.316 m/s), and CSB (0.587 m/s) on S2017, A2017, and S2019, respectively. In W2019 ice-covered conditions, the maximum depth-averaged velocities were approximately four times higher than during S2019 open-channel conditions (Tables 2 and 3), even though the discharge was ten times greater in S2019 (Table 1) than the preceding winter measurements. These indicate the pressurized conditions in ice-covered conditions (Table 3), as also calculated based on Equation (1).

Marco-Turbulence Patterns
Based on the preliminary data check and the analyses of depth-averaged velocity, the CS1 (inlet area), CSB (apex), and CS4 (outlet area) were selected for the most detailed analyses. These locations showed the most prominent differences in flow characteristics between measurement times. In addition, when all measurement locations were analyzed, macro-turbulence was formed predominantly (W2018, Table 1) in ice-covered conditions and when the open-channel conditions had the highest discharge (S2017, Table 1), i.e., when it differed the most from the ones analyzed in [9]. Note that the discharge of S2016 and S2018 would have been even higher, but there was no good quality data available from all measurement locations during those time steps. The figures of raw streamwise (i.e., along-stream) and vertical velocities of the selected measurement times and locations are included as the Supplementary Material (see Figures S3-S10). In addition, the macroturbulent flow analyses results (LISA) based on vertical flow velocity have been included as Supplementary Material (see Figures S11 and S12). All analyzed data sets reached a good signal-to-noise ratio (SNR) of at least 5 dB [9].

Spatial Variation of Macro-Turbulence under Ice-Covered Conditions
The W2018 results are presented as an example of ice-covered flow condition from all measurement locations. Overall, during ice-covered conditions, the high velocity core was located in the middle of the water column, and not within the surface layer ( Figures S3-S8). The streamwise (along-stream) velocity and vertical velocity data of W2018 revealed pulses of higher flow (Figures S3 and S4). These pulses reached through the whole water column, from surface to the near-bed layers. The inlet area had high streamwise velocities of 0.  Figure S4).
The statistically significant high and low flow velocity clusters were clearly apparent from the measurements of W2018 (Figure 3, Tables 5 and 6). Even though the streamwise strong (i.e., unbroken) low flow wedges (n = 5), reaching from surface layers to near-bed layers, outnumbered the strong high flow wedges (n = 4), the weaker wedges (which were broken by couple of cells in vertical direction) of high flow greatly outnumbered (n = 9) the weak low flow wedges (n = 3) at CS1 location. Both strong and weak macro-turbulent flow wedges were the most equally spaced at this inlet area during W2018 measurement time. The statistically significant flow wedges, which reached from the surface to riverbed, occurred in CSA also, but those were temporally shorter, being only 1-2 s long, than at the CS1, where 4-5 s long wedges occurred. Least strong wedges occurred at CSC and CSB at that time, and the statistically significant high and low flow clusters were only near the riverbed and ice-water interface layer. Thus, the statistically significant clusters did not reach clearly through the water column. Within the apex area (from CS2 to CSC), the detection of those wedges was overall more difficult as the statistically significant high and low flow areas varied there more randomly. At those measurement locations, there was often either high or low flow cluster at the surface layer, and then the opposite clusters (i.e., low or high flow cluster) at the near-bed layer. However, at the CS4 the macro-turbulence was again more visible, and temporally wider clusters of statistically significant high and low flow wedges appeared. The strong and weak high flow wedges outnumbered the low flow wedges at CS4. Thus, based on the streamwise flow, the macro-turbulence was the clearest at the inlet area, second clearest at the outlet area, and less clear in the apex measurement locations during the mid-winter ice-covered, mostly pressurized (Table 3), low flow conditions.  The statistically significant high and low flow wedges were also seen from the vertical velocity data of W2018 ( Figure S11, Table 6). Similar to the streamwise velocity, it was the clearest in the inlet area (CS1) during that measurement time, and the least visible at the apex area (from CS2 to CSB).   Tables 2 and 3). Note that the average values are calculated based on the measurement times of around equal length (i.e., others except S2018, W2019, and W2020). At CS4, there was spring measurement only from S2017, thus the average value was not calculated. The statistically significant high and low flow wedges were also seen from the vertical velocity data of W2018 ( Figure S11, Table 6). Similar to the streamwise velocity, it was the clearest in the inlet area (CS1) during that measurement time, and the least visible at the apex area (from CS2 to CSB). Different to the streamwise velocity is that the vertical significant flow clusters started to appear again right after the apex at CS3 and CSC. However, in measurement locations CSA-CS4 the statistically significant vertical high and low flow vertical wedges were shorter in time than at CS1 ( Figure S11).

Spatial Variation under Open-Channel Conditions
The S2017 measurements were used as an example of the open-channel flow conditions at each measurement location. The highest velocities of S2017 of each measurement location were closer to the water surface than the middle layers of the water column, which was the case in ice-covered flow conditions of W2018 (see Figures S3 and S5). When compared to the ice-covered flow conditions of W2018, the macro-turbulence in S2017 was more frequent and obvious in the open-channel flow conditions. The streamwise velocity values of CS1 showed that there was multiple high flow wedges of 0.75-0.8 m/s, which were separated with lower flow (~0.5-0.6 m/s) periods (see Figure S5). Note that when the streamwise velocity was high, often then the downward vertical velocities were high (see Figures S5 and S6). The vertical flows formed two layers on S2017 at the CS1, CSA, and CS2: there was clearly downward motion within the top section of the water column, and close to zero flow conditions or slightly upward flow within the bottom section of the water column. In winter (W2018) this phenomenon had been seen only at the CS1. Overall, the vertical velocities were slow in S2017 and the main direction of the flow was streamwise (Figures S5 and S6). In total, there were 25 statistically significant strong flow wedges at the inlet area (CS1) during 10 min-long measurements, out of which 13 were the high flow wedges (Figure 4, Tables 5 and 6).
Thus, on average there were 2-3 strong flow wedges per minute. This is a much higher frequency than encountered during any of the winter period measurements (Tables 5 and 6). In addition, at all other measurement locations, the occurrence of both strong and weak macro-turbulent wedges was more frequent on S2017 than W2018. Only at the CSA location was there one more strong high flow wedge on W2018 than on S2017, but the weak high flow wedges also outnumbered the winter conditions at that measurement location. However, when moving from the inlet to the apex area (CS2-CS3), similar changes in the characteristics of the statistically significant high and low flow wedges appeared in S2017 as in W2018 (Figures 3 and 4). The wedges became temporally narrower, (i.e., only a couple of seconds long in the apex area). In addition, the wedges became more broken. There was either high or low flow cluster at the surface layer, and then the opposite clusters (i.e., low or high flow cluster) at the near-bed layer in many occasions. This can be seen from the occurrence of the strong streamwise high and low flow wedges, whose numbers dropped, at the same time when the weak wedge numbers increased (Table 7). However, at the CSC and CS4, the streamwise flow wedges were again more obvious throughout the water column. The numbers of strong high flow wedges increased again closer to the frequencies of the inlet area's CS1 measurement location (CSC had 11 and CS4 had 8 strong high flow wedges; see Table 7). Thus, based on this, the macro-turbulence was the clearest at the inlet area, second clearest at the outlet area, and least clear in the apex measurement locations during the open-channel conditions of spring 2017.  Similar to the macro-turbulence observed based on the streamwise velocities, it was also more frequent on S2017 than on W2018 measurement time based on vertical velocities ( Figure S12, Table 7). The statistically significant high and low flow clusters did not vary much between the different parts of the bend. Therefore, the vertical zero, upward and downward flow formed statistically significant temporal clusters more consistently throughout the bend in open-channel conditions of S2017, than in ice-covered conditions of W2018.  . The LISA analysis results at each measurement location, calculated based on along-stream (i.e., streamwise) velocity of S2017. Measurement ID means the sample numbers (n), and therefore also represents seconds. The black dots indicate the occurrence times of statistically significant strong high and low flow wedges reaching from the surface to the riverbed. The grey dots indicate the times when there is weaker wedge (i.e., a cluster has broken with a few cells within the vertical water column).

Macro-Turbulence at the Meander Bend Inlet Area
The CS1, CSB, and CS4 were analyzed in more detail from all the measurement times, but only based on the streamwise velocities, as the W2018 and S2017 showed that the flow velocities were greatest in streamwise direction and the vertical flow was predominantly around 0 m/s. At the inlet (CS1), the statistically significant high and low flow wedges were less frequent during all other mid-winter measurement times than on W2018 (Figure 5). At the inlet area, least occurrences of macro-turbulence (from all measured ice-covered conditions) was observed on W2019, with only one weak high flow wedge and two weak low flow wedges. Thus, throughout the measurement, there was always at least a one or two cell gap in the wedges so that they did not reach throughout the water column. When detecting the raw velocity data sets ( Figure S7), the velocities were similarly high during other winters as on W2018. For example, on W2016, approximately 0.7-0.75 m/s streamwise velocities, which reached from the surface to the near-bed layers, occurred and lower velocities of 0.5-0.6 m/s appeared between these times of high flow wedges (see Figure S7).

Macro-Turbulence at the Meander Bend Inlet Area
The CS1, CSB, and CS4 were analyzed in more detail from all the measurement times, but only based on the streamwise velocities, as the W2018 and S2017 showed that the flow velocities were greatest in streamwise direction and the vertical flow was predominantly around 0 m/s. At the inlet (CS1), the statistically significant high and low flow wedges were less frequent during all other midwinter measurement times than on W2018 (Figure 5). At the inlet area, least occurrences of macroturbulence (from all measured ice-covered conditions) was observed on W2019, with only one weak high flow wedge and two weak low flow wedges. Thus, throughout the measurement, there was always at least a one or two cell gap in the wedges so that they did not reach throughout the water column. When detecting the raw velocity data sets ( Figure S7), the velocities were similarly high during other winters as on W2018. For example, on W2016, approximately 0.7-0.75 m/s streamwise velocities, which reached from the surface to the near-bed layers, occurred and lower velocities of 0.5-0.6 m/s appeared between these times of high flow wedges (see Figure S7).  Figure 3. Measurement ID means the sample numbers (n), and therefore also represents seconds. W2019 and W2020 measurements were shorter in time. The black dots indicate the occurrence times of statistically significant strong high and low flow wedges reaching from the surface to the riverbed. The grey dots indicate the times when there is weaker wedge (i.e., a cluster has broken with a few cells within the vertical water column). When open-channel conditions of CS1 were compared, the macro-turbulence was apparent during low flow conditions of A2016 (6.34 m 3 /s) and A2018 (3.69 m 3 /s), and during high discharge conditions of S2016 (Figure 6 and Figure S8). Note that in these A2016 and A2018 measurement time steps, both high and low flow wedges lasted longer in time, than in S2017, which had higher discharge condition at~11 m 3 /s. Thus, the overall macro-turbulent motion of the fluid took longer to pass by the measurement point, possibly due to the overall lower discharge and the flow depths of the autumn measurement times. In S2016, multiple turbulent wedges occurred consecutively, and resulted in temporally longer high flow periods. Despite the turbulent wedges taking longer in the shallower depth conditions, the highest velocities of these flow wedges were around 0.6-0.7 m/s on S2016, A2016, and A2018 ( Figure S8). Thus, these magnitudes were similar to the ice-covered conditions of W2018 and open-channel conditions of S2017, previously discussed. Thus, when the strongest statistically significant high flow wedges occurred at the inlet area of the meander bend, the streamwise flow velocity magnitudes were similar, despite the season in question. This verifies that the statistically significant flow variation throughout the vertical water column occurred at the inlet area during every season, whether being high spring or low autumn and winter flow, or whether being ice-covered pressurized (W2016, W2017, W2018: see Table 3) or non-pressurized (W2020, see Table 3) conditions.

Macro-Turbulence at the Meander Bend Apex
As the W2018 and S2017 data sets are already shown, the CSB measurement location representing the apex area had most often chaotic and unclear occurrence of the statistically significant high and low flow wedges. In addition to W2018 and S2017, there was also on A2016 and

Macro-Turbulence at the Meander Bend Apex
As the W2018 and S2017 data sets are already shown, the CSB measurement location representing the apex area had most often chaotic and unclear occurrence of the statistically significant high and low flow wedges. In addition to W2018 and S2017, there was also on A2016 and W2017 frequent statistically significant clusters of high and low flow either close to the riverbed or the water surface ( Figure 7). Thus, these clusters of CSB did not reach through the water column.
Water 2020, 12, 1874 23 of 33 Figure 7. The results of the LISA analysis at CSB during all measurement times based on streamwise velocity. The S2017 and W2018 results can be seen from Figures 3 and 4. Measurement ID means the sample numbers (n), and therefore also represents seconds. Note that the S2019 and W2020 measurements were shorter in time. The black dots indicate the occurrence times of statistically significant strong high and low flow wedges reaching from the surface to the riverbed, and the grey dots indicate the times when there is weaker wedge (i.e., cluster has broken with a few cells within the vertical water column).

Figure 7.
The results of the LISA analysis at CSB during all measurement times based on streamwise velocity. The S2017 and W2018 results can be seen from Figures 3 and 4. Measurement ID means the sample numbers (n), and therefore also represents seconds. Note that the S2019 and W2020 measurements were shorter in time. The black dots indicate the occurrence times of statistically significant strong high and low flow wedges reaching from the surface to the riverbed, and the grey dots indicate the times when there is weaker wedge (i.e., cluster has broken with a few cells within the vertical water column).
However, exceptions were W2016, W2019, S2019, and W2020, when strong wedges of either high or low streamwise velocities occurred throughout the water column (Figure 7 and Figure S9). Especially interesting is that on W2016 and W2019 the macro-turbulence was not frequent at the inlet areas (CS1 measurement location) but was more detectable at the apex area (CSB).
The streamwise velocities were high on W2016, A2016, and A2017, reaching 0.8 m/s during the high flow pulses ( Figure S9). During W2019, S2019, and W2020 measurement times, when the flow velocities were overall slow ( Figure S9), the statistically significant high and low flow clusters reaching from the surface to the riverbed were longer in temporal scale (Figure 7), than during the measurement times having overall higher flow magnitudes.

Macro-Turbulence at the Meander Bend Outlet Area
The statistically significant high and low flow wedges were apparent at the outlet area on A2016, W2017, and A2017 measurement times. During W2017 the strong high flow clusters occurred 10 times within the 10 min-long measurement of CS4 (Table 6) and were thus more frequent than in the inlet or apex areas. The velocities of the high flow wedges at CS4 during A2017 measurement time were around 0.5 m/s ( Figure S10). On A2016, the high flow pulses had 0.6-0.7 m/s streamwise velocities, which were similar to the velocities of the inlet areas, when the clearest macro-turbulent flow wedges occurred. Thus, CS4 also experienced these statistically significant high and low flow wedges during both open-channel and ice-covered seasons.
The flow structures of W2019 and W2020 differed from other measurements done at CS4, and also from other measurement locations (Figures 5-8). The streamwise flow was more chaotic on W2019, even though the statistically significant strong high flow wedges were more frequent than at the inlet area ( Table 6). The streamwise flow magnitudes were also slow throughout the measurement of W2019 (Figure 8 and Figure S10): the highest flow clusters were then only 0.3-0.4 m/s and slowest flow velocities were 0-0.2 m/s. These were reflected in the statistically significant autocorrelation areas. Similarly, the streamwise velocities of W2020 were between 0-0.2 m/s most of the time, and only a few times were close to 0.4 m/s ( Figure S10). These low flow magnitudes could have caused the fact that most of the cells had statistically significant autocorrelation with neighboring cells during that measurement time (Figure 8).

Measured Sediment Transport
The bedload was greatest at the CS4 measurement location on W2020, and second greatest at CS1. At CS4 the bedload was from 844.3 to 2697.4 mg/(m·s), depending on the measurement. At CS1 the bedload was from 76.8 to 391.1 mg/(m·s). Minimal transport was at CS2 and CS3, as their every sample was below 50 mg/(m·s) on 6 February 2020. Only one sample, CS2 HS1 was smaller (12.8 mg/(m·s)) than the one measured in 2017 February (22.8 mg/(m·s)). The sample "CS1 HS3" was medium sand, moderately well sorted, coarse skewed, and mesokurtic, according to the calculations done with Gradistat software [25]. The "CS4 HS1" sample was defined as coarse sand, moderately well sorted, symmetrical, and mesokurtic.
As mentioned in Section 3, each bedload sample was measured within 6 min, but not at the same time with the ADCP measurements to avoid disturbance to the flow measurement data. While it was not possible to define the effect of exact macro-turbulent flow wedges during separate bedload measurements, it was possible to show variability between the samples, which indicates differing flow conditions at each bedload measurement time.
Total suspended solids varied 0-0.5 mg/L on W2019, were 0 mg/L on A2018 and on W2018, varied 6.5-7.5 mg/L on 25 May 2018, and were again 0 mg/L on A2017 (Table 6). Therefore, the transport method was mainly bedload at the Pulmanki River during the ADCP measurement times, both during the ice-covered and open-channel conditions. Water color values and turbidity were also low during each measurement time ( Table 6). The color of the water was between 20 and 50 mg/L Pt, and turbidity 0-4 FTU. The samples represented the conditions described by the ADCP data and did not include spring snow-melt discharge peak conditions. Water 2020, 12, 1874 24 of 33 Figure 8. The results of the LISA analysis at CS4 during all measurement times based on streamwise velocity. The S2017 and W2018 results can be seen from Figures 3 and 4. Measurement ID means the sample numbers (n), and therefore also represents seconds. Note that the W2019 and W2020 measurements were shorter in time. The black dots indicate the occurrence times of statistically significant strong high and low flow wedges reaching from the surface to the riverbed, and the grey dots indicate the times when there is weaker wedge (i.e., cluster has broken with a few cells within the vertical water column). W2019 ice-water interface is marked in the figure.  Figures 3 and 4. Measurement ID means the sample numbers (n), and therefore also represents seconds. Note that the W2019 and W2020 measurements were shorter in time. The black dots indicate the occurrence times of statistically significant strong high and low flow wedges reaching from the surface to the riverbed, and the grey dots indicate the times when there is weaker wedge (i.e., cluster has broken with a few cells within the vertical water column). W2019 ice-water interface is marked in the figure.

Near-Bed Velocities: Exceedance of the Critical Velocity Threshold
The near-bed fluctuation at each measurement time was compared against the critical velocity threshold of sediment movement, and to the overall macro-turbulent flow characteristic of those measurement times. Based on the exceedance of the critical velocity threshold of the near-bed layer velocities, most prone conditions for riverbed erosion were at CS1 on S2017 (Figures 9-11), when also the occurrence of the strong macro-turbulent flow wedges was the highest at the inlet of the meander bend (Table 6). Second-and third-highest sensitivity was on W2018 and W2016, respectively, when also pulsating flow was clear and the occurrence of the strong macro-turbulent wedges was its highest of the winters ( Table 6). The streamwise velocity fluctuated much during W2016, S2016, A2016, and W2017, in that the velocities varied above and below the critical threshold of incipient motion of D90 grain sizes. Note that conditions at these times were either ice-covered low flow situations or autumn low flow situations. On W2020 there was great fluctuation above and under the threshold, with overall flow magnitudes lower than those observed on W2016 and W2017. During the winter measurements there were greater low velocity outliers within the boxplot calculations (Figures 10 and 11) than during the open-channel flow situations at CS1, indicating the importance of high flow wedges for the exceedance of the critical sediment transport thresholds. The streamwise velocity of the near-bed layer would not have exceeded the critical velocity to move the D90 particles on W2019 and S2019 flow conditions at CS1. The D50 particle sizes would not have started moving on W2019 and S2019, and the threshold of the erosion would have exceeded on W2020 only occasionally (Figures 9-11).
The streamwise near-bed velocities were also high at CSB, in particular on W2016 and A2017. The critical velocity threshold of D90 value of "CS4 HS1" sediment sample was exceeded within the 25-75 percentile ranges of those measurements ( Figure 11). On W2016 there were few measurement time-steps when D90 would not have moved. On A2016, S2017, and S2019 the D50 and D90 thresholds were exceeded during most of the samples. Fluctuation above and below the threshold was greater than on A2017 and W2016, as during A2017 and W2016 the near-bed velocities were mostly above the threshold. Most fluctuation of the velocities above and below the D50 threshold was on W2017, as the 25-75 percentile region was around the D50 threshold velocity values ( Figure 11). However, on that time, the occurrence of macro-turbulent wedges was low (Table 6). Thus, this indicates that the fluctuation of near-bed velocities at the apex was not that year related to the macro-turbulent flow wedges. The widest ranges of the outliers were on W2019 where the overall median velocity was below the D50 threshold. That same period, similar to W2017, the macro-turbulent wedges were low at this apex area ( Table 6). During that low flow period the fluctuation of the flow magnitude was the greatest for the initiation of the sediment transport. On W2020, the median velocities were below the critical thresholds for the incipient movement, and the outlier ranges were also smaller than on W2019 at CSB. Thus, least prone conditions for erosion were on W2020 at CSB.
At CS4, the critical velocity was exceeded most frequently during A2016, when a great amount of strong macro-turbulent flow wedges had occurred ( Table 5). The D50 threshold was frequently exceeded on W2016, W2017, W2018, S2017, and A2017. These all had multiple macro-turbulent flow wedges, except on W2016 (Table 6). Thus, during most of those periods the flow pulses caused the exceedance of the threshold. However, the velocities measured on W2019 and W2020 basically did not exceed the critical threshold of streamwise velocity and the statistically significant flow wedges had overall slower velocities than in other years. On W2020, when the "CS4 HS1" sample had been measured, the outlier negative (i.e., towards upstream), near-bed velocity values exceeded the D50 grain size threshold. At other measurement locations, there were several time steps (outlier values of the boxplots), where the upstream velocities were high enough that the movement of both D90 and D50 particles of bead load sample "CS4 HS1" would have started (Figures 10 and 11). the flow measurement samples had near-bed velocity values often close to zero. Thus, the sediment particles, whose motion had possibly started further upstream, would not have had velocities low enough for their settlement, as the velocities were predominantly within the transport zone defined by Hjulström [24].    In addition, the velocity threshold for deposition of the largest analyzed D90 grain size particles (0.07 m/s) was exceeded according to the median near-bed flow velocities during most of the measurement times at CS1, CSB, and CS4 locations (Figures 9-11). Only during W2019 and W2020 the flow measurement samples had near-bed velocity values often close to zero. Thus, the sediment particles, whose motion had possibly started further upstream, would not have had velocities low enough for their settlement, as the velocities were predominantly within the transport zone defined by Hjulström [24]. Figure 11. The boxplots of streamwise velocity (m/s) of the deepest cells at CS1, CSB, and CS4 measurement locations during each measurement time. The critical velocity thresholds for the movement of the largest measured D50 and D90 grain sizes are shown. Note that the deepest cells are at different depths depending on the measurement times. The boxplots show the median line, 1st and 3rd quartiles (25 and 75 percentile ranges) within the boxes, maximum observations below upper and lower fences (i.e., the fence is 1.5 times the difference between the 3rd and 1st quartiles), and outliers (dots). In particular, scattering occurred on W2019 at CSB, as there were wide ranges in the outlier values. Figure 11. The boxplots of streamwise velocity (m/s) of the deepest cells at CS1, CSB, and CS4 measurement locations during each measurement time. The critical velocity thresholds for the movement of the largest measured D50 and D90 grain sizes are shown. Note that the deepest cells are at different depths depending on the measurement times. The boxplots show the median line, 1st and 3rd quartiles (25 and 75 percentile ranges) within the boxes, maximum observations below upper and lower fences (i.e., the fence is 1.5 times the difference between the 3rd and 1st quartiles), and outliers (dots). In particular, scattering occurred on W2019 at CSB, as there were wide ranges in the outlier values.

Discussion
Similarly to Lotsari et al. [2] who demonstrated based on cross-sectional measurements conducted at the Pulmanki River in 2014-2015, the high velocity core was located in the middle of the water column in ice-covered flow conditions of 2016-2020, whereas in open-channel conditions the highest velocities are located closer to the water surface than the riverbed. The winter discharges of the Pulmanki River were mostly lower than those observed by Demers et al. [9], and the spring discharges were higher. However, the autumn discharges were in similar magnitude as in the study site of Demers et al. [9]. In comparison, the study site of Demers et al. [9] had ice roughness from low relief wavy forms to high amplitude undulations with scarp faces, thus the ice cover of the Pulmanki River was smoother by comparison. Despite these differences to Demers et al. [9] overall discharge and ice-cover conditions, macro-turbulence was observed at the Pulmanki River. The statistically significant clusters of high and low flow were observed during all measured seasons and locations within a symmetrical meander bend. However, there were also periods where the statistically significant high flow clusters were weaker and were "broken" in the middle of the vertical water column so that these wedges had a short gap between the surface and bottom layers of the water column. Particularly within the apex area, these clusters were sometimes difficult to define and some of the statistically significant flow areas were left outside of the frequency calculations if they occurred predominantly in the bottom half of the water column or top half of the water column. These types of clusters occurred both in open-channel and ice-covered conditions.
According to Sukhodolov et al. [11], periods of bursting flow (i.e., ejection events can be identified as periods when vertical velocity is greater than 0 and streamwise velocity is less than 0) and periods of sweep events occur when streamwise flow is greater than 0 and vertical velocity is less than 0. In addition, it has been observed that the vertical ejections would be able to reach throughout the water column, but sweeps are said to reach up to 0.7-0.8 times the river depth [9]. As the streamwise velocity was mostly positive (downstream) at the Pulmanki River, this indicates that the bursting ejections were not occurring (i.e., vertical flow was directed seldomly towards the surface at Pulmanki River) during the strong streamwise flow wedges. As the sweeps are defined to reach only up to 0.7-0.8 times the river depth [9], these bursting sweeps could influence the weak macro-turbulent flow structures. However, similar to Demers et al. [9], the length of the vertical stripes of coherent flow structures, reaching from the surface to the riverbed, were longer than the sampling frequency. Therefore, these strongest coherent flow structures occurred quasi-steadily and met the criteria of the macro-turbulent flow [9], and not only bursting ejections and sweeps at the channel bed or at the surface air-water/ice-water interfaces [11].
During ice-covered conditions, these statistically significant high and low flow wedges occurred most consistently in the upstream inlet area (CS1), and secondly at the outlet area (CS4). The flow pulsation was significant in the inlet area, where the water depth was also shallow under ice. There was more variability in flow velocities at the apex, and it was not possible to see clear macro-turbulent flow pulsation in these high velocity measurement locations close to the outer bank. Thus, the pattern of flow in the inlet area was similar to Demers et al. [9]. The riffle in the inlet area of the Pulmanki River could explain the macroturbulence, much like Marquis and Roy [7] stated. However, as we did not have measurements from the upstream riffle-pool section (i.e., further upstream of the CS1 measurement location), we cannot state this for sure. Note that during ice-covered flow, only at this inlet area was there a clear layer of downward flow, indicating that there was circulating cells in the vertical direction, in addition to the streamwise macro-turbulent pulsating flow to the downstream direction. This vertical circulation cell could indicate the conditions also observed by Shvidchenko and Pender [6] in their flume study. These authors observed an upward flow in the bottom layer of the water column and a downward flow in the top layer of the water column due to macro-turbulence.
Macro-turbulence was observed during the open-channel flow conditions in autumn and spring. Similar to the ice-covered conditions, the occurrence of the high and low flow pulses was the most consistent at the inlet (CS1) and then secondly at the outlet of a meander bend (CS4). More variable and chaotic nature of the flow velocities was observed in the apex area. Macro-turbulence occurred the most frequently at the inlet area in spring 2017, which had one of the highest discharges of the 5-10 min-long ADCP measurement times. These statistically significant strong (i.e., vertically unbroken) flow wedges reaching from the surface to the near-bed layers occurred 2-3 times per minute, in addition to the occurring weaker (i.e., broken within 5-20% of the depth of the measurement time) statistically significant flow clusters. This is far from the frequency observed by Buffin-Bélanger et al. [5], who observed the more than two second-long large-scale flow structures to occur up to nine times per minute. Similar to the ice-covered period, in spring open-channel high flow conditions (e.g., S2017), the vertical velocity was observed to be downwards at the surface layers within the inlet area (CS1) and at the two following cross-sections closer to the apex (CSA and CS2). Thus, the vertical flow circulation was similar to Shvidchenko and Pender [6] also during the spring open-channel discharge conditions of the Pulmanki River.
Despite the frequency of the strong flow wedges was not as consistent and frequent during other measurement times than on S2017 at the Pulmanki River, the clearest statistically significant high flow wedges had always streamwise velocities at least 0.6-0.7 m/s, independently of the season. One difference between the seasons was that the strong statistically significant flow wedges were temporally longer during the measurement times, when the overall discharge was slower. This was especially the case in the inlet and outlet areas of the meander bend. Thus, our results show similarities to a previous study [9] in the sense that the time and length scales were larger for the ice-covered flow situation than the open-channel flow conditions.
Buffin-Bélanger et al. [5] noticed that the macro-turbulent flow cells were tilted towards downstream, in that the average angle of the front of the wedges was 36 • . Despite the high velocity wedges occasionally being tilted at the Pulmanki River, they were mostly straight from surface to the bottom layers of the flow. Thus, we were not able to see this consistent tilting from the ice-covered and open-channel flow data of ours. Note that the depths of Buffin-Bélanger et al. [5] were mostly less than 30 or 20 cm (i.e., much smaller than the ones of this present study; Tables 2 and 3).
The apex area differed from the inlet and outlet areas in regards to the characteristics of the statistically significant flow clusters during ice-covered low flow (winter), open-channel high flow (spring), and open-channel low flow (autumn) conditions. The macro-turbulent flow wedges were the weakest and temporally shortest (i.e., only 1-2 s long) in the apex region (CS2-CSC). In the inlet areas the 4-5 s long wedges corresponded to the lengths of Demers et al. [9], who stated that the macro-turbulent flow wedges in streamwise velocity direction were in ice-covered conditions 0.4 times the depth. Thus, the apex area of Pulmanki River's bend had temporally shorter macro-turbulent flow structures than Demers at al. [9] observed. During S2017 measurement time, for example, when moving from inlet to the apex region, the numbers of the strong streamwise high and low flow wedges dropped, at the same time the weak wedge numbers increased. In addition, many statistically significant high and low flow clusters appeared often either near the riverbed or near the surface (ice-water/air-water interface) layer at the apex. There were also frequent high flow clusters at the surface layers, and low flow clusters simultaneously near the riverbed, or vice versa. Thus, the statistically significant clusters did not reach through the water column, and this turbulence was possibly caused by the friction effect of both the riverbed and ice-cover and not the actual macro-turbulence. It could have also been partially caused by the centrifugal forces and the friction caused by the outer bank. Unique to many other measurement times were W2016 and W2019 measurements, as during those times the statistically significant flow wedges reaching from the surface to the riverbed were not frequent at the inlet areas (CS1 measurement location), but were more detectable at the apex area (CSB). Due to the less frequent macro-turbulence at the inlet area, we suggest that the turbulence at the apex area was that time caused by local friction and centrifugal forces and not the macro-turbulence.
According to Shvidchenko and Pender [6], these depth-scale eddies are an important turbulence mechanism contributing to sediment transport, as they cause longitudinal troughs and ridges on bed and cause transport of particles along troughs. The data sets showed that the critical velocities for incipient motion were exceeded in both ice-covered and open-channel conditions. The most prone velocity situation for erosion, based on the exceedance of the critical velocity threshold of the near-bed layer streamwise velocities, was at CS1 on S2017 and second-and third-highest sensitivity were on W2018 and W2016, respectively, when variation of high and low near-bed streamwise flow velocities was also clear. The near-bed streamwise velocity fluctuated greatly during W2016, S2016, A2016, and W2017 (i.e., the velocities varied above and below the critical threshold of incipient motion of the sediment grain sizes of Pulmanki River). The macro-turbulence thus enabled the exceeding of the thresholds during these low flow open-channel and ice-covered conditions when the overall discharges were small. However, note that in addition to the existence of the macro-turbulent flow wedges, also local near-bed turbulence due to friction caused the threshold for incipient motion to be exceeded. There were spatial differences between the inlet, apex, and outlet areas, whether the threshold for motion was exceeded. The thresholds for D50 particle sizes to move were exceeded at CS1 during most of the measurement times, second often during different seasons at the apex (CSB), and most seldomly at the outlet area (CS4). However, at CSB the highest fluctuations of near-bed velocities occurred, when there was not as much macro-turbulent flow. This indicates that the small scale turbulence near the channel bed was the cause of the transport at the apex area. In addition, the highest near-bed velocities were not detected at the same measurement time at the inlet, apex, and outlet areas. This could be also due to the different depth and bedform conditions of the river channel related to the measurement times. In any case, when the macro-turbulent wedge frequencies, and the bedload amounts were compared together, the bedload transport was the greatest at CS4 and CS1, which also had the most macro-turbulent wedge occurrences (both in number, and per minute). Thus, similar to the results of previous studies on open-channel flow conditions [5][6][7][8], the macro-turbulent flow is an important contributor to sediment transport, such as that observed at the inlet and outlet of the meander bend at the Pulmanki River. However, our results show that it is especially important during low flow situations when the overall velocities and discharges are not high.

Conclusions
Macro-turbulent flow wedges were observed in every measurement location at the meandering Pulmanki River during ice-covered low flow conditions and open-channel flow conditions. Pressurized winter-time conditions were not a requirement for the development of the macro-turbulent flow.
These macro-turbulent flows occurred most consistently, based on the streamwise velocities, in the inlet area and second most often in the outlet area of the meander bend during both ice-covered and open-channel conditions. Centrifugal forces and local friction due to confining channel and river ice surface could have caused chaotic flow pattern at the apex. Thus, the macro-turbulent flow structures were not clear in the apex area of the analyzed meander bend.
In line with previous studies, the temporal length of the macro-turbulence was longer in conditions of low discharges and shallow depths, than in high discharge conditions. However, the longest statistically significant high flow wedges occurred at inlet and outlet areas and were temporally shorter around the apex. The frequency of the macro-turbulence was the greatest during the spring open-channel conditions. However, similar to all seasons, clearest and strongest macro-turbulent high flow wedges had velocities of at least 0.6-0.7 m/s. Despite near-bed velocity fluctuation and the macro-turbulent flow contributing to the sediment transport during the high discharge periods, during those measurement times the overall velocities were high, and therefore the macro-turbulence was not as critical to the sediment transport as during lower discharge conditions. There were spatial differences between the inlet, apex, and outlet areas, where the critical velocity threshold for sediment movement was exceeded by the near-bed streamwise velocities. The threshold was exceeded often at the inlet, where the statistically significant macro-turbulent flow wedges occurred the most consistently. Thus, the macro-turbulence was especially important for initiation of bed sediment movement during low flow situations (i.e., when the overall velocities were not high). The macro-turbulence enabled the flow velocities to exceed the sediment movement threshold more often than based on solely the turbulence caused by the friction of the riverbed.