Optical glass scratching can induce various types of cracks, among which median cracks are extremely detrimental and penetrate deeply under the surface. Due to deep-scratching process complexity, it is challenging to devise a method to predict median crack depth. Indentation testing has been examined comprehensively in prior research works. It has been found that using the correlation between scratch and indentation testing can simplify predictive method development. In this research, a numerical method based on indentation testing is proposed to determine median crack depth during deep scratching. In the first step, an FE model is configured to simulate the indentation testing process and the Cohesive Zone Method is applied to describe median crack behavior. The cohesive parameters calibrated through experimental indentation testing are implemented in the FE scratch model, and the results are compared with the experimental scratch test results. According to the results, the FE scratch model was enhanced by mode II fracture energy and the modeled friction coefficient. The indentation and scratch experiments were conducted with BK7, F2, Fused silica, K5, Pyrex, Quartz, SF6, and SF19. The experimental results prove that the nonlinearity of the median crack depth curve correlates with KIc. A comparison of the experimental and numerical results demonstrates the model is virtually functional for materials with KIc below 1000 kPa m1/2. Comparisons between the current findings and other studies infer the model and experimental results are accurate and reliable. © 2018 Elsevier Ltd and Techna Group S.r.l.