Sistem Sungai Purba Sundalandia: Merekonstruksi Jaringan Drainase Terendam dari Deglasiasi Terakhir

Penelitian oleh Dhani Irwanto, 13 Oktober 2025

Abstrak

Studi ini merekonstruksi sistem drainase purba Paparan Sunda—yang kini sebagian besar terendam di bawah Laut Jawa, Karimata, dan Cina Selatan—menggunakan data batimetri dan topografi resolusi tinggi. Integrasi GEBCO 2025 (15 detik-busur), SRTM v3 (1 detik-busur), dan model genangan deglasial mengungkap enam sistem sungai purba utama dan sebuah danau purba besar di Teluk Thailandia. Pemodelan Daerah Aliran Sungai (DAS) dilakukan dengan ambang batas geomorfik yang konsisten (DAS minimum ≥ 1.000 km²; panjang sungai ≥ 10 km) setelah mengoreksi artefak lintasan kapal dalam data GEBCO. Jaringan yang dihasilkan menggambarkan lanskap fluvial yang saling terhubung yang pernah menghubungkan daratan Sumatra, Jawa, Kalimantan, dan Semenanjung Malaya yang terpapar. Rekonstruksi ini menjelaskan arsitektur hidrologi purba yang menyusun koridor ekologi, transportasi sedimen, dan pergerakan manusia awal di Sundalandia Pleistosen Akhir.

Kata kunci: Sundalandia, sungai purba, deglasiasi, GEBCO 2025, pemodelan DAS, danau purba Teluk Thailandia, Sungai Molengraaff, hidrologi Pleistosen

1. Pendahuluan

Selama Glasial Maksimum Terakhir (LGM) dan deglasiasi berikutnya, Paparan Sunda merupakan salah satu dataran terpapar terbesar di dunia, menyatukan pulau Sumatra, Jawa, Kalimantan, dan Semenanjung Malaya. Naiknya permukaan laut lebih dari 120 m sejak ~21 ka BP secara progresif menenggelamkan platform benua ini, memecahnya menjadi kepulauan Indonesia dan Malaysia saat ini. Merekonstruksi sistem sungai purba sangat penting untuk memahami pola perutean air tawar dan sedimen, konektivitas ekologis dan biogeografis sebelum isolasi, dan respons sistem fluvial tropis terhadap transgresi pasca-glasial yang cepat. Karya-karya sebelumnya (Molengraaff 1921; Voris 2000; Sathiamurthy & Voris 2006) menguraikan peta drainase umum Sundalandia, tetapi mengandalkan data batimetri kasar. Dengan peningkatan terbaru dalam model elevasi digital, sekarang memungkinkan untuk menggambarkan saluran dan cekungan pada skala benua dengan realisme yang lebih besar. Makalah ini memperluas rekonstruksi sebelumnya tentang perubahan permukaan laut relatif (Irwanto 2025a), evolusi suhu permukaan laut (Irwanto 2025b), dan tingkat genangan deglasial (Irwanto 2025c) dengan memetakan jaringan hidrologi purba yang mengatur bekas daratan Sunda.

2. Data dan Metode

2.1 Sumber Data

Untuk mencapai rekonstruksi kerangka drainase purba yang realistis di seluruh Sundalandia, studi ini mengintegrasikan kumpulan data medan global beresolusi tertinggi yang tersedia untuk umum. Data batimetri dan topografi dipilih karena cakupan spasialnya yang saling melengkapi—bawah air dan terestrial—dan karena kompatibilitasnya dalam kerangka geodetik yang seragam. Sumber-sumbernya dirangkum di bawah ini.

1. Batimetri GEBCO 2025 (grid 15 detik-busur), mewakili data kedalaman laut global yang dihimpun dari survai multibeam dan altimetri satelit.
2. SRTM Versi 3 (grid 1 detik-busur), menyediakan elevasi daratan dengan akurasi tinggi yang diperoleh dari interferometri radar.
3. Dataset Genangan Deglasial Sundalandia (Irwanto 2025), sebelumnya diproduksi dari pemodelan permukaan laut, menyediakan permukaan referensi garis pantai untuk interpretasi hidrologi purba.

Semua kumpulan data diambil sampelnya kembali dan dimosaik menjadi model elevasi berkelanjutan yang direferensikan ke koordinat geografis WGS 84 untuk memastikan keselarasan vertikal dan horizontal yang konsisten.

2.2 Metodologi

Rekonstruksi sistem sungai purba Sundalandia mengikuti serangkaian analisis geomorfometri dan hidrologi dalam lingkungan SIG. Model elevasi digital (DEM) komposit dihasilkan melalui mosaik dan resampling batimetri GEBCO 2025 (15 detik-busur) dan topografi SRTM v3 (1 detik-busur) ke dalam grid geografis WGS 84 yang seragam. Permukaannya diratakan secara hidro untuk menghilangkan diskontinuitas di sepanjang garis pantai modern dan memastikan rute aliran yang konsisten melintasi antarmuka subaerial-submarine.

2.2.1 Minimalisasi Artefak

Anomali seperti parit batimetri—dikenal sebagai artefak lintasan kapal—diidentifikasi secara visual sebagai depresi linier yang sejajar dengan jalur survai di grid GEBCO. Anomali ini dikoreksi secara lokal melalui penyaringan median lingkungan dan penyuntingan manual simpul grid yang menyimpang menggunakan interpolasi bilinear. Tujuannya adalah untuk menekan jalur aliran silang buatan yang dapat mendistorsi konektivitas hidrologis sekaligus mempertahankan variabilitas topografi yang sebenarnya. Penghapusan artefak semacam itu secara menyeluruh tidak mungkin dilakukan tanpa data pengukuran multibeam asli; namun, pengaruh hidrologisnya diminimalkan hingga tingkat yang dapat diabaikan.

2.2.2 Pemodelan Hidrologi

Delineasi Daerah Aliran Sungai (DAS) dan ekstraksi aliran sungai dilakukan menggunakan algoritma akumulasi aliran standar dalam perangkat lunak SIG. Arah aliran diturunkan dari DEM yang telah dikoreksi menggunakan algoritma D8[1], diikuti dengan perhitungan akumulasi aliran dan orde aliran sungai. Ambang batas area kontribusi minimum sebesar 1.000 km² ditetapkan untuk aliran orde pertama, dan panjang saluran minimum 10 km diadopsi untuk mengecualikan drainase semu atau ephemeral. Jaringan aliran sungai yang dihasilkan kemudian divektorisasi dan divalidasi secara topologi untuk memastikan konektivitas dan hierarki drainase yang realistis.

2.2.3 Integrasi dengan Drainase Modern

Saluran purba yang dimodelkan diselaraskan secara spasial dengan muara sungai modern untuk menjaga kontinuitas genetik antara daerah aliran sungai subaerial dan bawah laut. Prosedur ini melibatkan penyesuaian jalur aliran terminal menuju muara dan muka delta yang ada, berdasarkan gradien hidrologi dan arah transpor sedimen yang disimpulkan dari analisis kemiringan dan kelengkungan. Hal ini memastikan bahwa sistem drainase purba yang dimodelkan tetap sesuai dengan muara sungai dan batas fisiografis yang ada.

2.2.4 Sintesis dan Visualisasi Hidrologi

Mosaik drainase akhir dikategorikan berdasarkan identitas cekungan dan diekspor sebagai vektor shapefile untuk visualisasi kartografi. Enam sistem skala regional didefinisikan melalui pengelompokan hierarkis zona akumulasi aliran, yang sesuai dengan sistem Laut Jawa, Laut Jawa Timur, Selat Karimata, Teluk Thailandia, Perpanjangan Mekong, dan Selat Malaka. Hasilnya dibandingkan dengan kontur batimetri regional dan rekonstruksi garis pantai deglasial (Irwanto 2025a; 2025c) untuk memvalidasi koherensi drainase di bawah permukaan laut -122 m (~22,5 ribu tahun lalu), yang sesuai dengan Maksimum Glasial Terakhir (LGM).

2.3 Keterbatasan

Sistem sungai purba yang direkonstruksi merepresentasikan model geomorfometri orde pertama yang dibatasi terutama oleh topografi dan batimetri. Beberapa proses geomorfik dan dinamis tidak diikutsertakan secara eksplisit karena kelangkaan dan inkonsistensi data regional. Akibatnya, hasilnya harus diinterpretasikan sebagai kerangka hidrologi umum, alih-alih geometri kanal purba yang eksak.

1. Dinamika sedimen — Proses-proses seperti progradasi delta, avulsi saluran, agradasi dataran banjir, dan pergeseran litoral tidak dimodelkan. Faktor-faktor ini dapat secara substansial mengubah morfologi lembah dan geometri muara seiring waktu, terutama selama transgresi tahap akhir.
2. Pengaruh bawah permukaan dan tektonik — Penurunan tektonik lokal, reaktivasi sesar, dan pengangkatan diferensial mungkin telah mengubah gradien drainase dan bentuk cekungan setelah pembentukan saluran awal. Pengaruh-pengaruh ini masih belum terkuantifikasi pada skala paparan.
3. Karst dan medan pelarutan — Di wilayah yang dilapisi oleh litologi karbonat (misalnya, Jawa utara, Kalimantan barat, dan batas Thailand–Melaya), drainase bawah permukaan dan pengembangan sinkhole mungkin telah memengaruhi konektivitas daerah tangkapan air dengan cara yang tidak tertangkap oleh model aliran berbasis permukaan.
4. Pemadatan sedimen dan penyesuaian isostatik — Penurunan pasca pengendapan dan pemulihan isostatik menyusul deglasiasi tidak diintegrasikan ke dalam koreksi DEM, sehingga menimbulkan ketidakpastian kecil pada elevasi absolut dan tingkat dasar relatif.
5. Kualitas data batimetri — Meskipun artefak telah diminimalkan, anomali lintasan kapal yang tersisa di grid GEBCO mungkin masih memengaruhi rute aliran lokal. Artefak ini—depresi seperti parit linear yang terbentuk selama grid data—telah dikurangi tetapi tidak dapat sepenuhnya dihilangkan tanpa akses ke data pengukuran kedalaman asli. Dampaknya diminimalkan pada skala regional tetapi mungkin tetap ada secara lokal.
6. Penyederhanaan temporal — Pemodelan ini mengasumsikan topografi kuasi-statis yang sesuai dengan satu permukaan laut acuan (−122 m RSL, ~22,5 ribu tahun lalu). Oleh karena itu, migrasi garis pantai yang progresif, redistribusi sedimen, dan reorganisasi hidrologis selama ribuan tahun berikutnya berada di luar cakupan rekonstruksi ini.

Secara keseluruhan, ketidakpastian di atas kemungkinan besar tidak akan mengubah konfigurasi umum enam sistem drainase utama yang diidentifikasi dalam studi ini, tetapi dapat memengaruhi posisi saluran lokal dan detail anak sungai. Penelitian selanjutnya yang mengintegrasikan stratigrafi seismik, inti sedimen, dan batimetri beresolusi lebih tinggi dapat semakin menyempurnakan realisme hidrologi purba rekonstruksi Sundalandia.

3. Hasil

3.1 Sistem Sungai Purba Utama

Pemodelan hidrologi menunjukkan jaringan koheren dari enam sistem drainase utama yang menempati Paparan Sunda sebelum penggenangan Holosen. Setiap sistem mengintegrasikan banyak anak sungai yang mengalir dari daratan Sumatra, Jawa, Kalimantan, dan Semenanjung Malaya. Besaran relatif dan wilayah penyumbangnya dirangkum dalam Tabel 1, sementara konfigurasi spasialnya diilustrasikan pada Gambar 1.

Gambar 1. Sistem sungai purba yang direkonstruksi dan daerah aliran sungai utama di Paparan Sunda. Arsir biru pucat menunjukkan luas genangan akibat glasial; biru tua menandai danau purba di Teluk Thailandia. Garis hitam menunjukkan sungai purba yang dimodelkan, garis abu-abu menggambarkan sungai modern. © 2025 Dhani Irwanto.

Tabel 1. Sistem Sungai Purba Utama di Sundalandia

Sistem	Wilayah Sumber Utama	Perkiraan Luas Daerah Aliran Sungai (km²)
Laut Jawa	Kalimantan Selatan, Jawa Utara, Sumatera Selatan	≈ 570.000
Laut Jawa Timur	Kalimantan Selatan (Barito, Kapuas-Murung, Kahayan)	≈ 180.000
Selat Karimata (Molengraaff)	Sumatera Timur, Kalimantan Barat	≈ 630.000
Teluk Thailandia	Semenanjung Melayu Timur, Cekungan Chao Phraya	≈ 1.020.000
Perpanjangan Mekong	Mekong Hilir dan tepi Laut Cina Selatan yang berdekatan	≈ 690.000
Selat Malaka	Sumatera Timur, Semenanjung Malaya Barat	≈ 260.000

Tabel 2. Konektivitas antara Sungai Modern dan Sistem Sungai Purba

Sistem Sungai Purba	Sungai Modern — Kalimantan	Sungai Modern — Sumatra	Sungai Modern — Jawa	Sungai Modern — Semenanjung Malaya/Daratan
Laut Jawa	Mendawai, Sampit, Pembuang	Tulang Bawang	Bengawan Solo, Serang, Cimanuk, Citarum	–
Laut Jawa Timur	Barito, Kapuas-Murung, Kahayan	–	–	–
Selat Karimata (Molengraaff)	Kapuas	Musi, Batanghari, Indragiri	–	–
Teluk Thailandia	–	–	–	Johor, Rompin, Endau, Kuantan, Kelantan, Tapi, Mae Klong, Chao Phraya, Bang Pakong
Perpanjangan Mekong	–	–	–	Sungai Mekong dan anak-anak sungainya
Selat Malaka	–	Kampar, Rokan, Barumun, Belawan	–	Malaka, Perak

Catatan: Tabel 2 menggambarkan kesinambungan antara muara sungai masa kini dan saluran purba yang dimodelkan, mendukung kesimpulan bahwa banyak muara modern berasal sebagai segmen terminal dari sistem purba ini.

3.2 Danau Purba di Teluk Thailandia

Depresi tertutup seluas kurang lebih 93.000 km² dengan ambang outlet mendekati −55 m dari permukaan laut relatif menunjukkan keberadaan danau purba yang luas di bagian tengah Teluk Thailandia. Morfologi cekungan menunjukkan retensi air tawar yang berkepanjangan selama tahap deglasial awal sebelum akhirnya meluap dan menerobos ke Laut Cina Selatan. Fase danau ini sesuai dengan catatan sedimen transisi yang didokumentasikan dalam studi inti regional (Horton dkk., 2005; Chabangborn dkk., 2020; Zhang dkk., 2022)[2] yang menunjukkan pergeseran dari fasies yang didominasi air tawar menuju lingkungan pengendapan yang lebih dipengaruhi estuari atau laut selama kenaikan permukaan laut.

Sebagai contoh, inti sedimen di sepanjang Teluk Thailandia bagian barat (misalnya, CP3, CP4, CP5) mencatat transisi stratigrafi yang konsisten dengan meningkatnya pengaruh laut sekitar 7,9 ribu tahun yang lalu, sebagaimana dibuktikan oleh rasio ukuran butir, kumpulan mikrofosil, dan masuknya serbuk sari bakau (Chabangborn dkk., 2020; Horton dkk., 2005). Temuan-temuan ini memberikan dukungan independen untuk transformasi air tawar-ke-muara yang luas yang sesuai dengan hidrologi danau purba yang dimodelkan dalam studi ini.

3.3 Efek Pengurangan Artefak

Penggambaran jalur aliran yang akurat melintasi landas kontinen dangkal memerlukan koreksi depresi buatan yang menyerupai parit yang dihasilkan melalui pemetaan data jalur kapal. Setelah penghalusan iteratif batimetri GEBCO 2025, beberapa perbaikan dicapai pada topologi drainase yang dimodelkan. Koreksi ini menghasilkan alinyemen lembah yang konsisten secara morfologis dan konektivitas yang lebih realistis antar cekungan yang berdekatan. Penyesuaian yang paling penting diuraikan di bawah ini:

1. Sistem Laut Jawa meluas ke arah barat, mengintegrasikan anak sungai Sumatera selatan yang sebelumnya salah arah menuju Selat Sunda.
2. Sistem Laut Jawa Timur yang berbeda muncul, mengisolasi daerah aliran sungai Barito, Kapuas-Murung, dan Kahayan dari cekungan Laut Jawa.
3. Saluran melintang linier yang sebelumnya dihasilkan oleh artefak jejak kapal telah dihilangkan, memulihkan drainase lengkung alami.

4. Diskusi

4.1 Signifikansi Paleogeografis

Jaringan yang direkonstruksi menunjukkan bahwa Paparan Sunda pernah berfungsi sebagai dataran fluvial yang bersebelahan. Zona konvergensi drainase di Laut Jawa, Selat Karimata, dan Teluk Thailandia selaras dengan deposentrum yang teridentifikasi dalam survai seismik dan analisis inti sedimen (Horton dkk., 2005; Chabangborn dkk., 2020). Hubungan ini memperjelas peran paparan sebagai penyerap sedimen sekaligus koridor pembuangan air tawar selama deglasiasi, memberikan konteks fisik bagi pola transgresi cepat dan fragmentasi garis pantai yang dijelaskan dalam Irwanto (2025c). Danau purba Teluk Thailandia dan transgresi lautnya selanjutnya merupakan contoh transisi dinamis dari lingkungan terestrial ke lingkungan laut, yang menggambarkan kontinuitas sedimen antara sistem fluvial purba dan cekungan paparan modern.

4.2 Bukti Biogeografis

Distribusi modern taksa air tawar dan estuari menyiratkan adanya kesinambungan historis di sepanjang jalur perairan purba ini. Ikan kurau air tawar (Polydactylus macrophthalmus), yang saat ini terbatas di Sungai Kapuas (Kalimantan) dan Musi–Batanghari (Sumatera), merupakan contoh pemisahan vikariat populasi yang pernah bergabung dengan sistem Sungai Molengraaff (Motomura dkk., 2001). Disjungsi yang serupa antara spesies mangrove dan moluska akuatik memperkuat hubungan paleohidrologi yang disimpulkan dari model ini.

4.3 Perbandingan dengan Model Sebelumnya

Dibandingkan dengan Voris (2000) dan Sathiamurthy & Voris (2006), rekonstruksi ini menawarkan fidelitas spasial yang lebih tinggi dan realisme hidrologis yang lebih baik. Resolusi GEBCO 2025 yang lebih halus menggambarkan meander dan kelengkungan anak sungai yang sebelumnya tidak terpecahkan, sementara koreksi artefak meningkatkan kontinuitas drainase. Sistem yang dihasilkan menunjukkan cekungan asimetris dan delta multi-cabang yang lebih konsisten dengan morfodinamika aluvial tropis.

4.4 Implikasi Hidrologi

Dataran rendah bergradien luas yang disimpulkan dari model menunjukkan sungai-sungai yang berkelok-kelok perlahan melintasi dataran banjir yang luas, mampu menopang lahan basah dan kompleks delta yang luas. Saluran-saluran ini kemungkinan mengangkut muatan sedimen yang besar menuju tepi paparan, memengaruhi dinamika nutrisi di dekat pantai dan akhirnya membentuk rangkaian punggungan terendam yang terlihat dalam batimetri saat ini.

4.5 Implikasi yang Lebih Luas bagi Sejarah Manusia dan Biotik

Selama penurunan permukaan laut, koridor sungai terpadu di Sundalandia menyediakan air tawar yang berkelanjutan, tanah yang subur, dan rute pelayaran melintasi paparan yang muncul. Koridor-koridor tersebut akan memfasilitasi penyebaran kelompok-kelompok manusia, memungkinkan pendudukan cekungan pedalaman dan tepi pantai jauh sebelum transgresi laut. Dataran-dataran sungai yang saling terhubung mungkin telah berfungsi sebagai urat nadi pertukaran budaya dan genetik di wilayah yang sekarang dikenal sebagai Kepulauan Asia Tenggara.

Jalur aluvial yang luas di sepanjang sistem Kapuas–Musi–Batanghari purba (Molengraaff) dan Teluk Thailandia memiliki kapasitas ekologis untuk menopang komunitas proto-pertanian. Lingkungan ini mencerminkan tatanan lingkungan peradaban sungai di tempat lain, menunjukkan bahwa Paparan Sunda menawarkan peluang serupa bagi perilaku produksi pangan dan permukiman awal, sebagaimana dibahas dalam bagian Peradaban Sungai di Irwanto (2015).

Setelah tergenang secara progresif, sungai-sungai induk berevolusi menjadi muara pesisir dan dataran delta, mempertahankan perannya sebagai poros komunikasi. Transformasi jaringan transportasi fluvial menjadi muara kemungkinan mendorong munculnya pengetahuan hidrolika dan navigasi, yang mendorong transisi dari budidaya di pedalaman ke eksploitasi sumber daya maritim.

Ketika penggenangan landas kontinen memutus rute benua, komunitas manusia beradaptasi dengan naiknya permukaan air dengan beralih ke mata pencaharian pesisir. Lembah-lembah sungai yang dulunya terlindung kini menjadi teluk dan selat yang terlindung—jalur alami bagi pelayaran awal. Paksaan lingkungan ini mungkin telah menumbuhkan orientasi maritim yang kemudian menjadi ciri budaya Austronesia dan budaya Asia Tenggara awal lainnya.

Perendaman Paparan Sunda memecah habitat yang dulunya berkesinambungan, mengisolasi spesies air tawar dan darat. Distribusi Polydactylus macrophthalmus yang terbagi di Sumatra dan Kalimantan (Motomura dkk., 2001) merupakan ciri vikariansi pasca-genangan. Proses serupa kemungkinan besar memengaruhi gajah di Kalimantan (Fernando dkk., 2003; Sharma dkk., 2018), kura-kura air tawar, dan vegetasi sungai, menghasilkan mosaik biogeografis yang terlihat saat ini.

4.6 Bukti Empiris Transisi Danau ke Muara di Teluk Thailandia

Berbagai studi inti sedimen dari Teluk Thailandia dan dataran pantai di sekitarnya memperkuat interpretasi tahap paleo-lacustrine yang diikuti oleh pengaruh laut progresif selama awal hingga pertengahan Holosen. Inti sedimen dari Teluk Thailandia bagian barat (CP3, CP4, CP5; Chabangborn dkk., 2020; Jiwarungrueangkul dkk., 2022) mengungkapkan suksesi stratigrafi di mana satuan lacustrine dan delta berbutir halus dilapisi oleh fasies estuari payau hingga laut, disertai peningkatan serbuk sari bakau, kelimpahan foraminifera, dan mikrofosil laut.

Demikian pula, Horton dkk. (2005) mendokumentasikan transisi lingkungan purba yang serupa di inti pesisir Semenanjung Melaya-Thailandia, dengan evolusi yang jelas dari rawa air tawar dan endapan fluvial menjadi sedimen dataran pasang surut dan estuari yang sinkron dengan kenaikan muka air laut pertengahan Holosen. Sintesis regional sedimentasi paparan (Zhang dkk., 2022) selanjutnya menunjukkan bahwa Paparan Sunda mengalami reorganisasi hidrologi yang meluas, di mana cekungan yang sebelumnya berada di bawah permukaan tanah berubah menjadi estuari yang tenggelam dan teluk laut dangkal seiring dengan kenaikan muka air laut yang cepat antara sekitar 10 dan 7 ribu tahun yang lalu.

Secara kolektif, kumpulan data ini memperkuat interpretasi bahwa depresi Teluk Thailandia awalnya berfungsi sebagai cekungan air tawar besar dan kemudian beralih menjadi teluk laut semi-tertutup—urutan yang konsisten dengan topografi yang dimodelkan dan rekonstruksi hidrologi yang disajikan dalam studi ini.

5. Kesimpulan

Analisis terpadu ini mengidentifikasi enam sistem sungai purba utama dan sebuah danau purba besar di Teluk Thailandia yang bersama-sama membentuk kerangka hidrologi Sundalandia yang muncul. Dengan menggabungkan batimetri GEBCO 2025, topografi SRTM v3, dan pemodelan hidrologi, studi ini menyempurnakan rekonstruksi sebelumnya dan menetapkan penggambaran arsitektur drainase paparan yang konsisten secara fisik. Selain geomorfologi, temuan ini menjelaskan bagaimana jaringan fluvial ini menyusun koridor ekologi dan jalur manusia sebelum transgresi Holosen, yang meletakkan dasar bagi penelitian interdisipliner di masa mendatang.

Referensi

Molengraaff, G.A.F. (1921). Modern Deep-Sea Research in the East Indian Archipelago.

Voris, H.K. (2000). Maps of Pleistocene Sea Levels in Southeast Asia. Journal of Biogeography, 27, 1153–1167.

Sathiamurthy, E., & Voris, H.K. (2006). Maps of Holocene Transgression and Pleistocene Coastlines, Southeast Asia.

Chabangborn, A., Phantuwongraj, S., Sinsakul, S., Choowong, M., & Nakagawa, T. (2020). Environmental changes on the west coast of the Gulf of Thailand during the Holocene. Quaternary International, 555, 3–16. https://doi.org/10.1016/j.quaint.2019.12.020

Horton, B. P., et al. (2005). Holocene sea levels and palaeoenvironments, Malay–Thai Peninsula. The Holocene, 15(8), 1189–1203. https://doi.org/10.1191/0959683605hl887rp

Zhang, H., Liu, S., Wu, K., Cao, P., Pan, H-J., Wang, H., … Shi, X. (2022). Evolution of sedimentary environment in the Gulf of Thailand since the last deglaciation. Quaternary International, 629, 36–43. https://doi.org/10.1016/j.quaint.2021.02.018

Jiwarungrueangkul, T., Jirapinyakul, A., Sompongchaiyakul, P., & Rattanakom, R. (2022). Response of sediment grain size to sea-level rise during the middle Holocene on the west coast of the Gulf of Thailand. Arabian Journal of Geosciences, 15, 167. https://doi.org/10.1007/s12517-022-09450-3

Motomura, H., et al. (2001). Redescription of a rare threadfin (Perciformes: Polynemidae), Polydactylus macrophthalmus (Bleeker, 1858), with designation of a lectotype and notes on distributional implications. Ichthyological Research, 48, 289–294.

Fernando, P., Vidya, T.N.C., Payne, J., Stuewe, M., Davison, G., Alfred, R.J., Andau, P., Bosi, E., Kilbourn, A., & Melnick, D.J. (2003). DNA analysis indicates that Asian elephants are native to Borneo and are therefore a high priority for conservation. PLoS Biology, 1(1), 110–115.

Sharma, R., Goossens, B., Heller, R., Rasteiro, R., Othman, N., Bruford, M.W., & Chikhi, L. (2018). Genetic analyses favour an ancient and natural origin of elephants on Borneo. Scientific Reports, 8, 880.

World Wildlife Fund (WWF). (n.d.). Borneo Pygmy Elephant. Retrieved from http://www.worldwildlife.org/species/borneo-pygmy-elephant

Irwanto, D. (2015). Sundaland: Tracing the Cradle of Civilizations. Sections “Riverine Civilizations” and “Kalimantan Elephants.”

Irwanto, D. (2025a). A Refined Relative Sea-Level Curve for Sundaland.

Irwanto, D. (2025b). Holocene and Deglacial Sea-Surface Temperatures in Sundaland.

Irwanto, D. (2025c). Deglacial Rapid Inundation and Land-Loss Rates of Sundaland.

Catatan Kaki

[1] Algoritma D8 (Deterministic Eight-node) merupakan model arah aliran standar dalam analisis SIG hidrologi. Algoritma ini menetapkan setiap sel raster satu arah kemiringan menurun menuju salah satu dari delapan sel tetangganya—utara, timur laut, timur, tenggara, selatan, barat daya, barat, atau barat laut—berdasarkan gradien penurunan paling curam. Pendekatan ini memungkinkan komputasi akumulasi aliran dan deliniasi DAS yang efisien di seluruh kumpulan data medan yang besar.

[2] Bukti inti langsung untuk sebuah “danau” air tawar yang berkesinambungan yang mencakup seluruh cekungan Teluk Thailandia adalah terbatas; oleh karena itu, interpretasi danau purba harus dianggap sebagai hipotesis geomorfologi kerja yang berasal dari topografi yang dimodelkan dan potensi hidrologi, dibatasi oleh analog sedimen regional.