Bối cảnh và vấn đề nghiên cứu

Lý do chọn đề tài và bối cảnh thực tiễn

Xói mòn đất trong sản xuất nông nghiệp là vấn đề môi trường và kinh tế cấp bách, ảnh hưởng trực tiếp đến độ phì nhiêu, chất lượng đất canh tác, làm suy giảm năng suất cây trồng và gây ô nhiễm môi trường thứ cấp. Xói mòn và suy thoái đất cũng gắn liền với các hiện tượng thiên tai và biến đổi khí hậu toàn cầu. Việc ứng dụng công nghệ thông tin và hệ thống thông tin địa lý (GIS) vào nghiên cứu xói mòn đất giúp lượng hóa mức độ mất đất, hỗ trợ xây dựng các giải pháp canh tác bền vững.

Trong các mô hình tính toán xói mòn hiện nay, phương trình mất đất phổ quát sửa đổi (RUSLE) được áp dụng rộng rãi nhưng vẫn bộc lộ một số hạn chế:

  • Tồn tại sai số đáng kể so với số liệu đo đạc thực tế, đặc biệt là tại các khu vực xảy ra đồng thời quá trình xói mòn và bồi tụ (như hệ thống kênh rạch, sông ngòi, vùng cửa sông).
  • Các công thức tính toán phức tạp đòi hỏi nhiều biến số đầu vào chi tiết, dẫn đến tốc độ xử lý dữ liệu và kết hợp mô hình bị chậm trên diện tích lớn.
  • Khó khăn trong việc chuẩn hóa và cấu hình các hệ số thành phần (factors) sao cho phù hợp với đặc thù khí hậu nhiệt đới gió mùa và hai dạng địa hình chính của Việt Nam là vùng đồi núi và đồng bằng.

Mục tiêu và nhiệm vụ nghiên cứu

Đề tài xác định 3 mục tiêu và nhiệm vụ cụ thể:

  1. Xây dựng bản đồ phân loại đất qua các năm: Phản ánh sự phân bố và quy luật chuyển dịch của các loại đất trên lãnh thổ, làm cơ sở nhận diện các khu vực có nguy cơ xói mòn.
  2. Thiết lập mô hình RUSLE tùy biến: Lấy hệ số xói mòn của đất (K-factor) làm trung tâm, kết hợp với các hệ số địa hình (LS-factor), lớp phủ (C-factor), biện pháp canh tác (P-factor) và lượng mưa (R-factor) để thử nghiệm các tổ hợp factor tối ưu, xác định chính xác các vùng bị xói mòn và bồi tụ.
  3. Cải tiến hiệu năng tính toán và giảm thiểu sai số: Rút ngắn thời gian xử lý của mô hình, khắc phục nhược điểm của các nghiên cứu trước ở khả năng mô phỏng xói mòn - bồi tụ tại khu vực mạng lưới sông ngòi, kênh rạch và các lưu vực đặc thù.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các loại đất và đặc tính thành phần lý hóa của đất; dữ liệu lượng mưa; độ dốc và chiều dài sườn dốc địa hình; các lớp dữ liệu bao phủ/sử dụng đất; các hệ số thành phần trong mô hình USLE, RUSLE, RUSLE2.
  • Đối tượng thụ hưởng/ứng dụng: Nông dân và ngành nông nghiệp; sinh viên, giảng viên và cán bộ nghiên cứu chuyên ngành Địa lý, Địa chất, Quản lý tài nguyên môi trường và Hệ thống thông tin địa lý.
  • Phạm vi không gian: Toàn bộ phạm vi đất liền của lãnh thổ Việt Nam (không bao gồm phần biển đảo trong giai đoạn tính toán mô hình thô; kế hoạch tích hợp không gian biển đảo thuộc giai đoạn hoàn thiện và triển khai hệ thống).
  • Phạm vi thời gian thực hiện: Từ ngày 10/03/2023 đến ngày 15/06/2023.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết và mô hình tính toán

Nghiên cứu dựa trên nền tảng lý thuyết của các phương trình xói mòn đất cơ bản:

  • Phương trình mất đất phổ quát USLE (Wischmeier & Smith, 1965/1978): Mô hình ước tính mất đất trung bình năm: $$A = R \times K \times L \times S \times C \times P$$
  • Phương trình sửa đổi RUSLE: Cải tiến từ USLE, tích hợp hệ số sinh thái sườn dốc $LS$ thay cho việc tách rời $L$ và $S$, đồng thời cho phép đánh giá trên nhiều loại bề mặt đất khác nhau (đất nông nghiệp, đất rừng xáo trộn, đất khai khoáng, công trường xây dựng): $$A = R \times K \times LS \times C \times P$$
  • Phương trình RUSLE2: Mở rộng khả năng tính toán trầm tích lắng đọng và bồi tụ theo ngày, nhưng yêu cầu tập dữ liệu đầu vào rất chi tiết và giải thuật phức tạp.
  • Phương trình MUSLE: Sử dụng các yếu tố dòng chảy mặt để cải thiện việc tính toán lượng mưa.
Mô hình Ưu điểm chính Nhược điểm chính
USLE Dễ sử dụng, dữ liệu đầu vào ít, áp dụng tốt cho đất trồng trọt. Số lượng yếu tố tính toán hạn chế, sai số cao khi áp dụng ngoài đất nông nghiệp chuẩn, giả định dòng chảy đồng đều.
RUSLE Ứng dụng cho nhiều loại hình sử dụng đất, xử lý được địa hình phân kỳ/hội tụ, tính toán tốt hơn tác động của thảm phủ và độ dốc. Tính toán phức tạp hơn USLE, tốc độ xử lý chậm, còn sai số ở vùng bồi tụ và kênh rạch nếu giữ nguyên công thức chuẩn.
RUSLE2 Kết quả có độ chính xác cao, tính toán chi tiết sự sản sinh và lắng đọng trầm tích theo chu kỳ ngày. Yêu cầu dữ liệu đầu vào khắt khe, khó thu thập đầy đủ ở quy mô lớn, thuật toán tính toán phức tạp.

Các hệ số trong mô hình gồm:

  • Hệ số $K$ (Soil Erodibility Factor): Tính toán theo công thức của Goldman et al. (1986) và Wischmeier et al. (1971), dựa trên tỷ lệ cát mịn, bùn (silt), sét (clay), phần trăm chất hữu cơ ($P_{om}$), chỉ số kết cấu đất ($S_{truc}$) và hệ số thấm mặt cắt ($f_{perm}$).
  • Hệ số $R$ (Rainfall Erosivity Factor): Sử dụng công thức Hurni (1985) nghiên cứu cho vùng có đặc thù khí hậu nhiệt đới tương đồng (Ethiopia), loại bỏ thành phần băng tuyết $R_g$: $R = 0.12 \times P$ (với $P$ là lượng mưa năm).
  • Hệ số $LS$ (Slope Length and Steepness Factor): Tính toán dựa trên lớp tích lũy dòng chảy (Flow Accumulation) và độ dốc trích xuất từ mô hình số độ cao (DEM) theo cải tiến của Jan D. Moore và John P.
  • Hệ số $C$ (Cover Management Factor) và Hệ số $P$ (Support Practice Factor): Xác định theo phân loại độ phủ đất và kỹ thuật canh tác (làm ruộng bậc thang, canh tác theo đường đồng mức, trồng xen băng) dựa trên các bảng định mức chuẩn (Shin, 1999).

Phương pháp nghiên cứu thực tế

  • Xử lý dữ liệu không gian Raster và Vector: Sử dụng cấu trúc dữ liệu lưới Raster (ma trận điểm ô vuông) để biểu diễn các biến số liên tục (DEM, lượng mưa, độ che phủ) và dữ liệu Vector (Polygon) để biểu diễn ranh giới khu vực hành chính, vùng đất.
  • Phân tích và mô hình hóa không gian trên hệ sinh thái ArcGIS: Sử dụng phần mềm ArcGIS Pro, công cụ ArcInfo để biên tập, tính toán ma trận raster; sử dụng hệ quản trị cơ sở dữ liệu địa lý (GeoDatabase) để tổ chức các bảng dữ liệu không gian (Spatial Data Tables) và phi không gian (Nonspatial Data Tables).
  • Phương pháp thực nghiệm tổ hợp Factor: Thử nghiệm và đánh giá so sánh 5 cấu hình mô hình tùy biến gồm: KLS, KLSP, KLSC, KLSPC và KC nhằm xác định mô hình tối ưu về tốc độ và khả năng phản ánh đúng thực tế xói mòn - bồi tụ tại Việt Nam.

Nguồn dữ liệu sử dụng

  1. FAO-DSMW (Digital Soil Map of the World): Bộ dữ liệu bản đồ đất kỹ thuật số toàn cầu của FAO, cung cấp thông tin phân loại đất theo chuẩn WRB và các thuộc tính lý hóa đất (độ phì, độ sâu, hàm lượng hữu cơ, kết cấu).
  2. GLC-SHARE (Global Land Cover-SHARE): Dữ liệu lớp phủ bề mặt toàn cầu của FAO với độ phân giải 30 giây cung, gồm 11 lớp phủ chính.
  3. HydroSHEDS Core: Dữ liệu thủy văn toàn cầu trích xuất từ dữ liệu radar địa hình SRTM (NASA/USGS), cung cấp mô hình DEM hiệu chỉnh, hướng dòng chảy và tích lũy dòng chảy (ACC grids cấp độ Continental và Global) ở các độ phân giải 3, 15, 30 giây cung.
  4. CRU-TS (Climatic Research Unit - University of East Anglia): Chuỗi dữ liệu khí hậu dạng lưới toàn cầu từ năm 1901 đến 2022 với độ phân giải $0.5^\circ \times 0.5^\circ$ (~50 km), cung cấp các thông số nhiệt độ, lượng mưa, áp suất hơi nước và tần suất ngày ẩm ướt.
Bộ dữ liệu Đơn vị phát triển / Nguồn Độ phân giải Thuộc tính trích xuất chính
FAO-DSMW FAO Toàn cầu Phân loại đất (WRB), thành phần cơ giới, chất hữu cơ
GLC-SHARE FAO 30 giây cung 11 lớp phủ bề mặt (rừng, đất trọt, đồng cỏ, đất trần...)
HydroSHEDS Core Univ. of California & WWF (từ SRTM - NASA/USGS) 3, 15, 30 giây cung DEM độ cao, Flow Accumulation (ACC), hướng dòng chảy
CRU-TS CRU - Univ. of East Anglia 0.5 độ (~50 km) Lượng mưa tháng/năm, nhiệt độ, độ ẩm, số ngày ẩm ướt

Nội dung chính theo từng chương

Chương 1. Mở đầu

Chương 1 trình bày tổng quan về tính cấp thiết của nghiên cứu xói mòn đất đối với sản xuất nông nghiệp, an ninh lương thực và môi trường sinh thái tại Việt Nam. Tác giả nêu bật những hạn chế hiện hữu của mô hình RUSLE truyền thống khi áp dụng vào thực tế địa phương, xác lập 3 mục tiêu nghiên cứu cụ thể, phạm vi không gian đất liền Việt Nam và khoảng thời gian thực hiện đề tài (10/03/2023 – 15/06/2023).

Chương 2. Các công trình liên quan

Tác giả phân tích các nghiên cứu điển hình trong và ngoài nước về ứng dụng mô hình xói mòn và GIS:

  • Nghiên cứu ứng dụng RUSLE2 tại lưu vực Suối Rạt, Bình Phước: Phân tích xói mòn - bồi tụ cho 56 mặt cắt sườn dốc; kết quả theo TCVN 5299:2009 chỉ ra xói mòn cấp IV chiếm 75.24% diện tích, cấp V chiếm 24.76%, cấp I–III chiếm tỷ lệ không đáng kể.
  • Nghiên cứu tích hợp RUSLE, Remote Sensing và GIS tại vùng El-Minya (Ai Cập): Đánh giá xói mòn vùng khô hạn, phân chia 4 cấp độ xói mòn và sử dụng chỉ số SWCR để chứng minh tác động của con người vượt trội biến đổi khí hậu.
  • Nghiên cứu so sánh USLE, MUSLE và RUSLE tại lưu vực Wadi Ghazouana (Algeria): Tốc độ xói mòn ước tính lần lượt đạt 9.65 t/ha/năm (USLE), 9.90 t/ha/năm (RUSLE) và 11.33 t/ha/năm (MUSLE), cho thấy MUSLE tính toán lượng mưa chi tiết hơn nhưng làm chậm tốc độ xử lý.

Từ đó, chương 2 tổng hợp các thách thức về quy mô không gian lớn, độ phức tạp địa hình và tốc độ xử lý, đề xuất phương hướng xây dựng mô hình tùy biến trên nền tảng ArcGIS.

Chương 3. Nghiên cứu và thực nghiệm

Nội dung chương tập trung vào cơ sở lý thuyết công nghệ và tiền xử lý các bộ dữ liệu:

  • Cơ chế dữ liệu ArcGIS: So sánh mô hình Raster và Vector. Phân tích vai trò của ArcSDE, ArcInfo và cấu trúc GeoDatabase trong quản trị dữ liệu không gian.
  • Tương đồng địa hình - khí hậu giữa Ethiopia và Việt Nam: Luận giải việc áp dụng các công thức thực nghiệm từ Ethiopia sang Việt Nam dựa trên sự tương đồng về hệ thống núi cao phức tạp (Hoàng Liên Sơn, Trường Sơn so với Simien, Bale), các cao nguyên lớn (Tây Nguyên so với các cao nguyên Ethiopia), các vùng đồng bằng châu thổ rộng lớn (Đồng bằng sông Hồng, Đồng bằng sông Cửu Long so với Rift Valley), cùng chế độ khí hậu nhiệt đới gió mùa phân hóa mùa mưa - mùa khô sâu sắc.
  • Tiền xử lý 4 bộ dữ liệu nền tảng:
    • GLC-SHARE: Phân tích tỷ lệ các lớp phủ đất, gồm Diện tích cây bao phủ (27.7%), Đất trần Bareland (15.0%), Đồng cỏ Grassland (13.0%), Đất trồng trọt Cropland (12.6%), Cây bụi Shrubs (9.5%), Tuyết/sông băng (9.7%), Thảm thực vật mỏng (7.0%), Vùng nước Water bodies (2.6%), Thực vật thân thảo (1.3%), Bề mặt nhân tạo (0.6%), Rừng ngập mặn Mangroves (0.7%).
    • CRU-TS: Trích xuất 13 thông số khí tượng (Pre, Tmp, Tmn, Tmx, Pet, Cld, Dtr, Frs, Rhm, Ssh, Vap, Wet, Wnd).
    • HydroSHEDS: Trích xuất lưới tích lũy dòng chảy ACC và dữ liệu độ dốc DTED từ radar SRTM.
  • Lý thuyết các hệ số và đề xuất mô hình RUSLE tùy biến: Phân tích khả năng tích hợp ảnh hưởng của lượng mưa ($R$) vào hệ số địa hình dòng chảy ($LS$) và hệ số canh tác ($P$). Đề xuất thử nghiệm 4 tổ hợp factor chính ($K \times LS$, $K \times LS \times P$, $K \times LS \times C$, $K \times LS \times C \times P$) cùng tổ hợp tinh gọn $K \times C$.
Lớp phủ bề mặt (GLC-SHARE) Tỷ lệ diện tích (%) Đặc điểm bảo vệ đất chống xói mòn
Vùng cây bao phủ (Tree Covered Area) 27.7% Che phủ lớn nhất, giữ ẩm đất, hấp thụ carbon, chống xói mòn cao nhất
Đất trần (Bareland) 15.0% Không có thực vật bao phủ, dễ bị rửa trôi và xói mòn nghiêm trọng nhất
Đồng cỏ (Grassland) 13.0% Hệ sinh thái thảo mộc rộng lớn, bảo vệ mặt đất ở mức trung bình
Đất trồng trọt (Cropland) 12.6% Đất canh tác nông nghiệp, mức độ xói mòn phụ thuộc biện pháp mùa vụ
Tuyết và sông băng (Snow & Glacier) 9.7% Tập trung ở vùng cực / núi cao quanh năm
Cây bụi (Shrubs Covered Area) 9.5% Thân thấp gần mặt đất, nhiều cành nhánh, chống xói mòn hiệu quả
Thảm thực vật mỏng (Sparse Vegetation) 7.0% Mật độ thực vật thấp do môi trường khắc nghiệt
Vùng nước (Water bodies) 2.6% Sông, hồ, kênh rạch chứa nước
Thực vật thân thảo (Herbaceous Vegetation) 1.3% Cây thân mềm, thảo dược ngắn ngày
Rừng ngập mặn (Mangroves) 0.7% Ven biển, cửa sông, đầm lầy, chống sạt lở bờ biển
Bề mặt nhân tạo (Artificial Surfaces) 0.6% Đô thị, đường sá, công trình xây dựng, khu công nghiệp

Chương 4. Thiết kế và lưu trữ dữ liệu địa lý

Chương 4 mô tả quy trình hiện thực hóa dữ liệu trên GeoDatabase:

  • Tổ chức GeoDatabase: Thiết kế hệ thống bảng dữ liệu không gian (Spatial Data Tables) lưu trữ tọa độ, hình học đối tượng (Geometry/Shape) của các lớp polygon phân loại đất và bảng dữ liệu phi không gian (Nonspatial Data Tables) quản lý các thuộc tính cơ lý hóa của đất.
  • Tính toán chi tiết các Factor trên lãnh thổ Việt Nam:
    • Tính $K$-factor: Áp dụng công thức giải tích cho từng loại đất từ FAO-DSMW, lập Scatter Plot CDOMSOIL - Shape và bản đồ phân loại đất liền.
    • Tính $R$-factor: Xử lý dữ liệu mưa CRU-TS giai đoạn 1991–2000 qua mô hình Hurni.
    • Tính $LS$-factor & $P$-factor: Xử lý DEM, làm mượt dữ liệu (DEM smoothing), tính toán ma trận dòng chảy tích lũy tại khu vực Dãy Trường Sơn, Tây Bắc, Đồng bằng sông Hồng và Đồng bằng sông Cửu Long.
    • Tính $C$-factor: Gán trọng số che phủ từ GLC-SHARE.
  • Xây dựng Custom RUSLE Model: Thực hiện kết hợp (combine) các lớp raster theo từng tổ hợp factor, phân tích độ xói mòn tại các lưu vực sông lớn (sông Hồng, sông Cửu Long) và các hồ chứa (Hồ Trị An, Hồ Dầu Tiếng).

Chương 5. Kết luận và hướng phát triển

Đánh giá tổng kết các kết quả đạt được về mặt học thuật và kỹ thuật mô hình hóa, đồng thời vạch rõ các hạn chế về dữ liệu và đề xuất định hướng phát triển trong tương lai.


Kết quả và đóng góp

Kết quả thực nghiệm chính

Quá trình thử nghiệm 5 tổ hợp mô hình trên toàn bộ dữ liệu đất liền Việt Nam mang lại các kết quả định lượng cụ thể:

Tên mô hình tổ hợp Các Factor thành phần Mức độ xói mòn tính toán Thời gian xử lý (kết hợp mô hình) Khả năng bao hàm vùng sông ngòi / bồi tụ
KLS $K, LS$ Phản ánh xói mòn địa hình thô Chậm hơn Kém, không thể hiện rõ bồi tụ
KLSP $K, LS, P$ Điều chỉnh theo tập quán canh tác Trung bình Hạn chế ở vùng hạ lưu
KLSC $K, LS, C$ Kết hợp độ che phủ thực vật Trung bình Thể hiện tốt vùng đồi núi
KLSPC $K, LS, P, C$ Mô phỏng đầy đủ 4 yếu tố Tương đối chậm Khá, nhưng tính toán cồng kềnh
KC Model (Tối ưu) $K, C$ 13.7 đến 674 tấn/ha/năm 1.32 giây Tốt nhất, bao hàm xói mòn và bồi tụ lưu vực
  • Mô hình KC cho kết quả tối ưu nhất với dải giá trị mất đất dao động từ 13.7 đến 674 tấn/ha/năm, tốc độ xử lý nhanh nhất đạt 1.32 giây.
  • Mô hình KC thể hiện năng lực vượt trội trong việc bao hàm và phản ánh đúng thực tế mức độ xói mòn, bồi tụ tại các khu vực phức tạp như hệ thống kênh rạch, sông ngòi và vùng đồng bằng ngập nước (Đồng bằng sông Cửu Long, Đồng bằng sông Hồng).

Giải pháp và kiến nghị

  • Sử dụng mô hình tổ hợp tùy biến tinh gọn ($KC$) cho các bài toán phân tích nhanh trên quy mô diện rộng nhằm giảm tải tính toán mà vẫn đảm bảo độ bao quát hiện tượng xói mòn - bồi tụ.
  • Khuyến nghị áp dụng đồng bộ các giải pháp quản lý đất đai dựa trên kết quả phân tích: tăng cường độ che phủ của rừng và thảm thực vật tại các sườn dốc vùng Tây Bắc và dãy Trường Sơn; duy trì hệ thống rừng ngập mặn và đê điều tại Đồng bằng sông Cửu Long để hạn chế xâm nhập mặn và sạt lở bờ sông, bờ biển.

Đóng góp của đề tài

  • Xây dựng thành công cơ sở dữ liệu địa lý GeoDatabase hoàn chỉnh tích hợp đa nguồn dữ liệu toàn cầu (FAO-DSMW, GLC-SHARE, HydroSHEDS, CRU-TS) phục vụ nghiên cứu xói mòn đất tại Việt Nam.
  • Khắc phục các nhược điểm cố hữu của mô hình RUSLE cũ về sai số vùng bồi tụ và thời gian xử lý chậm, cung cấp bộ tham số thực nghiệm phù hợp với hai dạng địa hình chính của Việt Nam.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của đề tài

  • Quá trình thu thập dữ liệu trên phạm vi toàn quốc gặp phải tình trạng dữ liệu thô bị trống (null) hoặc thiếu chính xác cục bộ ở một vài khu vực do độ phân giải của các bộ dữ liệu toàn cầu còn thô (như CRU-TS có độ phân giải 50 km).
  • Phạm vi nghiên cứu mới chỉ tập trung trên phần lãnh thổ đất liền, chưa xử lý và tích hợp đồng bộ dữ liệu không gian khu vực biển đảo và các hải đảo chủ quyền của Việt Nam.
  • Hệ thống đòi hỏi người dùng cuối phải có kiến thức nền tảng nhất định về GIS và phần mềm ArcGIS để vận hành và khai thác cơ sở dữ liệu.

Hướng phát triển tiếp theo

  • Triển khai (Deploy) cơ sở dữ liệu và mô hình lên hệ thống thông tin địa lý không gian ba chiều (3D GIS) để tăng tính trực quan và giúp người dùng phổ thông dễ dàng tiếp cận.
  • Mở rộng phạm vi mô hình hóa ra toàn bộ vùng biển đảo và hải đảo thuộc chủ quyền Việt Nam sau khi hoàn thiện giai đoạn triển khai.
  • Cập nhật các bộ dữ liệu viễn thám có độ phân giải cao hơn nhằm nâng cao độ chính xác tại các vi lưu vực có địa hình chia cắt mạnh.

Giá trị tham khảo

  • Đối tượng khai thác tài liệu:
    • Sinh viên, học viên chuyên ngành Công nghệ Thông tin, Hệ thống Thông tin Địa lý (GIS), Khoa học Môi trường, Quản lý Đất đai, Địa lý và Địa chất: Dùng làm tài liệu tham khảo về quy trình thu thập, chuyển đổi định dạng và tích hợp các nguồn dữ liệu không gian toàn cầu vào bài toán cụ thể tại Việt Nam.
    • Cán bộ nghiên cứu và quản lý nông nghiệp, môi trường: Tham khảo phương pháp luận tùy biến mô hình RUSLE và các bảng trọng số hệ số đất ($K$), độ phủ ($C$) phục vụ công tác quy hoạch canh tác và chống suy thoái đất.
  • Nội dung có giá trị tham khảo cao nhất:
    • Quy trình thiết kế, tổ chức cấu trúc bảng dữ liệu không gian và phi không gian trong GeoDatabase trên phần mềm ArcGIS Pro.
    • Phương pháp luận chuyển đổi và tính toán giải tích hệ số xói mòn đất ($K$-factor) theo Goldman et al. và công thức mưa ($R$-factor) theo Hurni áp dụng cho điều kiện Việt Nam.

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

1. Tại sao nghiên cứu lại lược bỏ thành phần $R_g$ trong công thức tính $R$-factor của Hurni (1985)?

Trong công thức gốc của Hurni (1985), chỉ số phụ $R_g$ được đưa vào để tính toán tác động xói mòn bổ sung do tuyết tan và nhiệt tan. Do đặc thù khí hậu Việt Nam là nhiệt đới gió mùa, không có hiện tượng tuyết phủ và băng tan diện rộng gây xói mòn, nên thành phần $R_g$ được lược bỏ nhằm tinh giản mô hình, công thức tính $R$-factor được rút gọn thành $R = 0.12 \times P$.

2. Mô hình tổ hợp nào cho kết quả tốt nhất trong các thực nghiệm và vì sao?

Mô hình tổ hợp $KC$ ($A = K \times C$) đạt kết quả tốt nhất với tốc độ tính toán nhanh nhất (1.32 giây), dải xói mòn mô phỏng đạt từ 13.7 đến 674 tấn/ha/năm và phản ánh tốt cả hiện tượng xói mòn lẫn bồi tụ tại khu vực sông ngòi, kênh rạch. Điều này đạt được do yếu tố mưa và dòng chảy đã được điều tiết đáng kể thông qua độ che phủ ($C$) và tính chất thấm giữ nước của đất ($K$).

3. Bộ dữ liệu khí hậu CRU-TS cung cấp những thông số nào và có độ phân giải bao nhiêu?

Bộ dữ liệu CRU-TS cung cấp chuỗi số liệu từ năm 1901 đến 2022 với độ phân giải không gian $0.5^\circ$ (~50 km), bao gồm 13 biến khí tượng: độ che phủ mây (Cld), khoảng nhiệt độ ngày (Dtr), tần suất ngày sương giá (Frs), khả năng thoát hơi nước (Pet), lượng mưa (Pre), độ ẩm tương đối (Rhm), thời lượng nắng (Ssh), nhiệt độ trung bình tháng (Tmp), nhiệt độ trung bình thấp nhất (Tmn), nhiệt độ trung bình cao nhất (Tmx), áp suất hơi (Vap), tần suất ngày ẩm ướt (Wet) và tốc độ gió (Wnd).

4. Tại sao tác giả lại so sánh đặc điểm địa hình và khí hậu của Việt Nam với Ethiopia khi thiết lập mô hình?

Tác giả so sánh với Ethiopia vì hai quốc gia có nhiều nét tương đồng về cấu trúc địa hình và khí hậu: cả hai đều có hệ thống dãy núi dài phức tạp xen lẫn các cao nguyên lớn và đồng bằng châu thổ rộng; đồng thời đều thuộc đới khí hậu nhiệt đới/cận nhiệt đới có sự phân chia mùa mưa - mùa khô rõ rệt. Sự tương đồng này cho phép kế thừa hợp lý các công thức thực nghiệm của Hurni (1985) đã được kiểm chứng tại Ethiopia.

5. Khóa luận sử dụng các lớp che phủ nào từ bộ dữ liệu GLC-SHARE và lớp nào chiếm tỷ lệ lớn nhất?

GLC-SHARE cung cấp 11 lớp phủ mặt đất: Bề mặt nhân tạo (0.6%), Đất trồng trọt (12.6%), Đồng cỏ (13.0%), Vùng cây bao phủ (27.7%), Cây bụi (9.5%), Thực vật thân thảo (1.3%), Rừng ngập mặn (0.7%), Thảm thực vật mỏng (7.0%), Đất trần (15.0%), Tuyết/sông băng (9.7%), và Vùng nước (2.6%). Trong đó, lớp Vùng cây bao phủ (Tree Covered Area) chiếm tỷ lệ lớn nhất với 27.7%.


Kết luận

Khóa luận đã nghiên cứu và ứng dụng thành công hệ thống thông tin địa lý ArcGIS cùng mô hình RUSLE tùy biến để phân loại đất và tính toán định lượng xói mòn đất trên quy mô toàn bộ đất liền Việt Nam. Bằng việc tích hợp đa nguồn dữ liệu toàn cầu chuẩn hóa (FAO-DSMW, GLC-SHARE, HydroSHEDS, CRU-TS) vào cơ sở dữ liệu GeoDatabase, tác giả đã chứng minh mô hình tổ hợp $KC$ mang lại hiệu năng xử lý tối ưu (1.32 giây) và giải quyết tốt bài toán mô phỏng xói mòn - bồi tụ tại các lưu vực sông ngòi. Đây là công trình có giá trị học thuật và thực tiễn cao, đóng góp cơ sở dữ liệu và phương pháp luận quan trọng cho công tác quản lý tài nguyên đất và quy hoạch nông nghiệp bền vững.