In this paper, mixed mode weight functions (MMWFs) for elliptical subsurface cracks under shear loadings are derived. Reference mixed mode stress intensity factors (MMSIFs), calculated by finite element analysis (FEA), under uniform shear loadings are used to derive MMWFs. The cracks have aspect ratios and crack depth to crack length ratios in the ranges of 0.2–1.0 and 0.05 to infinity, respectively. MMWFs are verified to calculate MMSIFs for any point of the crack front under linear and nonlinear two dimensional (2D) loadings. The derived MMWFs could be effectively used for evaluation of fatigue crack growth under complicated 2D stress distributions.