This study is aimed at collecting and processing digital waveforms and bulletin data from seismic networks in the Middle East and performing a joint inversion of five data types including absolute body wave travel times, differential body wave travel times, surface wave dispersion between station and event pairs, surface wave dispersion between station pairs from ambient noise cross-correlation, andreceiver functions for better constraint on the P and S velocity models for the Middle East. We have collected and processed waveform and bulletin data and constructed a phase picks database for the region. We have also processed ambient noise data and performed waveform cross-correlation to extract empirical Greens functions and measured Rayleigh wave dispersions. Surface wave dispersions between earthquakes and stations are also measured. Teleseismic waveforms are processed for receiver function analysis in the region. A joint tomographic inversion of the processed datasets is performed for high-resolution P and S velocity models for the Middle East. The obtained velocity models show complex features that are in general agreement with the regional tectonics. Resolution analysis indicates that data coverage is sufficient, and tomographic resolution is adequate for the study area. The event relocations from the inversion are also utilized to construct a GT5 database for the region.