Multifractal Analysis of Movement Behavior in Association Football

Research in football has been embracing the complex systems paradigm in order to identify different insights about key determinants of performance. The present study explored the multifractal properties of several football-related scenarios, as a candidate method to describe movement dynamics. The sample consisted of five footballers that were engaged in six different training situations (jogging, high intensity interval protocol, running circuit, 5 vs. 5, 8 vs. 8 and a 10 vs. 10 small-sided game). All kinematic measures were collected using a 100 Hz wireless and wearable inertial measurement unit (WIMUPRO©). Data were processed using a discrete wavelet leader transform in order to obtain a spectrum of singularities that could best describe the movement dynamics. The Holder exponent for each of all six conditions revealed mean values h < 0.5 indicating presence of long memory with anti-correlated behavior. A strong trend was found between the width of the multifractal spectrum and the type of task performed, with jogging showing the weakest multifractality ∆h = 0.215 ± 0.020, whereas, 10 vs. 10 small-sided game revealed the strongest ∆h = 0.992 ± 0.104. The Hausdorff dimension indicates that a maximal fluctuation rate occurs with a higher probability than that of the minimal fluctuation rate for all tasks, with the exception of the high intensity interval protocol. Moreover, the spectrum asymmetry values of jogging, running circuit, 5 vs. 5, 8 vs. 8 and 10 vs. 10 small-sided games reveal their multifractal structures are more sensitive to the local fluctuations with small magnitudes. The multifractal analysis has shown a potential to systematically elucidate the dynamics and variability structure over time for the training situations.


Introduction
Wearable technologies have widened the scope and the research questions about the key determinants of sport performance [1,2]. For example, technological development has provided the advent of new tracking devices and training-monitoring systems that allow sports scientists to explore interactions between tactical dimensions and physiological aspects of elite performance in team sports [3,4]. Furthermore, high-frequency inertial measurement units (IMUs) are being used to capture the players' locomotion and, thus, providing a valid measure to assist in the understanding of movement dynamics [5,6]. As an example, several studies have described the use of accelerometers to detect physical activities, movement patterns and to identify activity profiles in team sports players [5,7,8].
The way that each player moves on the pitch obeys the team's collective dynamics, thus depicting great variability of behavior with different magnitudes and temporal dynamics [9]. According to Sampaio and Maçãs [10], such behavior features highly nonlinear, deterministic and stochastic characteristics, furthermore, the measurement of complex biological systems such as human locomotion via IMUs, often generates extremely irregular signals that tend to exhibit intricate characteristics The global scale exponents α are calculated as the slope of the linear regression between the log (µ s ) versus the log (s). It follows that, if the slope of the log-log plot between the scale of observation and the descriptive statistics is equal 0 < α < 0.5, then the process has a memory, and it exhibits anti-correlations. Therefore, the successive random elements in a sequence present an increasing trend in the past that is likely to be followed by a decreasing trend in the future, and vice-versa. Conversely, if the slope is equal to 0.5 < α < 1, then the process has a memory, and it exhibits positive correlations. Therefore, the successive random elements in a sequence present an increasing trend in the past that is likely to be followed by an increasing trend in the future. If the slope is α = 0.5, then the process is indistinguishable from a random process with no memory (random walk). Finally, if the slope is 1 < α < 2, then the process is non-stationary. Most real-world phenomena yield signals that mix different patterns of correlations over time and in different time scales, most often presenting irregular behavior with interwoven periods of high and low correlation. In such cases the linear regression of log 2 µ(s) on log 2 s does not exhibit a straight-line characterizing mono-scaling behavior by a single Symmetry 2020, 12, 1287 3 of 12 scale exponent, conversely, additional exponents are necessary. In fact, signals with interwoven fractal subsets become more complex, hence requiring multiple scaling exponents for the full description of different parts of the signal. Consequently, in contrast to the monofractal formalism, where a global scale exponent α is defined for the entire signal, multifractal time series are defined by a set of local singularity exponents. According to Ihlen and Vereijken [19], signals present in biological systems exhibit multiscaling behavior, for example, signals of inter-stride intervals in human walking or inter-beat interval in electrocardiogram data. Moreover, almost all data points in the aforementioned signals can be represented as a singularity at time instant t 0 with a strength that is numerically defined by a singularity exponent α. These apparent random variations can be misinterpreted as pure noise originating from the system, however, studying the correlation structure of noise-like signals can elucidate important information about the process underpinning its generation.
Summing up all the previous considerations, there is a clear opportunity to explore how the multifractal analysis can reveal new insights about sports performance and discuss how its outcomes can extend current knowledge in the particular topic of movement analysis in football. Therefore, the aim of this study is to identify the multifractal spectrum of different football-related activities (jogging, high intensity interval protocol, circuit running, 5 vs. 5, 8 vs. 8 and 10 vs. 10 small-sided games).

Participants
Five young football players (16.6 ± 0.5 years, 1.76 ± 0.03 m height, and 68.4 ± 4.2 kg body mass) participated in this study. Criteria for inclusion were applied to ensure all players were engaged in amateur competitions and organized training schedules for, at least, 3 times a week. They also would have to be injury-free for at least 3 months. Participants, their legal guardians and their coaching staff, were fully informed about the purpose, benefits and risks of the study, and were provided with written informed consent before the study started. The study protocol conformed to the recommendations of the Declaration of Helsinki and was approved by the local institution (UID/DTP/04045/2019).

Data Collection
Inertial data of all participants were collected using individual WIMU units (RealTrack Systems, Almería, Spain). The validity and reliability of WIMU© system have been already reported previously and their operation procedures and handling are well documented elsewhere [20]. In order to decrease the measurement errors, the players used the same unit across all the testing situations. Accelerometer data were sampled at 100 Hz. After each session, data were downloaded using proprietary software (WIMU PRO, Real Track System, Almeria, Spain) and raw data were carried for further analysis with MATLAB ® software (MathWorks, Inc., Natick, MA, USA).
To get a comprehensive understanding of the fluctuation characteristics of player movement behavior, we investigated the multifractality of six accelerometry signals belonging to the following experimental conditions: jogging, high intensity interval protocol (HIIP), circuit running (CR), 5 vs. 5, 8 vs. 8 and 10 vs. 10 small-sided games (SSG), as illustrated in Figure 1. The conditions are described as follows:

•
Jogging: low intensity self-paced steady state running. • HIIP: consisting of 10 s of maximum effort of linear running interspersed with 20 s active recovery. • RC: comprised of walking, jogging, running and sprinting. To each element was assigned a number, 1 to 4, respectively, and a sequence of numbers was randomly generated using an Excel spreadsheet. The exercise intensity was changed every 5 s according to the sequence using a verbal command. spreadsheet. The exercise intensity was changed every 5 s according to the sequence using a verbal command.  The small-sided games were played in a pitch area proportional to the number of players on small (43 m × 64 m), medium (70 m × 64 m) and large (103 m × 64 m) fields, respectively. The activities of jogging, HIIP and RC were 5 min long. The small-sided games were 10 min long and the initial and final 2 min 30 s were rejected for analysis, avoiding adaptation to the task and fatigue issues, thus leaving 5 min for analysis. Prior to each data collection, participants received proper instruction regarding the execution of the task and subsequently performed a standardized warmup consisting of 5 min of calisthenics and active stretching, followed by 5 min of jogging.

Data Processing
All signals were analyzed in raw format and no pre-processing technique was employed. To investigate the multifractal features of the data, the Discrete Wavelet Leader analysis of the accelerometry signal was computed using MATLAB software (MathWorks© 2019b). A biorthogonal wavelet "bior1.5" was used as the mother wavelet function. The maximum decomposition level was set to log2 of the sample size (signal length). The regression type was set to "uniform", which means that the leader coefficients were uniformly weighted in all decomposition levels. The range of statistical moments (q-moment) used varied from −5 to +5. The small-sided games were played in a pitch area proportional to the number of players on small (43 m × 64 m), medium (70 m × 64 m) and large (103 m × 64 m) fields, respectively. The activities of jogging, HIIP and RC were 5 min long. The small-sided games were 10 min long and the initial and final 2 min 30 s were rejected for analysis, avoiding adaptation to the task and fatigue issues, thus leaving 5 min for analysis. Prior to each data collection, participants received proper instruction regarding the execution of the task and subsequently performed a standardized warmup consisting of 5 min of calisthenics and active stretching, followed by 5 min of jogging.

Data Processing
All signals were analyzed in raw format and no pre-processing technique was employed. To investigate the multifractal features of the data, the Discrete Wavelet Leader analysis of the accelerometry signal was computed using MATLAB software (MathWorks© 2019b). A biorthogonal wavelet "bior1.5" was used as the mother wavelet function. The maximum decomposition level was set to log 2 of the sample size (signal length). The regression type was set to "uniform", which means that the leader coefficients were uniformly weighted in all decomposition levels. The range of statistical moments (q-moment) used varied from −5 to +5.

Multifractal Discrete Wavelet Leader
The discrete wavelet leader (DWL) is an analytical technique rooted in signal processing that provides a unity of approach to the multifractal analysis complex signals [21]. The DWL deals well with nonstationary, heterogeneous and transient behavior, therefore, it can be useful for analyzing a range of complex signals. The DWL is a nonparametric method and requires no previous knowledge about the system that will be under study.
Essentially, the DWL aims to characterize the fluctuations of the time of the local singularity in a signal x[t] or function f (t) by quantifying the size of the sets of points in R d , sharing a same local singularity h. Local singularity can be quantified by pointwise exponents, the most common one being the Hölder exponent, h f (t 0 ), defined as follows: f belongs to C α (t), α ≥ 0, if there exists a polynomial P t with deg(P t ) < α and a constant C > 0 such that: f (t + a) − Pt(t + a) ≤ C|a| α when |a| → 0 . The Hölder exponent consists of the largest such α : Essentially, the larger h(t 0 ), the smoother f around t 0 , and conversely, the closer h(t 0 ) to 0, the more irregular f at t 0 . The wavelet leaders have shown to be, theoretically, the only multiscale measure able to characterize the Hölder exponent [22]. To describe the characteristics of a local singularity, DWL analysis relies upon the discrete wavelet transform (DWT) of a signal X = {x k , K ∈ Z}, that is given by: where ψ 0 is the "mother" wavelet with a compact time support, t is time script, j is the dyadic scale parameter, and k relates to the shift parameter. According to Ihlen and Vereijken [19], the DWT calculates the scale-dependent measure d x ( j, k) by a recursive high and low pass filtering of the signal X. The authors describe the process as follows: "the scale-dependent measure d j=1,k is the convolution product ⊗ of the signal x[t] and the high pass filter ψ hi t subsampled by 2. Then, the subsampled convolution product ⊗ of the time series x[t] and the low pass filter ψ lo t is used in another convolution and subsample procedure to define the scale-dependent measure d j=2,k . The subsampling by 2 leads to dyadic intervals λ j,k that halves the sample sizes j = 1/2 j for each new level j, but double the frequency band f = 2 j . The wavelet leaders are therefore defined as the local maxima of 3 subsequent intervals λ j ,k on scale j and for all j < j , as follows" [19]: The size and variability of the sets of points where the pointwise singularity exponent takes a given value h is mathematically formalized as the Hausdorff dimension D [23], which encapsulates the Hausdorff dimensions of the sets of points where a pointwise singularity exponent takes a given value h. The Legendre formalism describes the multifractal spectrum D(h) indirectly by the q-order statistical moment S( j, q) of the local scale-dependent measure L X ( j, k). Then, for a signal X, let S(q, j) denote the structure function and ζ(q) the scaling exponents, where q is the order (or moment) of multi-resolution. They are expressed as follows: where L X represents the largest wavelet coefficient calculated at all finer scales, n j is the number of leaders L X ( j, k) available at scale 2 j defined on each scale j. The q-order mass exponent ζ q is estimated by the linear regression of Equation (4) in log coordinates vs. log( j) for several q-orders.
The q-order scaling exponent ζ q is then Legendre-transformed into the corresponding dimension D(h) by the following equation: According to Ihlen and Vereijken [19], "the Legendre formalism assumes that ζ q is a convex function of q such that its tangent slope, the singularity exponent h q , is a monotonically decreasing function of q. The positive q-orders will amplify L X ( j, k) in Equation (4) for the large fluctuations in the signal whereas the negative q-orders will amplify L X ( j, k) for the small fluctuations" (seen in the red and black intervals in Figure 2, respectively). The q-order's singularity exponents h q for negative q represents the right tail of the singularity spectrum whereas positive q represents the left tail of the singularity spectrum. In the special case of the monofractal noise-like time series, ζ q is a linear function of q and, consequently, h q and D(h) are q-independent constants collapsing the singularity spectrum D h q into a very narrow support.
Symmetry 2020, 12, x FOR PEER REVIEW 6 of 14 According to Ihlen and Vereijken [19], "the Legendre formalism assumes that is a convex function of such that its tangent slope, the singularity exponent ℎ , is a monotonically decreasing function of . The positive -orders will amplify ( , ) in Equation (4) for the large fluctuations in the signal whereas the negative -orders will amplify ( , ) for the small fluctuations" (seen in the red and black intervals in Figure 2, respectively). The -order's singularity exponents ℎ for negative represents the right tail of the singularity spectrum whereas positive represents the left tail of the singularity spectrum. In the special case of the monofractal noise-like time series, is a linear function of and, consequently, ℎ and (ℎ) are -independent constants collapsing the singularity spectrum into a very narrow support.

Interpreting the Multifractal Singularity Spectrum
This section is devoted to the explanation and interpretation of the components of the multifractal spectrum. A narrow singularity spectrum represents relative systematic fluctuations in the signal that occur over similar amplitudes over time, resembling a mono-scaling structure. A broad singularity spectrum represents more intermittent and asystematic fluctuations in the movement that occur over a greater range of amplitudes in different time lengths, a benchmark of multifractal processes. Therefore, the range of the singularity spectrum is suggested to be a measure that captures the dynamics of the system [24,25].
The singularity exponents ℎ for the large fluctuations are small and located in the left tail of the spectrum, whereas ℎ for the small fluctuations are large and located in the right tail of the spectrum. The strength of multifractality is therefore described by the large deviation of the local singularity exponent ℎ from the central tendency ℎ(0). The monofractal signal is the case when ℎ is almost constant and presents very short support, in some cases, collapsing into a single point at ℎ [19].

Interpreting the Multifractal Singularity Spectrum
This section is devoted to the explanation and interpretation of the components of the multifractal spectrum. A narrow singularity spectrum represents relative systematic fluctuations in the signal that occur over similar amplitudes over time, resembling a mono-scaling structure. A broad singularity spectrum represents more intermittent and asystematic fluctuations in the movement that occur over a greater range of amplitudes in different time lengths, a benchmark of multifractal processes. Therefore, the range of the singularity spectrum is suggested to be a measure that captures the dynamics of the system [24,25].
The singularity exponents h for the large fluctuations are small and located in the left tail of the spectrum, whereas h for the small fluctuations are large and located in the right tail of the spectrum. The strength of multifractality is therefore described by the large deviation of the local singularity exponent h from the central tendency h(0). The monofractal signal is the case when h is almost constant and presents very short support, in some cases, collapsing into a single point at h [19]. The range h indicates the variety of singularity exponents that describe the dynamics of the system while the value of D(h) indicates the probability of the contribution of a respective h.
Following this, the range of a singularity spectrum is subdivided into a left limb (L) and a right limb (R) (Figure 3), related to the dynamics of the highest and lowest fluctuations, respectively. The asymmetry degree of the multifractal spectrum can be identified using ∆S. If ∆S < 0, the multifractal spectrum presents a long left tail, indicating that the structure in the time series is not sensitive to local fluctuations with small magnitudes. If ∆S > 0, the multifractal spectrum presents a long right tail, indicating that the structure in the signal is not sensitive to the local fluctuations with large magnitude [24]. In cases where high and low fluctuation components of a signal present a similar scaling complexity, the singularity spectrum will result in being approximately symmetrical and L = R [25]. In cases of more complexity in either of these components, the singularity spectrum becomes asymmetrical, resulting in one limb being broader than the other. The range ℎ indicates the variety of singularity exponents that describe the dynamics of the system while the value of (ℎ) indicates the probability of the contribution of a respective ℎ.
Following this, the range of a singularity spectrum is subdivided into a left limb ( ) and a right limb ( ) (Error! Reference source not found.3), related to the dynamics of the highest and lowest fluctuations, respectively. The asymmetry degree of the multifractal spectrum can be identified using Δ . If Δ < 0, the multifractal spectrum presents a long left tail, indicating that the structure in the time series is not sensitive to local fluctuations with small magnitudes. If Δ > 0, the multifractal spectrum presents a long right tail, indicating that the structure in the signal is not sensitive to the local fluctuations with large magnitude [24]. In cases where high and low fluctuation components of a signal present a similar scaling complexity, the singularity spectrum will result in being approximately symmetrical and = [25]. In cases of more complexity in either of these components, the singularity spectrum becomes asymmetrical, resulting in one limb being broader than the other.
The nature of the higher fluctuation rate occurring with a higher probability can be interpreted using Δ , calculated by the difference between the fractal dimensions of maximum and minimum probability subsets [24]. Therefore, if Δ > 0, "the maximal fluctuation rate occurs with a higher possibility than that of the minimal fluctuation rate, and vice versa" [24].

Results
The individual multifractal spectra corresponding to all conditions are presented in Figure 4. Average values for each condition can be found in Figure 5, with additional information related to those multifractal spectra presented in Table 1.
The mean generalized Holder exponent ℎ(0) for each of all six conditions reveled values less than 0.5. By definition, this means that all six conditions presented long memory with anti-correlated behavior.
Analysis among the width of the multifractal spectrum Δ ℎ revealed a trend between the condition performed and the width of the multifractal spectrum. Indeed, the 10 vs. 10 SSG had the widest (ℎ = 0.99) and jogging (Δℎ = 0.21) showed the narrowest.
The mean values of Δ for jogging, circuit running, 5 vs. 5, 8 vs. 8 and 10 vs. 10 were all negative, meaning that the maximal fluctuation rate occurred with a higher probability than the The nature of the higher fluctuation rate occurring with a higher probability can be interpreted using ∆D h , calculated by the difference between the fractal dimensions of maximum and minimum probability subsets [24]. Therefore, if ∆D h > 0, "the maximal fluctuation rate occurs with a higher possibility than that of the minimal fluctuation rate, and vice versa" [24].

Results
The individual multifractal spectra corresponding to all conditions are presented in Figure 4. Average values for each condition can be found in Figure 5, with additional information related to those multifractal spectra presented in Table 1.
The mean generalized Holder exponent h(0) for each of all six conditions reveled values less than 0.5. By definition, this means that all six conditions presented long memory with anti-correlated behavior.
Analysis among the width of the multifractal spectrum ∆h revealed a trend between the condition performed and the width of the multifractal spectrum. Indeed, the 10 vs. 10 SSG had the widest (∆h = 0.99) and jogging (∆h = 0.21) showed the narrowest.
The mean values of ∆D h for jogging, circuit running, 5 vs. 5, 8 vs. 8 and 10 vs. 10 were all negative, meaning that the maximal fluctuation rate occurred with a higher probability than the minimal fluctuation rate, however, not for all the players. Contrarily, the mean value of ∆D h for HIIP Symmetry 2020, 12, 1287 8 of 12 was positive, indicating that the minimal fluctuation rate occurred with a higher probability than the maximal fluctuation rate.

Discussion
This study aimed to explore the multifractal features of football players' movement behavior using accelerometry measures by means of the discrete wavelet leader method. Overall, our exploratory results revealed a strong trend between the task and the shape and size of the multifractal

Discussion
This study aimed to explore the multifractal features of football players' movement behavior using accelerometry measures by means of the discrete wavelet leader method. Overall, our exploratory results revealed a strong trend between the task and the shape and size of the multifractal The mean values of ∆S for the jogging, circuit running, 5 vs. 5, 8 vs. 8 and 10 vs. 10 were all positive, indicating that multifractal spectra had long right tails, and multifractal structures were more sensitive to local fluctuations with small magnitudes, if compared to the fluctuations with large magnitudes. Conversely, the mean value of ∆S for the HIIP is negative, indicating that its multifractal spectrum had a long left tail, and its multifractal structure was more sensitive to local fluctuations with large magnitudes.

Discussion
This study aimed to explore the multifractal features of football players' movement behavior using accelerometry measures by means of the discrete wavelet leader method. Overall, our exploratory results revealed a strong trend between the task and the shape and size of the multifractal spectrum. Indeed, more complex tasks led to wider multifractal spectrums while less complex tasks resulted in narrower multifractal spectrums. The increment in task complexity reshaped the multifractal spectrum most often systematically by stretching its right leg to the right assuming long-range positive correlations in lower fluctuations, whereas, the left end of the spectrum remained almost unchanged. This means that the analysis was sensible in capturing small fluctuations but ineffective in capturing larger fluctuations. Since the h min for jogging was already very small at the beginning h min = 0.069, it became harder to detect greater fluctuations in accelerometry. In conjunction, the high-pass wavelet filter appeared to have an inferior resolution than what was required to capture the high fluctuations, however, it performed considerably well on small fluctuations.
Comparing the values of the Hausdorf dimension of pointwise singularities D(h) for each condition demonstrates that the maximal fluctuations were associated with more probability in all tasks if compared to the small fluctuations, with the exception of HIIP. Indeed, the work-to-rest ratio (i.e., work 1:2 rest) used in the HIIP protocol clearly emphasized a longer duration for low fluctuations (rest) during the activity, thus, conferring more probability if compared to higher fluctuations.
Interestingly, the dynamics of SSG showed a similar pattern for the left leg of the spectrum, almost superimposing its components. This may suggest that the periods with large fluctuations behave similarly among these conditions. Contrarily, the right leg of the spectrum stretched to the right with an increase in the number of players and field dimension, which seems reasonable since players may experience more periods with lower intensity (walking or standing). The summation of this behavior turned out to make the signal uncorrelated, which corroborates with the increase of h(0) values from anti-correlated to no correlation.
The multifractal spectrum (width and shape) could reflect the adaptability of movement to external variation and shown the characteristic configuration for each activity. Adopting a similar approach Fuss, Subic [23] explored the Hausdorff dimension along with the average amplitude of the acceleration signal to identify wheelchair rugby accelerations among different activities. The results showed that Hausdorff dimension values of high speed pushing and high speed coasting were similar, and also show the average amplitude values of collisions and high speed pushing. In addition to this first exploration, it seems clear that the multifractal analysis of the movement behavior of the football players introduces a complementary dimension to the research of performance and coordination of human behavior in dynamic environments. In fact, the close relation between stochastic processes and movement behavior makes multifractal analyses an innovative and promising candidate tool for investigation.
Future studies might consider different parameters for the wavelet leader analysis, such as the wavelet type and the number of decomposition levels, as well as pre-processing techniques, such as filtering. It is important to highlight that the wavelet leader multifractal analysis is very sensitive to any of the above-mentioned conditions.
The jogging activity showed the narrowest multifractal spectrum, resembling a monofractal structure, with a h value centered at 0.15, which indicates an overall long memory anticorrelated behavior. Indeed, jogging is a stationary activity with regular fluctuations that might change due to inherent attributes of motor control (variability of movement) if compared to changes imposed by the task. The elevated value of D(h) in the left end of the spectrum indicates that the high fluctuations present more probability throughout the signal if compared to low fluctuations.
The HIIP activity showed a multifractal spectrum with an h value centered at 0.23, which indicates overall long memory anticorrelated behavior and a larger range of singularity exponents than the jogging activity. In HIIP, the players had to perform all-out accelerations causing tremendous rise in the amplitude of the accelerometry signal while maintaining jogging during resting. Consequently, the spectrum of singularities indicates higher fluctuations if compared to the jogging activity with a lower probability if compared to low fluctuations.
With the exception of player 5 who showed symmetry and almost equal probability of large and small fluctuations, all spectrums show a long left leg (R < L)-an aspect that can be attributed to the change in high intensity running and maintenance of small fluctuations D h (h max ) > D h (h min ). The fact that the low fluctuations were shifted a little to the right might indicate that adaptations in the rest period, to balance the effects of fatigue, made the fluctuations while jogging less anti-correlated, or less noisy.
The circuit activity showed a larger multifractal spectrum than jogging and HIIP activities. That was the first case an activity spectrum started to show positive correlations. This might have happened due to the inclusion of walking, which created a positive correlation in the spectrum. Clearly walking introduces a new pattern of behavior with very low fluctuation in its signal dynamics. Except for player 3, all spectrums showed a long right leg (R > L) and suppression of small fluctuations over large fluctuations D h (h max ) < D h (h min ), suggesting a higher probability of jogging, running and sprinting than walking.
The 5 vs. 5 SSG showed a multifractal spectrum with an h value centered at 0.40 which indicates an overall anticorrelated behavior very close to no correlation. Its range of singularity spectrums indicates richer fractality if compared to jogging, HIIP and circuit running, but it performed the poorest in the SSGs. The individual analysis demonstrates different behavior among players with some characterized by long memory with anti-correlation h < 0.5 (i.e., players 2 and 4) and others characterized with no correlation h = 0.5 (i.e., players 1, 3 and 4). With the exception of player 3, all spectrums show a long right leg (R > L) and suppression of small fluctuations over large fluctuations D_h (h_max ) < D_h (h_min).
The 8 vs. 8 SSG showed a multifractal spectrum with an h value centered at 0.41 which also indicates an overall anticorrelated behavior very close to no correlation. The individual analysis demonstrates different behaviors among players with some characterized by long memory with positive-correlation h < 0.5 (i.e., player 1), anti-correlation h < 0.5 (i.e., players 2 and 4) and no correlation h = 0.5 (i.e., players 3 and 5). All spectrums presented wide support. Players 2, 4 and 5 demonstrated long right leg (R > L) with suppression of small fluctuations over large fluctuations D h (h max ) < D h (h min ). Players 1 and 3 showed long left leg (R < L), and, while player 1 showed suppression of large fluctuations over small fluctuations D h (h max ) > D h (h min ), player 3 showed the opposite results.
The 10 vs. 10 SSG showed a multifractal spectrum with an h value centered at 0.44 which also indicates an overall anti-correlated behavior very close to no correlation. The singularity spectrums are characterized by the largest among all conditions. All spectrums show a long right leg (R > L) and suppression of large fluctuations over small fluctuations D h (h max ) < D h (h min ).

Conclusions
The wavelet leader multifractal analysis has shown to provide different insights into the analysis of complex signals in sports performance scenarios. However, more research is needed to explore improvements in the analysis and ways of implementation of this information in the training setting. Although the systematicity presented in the current study leveraged new and fruitful information for the analysis of movement behavior in football, the transfer to training is still a shortcoming. The use of asymmetry parameters, such as ∆D h , ∆h and ∆S, may be a good starting point to explore the demands of training sessions, allowing comparisons among players and providing substance to support different research questions. The current results are exploratory and might be limited by the sample expertise level, however, the experiment has the potential to be used in understanding several key performance questions, namely: (i) collective synergies, (ii) talent development and (iii) pre-to post-injury changes.