Nghiên cứu biến động lớp phủ rừng trong quá trình phát triển kinh tế - xã hội tỉnh Hòa Bình với sự trợ giúp của công nghệ viễn thám và GIS

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN NGUYỄN THỊ THÚY HẠNH NGHIÊN CỨU BIẾN ĐỘNG LỚP PHỦ RỪNG TRONG QUÁ TRÌNH PHÁT TRIỂN KINH TẾ - XÃ HỘI TỈNH HỊA BÌNH VỚI SỰ TRỢ GIÚP CỦA CÔNG NGHỆ VIỄN THÁM VÀ GIS LUẬN ÁN TIẾN SĨ ĐỊA LÝ HÀ NỘI, 2016 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN NGUYỄN THỊ THÚY HẠNH NGHIÊN CỨU BIẾN ĐỘNG LỚP PHỦ RỪNG TRONG QUÁ TRÌNH PHÁT TRIỂN KINH TẾ - XÃ HỘI TỈNH HỊA BÌNH VỚI SỰ TRỢ GIÚP CỦA CƠNG NGHỆ VIỄN THÁM VÀ GIS CHUYÊN NGÀNH: BẢN ĐỒ, VIỄN THÁM VÀ HỆ THÔNG TIN ĐỊA LÝ Mã số: 62 44 02 14 LUẬN ÁN TIẾN SĨ ĐỊA LÝ NGƢỜI HƢỚNG DẪN KHOA HỌC: PGS.TS.Nguyễn Ngọc Thạch TS Trần Vân Anh HÀ NỘI - 2016 LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu riêng Các số liệu, kết nêu luận án trung thực chƣa đƣợc cơng bố cơng trình khác Tác giả luận án Nguyễn Thị Thúy Hạnh Nguyễn Thị Thúy Hạnh iii LỜI CẢM ƠN Sau nhiều năm học tập nghiên cứu trƣờng Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, đến luận án đƣợc hoàn thành Tác giả xin trân trọng cảm ơn Đại học Quốc gia Hà Nội, trƣờng Đại học Khoa học Tự nhiên quan tâm, giúp đỡ nghiên cứu sinh suốt trình đào tạo Tác giả xin chân thành cảm ơn Ban lãnh đạo, thày cô giáo, cán Khoa Địa lý, Bộ môn Bản Đồ-Viễn thám GIS trực tiếp giảng dạy, trao đổi đóng góp nhiều ý kiến quý báu nhƣ động viên tinh thần để tác giả hoàn thành luận án Tác giả xin trân trọng cảm ơn Ban Lãnh đạo trƣờng Đại học Tài nguyên Môi trƣờng Hà Nội, Ban lãnh đạo khoa Trắc địa - Bản đồ - nơi tác giả cơng tác tồn thể đồng nghiệp, bạn bè tạo điều kiện giúp đỡ, chia sẻ công việc để tác giả hồn thành luận án Tác giả xin bày tỏ lòng biết ơn sâu sắc tới PGS.TS Nguyễn Ngọc Thạch TS.Trần Vân Anh tận tình hƣớng dẫn suốt trình nghiên cứu động viên, tạo điều kiện giúp đỡ để nghiên cứu sinh hoàn thành luận án Tác giả xin gửi lời cảm ơn tới PGS.TS Đinh Thị Bảo Hoa, PGS.TS Phạm Văn Cự, TS Hoàng Thị Thu Hƣơng - Khoa Địa lý- Trƣờng Đại học Khoa học Tự nhiên TS.Vũ Kim Chi - Viện Việt Nam học Khoa học phát triển - Đại học Quốc Gia Hà Nội nhiều Nhà khoa học khác - ngƣời tiên phong nghiên cứu lớp phủ, biến động lớp phủ để đến hôm tác giả đƣợc kế thừa tiếp tục phát triển hƣớng nghiên cứu, đặc biệt gợi ý từ khóa luận, luận án góp phần bổ sung, hồn thiện luận án Tác giả xin gửi lời cảm ơn tới nhà khoa học thuộc quan: Viện Địa lý, Viện Công nghệ Vũ trụ - Viện Hàn lâm Khoa học Công nghệ Việt Nam; Đại học Sƣ phạm Hà Nội; Đại học Mỏ - Địa chất Hà Nội, Đại học Lâm nghiệp,… đóng góp ý kiến từ nhiều khía cạnh để luận án trở nên hoàn thiện Tác giả xin gửi lời cảm ơn tới Sở Tài nguyên Mơi trƣờng Hịa Bình, Cục Thống kê Hịa Bình, Trung tâm Viễn thám Địa chất – Liên đoàn Bản đồ Địa chất miền Bắc, …đã cung cấp tƣ liệu cho luận án Cuối nhƣng không phần quan trọng, tác giả xin bày tỏ lòng biết ơn tới gia đình, đặc biệt ngƣời bạn đời ln sát cánh, chia sẻ động viên tinh thần suốt nhiều năm tác giả thực luận án Hà nội, ngày 10 tháng năm 2016 Tác giả Nguyễn Thị Thúy Hạnh iv MỤC LỤC MỤC LỤC v DANH MỤC HÌNH viii MỞ ĐẦU 1 Tính cấp thiết đề tài Mục tiêu nội dung nghiên cứu Phạm vi nghiên cứu Điểm luận án Luận điểm bảo vệ Cơ sở tài liệu thực đề tài Ý nghĩa khoa học thực tiễn Cấu trúc luận án CHƢƠNG 1: CƠ SỞ LÝ LUẬN VÀ PHƢƠNG PHÁP NGHIÊN CỨU LỚP PHỦ VÀ LỚP PHỦ RỪNG 1.1 Tổng quan cơng trình nghiên cứu có liên quan 1.1.1 Tổng quan tình hình sử dụng tƣ liệu viễn thám nghiên cứu lớp phủ lớp phủ rừng giới 1.1.2 Tình hình nghiên cứu lớp phủ lớp phủ rừng Việt Nam 11 1.1.3 Các cơng trình nghiên cứu liên quan đến lớp phủ lớp phủ rừng Hịa Bình 13 1.2 Cơ sở lý luận nghiên cứu lớp phủ rừng 17 1.2.1 Khái niệm lớp phủ rừng hệ thống phân loại lớp phủ 17 1.2.2 Các phƣơng pháp chiết tách thông tin từ ảnh viễn thám 21 1.2.3 Các phƣơng pháp nghiên cứu biến động viễn thám GIS 23 1.2.4 Các mô hình nghiên cứu biến động 27 1.3 Quan điểm, phƣơng pháp quy trình nghiên cứu 35 1.3.1 Quan điểm nghiên cứu 35 1.3.2 Phƣơng pháp nghiên cứu 36 1.3.3 Quy trình nghiên cứu 38 CHƢƠNG 2: HIỆN TRẠNG LỚP PHỦ RỪNG TỈNH HỊA BÌNH QUA CÁC GIAI ĐOẠN PHÁT TRIỂN KINH TẾ - XÃ HỘI 40 2.1 Các nhân tố ảnh hƣởng đến lớp phủ rừng 40 2.1.1 Vị trí địa lý 40 2.1.2 Điều kiện tự nhiên 41 v 2.1.3 Điều kiện kinh tế xã hội 49 2.2 Tích hợp viễn thám GIS nghiên cứu trạng lớp phủ rừng qua giai đoạn 67 2.2.1 Nghiên cứu thay đổi trạng thái lớp phủ rừng 68 2.2.2 Nghiên cứu phân bố không gian biến động lớp phủ rừng 74 CHƢƠNG ĐÁNH GIÁ BIẾN ĐỘNG LỚP PHỦ RỪNG TỈNH HÒA BÌNH VÀ ĐỀ XUẤT ĐỊNH HƢỚNG PHÁT TRIỂN-BẢO VỆ 87 3.1 Đặc điểm xu hƣớng biến động 87 3.2 Phân tích mối quan hệ biến động lớp phủ rừng yếu tố địa lý tự nhiên, kinh tế xã hội 91 3.2.1 Kết phân tích rừng suy thoái rừng 97 3.2.2 Kết phân tích tái sinh rừng tự nhiên 103 3.2.3 Kết phân tích mở rộng rừng trồng 108 3.3 Kiểm chứng thảo luận yếu tố gây biến động lớp phủ rừng qua hai giai đoạn nghiên cứu 113 3.3.1 Các yếu tố tác động trực tiếp đến biến động lớp phủ rừng qua hai giai đoạn nghiên cứu 113 3.3.2 Kết điều tra điểm chìa khóa 114 3.3.3 Biến động lớp phủ rừng sách phát triển kinh tế xã hội 120 3.4 Dự báo biến động lớp phủ rừng định hƣớng phát triển-bảo vệ lớp phủ rừng đến năm 2020 123 3.4.1 Dự báo biến động lớp phủ rừng đến năm 2020 123 3.4.2 Định hƣớng phát triển bảo vệ lớp phủ rừng tỉnh Hòa Bình đến năm 2020 128 KẾT LUẬN VÀ KIẾN NGHỊ 135 DANH MỤC CƠNG TRÌNH KHOA HỌC CỦA TÁC GIẢ 137 LIÊN QUAN ĐẾN LUẬN ÁN 137 TÀI LIỆU THAM KHẢO 139 PHỤ LỤC 151 Phụ lục 1: Tƣ liệu sử dụng 151 Phụ lục 2: Mã hóa biến giải thích 184 Phụ lục 3: Bản đồ tƣơng ứng với biến giải thích 185 Phụ lục 4: Kết chạy mơ hình hồi quy 186 vi DANH MỤC BẢNG Bảng 1 Những lớp đối tƣợng nghiên cứu luận án 21 Bảng Số liệu thống kê mật độ dân số tỉnh Hòa Bình 50 Bảng 2 Tăng trƣởng kinh tế tỉnh Hịa Bình………………………………….54 Bảng Cơ cấu giá trị sản xuất công nghiệp phân theo ngành…………… 55 Bảng Một số số tăng trƣởng kinh tế ngành nông nghiệp………… 57 Bảng Kết thực số số phát triển ngành lâm nghiệp……….58 Bảng Số lƣợng trƣờng học huyện 62 Bảng Số lƣợng học sinh giáo viên huyện 62 Bảng Thay đổi sử dụng đất qua năm 63 Bảng Chỉ số NDVI biến động số NDVI 70 Bảng 10 Thống kê diện tích loại lớp phủ 77 Bảng 11 Độ xác kết phân loại năm 1994 82 Bảng 12 Độ xác kết phân loại năm 2005 82 Bảng 13 Độ xác kết phân loại năm 2015 82 Bảng Gộp nhóm loại biến động………………………………… 92 Bảng Các biến phụ thuộc 95 Bảng 3 Các biến độc lập 95 Bảng Mô hình rừng suy thối rừng 1994-2005 97 Bảng So sánh đồ xác suất rừng suy thoái rừng với đồ biến động thực tế 99 Bảng Mơ hình rừng suy thối rừng giai đoạn 2005-2015 100 Bảng Mô hình tái sinh rừng tự nhiên giai đoạn 1994-2005 103 Bảng So sánh đồ xác suất tái sinh rừng tự nhiên đồ tái sinh rừng thực tế………………………………………………………………………….108 Bảng Mô hình tái sinh rừng tự nhiên giai đoạn 2005-2015 106 Bảng 10 Mơ hình mở rộng rừng trồng giai đoạn 1994-2005 108 Bảng 11 So sánh đồ xác suất mở rộng rừng trồng với BĐBĐ thực tế 110 Bảng 12 Mơ hình mở rộng rừng trồng giai đoạn 2005-2015 110 vii DANH MỤC HÌNH Hình 1 Phân tích vector biến động………………………………………… 25 Hình Nguyên tắc kỹ thuật trừ ảnh 25 Hình Diện tích dƣới đƣờng cong ROC 35 Hình Bản đồ hành tỉnh Hịa Bình 40 Hình 2 Mơ hình số độ cao tỉnh Hịa Bình 42 Hình Ảnh vệ tinh Landsat TM band 5, band tổ hợp band 432 42 Hình Bản đồ địa mạo tỉnh Hịa Bình 43 Hình Bản đồ sinh khí hậu tỉnh Hịa Bình 47 Hình Mạng lƣới sơng suối tỉnh Hịa Bình 46 Hình Bản đồ thổ nhƣỡng tỉnh Hịa Bình 48 Hình 2 Cơ cấu lao động theo ngành kinh tế năm 1994, 2001 2011……….56 Hình Cơ cấu ngành kinh tế năm 1994, 2000, 2010 2014 54 Hình 10 Hệ thống đƣờng giao thơng tỉnh Hịa Bình 60 Hình Phƣơng pháp trƣớc phân loại nghiên cứu biến động trạng thái lớp phủ rừng……………………………………………………………… .71 Hình 12 Ảnh NDVI qua năm 69 Hình 13 Histogram ảnh NDVI 1994-2005 (trái) 2005-2015 (phải) 70 Hình 14 Biểu đồ biến động trạng thái lớp phủ 1994-2005-2015 71 Hình 15 Biến động số NDVI giai đoạn 1994 - 2005 72 Hình 16 Biến động số NDVI giai đoạn 2005 - 2015 73 Hình 17 Nghiên cứu biến động lớp phủ rừng theo phƣơng pháp sau phân loại 74 Hình 18 Biểu đồ phân bố pixel không gian phổ 75 Hình 19 Kết phân loại ảnh 76 Hình 20 Biểu đồ cấu diện tích loại lớp phủ 76 Hình 21 Ảnh phân loại sau gộp nhóm 78 Hình 22 Biến động lớp phủ rừng khu Bảo tồn Phu Canh 78 Hình 23 Biến động lớp phủ rừng khu vực lịng hồ sơng Đà 79 Hình 24 Biến động lớp phủ rừng khu bảo tồn Thƣợng Tiến 79 Hình 25 Biến động lớp phủ rừng khu vực Lạc Thủy 80 viii Hình 26 Biến động lớp phủ rừng khu vực Lƣơng Sơn 80 Hình 27 Phân bố khơng gian lớp phủ rừng năm 1994 83 Hình 28 Phân bố khơng gian lớp phủ rừng năm 2005 84 Hình 29 Phân bố khơng gian lớp phủ rừng năm 2015 85 Hình Diện tích biến động theo khoảng cách đến đƣờng giao thơng giai đoạn 1994-2005 …………………………………………………………………87 Hình Diện tích biến động theo khoảng cách đến đƣờng giao thông phụ giai đoạn 1994-2005 87 Hình 3 Diện tích biến động theo khoảng cách đến sông suối lớn giai đoạn 1994-2005 88 Hình Diện tích biến động theo khoảng cách đến sông suối nhỏ giai đoạn 1994-2005 88 Hình Diện tích biến động theo độ cao địa hình giai đoạn 1994-2005 88 Hình Diện tích biến động theo độ dốc giai đoạn 1994-2005 89 Hình Diện tích biến động theo khoảng cách đến đƣờng giao thơng giai đoạn 2005-2015 89 Hình Diện tích biến động theo khoảng cách đến đƣờng giao thông phụ giai đoạn 2005-2015 89 Hình Diện tích biến động theo khoảng cách đến sông suối lớn giai đoạn 2005-2015 90 Hình 10 Diện tích biến động theo khoảng cách đến sông suối nhỏ giai đoạn 2005-2015 90 Hình 11 Diện tích biến động theo độ cao địa hình giai đoạn 2005-2015 90 Hình 12 Diện tích biến động theo độ dốc giai đoạn 2005-2015 91 Hình 13 Bản đồ biến động LPR 93 Hình 14 Bản đồ biến động LPR 94 Hình 15 Sơ đồ bƣớc phân tích mối quan hệ biến động lớp phủ rừng yếu tố địa lý 96 Hình 16 Xác suất rừng suy thoái rừng giai đoạn 1994-2005 98 Hình 17 Đƣờng cong ROC kiểm chứng mơ hình rừng suy thối rừng 99 Hình 18 Xác suất rừng suy thoái rừng giai đoạn 2005-2015 102 Hình 19 Xác suất tái sinh rừng tự nhiên giai đoạn 1994 - 2005 104 ix Hình 20 Đƣờng cong ROC kiểm chứng mơ hình tái sinh rừng tự nhiên 105 Hình 21 Xác suất tái sinh rừng tự nhiên giai đoạn 2005-2015 107 Hình 22 Xác suất mở rộng rừng trồng giai đoạn 1994-2005 109 Hình 23 Đƣờng cong ROC kiểm chứng mơ hình mở rộng rừng trồng 110 Hình 24 Xác suất mở rộng rừng trồng giai đoạn 2005-2015 112 Hình 25 Vị trí xã hộ gia đình đƣợc vấn 115 Hình 26 Vị trí xã hộ gia đình đƣợc vấn Huyện Đà Bắc 115 Hình 27 Rừng keo đƣợc trồng thơn Rổng Vịng, xã Lâm Sơn, H Lƣơng Sơn 116 Hình 28 Ngơi nhà đƣợc xây dựng khang trang tiền thu nhập từ trồng rừng thơn Rổng Vịng, xã Lâm Sơn, H Lƣơng Sơn 116 Hình 29 Ngơi nhà ngƣời Kinh thơn Bình Tiến, xã Dân Hạ, H Kỳ Sơn 117 Hình 30 Bên nhà ngƣời Dao xã Hiền Lƣơng, Đà Bắc 118 Hình 31 Cơ cấu số hộ gia đình theo nguồn thu nhập nhóm 119 Hình 32 Những ngơi nhà ngƣời Tày, xã Tân Minh, H Đà Bắc 120 Hình 33 Những đống củi đƣợc hộ gia đình ngƣời Tày, xã Tân Minh 120 Hình 34 Bản đồ dự báo xác suất rừng suy thối rừng đến năm 2020 125 Hình 35 Bản đồ dự báo xác suất tái sinh rừng tự nhiên đến năm 2020 126 Hình 36 Bản đồ dự báo xác suất mở rộng rừng trồng đến năm 2020 127 Hình Bản đồ phân vùng địa lý tự nhiên…………………………………131 x STT Tọa độ Huyện Xã Thơn, xóm Thơng tin lớp phủ Y (m) X (m) 160 2262668 569178 Yên Thủy Bảo Hiệu 161 2277169 547921 Lạc Sơn Q Hịa Xóm Vệ Đất trống, bụi 162 2266187 558284 Lạc Sơn Hình Hẻm Xóm Khƣơng RT chƣa khép tán 163 2261004 547646 Lạc Sơn Tân Mỹ Xóm Mặc Dân cƣ 164 2255410 551311 Lạc Sơn Ân Nghĩa Xóm Nghéo RT khép tán Đ 165 2255393 568043 n Thủy Phú Lai Xóm Rị Mặt nƣớc Đ 166 2268011 540044 Lạc Sơn Thƣợng Cốc Xóm Trang Nơng nghiệp Đ 167 2269337 577029 Lạc Thủy Liên Hồng Đồng Xuống RT chƣa khép tán Đ 168 2268430 559581 n Thủy Lạc Sỹ Xóm Sỹ Nơng nghiệp 169 2268318 559485 Yên Thủy Lạc Sỹ Xóm Sỹ RT khép tán 170 2267720 559471 Yên Thủy Lạc Sỹ Xóm Sỹ RTN nghèo 171 2266389 561536 Yên Thủy Lạc Lƣơng Đất Bãi Đất trống, bụi 172 2265380 568331 Yên Thủy Lạc Hƣng Áng RT chƣa khép tán 173 2265114 571352 Lạc Thủy Đồng Môn Cú Đẻ Dân cƣ 174 2273983 533843 Lạc Sơn Mỹ Thành Cỏ Giữa RT khép tán Đ 175 2272530 552416 Lạc Sơn Mỹ Thành Riệc RTN nghèo Đ 176 2259012 550894 Lạc Sơn Tân Mỹ Xóm Khí Đất trống, bụi 175 Thực địa APL RTN nghèo Đ Đ Đ STT Tọa độ Huyện Xã Thơn, xóm Thơng tin lớp phủ Y (m) X (m) 177 2261097 588514 Lạc Thủy Đồng Tâm Đồng Mối Đất trống, bụi 178 2268846 562616 Yên Thủy Lạc Sỹ Tháu Nông nghiệp 179 2268632 561876 Yên Thủy Lạc Sỹ Tháu Dân cƣ 180 2268635 560961 Yên Thủy Lạc Sỹ Tháu RT khép tán Đ 181 2266683 568318 Yên Thủy Lạc Hƣng Bông Bạc RTN nghèo Đ 182 2263588 548568 Lạc Sơn Liên Vũ Xóm Lả Đất trống, bụi Đ 183 2259342 553887 Lạc Sơn Ân Nghĩa Xóm Re Nơng nghiệp 184 2309853 537512 RT chƣa khép tán Đ 185 2308722 535814 Mặt nƣớc Đ 186 2308096 535860 Nông nghiệp 187 2306179 536325 Đất trống, bụi 188 2303684 537082 RT chƣa khép tán 189 2303346 537026 Nông nghiệp 190 2303252 536746 Dân cƣ 191 2305003 534484 RT khép tán 192 193 2305647 533712 176 Thực địa APL Đ Đ RTN nghèo Đ Đất trống, bụi Đ STT Tọa độ Y (m) X (m) Huyện Xã 194 Thơn, xóm Thơng tin lớp phủ Thực địa APL Dân cƣ 195 2305833 533157 RT khép tán Đ 196 2305809 532689 Mặt nƣớc Đ 197 Nông nghiệp 198 2306694 530083 Đất trống, bụi 199 2308012 528672 RT chƣa khép tán Đ 200 2308471 526756 Dân cƣ Đ 201 2309118 525435 RT khép tán Đ 202 2308482 523850 RTN nghèo 203 2307579 521114 Đất trống, bụi 204 2307281 518820 Nông nghiệp 206 2308815 516755 RTN nghèo 207 2310655 515982 RTN giàu S 208 2313186 515572 Đất trống, bụi Đ 209 2314275 513374 Nông nghiệp 210 2314539 511992 Dân cƣ Đ 211 2303448 534534 RT khép tán Đ 177 Đ STT Tọa độ Y (m) X (m) 212 2299637 213 Huyện Xã Thơn, xóm Thơng tin lớp phủ Thực địa APL 535290 Mặt nƣớc Đ 2298382 534382 RTN nghèo Đ 214 2297166 534276 Nông nghiệp 215 2295919 534609 RT chƣa khép tán 216 2294449 534516 RT khép tán 217 2292493 533446 RTN nghèo 218 2290634 533688 Đất trống, bụi Đ 219 2286499 532105 Dân cƣ Đ 220 2284880 530551 RT khép tán 221 2283245 529391 Đất trống, bụi 222 2281506 528703 Nông nghiệp 223 2279913 525927 RTN nghèo 224 2281644 522663 RT khép tán 225 2282028 518552 Mặt nƣớc Đ 226 2283778 515670 RTN giàu Đ 227 2284598 514804 Nông nghiệp 228 2285864 512827 Đất trống, bụi 178 Đ Đ Đ STT Tọa độ Y (m) X (m) 229 2286491 230 Huyện Xã Thơn, xóm Thơng tin lớp phủ Thực địa APL 508711 Nông nghiệp Đ 2290083 506741 Dân cƣ 231 2290522 502506 RT khép tán 232 2291205 499637 RTN nghèo 233 2291094 496620 RTN giàu 234 2291706 496131 RTN nghèo 235 2292287 494344 Đất trống, bụi Đ 236 2294728 491589 Nông nghiệp Đ 237 2293927 535060 RTN nghèo Đ 238 2293291 536773 RT khép tán S 239 2293347 539063 Mặt nƣớc Đ 240 2293513 541701 Nông nghiệp 241 2293421 543506 Nông nghiệp 242 2293314 546412 RTN nghèo 243 2292492 548104 Nông nghiệp 244 2289981 550047 RT chƣa khép tán 245 2289090 551430 RT khép tán 179 Đ Đ Đ STT Tọa độ Y (m) X (m) 246 2288586 247 Huyện Xã Thơn, xóm Thơng tin lớp phủ Thực địa APL 552672 RTN giàu Đ 2287037 554680 Đất trống, bụi 248 2285135 556582 Nông nghiệp 249 2284722 558188 Dân cƣ 250 2281033 560076 RT khép tán S 251 2278419 562324 Mặt nƣớc Đ 252 2276102 564560 Nông nghiệp 253 2276492 567956 Nông nghiệp 254 2277688 568855 RT chƣa khép tán 255 2277814 570067 RTN nghèo 256 2279564 569054 Đất trống, bụi Đ 257 2281222 568224 Nông nghiệp Đ 258 2282742 568492 RTN nghèo Đ 259 2286349 568842 RT khép tán S 260 2287766 568844 Dân cƣ 261 2298476 548426 Đất trống, bụi 262 2289021 551446 Dân cƣ 180 Đ Đ Đ STT Tọa độ Y (m) X (m) 263 2284805 264 Huyện Xã Thơn, xóm Thơng tin lớp phủ Thực địa APL 557949 RTN giàu Đ 2293838 543012 Nông nghiệp 265 2285493 558961 RTN nghèo 266 2282489 559208 RTN nghèo 267 2277281 570189 Nông nghiệp 268 2293458 544648 Nông nghiệp 269 2295800 540868 RT chƣa khép tán 181 S Thống kê điểm kiểm tra: 130 số 269 điểm thực địa đƣợc sử dụng để kiểm chứng kết phân loại ảnh Kết đƣợc tổng hợp bảng dƣới: Bảng Kết kiểm tra thực địa Lớp phủ Tổng số điểm kiểm tra Số điểm xác APL RTN giàu 11 9/11 RTN nghèo 18 15/18 RT khép tán 22 17/22 RT chƣa khép tán 23 19/23 Nông nghiệp 21 16/21 Đất trống, bụi 12 9/12 Dân cƣ 14 11/14 Mặt nƣớc 9/9 Tổng số điểm 130 xác 105/130 Ghi chú: Đ: đúng, S: sai APL: kết phân loại 1.3 Một số hình ảnh thực địa Phỏng vấn ngƣời dân Đà Bắc (6/2015) 182 Phỏng vấn ngƣời dân Thung Rếch (Kim Bắc – Kim Bôi) 3/2015 Hoạt động khai thác đất đá đoạn gần cầu Hòa Bình (6/2015) Thị Trấn Mai Châu (8/2014) Lạc Thủy 183 Phụ lục 2: Mã hóa biến giải thích Bảng Mã đồ thổ nhưỡng Đất phù sa Đất loại đá khác Bãi bồi cát Đất thƣờng xuyên đƣợc bồi trung tính Đất thƣờng xuyên đƣợc bồi chua Đất thƣờng xuyên đƣợc bồi glây Đất phù sa sông suối Đất dốc tụ Đất bạc màu 10 Đất đen sản phẩm phong hóa đá vơi 11 Đất feralit biến đổi trồng lúa 12 Đất feralit phù sa cổ 13 Đất feralit đá macma bazơ 14 Đất feralit đá vôi 15 Đất feralit đá sét, phiến sét 16 Đất Feralit đá sét xen đá vôi 17 Đất feralit đá cát 18 Đất feralit đá macma axit 19 Đất feralit đá biến chất 20 Đất feralit mùn đá macma bazơ 21 Đất feralit mùn đá vôi 22 Đất Feralit mùn đá sét xen đá vôi 23 Đất Feralit mùn đá macma axit 24 Núi đá vơi có lớp phủ rừng 25 Núi đá vơi Bảng 10 Mã đồ địa mạo Núi kiến trúc bóc mịn thể macma xâm nhập Núi kiến trúc bóc mịn Núi, sơn ngun khối tảng phận võng phức Dải núi thấp uốn nếp khối tảng Núi uốn nếp khối tảng phá hủy phận Núi kiến tạo xâm thực khối tảng Dải núi uốn nếp khối tảng Đồng đồi núi Đồng thung lũng trung bình - cao 184 10 Đồng Karst xâm thực 11 Đồng thung lũng cao 12 Đồng thấp ẩm ƣớt 13 Thung lũng hình chữ V 14 Thung lũng tích tụ 15 Thung lũng địa hào 16 Dải núi đá vôi võng phức 17 Các đồi khối tảng uốn nếp Bảng 11 Mã đồ khí hậu Mã Kí hiệu Nhiệt độ TB Lƣợng mƣa Độ dài mùa Độ dài mùa năm (0C) TB năm lạnh khô (tháng) (mm) IA1a ≥20 ≥2500 ≤3 ≤2 IB1b ≥20 2000-2500 ≤3 3-4 IB1c ≥20 2000-2500 ≤3 5-6 IC1c ≥20 1500-2000 ≤3 5-6 IIA1a 16- 20 ≥2500 ≤3 ≤2 IIA1b 16- 20 ≥2500 ≤3 3-4 IIB1b 16- 20 2000-2500 ≤3 3-4 IIIA2a 12-16 ≥2500 4-6 ≤2 IIIA2b 12-16 ≥2500 4-6 3-4 10 IIIB2b 12-16 2000-2500 4-6 3-4 11 IVB3b ≤12 2000-2500 7-9 3-4 Phụ lục 3: Bản đồ tƣơng ứng với biến giải thích 3.1 Các đồ kinh tế xã hội năm 2001 (Nguồn: Số liệu điều tra Nông nghiệp nông thôn cấp xã năm 2001) 185 Phụ lục 4: Kết chạy mơ hình hồi quy 4.1 Giai đoạn 1994-2005 4.1.1 Mất rừng suy thoái rừng Omnibus Tests of Model Coefficients Chi-square df Sig, 712,094 75 ,000 Block 712,094 75 ,000 Model 712,094 75 ,000 Step Step Model Summary Step -2 Log Cox & Snell R Nagelkerke R likelihood Square Square 5,526a ,072 ,152 a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, Classification Tablea Predicted Defor Observed Step Defor 1 Percentage Correct 8563 11 99,9 924 10 1,1 Overall Percentage 90,2 a, The cut value is ,500 4.1.2 Tái sinh rừng tự nhiên Omnibus Tests of Model Coefficients Chi-square df Sig, 1915,575 75 ,000 Block 1915,575 75 ,000 Model 1915,575 75 ,000 Step Step Model Summary Step -2Log likelihood Cox & Snell R Square Nagelkerke R Square 2,747a ,182 ,262 a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, 186 Classification Tablea Predicted Natu Percentage Correct 6251 572 91,6 1784 901 33,6 Observed Step Natu Overall Percentage 75,2 a, The cut value is ,500 4.1.3 Mở rộng rừng trồng Omnibus Tests of Model Coefficients Chi-square df Sig, 809,630 75 ,000 Block 809,630 75 ,000 Model 809,630 75 ,000 Step Step Model Summary Step -2 Log Cox & Snell R Nagelkerke R likelihood Square Square 6,093a ,082 ,169 a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, Classification Tablea Predicted Plan Percentage Correct 8514 24 99,7 961 ,9 Observed Step Plan Overall Percentage 89,6 a, The cut value is ,500 187 4.2 Giai đoạn 2005-2015 4.2.1 Mất rừng suy thoái rừng Omnibus Tests of Model Coefficients Chi-square df Sig, 110,902 76 ,006 Block 110,902 76 ,006 Model 110,902 76 ,006 Step Step Model Summary Step -2 Log Cox & Snell R Nagelkerke R likelihood Square Square 3,487a ,012 ,024 a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, Classification Tablea Predicted Defor Observed Step Defor 1 Percentage Correct 8569 100,0 939 ,1 Overall Percentage 90,1 a, The cut value is ,500 4.2.2 Tái sinh rừng tự nhiên Omnibus Tests of Model Coefficients Chi-square df Sig, 2285,650 76 ,000 Block 2285,650 76 ,000 Model 2285,650 76 ,000 Step Step Model Summary Step -2 Log Cox & Snell R Nagelkerke R likelihood Square Square 57,500a ,214 ,307 a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, 188 Classification Tablea Predicted Natu Percentage Correct 5896 915 86,6 1839 859 31,8 Observed Step Natu Overall Percentage 71,0 a, The cut value is ,500 4.2.3 Mở rộng rừng trồng Omnibus Tests of Model Coefficients Chi-square df Sig, 4094,718 74 ,000 Block 4094,718 74 ,000 Model 4094,718 74 ,000 Step Step Model Summary Step -2 Log Cox & Snell R Nagelkerke R likelihood Square Square 6,445a ,351 ,682 a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, Classification Tablea Predicted Plan Observed Step Plan 0 8149 390 Percentage Correct 201 97,6 720 64,9 Overall Percentage 93,8 a, The cut value is ,500 189