đánh giá nguy cơ xói mòn đất do tác động của biến đổi khí hậu khu vực các huyện bắc quang, quang bình, vị xuyên – tỉnh hà giang

BỘ GIÁO DỤC VÀ ĐÀO TẠO

BỘ NÔNG NGHIỆP VÀ PTNT

HỌC VIỆN NÔNG NGHIỆP VIỆT NAM

NGUYỄN VĂN ĐỨC

ĐÁNH GIÁ NGUY CƠ XÓI MÒN ĐẤT DO TÁC ĐỘNG CỦA
BIẾN ĐỔI KHÍ HẬU KHU VỰC CÁC HUYỆN BẮC QUANG,
QUANG BÌNH, VỊ XUYÊN – TỈNH HÀ GIANG

LUẬN VĂN THẠC SĨ

HÀ NỘI, NĂM 2015


BỘ GIÁO DỤC VÀ ĐÀO TẠO

BỘ NÔNG NGHIỆP VÀ PTNT

HỌC VIỆN NÔNG NGHIỆP VIỆT NAM

NGUYỄN VĂN ĐỨC

ĐÁNH GIÁ NGUY CƠ XÓI MÒN ĐẤT DO TÁC ĐỘNG CỦA
BIẾN ĐỔI KHÍ HẬU KHU VỰC CÁC HUYỆN BẮC QUANG,
QUANG BÌNH, VỊ XUYÊN – TỈNH HÀ GIANG

CHUYÊN NGÀNH: KHOA HỌC MÔI TRƯỜNG
MÃ SỐ: 60.44.03.01

NGƯỜI HƯỚNG DẪN KHOA HỌC
TS. NGUYỄN DUY BÌNH

HÀ NỘI, NĂM 2015


LỜI CAM ĐOAN
Tôi xin cam đoan rằng số liệu và kết quả nghiên cứu trong luận văn là
trung thực và chưa được sử dụng để bảo vệ một học vị nào.
Tôi xin cam đoan rằng các thông tin trích dẫn trong luận văn này đã được
ghi rõ nguồn gốc.
Hà Nội, ngày 24 tháng 10 năm 2015
Tác giả luận văn

NGUYỄN VĂN ĐỨC

ii


LỜI CẢM ƠN
Trong quá trình nghiên cứu thực hiện đề tài, tôi đã nhận được sự quan tâm
giúp đỡ của nhiều tập thể, cá nhân trong và ngoài học viện.
Trước hết, tôi xin bày tỏ lòng biết ơn sâu sắc đến các thầy cô giáo của
Khoa Môi Trường, Khoa Quản lý Đất đai - Học Viện Nông nghiệp Việt Nam đã
truyền đạt cho tôi những kiến thức bổ ích và đã tạo mọi điều kiện giúp đỡ tôi
hoàn thành luận văn này.
Tôi xin bày tỏ lòng biết ơn chân thành đến thầy giáo TS. Nguyễn Duy
Bình đã dành nhiều thời gian và tâm huyết, tận tình hướng dẫn, chỉ bảo tôi trong
quá trình thực hiện luận văn này.

Tôi xin chân thành cảm ơn Sở Tài nguyên và Môi trường tỉnh Hà Giang,
Trung tâm Tương tác Biển và Khí quyển, Trung tâm dự báo khí tượng thủy văn
tỉnh Hà Giang, Viện Quy hoạch và Thiết kế Nông nghiệp đã tạo điều kiện, cung
cấp những số liệu cần thiết và giúp đỡ tôi trong quá trình nghiên cứu tại địa bàn.
Cuối cùng, tôi xin bày tỏ lòng biết ơn sâu sắc tới gia đình, người thân, bạn
bè và đồng nghiệp đã luôn động viên khích lệ và giúp đỡ tôi hoàn thành tốt
nghiên cứu của mình.
Hà Nội, ngày 24 tháng 10 năm 2015
Tác giả luận văn

NGUYỄN VĂN ĐỨC

iii


MỤC LỤC
Lời cam đoan

ii

Lời cảm ơn

iii

Mục lục

iv

Danh mục bảng

vii

Danh mục hình

ix

Danh mục chữ viết tắt

xi

MỞ ĐẦU

1

1

Tính cấp thiết của đề tài

1

2

Mục tiêu nghiên cứu

2

3

Yêu cầu

2

Chương 1 TỔNG QUAN TÀI LIỆU

3

1.1

Xói mòn đất

3

1.1.1

Khái niệm và các nguyên nhân xói mòn đất

3

1.1.2

Các phương pháp ước lượng xói mòn đất

6

1.1.3

Công thức USLE và các yếu tố ảnh hưởng đến xói mòn đất

13

1.1.4

Tình hình nghiên cứu và hiện trạng xói mòn đất

21

1.1.5

Các nghiên cứu đánh giá tác động biến đổi khí hậu (BĐKH) đến
xói mòn đất

29

1.2

Biến đổi khí hậu

31

1.2.1

Khái niệm về biến đổi khí hậu

31

1.2.2

Các biểu hiện của biến đổi khí hậu

33

1.2.3

Kịch bản biến đổi khí hậu

37

Chương 2 NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU

46

2.1

Đối tượng nghiên cứu

46

2.2

Phạm vi nghiên cứu

46

2.3

Nội dung nghiên cứu

46

2.3.1

Đánh giá điều kiện tự nhiên - kinh tế xã hội khu vực nghiên cứu

46

2.3.2

Đánh giá các yếu tố xói mòn đất trong khu vực nghiên cứu

46

iv


2.3.3

Xây dựng dữ liệu bản đồ xói mòn đất do mưa năm 2010 và dự báo
năm 2020 với sự thay đổi lượng mưa theo kịch bản biến đổi khí hậu

46

2.3.4

Đánh giá nguy cơ xói mòn đất do biến đổi khí hậu trong khu vực

47

2.3.5

Đề xuất một số giải pháp giảm thiểu xói mòn

47

2.4

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

47

2.4.1

Phương pháp thu thập tài liệu thứ cấp

47

2.4.2

Phương pháp phân tích, tổng hợp dữ liệu

48

2.4.3

Phương pháp công thức tổng hợp xói mòn đất kết hợp GIS

48

2.4.4

Phương pháp thống kê

54

2.4.5

Phương pháp so sánh

54

Chương 3 KẾT QUẢ NGHIÊN CỨU

55

3.1

Điều kiện tự nhiên, kinh tế xã hội các huyện Bắc Quang, Quang
Bình và Vị Xuyên

55

3.1.1

Điều kiện tự nhiên

55

3.1.2

Tình hình phát triển kinh tế, xã hội

64

3.2

Đánh giá các yếu tố xói mòn đất trong khu vực nghiên cứu

68

3.2.1

Đánh giá yếu tố thổ nhưỡng

68

3.2.2

Đánh giá yếu tố lượng mưa và sự thay đổi do biến đổi khí hậu

71

3.2.3

Yếu tố địa hình

74

3.2.4

Đánh giá yếu tố che phủ

75

3.2.5

Đánh giá yếu tố sử dụng đất

77

3.3

Xây dựng dữ liệu bản đồ xói mòn đất do biến đổi lượng mưa theo
kịch bản biến đổi khí hậu

79

3.3.1

Xây dựng dữ liệu bản đồ hệ số xói mòn của đất (K)

79

3.3.2

Xây dựng dữ liệu hệ số xói mòn do mưa (R) năm 2010 và năm
2020 với sự thay đổi lượng mưa theo kịch bản biến đổi khí hậu

83

3.3.3

Xây dựng dữ liệu bản đồ hệ số LS

89

3.3.4

Xây dựng bản đồ hệ số về che phủ (C)

91

3.3.5

Xây dựng dữ liệu bản đồ hệ số P

94

v


3.3.6

Xây dựng dữ liệu bản đồ xói mòn đất năm 2010 và dự báo năm
2020 với sự thay đổi lượng mưa theo kịch bản biến đổi khí hậu

96

3.4

Đánh giá nguy cơ xói mòn đất trong khu vực nghiên cứu

100

3.4.1

Đánh giá mức độ xói mòn đất trong khu vực theo kịch bản BĐKH

100

3.4.2

Đánh giá xu hướng xói mòn đất theo kịch bản biến đổi khí hậu

104

3.5

Đề xuất biện pháp hạn chế xói mòn đất nhằm ứng phó với biến đổi
khí hậu

108

3.5.1

Các biện pháp tăng cường khả năng che phủ đất

108

3.5.2

Các biện pháp công trình

109

3.5.3

Các biện pháp canh tác

113

KẾT LUẬN VÀ KIẾN NGHỊ

114

1

Kết luận

114

2

Kiến nghị

115

TÀI LIỆU THAM KHẢO

116

vi

DANH MỤC BẢNG
STT

Tên bảng

Trang

1.1

Một số công thức tính hệ số R

16

1.2

Giá trị trị hệ số R tại một số khu vực phía Bắc Việt Nam

16

1.3

Hệ số xói mòn của các loại đất trên địa bàn tỉnh Hà Giang

18

1.4

Giá trị hệ số C một của một số loại lớp phủ trên địa bàn tỉnh Hà

Giang

1.5

21

Tỷ lệ đất dốc bị thoái hóa do ảnh hưởng của xói mòn phân theo các
vùng sinh thái ở Việt Nam

23

1.6

Lượng đất mất bình quân do xói mòn ở một số vùng đặc trưng

23

1.7

Tổng hợp mức độ xói mòn đất và diện tích khu vực nghiên cứu
năm 2008

1.8

28

Mức thay đổi lượng mưa năm (%) qua các thập kỷ của thế kỷ 21 so
với thời kỳ 1980 - 1999 ứng với kịch bản phát thải trung bình (B2)

44

3.1

Nhiệt độ trung bình năm của trạm Hà Giang và trạm Bắc Quang

58

3.2

Lượng mưa trung bình năm tại trạm Hà Giang và trạm Bắc Quang

59

3.3

Độ ẩm trung bình năm tại trạm Hà Giang và trạm Bắc Quang

59

3.4

Giờ nắng trong năm tại trạm Hà Giang và trạm Bắc Quang

60

3.5

Vận tốc gió trung bình tại trạm Hà Giang và trạm Bắc Quang

60

3.6

Thống kê diện tích một số loại đất trong khu vực nghiên cứu

61

3.7

Số lượng đàn gia súc trong khu vực nghiên cứu

66

3.8

Diện tích các loại đất khu vực nghiên cứu theo nhóm và loại đất

69

3.9

Mức tăng lượng mưa năm(%) khu vực nghiên cứu qua các thập kỷ
của thế kỷ 21 so với thời kỳ 1980 - 2000

3.10

72

Mức thay đổi lượng mưa mùa mưa và mùa khô qua các thập kỷ của
thế kỷ 21 so với thời kỳ 1980 - 1999 ứng với kịch bản phát thải

trung bình (B2) tại các trạm khí tượng trong khu vực

73

3.11

Thống kê các loại đối tượng che phủ trong khu vực

77

3.12

Diện tích các loại đất được áp dụng biện pháp canh tác trong khu vực

78

vii


3.13

Hệ số xói mòn của các loại đất trong khu vực nghiên cứu

80

3.14

Hệ số xói mòn của thảm thực vật trong khu vực nghiên cứu

91

3.15

Hệ số xói mòn của các loại hình canh tác và sử dụng đất trong khu vực

95

3.16

Mức độ xói mòn đất do mưa đến năm 2020 theo kịch bản BĐKH
khu vực nghiên cứu

3.17

100

Phân chia tỷ lệ (%) diện tích (ha) theo các cấp độ xói mòn và loại
hình canh tác trong khu vực nghiên cứu

3.18

Thống kê diện tích đất theo cấp độ xói mòn đến năm 2020 khu vực
huyện Bắc Quang.

3.19

103

Thống kê diện tích đất theo cấp độ xói mòn đến năm 2020 khu vực
huyện Quang Bình

3.20

102

104

Thống kê diện tích đất theo cấp độ xói mòn đến năm 2020 khu vực
huyện Vị Xuyên

104

3.21

So sánh thống kê số liệu xói mòn khu vực nghiên cứu tại các thời điểm 105

3.22

Thống kê số liệu xói mòn khu vực huyện Bắc Quang tại các thời điểm

3.23

Thống kê số liệu xói mòn khu vực huyện Quang Bình tại các thời điểm 107

3.24

Thống kê số liệu xói mòn khu vực huyện Vị Xuyên tại các thời điểm

viii

106
107


DANH MỤC HÌNH
STT

Tên hình

Trang

1.1

Biến đổi của nhiệt độ bề mặt Trái Đất theo thời gian

34

1.2

Xu hướng biến đổi một số khí nhà kính đến 2005

34

1.3

Kịch bản BĐKH đối với nhiệt độ cho 6 khu vực của châu Á

40

1.4

Kịch bản BĐKH đối với nhiệt độ cho khu vực Châu Á (IPCC/2007)

41

2.1

Bảng tra toán đồ hệ số K của Wischmeier và Smith (1978)

51

2.2

Quá trình tính giá trị hệ số LS

52

2.3

Quy trình mô phỏng mức độ xói mòn đất

54

3.1

Vị trí khu vực nghiên cứu

56

3.2

Mô hình số độ cao khu vực nghiên cứu

57

3.3

Dữ liệu bản đồ thổ nhưỡng khu vực nghiên cứu

70

3.4

Diễn biến tổng lượng mưa năm tại các trạm khí tượng khu vực

72

3.5

Mô hình số dữ liệu độ cao khu vực

74

3.6

Mô hình số dữ liệu độ dốc khu vực

75

3.7

Dữ liệu hiện trạng lớp phủ năm 2010 khu vực nghiên cứu

76

3.8

Dữ liệu hiện trạng các biện pháp canh tác năm 2010 trong khu vực
79

nghiên cứu
3.9

Lựa chọn các đối tượng theo loại đất (a) và gán giá trị hệ số K cho
loại đất được chọn (b)

3.10

81

Bảng thuộc tính (a) và dữ liệu bản đồ hệ số K (b) của khu vực
nghiên cứu

81

3.11

Chuyển đổi dữ liệu hệ số K sang dạng Raster

82

3.12

Dữ liệu bản đồ hệ số xói mòn của đất (K) trong khu vực nghiên cứu

82

3.13

Dữ liệu bản đồ phân bố lượng mưa năm 2010 khu vực nghiên cứu

83

3.14

Xử lý dữ liệu mức tăng lượng mưa bằng phần mềm ArcGIS

84

3.15

Kết quả xử lý dữ liệu phân bố lượng mưa trung bình (a) và dữ liệu
phân bố mức tăng lượng mưa (b) trong khu vực nghiên cứu

ix

85


3.16

Xử lý tính toán số liệu phân bố lượng mưa năm 2020 trên phần
mềm ArcGIS

86

3.17

Dữ liệu bản đồ phân bố lượng mưa năm 2020 khu vực nghiên cứu

86

3.18

Tính toán hệ số R từ dữ liệu phân bố lượng mưa trên phần mềm
ArcGIS

88

3.19

Dữ liệu hệ số xói mòn do mưa (R) năm 2010

88

3.20

Dữ liệu hệ số xói mòn do mưa (R) năm 2020 theo kịch bản BĐKH

89

3.21

Tách riêng dữ liệu hệ số LS khu vực nghiên cứu

90

3.22

Dữ liệu hệ số LS khu vực nghiên cứu

90

3.23

Lựa chọn các đối tượng theo loại che phủ rừng (a) và gán giá trị hệ
số C cho loại che phủ rừng được chọn (b)

3.24

Bảng thuộc tính (a) và dữ liệu bản đồ hệ số C (b) của khu vực
nghiên cứu

3.25

92
92

Dữ liệu bản đồ hệ số về che phủ đất (C) trong khu vực nghiên cứu
năm 2010

93

3.26

Tách riêng dữ liệu khu vực nghiên cứu từ dữ liệu thu thập

94

3.27

Dữ liệu bản đồ hệ số P khu vực nghiên cứu năm 2010

96

3.28

Tính toán mức độ xói mòn đất bằng phần mềm ArcGIS

97

3.29

Dữ liệu bản đồ xói mòn đất năm 2010 khu vực nghiên cứu

98

3.30

Dữ liệu bản đồ nguy cơ xói mòn đất năm 2020 theo kịch bản biến

đổi khí hậu khu vực nghiên cứu

3.31

99

Xây dựng rào chắn

110

(Nguyễn Viết Khoa và cs, 2008)

110

3.32

Bẫy đất

110

3.33

Đào rãnh theo đường đồng mức

111

3.34

Kích thước của rãnh

112

3.35

Tạo nền móng để xây dựng tường đá

112

x


DANH MỤC CHỮ VIẾT TẮT
BĐKH Biến đổi khí hậu
GIS

IPCC

Hệ thống thông tin địa lý (Geographic
information system)
Ban liên chính phủ về Biến đổi khí hậu (The
Intergovernmental Panel on Climate Change)

TCVN Tiêu chuẩn Việt Nam
TNMT Tài nguyên môi trường

xi


MỞ ĐẦU

1. Tính cấp thiết của đề tài
Hà Giang là tỉnh vùng núi cao ở địa đầu cực Bắc của Việt Nam. Đặc điểm
địa hình khu vực bị chia cắt với nhiều núi cao, vực sâu, nhiều dòng sông có độ
dốc lớn đã gây ra rất nhiều khó khăn cho hoạt động phát triển kinh tế - xã hội của
vùng đặc biệt là hoạt động nông nghiệp.
Trên địa bàn tỉnh Hà Giang hiện nay, phần lớn dân số làm nông nghiệp.
Tuy nhiên, do những khó khăn về địa hình, khí hậu nên năng suất nông nghiệp
của tỉnh thuộc loại thấp so với năng suất trung bình trong cả nước. Cả tỉnh chỉ có
khu vực các huyện Bắc Quang, Quang Bình, Vị Xuyên có năng suất cây trồng
cao do có điều kiện ít khó khăn hơn. Do đó, sự ưu tiên tập trung tăng năng suất
cây trồng của khu vực này nhằm bổ sung cho các khu vực khác là một trong
những định hướng trong chiến lược đảm bảo an toàn lương thực cho người dân
trong tỉnh.
Bên cạnh đó, theo kịch bản biến đổi khí hậu của Việt Nam cho thấy: Tại
khu vực phía Bắc khí hậu sẽ có những thay đổi lớn, điển hình là sự tăng lên về
nhiệt độ và lượng mưa trung bình hàng năm, mùa mưa trong năm rút ngắn lại,
xuất hiện tình trạng khô hạn kéo dài... Các trận mưa tập trung với cường độ lớn
làm tăng khả năng xói mòn của đất, chất dinh dưỡng cũng như các tính chất lý
hóa của đất bị thay đổi theo chiều hướng tiêu cực. Đây chính là một trong những
nguyên nhân dẫn đến suy thoái đất nông nghiệp. Hiện tượng này càng xảy ra
mạnh mẽ đối với các khu vực có lượng mưa lớn như khu vực ba huyện trên.
Trước các vấn đề cấp bách đó, đề tài: “Đánh giá nguy cơ xói mòn đất do
tác động của biến đổi khí hậu khu vực các huyện Bắc Quang, Quang Bình, Vị
Xuyên – tỉnh Hà Giang” được xây dựng nhằm dự báo, đánh giá về nguy cơ xói
mòn đất tại các vị trí trong khu vực để từ đó đưa ra các định hướng lâu dài về sử
dụng đất bền vững trong tương lai.

1

2. Mục tiêu nghiên cứu
- Đánh giá nguy cơ xói mòn đất do sự thay đổi lượng mưa theo kịch bản
biến đổi khí hậu đến năm 2020 khu vực ba huyện Bắc Quang, Quang Bình, Vị
Xuyên – tỉnh Hà Giang.
- Đề xuất các giải pháp trong sử dụng đất nhằm hạn chế các tác động gây
nên xói mòn đất trong khu vực nghiên cứu.
3. Yêu cầu
- Thu thập các thông tin về điều kiện tự nhiên, kinh tế - xã hội khu vực
nghiên cứu.
- Thu thập các dữ liệu liên quan làm đầu vào cho quá trình tính toán và mô
phỏng nguy cơ xói mòn đất.
- Thiết kế, xây dựng các lớp thông tin làm dữ liệu đầu vào cho quá trình
tính toán và mô phỏng mức độ xón mòn đất.
- Sử dụng phương pháp chồng ghép dữ liệu của công cụ GIS để tính toán
và mô phỏng mức độ xói mòn đất.

2


Chương 1. TỔNG QUAN TÀI LIỆU
1.1. Xói mòn đất
1.1.1. Khái niệm và các nguyên nhân xói mòn đất
1.1.1.1. Khái niệm về xói mòn đất
Sự tồn tại của cuộc sống con người phụ thuộc rất nhiều vào lớp đất trồng
trọt để sản xuất ra lương thực, thực phẩm và các nguyên liệu sản xuất công
nghiệp phục vụ cho đời sống của con nguời. Tuy nhiên sự tồn tại của lớp đất có
khả năng canh tác này lại luôn chịu những tác động mạnh mẽ của các điều kiện
tự nhiên và các hoạt động canh tác do con người có thể làm cho chúng bị thoái
hóa và dần mất đi khả năng sản xuất, trong đó một trong những nguyên nhân làm
cho đất bị thoái hóa mạnh nhất là do xói mòn.

Xói mòn đất là quá trình làm mất lớp đất trên mặt và phá huỷ các tầng đất
bên dưới do tác động của nước mưa, băng tuyết tan hoặc do gió. Ðối với đất sản
xuất nông nghiệp thì nước và gió là hai tác nhân quan trọng nhất gây ra xói mòn
và các tác nhân này có mức độ ảnh hưởng tăng giảm khác nhau theo các hoạt
động của con người đối với đất đai.
Hiện tượng mất đất do xói mòn mạnh hơn rất nhiều so với sự tạo thành đất
trong quá trình tự nhiên, một vài cm đất có thể bị mất đi chỉ trong một vài trận
mưa giông hoặc gió lốc trong khi để có được vài cm đất đó cần phải có thời gian
hàng trăm năm, thậm chí hàng ngàn năm mới tạo ra được. Trên thế giới hầu như
không có quốc gia nào là không chịu ảnh hưởng của xói mòn, nhất là ảnh hưởng
của xói mòn do nước và do gió (Trần Văn Chính và cs., 2006).
Có nhiều quan niệm về xói mòn khác nhau. Theo các tác giả nước ngoài,
xói mòn đất là những hiện tượng phá hủy và cuốn trôi theo đất cũng như quặng
xốp bằng dòng nước và gió thể hiện dưới nhiều hình thức và rất phổ biến. R.P.C
Morgan (2005) cũng cho rằng xói mòn đất là một quá trình gồm hai pha bao gồm
sự tách rời của các phần tử nhỏ từ mặt đất sau đó vận chuyển chúng dưới các tác
nhân gây xói như nước chảy và gió. Khi năng lượng không còn đủ để vận chuyển
các phần tử này, pha thứ ba – quá trình bồi lắng – sẽ xảy ra.
3


Theo một trong những cách tiếp cận khác khi nghiên cứu về lớp phủ thực
vật thì xói mòn là một quá trình động lực phá hủy độ màu mỡ của đất, làm mất
trạng thái cân bằng của cả vùng bị xói mòn lẫn vùng bị bồi tụ (Nguyễn Quang
Mỹ, 2005).
Như vậy, xói mòn đất được xem xét trên quan điểm là một quá trình động
lực, bao gồm sự phá hủy các lớp đất đá, mùn và vận chuyển chúng đi xa dưới tác
động của các nhân tố gây xói như gió, nước, băng, tuyết tan hoặc hoạt động của
sinh vật bao gồm cả các yếu tố nhân sinh.
Có thể nhận thấy đối với thực tiễn sản xuất nông nghiệp không có sự thoái

hóa đất nào mạnh và hiểm họa hơn xói mòn đất bởi nó liên quan đồng thời tới
các quá trình mất đất, mất chất dinh dưỡng và nước cho cây trồng đồng thời còn
gây ra các tác động xấu đến môi trường. Do đó việc nghiên cứu xói mòn là vô
cùng cần thiết cho mọi quốc gia, đặc biệt đối với nước ta là một nước nằm trong
vành đai nhiệt đới với 3/4 diện tích đất tự nhiên là đồi núi, thường xuyên phải
hứng chịu các hậu quả do xói mòn gây ra thì việc khống chế hiện tượng xói mòn
đất càng trở nên cực kỳ quan trọng để bảo vệ độ phì nhiêu của đất và bảo vệ môi
trường sinh thái (Trần Văn Chính và cs., 2006).
1.1.1.2. Các nguyên nhân dẫn đến xói mòn đất
Theo Giáo trình thổ nhưỡng Đại học Nông nghiệp Hà Nội (Trần Văn
Chính và cs., 2006), các nguyên nhân dẫn đến xói mòn đất bao gồm:
a, Xói mòn do nước
Kiểu xói mòn do nước gây ra do tác động của nước chảy tràn trên bề mặt
(nước mưa, băng tuyết tan hay tưới tràn). Hiện tượng xói mòn do nước gây ra
đối với đất sản xuất nông nghiệp thường mạnh nhất ở các bề mặt đất trống, sau
khi làm đất chuẩn bị gieo trồng. Ðể xảy ra xói mòn, nước cần có năng lượng để
tách các hạt đất, rồi sau đó vận chuyển chúng đi. Mưa và nước có thể tách được
các hạt đất song việc vận chuyển được chúng đi bao xa phải phụ thuộc vào
dòng chảy. Tác động của mưa gây ra xói mòn đối với đất gồm các tác động va
đập, phá vỡ, làm tách rời các hạt đất và sau đó vận chuyển các hạt đất bị phá
hủy theo các dòng chảy tràn trên mặt đất. Dòng chảy của nước có thể tạo ra các
4


rãnh xói, khe xói hoặc bị bóc theo từng lớp bao gồm kiểu xói mòn thẳng và
kiểu xói mòn phẳng.
Khi lớp đất trên bề mặt bị xói mòn thì rất khó khôi phục và những thiệt
hại của xói mòn có thể ảnh hưởng nghiêm trọng tới sức sản xuất của đất. Ở
vùng nhiệt đới có những nơi xói mòn làm mất 3 cm đất mặt hàng năm. Riêng
vùng đồi núi hàng năm bình quân mất đi khoảng 2 cm điều này làm cho đất ở

đây bị thoái hóa nhanh chóng. Trên những vùng đất cao, dốc, mưa lớn còn tạo
nên những dòng chảy cực đại trên sườn dốc và ngoài việc bào mòn lớp đất mặt
chúng còn có khả năng tạo ra những dòng xói hoặc rãnh xói.
Hiện tượng xói mòn do nước xảy ra ở khắp mọi nơi trên bề mặt trái đất,
song tập trung mạnh nhất ở những vùng khí hậu nhiệt đới ẩm nơi thường có
tổng lượng mưa hàng năm lớn, tập trung theo mùa với cường độ cao kết hợp
với đất có địa hình cao và dốc đã tạo ra những dòng chảy tràn lớn trên bề mặt
đất. Trong nhiều vùng đất dốc hiện tượng nước chảy tràn trên mặt không chỉ
làm mất đi một lượng nước lớn trong mùa mưa (khoảng 50 - 60% lượng mưa
hàng năm) mà kèm theo đó là hiện tượng đất bị mất do xói mòn mạnh và đây
mới chính là thiệt hại nghiêm trọng đối với đất canh tác hơn cả vì chúng làm
cho lớp đất mặt bị bào mòn dẫn đến hậu quả là đất bị mất dần và thậm chí
không còn khả năng sản xuất. Khi xói mòn xả y ra, thông thường những thành
phần hạt đất nhỏ, mịn trên cùng lớp đất mặt bị đẩy đi trước tiên và ở lớp đất này
thường tập trung độ phì nhiêu cao nhất do vậy hàm lượng các chất dinh dưỡng
bị mất đi trong đất do xói mòn cũng rất lớn.
Các dạng xói mòn do nước bao gồm xói mòn theo lớp và xói mòn theo
các khe, rãnh tạo ra các mương xói.
b, Xói mòn do gió
Kiểu xói mòn do gió là hiện tượng xói mòn gây ra bởi sức gió. Ðây là
hiện tượng xói mòn có thể xảy ra ở bất kì nơi nào khi có những điều kiện thuận
lợi sau đây:
- Ðất khô, tơi và bị tách nhỏ đến mức độ gió có thể cuốn đi;
- Mặt đất phẳng, có ít thực vật che phủ thuận lợi cho việc di chuyển của gió;
5


- Diện tích đất đủ rộng và tốc độ gió đủ mạnh để mang được các hạt đất đi.
Thông thường đất cát là loại đất rất dễ bị xói mòn do gió vì sự liên kết
giữa các hạt cát là rất nhỏ, đất lại bị khô nhanh. Ảnh hưởng của xói mòn do gió

gây ra khá phức tạp, đất bị chuyển dịch đi có thể dưới các hình thức nhảy cóc,
trườn trên bề mặt hoặc bay lơ lửng.
1.1.2. Các phương pháp ước lượng xói mòn đất
1.1.2.1. Phương pháp thực nghiệm
a, Định lượng xói mòn bằng thực nghiệm
Phương pháp đóng cọc
Trên sườn dốc bố trí cọc ở 3 vị trí đỉnh, sườn dốc và chân dốc. Điểm đặt
cọc ở chân dốc phải ở vị trí không bị xói mòn và không bồi tụ (điểm chuyển
tiếp). Sau từng thời gian nhất định, do độ sâu bị bào mòn (so với thời điểm đóng
cọc) được kết quả cọc ở đỉnh có độ sâu bào mòn là hd, cọc ở sườn là hs.Ta có độ
sâu bào mòn trung bình được tính theo công thức:
0,5 (hs+hd)
Hạn chế của phương pháp này là dễ mất cọc, thời gian theo dõi lâu, cọc
đóng ảnh hưởng đến dòng chảy, phương pháp xác định độ sâu bào mòn bằng mắt
thông thường v.v... Do vậy độ chính xác của lượng đất mất xác định được là
không cao. Tuy nhiên trong thực tế vẫn có thể sử dụng được để đánh giá chung
về hiện tượng xói mòn đất.
Phương pháp phẫu diện
Phương pháp này do Zakharov đề xướng trên cơ sở các khu vực có hình
thái địa lý khác nhau, nếu không bị bào mòn thì sẽ có độ dầy các tầng đất tương
đối ổn định. Nếu bị xói mòn thì độ dầy tầng A sẽ bị giảm, quá trình xói mòn kéo
dài thì dẫn sang tầng B và tầng C, v.v. cũng bị xói mòn.
Việc xác định mức độ xói mòn được tiến hành bằng cách so sánh các vị trí
đất chưa bị xói mòn với các vị trí bị xói mòn của cùng một loại đất trong cùng
các điều kiện đồng nhất.
Phương pháp này nói chung độ chính xác không cao, cần phải có kinh
nghiệm nhất định mới xác định được vị trí không bị xói mòn, không bị bồi lắng,
6

đồng thời phải xác định được thời điểm bắt đầu trồng trọt. Điều này trong thực tế
là hết sức khó khăn. Để khắc phục nhược điểm này, khi áp dụng đối với khu vực
làm đất, người ta tính từ thời điểm làm đất sau cùng và xác định độ bào mòn của
tầng đất cày.
Tuy nhiên phương pháp Zakharov có thể đánh giá tổng quát về hiện tượng
xói mòn trên khu vực nghiên cứu, từ đó xác định được phương pháp nghiên cứu
đánh giá xói mòn đất cụ thể, chính xác hơn (Hồ Kiệt, 2002).
Phương pháp đo bộ rễ thực vật
Phương pháp này được nhiều nhà khoa học Châu Âu, Châu Phi áp dụng ở
những khu vực trồng cỏ lâu năm. Cơ sở của phương pháp này là khi bị bào mòn
cổ rễ thực vật vẫn ở độ cao nhất định và như vậy cổ rễ nổi cao lên trên mặt đất.
Xác định thời điểm trồng đến thời điểm nghiên cứu ta xác định được mức độ xói
mòn trong khoảng thời gian đó và tính được lượng đất bị bào mòn hàng năm.
Thực chất độ chính xác của phương pháp này không cao song cũng như phương
pháp đóng cọc, cho phép đánh giá xói mòn đất ở khu vực, mặt khác để áp dụng,
không đòi hỏi phương tiện nghiên cứu phức tạp.
Phương pháp thu hứng, thu lượng đất bị xói mòn
Phương pháp này thực hiện như sau: trên diện tích nhất định được cô lập
khỏi sự ảnh hưởng dòng chảy và lượng đất xói mòn từ các diện tích khác (dùng
bờ, rãnh ngăn...), phía cuối ô đất, đặt bộ phận thu đất. Sau từng khoảng thời gian,
cân lượng đất bồi lắng thu được và tính toán sẽ có lượng đất bị xói mòn trên một
đơn vị diện tích. Có hai phương pháp như sau:
- Phương pháp thu toàn bộ: suốt cả bề rộng của khoảng nghiên cứu, đào
một đoạn rãnh, lót ni lông, phía dưới có đóng các cọc rỗng để thoát nước. Đem
lượng đất thu được cân toàn bộ.
- Phương pháp thu một bộ phận đất bị xói mòn: Lượng nước (bao gồm cả
đất bị xói mòn) được trộn đều và dẫn đến dụng cụ phía trước (Runoff Plots:Vnotch hoặc H-flume) thành các phần bằng nhau, sao cho đại diện cho lượng nước
toàn bộ và đơn giản cho quá trình tính toán rồi đưa vào bộ phận thu hứng
nước.Lượng nước thu được để lắng, lấy phần lắng đọng và xác định khối lượng.
7

Từ đó tính toán được lượng đất bị xói mòn trên một đơn vị diện tích trong thời
gian nghiên cứu ở khu vực.
Phương pháp xác định xói mòn bằng cách thu đất bị rửa trôi cho kết quả
chính xác, đánh giá được lượng đất mất do tất cả các dạng xói mòn nước. Tuy
nhiên nó đòi hỏi thi công, thí nghiệm phức tạp và trang bị đắt tiền, do vậy
phương pháp này thường được sử dụng để phục vụ cho việc nghiên cứu, thiết lập
các tham số trong các phương trình hoặc các mô hình chuẩn đoán xói mòn đất
(Nguyễn Duy Bình, 2012).
b, Chất đánh dấu
Để hiểu rõ về cách xác định xói mòn đất cần những thông tin về sự phân
bố không gian và mức độ xói mòn và lắng đọng, các nhà khoa học đã đưa ra các
phương pháp như phương pháp đo từ khoáng vật và phương pháp sử dụng
nguyên tử đánh dấu:
Phương pháp đo từ khoáng vật
Phương pháp này sử dụng nguyên tử nuclit phóng xạ để ước tính sự tích
lũy hay rửa trôi của đất khi các nuclit phóng xạ hấp thụ. Dearing et al. (1986)
cũng sử dụng phương pháp này để đo từ tính vết tích sự di chuyển tầng mặt (tầng
đế cày) trên sườn dốc ở Oxfordshire (Mỹ). Cơ sở cho việc xác định vết tích này
đó là sự khác biệt từ tính của tầng mặt và tầng phụ thông qua quá trình hình
thành thổ nhưỡng được phân tích để quyết định mức độ xói mòn và lắng đọng
trên một lưu vực (Thompson và Oldfeld, 1986). Các chỉ thị từ tính của xói mòn
đất kết hợp với quá trình hình thành đất được công bố và dữ liệu các chất lơ lửng
lặng đọng từ khu vực thí nghiêm trên 1 lưu vực để xác định quá trình di chuyển
trầm tích, lắng đọng. Tuy nhiên sử dụng rộng rãi phương pháp này cũng gặp
nhiều khó khăn ở mức độ không gian do số lượng mẫu lớn, giá thành phân tích
cao và đòi hỏi nhiều lao động. Hơn nữa phạm vi hạn chế của những nghiên cứu
tiên phong cần thiết phải thúc đẩy tiếp tục việc nghiên cứu về công nghệ từ tính
khoáng vật để hiểu rõ hơn về phân bố lại đất và dòng chảy bùn cát trên sườn dốc

(Dearing et al., 1986).

8


Phương pháp 137Cs
Phương pháp này đã được sử dụng ở cả các nước phát triển cũng như các
nước đang phát triển (Ritchie, 2001).137Cs là chất phóng xạ nhân tạo với chu kỳ
phân rã là 30,12 năm được tạo ra bởi sự phân rã hạt nhân. Giá trị 137Cs khi nó là
chất phóng xạ đánh dấu được hấp thụ nhanh chóng và mạnh mẽ bởi các hạt đất
min như là khoáng vật sét và hợp chất mùn (Tamura, 1964 và Bachhuber et al.,
1982). Trong hệ sinh thái nông nghiệp, sự phân phối này được thừa nhận là kết
quả trực tiếp của hiện tượng xói mòn đất, vận chuyển và phân hủy các hạt đất
xuất hiện trong suốt giai đoạn chính của sự lắng đọng ở trong bầu khí quyển tới
thời điểm lấy mẫu (Walling và Quine, 1995). Để định lượng ước tính tốc độ xói
mòn và lắng đọng từ việc đo 137Cs thì cần thiết phải thiết lập mối quan hệ giữa độ
lớn của độ lệch từ bảng kiểm kê tham chiếu và khu vực mất đất hay tăng thêm.
c, Viễn thám
Viễn thám đã và đang được sử dụng cho những nghiên cứu về xói mòn đất
thông qua việc giải đoán ảnh vệ tinh bao gồm giải đoán những hình ảnh xói mòn
cũng như cung cấp dữ liệu cho mô hình đầu vào. Bắt đầu từ năm 1972 với sự
khởi đầu của vệ tinh Landsat1, ảnh vệ tinh đã được phổ biến và sử dụng rộng rãi
cho cộng đồng các nhà làm khoa học. Trong suốt 30 năm, có rất nhiều nghiên
cứu đã được công bố về việc áp dụng ảnh vệ tinh trong việc đánh giá xói mòn đất
trong nhiều cách thức khác nhau, các dữ liệu vệ tinh chủ yếu được lấy từ vệ tinh
Landsat, SPORT, ASTER, ENVISAT ASAR, IKONOS, QickBird.
Hơn nữa, ngoài việc dò tìm phát hiện ra xói mòn thì dữ liệu vệ tinh còn
được ứng dụng cho việc đánh giá ảnh hưởng xói mòn đất. Nhiều nghiên cứu đã
chỉ ra mối tương tác có ý nghĩa giữa sự tập trung các chất lắng đọng trong một
thể tích nước nhất định và phản xạ phổ trên bề mặt khí quyển được lấy từ dữ liệu

viễn thám như Landsat (Nellis et al., 1998), IRS LISS-1 (Choubey, 1998) và đa
phổ SPOT HRV (Chacontorres et al., 1992). Sự kết hợp giữa phương trình hàm
lắng đọng lơ lửng và dữ liệu viễn thám là cơ sở cho việc đánh giá những vùng có
tốc độ xói mòn nhanh.

9


1.1.2.2. Phương pháp mô hình
a, Mô hình công thức thực nghiệm/ thống kê
Công thức USLE
Sau rất nhiều năm nghiên cứu kết hợp với kinh nghiệm thực tiễn, các nhà
khoa học Weischmaier và Smith đã xác định được phương trình dự tính lượng
đất xói mòn do nước gây ra, thường được gọi là phương trình mất đất phổ dụng
có công thức sau:
A=R.K.L.S.C.P
Trong đó:
A - Lượng đất mất bình quân trong năm (tấn/ha/ năm);
R - Yếu tố mưa và dòng chảy;
K - Hệ số bào mòn của đất (tấn/ha/ đơn vị chỉ số xói mòn);
L - Yếu tố chiều dài của sườn dốc;
S - Yếu tố độ dốc;
C - Yếu tố che phủ và quản lý đất;
P - Yếu tố hoạt động điều tiết chống xói mòn;
Việc hiểu rõ các yếu tố này sẽ giúp cho xác định hướng khống chế xói
mòn đối với đất (Trần Văn Chính và cs., 2006). Chi tiết phương pháp xác định
giá trị các yếu tố trên được trình bày riêng trong ở mục 1.1.3.
Mô hình MMF
Mô hình Morgan/Morgan/Finney model (Morgan, 1986) được phát triển
nhằm dự báo lượng đất bị xói mòn hàng năm cho diện tích cỡ thửa ruộng trên

đất dốc. Mô hình này bao gồm phần mô phỏng quá trình dòng chảy và phần mô
phỏng quá trình bồi lắng. Nó được coi là có tính chất động thái nhiều hơn so
với mô hình USLE.
Trong mô hình MMF, trước hết mô hình tính toán động năng (đơn vị J/M)
của dòng chảy với công thức E = R (11,9 + 8,7 log10 I) trong đó lượng mưa và
cường độ gây xói của trận mưa là dữ liệu đầu vào. Sau đó lượng hạt đất bị xói
mòn (F, kg/m2) được ước lượng theo công thức, F = K × E(0,05P) × 10-3 , trong đó
P = phần trăm lượng mưa rơi xuống tán lá và thân cành cây cối để sau đó bị bốc
10


hơi hay tạo dòng chảy theo thân cây, phụ thuộc vào loại che phủ thực vật; và hệ
số K (chỉ số đất bị tách rời) là khối lượng hạt đất bị tách rời trên đơn vị năng
lượng mưa rơi (gram/Joule).
Lưu lượng dòng chảy tràn trên mặt đất (Q) là hàm số của lượng mưa và
sức giữ ẩm đồng ruộng của đất (soil moisture capacity) tương ứng với hiện
trạng lớp phủ thực vật. Dữ liệu đầu vào của mô hình bao gồm lượng mưa trung
bình ngày (Ro) và sức giữ ẩm đồng ruộng của đất ứng với hiện trạng lớp phủ
thực vật thực tế (Rc). Rc được tính theo công thức
Rc =1000* MS × BD × RD × (Et/Eo)0,5
Trong đó MS = độ ẩm đồng ruộng (%); BD = dung trọng đất (bulk
density) của tầng đất mặt (Mg/m3); RD = độ sâu tầng rễ cây hoạt động (m);
Et/Eo = tỷ lệ giữa lượng bốc thoát nước thực tế và tiềm năng.
Sau đó lưu lượng dòng chảy tràn (Q) sẽ bằng Q = R (-Rc/Ro). Khả năng
vận chuyển đất của dòng chảy (T) phụ thuộc vào lưu lượng dòng chảy tràn
(Q), độ dốc địa hình (S) và độ che phủ của cây trồng (C), theo phương trình
T = C Q2 × sin S × 10-3 (kg/m2). Cuối cùng giá trị bé nhất giữa khả năng của
dòng chảy tràn về vận chuyển đất và tách rời đất được coi như lượng đất bị
xói mòn (Morgan, 1995).
b, Mô hình động thái trên cơ sở vật lý

Những mô hình này dựa trên việc áp dụng định luật bảo toàn vật chất và
năng lượng. Phần lớn chúng sử dụng các phương trình vi phân về tính liên tục bảo
toàn vật chất khi vật chất biến động theo không gian và thời gian (Morgan, 1986).
Mô hình EUROSEM
EUROSEM (EURopean Soil Erosion Model, Mô hình xói mòn đất của
Châu Âu) là mô hình dự báo xói mòn đất dựa trên mô phỏng từng quá trình vật
lý trong cơ chế xói mòn đất. Nó được thiết kế với mục đích mô phỏng cho từng
trận mưa riêng biệt và có thể sử dụng để đánh giá tác động của các biện pháp
bảo vệ đất (Morgan, 1995).
Mô hình KINEROS
KINEROS (KInematic EROsion Simulation Model, Mô hình mô phỏng
11


động học xói mòn) diễn tả toán học các quá trình tạo nên dòng chảy tràn trên
mặt đất. Mô hình KINEROS thuộc loại mô hình dựa trên nguyên lý các quá
trình vật lý, mô phỏng quá trình xói mòn đất trên diện tích lưu vực nhỏ khi có
trận mưa xẩy ra. Lượng mưa hữu hiệu được tích tụ lại trên mặt đất cho đến khi
độ sâu nước đủ lớn để tạo nên dòng chảy tràn.
Sau đó dòng chảy tràn được mô phỏng trên cơ sở tìm nghiệm phương
trình liên tục với giả thiết xấp xỉ cơ học. Vùng địa hình được biểu thị bằng kết
nối nhiều phần tử diện tích và kênh nối. Như vậy dòng chảy được tính truyền lần
lượt theo thứ tự các diện tích và kênh dẫn phần tử cho đến cửa ra của lưu vực.
Lượng đất xói mòn được ước tính cho từng diện tích phần tử và được
tính truyền hay cho bồi lắng theo thứ tự liên kết các phần tử. KINEROS mô
phỏng lượng vận chuyển đất bị xói mòn qua quá trình tách rời của giọt mưa rơi,
vận chuyển của dòng chảy tập trung. Mô hình sử dụng phương pháp cân bằng
vật chất để tính toán dòng bùn cát (Doe et al., 1999).
Mô hình WEPP
Mô hình dự án dự báo xói mòn đất do nước WEPP (Water erosion

prediction Project) được coi là một công nghệ mới về dự báo xói mòn đất.
WEPP được phát triển với mục đích thay thế công thức phổ cập xói mòn đất
USLE, MUSLE, và RUSLE, và tăng cường khả năng mô phỏng và dự báo đối
với nhiều địa hình địa mạo khác nhau. Mô hình có thể ứng dụng được cho cả
vùng đất canh tác, đồng cỏ chăn nuôi, đất rừng, đất xây dựng và đất vùng khai
thác khoáng sản. WEPP thuộc loại mô hình dựa trên nguyên lý các quá trình vật
lý, mô phỏng liên tục theo thời gian và không gian quá trình xói mòn đất. Nó
được phát triển dựa trên nguyên lý cơ bản về mô phỏng xác suất quá trình khí
tượng, lý thuyết dòng thấm, thủy văn, vật lý đất, và khoa học cây trồng, thủy
lực và động học xói mòn đất. Những điểm nổi bật của WEPP bao gồm khả năng
dự báo phân bổ lượng đất xói theo thời gian và không gian (tổng lượng đất bị
xói trên toàn bộ sườn dốc hay tại từng vị trí, ở mọi thời điểm theo thời gian như
từng ngày, tháng, hay trung bình năm). Đồng thời mô hình cho phép ngoại suy
kết quả mô hình cho rất nhiều điều kiện khác nhau. Khi ứng dụng cho một lưu
12


vực sông thì mô hình còn ước tính lưu lượng dòng chảy bùn cát sinh ra từ toán
bộ diện tích mô phỏng.
WEPP bao gồm hai phiên bản, một cho diện tích sườn dốc và một cho
toàn bộ lưu vực. Phiên bản diện tích sườn dốc được tích hợp trong phiên bản
lưu vực, nhằm tích toán lượng đất bị xói mòn và chuyển vận trên từng sườn dốc
xuống dòng nước sông suối.
Trong những năm gần đây, xu hướng sử dụng mô hình trên thế giới ngày
càng hướng về các mô hình động thái quá trình vật lý (Siepel et al., 2002).
Những mô hình này dựa trên việc giải quyết các phương trình toán học biểu
diễn bản chất từng quá trình thành phần của hiện tượng xói mòn đất, trong khi
tuân thủ định luật bảo toàn vật chất và năng lượng (Morgan, 1995). Một số mô
hình động thái vật lý bị đánh giá là quá phức tạp, với rất nhiều thông số và đòi
hỏi khối lượng lớn dữ liệu đầu vào. Tuy vậy người ta cho rằng hiện nay những

mô hình này có thể ước tính được lượng đất xói mòn chính xác hơn trong khi
cần ít công sức hiệu chỉnh và kiểm chứng mô hình hơn (Siepel et al., 2002).
Mô hình SWAT
Mô hình SWAT mô phỏng dòng chảy tràn mặt đất theo phương pháp
đường cong SCS; sử dụng phương trình MUSLE để tính lượng đất bị xói mòn.
Lượng bồi lắng dựa trên ước tính về lực cắt, vận tốc dòng chảy và thành phần cơ
giới của đất. Công cụ đánh giá đất và nước SWAT được phát triển từ mô hình
SWRRB nhưng nhằm áp dụng cho lưu vực lớn và địa hình phức tạp. Nó chia
diện tích thành các lưới ô để có thể biểu thị được những biến động theo không
gian của các thông số mô hình. Giống như SWRRB và SPUR models, nó cũng sử
dụng công thức phổ cập cải tiến (Modified USLE) để tính toán dòng bùn cát
(Doe và cộng sự, 1999).
1.1.3. Công thức USLE và các yếu tố ảnh hưởng đến xói mòn đất
Phương trình phổ cập hay tổng hợp xói mòn đất (Universal Soil Loss
Equation, USLE) là một mô hình thực nghiệm đơn giản do Bộ nông nghiệp Mỹ
(USDA) phát triển và liên tục cải tiến từ những năm 1940, có mục đích ước tính
lượng đất xói mòn cho từng ô ruộng mà sau này được cải tiến thành công thức
13