Telah dilakukan studi analisis hazard kegempaan dengan menggunakan metode Gaussian dan nilai PGA yang berpengaruh di wilayah Yogyakarta. Penulisan ini bertujuan untuk melakukan analisis hazard gempa  dengan  pemodelan  3-dimensi  untuk  periode  ulang  500  tahun.  Analisis  probabilitas  gempa dilakukan pada T=0 detik, T=0,2 detik (periode pendek) dan T=1 detik (periode panjang) yang dibagi per-grid. Model katalog gempa yang digunakan dalam analisis ini adalah data gempa dari tahun tahun 1937 sampai tahun 2011. Skala magnitudo yang digunakan magnitudo ⥠5, dan kedalaman maksimum 60 km, dengan koordinat 110°10â -110°80â BT dan 7,10°-7,25° LS yang diambil dari BMKG. Hasil penelitian menunjukkan bahwa untuk probabilitas Yogyakarta menghasilkan nilai PGA pada grid 1 antara 0,00848g - 0,00142g, pada grid 2 antara 0,01006g - 0,00136g, pada grid 3 antara 0,07677g - 0,01049g dan pada grid 4 berkisar 0,16825g - 0,01102g. untuk probabilitas terlampaui 10% dalam 50 tahun umur bangunan atau setara dengan periode ulang gempa 500 tahun. Kata kunci : Seismik hazard, PGA, metode Gaussian, Magnitudo