<p>Convergence-related shortening gets primarily accommodated in faults, fault-related folds and penetrative strain in fold thrust belts (FTB). For example, in the Himalayan FTB, ~477-919 km minimum orogen-scale shortening is accommodated by a series of folded, south vergent thrust systems that vary laterally in their geometry resulting in laterally varying shortening distribution. From hinterland to foreland, these major thrust faults are the Main Central thrust, the Pelling-Munsiari thrust, the Lesser Himalayan duplex, the Main Boundary thrust, and the Main Frontal thrust. In the Sikkim Himalayan FTB, the structural geometry of these thrust sheets laterally varies over ~15 km. Based on two regional, transport-parallel balanced cross-sections, ~542-589 km minimum wedge-scale shortening has been estimated. To quantify grain-scale shortening, we analyzed 201 thin-sections cut from 96 quartz-rich samples (sandstone, quartzite, phyllite, schist, and gneiss) and calculated penetrative strain from them. Penetrative strain results indicate that ~25-26% of total Himalayan shortening is recorded at the grain-scale in this section of the eastern Himalaya.</p><p>In the internal thrust sheets, the strain magnitude (R<sub>S</sub>) remains higher (~1.4-2.43 ), and it progressively decreases in the frontal thrust sheets (~1.08-1.51). The normalized Fry and the R<sub>f</sub>-&#966; are the two most commonly used graphical methods to estimate best-fit strain ellipse parameters, i.e., R<sub>S</sub> and &#966; (long-axis orientation). However, in thrust sheets with less deformed sandstones, where initial grain shapes were not spherical, these graphical methods do not accurately estimate the best-fit strain ellipse parameters. The central vacancy in the Fry plot was objectively fitted using the enhanced normalized Fry (ENFRY), the point-count density (PCD), the continuous function method (CFM), and weighted least square (WLS) methods. From the R<sub>f</sub>-&#966; data, we calculated the best-fit strain ellipse using the shape matrix eigenvector (SME), centroids of the hyperbolic plot (HP), Elliot&#8217;s polar graph (EPG), and R<sub>f</sub>-&#966; graph, harmonic mean (HM) and vector mean (VM) methods. In this study, we calculate the accuracy of these strain methods as a function of the strain magnitude and structural position within the orogenic wedge. The SME and HP methods record the lowest bootstrap errors in the strain parameters in the internal thrust sheets. In contrast, R<sub>S</sub> and &#966; values estimated by the WLS method records the lowest bootstrap error in the frontal thrust sheets, followed by the SME, HP, and EPG methods. We also created six synthetic aggregates containing 150-170 random elliptical grains with random long-axis orientations. We deformed these aggregates under pure-shear, simple-shear, and general-shear conditions at various strain increments. We have generated 7560 strain data. To understand the accuracy of these strain methods in estimating penetrative strain, we calculated the Root Mean Square Error (RMSE) and Mean Absolute Error (MAE) for every strain method and every type of deformation. Experimental results indicate that the SME and HP methods record the lowest errors in the R<sub>S</sub> and &#966; values. In low strain conditions (R<sub>S</sub><1.5), the SME, HP, and EPG methods record lower errors in the strain parameters. Therefore, this study shows that the SME and HP methods overall yield a better penetrative strain estimate.</p>