## 2 Materials and Methods
### Data Processing  
The marker data were digitized and labelled via QTM [15], and the recorded data were exported to Visual3D biomechanical software (v2023.09.3, C-Motion, Inc., Germantown, MD, USA) for the further computation of kinetic variables. A seven-segment rigid-line body model was built to represent the trunk, pelvis, left- and right-side thigh, shank, and foot segments. Then, this model was used to compute the motion of each segment and joint using a Cardan angle sequence x-y-z [10], where x represents the medial-lateral axis, y represents the anterior-posterior axis, and z represents the longitudinal axis. The raw GRF and marker data were smoothed via a low-pass Butterworth filter with 50 Hz and 12 Hz cut-off frequencies, respectively [10]. The sagittal plane joint power was calculated using Newtonian inverse dynamics, which is the product of the net joint moment and the corresponding joint angular velocity [10; 12]. The GRF and joint power that had been processed in Visual3D were then exported to MATLAB (R2022a; The MathWorks, Inc., Natick, MA, USA) for further calculation.  
The GRF was down-sampled five times to correspond to the motion data at 200 Hz [12] to calculate time-synchronized vertical velocity and displacement that acted on the participants' COM consistent with the previous literature [10; 26]. The braking phase of the CMJ was defined as the time between the local peak negative COM velocity to the COM velocity reaching zero [26]. A time interval between the COM velocity reaching zero and the take-off instant (the net GRF being lower than five times the standard deviation of GRF extracted from only the middle part of the flight phase) was defined as the propulsive phase of the CMJ [26; 27]. The touchdown of the CMJ was defined as the GRF value being larger than five times the standard deviation of GRF extracted from the flight phase [27]. The first touchdown during the DJ was identified as the GRF being larger than five times the standard deviation of vertical GRF extracted from the standing still period (i.e., participants stood on the drop box) [24]. Subsequently, the time between the local peak negative COM velocity and the local lowest COM vertical position was defined as the braking phase of the DJ [24]. The propulsive phase of the DJ was determined as the time between the local lowest COM vertical position and take-off instant (using the same force threshold as the CMJ) [24]. The method to define the braking and propulsive phases in the CMRJ1 was identical to CMJ, while the method to define these phases in CMRJ2 was identical to DJ [7]. Furthermore, the jump height of three jumps in the present study was calculated via the impulse-momentum methods from the COM height at the take-off instant to the height point of the flight phase [5]; the strategy-based metrics were calculated in accordance with methods described by [PERSON] et al. [1].  
The mean negative and positive net joint works were calculated by integrating (through trapezoid rule) the negative and positive portions of the power-time curve during braking and propulsive phases of each jump action [10; 16; 28]. The negative work reflects energy absorption and positive work reflects energy generation [16]. The negative and positive hip, knee, and ankle joint work values were summed first to obtain the total negative and positive work, respectively [10]. Finally, the negative and positive work of each joint was divided by the total negative and positive work values to calculate the individual joint contribution to the total negative and positive lower extremity work [11]. All kinetic variables (joint work and contributions) were normalized to participants' body mass to reduce its influence on statistical comparison.