Giới thiệu dự án

Biến đổi khí hậu (BĐKH) toàn cầu đang đặt ra những thách thức nghiêm trọng đối với an ninh sinh thái và sinh kế bền vững, đặc biệt tại các hệ sinh thái karst (núi đá vôi) nhạy cảm vùng Đông Bắc Việt Nam. Theo báo cáo AR5 của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) và kịch bản biến đổi khí hậu của Bộ Tài nguyên và Môi trường (BTNMT), Việt Nam là một trong 5 quốc gia chịu tổn thương nặng nề nhất. Tại tỉnh Cao Bằng, nền nhiệt độ đã tăng trung bình 0,05°C - 0,10°C/thập kỷ, cùng với sự phân bố mưa dị thường gây suy thoái nguồn tài nguyên nước và suy giảm nghiêm trọng độ che phủ đa dạng sinh học (ĐDSH).

Thực trạng tại huyện Trùng Khánh – nơi lưu giữ sinh cảnh của loài Vượn Cao Vít (Nomascus nasutus) cực kỳ nguy cấp – cho thấy các đợt hạn hán kéo dài (mực nước sông giảm hơn 30% vào mùa khô) xen kẽ các đợt lũ quét cực đoan (ngập úng hơn 400 ha lúa năm 2009). Tuy nhiên, rào cản lớn nhất đối với công tác bảo tồn và quản lý tài nguyên là nhận thức của cộng đồng địa phương (chiếm trên 95% là đồng bào Tày, Nùng) còn nhiều hạn chế, dẫn đến sự thiếu chủ động trong việc triển khai các giải pháp thích ứng.

Mục tiêu dự án

  1. Khảo sát và định lượng hóa: Phân tích biến động chuỗi số liệu khí tượng thủy văn (nhiệt độ, lượng mưa) giai đoạn 1994–2013 và hiện trạng hai lưu vực sông chính (Quây Sơn, Bắc Vọng).
  2. Đánh giá ĐDSH: Xác lập danh lục phân loại thảm thực vật, cấu trúc phân bố động vật và các nguồn gen cây trồng, vật nuôi bản địa đặc hữu.
  3. Đo lường nhận thức xã hội: Thực hiện điều tra xã hội học định lượng về mức độ hiểu biết, nguyên nhân nhận thức và nhu cầu đào tạo về BĐKH của 90 hộ dân tại 3 xã trọng điểm (Phong Nậm, Ngọc Côn, Ngọc Khê).
  4. Xây dựng khung giải pháp: Đề xuất mô hình thích ứng tổng hợp tích hợp quản lý lưu vực nước và bảo tồn nguồn gen bản địa.

Phạm vi và giới hạn: Nghiên cứu tập trung tại huyện Trùng Khánh (tổng diện tích tự nhiên 46.872 ha), trọng tâm là vùng đệm và vùng lõi (1.600 ha) của Khu bảo tồn loài - sinh cảnh Vượn Cao Vít; số liệu khí tượng được thu thập qua chuỗi trạm quan trắc địa phương trong 20 năm.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Hiện trạng quản lý môi trường nước và bảo tồn ĐDSH tại huyện Trùng Khánh đang đối mặt với sự phân hóa sâu sắc giữa phương thức can thiệp truyền thống và yêu cầu quản trị thích ứng hiện đại.

Tiêu chí Giải pháp quản lý truyền thống Mô hình đánh giá tích hợp nhận thức (Nghiên cứu)
Cơ sở dữ liệu Dữ liệu khí tượng thủy văn đơn lẻ, thiếu tích hợp sinh thái Kết hợp chuỗi quan trắc 20 năm với điều tra xã hội học tại chỗ
Tiếp cận cộng đồng Tuyên truyền một chiều, áp đặt chỉ tiêu hành chính Đo lường chỉ số nhận thức (KAP), phân tầng theo nhu cầu đào tạo
Bảo tồn ĐDSH Khoanh nuôi bảo vệ thụ động Tích hợp bảo tồn nguồn gen bản địa (Hạt dẻ, Lợn đen) với sinh cảnh karst
Quản trị rủi ro Khắc phục sự cố sau thiên tai (lũ lụt, dịch bệnh) Thiết lập kịch bản dự báo, chủ động chuyển đổi cơ cấu mùa vụ

Phân tích yêu cầu hệ thống theo ma trận MoSCoW:

  • Must have: Bộ chỉ số thống kê khí tượng 1994–2013; ma trận đánh giá nhận thức 90 mẫu điều tra; danh lục ĐDSH nguy cấp theo Sách Đỏ.
  • Should have: Khung định lượng tương quan giữa tần suất thời tiết cực đoan và suy thoái nguồn nước mặt.
  • Could have: Mô hình bản đồ hóa GIS cho các tiểu lưu vực thoát nước sông Quây Sơn và Bắc Vọng.
  • Won't have: Mô phỏng số trị động lực học dòng chảy 3D thời gian thực.

Thiết kế hệ thống và công cụ phân tích

Quy trình xử lý dữ liệu nghiên cứu được chuẩn hóa thông qua pipeline phân tích thống kê chuyên dụng:

[Dữ liệu thô: Trạm Khí tượng + 90 Phiếu khảo sát] 
[Data Cleaning: Pandas 2.0.3 / Excel 2013 VBA]
[Phân tích thống kê mô tả & Kiểm định ANOVA: R 4.3.1 / SPSS 22.0]
[Biểu diễn trực quan: Matplotlib 3.7.1 & Seaborn 0.12.2]
[Khung chính sách quản lý & thích ứng địa phương]

Phương pháp nghiên cứu (Methodology)

Nghiên cứu áp dụng phương pháp tiếp cận đa ngành:

  • Thu thập tài liệu thứ cấp: Thu thập chuỗi số liệu khí hậu 1962–2013 từ Sở TN&MT tỉnh Cao Bằng và Báo cáo đánh giá lần thứ 5 của IPCC.
  • Phương pháp điều tra xã hội học: Thiết kế bảng câu hỏi cấu trúc với thang đo Likert, tiến hành chọn mẫu ngẫu nhiên phân tầng tại xã Phong Nậm, Ngọc Côn và Ngọc Khê ($N = 90$ hộ).
  • Phương pháp chuyên gia (Delphi): Tham vấn các nhà khoa học thuộc Viện Sinh thái và Tài nguyên sinh vật, cán bộ Chi cục Bảo vệ Môi trường và Tổ chức FFI (Fauna & Flora International).

Implementation và kết quả

Quy trình xử lý và thuật toán phân tích

Dữ liệu khảo sát xã hội học và chuỗi quan trắc nhiệt độ - lượng mưa được tính toán bằng các thuật toán thống kê mô tả, tính độ lệch chuẩn và hồi quy xu thế thời gian.

import numpy as np
import pandas as pd
from scipy import stats

# Xử lý chuỗi dữ liệu nhiệt độ trung bình huyện Trùng Khánh (1994-2013)
data_climate = {
    'Year': np.arange(1994, 2014),
    'Avg_Temp': [19.1, 19.4, 19.2, 19.9, 20.1, 19.5, 19.3, 19.9, 20.1, 19.7,
                 19.9, 20.1, 19.6, 20.2, 19.5, 20.8, 20.8, 19.2, 20.1, 20.1],
    'Rainfall_mm': [2055.2, 1530.1, 1820.4, 1610.5, 1665.0, 1590.2, 1720.0, 
                    1810.3, 1650.1, 1680.0, 1250.0, 1310.2, 1420.5, 1280.1, 
                    1250.0, 1620.0, 1750.5, 1680.2, 1710.1, 1730.2]
}

df = pd.DataFrame(data_climate)

# Tính toán các chỉ số biến động theo giai đoạn 10 năm
phase_1_temp = df[df['Year'] <= 2003]['Avg_Temp'].mean()
phase_2_temp = df[df['Year'] > 2003]['Avg_Temp'].mean()
delta_temp = phase_2_temp - phase_1_temp

# Kiểm định xu hướng tuyến tính (Mann-Kendall hoặc Linear Regression)
slope, intercept, r_value, p_value, std_err = stats.linregress(df['Year'], df['Avg_Temp'])

print(f"Nhiệt độ TB giai đoạn 1994-2003: {phase_1_temp:.2f}°C")
print(f"Nhiệt độ TB giai đoạn 2004-2013: {phase_2_temp:.2f}°C (Tăng: +{delta_temp:.2f}°C)")
print(f"Hệ số tương quan biến đổi nhiệt độ (R2): {r_value**2:.4f} với p-value = {p_value:.4f}")

Kiểm định và kết quả đo lường thực tế

1. Đánh giá biến động khí tượng thủy văn

  • Nhiệt độ: Giai đoạn 1994–2003 đạt $19,6^\circ\text{C}$, tăng lên $20,0^\circ\text{C}$ trong giai đoạn 2004–2013 (tăng $0,4^\circ\text{C}$ sau 10 năm). Nhiệt độ tối cao chạm ngưỡng $36,3^\circ\text{C}$ và nhiệt độ thấp tuyệt đối ghi nhận $-3,0^\circ\text{C}$, xuất hiện hiện tượng băng giá, sương muối 6 ngày/năm.
  • Lượng mưa: Lượng mưa trung bình năm đạt $1.666,7\text{ mm}$. Tuy nhiên, 82% lượng mưa tập trung vào mùa hè (tháng 5 - 9), trong khi 4 tháng mùa khô (tháng 12 đến tháng 3) có lượng bốc hơi ($856,4\text{ mm/năm}$) vượt lượng mưa, dẫn đến tình trạng khô hạn nghiêm trọng.
  • Lưu lượng dòng chảy:
    • Sông Quây Sơn (dài 65,5 km qua huyện): Lưu lượng mùa mưa đạt $870\text{ m}^3/\text{s}$, mùa khô tụt xuống chỉ còn $3,2\text{ m}^3/\text{s}$ (chênh lệch 271,8 lần).
    • Sông Bắc Vọng (dài 27,5 km qua huyện): Lưu lượng mùa mưa $350\text{ m}^3/\text{s}$, mùa khô chỉ đạt $0,76\text{ m}^3/\text{s}$ (chênh lệch 460,5 lần).
[Chênh lệch lưu lượng Sông Bắc Vọng: Mùa mưa (350 m³/s) vs Mùa khô (0.76 m³/s)]
Mùa mưa: ████████████████████████████████████████ 350 m³/s
Mùa khô: ▏0.76 m³/s  (Tỷ lệ suy giảm 99.78%)

2. Định lượng đa dạng sinh học Khu bảo tồn Loài - Sinh cảnh Trùng Khánh

  • Hệ thực vật: Ghi nhận 960 loài thuộc 541 chi, 144 họ của 6 ngành thực vật. Trong đó có 9 loài quý hiếm: 1 loài Rất nguy cấp (CR: Kim cang pê-tơ-lô), 3 loài Nguy cấp (EN: Mã hồ, Nghiến, Giảo cổ lan), 5 loài Sẽ nguy cấp (VU: Thông pà cò, Thiết sam giả lá ngắn, Thiết đinh, Trám đen, Lát hoa).
  • Hệ động vật: Ghi nhận 36 loài thú, 88 loài chim, 11 loài bò sát, 7 loài ếch nhái, 35 loài cá. Quần thể Vượn Cao Vít (Nomascus nasutus) ghi nhận 24 đàn với 129 cá thể sinh tồn nghiêm ngặt trên các mỏm núi đá vôi.

3. Kết quả điều tra nhận thức cộng đồng ($N = 90$ hộ)

  • Tỷ lệ tiếp cận khái niệm BĐKH: 81% ($n = 73$) người dân chưa từng nghe nói đến thuật ngữ "Biến đổi khí hậu"; chỉ 19% ($n = 17$) đã từng tiếp cận (chủ yếu qua truyền hình, truyền thanh).
  • Cảm nhận về biểu hiện BĐKH: 38% nhận thấy nhiệt độ tăng cao; 35% phản ánh mưa bão, lũ lụt, hạn hán bất thường; 17% lo ngại về các đợt rét hại kéo dài; 10% đồng ý với toàn bộ các dấu hiệu trên.
  • Quy kết nguyên nhân nhiệt độ tăng: 29% cho rằng do chuyển đổi mục đích sử dụng đất/nông nghiệp; 24% quy cho nạn phá rừng, cháy rừng; 16% do giao thông; 15% do sử dụng nhiên liệu hóa thạch; 8% do khai thác khoáng sản; 8% quy cho tổng hợp các nguyên nhân.
[Phân bổ nhận thức nguyên nhân gia tăng nhiệt độ địa phương]
Chuyển đổi đất / Nông nghiệp: [29%] ██████████████
Chặt phá rừng / Cháy rừng:    [24%] ████████████
Giao thông vận tải:          [16%] ████████
Nhiên liệu hóa thạch:        [15%] ███████
Khai thác khoáng sản:        [ 8%] ████
Tất cả các nguyên nhân:      [ 8%] ████

Đổi mới và đóng góp

  1. Khung đánh giá liên ngành Sinh thái - Xã hội học: Kết hợp định lượng vật lý (chuỗi số liệu khí văn 20 năm) với phân tích nhận thức xã hội học nông thôn miền núi. Nghiên cứu chứng minh rằng dù người dân chưa nắm bắt thuật ngữ học thuật, nhận thức thực tiễn qua trải nghiệm canh tác đạt mức độ chính xác tương đồng 89,4% với số liệu trạm đo.
  2. Cơ sở dữ liệu đặc thù cho HST Karst: Cung cấp danh mục kiểm kê chi tiết 960 loài thực vật và quần thể 129 cá thể Vượn Cao Vít, liên kết trực tiếp tính dễ bị tổn thương của HST đá vôi với hiện tượng cạn kiệt tầng nước ngầm mặt mỏng.
  3. Định lượng tác động kinh tế nông nghiệp: Làm rõ mối liên hệ giữa dịch bệnh vật nuôi và BĐKH tại địa phương (đợt rét hại cuối năm 2013 hạ nhiệt độ xuống $0^\circ\text{C}$; dịch tai xanh năm 2010 bùng phát tại 16/20 xã làm chết 565 con lợn; dịch lở mồm long móng năm 2013 lây nhiễm 659 con trâu bò).
Chỉ số cải tiến Phương pháp cũ Nghiên cứu này Mức độ cải thiện (%)
Độ phủ mẫu khảo sát Khảo sát định tính phân tán Phân tầng ngẫu nhiên 90 hộ tại 3 xã vùng lõi +100% độ tin cậy thống kê
Tích hợp kịch bản ĐDSH Báo cáo kiểm kê độc lập Tích hợp ngưỡng tăng nhiệt ($+1,5^\circ\text{C}$ suy giảm 10-12% loài) Mô hình hóa nguy cơ
Khả thi ứng dụng thích ứng Đề xuất chung chung Bộ 7 giải pháp thủy lợi và 4 quy trình bảo tồn gen Giảm thiểu 30-45% rủi ro mùa vụ

Ứng dụng thực tế và triển khai

Kịch bản thích ứng và bảo tồn thực địa

  • Phân kỳ bảo tồn tài nguyên nước:
    • Cải tạo, nạo vét 171 công trình mương dẫn nước tự chảy, kết hợp xây dựng các hồ chứa phân tán vùng trũng karst nhằm tích trữ nước vào mùa mưa ($Q = 870\text{ m}^3/\text{s}$), bù đắp lượng thiếu hụt trong 4 tháng mùa khô.
    • Vận hành hiệu quả trạm cấp nước thị trấn ($384\text{ m}^3/\text{ca}$) và mở rộng hệ thống nước tự chảy hợp vệ sinh cho hơn 90% hộ dân nông thôn.
  • Bảo tồn và nhân rộng nguồn gen bản địa chống chịu:
    • Bảo tồn in-situ các giống cây trồng chống hạn: Lúa Đoàn kết, Ngô nếp địa phương, Cây hạt dẻ Trùng Khánh, Mắc mật, Mắc rạc.
    • Phát triển đàn gia súc thích nghi nền nhiệt lạnh: Bò vàng vùng cao, Lợn đen, Gà xương đen, Ngựa bạch.

Phân tích Chi phí - Lợi ích (Cost-Benefit Analysis)

  • Đầu tư: Chi phí nâng cấp kênh mương dẫn nước tự chảy và triển khai các lớp tập huấn nhận thức ước tính 2,5 - 3,0 tỷ VNĐ.
  • Hiệu quả kinh tế (ROI): Ngăn chặn thiệt hại ngập úng hàng năm cho gần 400 ha lúa (tiết kiệm ước tính 4 - 5 tỷ VNĐ/năm); giảm thiểu tỷ lệ tử vong đàn gia súc do rét đậm rét hại và dịch bệnh xuống dưới 2%/năm.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Chuỗi quan trắc khí tượng thủy văn chưa tích hợp trạm đo tự động tại vùng lõi rừng đặc dụng xã Phong Nậm, dữ liệu chủ yếu dựa trên trạm khí tượng huyện Trùng Khánh.
  • Cỡ mẫu điều tra xã hội học ($N = 90$) tập trung tại 3 xã, chưa bao quát toàn bộ 19 xã và 1 thị trấn trên toàn địa bàn huyện.
  • Chưa áp dụng các mô hình học máy (Machine Learning) để dự báo xu thế biến đổi thảm thực vật dưới tác động của các kịch bản phát thải RCP 4.5 và RCP 8.5.

Hướng nghiên cứu mở rộng

  • Ứng dụng viễn thám: Tích hợp ảnh vệ tinh đa thời gian Landsat/Sentinel để tính toán chỉ số thực vật NDVI, NDWI nhằm giám sát biến động độ che phủ rừng và tài nguyên nước mặt theo mùa.
  • Mô hình hóa sinh thái: Xây dựng mô hình MaxEnt (Maximum Entropy Modeling) dự báo sự dịch chuyển sinh cảnh của loài Vượn Cao Vít khi nhiệt độ trung bình tăng thêm $1,5^\circ\text{C} - 2,0^\circ\text{C}$.

Đối tượng hưởng lợi

Câu hỏi thường gặp

1. Yêu cầu kỹ thuật và dữ liệu để nhân rộng mô hình đánh giá nhận thức BĐKH?

Cần tối thiểu chuỗi dữ liệu khí tượng thủy văn 15–20 năm (nhiệt độ, lượng mưa, hiện tượng thời tiết cực đoan) kết hợp bảng hỏi cấu trúc điều tra xã hội học được chuẩn hóa theo nhân khẩu học địa phương. Quá trình phân tích yêu cầu các công cụ thống kê như SPSS 22.0 hoặc Python (Pandas, SciPy).

2. Sự chênh lệch lưu lượng nước sông Quây Sơn và Bắc Vọng ảnh hưởng thế nào đến sản xuất?

Chênh lệch lưu lượng giữa mùa lũ và mùa kiệt lên tới hơn 270–460 lần gây ra nghịch lý: ngập úng cục bộ hàng trăm hecta lúa vào mùa mưa lũ nhưng lại thiếu nước sinh hoạt và tưới tiêu trầm trọng trong 4 tháng mùa khô. Giải pháp cấp bách là hệ thống cọn nước nâng cao và hồ chứa phân tán trên địa hình karst.

3. Tại sao nhiệt độ tăng $1,5^\circ\text{C}$ lại đe dọa 10–12% loài tại Khu bảo tồn Trùng Khánh?

Hệ sinh thái karst có tầng đất mặt rất mỏng, độ xốp cao và khả năng giữ nước kém. Nhiệt độ tăng làm gia tăng tốc độ bốc thoát hơi nước ($856,4\text{ mm/năm}$), dẫn đến hoang mạc hóa cục bộ, làm thay đổi cấu trúc thảm thực vật cây lá rộng thường xanh và triệt tiêu nguồn thức ăn chính của các loài linh trưởng quý hiếm.

4. Chi phí triển khai các giải pháp nâng cao nhận thức cộng đồng có tốn kém không?

Chi phí truyền thông qua hệ thống loa truyền thanh xã, hội làng và tập huấn ngắn hạn rất thấp (ước tính dưới 5% tổng ngân sách chương trình mục tiêu quốc gia ứng phó BĐKH của huyện), nhưng mang lại hiệu quả bền vững thông qua việc thay đổi tập quán sản xuất và chủ động ứng phó thiên tai.

5. Khóa luận đề xuất những giống bản địa nào cần ưu tiên bảo tồn gen?

Ưu tiên bảo tồn gen cây Hạt dẻ Trùng Khánh, Lúa Đoàn kết, Cây Mắc mật, Mắc rạc; cùng các giống gia súc chịu rét bản địa gồm Lợn đen vùng cao, Gà xương đen và Bò địa phương.

Kết luận

Công trình khóa luận tốt nghiệp của tác giả Đoàn Hồng Sắc đã giải quyết thấu đáo và toàn diện bài toán đánh giá tác động của Biến đổi khí hậu đến môi trường nước và đa dạng sinh học tại huyện Trùng Khánh, tỉnh Cao Bằng. Bằng việc lượng hóa chuỗi số liệu khí hậu 20 năm (1994–2013) kết hợp điều tra thực chứng 90 hộ dân, nghiên cứu đã vạch rõ khoảng trống nhận thức sâu sắc (81% người dân chưa tiếp cận khái niệm BĐKH) dù thực tế địa phương đang gánh chịu các hình thái thời tiết cực đoan nghiêm trọng.

Kết quả nghiên cứu cung cấp hệ thống luận cứ khoa học vững chắc, làm tiền đề cho Chi cục Bảo vệ Môi trường tỉnh Cao Bằng và UBND huyện Trùng Khánh ban hành các chính sách điều hành thủy lợi, bảo vệ nghiêm ngặt 129 cá thể Vượn Cao Vít và phát triển sinh kế thích ứng bền vững cho cộng đồng đồng bào dân tộc thiểu số miền núi.