Non-linear acoustic technique is an attractive approach in evaluating early fatigue as well as cracks in material. However, its accuracy is greatly restricted by external non-linearities of ultra-sonic measurement systems. In this work, an acoustical data-driven deviation detection method, called the consensus self-organizing models (COSMO) based on statistical probability models, was introduced to study the evolution of localized crack growth. By using pitch-catch technique, frequency spectra of acoustic echoes collected from different locations of a specimen were compared, resulting in a Hellinger distance matrix to construct statistical parameters such as z-score, p-value and T-value. It is shown that statistical significance p-value of COSMO method has a strong relationship with the crack growth. Particularly, T-values, logarithm transformed p-value, increases proportionally with the growth of cracks, which thus can be applied to locate the position of cracks and monitor the deterioration of materials.