<p class="MsoNormal" style="text-align: justify; text-justify: inter-ideograph; line-height: 120%; layout-grid-mode: char; mso-layout-grid-align: none; margin: 12.0pt 0cm 6.0pt 0cm;"><span lang="EN-US" style="font-size: 10.0pt; mso-bidi-font-size: 11.0pt; line-height: 120%; font-family: 'Times New Roman',serif;">Variability characterization is a key component in the reliability assessment of geotechnical systems, particularly in scenarios involving tailings deposits, where spatial heterogeneity can critically impact design safety. Despite advances in modeling spatial randomness, many reported cases still rely on theoretical assumptions to define appropriate statistical characterizations&mdash;such as theoretical probability density functions (PDFs) and correlation structures&mdash;which may misrepresent site-specific conditions. To bridge this gap, piezocone penetration tests (CPTu) stand out as promising tools for providing continuous measurements along a vertical profile that can be used to define statistical behavior and avoid bias. This paper presents a spatial variability characterization of a bauxite tailings deposit based on mechanical parameters derived from CPTu data. The study includes basic statistical analysis&mdash;mean (<em>&mu;</em>), standard deviation (</span><span lang="EN-US" style="font-size: 10.0pt; mso-bidi-font-size: 11.0pt; line-height: 120%; font-family: Symbol;">s</span><span lang="EN-US" style="font-size: 10.0pt; mso-bidi-font-size: 11.0pt; line-height: 120%; font-family: 'Times New Roman',serif;">), coefficient of variation (CV)&mdash;alongside a comparison with theoretical PDFs. Subsequently, spatial correlation is evaluated through covariance analysis and estimation of the vertical scale of fluctuation (</span><em><span lang="EN-US" style="font-size: 10.0pt; line-height: 120%; font-family: 'Times New Roman',serif;">&delta;</span></em><span lang="EN-US" style="font-size: 10.0pt; mso-bidi-font-size: 11.0pt; line-height: 120%; font-family: 'Times New Roman',serif;">), using a dedicated subroutine that fits theoretical autocorrelation models (TAMs). The deposit is classified as highly variable according to the I<sub>cRW</sub> index. The normal and Weibull PDFs best represent the data distributions. The vertical scales of fluctuation vary significantly: 0.01 m to 4.43 m for cone resistance (<em>q</em><sub>c</sub>), 0.01 m to 4.36 m for sleeve friction (<em>f</em><sub>s</sub>), and 0.01 m to 5.00 m for pore water pressure (<em>u</em><sub>2</sub>). These findings offer valuable input for probabilistic stability and serviceability analyses, contributing to safer and more informed geotechnical designs involving mine tailings. </span></p>