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 ận Lu án NGHIÊN CỨU BIẾN ĐỘNG LỚP PHỦ RỪNG tiê TRONG QUÁ TRÌNH PHÁT TRIỂN KINH TẾ - XÃ HỘI n sĩ TỈNH HỊA BÌNH VỚI SỰ TRỢ GIÚP CỦA CƠNG NGHỆ ịa Đ VIỄN THÁM VÀ GIS lý 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 ận Lu NGHIÊN CỨU BIẾN ĐỘNG LỚP PHỦ RỪNG án TRONG QUÁ TRÌNH PHÁT TRIỂN KINH TẾ - XÃ HỘI tiê TỈNH HỊA BÌNH VỚI SỰ TRỢ GIÚP CỦA CÔNG NGHỆ VIỄN THÁM VÀ GIS n sĩ CHUYÊN NGÀNH: BẢN ĐỒ, VIỄN THÁM VÀ HỆ THÔNG TIN ĐỊA LÝ ịa Đ Mã số: 62 44 02 14 lý 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 ận Lu Nguyễn Thị Thúy Hạnh án n tiê Nguyễn Thị Thúy Hạnh sĩ ịa Đ lý iii ận Lu 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 toà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, hoà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 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ả án n tiê sĩ ịa Đ lý 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 Lu Ý nghĩa khoa học thực tiễn ận Cấu trúc luận á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 n tiê 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 sĩ 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 ịa Đ 1.2 Cơ sở lý luận nghiên cứu lớp phủ rừng 17 lý 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 Lu 3.2.1 Kết phân tích rừng suy thối rừng 97 3.2.2 Kết phân tích tái sinh rừng tự nhiên 103 ận 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 án giai đoạn nghiên cứu 113 n tiê 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 sĩ ịa Đ 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 lý 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 Lu Bảng Chỉ số NDVI biến động số NDVI 70 ận 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 án 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 tiê Bảng Gộp nhóm loại biến động………………………………… 92 n Bảng Các biến phụ thuộc 95 sĩ 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 ịa Đ 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 lý Bảng Mơ hình rừng suy thoá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 Lu Hình Mạng lƣới sơng suối tỉnh Hịa Bình 46 ận 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 án 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 n tiê 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 sĩ Hình 13 Histogram ảnh NDVI 1994-2005 (trái) 2005-2015 (phải) 70 ịa Đ 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 lý 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 Lu 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 ận 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 án 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 tiê n 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 sĩ 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 ịa Đ 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 lý 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 thố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 Lu Hình 29 Ngơi nhà ngƣời Kinh thơn Bình Tiến, xã Dân Hạ, H Kỳ Sơn 117 ận 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 án 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 tiê Hình 34 Bản đồ dự báo xác suất rừng suy thoá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 n Hình 36 Bản đồ dự báo xác suất mở rộng rừng trồng đến năm 2020 127 sĩ Hình Bản đồ phân vùng địa lý tự nhiên…………………………………131 ịa Đ lý x Tọa độ Xã Thơng tin lớp phủ Thơn, xóm 160 2262668 569178 Yên Thủy Bảo Hiệu 161 2277169 547921 Lạc Sơn Quý 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 Yê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 Yên Thủy Lạc Sỹ Xóm Sỹ Nơng nghiệp 169 2268318 559485 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 án n X (m) ận Lu Y (m) ịa Huyện tiê STT APL RTN nghèo Đ sĩ lý 175 Thực địa Đ Đ Đ STT Tọa độ Huyện Xã Thơn, xóm Thơng tin lớp phủ 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 185 2308722 535814 186 2308096 535860 187 2306179 536325 188 2303684 537082 189 2303346 537026 190 2303252 536746 191 2305003 534484 ận Lu Y (m) án Đ RT chƣa khép tán Đ Mặt nƣớc Đ n tiê Nông nghiệp Đ sĩ Đất trống, bụi RT chƣa khép tán ịa Đ Nông nghiệp lý 2305647 APL Dân cƣ RT khép tán 192 193 Thực địa 533712 176 RTN nghèo Đ Đất trống, bụi Đ STT Tọa độ Y (m) X (m) Huyện Xã Thông tin lớp phủ Thơn, xóm Thực địa 194 APL Dân cƣ 2305833 533157 RT khép tán Đ 196 2305809 532689 Mặt nƣớc Đ 197 ận Lu 195 Nông nghiệp 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 203 2307579 521114 204 2307281 518820 206 2308815 516755 207 2310655 515982 208 2313186 515572 209 2314275 513374 Nông nghiệp 210 2314539 511992 Dân cƣ Đ 211 2303448 534534 RT khép tán Đ án 198 tiê n RTN nghèo Đất trống, bụi Đ sĩ Nông nghiệp RTN nghèo ịa RTN giàu S Đất trống, bụi Đ lý 177 Đ STT Tọa độ Y (m) X (m) 212 2299637 213 Huyện Xã Thơng tin lớp phủ Thơn, xóm 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 219 2286499 532105 220 2284880 530551 221 2283245 529391 222 2281506 528703 223 2279913 525927 224 2281644 522663 225 2282028 518552 226 2283778 227 228 ận Lu Thực địa Đ án Đ Dân cƣ Đ tiê Đất trống, bụi n RT khép tán Đất trống, bụi Đ sĩ Nông nghiệp RTN nghèo Đ ịa RT khép tán Đ 515670 RTN giàu Đ 2284598 514804 Nông nghiệp 2285864 512827 Đất trống, bụi lý Mặt nƣớc 178 Đ STT Tọa độ Y (m) X (m) 229 2286491 230 Huyện Xã Thông tin lớp phủ Thơn, xóm 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 236 2294728 491589 237 2293927 535060 238 2293291 536773 239 2293347 539063 240 2293513 541701 241 2293421 543506 242 2293314 546412 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 ận Lu Thực địa án Đ Đất trống, bụi Đ Nông nghiệp Đ RTN nghèo Đ RT khép tán S Mặt nƣớc Đ n tiê Đ sĩ Nông nghiệp ịa Nông nghiệp lý 179 Đ RTN nghèo Đ STT Tọa độ Y (m) X (m) 246 2288586 247 Huyện Xã Thông tin lớp phủ Thơn, xóm 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 253 2276492 567956 254 2277688 568855 255 2277814 570067 256 2279564 569054 257 2281222 568224 258 2282742 568492 259 2286349 568842 260 2287766 568844 Dân cƣ 261 2298476 548426 Đất trống, bụi 262 2289021 551446 Dân cƣ ận Lu Thực địa Đ án Nông nghiệp tiê Nông nghiệp n RT chƣa khép tán RTN nghèo Đ sĩ Đất trống, bụi Đ Nông nghiệp Đ RTN nghèo Đ RT khép tán S ịa lý 180 Đ Đ STT Tọa độ Y (m) X (m) 263 2284805 264 Huyện Xã Thơng tin lớp phủ Thơn, xóm 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 ận Lu Thực địa án RT chƣa khép tán n tiê Đ sĩ ịa lý 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 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 ận Lu RTN giàu Tổng số điểm 130 xác án 1.3 Một số hình ảnh thực địa n tiê Ghi chú: Đ: đúng, S: sai APL: kết phân loại 105/130 sĩ ịa Đ lý 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 ận Lu án tiê Hoạt động khai thác đất đá đoạn gần cầu Hịa Bình (6/2015) n sĩ ịa Đ lý 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 14 16 Đất feralit đá macma bazơ Đất feralit đá vôi Đất feralit đá sét, phiến sét án 15 Đất feralit phù sa cổ ận 13 Lu Đất Feralit đá sét xen đá vôi Đấ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 n tiê 17 sĩ ịa Đ lý 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 ≥20 2000-2500 ≤3 5-6 ≥20 1500-2000 ≤3 5-6 16- 20 ≤3 ≤2 IB1c IC1c IIA1a ≥2500 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 n tiê án ận Lu sĩ ịa Đ 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) lý 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 ,152 Lu ,072 ận a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, án Classification Tablea Predicted Defor 8563 924 99,9 10 1,1 90,2 ịa Đ Overall Percentage 11 sĩ n Step Defor tiê Observed a, The cut value is ,500 4.1.2 Tái sinh rừng tự nhiên lý 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 Percentage Correct 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 Sig, 809,630 75 ,000 Block 809,630 75 ,000 Model 809,630 75 ,000 ận Step Step Lu Chi-square df án Model Summary tiê 6,093a ,082 ,169 n Step -2 Log Cox & Snell R Nagelkerke R likelihood Square Square sĩ ịa Đ a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, lý 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 ận Lu a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, Classification Tablea án Predicted Defor 8569 100,0 939 ,1 n Step Defor tiê Observed Percentage Correct sĩ Overall Percentage 90,1 ịa Đ a, The cut value is ,500 4.2.2 Tái sinh rừng tự nhiên Omnibus Tests of Model Coefficients Sig, 2285,650 76 ,000 Block 2285,650 76 ,000 Model 2285,650 76 ,000 Step Step lý Chi-square df 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 Sig, 4094,718 74 ,000 Block 4094,718 74 ,000 Step Step ận Lu Chi-square df Model 4094,718 ,000 án Model Summary 74 6,445a tiê Step -2 Log Cox & Snell R Nagelkerke R likelihood Square Square ,351 ,682 n a, Estimation terminated at iteration number 20 because maximum iterations has been reached, Final solution cannot be found, sĩ ịa Đ Classification Tablea Predicted Plan Step Plan 0 8149 390 201 97,6 720 64,9 lý Observed Percentage Correct Overall Percentage 93,8 a, The cut value is ,500 189