BAB 1 PENDAHULUAN

1.1 Latar Belakang

Secara geotektonik, Indonesia berada di jalur pertemuan tiga lempeng besar Indo-Australia, Eurasia, dan Pasifik yang menjadikannya salah satu negara dengan tingkat aktivitas gempa tertinggi di dunia. Kondisi tersebut menyebabkan Pulau Sumatera, khususnya Provinsi Aceh, menjadi wilayah yang sangat aktif secara seismik. Aceh terletak di dekat zona subduksi Sunda Megathrust dan juga dilalui oleh segmen Aceh dari Sesar Sumatera, dua sumber utama yang sering memicu gempa dengan magnitudo besar. Kajian seismisitas yang dilakukan oleh BMKG menunjukkan bahwa wilayah Aceh memiliki nilai parameter a dan b pada persamaan Gutenberg–Richter yang menunjukkan frekuensi kejadian gempa relatif tinggi dan berpotensi signifikan terhadap bahaya geotektonik di masa depan(Lusiani dkk, 2020).

Selain dipengaruhi oleh aktivitas subduksi, struktur sesar aktif di daratan Aceh juga memberikan kontribusi besar terhadap karakteristik gempa lokal. Penelitian geologi terbaru menunjukkan bahwa segmen Aceh pada Great Sumatran Fault (GSF) terdiri dari beberapa subsegmen yang memiliki orientasi berbeda, menghasilkan pola deformasi kompleks dan potensi guncangan kuat di sekitar Banda Aceh hingga Lhokseumawe. Analisis distribusi spasial gempa di sepanjang jalur ini penting untuk mengenali zona-zona dengan konsentrasi kejadian tinggi serta memetakan pola klaster seismik yang dapat digunakan dalam mitigasi risiko. Studi klasterisasi gempa di Pulau Sumatera menunjukkan bahwa wilayah Aceh termasuk dalam zona dengan kepadatan epicenter paling tinggi di Indonesia bagian barat. (Hady & Marliyani,2022).

Analisis spasial menjadi pendekatan yang efektif untuk menggambarkan persebaran dan pola kejadian gempa bumi, data spasial umumnya terbagi menjadi dua jenis utama, yaitu vektor (titik, garis, dan poligon) dan raster, di mana data titik digunakan untuk merepresentasikan lokasi peristiwa seperti titik pusat gempa (Modul Statistika Spasial, Pertemuan 1, 2025). Dengan perangkat lunak R dengan paket seperti sf, ggplot2, dan rnaturalearth, data koordinat gempa dapat diolah menjadi objek spasial dan divisualisasikan ke dalam peta informatif. Visualisasi ini membantu memahami hubungan antara lokasi epicenter, jalur sesar, serta kondisi geologi regional, sehingga menjadi dasar dalam perencanaan mitigasi bencana. Praktikum ini bertujuan melatih mahasiswa dalam pengolahan data spasial dan pembuatan peta sebaran gempa bumi di Aceh, sekaligus memperkuat pemahaman tentang penerapan sistem informasi geografis (SIG) dalam studi geofisika dan kebencanaan(Bivand dkk, 2013).

1.2 Tujuan

Mahasiswa dapat membuat peta dengan paket pemrograman R
Mahasiswa dapat membuat peta dengan data titik spasial pada R

1.3 Rumusan Masalah

Bagaimana membuat peta dengan paket pemrograman R
Bagaiman membuat peta dengan data titik spasial pada R

1.4 Batasan Masalah

Buatlah peta epicenter gempa bumi di Indonesia pada tahun 2024 (gunakan packages R yang dianggap gampang, setiap praktikan menggambarkan peta provinsi yang berbeda).
Buatlah peta semua provinsi di Indonesia dengan memberikan warna yang berbeda pada setiap kab/kota (setiap praktikan membuat peta provinsi yang berbeda).

BAB 2 PEMBAHASAN

2.1. Peta tanpa file SHP

Download package

library(sf)

## Warning: package 'sf' was built under R version 4.4.3

## Linking to GEOS 3.13.0, GDAL 3.10.1, PROJ 9.5.1; sf_use_s2() is TRUE

library(ggplot2)
library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(rnaturalearth)

## Warning: package 'rnaturalearth' was built under R version 4.4.3

library(rnaturalearthdata)

## Warning: package 'rnaturalearthdata' was built under R version 4.4.3

## 
## Attaching package: 'rnaturalearthdata'

## The following object is masked from 'package:rnaturalearth':
## 
##     countries110

library(viridis)

## Loading required package: viridisLite

Input data gempa

setwd("F:/Spasial ")
data <- read.csv("Dataaceh.csv")

Konversi menjadi data spasial (sf)

gempa_sf <- st_as_sf(data, coords = c("longitude", "latitude"), crs = 4326)

Filter peta shp menjadi provinsi Aceh

indo <- ne_states(country = "indonesia", returnclass = "sf")
aceh <- indo %>% filter(name == "Aceh")

buat plot peta sebaran gempa tanpa shp hanya packages dari sf

ggplot() +
  geom_sf(data = aceh, fill = "gray", color = "black") +
  geom_sf(data = gempa_sf, aes(size = mag, color = depth), alpha = 0.7) +
  scale_color_viridis(option = "plasma", name = "Kedalaman (km)") +
  scale_size_continuous(name = "Magnitudo", range = c(2, 6)) +
  labs(title = "Sebaran Gempa di Provinsi Aceh tahun 2024",
       subtitle = "Sumber: USGS",
       x = "Longitude", y = "Latitude",
       caption = "Visualisasi: ggplot2 + sf") +
  theme_minimal() +
  theme(plot.title = element_text(face = "bold", size = 14),
        legend.position = "right")

2.2 Peta sebaran gempa dengan SHP

Input data

setwd("F:/Spasial")
indo_shp<-st_read("gadm41_IDN_2.shp")

## Reading layer `gadm41_IDN_2' from data source `F:\Spasial\gadm41_IDN_2.shp' using driver `ESRI Shapefile'
## Simple feature collection with 502 features and 13 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: 95.00971 ymin: -11.00761 xmax: 141.0194 ymax: 6.076941
## Geodetic CRS:  WGS 84

Filter peta shp menjadi provinsi Aceh

aceh_shp<-indo_shp %>% filter(NAME_1 == "Aceh")
aceh_shp

Plot hasil sebaran gempa dengan SHP

ggplot(aceh_shp) +
  geom_sf(aes(fill = NAME_2), color = "black", size = 0.3) +
  geom_sf(data = gempa_sf, aes(size = mag, color = depth), alpha = 0.7) +
  guides(fill = "none") +
  scale_color_viridis(option = "plasma", name = "Kedalaman (km)") +
  scale_size_continuous(name = "Magnitudo", range = c(2, 6)) +
  labs(title = "Sebaran Gempa di Provinsi Aceh Tahun 2024",
       subtitle = "Sumber: USGS",
       x = "Longitude", y = "Latitude",
       caption = "Visualisasi: ggplot2 + sf") +
  theme_minimal() +
  theme(plot.title = element_text(face = "bold", size = 14),
        legend.position = "right")

2.3 Pembahasan

Pada plot 1 menunjukkan bahwa aktivitas gempa bumi terjadi di berbagai wilayah, terutama di bagian barat dan selatan Aceh, termasuk di sekitar Pulau Simeulue yang berada di lepas pantai Samudra Hindia. Sebagian besar titik gempa memiliki kedalaman kurang dari 75 km, yang ditunjukkan oleh warna ungu hingga biru, menandakan bahwa gempa tersebut tergolong dangkal dan berpotensi menimbulkan dampak signifikan di permukaan, terutama jika terjadi di wilayah laut. Magnitudo gempa yang terekam berkisar antara 4,0 hingga 5,6 dengan mayoritas berada pada kisaran 4,0–5,0, sehingga dikategorikan sebagai gempa ringan hingga sedang. Sebaran titik gempa memperlihatkan bahwa aktivitas tektonik lebih dominan di sepanjang zona subduksi lempeng Indo-Australia dan Eurasia yang melintas di barat Aceh. Kondisi ini menunjukkan bahwa wilayah pesisir barat dan selatan Aceh memiliki tingkat kerawanan gempa yang tinggi, sehingga perlu tetap waspada terhadap potensi gempa bumi dan tsunami. Pada plot 2 menunjukkan bahwa aktivitas seismik tersebar merata di sepanjang pantai barat dan selatan, terutama di wilayah Samudra Hindia bagian barat Sumatra yang berdekatan dengan zona subduksi lempeng Indo-Australia dan Eurasia. Sebagian besar gempa memiliki kedalaman dangkal (kurang dari 50 km) yang ditunjukkan oleh warna ungu hingga biru, sehingga berpotensi menimbulkan dampak signifikan di permukaan dan kemungkinan tsunami apabila terjadi di laut. Magnitudo gempa yang teramati berkisar antara 4,0 hingga 5,6, dengan dominasi gempa ringan hingga sedang serta beberapa kejadian lebih kuat di sekitar wilayah barat daya Aceh dan Pulau Simeulue. Pola ini mengindikasikan bahwa sepanjang tahun 2024, Aceh mengalami aktivitas tektonik yang cukup tinggi di zona pertemuan lempeng aktif, sehingga wilayah pesisir barat dan selatan perlu tetap waspada terhadap potensi gempa bumi dan tsunami.

BAB 3 KESIMPULAN

Berdasarkan hasil praktikum analisis peta sebaran gempa di Provinsi Aceh tahun 2024, dapat disimpulkan bahwa wilayah Aceh termasuk daerah dengan aktivitas seismik yang tinggi. Sebaran gempa didominasi di bagian barat dan selatan Aceh, terutama di sekitar perairan Pulau Simeulue yang berada di zona subduksi antara lempeng Indo-Australia dan Eurasia. Mayoritas gempa yang terjadi memiliki magnitudo antara 4,0 hingga 5,6 dengan kedalaman kurang dari 75 km, yang menunjukkan bahwa sebagian besar gempa bersifat dangkal dan berpotensi menimbulkan dampak di permukaan serta ancaman tsunami apabila berpusat di laut. Hasil ini menggambarkan bahwa aktivitas tektonik di Aceh masih sangat aktif sehingga diperlukan pemantauan seismik secara berkala serta peningkatan kesadaran dan kesiapsiagaan masyarakat terhadap potensi bencana gempa bumi dan tsunami.

DAFTAR PUSTAKA

Modul Statistika Spasial Pertemuan 1. (2025). Visualisasi data spasial dengan R. Program Studi Statistik.
Bivand, R. S., Pebesma, E., & Gómez-Rubio, V. (2013). Applied spatial data analysis with R (2nd ed.). Springer.
Anggraini, D., Indratno, S. W., Mukhaiyar, U., & Puspito, A. N. T. (2023). Cluster Analysis of Sumatra Island Earthquake Distribution. Desimal: Jurnal Matematika.

Laporan 1 Pengantar Statistika Spasial

Silvi Anggraini (F1F022043)

2025-11-04