Đánh giá khả năng xói mòn đất ở huyện ĐakRong, tỉnh Quảng Trị mô hình RMMF

Header of 133 TẠPPage CHÍ1KHOA HỌC, Đại học Huế, Tập 74A, Số 5, (2012), 173-184 ĐÁNH GIÁ KHẢ NĂNG XÓI MÒN ĐẤT Ở HUYỆN ĐAKRÔNG, TỈNH QUẢNG TRỊ BẰNG MÔ HÌNH RMMF (REVISED MORGAN-MORGAN-FINNEY) Trương Đình Trọng, Nguyễn Quang Việt, Đỗ Thị Việt Hương Trường Đại học Khoa học, Đại học Huế Tóm tắt Đakrông huyện miền núi tỉnh Quảng Trị có địa hình phân hóa phức tạp; lượng mưa lớn tập trung; thảm phủ thực vật nghèo dần khai thác rừng tập quán đốt nương làm rẫy người dân khu vực Do đó, khả đất bị xói mòn xảy lớn Với trợ giúp công nghệ GIS, tác giả sử dụng mô hình RMMF để tính toán khả xói mòn đất Dữ liệu đầu vào mô hình gồm 15 thông số liên quan đến địa hình, thổ nhưỡng, khí hậu thảm phủ Kết tính toán cho thấy lượng đất xói mòn biến thiên từ 0-957(tấn/ha.năm) chia thành cấp xói mòn Xói mòn yếu chiếm phần lớn diện tích (47,28%), trung bình chiếm 5,43%, cấp xói mòn mạnh chiếm diện tích nhỏ Qua việc so sánh mối tương quan lượng đất xói mòn với nhân tố ảnh hưởng, nhận thấy bên cạnh nhân tố địa hình, đặc điểm thổ nhưỡng khí hậu thảm phủ thực vật đóng vai trò lớn, định đến lượng đất xói mòn lãnh thổ nghiên cứu Vì vậy, việc bảo vệ lớp phủ thực vật khu vực địa hình dốc biện pháp hữu hiệu để hạn chế xói mòn đất Từ khóa: xói mòn đất, xói mòn, RMMF, mô hình RMMF, Đakrông Đặt vấn đề Hiện nay, tác động trình tự nhiên người ngày mạnh nên trình thoái hóa đất diễn ngày trầm trọng, đe dọa đến nông nghiệp nhiều nước giới Việt Nam Trong đó, sa mạc hóa xói mòn đất xem hai trình ảnh hưởng đến suy thoái tài nguyên đất Theo kết công bố H.Eswaran cộng (2001), khả sản xuất số khu vực giới bị giảm đến 50% xói mòn sa mạc hóa; Ở khu vực Nam Á sản xuất ngũ cốc giảm khoảng 36 triệu tấn/năm xói mòn nước tương đương với 5,4 tỉ USD; Ở mức độ toàn cầu, hàng năm bề mặt trái đất 75 tỉ đất tương đương với giá trị kinh tế 400 tỉ USD, tính theo đầu người công dân trái đất năm thiệt hại xấp xỉ 70 USD Huyện Đakrông có địa hình phần lớn đồi núi; khí hậu mang tính chuyển tiếp Đông Tây Trường Sơn với lượng mưa trung bình năm 2000 mm; thảm thực vật tự nhiên bị giảm sút khai thác rừng tập quán đốt nương làm rẫy Vì đất có khả bị xói mòn lớn Do đó, việc đánh giá, dự báo định lượng xói mòn 173 Footer Page of 133 Header 174 Page of 133 Đánh giá khả xói mòn đất huyện Đakrông… làm sở khoa học cho việc lập kế hoạch quản lý, sử dụng tài nguyên đất cách hữu hiệu cần thiết Hiện nay, số công trình tiến hành nghiên cứu để xác định nguy xói mòn đất bao gồm nghiên cứu định tính lẫn định lượng Trong đó, mô hình RMMF [1] phát triển gần dùng cho tính toán dự báo xói mòn áp dụng cho nhiều vùng có điều kiện tự nhiên khác giới cho kết đánh giá tin cậy Nội dung nghiên cứu 2.1 Sử dụng công cụ GIS cho mô hình RMMF Cơ sở liệu cần thiết cho mô hình RMMF bao gồm: - Bản đồ Thổ nhưỡng tỷ lệ 1/50.000 số liệu phân tích đất tỉnh Quảng Trị - Bản đồ Hiện trạng thảm thực vật tỉnh Quảng Trị tỷ lệ 1/50.000 - Bản đồ Địa hình tỉnh Quảng Trị tỷ lệ 1/50.000 - Bản đồ Phân vùng khí hậu tỉnh Quảng Trị tỷ lệ 1/50.000 - Tài liệu khí hậu số liệu quan trắc trạm khí tượng tỉnh Quảng Trị Tất liệu đồ, số liệu phân tích mã hóa nhờ trợ giúp công cụ GIS theo 03 bước: (1) số hóa tạo sở liệu, (2) chuyển liệu sang ArcGIS raster hóa (3) chạy công thức tính toán mô hình Dữ liệu mưa bao gồm: Lượng mưa trung bình năm số hóa từ đồ phân bố mưa năm tỉnh Quảng Trị Mapinfo chuyển sang định dạng ArcGIS; số ngày mưa dựa vào trạm quan trắc: Trạm Đông Hà (thuộc Đông Trường Sơn), trạm Khe Sanh (thuộc Tây Trường Sơn) Trên sở số liệu quan trắc trạm đồ khí hậu tỉnh để thành lập đồ phân hóa số ngày mưa Độ dốc nội suy từ mô hình số độ cao (DEM) công cụ Interpolation ArcGIS 9.3 Trong đó, DEM xây dựng từ điểm độ cao đồ địa hình Tất liệu đồ chuyển sang định dạng ArcGIS từ liệu Mapinfo, sau raster hóa thành đồ thành phần (với kích thước pixel 15m) 2.2 Ứng dụng mô hình RMMF đánh giá xói mòn đất Mô hình RMMF dùng để tính lượng đất xói mòn cho khu vực đồi núi Mô hình RMMF đòi hỏi 15 thông số đầu vào thể bảng Footer Page of 133 Header Page ofTRƯƠNG 133 ĐÌNH TRỌNG, NGUYỄN QUANG VIỆT, ĐỖ THỊ VIỆT HƯƠNG 175 Bảng Các thông số đầu vào mô hình RMMF (Morgan 2000) [7] Nhân tố xói mòn Thông số tính toán Mưa Thổ nhưỡng Giải thích thông số Nguồn liệu xác định giá trị cho thông số R Lượng mưa TB năm (mm) Số liệu khí tượng/Bản đồ phân bố mưa năm Rn Số ngày mưa năm Khí hậu Quảng Trị I Cường độ mưa (mm/h) Tài liệu tham khảo MS Khả trữ ẩm đất (%w/w) Tài liệu theo thành phần giới BD Dung trọng lớp đất bề mặt Tài liệu theo thành phần (mg/m3) giới EHD Độ sâu thủy học đất (m) Tài liệu tham khảo theo thảm phủ K Chỉ số xói mòn đất (g/J) Số liệu phân tích Bản đồ đất Tính dính kết bề mặt (kPa) Tài liệu, theo loại đất S Độ dốc (0) Nội suy từ DEM (15m) A Tỷ lệ (0 - 1) lượng mưa bị cản Tài liệu tham khảo theo thảm phủ thảm phủ Et/Eo Tỷ lệ bốc thoát nước thực tế Tài liệu tham khảo theo (Et) tiềm (Eo) thảm phủ C Nhân tố quản lý thảm phủ (kết Tài liệu tham khảo theo hợp C P USLE) thảm phủ COH Địa hình Thảm phủ Footer Page of 133 CC Độ che phủ tán (0 - 1) Quan trắc Bản đồ trạng thảm phủ GC Độ che phủ bề mặt đất (0 - 1) Quan trắc Bản đồ trạng thảm phủ PH Chiều cao thảm phủ Quan trắc Bản đồ trạng thảm phủ Header 176 Page of 133 Đánh giá khả xói mòn đất huyện Đakrông… Thảm phủ Mưa Thổ nhưỡng Các thông số thảm phủ Địa hình Các thông số thổ nhưỡng Mưa hữu hiệu (ER) Mưa qua tán (LD) DEM Mưa trực tiếp (DT) Năng lượng KE (LD) Năng lượng KE (DT) Tổng lượng mưa (KE) Tổng dòng chảy (Q) Lượng đất tách mưa (F) Lượng đất tách dòng chảy (H) Tổng lượng đất tách (J) Khả vận chuyển (G) Bản đồ độ dốc Min (J;G) Lượng đất Hình Sơ đồ minh hoạ ứng dụng mô hình RMMF đánh giá xói mòn đất [2] Bảng Giá trị thông số đầu vào thuộc nhóm thảm phủ mô hình RMMF CC (%) GC (%) Rừng rậm thường xanh bị 16.00 tác động 85 Rừng rậm thường xanh bị 12.00 tác động mạnh Trảng bụi thứ sinh có Loại thảm phủ Footer Page of 133 PH (m) 4.50 Et/Eo A (%) EHD (m) C 98 0.95 30 0.20 0.007 70 80 0.90 27 0.20 0.090 65 72 0.90 25 0.20 0.600 Header Page ofTRƯƠNG 133 ĐÌNH TRỌNG, NGUYỄN QUANG VIỆT, ĐỖ THỊ VIỆT HƯƠNG 177 gỗ rải rác Trảng bụi thứ sinh gỗ 3.50 60 68 0.90 23 0.20 0.700 Trảng cỏ thứ sinh 0.70 35 70 0.86 20 0.14 0.830 Rừng trồng 4.50 50 45 0.75 25 0.20 0.150 Cây lâu năm 4.20 35 30 0.70 21 0.20 0.170 Cây hàng năm (ngô, sắn) 1.10 25 25 0.68 20 0.12 0.280 Nương rẫy 1.00 20 20 0.68 17 0.12 0.600 Lúa 0.60 50 50 1.35 20 0.12 0.060 Ghi chú: A, CC, GC, PH đo đạc thực địa; giá trị lại tham khảo giá trị từ tài liệu trích dẫn [1], [7], [8] * Tính toán lượng mưa: - Lượng mưa hữu hiệu: ER = R*(1 - A) (1) ER chia thành phần: Lượng mưa trực tiếp (DT), phần qua tán đến bề mặt (LD): LD = ER*CC (2) DT = ER – LD (3) Năng lượng động học DT (KE(DT); J/m2 ) (Kinetic Energy of Direct Throughfall) xác định hàm cường độ mưa (I; mm/h) Công thức sau phát triển Hudson (1965) Zimbabwe có khí hậu nhiệt đới [8]: KE(DT) = DT* (29.8 - (127.5/I)) (4) Năng lượng động học LD (KE(LD); J/m2) (Kinetic Energy of Leaf Drainage) phụ thuộc vào chiều cao tán (PH; m): KE(LD) = LD* ((15.8*PH0.5) - 5.87) (5) Tổng lượng động học lượng mưa hữu hiệu (KE;J/m2): KE = KE(DT) + KE(LD) (6) * Tính toán dòng chảy mặt: Công thức tính toán dòng chảy đề xuất Kirkby (1976), dòng chảy mặt phát sinh lượng mưa ngày vượt khả trữ ẩm đất (Rc; mm) Q = R*exp(-Rc/Ro) (7) Trong đó: Ro = R/Rn; Ro lượng mưa trung bình ngày (mm) Khả trữ ẩm đất tính toán: Rc = 1000*MS*BD*EHD*(Et/Eo) (8) Footer Page of 133 Header 178 Page of 133 Đánh giá khả xói mòn đất huyện Đakrông… Bảng Hệ số K Bảng Các giá trị hướng dẫn thổ nhưỡng cho RMMF (Morgan 2000; Morgan Duzant 2007) [7] Ghi chú: Hệ số K xác định phân tích loại đất tham khảo số liệu [5] *Tính toán lượng đất tách ra: Công thức mô hình RMMF để ước tính khả tách phần tử đất ảnh hưởng hạt mưa: F = K*KE*10-3 (9) Hệ số ứng chịu xói mòn đất xác định dựa vào mẫu đất phân tích thành phần cấp hạt, độ mùn, tính thấm cấu trúc đất Từ đó, xác định K dựa vào toán đồ USDA Công thức ước tính lượng đất tách dòng chảy mặt dựa vào công thức thực nghiệm Quansah (1982): H = Z*Q1.5* sinS*(1 - GC)*10-3 (10) Trong đó, kháng cự đất Z = 1/0.5*COH (11) Tổng lượng đất tách tính: J = F + H (12) * Tính toán khả vận chuyển dòng chảy: G = C*Q2*sinS*10-3 (S tính độ) (13) * Tính toán lượng đất xói mòn: Tính toán tổng lượng đất tách tác động mưa dòng chảy, sau so sánh với khả vận chuyển dòng chảy mặt Giá trị lượng đất xói mòn năm (Meyer Wischmeier, 1969) Lượng đất xói mòn = Min (J;G) (14) Footer Page of 133 Header Page ofTRƯƠNG 133 ĐÌNH TRỌNG, NGUYỄN QUANG VIỆT, ĐỖ THỊ VIỆT HƯƠNG 179 2.3 Kết đánh giá khả xói mòn đất mô hình RMMF 2.3.1 Xác định thông số cho mô hình a Năng lượng công kích mưa (KE) Năng lượng mưa lực hạt mưa tác động lên bề mặt đất Trong mô hình RMMF, KE xác định R, A, CC PH Kết tính toán KE theo công thức (6) thể hình b Lưu lượng dòng chảy mặt (Q) Dòng chảy mặt phát sinh mưa độ trữ ẩm khả thấm bị bão hòa Q phụ thuộc vào MS, BD, EHDvà Et/Eo MS BD định loại đất, EHD Et/Eo quan hệ chặt chẽ với thảm phủ thực vật Kết tính toán theo (7) thể hình Hình Sơ đồ lượng mưa Hình Lưu lượng dòng chảy mặt c Tính toán thông số tách phần tử đất * Tác động hạt mưa (F): Các phần tử đất tách mưa phụ thuộc vào K KE Chỉ số K thể kháng cự trình tách vận chuyển phần tử đất Giá trị K cao dễ dàng bị tách Kết tính toán theo (9) trình bày hình * Tác động dòng chảy mặt (H): Dòng chảy mặt tách phần tử đất vận tốc di chuyển đủ lớn để tách Mô hình sử dụng (10) để xác định lượng đất tách bao gồm: COH, Q, GC độ dốc Kết tính toán thể hình Footer Page of 133 Header 180 Page of 133 Đánh giá khả xói mòn đất huyện Đakrông… Hình Khả tách mưa Hình Khả tách dòng chảy mặt Hình Tổng lượng đất bị tách Hình Khả vận chuyển dòng chảy * Tổng lượng đất tách (J): Tổng lượng đất tách kết lượng đất tách mưa (F) dòng chảy mặt (H) (Hình 6) Theo hình 6, J có giá trị dao động từ 4,5 đến 95,7 kg/m2 Trong đó, diện tích Footer Page of 133 Header Page ofTRƯƠNG 133 ĐÌNH TRỌNG, NGUYỄN QUANG VIỆT, ĐỖ THỊ VIỆT HƯƠNG 181 loại hàng năm, nương rẫy thảm bụi có giá trị cao; mặt khác số khu vực diện tích rừng lượng mưa qua tán có giá trị cao nên khả tách phần tử đất lớn c Khả vận chuyển dòng chảy mặt (G) Sức tải dòng chảy mặt khả dòng chảy vận chuyển vật liệu từ nơi đến nơi khác G phụ thuộc vào C, Q độ dốc Kết tính toán theo (13) thể hình Đại lượng G biến thiên khoảng giá trị từ đến 4.471,21 kg/m2 có phân hóa rõ loại thảm phủ Khả vận chuyển cao tập trung chủ yếu diện tích hàng năm, lâu năm bụi thứ sinh có hệ số thảm phủ C cao Rừng tự nhiên có hệ số thảm phủ C thấp sinh khối cao nên hạn chế nhiều khả vận chuyển vật liệu đất bề mặt d Tính toán lượng đất xói mòn: Theo tính toán, lượng đất xói mòn dao động từ đến 957 tấn/ha/năm, chủ yếu phổ biến dãy giá trị từ đến 100 tấn/ha/năm, pixel có giá trị 100 tấn/ha/năm không đáng kể Do đó, lượng đất trung bình toàn lãnh thổ 13,42 tấn/ha/năm, tổng lượng đất xói mòn hàng năm vào khoảng 1,67 triệu Phần lớn diện tích lãnh thổ có lượng đất xói mòn thấp, khu vực xói mòn cao chủ yếu tập trung dọc thung lũng sông Quảng Trị Đường (từ xã Ba Lòng đến xã Hướng Hiệp); dọc đường Hồ Chí Minh nơi có độ che phủ thấp, chủ yếu hàng năm, nương rẫy bụi thứ sinh có độ che phủ thấp Những diện tích che phủ rừng tự nhiên có lượng xói mòn thấp 2.3.2 Phân cấp xói mòn Hiện nay, có nhiều thang phân cấp xói mòn khác cho vùng, chưa có số liệu cụ thể phân cấp xói mòn chung cho tất vùng Do đó, tác giả dựa vào đặc điểm lãnh thổ nghiên cứu tạm thời chia thành cấp xói mòn sau: Bảng Phân hạng cấp xói mòn đất huyện Đakrông Lượng đất (tấn/ha/năm) Diện tích (ha) Tỷ lệ % < 20 57.593,05 47,28% 20 – 40 6.057,77 5,43% Cấp – Mạnh 40.1 – 60 5.999,74 5,88% Cấp – Rất mạnh 60.1 – 80 8.356,22 9,10% > 80 43.795,41 35,96% Cấp xói mòn Cấp – Yếu Cấp – Trung bình Cấp – Nguy hiểm So sánh kết phân cấp với đồ thành phần (các hệ số), nhận thấy phần lớn lãnh thổ có khả xói mòn yếu, điều phù hợp với thảm phủ rừng tương đối Footer Page of 133 Header 182 Page 10 of 133 Đánh giá khả xói mòn đất huyện Đakrông… nhiều hạn chế ảnh hưởng mưa độ dốc địa hình Xói mòn trung bình chiếm tỉ lệ 5,43%, nơi nằm khu vực có độ dốc lượng mưa lớn, thảm phủ chủ yếu rừng trồng, thảm bụi thứ sinh Các cấp xói mạnh mạnh chiếm diện tích gần 15%, chủ yếu xảy khu vực canh tác nương rẫy trồng loại (ngô, sắn) có độ che phủ bề mặt thấp phân bố phần diện tích dốc Cấp nguy hiểm chiếm diện tích lớn 35,96% tập trung hàng năm nương rẫy khu vực đất dốc, đất có hệ số xói mòn cao Hình Lượng đất xói mòn Hình Phân cấp xói mòn đất Kết luận 3.1 Mô hình RMMF công cụ hữu hiệu để tính toán xói mòn khu vực đồi núi Mô hình thiết kế dựa vào động lực gây xói mòn, thông số đòi hỏi phải nội suy quan trắc thực địa cách tỉ mỉ Các yếu tố đầu vào tính toán dễ dàng nhờ trợ giúp GIS Mối quan hệ lượng đất xói mòn nhân tố ảnh hưởng nhận thấy cách dễ dàng 3.2 So sánh mối tương quan lượng đất với nhân tố gây xói mòn cho thấy bên cạnh độ dốc lượng mưa thảm phủ thực vật đóng vai trò định lớn đến lượng đất xói mòn lãnh thổ nghiên cứu 3.3 Trong trình xói mòn, khả tách phần tử đất mưa đóng vai trò chủ đạo so với dòng chảy Điều phù hợp với lượng mưa lớn lớp phủ thực vật có tầng tán bề mặt đất che phủ tốt 3.4 Kết đánh giá khả xói mòn đất mô hình RMMF lãnh thổ huyện Đakrông cho thấy diện tích cấp xói mòn có phân hóa lớn Mức độ xói Footer Page 10 of 133 Header Page 11 TRƯƠNG of 133 ĐÌNH TRỌNG, NGUYỄN QUANG VIỆT, ĐỖ THỊ VIỆT HƯƠNG 183 mòn yếu nguy hiểm chiếm phần lớn diện tích tự nhiên, cấp xói mòn trung bình cấp xói mòn mạnh chiếm phần nhỏ diện tích bề mặt lãnh thổ 3.5 Phần lớn diện tích có lượng đất xói mòn lớn tập trung diện tích hàng năm, nương rẫy bụi thứ sinh Điều cho thấy, hoạt động nhân tác đóng vai trò lớn việc thúc đẩy khả xói mòn lãnh thổ nghiên cứu, cần phải có biện pháp sử dụng bảo vệ đất cách hợp lý TÀI LIỆU THAM KHẢO [1] Morgan R.P.C and Duzant J.H, Modified MMF (Morgan-Morgan-Finney) model for evaluating effects of crops and vegetation cover on soil erosion, Journal of Earth surface processes and Landfoms 32, (2008), 90-106 [2] Khatereh Polous, Effect of spatial resolution on erosion assessment in Namchun watershed, Thailand, Facuty of Geo-Information science and Earth observation university of Twente, Enschede, The Netherlands, 2010 [3] Sở Khoa học công nghệ Môi trường tỉnh Quảng Trị, Báo cáo chuyên đề thuyết minh đồ đất tỉnh Quảng Trị tỷ lệ 1:50.000, Đông Hà, 2000 [4] Sở Khoa học công nghệ Môi trường tỉnh Quảng Trị, Số liệu kết phân tích, tả tả phân tích đất, Đông Hà, 2000 [5] Sở Tài nguyên & Môi trường tỉnh Quảng Trị, Đánh giá tình hình xói lở bồi lắng dòng sông hệ thống sông Thạch Hãn tỉnh Quảng Trị, Quảng Trị, 2010 [6] Andi Sukman, Assessing erosion hazard using Revised Morgan Morgan and Finney (MMF) erosion model and microtopography features; A case study in river Oyo subcatchment, Master of Science in Geo-Information for Spatial Planning and Risk management, Gadjah Mada university, International Institute for Geo-Information Science and Earth observation, 2009 [7] Ugyen Thinley, Spatial Modeling for Soil erosion assessment in upper Lam Phra Phloeng watershed, Nakhon Ratchasima, Thailand, 2008 [8] Vũ Anh Tuấn, Nghiên cứu biến động trạng thảm thực vật ảnh hưởng tới trình xói mòn lưu vực sông Trà Khúc phương pháp viễn thám hệ thông tin địa lý, Luận án Tiến sỹ (Mã số: 1.07.14), Hà Nội, 2000 Footer Page 11 of 133 Header 184 Page 12 of 133 Đánh giá khả xói mòn đất huyện Đakrông… ASSESSMENT OF POTENTIAL SOIL EROSION BY RMMF MODEL (REVISED MORGAN – MORGAN – FINNEY) IN DAKRONG DISTRICT, QUANG TRI PROVINCE Truong Dinh Trong, Nguyen Quang Viet, Do Thi Viet Huong College of Sciences, Hue University Abstract Dakrong, a mountainous district in Quang Tri Province, is characterized by a diversified topography with a large variety of height, high rainfall, and deteriorating land cover due to exploitation of forest and the practice of burning vegetation to make land for cutivation Thus, there is a high possibility of soil erosion With the support of GIS technology, the authors used RMMF model to calculate potential soil erosion The input data of the model include 15 coeffcients relating to topography, soil properties, climate and land cover Simulations of RMMF include estimates of rainfall energy, runoff, soil particle detachment by raindrop impact, soil particle detachment by runoff, transport capacity of runoff and soil loss Results showed that the anual amount of soil loss range from a minimum of tons/ha to a maximum of 957 tons/ha and is divided into classes of erosion Weak class erosion covers the largest area of the region researched (47.28% of total area), moderate class occupies 5.43% of total area, while strong classes only covers a small area Through studying the relationships between soil loss mass and erosion factors, we found that vegetational cover plays a more significant role in determining soil loss than topography, climate and properties of soil Therefore, the protection of forest floor in sloping area is one of the most effective methods to reduce soil erosion Keywords: Soil erosion, Erosion, RMMF model, Dakrong Footer Page 12 of 133  ... rẫy khu vực đất dốc, đất có hệ số xói mòn cao Hình Lượng đất xói mòn Hình Phân cấp xói mòn đất Kết luận 3.1 Mô hình RMMF công cụ hữu hiệu để tính toán xói mòn khu vực đồi núi Mô hình thiết kế... 179 2.3 Kết đánh giá khả xói mòn đất mô hình RMMF 2.3.1 Xác định thông số cho mô hình a Năng lượng công kích mưa (KE) Năng lượng mưa lực hạt mưa tác động lên bề mặt đất Trong mô hình RMMF, KE xác... phần (với kích thước pixel 15m) 2.2 Ứng dụng mô hình RMMF đánh giá xói mòn đất Mô hình RMMF dùng để tính lượng đất xói mòn cho khu vực đồi núi Mô hình RMMF đòi hỏi 15 thông số đầu vào thể bảng Footer