103 4 Khu du lịch sinh thái Gáo Giồng 169 582 34,5 5 Khu sinh thái Đồng Sen Gò Tháp 75 600 6,3 6 Làng nổi Tân Lập 90 500 121,5 7 Khu bảo tồn sinh thái Đồng Tháp Mười 8 Khu bảo tồn Cây dược liệu Đồng T[.]
4 Khu du lịch sinh thái Gáo Giồng Khu sinh thái Đồng Sen Gò Tháp Làng Tân Lập Khu bảo tồn sinh thái Đồng Tháp Mười Khu bảo tồn Cây dược liệu Đồng Tháp Mười 169.582 75.600 90.500 - 34,5 6,3 121,5 - 3.4.2 Xây dựng lớp liệu áp lực rủi ro lên hệ sinh thái đất ngập nước Đồng Tháp Mười 3.4.2.1 Xây dựng lớp liệu áp lực rủi ro lên khả cung cấp bãi đẻ trứng, vườn ươm, môi trường sống (mức độ đa dạng sinh học) - Xây dựng liệu khơng gian số lồi ngoại lai xâm hại Các số liệu điều tra, khảo sát, thu thập tổng hợp qua tài liệu [3], [7], [72], [84], [126], [127], [134] số loài ngoại lai xâm hại khu vực vườn quốc gia Tràm Chim 12, khu bảo tồn đất ngập nước Láng Sen 03, Xẻo Quýt 06 Các khu vực cịn lại chưa có số liệu cụ thể số loài ngoại lai diện khu vực Dữ liệu số lượng loài ngoại lai xâm hại khơng gian hố cách gán vào lớp liệu hệ sinh thái đất ngập nước sử dụng công cụ Editor ArcGIS Pro 3.0 Kết thể Hình 3.16 Hình 3.16 Dữ liệu mức độ đa dạng sinh học giá trị dịch vụ du lịch đất ngập nước Đồng Tháp Mười năm 2018 103 - Xây dựng liệu độ mặn nước Độ mặn nước áp lực ảnh hưởng tới mục tiêu bảo vệ tính đa dạng HST ĐNN dịch vụ du lịch, giải trí, nghiên cứu khoa học, đại diện độ mặn nước Dữ liệu độ mặn nước trích xuất từ liệu xâm nhập mặn theo kịch biến đổi khí hậu nghiên cứu xây dựng định dạng MapInfo đề tài cấp nhà nước [75] Việc chuyển đổi định dạng, cắt theo ranh giới gán giá trị mức độ xâm nhập mặn vào lớp hệ sinh thái thực phần mềm ArcGIS Pro 3.0 Kết thể Hình 3.17 cho thấy giới hạn xâm nhập mặn ảnh hưởng tới khu bảo tồn sinh thái Đồng Tháp Mười với mức 0,25 g/l, khu vực khác không chịu tác động yếu tố Hình 3.17 Khu vực bị xâm nhập mặn theo kịch biến đổi khí hậu - Xây dựng lớp liệu nhạy cảm cháy Mức độ nhạy cảm cháy vùng đất ngập nước Đồng Tháp Mười đại diện thông qua số lượng điểm cháy diễn Số lượng điểm cháy nhiều mức độ nhạy cảm cháy cao ngược lại Dữ liệu điểm cháy thu thập từ năm 2000-2018 hệ tọa độ WGS84, định dạng *shp, đó, bước việc xây dựng lớp liệu chuyển đổi hệ tọa độ VN2000 Bước thực chức Project ArcGIS Pro 3.0 (Analysis → Tools → Geoprocessing → Data Management Tools → Projections and Transformations → Project) 104 Dữ liệu MODIS VIIRS cung cấp thông tin chất lượng xử lý điểm cháy trường thông tin "Confidence" Dữ liệu điểm cháy xử lý từ ảnh MODIS đánh giá mức độ tin cậy thông qua tỷ lệ %; VIIRS đánh giá thơng qua mức thấp (l), bình thường (n) cao (h) Trong luận án này, điểm có mức độ tin cậy thấp (l) 50% loại bỏ khỏi số liệu trước đưa vào tính tốn chức Select by Attribute ArcGIS Pro 3.0 Hình 3.18 Phân bố điểm cháy khu vực Đồng Tháp Mười giai đoạn 2000 - 2018 Sau loại bỏ điểm mức độ tin cậy thấp, liệu điểm cháy gán vào khu vực đất ngập nước chức Spatial Join, Select by Attribute, Field Calculator Số liệu gán cho thấy, tính từ năm 2000 đến 2018, vườn quốc gia Tràm Chim có 77 điểm cháy (cao nhất), khu bảo tồn đất ngập nước Láng Sen có 21 điểm cháy (cao thứ 2), khu bảo tồn dược liệu Đồng Tháp Mười có 19 điểm cháy (xếp thứ 3), thứ tư Đồng Sen Gị Tháp có 10 điểm cháy, Gáo Giồng có điểm cháy, Tân Lập có điểm cháy (thứ 6), thứ Xẻo Quýt với điểm cháy cuối khu bảo tồn sinh thái Đồng Tháp Mười với số điểm cháy 3.4.2.2 Xây dựng lớp liệu đánh giá áp lực rủi ro lên du lịch giải trí Biên tập, tách lớp liệu đường giao thông, khu công nghiệp, khu dân cư từ sở liệu địa lý tỷ lệ 1:2.000, 1:5.000; khu nuôi trồng thuỷ sản từ đồ trạng sử dụng đất công cụ biên tập; xây dựng lớp liệu khoảng cách từ lớp liệu công cụ Euclidean Distance ArcGIS Pro 3.0 105 Lớp raster khoảng cách từ khu vực chăn nuôi (nuôi trồng thuỷ sản) Lớp raster khoảng cách từ khu vực công nghiệp Lớp raster khoảng cách từ đường giao thông Lớp raster khoảng cách từ khu vực dân cư Hình 3.19 Các lớp raster khoảng cách từ khu vực chăn nuôi, khu vực công nghiệp, khu vực dân cư, đường giao thông ảnh hưởng tới dịch vụ du lịch 3.4.3 Xây dựng lớp liệu tình trạng suy thoái hệ sinh thái đất ngập nước Đồng Tháp Mười 3.4.3.1 Xây dựng lớp liệu mức độ giảm diện tích lớp phủ thực vật Dữ liệu đánh giá mức độ suy giảm diện tích lớp phủ thành lập thông qua việc phân loại lớp phủ ảnh vệ tinh thời điểm 2010 2018 thông qua số NDVI Chồng xếp lớp phủ thực vật phân loại để xác định diện tích suy giảm, theo đó, diện tích thực phủ vườn quốc gia Tràm Chim giảm -383,75 ha, khu di tích Xẻo Quýt giảm -24,30 ha, khu vực khác có diện tích biến động nên xem khơng thay đổi Các bước tính tốn biến động gán số liệu mức độ suy giảm diện tích lớp phủ thực vật vào lớp hệ sinh thái đất ngập nước cho khu vực phần mềm ArcGIS Pro 3.0 3.4.3.2 Xây dựng liệu mức độ giảm số lượng động thực vật nguy cấp, quý Thông tin suy giảm số lượng động thực vật nguy cấp, quý khu vực nghiên cứu thu thập từ tài liệu [3], [7], [72], [84], [126], [127], [134] kết hợp với điều tra, khảo sát người dân, nhà quản lý khu vực Tuy nhiên, số liệu có bị hạn chế, có thông tin mức độ suy giảm số lượng 106 Sếu đầu đỏ khu vực vườn quốc gia Tràm Chim khu bảo tồn Láng Sen Các thông tin gán vào lớp liệu hệ sinh thái phần mềm ArcGIS Pro 3.0 Diễn biến số lượng Sếu đầu đỏ vườn quốc gia Tràm chim từ 1986-2018 thể Hình 3.20 Số cá thể 1200 1052 1000 814 800 741 665 631 605 600 511 503 479 469 400 167 200 68 64 113 128 159 93 89 48 125 126 84 85 94 52 14 21 21 14 11 2018 302 2013 271 2017 2016 2015 2014 2012 2011 2010 2009 2008 2007 2006 2005 2004 2003 2002 2001 2000 1999 1998 1997 1996 1995 1994 1993 1992 1991 1990 1989 1988 1987 1986 Hình 3.20 Diễn biến số lượng Sếu đầu đỏ vườn quốc gia Tràm Chim (nguồn: Vườn quốc gia Tràm Chim, 2018) 3.4.3.3 Xây dựng lớp liệu mức độ ô nhiễm môi trường (chất lượng nước mặt) Dữ liệu mức độ ô nhiễm môi trường nước mặt kế thừa kết nghiên cứu đề tài [113] Dữ liệu thành lập dựa số liệu điểm quan trắc môi trường nước mặt địa bàn ba tỉnh Đồng Tháp, Long An, Tiền Giang vào tháng 11 năm 2018; ảnh vệ tinh Landsat-8 OLI chụp ngày 15/11/2018 [92-94] Hình 3.21 Chất lượng nước mặt khu vực Đồng Tháp Mười tháng 11 năm 2018 [113] 107 3.4.4 Tích hợp liệu phục vụ đánh giá mức độ ưu tiên bảo vệ, phục hồi hệ sinh thái đất ngập nước khu vực Đồng Tháp Mười Do lớp liệu tổng số loài, số loài nguy cấp quý hiếm, số lượng loài ngoại lai xâm hại, độ mặn nước, nhạy cảm cháy, số lượng khách doanh thu từ du lịch trung bình năm, mức độ giảm diện tích lớp phủ thực vật, mức độ suy giảm số lượng động thực vật nguy cấp, quý gán trực tiếp vào lớp liệu hệ sinh thái nên cần chồng xếp lớp lớp liệu khoảng cách mức độ nhiễm nước để tích hợp liệu Sử dụng công cụ Zonal Statistics, Spatial Join, Join Data để gán thuộc tính khoảng cách chất lượng nước vào hệ sinh thái đất ngập nước thuộc khu vực bảo tồn lựa chọn nghiên cứu; sử dụng công cụ xuất liệu thuộc tính ArcGIS Pro 3.0 định dạng Excel phục vụ tính tốn, xếp hạng mức độ ưu tiên bảo vệ, phục hồi hệ sinh thái đất ngập nước Số liệu đánh giá mức độ ưu tiên bảo vệ, phục hồi hệ sinh thái đất ngập nước trình bày Bảng 3.7 3.8 Bảng 3.7 Số liệu phục vụ đánh giá mức độ ưu tiên bảo vệ HST ĐNN Đồng Tháp Mười Khu vực TC LS DTME XQ GG DTMM TL DSGT BV1 667 382 337 334 389 91 0 BV2 57 20 15 10 0 BV3 12 0 0 BV4 BV5 2280,9 2036,2 1942,2 3037,4 2870,9 2361,9 2439,7 3380,1 1066,7 1149,6 1975,6 2146,9 1337,9 1094,2 1174,8 1821,1 Chỉ tiêu BV6 BV7 BV8 BV9 937,9 774,5 655,7 313,1 1368,7 810,6 694,7 1812,4 0 0,25 0 0 77.0 21.0 1.0 5.0 19.0 2.0 10.0 210,9 127,7 163,9 84,6 124,1 909,4 620,5 186,7 BV10 66.000 30.200 72.000 169.582 90.500 75.600 BV11 6,7 2,833 16,92 34,5 121,5 6,3 Bảng 3.8 Số liệu đánh giá mức độ ưu tiên phục hồi HST ĐNN Đồng Tháp Mười Khu vực TC LS DTME XQ GG DTMM TL DSGT PH1 -383,8 0,0 0,0 -24,3 0,0 0,0 0,0 0,0 Chỉ tiêu PH2 1041 258 0 0 0 PH3 5 4 Lớp liệu đồ đánh giá HST ĐNN chứa thông tin trường liệu tương ứng với tiêu chí đánh giá mục tiêu ưu tiên bảo vệ, phục hồi Kết thể hình 108 Hình 3.22 Dữ liệu tích hợp tiêu chí đánh giá mức độ ưu tiên bảo vệ, phục hồi HST ĐNN Đồng Tháp Mười 3.5 Phân vùng ưu tiên bảo vệ, phục hồi hệ sinh thái đất ngập nước khu vực Đồng Tháp Mười 3.5.1 Phân vùng ưu tiên bảo vệ hệ sinh thái đất ngập nước 3.5.1.1 Tính tốn trọng số tiêu chí đánh giá AHP Mức độ quan trọng theo cặp 11 tiêu xác định sở xin ý kiến chuyên gia liên quan đến lĩnh vực bảo tồn thiên nhiên đa dạng sinh học (03) tài nguyên môi trường (02), quản lý đất đai, kinh tế, xã hội (02) Việc xin ý kiến thực trực tuyến thông qua trang web https://bpmsg.com/ Đây ứng dụng trực tuyến miễn phí cho phép chuyên gia xác định mức độ quan trọng tương đối cặp tiêu tự động tính tốn xác định tính qn việc cho điểm Các chuyên gia thực điều chỉnh thang cho điểm tiêu chí cho tỷ số quán (CR) nhỏ 10%, lúc thang cho điểm thoả mãn Sau tổng hợp ý kiến chuyên gia, luận án kiểm tra lại mức độ quán địa trên, chi tiết thể Phụ lục Trong đó, Số lồi nguy cấp, q (BV2) có mức độ ưu tiên 22,8%, Tổng số loài (BV1) 18,7%, Số loài ngoại lai xâm hại (BV3) 13,4% Ma trận so sánh cặp ma trận chuẩn tiêu ưu tiên bảo vệ thể Bảng 3.9 Bảng 3.10 109 Bảng 3.9 Ma trận so sánh cặp tiêu ưu tiên bảo vệ HST ĐNN Đồng Tháp Mười Tên tiêu Tổng số loài Số loài nguy cấp, quý Số lượng loài ngoại lai xâm hại Khoảng cách từ khu vực chăn nuôi đến khu vực ĐNN Khoảng cách từ khu vực công nghiệp đến khu vực ĐNN Khoảng cách từ khu vực dân cư đến khu vực ĐNN Khoảng cách từ đường giao thông đến khu vực ĐNN Độ mặn nước Nhạy cảm cháy rừng Số lượng khách du lịch trung bình năm Doanh thu từ du lịch trung bình năm Tổng Chỉ tiêu BV1 BV1 BV2 BV3 BV4 BV5 BV6 BV7 BV8 BV9 BV10 BV11 1/2 3 BV2 3 3 BV3 1/2 1/3 2 BV4 1/3 1/3 1/2 2 BV5 1/3 1/3 1/2 1/2 2 2 BV6 1/3 1/3 1/2 1/2 1/2 3 BV7 1/4 1/3 1/3 1/3 1/2 1/2 2 BV8 1/5 1/4 1/5 1/4 1/2 1/3 1/2 5 BV9 1/7 1/9 1/6 1/5 1/2 1/3 1/2 1/5 BV10 1/4 1/5 1/4 1/3 1/3 1/3 1/3 1/2 1/2 BV11 1/9 1/8 1/7 1/6 1/6 1/6 1/6 1/5 1/3 1/7 5,45 3,85 8,59 11,28 13,50 14,67 19,50 26,90 40,83 30,14 64,00 Bảng 3.10 Ma trận chuẩn so sánh cặp tiêu ưu tiên bảo vệ HST ĐNN Đồng Tháp Mười Tên tiêu Tổng số loài Số loài nguy cấp, quý Số lượng loài ngoại lai xâm hại Khoảng cách từ khu vực chăn nuôi đến khu vực ĐNN Khoảng cách từ khu vực công nghiệp đến khu vực ĐNN Khoảng cách từ khu vực dân cư đến khu vực ĐNN Khoảng cách từ đường giao thông đến khu Chỉ tiêu BV1 BV1 BV2 BV3 BV4 BV5 BV6 BV7 BV8 BV9 BV10 BV11 0,18 0,13 0,23 0,27 0,22 0,20 0,21 0,19 0,17 0,13 0,14 Trọng số 0,19 BV2 0,37 0,26 0,35 0,27 0,22 0,20 0,15 0,15 0,22 0,17 0,13 0,23 BV3 0,09 0,09 0,12 0,18 0,15 0,14 0,15 0,19 0,15 0,13 0,11 0,14 BV4 0,06 0,09 0,06 0,09 0,15 0,14 0,15 0,15 0,12 0,10 0,09 0,11 BV5 0,06 0,09 0,06 0,04 0,07 0,14 0,10 0,07 0,05 0,10 0,09 0,08 BV6 0,06 0,09 0,06 0,04 0,04 0,07 0,10 0,11 0,07 0,10 0,09 0,08 BV7 0,05 0,09 0,04 0,03 0,04 0,03 0,05 0,07 0,05 0,10 0,09 0,06 110 vực ĐNN Độ mặn nước BV8 0,04 0,06 0,02 0,02 0,04 0,02 0,03 0,04 0,12 0,07 0,08 0,05 Nhạy cảm cháy rừng Số lượng khách du lịch trung bình năm Doanh thu từ du lịch trung bình năm BV9 0,03 0,03 0,02 0,02 0,04 0,02 0,03 0,01 0,02 0,07 0,05 0,03 BV10 0,05 0,05 0,03 0,03 0,02 0,02 0,02 0,02 0,01 0,03 0,11 0,04 BV11 0,02 0,03 0,02 0,01 0,01 0,01 0,01 0,01 0,01 0,00 0,02 0,01 Với 11 tiêu, tính tốn giá trị riêng ma trận (λmax) = 11,93, số quán (CI = 0,09), tra bảng số ngẫu nhiên (RI = 1,51), tính tốn tỷ số qn CR = 0,06 < 0,1 cho thấy tính quán cho điểm so sánh cặp trọng số đáp ứng yêu cầu chấp nhận sử dụng tính tốn bước 3.5.1.2 Xác định thứ bậc ưu tiên bảo vệ TOPSIS - Xây dựng ma trận định ma trận định chuẩn hoá Ma trận định thể giá trị tiêu sử dụng để đánh giá mức độ ưu tiên Đây ma trận với đại lượng có thứ nguyên, để so sánh, đánh giá cần đưa giá trị đại lượng khơng có thứ ngun cách tính tốn ma trận định chuẩn hố Kết thể Bảng 3.11 3.12 Bảng 3.11 Ma trận định ưu tiên bảo vệ HST ĐNN Đồng Tháp Mười CT KV TC LS DTME XQ GG DTMM TL DSGT BV1 BV2 BV3 BV4 BV5 BV6 BV7 BV8 BV9 BV10 BV11 667 382 337 334 389 91 0 987,7 57 20 15 10 0 63,9 12 0 0 13,7 2280,9 2036,2 1942,2 3037,4 2870,9 2361,9 2439,7 3380,1 7.317,4 1066,7 1149,6 1975,6 2146,9 1337,9 1094,2 1174,8 1821,1 4.319,2 937,9 774,5 655,7 313,1 1368,7 810,6 694,7 1812,4 2.881,9 210,9 127,7 163,9 84,6 124,1 909,4 620,5 186,7 1.164,9 0 0,25 0 0 0,3 77,0 21,0 1,0 5,0 19,0 2,0 10,0 82,8 66000 30200 72000 169582 90500 75600 230.468,4 6,7 2,833 16,92 34,5 121,5 6,3 127,8 Bảng 3.12 Ma trận định chuẩn hóa ưu tiên bảo vệ HST ĐNN Đồng Tháp Mười CT KV TC LS DTME XQ GG DTMM TL DSGT BV1 BV2 BV3 BV4 BV5 BV6 BV7 BV8 BV9 BV10 BV11 0,68 0,39 0,34 0,34 0,39 0,09 - 0,89 0,31 0,13 0,23 0,11 0,16 - 0,87 0,22 0,44 - 0,7 0,7 0,7 0,6 0,6 0,7 0,7 0,5 0,8 0,7 0,5 0,5 0,7 0,7 0,7 0,6 0,7 0,7 0,8 0,9 0,5 0,7 0,8 0,4 0,8 0,9 0,9 0,9 0,9 0,2 0,5 0,8 1,0 1,0 1,0 1,0 1,0 1,0 1,0 0,93 0,25 0,01 0,06 0,23 0,02 0,12 0,29 0,13 0,31 0,74 0,39 0,33 0,05 0,02 0,13 0,27 0,95 0,05 - Tính tốn ma trận định chuẩn hố trọng số Đây bước tính tốn lại ma trận định dựa ma trận định 111 chuẩn hố (quy đại lượng khơng có thứ ngun) ma trận trọng số tiêu chí tính tốn AHP Kết thể Bảng 3.13 Bảng 3.13 Ma trận định chuẩn hóa trọng số ưu tiên bảo vệ HST ĐNN Đồng Tháp Mười W 0,19 0,23 0,14 0,11 0,08 0,08 0,06 0,05 0,03 0,04 0,01 BV1 BV2 BV3 BV4 BV5 BV6 BV7 BV8 BV9 BV10 BV11 TC 0,1273 0,2012 0,1178 0,0749 0,0602 0,0513 0,0476 0,0488 0,0273 0,0103 0,0007 LS 0,0729 0,0706 0,0295 0,0786 0,0587 0,0556 0,0518 0,0488 0,0074 0,0047 0,0003 DTME 0,0643 0,0282 0,0000 0,0800 0,0434 0,0587 0,0500 0,0000 0,0000 0,0000 0,0000 XQ 0,0638 0,0529 0,0589 0,0637 0,0402 0,0678 0,0539 0,0488 0,0004 0,0112 0,0018 GG 0,0743 0,0247 0,0000 0,0661 0,0552 0,0399 0,0520 0,0488 0,0018 0,0264 0,0037 DTMM 0,0174 0,0353 0,0000 0,0737 0,0597 0,0546 0,0128 0,0488 0,0067 0,0000 0,0000 TL 0,0000 0,0000 0,0000 0,0726 0,0582 0,0577 0,0272 0,0488 0,0007 0,0141 0,0132 DSGT 0,0000 0,0000 0,0000 0,0586 0,0463 0,0282 0,0488 0,0488 0,0035 0,0118 0,0007 CT KV - Xác định lời giải lý tưởng dương, lý tưởng âm Lời giải lý tưởng dương (V+j), lý tưởng âm (V-j) xác định theo công thức 2.5, 2.6 cho tiêu chí lợi ích tiêu chí tiêu cực BV1, BV2, BV3, BV10, BV9, BV11 tiêu đại diện cho tính lợi ích ưu tiên bảo vệ (giá trị cao mức độ ưu tiên lớn); BV4, BV5, BV6, BV7, BV8 tiêu mang tính tiêu cực theo ưu tiên bảo vệ (giá trị thấp mức độ ưu tiên cao ngược lại giá trị cao mức độ ưu tiên thấp) Kết thể Bảng 3.14 Bảng 3.14 Xác định lời giải lý tưởng dương, lý tưởng âm tiêu chí đánh giá ưu tiên bảo vệ HST ĐNN Đồng Tháp Mười CT LG V+j= V-j= BV1 BV2 BV3 BV4 BV5 BV6 BV7 BV8 BV9 BV10 BV11 0,1273 0,0000 0,2012 0,0000 0,1178 0,0000 0,0586 0,0800 0,0402 0,0602 0,0282 0,0678 0,0128 0,0539 0,0000 0,0488 0,0273 0,0000 0,0264 0,0000 0,0132 0,0000 - Xác định khoảng cách tổng hợp từ khu vực vực đất ngập nước lựa chọn tới lời giải lý tưởng dương, lý tưởng âm Việc tính khoảng cách tổng hợp theo tiêu chí khu vực đất ngập nước đến lời giải lý tưởng dương (S+i), lý tưởng âm (S-i) tương ứng (theo công thức 2.7 2.8) nhằm mục đích xác định mức độ gần gũi lựa chọn tới lời giải lý tưởng dương Kết tính toán thể Bảng 3.15 112