The objective of this study is to characterize the recent stress distribution on the Main Marmara Fault (MMF) through quasi-static modeling of inter-seismic strain accumulation by considering past earthquakes and het- erogeneous interseismic coupling. Since the MMF is prone to creating large earthquakes, it is crucial to infer the state of stress on the fault. The obtained pre-stress distribution will serve as a basis for simulating dynamic earthquake ruptures via derived stress changes and calculating a realistic map of the peak ground velocity. Due to lack of observed data for future events, International Monitoring System seismic stations may be used to verify future ground motion simulations for destructive earthquakes. To calculate accumulated stress distri- bution on the locked part of the MMF, interpolated GPS velocities and slip-rates along the unlocked portions of the fault interface are used as boundary conditions. Previous studies of interseismic coupling and seismic- ity studies including the repeating earthquakes are also considered. A three-dimensional-FEM mesh for the region is built using a mesh size of ~500 m. Consequently, heterogeneous stress distributions are obtained via elaborative usage of recent geodetical investigations. Because the magnitude of stress can’t be measured within the earth crust, such heterogeneous stress change calculations are crucial.