The P and S wave receiver functions and their joint inversions are used to study the lithosphere beneath the Lake Van region using approximately 600 teleseismic earthquake data (Mw ≥ 5.8) at different azimuths collected from 10 permanent broadband stations operated by Kandilli Observatory and Earthquake Research Institute and Disaster and Emergency Management Authority in the region. The dataset is taken from the European Integrated Data Archive. The simulated annealing method is used for the joint inversion of P-wave and S-wave receiver functions, since the P410s phase cannot be detected reliably in the depth stacks. This may be due to the low olivine content and high basalt content at this depth. The inversion process is therefore performed without travel time residuals. The crustal thickness is observed at nearly 45 km in all the velocity models obtained from the inversion. The relatively low-velocity layer (Vs = 3.4 km/s) at depths of the middle crust may be associated with volcanic centers near Nemrut, in the west of the region. Additionally, another low-velocity layer (Vs = ~ 3.0 km/s) is observed in the upper crust around the Suphan Volcano. Also, Vp/Vs and Poisson’s ratios are calculated for the study area. Their high values (Vp/Vs ≥ 1.85, and σ ≥ 0.285) correspond to the partial melting of the lower crust in the region. According to velocity models obtained from the inversion results, the most important outcome is that the average variance of the Moho discontinuity (~ 45 km) and lithosphere–asthenosphere boundary (~ 90 km) under the seismograph stations can be interpreted as a likely part of the Arabian oceanic plate in relation to the collision zone.