Tổng quan nghiên cứu
Khu vực lưu vực sông Đà thuộc địa phận tỉnh Sơn La có diện tích tự nhiên khoảng 9.005,2 km2, chiếm 64,1% tổng diện tích toàn tỉnh (14.174,44 km2). Đây là địa bàn sinh sống của hơn 1.092.700 người thuộc 12 dân tộc anh em, đồng thời là trung tâm năng lượng chiến lược quốc gia. Việc triển khai Dự án Thủy điện Sơn La theo Quyết định số 92/QĐ-TTg của Thủ tướng Chính phủ với công suất thiết kế 2.400 MW, dung tích hồ chứa đạt 9,26 tỷ m3 nước và diện tích lưu vực 43.760 km2 đã tạo ra những tác động mang tính lịch sử đối với không gian sinh thái lãnh thổ. Quá trình ngăn dòng, tích nước và di dời 18.968 hộ dân (tương đương khoảng 90.000 nhân khẩu) đã thúc đẩy làn sóng chuyển dịch cơ cấu sử dụng đất và lớp phủ mặt đất với quy mô chưa từng có tiền lệ.
Mục tiêu cốt lõi của nghiên cứu là ứng dụng tích hợp công nghệ viễn thám (Remote Sensing - RS) và Hệ thông tin địa lý (GIS) nhằm lượng hóa chi tiết biến động sử dụng đất và lớp phủ mặt đất qua 3 mốc thời gian: năm 2003 (trước khởi công), năm 2009 (giai đoạn thi công cao điểm) và năm 2014 (giai đoạn vận hành ổn định). Phạm vi nghiên cứu được phân định thành 2 vùng sinh thái đặc trưng: khu vực thượng lưu trên đập có diện tích 217.307,9 ha (chiếm 15,3% diện tích tự nhiên tỉnh) và khu vực hạ lưu dưới đập có diện tích 683.209,6 ha (chiếm 48,2% diện tích tự nhiên tỉnh). Kết quả nghiên cứu cung cấp hệ thống cơ sở dữ liệu số hóa có độ chuẩn xác cao, hỗ trợ đắc lực cho công tác quản trị tài nguyên, tối ưu hóa quy hoạch đất đai và bảo đảm sinh kế bền vững cho đồng bào vùng tái định cư Tây Bắc.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết áp dụng
Nghiên cứu được xây dựng trên nền tảng lý thuyết phân loại lớp phủ bề mặt phân cấp theo Anderson (1976) kết hợp hệ thống phân loại thích ứng với điều kiện Việt Nam của Nguyễn Ngọc Thạch (2005). Cấu trúc phân loại gồm 7 nhóm cấp 1 (đô thị, đất nông nghiệp, đất bỏ hoang, đất rừng, mặt nước, đất ướt, đất hoang) và 27 phân lớp cấp 2, phản ánh đầy đủ đặc trưng thảm phủ vùng núi cao.
Về nguyên lý viễn thám, nghiên cứu ứng dụng quy luật phản xạ phổ điện từ của các vật thể tự nhiên và chỉ số khác biệt thực vật chuẩn hóa NDVI (Normalized Difference Vegetation Index) theo công thức toán học: NDVI = (NIR - RED) / (NIR + RED). Trong đó, NIR là độ phản xạ phổ hồng ngoại gần (bước sóng 0,75 - 1,0 µm) và RED là độ phản xạ dải sóng đỏ (bước sóng 0,55 - 0,7 µm). Chỉ số này cho phép bóc tách chính xác sinh khối thực vật dày hay thưa và phân định ranh giới giữa thảm rừng, nương rẫy với bề mặt đất trống, mặt nước. Bên cạnh đó, luận văn vận dụng lý thuyết động lực cảnh quan sinh thái lưu vực để phân tích mối tương quan giữa sự xuất hiện của công trình nhân tạo quy mô lớn và quá trình suy thoái hoặc phục hồi tài nguyên đất đai.
Ba khái niệm trung tâm được chuẩn hóa xuyên suốt gồm: Lớp phủ mặt đất (Land cover - các đối tượng vật lý quan sát được trên bề mặt), Hiện trạng sử dụng đất (Land use - mục đích khai thác kinh tế - xã hội của con người) và Biến động sử dụng đất (sự chuyển hóa không gian và chức năng của các khoanh đất theo thời gian).
Phương pháp nghiên cứu
Nguồn dữ liệu đầu vào bao gồm chuỗi ảnh vệ tinh Landsat đa thời gian (Landsat 5 TM năm 2003, Landsat 7 ETM+ năm 2009 và Landsat 8 OLI năm 2014) kết hợp bản đồ địa hình tỷ lệ 1:25.000, 1:50.000, mô hình số bề mặt DSM/DEM tỷ lệ 1:25.000 và hệ thống dữ liệu nền địa lý tỉnh Sơn La.
Để đảm bảo độ tin cậy thống kê, nghiên cứu thiết lập cỡ mẫu gồm 350 điểm kiểm chứng thực địa phân bố trên toàn bộ 900.520 ha của lưu vực. Phương pháp chọn mẫu ngẫu nhiên phân tầng (Stratified Random Sampling) được áp dụng để đảm bảo mỗi lớp phủ mặt đất đều có số lượng điểm đại diện tương ứng với tỷ lệ diện tích thực tế.
Thuật toán phân loại có kiểm định Maximum Likelihood trên phần mềm ENVI kết hợp phương pháp so sánh sau phân loại (phương pháp bản đồ - bản đồ) trong môi trường ArcGIS Desktop (Spatial Analyst Tools) được lựa chọn là giải pháp kỹ thuật tối ưu. Lý do lựa chọn là thuật toán Maximum Likelihood khai thác triệt để xác suất thống kê Bayes dựa trên các vector đặc trưng phổ đa chiều, hạn chế tối đa hiện tượng nhầm lẫn phổ do bóng núi và sương mù tại địa hình chia cắt phức tạp. Phương pháp bản đồ - bản đồ cho phép xuất trực tiếp ma trận biến động không gian (From-To matrix), chỉ rõ nguồn gốc và đích đến của từng hecta đất bị chuyển đổi giữa các chu kỳ nghiên cứu.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Thứ nhất, tại khu vực trên đập (217.307,9 ha), diện tích mặt nước chuyên dùng và lòng hồ chứa nước đã mở rộng vượt bậc trong giai đoạn 2009 - 2014 khi hồ thủy điện hoàn thành tích nước lên cao trình 215 m. Biến động này làm chìm ngập hoàn toàn hàng nghìn hecta đất lúa màu màu mỡ ven sông suối thuộc huyện Quỳnh Nhai và Mường La, biến các thung lũng nông nghiệp trũng thấp thành thủy vực rộng lớn với dung tích 9,26 tỷ m3 nước.
Thứ hai, không gian cư trú và canh tác nông nghiệp có sự dịch chuyển mạnh mẽ lên các sườn dốc cao hơn. Việc di chuyển 18.968 hộ dân tái định cư đã thúc đẩy diện tích đất xây dựng công trình hạ tầng và điểm dân cư nông thôn mới tăng trên 45% tại các huyện tiếp nhận, kéo theo sự thu hẹp của diện tích đất nương đồi bỏ hoang gần các trục lộ giao thông mới dài 50 km.
Thứ ba, cơ cấu thảm thực vật có sự chuyển biến tích cực tại các khu vực bảo vệ nghiêm ngặt. Tổng diện tích 3 loại rừng trên địa bàn tỉnh đạt 624.380,8 ha vào năm 2013, trong đó đất rừng phòng hộ chiếm tỷ trọng lớn nhất với 397.345,4 ha (chiếm 63,6%), đất rừng sản xuất đạt 178.920,4 ha (chiếm 28,7%) và đất rừng đặc dụng đạt 48.115,0 ha (chiếm 7,7%), nâng tỷ lệ che phủ rừng toàn tỉnh lên mức 50,0%.
Thứ tư, tại khu vực dưới đập (683.209,6 ha), áp lực canh tác nông nghiệp trên đất dốc gia tăng đáng kể. Thống kê theo độ dốc địa hình năm 2014 cho thấy diện tích canh tác trên các sườn núi có độ dốc từ 15 đến trên 20 độ chiếm tỷ lệ hơn 60% tổng diện tích đất nông nghiệp toàn vùng, làm gia tăng nguy cơ thoái hóa đất cục bộ tại Bắc Yên và Phù Yên.
Thảo luận kết quả
Nguyên nhân sâu xa của các biến động nêu trên bắt nguồn từ sự tương tác giữa quyết định đầu tư hạ tầng năng lượng tầm vóc quốc gia và quá trình thích ứng sinh kế của cộng đồng dân cư bản địa. So sánh với các nghiên cứu tương tự tại các công trình thủy điện lớn ở Đông Nam Á, tốc độ biến đổi cảnh quan tại lưu vực sông Đà diễn ra tập trung và triệt để hơn nhờ chiến dịch di dân tái định cư thần tốc hoàn thành vào tháng 4 năm 2010.
Về phương diện trực quan hóa dữ liệu, các diễn biến phức tạp của tài nguyên đất được minh họa trực quan thông qua:
- Hệ thống ma trận chuyển dịch sử dụng đất (Land transition matrix) thể hiện cụ thể diện tích luân chuyển giữa 7 loại đất qua các giai đoạn 2003 - 2009 và 2009 - 2014.
- Chuỗi biểu đồ cột so sánh tương quan tăng giảm diện tích từng lớp phủ giữa hai vùng trên đập và dưới đập.
- Bản đồ chuyên đề phân vùng độ dốc kết hợp mô hình DEM 1:25.000, giúp chỉ rõ các điểm nóng có nguy cơ sạt lở và xói mòn đất khi thảm thực vật bị suy giảm.
Đề xuất và khuyến nghị
Thứ nhất, quy hoạch và khai thác hiệu quả tài nguyên đất bán ngập vùng lòng hồ thủy điện:
- Hành động: Tổ chức cắm mốc ranh giới, thiết lập mô hình luân canh cây ngắn ngày (ngô lai, rau màu chịu ngập, cây họ đậu) xen kẽ nuôi trồng thủy sản trên diện tích ước tính khoảng 2.500 ha đất bán ngập khi mực nước hồ dao động từ cao trình 175 m đến 215 m.
- Mục tiêu: Nâng cao giá trị thu nhập đạt trên 60 triệu đồng/ha/năm, ổn định đời sống cho hơn 5.000 hộ dân lòng hồ.
- Lộ trình: Triển khai từ năm 2024 đến năm 2026.
- Chủ thể thực hiện: Sở Nông nghiệp và Phát triển nông thôn tỉnh Sơn La chủ trì, phối hợp với Ủy ban Nhân dân các huyện Quỳnh Nhai, Mường La, Thuận Châu.
Thứ hai, nhân rộng mô hình canh tác nông lâm kết hợp trên đất dốc (SALT - Sloping Agricultural Land Technology):
- Hành động: Xây dựng các băng rào chắn sinh học bằng cỏ Vetiver, cây cốt khí hoặc keo dậu xen kẽ các hàng cây ăn quả (xoài, nhãn ghép, mận hậu) trên các sườn đồi có độ dốc từ 15 đến 25 độ.
- Mục tiêu: Giảm thiểu 70% lượng đất xói mòn rửa trôi, cải thiện độ phì nhiêu tầng mặt trên diện tích 20.000 ha đất nương rẫy thoái hóa.
- Lộ trình: Giai đoạn 2024 - 2028.
- Chủ thể thực hiện: Chi cục Trồng trọt và Bảo vệ Thực vật tỉnh Sơn La kết hợp với Trung tâm Khuyến nông địa phương.
Thứ ba, thiết lập hệ thống WebGIS giám sát viễn thám định kỳ tài nguyên lưu vực:
- Hành động: Ứng dụng ảnh vệ tinh độ phân giải cao kết hợp thuật toán học máy để cập nhật biến động lớp phủ mặt đất định kỳ 6 tháng một lần, tự động cảnh báo các điểm phá rừng phòng hộ và lấn chiếm hành lang thoát lũ.
- Mục tiêu: Rút ngắn 80% thời gian phát hiện vi phạm quy hoạch đất đai, bảo vệ an toàn cho 397.345,4 ha rừng phòng hộ đầu nguồn.
- Lộ trình: Vận hành từ quý 1 năm 2025.
- Chủ thể thực hiện: Sở Tài nguyên và Môi trường tỉnh Sơn La điều hành, phối hợp Công ty Thủy điện Sơn La.
Thứ tư, tái cơ cấu sinh kế hậu tái định cư thông qua nuôi trồng thủy sản lồng bè và du lịch sinh thái:
- Hành động: Hỗ trợ vốn vay ưu đãi, tập huấn kỹ thuật nuôi cá lồng công nghệ cao trên mặt nước hồ chứa dung tích 9,26 tỷ m3 kết hợp phát triển du lịch trải nghiệm văn hóa lòng hồ sông Đà.
- Mục tiêu: Đưa sản lượng thủy sản nuôi trồng tăng trưởng trên 15%/năm, tạo việc làm ổn định cho hơn 10.000 lao động nông thôn.
- Lộ trình: Thực hiện từ năm 2024 đến năm 2030.
- Chủ thể thực hiện: Ủy ban Nhân dân tỉnh Sơn La phối hợp Tập đoàn Điện lực Việt Nam (EVN).
Đối tượng nên tham khảo luận văn
Nhóm 1: Cơ quan quản lý nhà nước về tài nguyên, môi trường và nông nghiệp (Sở TN&MT, Sở NN&PTNT, Chi cục Lâm nghiệp). Nhóm này có thể khai thác hệ thống bản đồ số hóa tỷ lệ 1:25.000 và quy trình giám sát lớp phủ tự động để phục vụ công tác lập quy hoạch sử dụng đất cấp tỉnh, thẩm định dự án giao đất lâm nghiệp và quản lý an toàn hồ đập.
Nhóm 2: Đơn vị tư vấn đánh giá tác động môi trường (ĐTM) và các chủ đầu tư dự án năng lượng thủy điện quy mô lớn. Nhóm này có thể kế thừa khung phương pháp luận lượng hóa tác động không gian đa thời gian của các công trình công suất trên 1.000 MW đến cấu trúc thảm phủ sinh thái và dân sinh kinh tế.
Nhóm 3: Giảng viên, nhà nghiên cứu, nghiên cứu sinh và học viên cao học chuyên ngành Bản đồ, Viễn thám, GIS, Địa lý tự nhiên và Quản lý đất đai. Tài liệu cung cấp quy trình chuẩn mực từ tiền xử lý ảnh Landsat, hiệu chỉnh khí quyển, chọn mẫu kiểm định thực địa đến thuật toán Maximum Likelihood và đánh giá ma trận nhầm lẫn.
Nhóm 4: Chính quyền địa phương cấp huyện, xã và các tổ chức phi chính phủ phát triển cộng đồng tại vùng Tây Bắc. Nguồn tư liệu này cung cấp các luận cứ khoa học thực tiễn để định hướng phân bổ nguồn vốn hỗ trợ sinh kế, chuyển giao mô hình canh tác SALT và bảo tồn vốn rừng đầu nguồn.
Câu hỏi thường gặp
Công nghệ viễn thám và GIS mang lại ưu thế vượt trội gì so với phương pháp điều tra mặt đất truyền thống trong nghiên cứu biến động đất đai vùng núi?
Phương pháp viễn thám và GIS cho phép bao quát toàn diện diện tích lưu vực rộng lớn 9.005,2 km2 có địa hình đồi núi dốc đứng hiểm trở. Ứng dụng này giúp rút ngắn chu kỳ kiểm kê đất đai từ 18 tháng xuống còn vài tuần xử lý ảnh số, tiết kiệm trên 60% kinh phí điều tra và đảm bảo tính đồng nhất cao về mặt không gian và thời gian.
Tại sao luận văn lại tách biệt phạm vi nghiên cứu thành hai khu vực trên đập và dưới đập thủy điện Sơn La?
Việc phân chia thành khu vực trên đập (217.307,9 ha) và dưới đập (683.209,6 ha) nhằm bóc tách cơ chế tác động chuyên biệt. Vùng trên đập chịu ảnh hưởng trực tiếp của việc ngập nước lòng hồ 9,26 tỷ m3 và di dời 18.968 hộ dân, trong khi vùng dưới đập phản ánh tác động gián tiếp từ sự biến đổi dòng chảy và áp lực gia tăng canh tác đất dốc.
Độ chính xác của các bản đồ hiện trạng và biến động sử dụng đất trong nghiên cứu được kiểm soát như thế nào?
Độ chính xác được kiểm soát chặt chẽ qua việc nắn chỉnh hình học theo bản đồ địa hình 1:25.000, hiệu chỉnh sai số bức xạ khí quyển và đối chiếu với 350 điểm kiểm tra thực địa độc lập. Ma trận nhầm lẫn cho thấy độ chính xác tổng thể của các bản đồ phân loại đạt trên 85% và hệ số Kappa vượt mức 0,82, đáp ứng đầy đủ tiêu chuẩn đo đạc bản đồ quốc gia.
Mô hình canh tác trên đất dốc SALT giải quyết thách thức môi trường nào tại lưu vực sông Đà?
Mô hình SALT giải quyết triệt để tình trạng xói mòn và rửa trôi đất canh tác nương rẫy trên các sườn dốc từ 15 đến 25 độ tại Sơn La. Bằng cách thiết lập các dải băng cây xanh cố định đạm chắn ngang sườn dốc, mô hình giúp giữ lại hơn 70% lượng phù sa tầng mặt, duy trì độ phì nhiêu và hạn chế bồi lắng lòng hồ thủy điện.
Dữ liệu không gian của luận văn có khả năng tích hợp vào các hệ sinh thái đô thị thông minh và hạ tầng dữ liệu không gian quốc gia (NSDI) không?
Toàn bộ cơ sở dữ liệu địa lý trong nghiên cứu được xây dựng đồng bộ trên hệ quy chiếu quốc gia VN-2000, múi chiếu 3 độ với cấu trúc phân lớp chuẩn mực. Dữ liệu này hoàn toàn tương thích để tích hợp trực tiếp vào hạ tầng WebGIS của tỉnh Sơn La và cổng dữ liệu không gian quốc gia phục vụ giám sát trực tuyến theo thời gian thực.
Kết luận
- Luận văn đã xây dựng thành công quy trình công nghệ tích hợp viễn thám và GIS để thành lập bộ bản đồ hiện trạng và biến động sử dụng đất tỷ lệ 1:25.000 qua 3 mốc lịch sử (2003, 2009, 2014) trên diện tích lưu vực sông Đà rộng 9.005,2 km2.
- Nghiên cứu đã định lượng hóa sự chuyển dịch không gian sâu sắc dưới tác động của Công trình Thủy điện Sơn La công suất 2.400 MW, đặc biệt là sự mở rộng đột biến của mặt nước hồ chứa 9,26 tỷ m3 và sự hình thành các điểm dân cư tái định cư cho 18.968 hộ dân.
- Phân tích làm rõ sự phân hóa sinh thái rõ rệt giữa vùng thượng lưu trên đập (217.307,9 ha) chịu tác động dâng ngập và vùng hạ lưu dưới đập (683.209,6 ha) chịu áp lực canh tác nông nghiệp trên đất dốc.
- Đề xuất hệ thống giải pháp kỹ thuật có tính ứng dụng cao, bao gồm khai thác quỹ đất bán ngập, nhân rộng mô hình nông lâm kết hợp SALT và thiết lập hệ thống giám sát vệ tinh định kỳ để bảo vệ 397.345,4 ha rừng phòng hộ.
- Công trình khẳng định giá trị phương pháp luận khoa học vững chắc, cung cấp nguồn tư liệu tham khảo thiết thực cho các nhà hoạch định chiến lược trong quản lý lưu vực sông và quy hoạch thủy điện bền vững.
Đóng góp lớn nhất của luận văn là đã chuẩn hóa phương pháp tiếp cận công nghệ không gian hiện đại để giải quyết bài toán quản trị tài nguyên phức tạp tại các vùng miền núi nhạy cảm sinh thái. Bước tiếp theo trong giai đoạn 2025 - 2030, các cơ quan chức năng cần tích hợp mô hình phân tích không gian này vào hệ thống WebGIS giám sát thời gian thực của tỉnh. Hãy liên hệ ngay với tác giả hoặc cơ sở đào tạo để tiếp cận toàn văn cơ sở dữ liệu số và bản đồ chi tiết phục vụ cho các dự án phát triển bền vững tại địa phương.