In Chile, sand tailings dams represent the most common deposits of mining residues. These structures present a potential risk in terms ofmechanical instability due to their potential susceptibility to seismic liquefaction. In order to manage these risks, it is necessary to take aprobabilistic approach, thus accounting for inherent variability of material properties. However, in practice, implementing such an approach isimpeded by the difficulty of acquiring and managing the data to be used in the reliability calculations and is conditioned by the relevance of theprobabilistic models chosen to represent this variability. This paper proposes a method for onsite determination of the tailings relative density(DR%), and its variability, using dynamic penetration tests. This method was applied to typical Chilean sand tailings dams, and proposes a singlemodel for all such tailings dams by associating a probability model to the variation of DR%. Finally, the validity of this approach is demonstrated byperforming a reliability calculation of liquefaction potential (which is the main cause for the failure of this type of structure in this country) for aparticular sands tailing dam.En Chile los tranques de arenas de relaves representan los depósitos de residuos mineros más comunes. Estas estructuras presentan un potencialriesgo en términos de inestabilidad mecánica debido al fenómeno de la licuefacción sísmicamente inducida. Para conocer este riesgo, es necesarioconsiderar una aproximación probabilística. Sin embargo, en la práctica, la implementación de estas aproximaciones es muy limitada, debido ladificultad obtener los datos necesarios para los cálculos de fiabilidad, y a la selección de los modelos probabilísticos necesarios para representar suvariabilidad. Este artículo propone un método para estimar in-situ la densidad relativa (DR%) en arena de relaves y su variabilidad, a partir deensayos de resistencia a la penetración dinámica. El método se aplica en tranques de arena de relaves chilenos, proponiendo un único modelo parala función de probabilidad del parámetro DR%. Finalmente, el interés de esta aproximación se demuestra mediante un cálculo de fiabilidad delpotencial riesgo de licuefacción, correspondiente a la principal causa de falla en países sísmicamente activos como Chile.