Physiological Responses to Rifle Carriage During Roller-Skiing in Elite Biathletes

Purpose: This study aimed to investigate the physiological factors affected by rifle carriage during biathlon skiing performance, as well as the sex differences associated with rifle carriage. Methods: Seventeen national- and international-level biathletes (nine females and eight males; age 23.0 ± 3.3 years, V.O2max 59.4 ± 7.6 mL.kg–1.min–1) performed a submaximal incremental test and a maximal time-trial (TT) using treadmill roller-skiing (gear 3, skating technique) on two occasions separated by at least 48 h. One condition involved carrying the rifle on the back (WR) and the other condition no rifle (NR) and the tests were randomized. Submaximal V.O2, skiing speed at 4 mmol.L–1 of blood lactate (speed@4mmol), gross efficiency (GE), aerobic (MRae), and anaerobic (MRan) metabolic rates, and V.O2max were determined. Results: Submaximal V.O2 (at all intensities) and GE (16.7 ± 0.9 vs. 16.5 ± 1.1%) were higher for WR compared to NR (p < 0.05), while speed@4mmol was lower (3.1 ± 0.4 vs. 3.3 ± 0.5 m.s–1, p = 0.040). TT performance was improved (4.6 ± 0.4 vs. 4.3 ± 0.4 m.s–1, p < 0.001) and MRan was higher (31.3 ± 8.0 vs. 27.5 ± 6.5 kJ.min–1, p < 0.01) for NR compared to WR, with no difference in V.O2max or MRae. For skiing WR, TT performance was correlated to speed@4mmol (r = 0.81, p < 0.001), MRan (r = 0.65, p < 0.01), V.O2max (r = 0.51, p < 0.05), and relative muscle (r = 0.67, p < 0.01) and fat (r = −0.67, p < 0.01) masses. Speed@4mmol together with MRan explained more than 80% of the variation in TT performance (WR 84%, NR 81%). Despite a higher relative mass of the rifle in females compared with males (5.6 ± 0.4 vs. 5.0 ± 0.4% of body mass, p = 0.012), there were no sex differences associated with rifle carriage measured as absolute or relative differences. Conclusion: Rifle carriage in biathlon skiing led to significantly higher physiological demands during submaximal exercise and reduced performance during maximal treadmill roller-skiing compared to NR for both sexes. The most important variables for performance in biathlon treadmill skiing seem to be speed@4mmol combined with MRan, both of which were lower for WR compared to NR. To improve skiing performance in biathlon, improving speed at 4 mmol.L–1 of blood lactate and anaerobic energy delivery while carrying the rifle are recommended.


INTRODUCTION
Biathlon is an intermittent endurance sport that combines rifle shooting and cross-country (XC) skiing while carrying a rifle. A biathlon competition consists of 6-20 km of XC skiing divided into three or five laps each lasting ∼3-9 min, depending on the competition discipline, and between each lap athletes shoot five shots in a prone or standing position. While few scientific studies have investigated the physiological, biomechanical, technical and tactical demands of biathlon (reviewed in Laaksonen et al., 2018), more have examined XC skiing events (i.e., competing without a rifle). Regarding the physiological demands, biathlon has multiple similarities to XC skiing, especially sprint XC skiing, which also involves repeated bouts of skiing (up to 4 × 2-4 min). However, the rest periods between the bouts are significantly longer in sprint XC races (i.e., 20-45 min) (FIS, 2018). Thus, reported findings on physiological responses in sprint XC skiing studies may not be directly transferrable to skiing performance in biathlon.
As in other endurance sports, maximal oxygen uptake (VO 2max ), lactate threshold, efficiency of movement and anaerobic energy delivery all have an important role in biathlon skiing performance (Laaksonen et al., 2018). Indeed,VO 2max values exceeding 80 mL.kg.min −1 have been reported for male biathletes, with approximately 15% lower values for females (Tønnessen et al., 2015). Moreover, bothVO 2max and lactate threshold have been shown to correlate with biathlon skiing performance (Rundell, 1995;Rundell and Bacharach, 1995). While gross efficiency (GE) has been shown to correlate well with XC skiing performance (Sandbakk et al., 2013;Andersson et al., 2017), no previous studies have investigated the importance of GE in biathlon when skiing with a rifle. In addition, although performance in biathlon is predominantly dependent on large contributions from aerobic energy delivery, anaerobic energy contributions are also likely to play an important role. This is the case for XC skiing (Stöggl et al., 2007;Andersson et al., 2019), where it has been estimated that ∼ 17-26% of the total energy delivery during maximal XC skiing lasting ∼ 2-4 min is provided by anaerobically-derived energy (Losnegard et al., 2012;McGawley and Holmberg, 2014;Andersson et al., 2016Andersson et al., , 2017. Unlike in XC ski racing, biathletes carry a rifle weighing >3.5 kg on their back while skiing. It has been estimated that the energy cost could be ∼7% higher while carrying an extra load of ∼4 kg while skiing (Fredrick, 1987). Previous studies investigating physiological responses in biathlon have reported a 2-8% higher oxygen uptake (VO 2 ) and a 6-11% higher ventilation rate (V E ) among male and female biathletes when carrying a rifle during submaximal roller-skiing compared to no rifle (Rundell and Szmedra, 1998;Stöggl et al., 2015). These differences were also reported to be greater for the female participants compared to the males, which may be due to the higher additional load relative to body mass in females. The two aforementioned previous studies used skiing speeds of between 2.2 and 3.0 m.s −1 . However, mean skiing speed for the ten best athletes during a biathlon World Cup sprint competition is currently around 6.3 and 7.2 m.s −1 for females and males, respectively (including both uphill, downhill, and flat terrain), and this speed has increased by ∼7% over the last 15 years (Laaksonen et al., 2018).
The complex nature of intermittent XC skiing and rifle carriage in biathlon makes the sport unique in its specific demands. With only two studies having investigated the effects of rifle carriage on a limited number of physiological responses in biathlon (Rundell and Szmedra, 1998;Stöggl et al., 2015), more research is needed to fully understand the physiological determinants of biathlon in order to improve sport-specific training and preparation for competition. Therefore, this study aimed to investigate the physiological factors affected by rifle carriage during biathlon skiing performance, as well as the sex differences associated with rifle carriage. This knowledge will help biathletes and coaches to plan and monitor training. It was hypothesized that the physiological responses would be increased when carrying a rifle, and that this effect would be larger in female biathletes compared to males.

Participants
Seventeen well-trained Swedish biathletes (nine females and eight males), competing at national and international levels, volunteered to participate in the study ( Table 1). All participants were fully informed of the nature of the study before providing written consent to participate. The study was pre-approved by the Regional Ethical Review Board in Umeå, Sweden (#2018-122-31M) and was performed in accordance with the Declaration of Helsinki.

Study Overview
Participants visited the laboratory for testing on two separate occasions during the preparation phase of the season (October), with 2-6 days (and a minimum of 48 h) separating the visits. The two testing sessions were performed in a randomized order and involved the same submaximal and maximal exercise protocol, one test skate roller-skiing while carrying the rifle on the back (WR) and one test without the rifle (NR). All participants used the same rifle during testing, which weighed 3.5 kg, and their own biathlon harness weighing 0.2 kg. The testing was performed on a treadmill using gear 3, which all participants were highly familiar with from previous testing and training sessions.

Testing Procedures
The skiing tests were carried out on a motor-driven treadmill (belt dimensions 3.3 × 2.5 m; Rodby Innovation AB, Vänge, Sweden). Participants wore a safety harness around their waist, which was suspended from the ceiling and connected to an emergency brake stopping the treadmill within 1 s in the case of a fall. All participants used Pro-Ski S2 roller skis (Sterners, Dala-Järna, Sweden) equipped with NNN (Rottefella, Klockarstua, Norway, rolling resistance 0.0224) or SNS (Salomon, Annecy, France, rolling resistance 0.0216) bindings. The roller-skis were pre-warmed for at least 60 min in a heating box before all tests to minimize the variation in rolling resistance (µ R ), which was determined as described previously by Ainegren et al. (2008). Expired air was sampled at 10-s intervals during submaximal and maximal exercise using an AMIS 2001 metabolic system (model C, Innovision A/S, Odense, Denmark), which was calibrated before each testing session using a 3-L syringe and a calibration gas with a known mixture of 16% O 2 and 4.5% CO 2 . On arrival at the laboratory participants rested in a supine position for 10 min before a capillary blood sample was collected from the fingertip for the subsequent analysis of resting blood lactate concentration (Biosen 5140, EKF Diagnostic GmbH, Magdeburg, Germany). Body mass and height were also measured (Seca 764, Hamburg, Germany). Thereafter, participants performed a standardized warm-up by roller-skiing for 6 min at the same workload as for the first submaximal exercise level during the subsequent test. The submaximal test consisted of 3-5 levels each lasting 4 min where the gradient was constant (3.5 • for the females, 4.5 • for males) and the speed increased by 0.56 m.s −1 (2 km.h −1 ) for every workload. The female skiers commenced at a velocity of 1.94 m.s −1 (7 km.h −1 ), while the males commenced at 2.22 m.s −1 (8 km.h −1 ). Expired air was analyzed from measurements taken during the last 30 s of each 4-min level and heart rate (HR) was recorded throughout the submaximal test (Polar S810, Polar Electro Oy, Kempele, Finland). Between each level, the treadmill was stopped for 60 s for fingertip blood sampling and to record the rating of perceived exertion (RPE; Borg-scale 6-20) for breathing, arms (RPE arm ) and legs (RPE leg ). The submaximal test was terminated when the respiratory exchange ratio (RER) was >1.00, theV E/V O 2 ratio was >30 and HR was >90% of the known maximal HR reported by the participant.
After the submaximal test participants recovered actively for 4 min (roller-skiing at the first submaximal workload) then rested passively for 5 min. Participants then completed 5 min of treadmill roller-skiing as a re-warm-up, including three sprints lasting 15 s, before commencing a maximal time trial (TT). The TT simulated competition skiing conditions, whereby the females completed 900 m at a gradient of 3.5 • and the males completed 1000 m at 4.5 • as fast as possible. The participants were able to adjust the speed by moving forward and backward on the treadmill, with a system consisting of two laser beams detecting the position of the participant (Swarén et al., 2013). The speed of the treadmill increased by 0.19 m.s −1 .s −1 or decreased by 0.11 m.s −1 .s −1 as the participant moved to the front or rear of the treadmill, respectively. During the TT, expired air and HR were measured throughout and a fingertip blood sample was collected 2 min after the test for the analysis of maximal blood lactate concentration. Performance was defined as mean speed during the TT.
Within 1-3 days of the roller-skiing tests body composition was assessed using iDXA (Encore 2007, Version 11.4, General Electric Company, Madison, WI, United States) in order to measure total body mass, as well as relative and absolute muscle and fat masses. The participants were tested in a fasted state (in the morning before breakfast) and lay in a supine position during the scan, which lasted 7-8 min. When analyzing the iDXA data, the body was also divided into upper-and lower-body sections. Arms, legs, and trunk were also analyzed separately for body composition and the length of the different body segments were measured. The length of the upper-body was measured without the head.

Calculations
Respiratory values (VO 2 ,VCO 2 ,V E , and RER) and HR were calculated as the mean of the last 30 s of each workload during the submaximal test. Speed, HR andVO 2 corresponding to a blood lactate concentration of 4 mmol.L −1 (speed @4mmol , HR @4mmol , andVO 2@4mmol , respectively) were calculated from the individual exponential relationships between blood lactate concentration and speed,VO 2 and HR. Heart rate,VO 2 ,VCO 2 , andV E were recorded throughout the TT and the highest consecutive 30-s values were reported as maximal values. Power output (PO) during roller-skiing was calculated as the power exerted to elevate the total mass against gravity and to overcome the rolling resistance: where M tot = body mass (kg) including equipment (also the rifle, if used), g = gravitational acceleration (m.s −2 ), v = velocity (m.s −1 ), µ R = rolling resistance and α = treadmill incline ( • ) (Andersson et al., 2016). The aerobic metabolic rate (MR ae ) was determined from the meanVO 2 and RER during the last minute of the highest submaximal workload with a RER < 1.00: where E GROSS = gross energy expenditure, given by: (3) Gross efficiency (GE) was calculated from the submaximal workload closest to, but not above, RER = 1.00: The total metabolic demand (MR DEMAND ) of the TT was calculated using the GE method: where PO TT was calculated using Eq. 1 with the mean speed for the TT. The MR ae for the TT was calculated using Eqs. 2 and 3 and the meanVO 2 for the entire TT, assuming RER = 1.0. Anaerobic metabolic rate (MR an ) was then calculated as follows: E GROSS was calculated for the first three submaximal levels using values forVO 2 (L.min −1 ) and RER values ≤ 1.00 using Eq. 3.

Statistical Analyses
All statistical tests were processed using SPSS 24.0 Software (SPSS, Inc., Chicago, IL, United States). Shapiro-Wilk tests showed a normal distribution for all variables for the group in total. For analyses separated for sex, non-parametric tests were used due to the small sample sizes. For the submaximal test, a two-way ANOVA with repeated measures (skiing condition [WR/NR] × submaximal level) was used to investigate differences in respiratory variables, HR and energy system contributions. If a significant difference was observed, repeated measures one-way ANOVA with a Bonferroni post hoc test was used to identify specific differences for the group in total, while a Friedman's test followed by a post hoc Wilcoxon signed rank test was used on RPE and when the sexes were tested separately. Paired samples T-tests (for the group in total) and Wilcoxon Signed Rank tests (for females and males separately) were used to determine the differences between WR and NR for respiratory variables (VO 2, RER, andV E ), HR and blood lactate during the maximal TT. Responses to rifle carriage were calculated as the absolute and relative differences between WR and NR and were compared between the sexes using a Mann-Whitney U test. Effect sizes (ES) were calculated using Cohen's d with the criteria 0.2 = small effect, 0.5 = moderate effect, and 0.8 = large effect. To analyze the relationship between TT performance (dependent variable) and the independent variables for endurance performance (speed @4mmol ,VO 2max , GE, MR an , MR ae , muscle mass, and fat mass), Pearson's product moment correlation tests were conducted for WR. When separating the group by sex, Spearman's Rho was used due to the small sample size. In addition, a multiple regression analysis was conducted with speed @4mmol ,VO 2max , GE, MR an and muscle mass as independent variables to determine the explained variation in TT performance. The TT was also divided into four quarters according to distance covered (i.e., 225 and 250 m for the females and males, respectively) in order to analyze the relationship between aerobic and anaerobic energy contributions during different sections of the TT. Mean PO TT , speed, MR ae and MR an were compared using a two-way repeated measures ANOVA (skiing condition [WR/NR] × quarter), followed by a oneway repeated measures ANOVA with a Bonferroni post hoc test (for the group in total) or a Friedman's test followed by a post hoc Wilcoxon signed rank test (for the females and males separately). The level of statistical significance was set at α ≤ 0.05 and data are presented as mean ± standard deviation (SD), except for RPE, where data are presented as median and interquartile range (IQR).

Submaximal Test
Responses during submaximal exercise for WR and NR are presented in Table 2. All participants (n = 17) performed the first three levels of the submaximal test for both WR and NR, while only 14 participants (eight females, six males) performed four submaximal levels for NR and 12 of those (seven females, five males) performed four levels for WR.
TheVO 2 was 3-5% higher for WR than for NR during submaximal levels 1-4. TheVCO 2 was 4-6% higher andV E was 7-10% higher, while there were no differences in HR or RER. The blood lactate concentration was higher for WR compared to NR only for the fourth submaximal level. The only difference in RPE between WR and NR was for submaximal levels 2 (RPE arm 13 (12-15) vs. 12 (11-13), p < 0.05) and 3 (RPE leg 15 (14-15) vs. 14 (12-15), p < 0.05). The individual values for submaximalVO 2 , HR and blood lactate for WR and NR are shown in Figures 1A-F. For the group in total (both sexes), E GROSS was higher during WR compared to NR for all submaximal levels. GE was also higher during WR compared to NR for the group in total (16.7 ± 0.9 vs. 16.5 ± 1.1%, p < 0.05, ES = 0.54). No interaction effects were found between skiing condition (WR vs. NR) and submaximal level for any of the variables.
While the mass of the rifle in relation to body mass was greater for the females than the males (5.6 ± 0.4 vs. 5.0 ± 0.4% of body mass, respectively; p = 0.012), there were no sex differences in the responses to rifle carriage (i.e., absolute or relative differences between WR and NR) for HR,VO 2,VE , RER,VO 2@4mmol , HR @4mmol , speed @4mmol , or blood lactate during the submaximal test. The males had a higher absolutė VO 2@4mmol than the females for both WR and NR (both p = 0.002), and a higher GE during NR (17.0 ± 1.3% for males TABLE 2 | Responses to skiing with (WR) and without (NR) the rifle during four submaximal levels (Sub 1-4).
Frontiers in Physiology | www.frontiersin.org vs. 16.0 ± 0.5% for females, p < 0.05), but there was no sex difference in GE for WR.

Maximal TT
Responses during the maximal TT for WR and NR are presented in Table 3. Time to complete the TT was 14.3 ± 9.4 s longer for WR compared to NR (p < 0.001; n = 17). No differences were found inVO 2max between WR and NR, whereas MR an (both absolute and relative) was higher for NR (Table 3). No differences in maximal RPE (for breathing, arms, or legs) were observed between WR and NR. Compared to the female athletes, the males had a higherVO 2max (absolute and relative),V E and absolute MR ae during the TT for both WR and NR, and a higher absolute MR an for WR ( Table 3). The differences in mean speed and mean power output between the four quarters and the two skiing conditions (WR/NR) during the TT for the group in total are illustrated in Figure 2. There was no difference in MR ae between WR and NR during the TT, while Data are presented as mean ± SD. Different from WR: * p < 0.05, * * p < 0.01, * * * p < 0.001. Different from females for the same condition (WR/NR): ‡ ‡ p < 0.001, ‡ ‡ ‡ p < 0.001.VO 2 , oxygen uptake;V E , ventilation rate; HR, heart rate; MR an , Anaerobic metabolic rate; MR ae , Aerobic metabolic rate; ES, effect size.
the absolute MR an was higher during the first half of the TT for NR (Figure 3).

DISCUSSION
This study aimed to investigate the physiological factors affected by rifle carriage during biathlon skiing performance, as well as the sex differences associated with rifle carriage. The main findings showed that: (1) rifle carriage affects respiratory responses, GE and speed @4mmol during submaximal treadmill roller-skiing, as well as the anaerobic energy contribution and performance during a maximal roller-skiing TT; (2) speed @4mmol and anaerobic energy contribution seem to explain the majority of the variation in maximal TT performance lasting ∼ 3-4 min; and (3) although the mass of the rifle was greater relative to body mass for the female athletes, there were no sex differences in the physiological responses to rifle carriage. Previous studies investigating the physiological response to rifle carriage have reported average treadmill roller-skiing speeds of 2.22-2.97 m.s −1 (Rundell and Szmedra, 1998;Stöggl et al., 2015), while average speeds in the present study were 1.94-4.19 m.s −1 for females and 2.22-4.43 m.s −1 for males during the submaximal and maximal tests. The difference in top average speeds reported in the present study compared to the study by Stöggl et al. (2015), both of which involved skiing at competition speeds, can be at least partly explained by the steeper treadmill incline used in the study by Stöggl et al. (2015), which was 5 • compared to 3.5-4.5 • in the present study. While the average TT speeds in the present study were more closely matched to real-world competition speeds than those reported in previous studies, they are still much lower than the mean speeds of 6.3 and 7.2 m.s −1 for females and males, respectively, reported by Laaksonen et al. (2018) for World Cup competitions. However, it is important to recognize that these speeds were derived from the top ten athletes in the World Cup, and that World Cup races include all types of terrain (i.e., flat, uphill, and downhill).
During submaximal exerciseVO 2 ,VCO 2 ,V E , and E GROSS were higher for WR compared to NR, which is consistent with findings reported by Stöggl et al. (2015) and Rundell and Szmedra (1998). In addition, blood lactate concentration was higher for WR compared to NR following the fourth submaximal level for the group in total as well as for the males, suggesting that the aerobic energy system is able to meet the energy demands at lower intensities for WR, while this is not possible at higher workloads. This is consistent with previous studies, where differences in blood lactate concentrations were demonstrated for workloads around or  above the lactate threshold (Rundell and Szmedra, 1998;Stöggl et al., 2015). However, it is relevant to note that the interindividual differences in blood lactate responses when skiing WR and NR were large (see Figure 1) and it may therefore be beneficial to individualize training while carrying the rifle depending on the blood lactate responses of each athlete.
Biathletes and XC skiers often use HR monitoring to control exercise intensities when training. In the present study, no differences in HR were observed between WR and NR, although the respiratory responses (VO 2 ,VCO 2 ,V E ) were higher during WR. Day-to-day variations in HR of up to 4.1% have been reported during submaximal exercise (Taylor, 1944). Therefore, both HR and blood lactate concentrations are likely not sensitive enough to detect the differences in respiratory responses between WR and NR. This is important to consider when prescribing and analyzing training loads, since the higher metabolic demands of skiing WR will probably not be reflected by HR recordings or sampling of blood lactate concentrations following a training or racing bout. Thus, more accurate methods of monitoring and quantifying training intensities in biathlon, particularly when skiing WR, are required. Methods for calculating power output are currently being investigated within XC skiing, which could improve training-load quantification. Such possibilities need further scientific investigation within biathlon.
Previous research has shown that around 60% of performance in biathlon sprint competitions can be explained by skiing speed (Luchsinger et al., 2018;Dzhilkibaeva et al., 2019), with the remaining 40% explained by the shooting time and accuracy. Speed at the lactate threshold,VO 2max , GE, and anaerobic energy contribution are suggested to be the most important factors for performance in endurance-based sports (Joyner and Coyle, 2008), including XC and biathlon skiing (Rundell, 1995;Rundell and Bacharach, 1995;Larsson et al., 2002;Carlsson et al., 2012Carlsson et al., , 2016. Therefore, these variables together with muscle mass were used as independent variables in the multiple regression analyses in the present study. Results showed that the most important factors for skiing performance in a TT lasting ∼ 3-4 min were the speed at 4 mmol·L −1 of blood lactate and the anaerobic energy contribution, which is supported by previous research (Rundell, 1995;Rundell and Bacharach, 1995;Andersson et al., 2017).
The reported training load performed with the rifle was approximately 4 h.week −1 (i.e., ∼ 30% of the endurance-based training) for the participants in the present study, but the individual variation was large (ranging from 1-10 h.week −1 ). The speed @4mmol and the anaerobic energy contributions were lower during WR compared to NR, which may be a result of the lower training load performed while carrying the rifle on the back. Although not investigated in the present study, anecdotal reports indicate that biathletes often perform highintensity training without carrying the rifle, which may therefore affect developments in anaerobic energy contribution. These findings suggest that training to improve speed @4mmol and anaerobic energy contribution while skiing with the rifle may help to improve skiing performance while carrying the rifle. While speculative, it may be hypothesized that increasing the relative training load with the rifle may be required by biathletes in order to improve biathlon skiing performance.
The duration of the TT in the present study corresponds to the ski-lap durations in a biathlon single mixed-relay (∼ 3-3.5 min per lap). Although the skiing laps are longer in other biathlon disciplines, the importance of the speed @4mmol and anaerobic energy contribution for performance may still play an important role due to the variation in terrain and the requirement to increase speed rapidly (e.g., sprinting to the finish line). Given the intermittent nature of biathlon skiing, the variation in terrain (which changes continuously between uphill, downhill, and flat sections) and the variation in competition duration and formats, more research is needed to understand the specific demands and energy system contributions during biathlon competitions, particularly with respect to the unique requirement for rifle carriage (compared to XC skiing, for example).
In biathlon races, the sex difference in skiing speed is 10-15% and this difference explains more than 90% of the difference in overall performance between male and female athletes (Laaksonen et al., 2018;Luchsinger et al., 2018Luchsinger et al., , 2019. In the present study, despite the greater relative weight of the rifle for the females compared to the males (5.6 vs. 5.0% of body weight), no sex differences in absolute or relative responses to rifle carriage were identified. This is inconsistent with previous research showing that female biathletes had greater differences than males inV E (Stöggl et al., 2015) andVO 2 (Rundell and Szmedra, 1998) between WR and NR. This difference in findings may be explained by the greater difference in body mass between the male and female biathletes in the study by Rundell and Szmedra (1998) compared to the present study (∼ 16 kg vs. ∼ 7 kg).

CONCLUSION
As hypothesized, rifle carriage increased the respiratory responses during submaximal roller-skiing, without any differences during maximal exercise. However, in contrast to our hypothesis, no sex differences were observed. In general, rifle carriage led to a decrease in the skiing speed at a blood lactate concentration of 4 mmol.L −1 , as well as decreases in the anaerobic energy contribution and performance during a maximal TT lasting ∼ 3-4 min. The most important variables for predicting skiing speed in biathlon identified in the present study were the speed at 4 mmol.L −1 of blood lactate concentration and the anaerobic metabolic rate, both of which were lower when skiing with the rifle compared to without. Thus, to improve skiing performance in biathlon, improving skiing speed at 4 mmol.L −1 of blood lactate concentration and anaerobic energy delivery while carrying the rifle are recommended.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by the Regional Ethical Review Board in Umeå, Sweden. The participants provided their written informed consent to participate in this study.

AUTHOR CONTRIBUTIONS
MJ initiated the study, collected and analyzed the data, and wrote the first draft of the manuscript. KM and ML contributed to the study design, data analysis, and writing of the manuscript. MJ, KM, and ML approved the final version to be published and agreed to be accountable for all aspects of the work.

FUNDING
This study was part-funded by Östersund City Council.