1. Introduction
Alpine ski racing, with its varying terrain and slope, gate set-up, and snow conditions, is extraordinarily challenging, requiring a high degree of technical skill and motor control [
1]. These demands are compounded by the high average speeds involved: 54 km/h for Slalom, 65 km/h for Giant Slalom, 86 km/h for Super-G, and 94 km/h for Downhill [
2]. Despite this complexity, which demands sophisticated biomechanical analysis [
3], the only two forces acting that can slow down skiers are friction between the skis and the snow and aerodynamic drag [
4,
5,
6,
7] and skiing performance is improved by virtually any lessening of the former [
8].
The skier’s posture, i.e., his/her area of frontal exposure, is the major factor determinant of aerodynamic drag, with less exposure resulting in higher speed and reduced overall race time [
3]. Indeed, in connection with downhill competitions, differences in aerodynamic drag account for almost half of the difference in the finishing times of slower and faster skiers [
9]. However, in giant slaloms, aerodynamic drag is not thought to be a major concern, accounting for only 15% of the total energy loss per turn [
6]. At the same time, Meyer and colleagues [
5] found that a giant slalom skier can reduce his/her energy dissipation due to aerodynamic drag by altering posture from dynamic to compact and, indeed, in this discipline aerodynamic drag accounts for 5–28% of all energy losses [
6]. Clearly, the downhill and super-giant slalom, with their considerably higher speeds, involve more aerodynamic drag than slalom and giant slalom [
10,
11]. In light of the fact that the difference between the finishing times of elite alpine skiers competing for medals is often on the order of hundredths of a second [
12], even very minor optimization of aerodynamic drag can be decisive.
To measure aerodynamic drag in connection with alpine ski racing, several experimental approaches have been developed for both speed [
13,
14,
15] and technical disciplines [
5,
6,
16]. Assessment of the coefficient of aerodynamic drag (C
d) and cross-sectional area (S), as well as C
d∙S, the product of these two, in wind tunnels have attempted to optimize body posture, showing that height, joint angles, and positioning of the arms and head all influence the aerodynamic drag [
5,
6,
16]. Moreover, C
d∙S for the technical disciplines ranged from 0.63–0.66 in the high, 0.51–0.55 in the middle, and 0.23–0.24 in the tuck positions, which is generally lower than the corresponding values reported for the speed disciplines (0.15–0.4) [
8,
13].
Investigation of different body postures in a wind tunnel demonstrated that aerodynamic drag was 41% less with the arms in a tuck position than stretched out (as when jumping) [
17]. Although the ski poles are positioned more or less transverse to the wind flow during the performance of all of the technical disciplines, including during the intensive turning involved in the speed disciplines (see, for example,
Figure 1) and the poles (typically 120–135 cm in length for male skiers) are considerably longer than the arms and have baskets, their impact on aerodynamic drag remains largely unexplored. Accordingly, the contribution of the ski poles to total aerodynamic drag in connection with various postures during alpine skiing at different speeds was assessed here in a wind tunnel. More specifically, we evaluated (i) the contribution of the cross-sectional area of the poles (S
p) to the skier’s total cross-sectional frontal area (S); (ii) intra-subject variation in this respect; and (iii) the influence of S
p on the coefficient of aerodynamic drag (C
d) and the product C
d∙S.
2. Materials and Methods
The data analyzed in the present study were taken from two previous investigations, where the collection of these data is described in detail [
6,
16]. Therefore, below we describe only essential aspects of our methodology, concentrating on details of data processing that differ from those of earlier reports.
2.1. Participants
The 10 male skiers (body mass: 80.8 ± 7.0 kg; height: 1.81 ± 0.05 m, means ± standard deviations) involved in this study were all members of the Swedish alpine ski team and signed their written informed consent prior to participation. The current study involving the reuse of some data was approved by the Commission for Ethical Issues in Sports in Ljubljana, Slovenia (Ref. No. 4/2023).
2.2. Instruments, the Experimental Design and Data Collection
Data were collected while the subjects were standing inside a closed-loop wind tunnel (12.5 m long, 2.7 m wide, 1.8 m high, contraction ratio 1:4.23) equipped with a 220 kW centrifugal blower that can produce wind speeds up to 30 m/s. The standard deviations associated with the measurements were 0.14 m/s and 1.8 N for wind speed and aerodynamic drag force, respectively. Aerodynamic drag was measured with a force platform (model 9286AA, Kistler Instrument Corp., Winterthur, Switzerland), which was set to zero prior to each test under windless conditions. The cross-sectional area of the poles (S
p) and the total cross-sectional area of the skier (S) were assessed by application of click-based interactive segmentation [
18] to the images taken in the coronal plane from behind with an HDR-HC7 camcorder purchased from the Sony Corporation (Tokyo, Japan) (
Figure 1). Fifteen images (5 in each of the three positions) were re-processed to determine the reproducibility of the measurement of cross-sectional area, resulting in a difference of 0.31 ± 0.24% (mean ± SD).
During the wind tunnel test, the participants wore their own suitably sized giant slalom alpine ski racing gear, i.e., racing suit, gloves, helmet, goggles, ski boots, and poles, enabling them to familiarize themselves with the positioning. The conical poles were round in cross-section and had round baskets. The stance width (i.e., the fixed distance between the ski bindings) was adjusted individually. Upon receiving a given visual signal, participants assumed one of three different typical positions (tuck, middle, high) five times in a randomized order at the respective wind speeds: 40, 60, and 80 km/h. Prior to this testing, these three typical positions were recorded while two of our subjects adopted them when skiing down a training course for Giant Slalom. The aerodynamic drag (Fd) was considered to be the 5-s average of the force in the direction in which the wind was blowing, once this signal had stabilized completely.
2.3. Calculations
As is standard practice in studies involving wind tunnels, the skier’s cross-sectional area (S) multiplied by the coefficient of aerodynamic drag (C
d) was obtained from the Rayleigh drag equation C
d∙S = 2∙
Fd/(
ρ·
V2). In this equation,
Fd represents the aerodynamic drag,
V the wind speed, and
ρ the density of the air. The cross-sectional areas, both total (S) and of the ski poles (S
p) alone, were derived from image processing as described above. The value of S obtained by interactive segmentation was in pixels and special calibration marks in the wind tunnel were utilized to convert these pixels into square meters. Division of C
d∙S by S provided C
d. Since in the combined tuck and high positions, the skier obstructed roughly 6–12% of the airflow in the wind tunnel, C
d was corrected for blockage as proposed by Elfmark and colleagues [
19].
2.4. Statistical Analyses
The cross-sectional areas S and Sp are presented as means and standard deviations, as well as the percentage of S accounted for by Sp (i.e., (Sp/S)∙100). The individual variation in Sp for each skier was expressed as the difference between the maximal and minimal measured values of Sp and Sp/S (min–max range), together with the coefficient of variance (CV) for Sp.
To compare the cross-sectional area with the different positions and speeds, three Bayesian linear multilevel models were applied to S, Sp and Sp/S, respectively. In these models, the speed and position were fixed factors, and a random intercept was utilized for each athlete. To examine the effect of Sp on Cd and Cd∙S, linear Bayesian multilevel models were fitted for each position, Sp, and body cross-sectional area, again employing a random intercept for each skier. Weakly informative priors were included on standardized values in each model, with intercept and coefficients of ~N(0, 2.5) and residuals of ~exp(1). With each model, 4 chains of 4000 warm-up and 4000 sampling iterations were performed, for a total of 16,000 sampling iterations. There were no diverging transitions, the minimal effective sample size to total sample size ratio of >0.26, and maximal Rhat value < 1.0016.
The highest-density intervals were utilized to calculate the confidence intervals (CI). In addition, the probabilities of direction (pd) (values from 0.5–1), which indicate how probable it is that a positive effect is actually positive and a negative effect actually negative, were determined. The pd can be interpreted as pd < 0.95 unclear, pd > 0.95 possibly existing, pd > 0.97 likely existing, pd > 0.99 probably existing, and pd > 0.999 certainly [
20]. The R 4.2.2 software [
21] was employed to carry out all analyses and rstanarm 2.21.3 [
22] for fitting models.
3. Results
The cross-sectional areas of the skiers’ bodies (S) and of their poles (S
p) in the different positions and at different wind speeds are shown in
Table 1, together with the intra-skier variation between trials.
Figure 2 illustrates the cross-sectional area of the poles for each skier during all of the trials and the images of Skier 8 presented in
Figure 3 highlight the extent to which the positioning of the poles and corresponding S
p obtained at 40, 60, and 80 km/h varied. Also noteworthy is the asymmetric positioning of the left and right poles (which was common for this particular skier), as well as the variation in the height of the shoulders and head.
3.1. The Cross-Sectional Area of the Ski Poles
The skiers demonstrated a somewhat smaller Sp in the high and middle positions compared to the tuck position, as well as in the middle versus high position (High–Middle = −0.0009 m2, 95% CI [−0.0017, −0.0001 m2], pd = 0.987; High–Tuck = 0.0186 m2, 95% CI [0.0178, 0.0194 m2], pd = 1.000; Middle–Tuck = 0.0195 m2, 95% CI [0.0186, 0.0202 m2], pd = 1.000). The percentage of the total cross-sectional area occupied by the poles differed between the positions, being greatest in the middle position, followed by the high position (High–Middle = −0.5%, 95% CI [−0.7, −0.3%], pd = 1.000; High–Tuck = 2.1%, 95% CI [1.9, 2.3%], pd = 1.000; Middle–Tuck = 2.6%, 95% CI [2.4, 2.7%], pd = 1.000).
In the high and middle positions, Sp at a wind speed of 40 km/h was probably higher than at 60 km/h and was certainly higher than at 80 km/h (40–60 km/h = 0.0013 m2, 95% CI [0.0003, 0.0022 m2], pd = 0.994; 40–80 km/h = 0.0018 m2, 95% CI [0.0008, 0.0028 m2], pd = 1.000; 60–80 km/h = 0.0005 m2, 95% CI [−0.0005, 0.0015 m2], pd = 0.849). At the same time, in the tuck position, it was unclear whether wind speed influenced Sp (40–60 km/h = −0.0001 m2, 95% CI [−0.0015, 0.0013 m2], pd = 0.561; 40–80 km/h = −0.0008 m2, 95% CI [−0.0021, 0.0006 m2], pd = 0.863; 60–80 km/h = −0.0007 m2, 95% CI [−0.0020, 0.0007 m2], pd = 0.830).
In the high and middle positions, the cross-sectional area of the poles comprised a somewhat higher percentage of S at a wind speed of 40 km/h than at the two higher speeds (40–60 km/h = 0.3%, 95% CI [0.1, 0.5%], pd = 0.999; 40–80 km/h = 0.4%, 95% CI [0.2, 0.6%], pd = 1.000; 60–80 km/h = 0.0%, 95% CI [−0.2, 0.2%], pd = 0.667). At the same time, the influence of wind speed on the Sp in the tuck position was again unclear (40–60 km/h = 0.0%, 95% CI [−0.3, 0.2%], pd = 0.587; 40–80 km/h = −0.2%, 95% CI [−0.5, 0.1%], pd = 0.928; 60–80 km/h = −0.2%, 95% CI [−0.5, 0.1%], pd = 0.894).
3.2. Effects of the Cross-Sectional Area of the Poles on Cd and Cd·S
The intra-athlete influence of a 0.01 m
2 increase in the cross-sectional area of the ski poles on the C
d and C
d · S with all wind speeds (controlling for S) can be seen in
Table 2. In the case of the high position, S
p had an unclear effect on C
d, while exerting a likely positive effect on C
d·S. With the middle position, there were likely positive effects on both C
d and C
d·S. Finally, in the tuck position, there was a probable negative effect on C
d and possible negative effect on C
d·S. The general relationship between the cross-sectional area of the ski poles and C
d and C
d·S is illustrated in
Figure 4.
Table 3 documents the effects of increasing S
p by the average range for all the skiers combined at each position on C
d and C
d·S (m
2).
4. Discussion
The major findings from this investigation on the contribution of the ski poles to total aerodynamic drag during alpine skiing were as follows: (i) the average cross-sectional area of the ski poles (Sp) ranged from 0.0067 in the tuck position to 0.0262 m2 in the middle position; (ii) these areas accounted for 2.2–4.8% of the total cross-sectional area; (iii) the contribution of the ski poles to the total cross-sectional area was dependent on wind speed in the high and middle positions; (iv) at 80 km/h, the inter-subject variations in this contribution ranged from 0.0018 m2 (27.6%) in the tuck position to 0.0072 m2 (30.5%) in the high position; and, finally, (v) Sp exerted a likely effect on the coefficient of drag (Cd) and Cd∙S (product of coefficient of drag and total cross-sectional area).
This study involved 10 elite skiers who assumed three different standard positions five times at each of three different wind speeds in a wind tunnel. We did not standardize these conditions, as was done in certain previous studies [
8,
13,
15], because the variations in the positions of the skiers and their poles, as well as in the wind speed, allowed us to examine the influence of these different conditions on the impact of the poles of aerodynamic drag. Moreover, the posture of a skier during a race varies considerably, often deviating from the basic positions, as he/she adapts to the varying terrain, snow conditions, and gate positions [
23].
Our current findings reveal that the S
p could account for as much as 5.0 ± 0.8% of the total frontal cross-sectional area (
Table 1), making the greatest absolute contribution in the high position and the most pronounced relative contribution in the middle position. Furthermore, the intra-skier variations in S
p were relatively extensive (17.1–34.8%), more so apparently at the two highest speeds examined. This variation for all of the skiers combined was even more pronounced (35.6–63.0%), as were the absolute differences between skiers (
Figure 2). Knowing that aerodynamic drag is one of the two resistive forces in alpine skiing [
4,
7,
24,
25] contributing to energy dissipation [
6,
26], these observations indicate that the ski poles may be responsible for variations in total energy loss as great as 1% (considering that maximally 28% of the energy loss during giant slalom is due to aerodynamic drag [Supej et al., 2013] [
6], a 5.8% relative cross-sectional area of the poles, and 63% variation), with an average of approximately 0.3% (considering that, on average, 15% of the energy loss during giant slalom is due to aerodynamic drag [Supej et al., 2013] [
6], a 4.5% relative cross-sectional area of the poles and 40% variation (the latter two values being approximately the mean values for the middle and high positions). No more than a fraction of a second often separates the finishing times of the gold and silver medallists [
12], e.g., 0.01 s or 0.015% in the case of the Super Giant Slalom for men in the recent 2023 Alpine Skiing World Championships.
As illustrated in
Figure 3, the position in which the skiers held their poles varied considerably, which, since there is an obvious geometrical relationship between the cross-sectional area and inclination of the poles in the airstream, influenced the values of the cross-sectional area obtained. It should be noted that in the example shown in this figure the S
p for the poles declined with increasing wind speed, whereas the graph in
Figure 2 shows that the opposite was more often observed for this particular Skier 8. The variation in the positioning of the poles between skiers was even greater than this variation for one and the same subject (
Table 1) and it would be worthwhile investigating further the pole position (for example, more or less parallel to the direction of skiing, i.e., to the wind flow) that minimizes drag optimally. However, this potential improvement in performance indicated here is considerably lower than the ~40% reduction in aerodynamic drag that could be achieved by training the arms according to previous reports [
17,
24] and future investigation should take both the arms and poles into consideration.
The smallest relative contribution of S
p to aerodynamic drag was associated with the tuck position, even though this is the position where the total frontal cross-sectional area of the skier is by far the smallest (
Table 1). Perhaps the fact that the poles are held much more parallel to the ground in the tuck position reduces their cross-sectional area to approximately one-quarter of that in the high and middle positions. To further reduce drag, some skiers who specialize in the speed disciplines use folded poles that allow them to align the shape of their poles to that of their bodies. However, even in such a situation, the poles can give rise to strong vortices that increase aerodynamic drag [
27].
As might be expected, we observed a probable positive effect of S
p on C
d in both the middle and high positions and on C
d∙S in the high position, such that a change in S
p of 0.01 m
2 might alter C
d and C
d∙S by as much as 0.02 m
2 (
Table 2). More surprising were the negative effects of S
p on C
d and C
d∙S in the tuck position, implying that these values both decrease with increasing S
p. The finding is probably due to the fact that the less favorable positioning of the poles (larger S
p) allows more favorable aerodynamic positioning of the skier. For example, if the skier inclines the poles slightly more forwards (tips higher than the grips), this would allow a lower, more aerodynamic positioning of the head and torso, which contribute strongly to the total aerodynamic drag [
8,
28].
Another interesting observation here was that both Sp and its relative contribution to the total aerodynamic drag were on average higher in the middle and high positions at 40 km/h than at 60 and 80 km/h, even though the skiers were instructed to assume the same position at all three wind speeds. Whether the skiers instinctively changed their position somewhat to reduce the aerodynamic drag (i.e., by positioning the poles in a manner that reduced their cross-sectional area) and/or whether this was a consequence of the drag on the poles forcing them into a more inclined position is not possible to know and it is difficult to predict whether similar behavior would occur during actual skiing.
The limitations of the current investigation are similar to those associated with a previous study based on the same measurements [
16], as well as other experiments conducted in a wind tunnel [
5,
8,
13]. In short, despite their familiarity with skiing positions which they assume almost automatically, elite skiers may not adopt exactly the same positions in a wind tunnel as on snow. Moreover, to maintain balance and/or adapt to the prevailing terrain and snow conditions, elite skiers switch between many different positions during a race, whereas here only the three most typical positions were considered. The lower shell of the ski boots was elevated approximately 30 mm by the ski binding and was probably affected by the boundary layer of air near the surfaces of the tunnel, due to friction between the moving air stream and solid surfaces. However, most importantly, the poles with their baskets were high enough to assume they were in the free air stream as shown by an additional 2D kinematic analysis of the positions of the pole baskets, which revealed that their heights (left and right combined) were 28.4 ± 13.4, 17.1 ± 11.2, and 93.0 ± 12.1 cm (mean ± SD) for the high, middle, and tuck positions, respectively. To minimize the error of the drag coefficient, we employed a blocking correction with an error of 2.2%, obtained for alpine skiers tested in the same wind tunnel [
19].
Last, but not least, derivation of the cross-sectional areas of the skier and poles through click-based interactive segmentation of the video recordings may not be entirely accurate. In the case of the poles in particular, the pixel size of the video recordings plays an important role. Thus, since the diameters of the conical poles ranged from ~10 to ~19 mm, the use of a pixel size of appropriately 1 mm could introduce an error of ~10%, even though the reproducibility of the determination of the cross-sectional area by image processing was only 0.31%. Since this error can be assumed to be of the same proportion in connection with all measurements, measurements of the absolute values of the cross-sectional area of the poles and their contribution to the total cross-sectional area will be affected most.