GLOF: Glacial Lake Outburst Flood พัฒนา ราชวงศ์ อาศรมภูมิวิทยาศาสตร์ สาขาวิชาภูมิศาสตร์ มหาวิทยาลัยนเรศวร เหตุการณ์ภัยพิบัติน้ำท่วมฉับพลันที่เพิ่งเกิดขึ้นในเนปาลเมื่อวานนี้ คือ อุทาหรณ์ที่สะท้อนถึงวิกฤตความสุ่มเสี่ยงของภูมิภาคหิมาลัยได้อย่างชัดเจนที่สุด เมื่อมวลน้ำมหาศาลจากการทะลักของทะเลสาบธารน้ำแข็ง (GLOF: Glacial Lake Outburst Flood) ไหลซัดทำลายล้างโครงสร้างพื้นฐานและชุมชนท้ายน้ำในพริบตา ปรากฏการณ์นี้ไม่ใช่เพียงอุบัติเหตุทางธรรมชาติที่เกิดขึ้นเป็นครั้งคราว หากแต่เป็นผลลัพธ์เชิงระบบของการเปลี่ยนแปลงสภาพภูมิอากาศโลก ที่เร่งให้อ่างกักเก็บน้ำแข็งบนยอดเขาสูงละลายและเกินจุดรับไหว (tipping point) การเข้าใจภัยพิบัติดังกล่าวผ่านกรอบความคิดเชิงระบบ (systems thinking) จึงเป็นหัวใจสำคัญในการมองเห็นความเชื่อมโยงระหว่างการเพิ่มขึ้นของอุณหภูมิโลก การเปลี่ยนรูปของธารน้ำแข็ง และความมั่นคงของมนุษย์ในพื้นที่เปราะบางข้ามพรมแดน ทั้งนี้เมื่อมองปรากฏการณ์ภัยพิบัตินี้ผ่านสายตาของนักคิดเชิงระบบที่เห็นเทือกเขาหิมาลัยเป็นเหมือน มหาเครื่องจักรแห่งธรรมชาติ ที่กำลังสูญเสียสมดุล โดยสภาวะดั้งเดิม ธารน้ำแข็งทำหน้าที่เป็นอ่างกักเก็บ (reservoir) ขนาดใหญ่ที่สุด มันคอยกักเก็บน้ำในรูปของแข็งและค่อยๆ ปลดปล่อยกระแสการไหล (flux) ออกมาเลี้ยงผู้คนนับพันล้านอย่างสม่ำเสมอตามฤดูกาล แต่เมื่อกิจกรรมมนุษย์อัดฉีดก๊าซเรือนกระจกขึ้นสู่ชั้นบรรยากาศ อุณหภูมิโลกที่สูงขึ้นได้เข้าไปเร่งอัตราการละลาย (input flux) ให้เร็วกว่าที่ธรรมชาติจะสร้างหิมะใหม่มาทดแทนได้ น้ำที่ละลายออกมาอย่างรวดเร็ว เริ่มมองหาที่อยู่ใหม่ มันไหลไปสะสมรวมกันเกิดเป็นทะเลสาบธารน้ำแข็ง (glacial lakes) ซึ่งกลายเป็น อ่างกักเก็บใหม่ ที่ตั้งอยู่อย่างหมิ่นเหม่บนมือของเปรตภูเขา โดยมีเพียงเขื่อนหินและน้ำแข็งธรรมชาติคอยกั้นไว้ ตรงนี้เองที่ลูปป้อนกลับเชิงบวก (positive feedback loop) เริ่มทำงานและเร่งวิกฤตให้รุนแรงขึ้น เมื่อธารน้ำแข็งหดตัว พื้นผิวสีขาวที่เคยสะท้อนแสงแดดกลับหายไป เปิดให้เห็นผืนน้ำและแผ่นดินสีเข้มที่ดูดซับความร้อนได้ดีกว่าเดิม ยิ่งเร่งให้ธารน้ำแข็งรอบข้างละลายไวขึ้น และเมื่อมวลน้ำในทะเลสาบขยายตัวและเชื่อมรวมกันเป็นแอ่งใหญ่ขึ้น ผืนน้ำกว้างนี้จะสะสมความร้อนไว้ และกัดเซาะฐานธารน้ำแข็งที่อยู่ติดกัน ดันให้ธารน้ำแข็งถอยร่นและละลายหนักขึ้นไปอีก ระบบถูกบีบให้สุ่มเสี่ยงเดินหน้าเข้าหาจุดเปลี่ยนผ่านวิกฤตตลอดเวลา มวลน้ำและแรงดันมหาศาลสะสมขึ้นเรื่อยๆ จนกระทั่งถึงจุดที่ปัจจัยภายนอกเพียงเล็กน้อย ไม่ว่าจะเป็นแผ่นดินไหว หรือก้อนหิมะถล่มลงมาสะกิด เขื่อนธรรมชาติกักเก็บน้ำก็พังทลายลงทันที ในเสี้ยววินาทีนั้น ระบบได้เปลี่ยนสถานะจากการกักเก็บไปสู่การปลดปล่อยฉับพลัน เกิดเป็น output flux ขนาดมหึมา คือ GLOF มวลน้ำซัดทะลักลงตามไหล่เขา กวาดทำลายทุกสิ่งในชุมชนท้ายน้ำ และทิ้งไว้เพียงร่องรอยของระบบธรรมชาติที่สมดุลพังทลายลงเพราะความร้อนที่มนุษย์เป็นผู้จุดขึ้น ปรากฏการณ์นี้เกิดขึ้น บ่อยครั้งและเป็นภัยพิบัติประจำภูมิภาค สำหรับประเทศแถบเทือกเขาหิมาลัย ทั้งเนปาล ทิเบต และภูฏาน เนื่องจากทั้งสามพื้นที่มีธรณีแปรโครงสร้างและลักษณะทางกายภาพเหมือนกัน คือ มีธารน้ำแข็งจำนวนมากที่กำลังหดตัวอย่างรวดเร็ว ที่ทิเบต ซึ่งเป็นพื้นที่ที่มีการบันทึกเหตุการณ์ GLOF มากที่สุดในภูมิภาคนี้ โดยตั้งแต่ปี 1935 เป็นต้นมา เกิดเหตุการณ์ GLOF ขนาดใหญ่ข้ามระดับวิกฤตไปแล้ว มากกว่า 40 ครั้ง เนปาลมีเหตุการุการณ์ใหญ่เคยเกิดขึ้นเมื่อปี 1981 ที่ทะเลสาบ Cirenmaco เกิดน้ำท่วมฉับพลันทะลักจากฝั่งทิเบตข้ามชายแดนเข้ามายังเนปาล ทำลายล้างโรงไฟฟ้าพาวเวอร์สเตชัน สะพาน และคร่าชีวิตผู้คนไปมากกว่า 200 คน ในเนปาล เคยเกิดเหตุการณ์ GLOF ใหญ่ที่บันทึกไว้ชัดเจนอย่างน้อย 26 ครั้ง โดยในจำนวนนี้ มี 10 ครั้งที่จุดกำเนิดน้ำท่วมเกิดในฝั่งทิเบต แล้วทะลักข้ามแดนลงมาทำลายล้างชุมชนในเนปาล เหตุการุการณ์ใหญ่เคยเกิดขึ้นเมื่อปี 1994 เมื่อเขื่อนธรรมชาติของทะเลสาบ Luggye Tsho พังทลาย ปล่อยมวลน้ำท่วมเมืองโบราณ Punakha ส่งผลให้มีผู้เสียชีวิต 21 คน และทำลายพื้นที่เกษตรกรรมและโบราณสถานสำคัญเสียหายหนัก ขณะที่ที่ภูฏาน มีการบันทึกเหตุการณ์ GLOF ครั้งใหญ่มาแล้วอย่างน้อย 21 ครั้ง โดยมีทะเลสาบธารน้ำแข็งที่ถูกจัดให้อยู่ในกลุ่มอันตรายสูง ณ ปัจจุบันอีกอย่างน้อย 17 แห่ง เหตุการณ์สำคัญที่เคยเกิดขึ้นจริงเมื่อปี2016สิบปีที่แล้วคือ 2016 ที่บริเวณทะเลสาบ Bhote Koshi / Gonggha เกิด GLOF จากฝั่งทิเบต แล้วไหลทะลักเข้าเนปาล ทำลายถนนสายหลักที่เชื่อมการค้าเนปาล-จีน และสะพานสำคัญหลายแห่ง สำหรับเหตุการณ์ภัยพิบัติน้ำท่วมฉับพลันที่เพิ่งเกิดขึ้นในเนปาลเมื่อวานนี้ คือ อุทาหรณ์ที่สะท้อนถึงวิกฤตความสุ่มเสี่ยงของภูมิภาคหิมาลัยได้อย่างชัดเจนที่สุด เมื่อมวลน้ำมหาศาลจาก GLOF ไหลซัดทำลายล้างโครงสร้างพื้นฐานและชุมชนท้ายน้ำในพริบตา ปรากฏการณ์นี้ไม่ใช่เพียงอุบัติเหตุทางธรรมชาติที่เกิดขึ้นเป็นครั้งคราว หากแต่เป็นผลลัพธ์เชิงระบบของการเปลี่ยนแปลงสภาพภูมิอากาศโลก ที่เร่งให้อ่างกักเก็บน้ำแข็งบนยอดเขาสูงละลายและเกินจุดรับไหว การเข้าใจภัยพิบัติดังกล่าวผ่านกรอบความคิดเชิงระบบจึงเป็นหัวใจสำคัญในการมองเห็นความเชื่อมโยงระหว่างการเพิ่มขึ้นของอุณหภูมิโลก การเปลี่ยนรูปของธารน้ำแข็ง และความมั่นคงของมนุษย์ในพื้นที่เปราะบางข้ามพรมแดน ในพื้นที่ทั้งสามที่กล่าวถึงนี้ มีต้นเหตุที่ทำให้เกิด GLOF หลักๆ เพียงแค่สองอย่าง ประการแรกเลย คือ อัตราการละลายสูงที่สุดในโลก ทั้งนี้เทือกเขาหิมาลัยได้ชื่อว่าเป็นขั้วโลกที่สาม (The Third Pole) ซึ่งกำลังอุ่นขึ้นเรื่อยๆ ในอัตราที่เร็วกว่าค่าเฉลี่ยโลก ทำให้เกิดทะเลสาบธารน้ำแข็งใหม่ๆ และขยายขนาดขึ้นอย่างรวดเร็ว อีกประการหนึ่ง ลักษณะธรณีวิทยาที่ไม่มั่นคง ทำให้เขื่อนธรรมชาติที่กั้นน้ำไว้ ที่สร้างขึ้นมาด้วยกองหินและดินที่ธารน้ำแข็งพัดมา (moraine) ซึ่งไม่มีความคงทนแข็งแรง เมื่อเจอกับแผ่นดินไหวบ่อยครั้ง หรือก้อนหิมะ/หินถล่ม (avalanche) ตกลงไปกระแทก น้ำก็พร้อมจะทะลักออกมาทันที เมื่อพิจารณาปัจจัยทั้งสองประการผ่านกรอบความคิดเชิงระบบ ทั้งจากการละลายอย่างรวดเร็วของขั้วโลกที่สาม และความยืดหยุ่นต่ำของคันดินหินธรรมชาติ จะพบว่า ปัญหาไม่ได้เกิดจากปัจจัยใดปัจจัยหนึ่งโดยลำพัง หากแต่เกิดจากปฏิสัมพันธ์เชิงบวกเร่งความเร็ว (reinforcing feedback loop) และการข้ามจุดวิกฤตของระบบธรรมชาติ การเร่งตัวของอุณหภูมิโลกทำให้อัตราการไหลเข้าของน้ำ (input flux) สูงเกินกว่าความจุเชิงโครงสร้างของคันกั้นธรรมชาติที่ไม่มั่นคง เมื่อคันกั้นที่เปราะบางนี้ถูกกระตุ้นด้วยสิ่งรบกวนภายนอกอย่างแผ่นดินไหวหรือหิมะถล่ม ระบบจึงเกิดการพังทลายเชิงโครงสร้าง และปลดปล่อยพลังงานออกมาอย่างฉับพลัน (output flux) การแก้ปัญหานี้จึงไม่สามารถใช้เพียงวิธีทางวิศวกรรมเฉพาะจุด แต่ต้องแก้ที่ จุดคานงัดเชิงระบบ (leverage points) เพื่อลดความเสี่ยงทั้งในระดับโลกและระดับท้องถิ่น ทั้งหลายทั้งปวง มีข้อเสนอเป็นทางออกเชิงระบบ (Systemic Interventions & Leverage Points) เพื่อให้เกิดความปลอดภัยในการอยยู่ร่วมกับธรรมชาติ และมีความสามารถในการรับมือกับสิ่งที่เกิดขึ้นอีกในอนาคต ดังนี้ 1. การบริหารจัดการอ่างกักเก็บและกระแสการไหล (managing reservoirs & fluxes) ทางออกทางกายภาพ (engineering & nature-based solutions): การลดระดับน้ำในทะเลสาบเสี่ยงสูงเชิงรุก (siphon systems & controlled drainage) เพื่อควบคุมปริมาณการกักเก็บ (reservoir volume) ไม่ให้ชนขีดจำกัดวิกฤต พร้อมปรับปรุงคันกั้นธรรมชาติด้วยโครงสร้างวิศวกรรมผสมผสาน เพื่อเพิ่มความต้านทานต่อแรงกระแทกจากหิมะถล่ม 2. การตัดลูปป้อนกลับรุนแรงและรับมือสิ่งรบกวน (disrupting feedback & disturbance response) ระบบเตือนภัยล่วงหน้าที่เชื่อมโยงกัน (integrated early warning systems): ติดตั้งเซนเซอร์ตรวจจับแรงสั่นสะเทือน การเคลื่อนตัวของหิมะ และระดับน้ำ real-time เพื่อเปลี่ยนสัญญาณการกระตุ้นภายนอกให้กลายเป็นระบบเตือนภัยทันที ซึ่งช่วยซื้อเวลาให้ชุมชนท้ายน้ำอพยพได้ทัน 3. การกำกับดูแลข้ามพรมแดน (transboundary systemic governance) ความร่วมมือระดับภูมิภาค (regional data sharing): ภัยพิบัติจาก GLOF ไม่เคยหยุดที่พรมแดนการเมือง ประเทศในแถบหิมาลัย ต้องสร้างข้อตกลงแลกเปลี่ยนข้อมูลดาวเทียมและสภาพอากาศแบบเรียลไทม์ เพื่อวางแผนรับมือมวลน้ำหลากข้ามชายแดนร่วมกัน 4. การปรับโครงสร้างระดับรากเหง้า (structural leverage point) การกู้วิกฤตสภาพภูมิอากาศโลก (global climate mitigation): ในระดับโครงสร้างใหญ่ที่สุด การชะลออัตราการ warming ของขั้วโลกที่สามจำต้องอาศัยการลดการปล่อยก๊าซเรือนกระจกในระดับโลก เพราะหากไม่สามารถลดอุณหภูมิโลกได้ อัตราการละลายของธารน้ำแข็งจะยังคงสร้างทะเลสาบแห่งใหม่ขึ้นมารวดเร็วกว่าที่มนุษย์จะตามไปจัดการได้ทั้งหมด ที่มา - วิเคราะห์ด้วยกรอบ Cyclic & System Thinking จากบทความ Robin George Andrews (2019) Flash floods threaten more people than thought, as ice melts. National Geographic., Published April 18. ใน https://www.nationalgeographic.com/.../himalayan-glacial...

เชื้อเพลิงฟอสซิล (Fossil Fuels)

พัฒนา ราชวงศ์ อาศรมภูมิวิทยาศาสตร์

สาขาวิชาภูมิศาสตร์ มหาวิทยาลัยนเรศวร

การนำเชื้อเพลิงฟอสซิล (Fossil Fuels) ขึ้นมาใช้เป็นการแทรกแซงวัฏจักรคาร์บอนของมนุษย์ที่สร้างผลกระทบต่อ Earth System Dynamics ในระดับโครงสร้างลึก เนื่องจากเป็นการเปลี่ยนสถานะของสสารข้ามมิติเวลา (Temporal Scales) และเปลี่ยนสมดุลของระบบโลกอย่างรวดเร็ว

พิจารณาเรื่องนี้ภายใต้กรอบ cyclic and system thinking จะทำให้้ห๋นรายละเอียดเชิงลึกขององค์ประกอบ 3 ส่วน ดังนี้

1. แหล่งกักเก็บ (Reservoirs)

เมื่อพิจารณาสภาวะดั้งเดิมก่อนยุคอุตสาหกรรม คาร์บอนในโลกถูกแบ่งออกเป็น 2 ระบบหลักตามระยะเวลากักเก็บ (Residence Time) คือ

แหล่งกักเก็บฟอสซิลในชั้นหิน (Geological / Lithospheric Fossil Reservoir) เป็นคาร์บอนอินทรีย์ที่ถูกทับถม ควบแน่น และกักเก็บไว้ใต้ความร้อนและความดันสูงนับหลายร้อยล้านปี (ยุค Carboniferous ถึง Mesozoic) ได้แก่ ถ่านหิน น้ำมันดิบ และก๊าซธรรมชาติ มี Residence Time > 100,000,000 ปี (อยู่ในระบบ Slow Carbon Cycle) โดยมีสถานะดั้งเดิมเป็นระบบปิด (Locked Pool) ที่แทบไม่มีการแลกเปลี่ยนมวลสารกับผิวโลก ยกเว้นการรั่วซึมตามธรรมชาติหรือภูเขาไฟระเบิดในปริมาณที่น้อยมาก

แหล่งกักเก็บหมุนเวียนเร็วหน้าผิวโลก (Active Surface Carbon Reservoirs) ได้แก่ บรรยากาศ (Atmosphere) มหาสมุทร (Hydrosphere) ชีวภาค (Biosphere) และชั้นดิน (Pedosphere) มี Residence Time ตั้งแต่ไม่กี่วัน (การหายใจ) ถึงหลักสิบ/หลักร้อยปี (การเติบโตของป่าไม้) ซึ่งเรียกว่า Fast Carbon Cycle

การขุดเจาะฟอสซิล (Oil & Gas Extraction / Coal Mining) ทำหน้าที่เป็น "Short-circuit Mechanism" หรือกลไกลัดวงจรที่ดึงคาร์บอนจาก Slow Carbon Cycle เข้าสู่ Fast Carbon Cycle ด้วยอัตราเร็วที่มากกว่ากระบวนการคืนกลับตามธรรมชาติถึง 100,000 เท่า ส่งผลให้ปริมาณคาร์บอนรวมในระบบหมุนเวียนหน้าผิวโลกอิ่มตัวเกินขีดความสามารถในการรองรับ (Carrying Capacity)

2. การปฏิสัมพันธ์ข้ามวัฏจักร

การเผาไหม้เชื้อเพลิงฟอสซิล (C_fossil + O2 —> CO2 + Energy) ปล่อยก๊าซคาร์บอนไดออกไซด์มหาศาล ซึ่งแพร่กระจายและทำปฏิกิริยาข้ามวัฏจักรต่างๆ ดังนี้

A. ภาวะมหาสมุทรเป็นกรด 

Geosphere —> Atmosphere —> Hydrosphere

เมื่อ CO2 ในบรรยากาศเพิ่มขึ้น มหาสมุทรจะดูดซับ CO2 ส่วนเกินเพื่อรักษาดุลยภาพความดันก๊าซ ส่งผลให้เกิดปฏิกิริยาลูกโซ่ทางเคมี

1. เกิดกรดคาร์บอนิก (H2CO3) ทำให้ pH ของน้ำทะเลลดลง (เกิด Ocean Acidification)

2. ไฮโดรเจนไอออน (H+) ที่เพิ่มขึ้น เข้าไปจับกับ คาร์บอเนตไอออน (CO3^2-) ทำให้ปริมาณคาร์บอเนตที่สิ่งมีชีวิตใช้สร้างเปลือก/ปะการังลดลง

3. กระบวนการนี้ย้อนกลับไปทำลายการตกตะกอนของหินปูน (Carbonate Sedimentation) ใน Rock Cycle ของท้องทะเล

B. วัฏจักรน้ำและพลังงานบรรยากาศ

Atmosphere —> Water Cycle

CO2 จากฟอสซิลกักเก็บความร้อนในบรรยากาศ ทำให้อุณหภูมิเฉลี่ยของโลกสูงขึ้น ซึ่งตามสมการ Clausius-Clapeyron ทุกๆ อุณหภูมิที่เพิ่มขึ้น 1 เซลเซียส บรรยากาศจะกักเก็บไอน้ำได้เพิ่มขึ้นประมาณ 7% ส่งผลให้อัตราการระเหยของน้ำ (Evaporation) จากมหาสมุทรและผืนดินเร่งตัวขึ้น นำไปสู่การเกิดคลื่นความร้อน ภัยแล้งรุนแรงในบางพื้นที่ และเกิดฝนตกหนักสุดขั้ว (Extreme Precipitation) ในอีกพื้นที่หนึ่ง

C. ชีวภาคและชั้นดิน

Atmosphere —> Biosphere & Pedosphere

สภาพอากาศที่อุ่นขึ้นและแปรปรวน เร่งอัตราการหายใจของจุลินทรีย์ในดิน (Soil Microbial Respiration) จุลินทรีย์จะย่อยสลายสารอินทรีย์วัตถุในดิน (Soil Organic Matter) ได้เร็วขึ้น เปลี่ยนคาร์บอนในดินให้กลายเป็น CO2 ปล่อยกลับสู่บรรยากาศ ซ้ำเติมปริมาณคาร์บอนในอากาศให้สูงขึ้นไปอีก

3. วงจรป้อนกลับเชิงระบบ

การนำเชื้อเพลิงฟอสซิลมาใช้กระตุ้นให้เกิดทั้งแรงเสริมแรงปั่นป่วน (Positive Feedback Loops) และแรงต้านการเปลี่ยนแปลง (Negative Feedback Loops) ดังนี้

3.1 วงจรป้อนกลับแบบบวก - ขยายผลกระทบจนระบบเสียเสถียร (Positive Feedback Loops)

1. วงจรการลดความสามารถในการละลายของมหาสมุทร (Ocean Solubility Feedback)

2. วงจรการละลายของดินแข็งอาร์กติก (Permafrost-Methane / Carbon Feedback)

CO2 จากการเผาฟอสซิลทำให้เขตร้อนและขั้วโลกอุ่นขึ้น —> ชั้นดินเฮือกแข็งคงตัว (Permafrost) ละลาย จุลินทรีย์ตื่นขึ้นมาย่อยสลายซากสิ่งมีชีวิตโบราณ ปลดปล่อยก๊าซมีเทน (CH4) และ CO2 ปริมาณมหาศาลเข้าสู่บรรยากาศ

3. วงจรไฟป่าและการเสื่อมโทรมของผืนป่า (Wildfire & Biomass Decay Feedback)

โลกร้อนขึ้น —> อากาศแห้งแล้ง เกิดไฟป่าบ่อยครั้งและรุนแรงขึ้น ป่าไม้ซึ่งเป็นแหล่งกักเก็บคาร์บอนหลัก (Carbon Sink) ถูกเผาทำลาย เปลี่ยนสถานะกลายเป็นแหล่งปลดปล่อยคาร์บอน (Carbon Source) สู่บรรยากาศ

3.2 วงจรป้อนกลับแบบลบ - กลไกหน่วง/สร้างดุลยภาพตามธรรมชาติ (Negative Feedback Loops)

1. ปรากฏการณ์เร่งการโตของพืช (CO2 Fertilization Effect)

ความเข้มข้นของ CO2 ในอากาศที่สูงขึ้น กระตุ้นให้พืชบนบกบางชนิดเพิ่มอัตราการสังเคราะห์ด้วยแสงและการปิดปากใบเพื่อลดการสูญเสียน้ำ พืชจึงช่วยดูดซับ CO2 กลับมาเก็บไว้ในมวลชีวภาพมากขึ้นชั่วครา กลไกนี้ถูกจำกัดด้วยปริมาณธาตุอาหารในดิน (เช่น ไนโตรเจนและฟอสฟอรัส) และความแห้งแล้งสะสม

2. การกัดกร่อนหินซิลิเกตในสเกลยาวนาน (Enhanced Silicate Weathering Feedback)

CO2 สูงขึ้น —> อุณหภูมิและปริมาณฝนเพิ่มขึ้น —> เกิดกรดคาร์บอนิกชะล้างหินซิลิเกตบนเปลือกโลก ไอออนแคลเซียมและไบคาร์บอเนตจะถูกพัดลงทะเล และตกตะกอนเป็นหินปูน ตรึงคาร์บอนกลับสู่ธรณีภาค กระบวนการนี้ใช้เวลา หลักหมื่นถึงหลักแสนปี จึงไม่สามารถยับยั้งวิกฤตการปล่อยฟอสซิลของมนุษย์ที่เกิดขึ้นในระดับรวดเร็วเป็นหลักสิบปีได้ทัน

แผนที่ธรณีวิทยา

พัฒนา ราชวงศ์ อาศรมภูมิวิทยาศาสตร์

สาขาวิชาภูมิศาสตร์ มหาวิทยาลัยนเรศวร

 

การทำแผนที่ธรณีวิทยา (Geological Mapping) คือกระบวนการสำรวจ บันทึก และถ่ายทอดข้อมูลลักษณะทางกายภาพของเปลือกโลก ลงบนแผนผังย่อส่วนที่มีระบบพิกัดอ้างอิงทางภูมิศาสตร์ หลักการสำคัญไม่ได้เป็นเพียงการบันทึกตำแหน่งหินที่พบบนผิวดิน แต่เป็นการรวบรวมหลักฐานเพื่อตีความโครงสร้างสามมิติใต้ผิวโลกและวิวัฒนาการทางเวลากว่าหลายร้อยล้านปี กระบวนการทำแผนที่ธรณีวิทยาอาศัยหลักการและขั้นตอนพื้นฐาน ดังนี้

 

·  การสำรวจโผล่หิน (Outcrop Mapping) เป็นการเดินสำรวจภาคสนามเพื่อบันทึกพิกัด ชนิดหิน (Lithology) เนื้อหิน และองค์ประกอบทางแร่ บริเวณที่หินแข็งโผล่พ้นดินตามลำน้ำ หน้าผา หรือการขุดเจาะ

·  การวัดทิศทางเชิงโครงสร้าง (Structural Measurements) ด้วยการใช้เข็มทิศธรณี (Brunton Compass) หรืออุปกรณ์ดิจิทัล วัดค่าทิศทางแนวรอยชั้น (Strike) และมุมเอียงเท (Dip) ของชั้นหิน รวมถึงแนวรอยแตก (Joints) และรอยเลื่อน (Faults)

·  การจัดลำดับชั้นหิน (Stratigraphic Succession) ตามการใช้กฎการสืบชั้น (Law of Superposition) และซากดึกดำบรรพ์ดรรชนี (Index Fossils) เพื่อจัดลำดับอายุของหินจากแก่ที่สุดไปหาอ่อนที่สุด

·  การเปลี่ยนผ่านสู่เทคโนโลยีดิจิทัล (Digital Transition) ปัจจุบันการทำแผนที่ภาคสนามพัฒนาไปสู่การใช้แท็บเล็ตภาคสนาม (Digital Field Mapping) ภาพถ่ายโดรนความละเอียดสูง (UAV Photogrammetry) และการสแกนด้วยเลเซอร์ (LiDAR) เพื่อสกัดข้อมูลเชิงโครงสร้างในพื้นที่เข้าถึงยากได้อย่างแม่นยำ

 

แผนที่ธรณีวิทยา

 

แผนที่ธรณีวิทยา คือ แผนที่แสดงการกระจายตัวของชนิดหิน หมวดหิน และวัตถุปกคลุมผิวดิน พร้อมทั้งแสดงเส้นโครงสร้างทางธรณีวิทยาที่ตัดผ่านพื้นที่นั้นๆ

 

องค์ประกอบมาตรฐานบนแผ่นระวาง

 

1.   ส่วนรายละเอียดเชิงพื้นที่ (Map Face) แสดงขอบเขตการปรากฏของกลุ่มหิน/หมวดหินด้วยรหัสสีสากล พร้อมสัญลักษณ์กำกับโครงสร้าง เช่น เส้นทึบแสดงรอยเลื่อนที่ชัดเจน เส้นประแสดงรอยเลื่อนคาดคะเน และสัญลักษณ์แกนรอยคดโค้ง (Anticline/Syncline)

2.   แท่งลำดับชั้นหินและคำอธิบายสัญลักษณ์ (Stratigraphic Column & Legend) จัดเรียงหมวดหินตามอายุกาลวิทยา โดยให้หินที่มีอายุแก่ที่สุดอยู่ด้านล่างสุด และหินอ่อนที่สุดอยู่ด้านบนสุด พร้อมระบุชนิดหินหลักและสภาพแวดล้อมการสะสมตัว

3.   รูปตัดขวางทางธรณีวิทยา (Geological Cross-section) ภาพจำลองชั้นหินใต้ผิวดินตามแนวตัดเฉพาะ (Cross-section Line) ช่วยให้ผู้ใช้มองเห็นความหนา มุมเอียงเท และรูปแบบการถูกตัดขาดของชั้นหินใต้ดิน

4.   แผนผังตำแหน่งระวาง (Index Map) และพิกัด (Coordination) เป็นการระบุระบบพิกัด (เช่น UTM) ทิศเหนือจริง/ทิศเหนือแม่เหล็ก และมาตราส่วนแผนที่

 

กรณีศึกษา: แผนที่ธรณีวิทยาระวาง NE 47-14 (เมาะลำเลิง / Moulmein) ระวาง NE 47-14 หรือ ระวางเมาะลำเลิง เป็นหนึ่งในระวางที่มีความซับซ้อนทางธรณีวิทยาและโครงสร้างเทกโทนิกส์สูงที่สุดในประเทศไทย ครอบคลุมพื้นที่ชายแดนไทย-เมียนมา บริเวณจังหวัดตาก (อำเภอแม่สอด พบพระ แม่ระมาด) อำเภออมก๋อย จังหวัดเชียงใหม่ และรอยต่อลุ่มน้ำปิง

 

แผนที่จำลองโครงสร้างทางธรณีวิทยาของระวาง NE 47-14 แสดงแนว กลุ่มรอยเลื่อนเมย / แม่สอด (Moei / Mae Ping Fault Zone) ซึ่งเป็นรอยเลื่อนตามแนวระดับขนาดใหญ่ (Strike-slip Fault) วางตัวในทิศNW-SE โครงสร้างนี้ตัดผ่านมวลหินแปรเกรดสูงยุคก่อนแคมเบรียน (Precambrian) และกลุ่มหินตะกอนยุคพาลีโอโซอิก-เมโซโซอิก เกิดเป็นแอ่งเทอร์เชียรีแม่สอด (Mae Sot Tertiary Basin) ซึ่งเป็นแหล่งสะสมตัวของหินดินดานน้ำมัน (Oil Shale) และหินตะกอนยุคใหม่

 

ลักษณะทางธรณีวิทยาและทรัพยากรในระวาง NE 47-14 แสดงกลุ่มหินตะกอนและหินแปร พบหินปูน หินดินดาน และหินทราย ยุคคาร์บอนิเฟอรัส–เพอร์เมียน รวมถึงหินแปรเกรดสูง (Gneiss/Schist) บริเวณเทือกเขาถนนธงชัย มวลหินอัคนีแทรกซอน แนวหินแกรนิต (Granite Plutons) วางตัวในทิศเหนือ-ใต้ ซึ่งสัมพันธ์กับการเกิดแหล่งแร่โลหะ (เช่น สังกะสี ดีบุก วุลแฟรม) และความสำคัญเชิงภัยพิบัติ การเป็นที่ตั้งของกลุ่มรอยเลื่อนมีชีวิตเมย ทำให้ระวางนี้เป็นพื้นที่เฝ้าระวังภัยพิบัติแผ่นดินไหวและดินถล่มระดับภูมิภาค

 

สำหรับการสังเคราะห์วรรณกรรมวิชาการระดับสากลจากวารสารชั้นนำ เช่น Journal of Structural Geology, Geomorphology, Ore Geology Reviews, และ Earth-Science Reviews สามารถแบ่งประเภทแผนที่ตามแนวทางการวิจัยของนักธรณีวิทยาได้เป็น 4 กลุ่มหลัก

 

แผนที่ธรณีวิทยาเชิงโครงสร้างและดิจิทัล

 

งานวิจัยของ Cawood et al. (2017) และ Tavani et al. (2020) เน้นการวิเคราะห์เรขาคณิตของมวลหิน แรงเทกโทนิกส์ และสถาปัตยกรรมรอยเลื่อน โดยได้พัฒนากระบวนการทำแผนที่จากแบบจำลองหินดิจิทัล (Digital Outcrop Models: DOMs) โดยใช้โดรน เพื่อเก็บข้อมูล Strike/Dip อัตโนมัติ ส่วน Sylvester (1988) กำหนดกรอบการทำแผนที่กลุ่มรอยเลื่อนตามแนวระดับ (Strike-slip Fault Zones)

 

แผนที่ธรณีสัณฐานวิทยาและภูมิประเทศรายละเอียดสูง

 

งานวิจัยของ Tarolli (2014) เป็นการศึกษารูปร่างภูมิประเทศ กระบวนการเปลี่ยนแปลงผิวดิน และระบบทางน้ำ สรุปการใช้ข้อมูล High-Resolution Topography (LiDAR) สกัดแผนที่ธรณีสัณฐานระดับเล็ก Brierley & Fryirs (2000) นำเสนอกรอบ River Styles Framework เพื่อจำแนกประเภทแม่น้ำระดับลุ่มน้ำจากโครงสร้างธรณีวิทยา และ Bishop et al. (2012) นำเสนอการทำแผนที่สี่มิติ (4D Spatiotemporal Mapping)

 

แผนที่สำรวจแร่ด้วยการรีโมตเซนซิง

 

งานวิจัยของ van der Meer et al. (2012) และ Cudahy & Ramanaidou (1997) ใช้ข้อมูลภาพถ่ายดาวเทียม Multispectral/Hyperspectral (เช่น ASTER) เน้นการใช้คลื่นแม่เหล็กไฟฟ้าตรวจจับโซนการแปรสภาพของหินเพื่อตีกรอบพื้นที่เป้าหมาย ทำแผนที่กลุ่มแร่แปรสภาพทางน้ำร้อน (Hydrothermal Alteration Minerals) เช่น Alunite และ White Micas

 

แผนที่และแบบจำลองธรณีวิทยา 3 มิติ/4 มิติ

 

งานวิจัยของ Wu et al. (2005) เน้นการเปลี่ยนผ่านจากการนำเสนอภาพ 2 มิติไปสู่ปริมาตรมวลหิน 3 มิติใต้ผิวดิน สรุปเป็นโครงสร้างข้อมูลเชิงพื้นที่แบบ Voxel, Wireframe และ B-Rep เพื่อสร้าง 3D Geological Model ส่วน Scotese (2004) สังเคราะห์ แผนที่ธรณีวิทยาบรรพกาล (Paleogeographic Maps) แสดงการเคลื่อนตัวของแผ่นทวีปตั้งแต่มหาทวีปพันเจีย (Pangea)

 

แผนที่เฉพาะเรื่องทางธรณีวิทยา

 

แผนที่เฉพาะเรื่อง (Thematic Geological Maps) คือ การนำข้อมูลธรณีวิทยาพื้นฐานมาซ้อนทับและประมวลผลร่วมกับปัจจัยทางกายภาพอื่นๆ บนระบบ GIS เพื่อตอบโจทย์การใช้งานเฉพาะด้าน

 

ประเภทแผนที่เฉพาะเรื่อง	ข้อมูลนำเข้าหลัก	การนำไปประยุกต์ใช้	งานวิจัยอ้างอิง
แผนที่รอยเลื่อนมีชีวิต (Active Fault Map)	แนวรอยเลื่อน, ร่องสำรวจ (Trenching), อัตราการเลื่อนตัว (Slip Rate)	ประเมินพิบัติภัยแผ่นดินไหว และกำหนดเขตพื้นที่เสี่ยง	Yeats & Prentice (1996); El Hamdouni et al. (2008)
แผนที่เสี่ยงภัยดินถล่ม (Landslide Susceptibility Map)	ชนิดหิน, รอยเลื่อน, ความชัน, ปริมาณน้ำฝน, การใช้ประโยชน์ที่ดิน	วางผังเมือง, จัดทำระบบเตือนภัยพิบัติล่วงหน้า	Guzzetti et al. (1999); Reichenbach et al. (2018)
แผนที่ศักยภาพแหล่งแร่ (Mineral Prospectivity Map: MPM)	แผนที่หิน, รอยเลื่อน, ข้อมูลเคมีธรณี, โซน Alteration	ตีกรอบพื้นที่สำรวจแร่เศรษฐกิจด้วย Machine Learning	Hronsky & Groves (2008); Zuo & Xiong (2018)
แผนที่ธรณีวิศวกรรม (Engineering/Geotechnical Map)	กำลังรับน้ำหนักของหิน/ดิน, คุณสมบัติทางกลศาสตร์, ระดับน้ำบาดาล	วางแผนฐานรากอาคาร อุโมงค์ สะพาน และเขื่อน	-

 

แผนธรณีวิทยาประเทศไทย

 

การทำแผนที่ธรณีวิทยาของประเทศไทยอยู่ในความรับผิดชอบหลักของ กรมทรัพยากรธรณี กระทรวงทรัพยากรธรรมชาติและสิ่งแวดล้อม

 

แผนที่ธรณีวิทยา มาตราส่วน 1:250,000

 

กรมทรัพยากรธรณีได้จัดทำแผนที่ธรณีวิทยาครอบคลุมพื้นที่ทั่วประเทศไทยรวมทั้งสิ้น 47 ระวาง โดยใช้อ้างอิงระบบพิกัด UTM

 

• ระดับรายละเอียด: จัดเป็นแผนที่ระดับภูมิภาค (Regional Scale) หรือการสำรวจสังเขป เน้นการจัดกลุ่มหินออกเป็น หมวดหิน (Formation) หรือ กลุ่มหิน (Group)
• การนำไปใช้: เหมาะสำหรับการวางแผนเชิงนโยบายระดับประเทศ การบริหารจัดการลุ่มน้ำ และการประเมินศักยภาพทรัพยากรธรรมชาติ

 

กรณีศึกษา ระวาง NE 47-14 (เมาะลำเลิง / Moulmein)

 

เป็นระวางทางตะวันตกเฉียงเหนือ ครอบคลุมพื้นที่อำเภอแม่สอด พบพระ แม่ระมาด (จังหวัดตาก) อำเภออมก๋อย (จังหวัดเชียงใหม่) และรอยต่อประเทศเมียนมา

• ลักษณะเด่น: มี แนวรอยเลื่อนเมย / แม่สอด (Moei / Mae Ping Fault Zone) ซึ่งเป็นรอยเลื่อนตามแนวระดับขนาดใหญ่ ตัดผ่านมวลหินแปรเกรดสูงยุคก่อนแคมเบรียน และทำให้เกิด แอ่งเทอร์เชียรีแม่สอด ซึ่งเป็นแหล่งสะสมตัวของหินดินดานน้ำมัน (Oil Shale)

 

ช่องทางการเข้าถึงข้อมูลแผนที่

 

• ไฟล์ระวาง PDF: ดาวน์โหลดผ่านระบบบริการแผนที่ของกรมทรัพยากรธรณี
• ระบบ GIS & Web Mapping: เข้าใช้งานผ่าน DMR GIS Portal และคลังข้อมูลเปิดภาครัฐ (Data.go.th / Data.dmr.go.th) ซึ่งให้บริการในรูปแบบ Shapefile และ ArcGIS REST Service API

 

แผนที่ธรณีวิทยาสามมิติ

แผนที่และแบบจำลองธรณีวิทยา 3 มิติ (3-D Geological Maps & Models) เป็นการพัฒนาก้าวสำคัญจากการนำเสนอภาพ 2 มิติบนระวางกระดาษ ไปสู่การสร้าง "แบบจำลองปริมาตรมวลหินใต้ผิวดิน (Volume Modeling)" ในสามมิติจริง (X, Y, Z) ซึ่งช่วยให้เห็นโครงสร้าง ความหนา และการต่อเนื่องของชั้นหินใต้พื้นผิวโลกได้อย่างแม่นยำ

 

หลักการของ 3-D Geological Maps

 

หลักการสำคัญของการทำแผนที่ธรณีวิทยา 3 มิติ คือการนำข้อมูลจากต่างที่มาและต่างมิติ มาเชื่อมโยงและคำนวณผ่านอัลกอริทึมเชิงพื้นที่ (Spatial Interpolation & Geostatistics) โดยมีหลักการพื้นฐานดังนี้

·       Boundary Representation (B-Rep) & Surface Modeling: การสร้างพื้นผิวรอยต่อระหว่างหมวดหิน (Geological Interfaces) และระนาบรอยเลื่อน (Fault Planes) ให้เป็นระนาบ 3 มิติในรูปแบบตาข่ายสามเหลี่ยม (TIN - Triangulated Irregular Network) หรือพื้นผิวโค้งระดับสูง (NURBS)

·       Volumetric Modeling (Voxel / Mesh): การแบ่งปริมาตรมวลหินใต้ดินออกเป็นลูกบาศก์เล็กๆ (Voxels) หรือตาข่ายองค์ประกอบย่อย (3D Unstructured Grids/Mesh) เพื่อให้สามารถบรรจุคุณสมบัติทางกายภาพของหิน เช่น ความพรุน (Porosity), ความหนาแน่น (Density), หรือความเร็วคลื่นไหวสะเทือน (Seismic Velocity) ลงไปในทุกๆ จุดของปริมาตรได้

·       Geological Logic Rules: การกำหนดกฎเกณฑ์ทางธรณีวิทยาบังคับกระบวนการประมาณค่า เช่น กฎการสืบชั้น (Stratigraphic Age Rules), ความสัมพันธ์การตัดขวางของรอยเลื่อน (Fault Cross-cutting Relationships), และมุมเอียงเท (Strike/Dip) เพื่อไม่ให้โมเดลขัดแย้งกับความเป็นจริงทางธรรมชาติ

 

ขั้นตอนและกระบวนการจัดทำ (Workflow)

 

1. การรวบรวมและเตรียมข้อมูล (Data Integration)

 

รวบรวมข้อมูลหลากมิติเข้าสู่ระบบ GIS/Modelingข้อมูลผิวดิน: แบบจำลองความสูงเชิงเลข (DEM/LiDAR), แผนที่ธรณีวิทยา 2D, แนวรอยเลื่อน, และจุดวัด Strike/Dipข้อมูลใต้ผิวดิน: ข้อมูลหลุมเจาะ (Borehole Log), ภาพถ่ายคลื่นไหวสะเทือน (Seismic Reflection 2D/3D), และข้อมูลธรณีฟิสิกส์ (Gravity/Magnetic)

 

2. การกำหนดกรอบเรขาคณิตรอยเลื่อน (Fault Modeling)

 

สร้างโครงข่ายรอยเลื่อน 3 มิติเป็นอันดับแรกกำหนดระนาบรอยเลื่อน (Fault Surfaces) ในระบบ 3 มิติ เนื่องจากรอยเลื่อนเป็นตัวตัดขาดและเลื่อนชั้นหินออกจากกัน การสร้างรอยเลื่อนที่ถูกต้องก่อนจะช่วยกำหนดขอบเขต (Fault Blocks) สำหรับชั้นหินได้อย่างแม่นยำ

 

3. การสร้างพื้นผิวชั้นหิน (Stratigraphic Surface Modeling)

 

เชื่อมโยงขอบเขตชั้นหินระหว่างหลุมเจาะและผิวดินนำจุดตัดชั้นหินจากหลุมเจาะ (Well Tops) และแนวขอบเขตหินบนผิวดินมาคำนวณสร้างพื้นผิวของแต่ละหมวดหิน โดยใช้วิธีการคณิตศาสตร์ เช่น Kriging, Radial Basis Functions (RBF), หรือ Implicit Surface Modeling

 

4. การสร้างปริมาตรและประมาณค่าคุณสมบัติ (Volumetric & Property Modeling)

 

เติมเต็มข้อมูลเชิงคุณสมบัติลงในมวลหินสร้างโครงข่ายปริมาตร 3D Grid และใช้เทคนิคสถิติภูมิศาสตร์ (Geostatistics) เพื่อประมาณค่าคุณสมบัติทางธรณีวิศวกรรม ธรณีเคมี หรือปิโตรฟิสิกส์ลงไปในเนื้อหิน

 

5. การตรวจสอบและปรับแก้แบบจำลอง (Validation & Model Refinement)

 

ตรวจสอบความสมเหตุสมผลทางธรณีวิทยาตัดรูปตัดขวาง (Cross-section) จากโมเดล 3D มาเปรียบเทียบกับภาพจริงทางภาคสนาม หรือทำการคำนวณย้อนกลับทางธรณีฟิสิกส์ (Forward/Inverse Geophysical Modeling) เพื่อปรับแก้โมเดลให้ตรงกับค่าสนามแรงดึงดูด/แม่เหล็กจริง

ประโยชน์ของ 3-D Geological Maps

การเปลี่ยนจากแผนที่ 2D มาเป็นแบบจำลอง 3D เพิ่มประสิทธิภาพอย่างมากในหลายภาคส่วน:

ภาคส่วน / การประยุกต์ใช้	ประโยชน์หลักของแบบจำลอง 3 มิติ
1. การสำรวจแร่และพลังงาน	ตีกรอบแหล่งแร่ สภาพกับเก็บปิโตรเลียม (Hydrocarbon Reservoirs) และแหล่งพลังงานความร้อนใต้พิภพ ช่วยคำนวณปริมาตรสำรอง (Reserve Estimation) และลดความเสี่ยงในการขุดเจาะหลุมสำรวจราคาแพง
2. วิศวกรรมฐานรากและงานอุโมงค์	ช่วยวิเคราะห์แนวขุดเจาะอุโมงค์ (Tunneling) เขื่อน หรือรถไฟใต้ดิน โดยแสดงจุดรอยเลื่อน มวลหินผุ หรือโพรงหินปูน (Karst) ใต้ดินล่วงหน้า เพื่อป้องกันอุบัติเหตุระหว่างก่อสร้าง
3. ทรัพยากรน้ำบาดาล (Hydrogeology)	สร้างแบบจำลองการไหลของน้ำบาดาล (Groundwater Flow Modeling) แสดงขอบเขตชั้นน้ำบาดาล (Aquifers) และทิศทางการแพร่กระจายของสารปนเปื้อนในระดับลึก
4. การประเมินภัยพิบัติทางธรณีวิทยา	วิเคราะห์โครงสร้างรอยเลื่อนมีชีวิตใต้ดินแบบสามมิติ ประเมินระนาบการการพังทลายของชั้นดิน/หินในการเกิดดินถล่มขนาดใหญ่ (Landslide Slip Surfaces)
5. การสื่อสารและการตัดสินใจเชิงนโยบาย	ช่วยให้ผู้เชี่ยวชาญ สถาปนิก วิศวกร ตลอดจนผู้กำหนดนโยบายและประชาชนทั่วไป เข้าใจสภาพโครงสร้างใต้ดินได้ง่ายกว่าการอ่านแผนที่ 2D

สรุป

 

แผนที่ธรณีวิทยาถือเป็นโครงสร้างพื้นฐานทางปัญญา (Intellectual Infrastructure) ที่มีความสำคัญต่อการพัฒนาประเทศในหลากหลายมิติ

 

1.        การบริหารจัดการทรัพยากรธรรมชาติ: เป็นเข็มทิศในการสำรวจแหล่งแร่เศรษฐกิจ แหล่งพลังงาน (ปิโตรเลียม หินดินดานน้ำมัน พลังงานความร้อนใต้พิภพ) และการประเมินศักยภาพแหล่งน้ำบาดาล

2.      การบรรเทาและรับมือภัยพิบัติธรรมชาติ: เป็นข้อมูลฐานในการประเมินพื้นที่เสี่ยงภัยแผ่นดินไหวจากรอยเลื่อนมีชีวิต ดินถล่ม รูยุบในพื้นที่หินปูน และการทรุดตัวของแผ่นดิน

3.      งานวิศวกรรมและโครงสร้างพื้นฐาน: ช่วยในการเลือกตำแหน่งที่ตั้งและวางแผนออกแบบฐานรากสำหรับเขื่อน อุโมงค์ ทางหลวง รถไฟความเร็วสูง และสิ่งก่อสร้างขนาดใหญ่

4.      การวางผังเมืองและการใช้ประโยชน์ที่ดิน: ป้องกันการก่อสร้างทับซ้อนบนแนวรอยเลื่อนเสี่ยงภัย หรือพื้นที่แหล่งแร่สำคัญ และช่วยจัดสรรพื้นที่เกษตรกรรมตามคุณสมบัติทางธรณีวิทยา

5.      การศึกษาวิจัยและสิ่งแวดล้อม: ช่วยในการทำความเข้าใจประวัติศาสตร์วิวัฒนาการของโลก และการประเมินการเปลี่ยนแปลงของสภาพแวดล้อมบรรพกาล

 

บรรณานุกรม

 

กรมทรัพยากรธรณี. (2544). แผนที่ธรณีวิทยาประเทศไทย มาตราส่วน 1:250,000 ระวางเมาะลำเลิง (NE 47-14). กองธรณีวิทยา กรมทรัพยากรธรณี กระทรวงทรัพยากรธรรมชาติและสิ่งแวดล้อม.

กรมทรัพยากรธรณี. (2567). ระบบสารสนเทศภูมิศาสตร์และบริการข้อมูลแผนที่ธรณีวิทยา (DMR GIS Portal). สืบค้นจาก https://gisportal.dmr.go.th/

Bishop, M. P., James, L. A., Shroder, J. F., & Walsh, S. J. (2012). Remote sensing and GIScience in geomorphology: Introduction to the special issue. Geomorphology, 137(1), 1–12.

Brierley, G. J., & Fryirs, K. A. (2000). River Styles, a geomorphic approach to catchment characterization: Implications for river management. Geomorphology, 32(3–4), 361–385.

Cawood, A. J., Bond, C. E., Inkpen, R. J., & Butler, R. W. H. (2017). Digital geological mapping and 3D visualisation: A workflow for structural geology field data. Journal of Structural Geology, 103, 117–129.

Cudahy, T. J., & Ramanaidou, E. R. (1997). Hyperspectral remote sensing of white micas: Theory, methods and applications in mineral exploration and regional geology. Earth-Science Reviews, 42(1–2), 19–35.

El Hamdouni, R., Irigaray, C., Fernández, T., Chacón, J., & Keller, E. A. (2008). Assessment of active tectonics based on geomorphic indices, Southern Spain. Geomorphology, 96(1–2), 150–173.

Guzzetti, F., Carrara, A., Cardinali, M., & Reichenbach, P. (1999). Landslide hazard and risk mapping: A review of standard approaches and current trends. Earth-Science Reviews, 45(3–4), 181–216.

Hronsky, J. M. A., & Groves, D. I. (2008). Science of targeting: Definition, applied methodology and case studies. Ore Geology Reviews, 33(3–4), 479–491.

Reichenbach, P., Rossi, M., Malamud, B. D., Mihir, M., & Guzzetti, F. (2018). A review of statistically-based landslide susceptibility models. Geomorphology, 315, 60–78.

Scotese, C. R. (2004). Pangea to the modern world: Tectonic controls on regional geology, sedimentation, and climate. Earth-Science Reviews, 67(1–2), 33–55.

Sylvester, A. G. (1988). Strike-slip faults. Journal of Structural Geology, 10(6), 665–689.

Tavani, S., Granado, P., Corradetti, A., Señoráns, R., & Muñoz, J. A. (2020). Building virtual outcrops, extracted orientations, and structural mapping with drones. Journal of Structural Geology, 137, 104082.

Tarolli, P. (2014). High-resolution topography for understanding earth surface processes: Research frontiers. Geomorphology, 216, 295–312.

van der Meer, F. D., van der Werff, H. M., van Ruitenbeek, F. J., Hecker, C. A., Bakker, W. H., Noomen, M. F., ... & Woldai, T. (2012). Multi- and hyperspectral geologic remote sensing: A review. Ore Geology Reviews, 44, 112–131.

Wu, Q., Xu, H., & Zou, X. (2005). Three-dimensional geological modelling: Methods, technical challenges and applications. Earth-Science Reviews, 71(1–2), 85–103.

Yeats, R. S., & Prentice, C. S. (1996). Active faulting and paleoseismology: History, progress, and present status. Earth-Science Reviews, 40(1–2), 71–89.

Zuo, R., & Xiong, Y. (2018). Big data analytics of geochemical data for mineral prospectivity mapping. Ore Geology Reviews, 94, 311–320.