ĐẠI HỌC HUẾ TRƯỜNG ĐẠI HỌC SƯ PHẠM CAO THỊ THANH NGUYÊN NGHIÊN CỨU MẬT ĐỘ, BIẾN ĐỘNG SỐ LƯỢNG VÀ XÁC SUẤT PHÁT HIỆN LỒI NHƠNG CÁT SỌC (Leiolepis guentherpetersi) Ở VÙNG CÁT VEN BIỂN HUYỆN PHÚ LỘC, TỈNH THỪA THIÊN HUẾ LUẬN VĂN THẠC SĨ SINH HỌC THEO ĐỊNH HƯỚNG NGHIÊN CỨU Thừa Thiên Huế, năm 2018 ĐẠI HỌC HUẾ TRƯỜNG ĐẠI HỌC SƯ PHẠM CAO THỊ THANH NGUYÊN NGHIÊN CỨU MẬT ĐỘ, BIẾN ĐỘNG SỐ LƯỢNG VÀ XÁC SUẤT PHÁT HIỆN LỒI NHƠNG CÁT SỌC (Leiolepis guentherpetersi) Ở VÙNG CÁT VEN BIỂN HUYỆN PHÚ LỘC, TỈNH THỪA THIÊN HUẾ Chuyên ngành: Động vật học Mã số: 60420103 Giảng viên hướng dẫn khoa học: TS NGƠ VĂN BÌNH GS.TS NGÔ ĐẮC CHỨNG LUẬN VĂN THẠC SĨ SINH HỌC THEO ĐỊNH HƯỚNG NGHIÊN CỨU Thừa Thiên Huế, năm 2018 LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu riêng tơi hướng dẫn khoa học TS Ngơ Văn Bình GS.TS Ngơ Đắc Chứng Các số liệu kết nghiên cứu ghi luận văn trung thực, đồng tác giả cho phép sử dụng chưa cơng bố cơng trình khác Tác giả luận văn Cao Thị Thanh Nguyên LỜI CẢM ƠN Đầu tiên, tơi xin bày tỏ lịng biết ơn sâu sắc đến thầy giáo GS.TS Ngô Đắc Chứng thầy giáo TS Ngơ Văn Bình tận tình hướng dẫn, trang bị cho kiến thức, kỹ cần thiết giúp đỡ suốt thời gian thực đề tài Tôi xin chân thành cảm ơn Ban Chủ nhiệm khoa Sinh học; thầy giáo, cô giáo Bộ môn Động vật học quý thầy cô khoa Sinh học phòng Đào tạo Sau Đại học Trường Đại học Sư phạm - Đại học Huế Xin chân thành cảm ơn lãnh đạo Sở Giáo dục Đào tạo tỉnh Kon Tum; Ban Giám hiệu quý thầy tổ Hóa - Sinh trường PT DTNT Tu Mơ Rông tỉnh Kon Tum tạo điều kiện thuận lợi cho tơi q trình học tập thực đề tài Trong trình thực luận văn, nhận giúp đỡ bà nhân dân sống huyện Phú Lộc - nơi thực đề tài mặt thông tin điều kiện để thực điều tra, khảo sát; bạn lớp Động vật học Khóa 25 hỗ trợ, giúp đỡ tơi q trình thu mẫu Tơi xin trân trọng cảm ơn giúp đỡ quý báu Cuối cùng, xin gửi lời cảm ơn đến bạn bè, đồng nghiệp đặc biệt gia đình ln quan tâm, chia sẻ khó khăn, động viên tạo điều kiện cho yên tâm học tập nghiên cứu để hoàn thành luận văn Huế, tháng 10 năm 2018 Tác giả luận văn Cao Thị Thanh Nguyên MỤC LỤC Trang MỞ ĐẦU Lí chọn đề tài Mục đích nghiên cứu Nội dung nghiên cứu Ý nghĩa đề tài CHƯƠNG TỔNG QUAN TÀI LIỆU 1.1 Tổng quan vấn đề nghiên cứu 1.1.1 Những vấn đề nghiên cứu giống Leiolepis giới 1.1.2 Tình hình nghiên cứu giống Leiolepis Việt Nam 13 1.1.3 Những nghiên cứu liên quan đến tỷ suất chiếm điểm chương trình giám sát động vật quy mơ lớn 17 1.2 Điều kiện tự nhiên vùng nghiên cứu 20 1.2.1 Giới hạn khu vực nghiên cứu 20 1.2.2 Đặc điểm khu vực nghiên cứu 21 CHƯƠNG VẬT LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU 25 2.1 Vật liệu nghiên cứu 25 2.2 Phương pháp nghiên cứu 25 2.2.1 Thời gian địa điểm nghiên cứu 25 2.2.2 Đối tượng nghiên cứu 27 2.2.3 Phương pháp khảo sát thực địa 27 2.2.4 Phương pháp xác định đặc điểm hình thái lồi 28 2.2.5 Phương pháp nghiên cứu việc sử dụng vi môi trường sống 29 2.2.6 Phương pháp nghiên cứu phương thức hoạt động 29 2.2.7 Phương pháp ước tính mật độ quần thể 29 2.2.8 Phương pháp nghiên cứu xác suất phát tỷ suất chiếm điểm 30 2.2.9 Phương pháp điều tra, vấn 34 2.2.10 Xử lý số liệu phân tích thống kê .34 CHƯƠNG KẾT QUẢ VÀ THẢO LUẬN 35 3.1 Đặc điểm hình thái lồi Nhơng cát sọc Leiolepis guentherpetersi 35 3.1.1 Đặc điểm hình thái 35 3.1.2 Mối quan hệ số số đo hình thái 37 3.2 Vi môi trường sống phương thức hoạt động Nhông cát sọc .40 3.2.1 Sử dụng vi môi trường sống 40 3.2.2 Phương thức hoạt động 43 3.3 Ước tính mật độ quần thể, biến động số lượng cá thể Nhơng cát sọc 45 3.3.1 Ước tính mật độ quần thể 45 3.3.1.1 Ước tính mật độ quần thể theo ô tiêu chuẩn 45 3.3.1.2 Ước tính mật độ hang 47 3.3.1 Ước tính mật độ quần thể 458 3.4 Xác suất phát tỷ suất chiếm điểm 50 3.4.1 Xác suất phát tỷ suất chiếm điểm lồi hai mơi trường sống mùa mưa 50 3.4.2 Xác suất phát tỷ suất chiếm điểm hai môi trường sống mùa khô 55 3.5 Sử dụng bền vững Nhông cát sọc 59 KẾT LUẬN VÀ KIẾN NGHỊ 62 Kết luận .62 Kiến nghị 63 TÀI LIỆU THAM KHẢO 64 PHỤ LỤC DANH MỤC CHỮ VIẾT TẮT BM : Khối lượng thể SVL : Chiều dài thân TL : Chiều dài đuôi HL : Dài đầu HW : Rộng đầu HH : Cao đầu HE : Chiều cao lỗ tai AG : Khoảng cách nách - háng DN : Khoảng cách hai mũi FL : Dài chi trước HB : Dài chi sau TIB : Dài cẳng chân FA : Dài cẳng tay DE : Đường kính lỗ tai ngồi DiE : Đường kính mắt MW : Rộng miệng AICc : Tiêu chuẩn thông tin Akaike kích thước mẫu nhỏ ∆AICc AIC weight ML K : Sự khác biệt tuyệt đối giá trị AIC so với mơ hình tốp : Tầm ảnh hưởng mơ hình : Model Likelihood = Khả xảy mơ hình : Thơng số mơ hình suy luận GB : Các điểm nghiên cứu thuộc môi trường sống gần biển XB : Các điểm nghiên cứu thuộc môi trường sống xa biển N : Nhiệt độ DA : Độ ẩm M : Tình hình nắng, mưa nắng mưa không xác định M1 : Nắng M2 : Mưa M3 : Tình hình nắng mưa khơng xác định DANH MỤC BẢNG Bảng 2.1 Các số đo hình thái Nhơng cát sọc 28 Bảng 2.2 Dự đoán chi tiết đánh giá tỷ suất chiếm điểm lồi 32 Bảng 3.1 Đặc điểm hình thái lồi Nhông cát sọc 36 Bảng 3.2 Nhiệt độ độ ẩm vi môi trường sống nơi phát lồi Nhơng cát sọc Leiolepis guentherpetersi vùng nghiên cứu 41 Bảng 3.3 Ước tính mật độ quần thể trung bình qua 18 khảo sát Lộc Vĩnh Vinh Hiền 46 Bảng 3.4 Ước tính mật độ hang quần thể Nhơng cát sọc dựa vào số lượng hang ô tiêu chuẩn .47 Bảng 3.5 Tóm tắt hai mơ hình để kiểm tra mức ý nghĩa thống kê khả phát lồi Nhơng cát sọc mùa mưa vùng nghiên cứu 51 Bảng 3.6 Tóm tắt q trình chọn lọc mơ hình AIC lồi Nhơng cát sọc thơng qua 50 điểm với khảo sát thuộc mùa mưa vùng nghiên cứu 53 Bảng 3.7 Tóm tắt mơ hình ứng cử viên để suy luận mức độ ảnh hưởng yếu tố thời tiết đến khả phát loài L guentherpetersi vùng nghiên cứu 54 Bảng 3.8 Tóm tắt q trình chon lọc mơ hình AIC lồi Nhơng cát sọc thơng qua 50 điểm với khảo sát thuộc mùa khô vùng nghiên cứu 57 Bảng 3.9 Tóm tắt mơ hình ứng cử viên để suy luận mức độ ảnh hưởng yếu tố nhiệt độ thời tiết đến khả phát loài L guentherpetersi vùng nghiên cứu 58 DANH MỤC HÌNH Hình 2.1 Bản đồ khu vực nghiên cứu, thu mẫu, giám sát .26 Hình 2.2 Nhơng cát sọc (Leiolepis guentherpetersi) vùng cát ven biển huyện Phú Lộc, tỉnh Thừa Thiên Huế .27 Hình 2.3 Mơ tả cách quan sát giám sát điểm 31 Hình 3.1 Nhơng cát sọc (Leiolepis guentherpetersi) vùng cát ven biển huyện Phú Lộc, tỉnh Thừa Thiên Huế .35 Hình 3.2 Hồi quy khơng tuyến tính chiều dài thân với khối lượng thể .37 Hình 3.3 Hồi quy khơng tuyến tính chiều dài với khối lượng thể .38 Hình 3.4 Hồi quy tuyến tính chiều dài thân với dài đầu 38 Hình 3.5 Hồi quy tuyến tính chiều dài thân với rộng đầu .39 Hình 3.6 Hồi quy tuyến tính chiều dài thân với khoảng cách nách - háng .40 Hình 3.7 Tỷ lệ phần trăm Nhơng cát sọc vi môi trường sống nhiệt độ tương ứng vi mơi trường sống 42 Hình 3.8 Tỷ lệ phần trăm Nhông cát sọc vi môi trường sống độ ẩm tương ứng vi mơi trường sống 42 Hình 3.9 Mối quan hệ nhiệt độ khơng khí với số lượng cá thể hoạt động .43 Hình 3.10 Thời gian hoạt động Nhơng cát sọc 44 Hình 3.11 Phương thức hoạt động L guentherpetersi .44 Hình 3.12 Mật độ quần thể Nhơng cát sọc vùng Lộc Vĩnh vùng Vinh Hiền hai mùa nghiên cứu 46 Hình 3.13 Số lượng cá thể Nhông cát sọc L guentherpetersi theo thời gian ngày môi trường sống 48 Hình 3.14 Số lượng cá thể Nhông cát sọc theo thời gian ngày mùa mưa mùa khô .49 Hình 3.15 Xác suất phát lồi Nhơng cát sọc (Leiolepis guentherpetersi) từ mơ hình [ψ(.),p(survey)] vùng nghiên cứu mùa mưa 52 Hình 3.16 Xác suất phát lồi Nhơng cát sọc (Leiolepis guentherpetersi) từ mơ hình [ψ(.),p(survey)] vùng nghiên cứu mùa khô 56 Hình 3.17 Hoạt động đốt khai thác gỗ phi lao vùng cát ven biển huyện Phú Lộc (A, B); đào hang bắt Nhông cát sọc (C); xây dựng khu nghỉ dưỡng sinh thái ven biển huyện Phú Lộc (D) 60 MỞ ĐẦU Lí chọn đề tài Nhơng cát sọc (Leiolepis guentherpetersi Darevsky & Kupriyanova, 1993) loài đặc hữu Việt Nam, có vùng phân bố hẹp từ tỉnh Thừa Thiên Huế đến tỉnh Quảng Ngãi (Nguyen et al., 2009; Grismer et al., 2014) Chúng chủ yếu phân bố dải cát ven biển Lồi có phương thức sinh sản trinh sinh (parthenogenetic species), tất cá thể quần thể (Darevsky & Kupriyanova, 1993; Grismer & Grismer, 2010) Nhơng cát sọc có giá trị thực phẩm dược liệu Chúng có kích thước tương đối lớn, thịt thơm ngon bổ dưỡng Đối với trẻ nhỏ, người lớn tuổi, người đau ốm dùng thịt Nhông cát sọc để nâng cao thể lực Đồng thời, thịt Nhông cát sọc chế biến thành ăn đặc sản cho người vùng du khách số tỉnh miền Trung Người ta cịn dùng rượu ngâm Nhơng cát sọc để hỗ trợ chữa bệnh tê liệt, đau nhức xương khớp bắp, hen suyễn (Nguyễn Lân Hùng, 2009) Mặt khác, nhơng cát sọc ăn nhiều loại trùng có hại cho hệ sinh thái nông nghiệp cào cào, châu chấu, bọ xít, ruồi, bướm,… Do chúng mắt xích quan trọng chuỗi thức ăn đóng vai trị định việc cân hệ sinh thái vùng cát ven biển (Trần Quốc Dung Ngơ Quốc Trí, 2012) Hiện nay, việc mở rộng khu du lịch ven biển, xây dựng khu nghỉ dưỡng sinh thái, mở rộng vùng đất canh tác, nuôi trồng thủy sản,…đang diễn quy mô lớn vùng cát ven biển huyện Phú Lộc, tỉnh Thừa Thiên Huế làm môi trường sống tự nhiên lồi nhơng cát Hơn nữa, việc người thường xuyên săn bắt khiến loài bị giảm mạnh số lượng cá thể quần thể, làm giảm đa dạng sinh học, ảnh hưởng đến cân sinh học hệ sinh thái Các nghiên cứu Nhông cát sọc (L guentherpetersi) giới nước chủ yếu tập trung vào phân loại, ghi nhận phân bố, mơ tả đặc điểm hình thái, sinh thái điều kiện ni phân tích kiểu nhân (Darevsky & Kupriyanova, 1993; Lê Thị Nga Ngô Đắc Chứng, 2009; Nguyen et al., 2009; Trần Quốc Dung cộng sự, 2009; Grismer & Grismer, 2010; Trần Quốc Dung psi site6 : 0.1326 0.3146 0.0007 - 0.9702 psi site7 : 0.1326 0.3146 0.0007 - 0.9702 psi site8 : 0.1326 0.3146 0.0007 - 0.9702 psi site9 : 0.1326 0.3146 0.0007 - 0.9702 psi 10 site10 : 0.1326 0.3146 0.0007 - 0.9702 psi 11 site11 : 0.1326 0.3146 0.0007 - 0.9702 psi 12 site12 : 0.1326 0.3146 0.0007 - 0.9702 psi 13 site13 : 0.1326 0.3146 0.0007 - 0.9702 psi 14 site14 : 0.1326 0.3146 0.0007 - 0.9702 psi 15 site15 : 0.1326 0.3146 0.0007 - 0.9702 psi 16 site16 : 0.6383 1.0093 0.0003 - 0.9999 psi 17 site17 : 0.6383 1.0093 0.0003 - 0.9999 psi 18 site18 : 0.6383 1.0093 0.0003 - 0.9999 psi 19 site19 : 0.6383 1.0093 0.0003 - 0.9999 psi 20 site20 : 0.6383 1.0093 0.0003 - 0.9999 psi 21 site21 : 0.6383 1.0093 0.0003 - 0.9999 psi 22 site22 : 0.6383 1.0093 0.0003 - 0.9999 psi 23 site23 : 0.6383 1.0093 0.0003 - 0.9999 psi 24 site24 : 0.6383 1.0093 0.0003 - 0.9999 psi 25 site25 : 0.6383 1.0093 0.0003 - 0.9999 psi 26 site26 : 0.6383 1.0093 0.0003 - 0.9999 psi 27 site27 : 0.6383 1.0093 0.0003 - 0.9999 psi 28 site28 : 0.6383 1.0093 0.0003 - 0.9999 psi 29 site29 : 0.6383 1.0093 0.0003 - 0.9999 psi 30 site30 : 0.6383 1.0093 0.0003 - 0.9999 psi 31 site31 : 0.6383 1.0093 0.0003 - 0.9999 psi 32 site32 : 0.6383 1.0093 0.0003 - 0.9999 psi 33 site33 : 0.6383 1.0093 0.0003 - 0.9999 psi 34 site34 : 0.6383 1.0093 0.0003 - 0.9999 psi 35 site35 : 0.6383 1.0093 0.0003 - 0.9999 psi 36 site36 : 0.6383 1.0093 0.0003 - 0.9999 psi 37 site37 : 0.6383 1.0093 0.0003 - 0.9999 psi 38 site38 : 0.6383 1.0093 0.0003 - 0.9999 psi 39 site39 : 0.6383 1.0093 0.0003 - 0.9999 psi 40 site40 : 0.6383 1.0093 0.0003 - 0.9999 psi 41 site41 : 0.1326 0.3146 0.0007 - 0.9702 psi 42 site42 : 0.1326 0.3146 0.0007 - 0.9702 psi 43 site43 : 0.1326 0.3146 0.0007 - 0.9702 psi 44 site44 : 0.1326 0.3146 0.0007 - 0.9702 psi 45 site45 : 0.1326 0.3146 0.0007 - 0.9702 psi 46 site46 : 0.1326 0.3146 0.0007 - 0.9702 psi 47 site47 : 0.1326 0.3146 0.0007 - 0.9702 psi 48 site48 : 0.1326 0.3146 0.0007 - 0.9702 psi 49 site49 : 0.1326 0.3146 0.0007 - 0.9702 psi 50 site50 : 0.1326 0.3146 0.0007 - 0.9702 =========================================================== DERIVED parameter - Psi-conditional = [Pr(occ | detection history)] Site psi-cond Std.err 95% conf interval psi-cond site1 : 0.2891 0.1551 0.0848 - 0.6409 P 13 psi-cond site2 psi-cond site3 psi-cond site4 psi-cond site5 psi-cond site6 psi-cond site7 psi-cond site8 psi-cond site9 psi-cond 10 site10 psi-cond 11 site11 psi-cond 12 site12 psi-cond 13 site13 psi-cond 14 site14 psi-cond 15 site15 psi-cond 16 site16 psi-cond 17 site17 psi-cond 18 site18 psi-cond 19 site19 psi-cond 20 site20 psi-cond 21 site21 psi-cond 22 site22 psi-cond 23 site23 psi-cond 24 site24 psi-cond 25 site25 psi-cond 26 site26 psi-cond 27 site27 psi-cond 28 site28 psi-cond 29 site29 psi-cond 30 site30 psi-cond 31 site31 psi-cond 32 site32 psi-cond 33 site33 psi-cond 34 site34 psi-cond 35 site35 psi-cond 36 site36 psi-cond 37 site37 psi-cond 38 site38 psi-cond 39 site39 psi-cond 40 site40 psi-cond 41 site41 psi-cond 42 site42 psi-cond 43 site43 psi-cond 44 site44 psi-cond 45 site45 psi-cond 46 site46 psi-cond 47 site47 psi-cond 48 site48 psi-cond 49 site49 psi-cond 50 site50 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 0.2886 0.2941 0.2905 0.2976 0.2821 1.0000 0.2872 0.2808 0.2865 1.0000 0.2826 0.2862 1.0000 0.2193 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.1379 1.0000 0.1510 0.2140 0.1865 0.2122 0.1808 0.1775 0.1626 0.1549 0.1562 0.1554 0.1570 0.1533 0.0000 0.1548 0.1532 0.1545 0.0000 0.1538 0.1546 0.0000 0.1348 0.0001 0.0001 0.0000 0.0000 0.0000 0.0001 0.0000 0.0001 0.0001 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0000 0.0001 0.0000 0.0001 0.0000 0.0000 0.0001 0.0000 0.0000 0.0000 0.1016 0.0000 0.1076 0.1344 0.1235 0.1338 0.1216 0.1202 0.1135 0.0846 - 0.6403 0.0871 - 0.6455 0.0855 - 0.6421 0.0886 - 0.6486 0.0818 - 0.6341 1.0000 - 1.0000 0.0839 - 0.6394 0.0811 - 0.6333 0.0837 - 0.6386 1.0000 - 1.0000 0.0818 - 0.6354 0.0834 - 0.6386 1.0000 - 1.0000 0.0566 - 0.5679 0.0000 - 1.0000 0.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 0.0000 - 1.0000 0.0000 - 1.0000 0.0000 - 1.0000 0.0000 - 1.0000 0.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 0.0000 - 1.0000 1.0000 - 1.0000 0.0000 - 1.0000 0.0000 - 1.0000 0.0000 - 1.0000 0.0000 - 1.0000 1.0000 - 1.0000 0.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 0.0291 - 0.4605 1.0000 - 1.0000 0.0332 - 0.4796 0.0538 - 0.5658 0.0444 - 0.5307 0.0531 - 0.5639 0.0423 - 0.5245 0.0412 - 0.5202 0.0365 - 0.4987 P 14 =========================================================== Assessing Model Fit for Single-season model: History(cohort) Observed Expected Chi-square 000000000( 0) 33.0000 31.397298074 0.08 000000001( 6) 6.0000 2.174268369 6.73 000000010( 10) 2.0000 2.309931757 0.04 001000011( 17) 1.0000 0.038077972 24.30 100000001( 18) 1.0000 1.402289526 0.12 100100000( 26) 1.0000 0.346489765 1.23 000100010( 27) 1.0000 0.103634892 7.75 001100000( 28) 1.0000 0.079373391 10.68 100000000( 29) 1.0000 3.610504111 1.89 010000010( 36) 1.0000 0.236996288 2.46 010010001( 38) 1.0000 0.013195308 73.80 010000000( 39) 1.0000 1.552568753 0.20 Test Statistic = 136.01 min(expect)=1.319531e-002 -Test Statistic (data) = 136.01 From 100 parametric bootstraps Probability of test statistic >= observed = 0.1485 Lowest simulated Test Stat = 13.3293 Average simulated Test Stat = 128.6596 Median simulated Test Stat = 47.3123 Highest simulated Test Stat = 2495.8243 Estimate of c-hat = 1.0531 (=TestStat/AvgTestStat) Distribution of simulated test statistics -TestStat Prop 187.1868 0.8400:******************************************************************* 311.9780 0.0900:******** 436.7693 0.0100:* 561.5605 0.0100:* 686.3517 0.0100:* 811.1429 0.0100:* 935.9341 0.0100:* 1060.7253 0.0100:* -avg param values for 100 bootstrap iterations -0.357021 26.298602 0.226425 -9.005072 -51.031347 -59.618697 0.708071 Bootstrap Variance-Covariance Matrix intercept Psi-Bar 13.0926 0.0001 0.0259 -1.3523 -0.0000 0.0000 0.2111 P 15 0.0001 0.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 0.0259 0.0000 0.0018 -0.0659 -0.0000 0.0000 0.0016 -1.3523 -0.0000 -0.0659 2.4821 -0.0000 0.0000 -0.0725 -0.0000 0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 0.0000 -0.0000 0.0000 0.0000 0.2111 0.0000 0.0016 -0.0725 -0.0000 0.0000 0.0107 Bootstrap estimate of SE for overall proportion of sites occupied = 0.1035 CPU time= 11 seconds (0.18 min) P 16 PHỤ LỤC 4b: KẾT QUẢ CỦA MƠ HÌNH CĨ NHIỀU THƠNG SỐ NHÂT [ψ(GB),p(N,DA,M)] PRESENCE - Presence/Absence-Site Occupancy data analysis Wed Aug 15 14:32:36 2018, Version 2.12.10 ==>C:\Progra~1\Presence\presence.exe ==>i=simdata.pao ==>l=pres_psi(GB)_p(N_DA_M1_M2_M3).out ==>name=psi(GB),p(N,DA,M1,M2,M3) ==>j=simdata.dm ==>bothvc ==>nose ==>nboot=100 ==>boot2=100 ==>link=cloglog ==>amoeba ==>echo varcov: nsig=6 eps=1.000000e-002 no model name N,T >50,9 1:000000000 2:000000000 3:000000000 4:000000000 5:000000000 6:000000000 7:000010000 8:000000000 9:000000000 10:000000000 11:001000000 12:000000000 13:000000000 14:000000000 15:000000000 16:000000000 17:000000000 18:100100000 19:001001000 20:001000010 21:000000000 22:000000000 23:000000000 24:000000000 25:000000000 26:010001100 27:000000100 P 17 28:010000000 29:000001101 30:000000011 31:001101000 32:001100000 33:001111101 34:000000000 35:000000000 36:000000000 37:100110001 38:100110011 39:001111110 40:000000010 41:000000000 42:000000000 43:000010110 44:000000000 45:000000000 46:000000000 47:000000000 48:000000000 49:000000000 50:000000000 ********* Input Data summary ******* Number of sites = 50 Number of sampling occasions = Number of states =2 Number of missing observations = Data checksum = 47221 NSiteCovs >2 site_covname[0]=GB site_covname[1]=XB NSampCovs >5 samp_covname[0]=N samp_covname[1]=DA samp_covname[2]=M1 samp_covname[3]=M2 samp_covname[4]=M3 Primary periods=1 Secondary periods: Naive occupancy estimate = 0.3600 MUA KHO-Presence 12.10 N=50 T=9 Groups=1 bootstraps=100 >1-9 Matrix 1: rows=2, cols=3 -,a1,a2, psi GB ======================== P 18 Matrix 2: rows=10, cols=6 -,b1 , b2, b3, b4, b5, P[1] N DA M1 M2 M3 P[2] N DA M1 M2 M3 P[3] N DA M1 M2 M3 P[4] N DA M1 M2 M3 P[5] N DA M1 M2 M3 P[6] N DA M1 M2 M3 P[7] N DA M1 M2 M3 P[8] N DA M1 M2 M3 P[9] N DA M1 M2 M3 ======================== Matrix 3: rows=0, cols=0 ======================== Matrix 4: rows=0, cols=0 ======================== Matrix 5: rows=0, cols=0 ======================= Matrix 6: rows=0, cols=0 ======================== modtype=1 Custom Model: psi(GB),p(N,DA,M1,M2,M3) attempt 0/0 (iv= 0.00 0.00 0.00 ), optmiz done -2loglike= 218.5246 random initial value vectors attempted: Initial values: 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 Number of significant digits = 7.3 Model has been fit using the logistic link for p, cloglog for psi Number of parameters =7 Number of function calls = 256 -2log(likelihood) = 218.5246 AIC = 232.5246 Untransformed Estimates of coefficients for covariates (Beta's) =========================================================== estimate std.error A1 psi.a1 : -1.650895 0.000000 A2 psi.GB : 1.960749 0.000000 B1 P[1].N : 0.312714 0.000000 B2 P[1].DA : 0.046383 0.000000 B3 P[1].M1 : -14.383210 0.000000 B4 P[1].M2 : 23.908783 0.000000 B5 P[1].M3 : -13.668716 0.000000 =========================================================== Individual Site estimates of Site estimate Std.err 95% conf interval psi site1 : 0.1746 0.0000 -4.2116 - -4.2116 psi site2 : 0.1746 0.0000 -4.2116 - -4.2116 psi site3 : 0.1746 0.0000 -4.2116 - -4.2116 psi site4 : 0.1746 0.0000 -4.2116 - -4.2116 P 19 psi site5 : 0.1746 0.0000 -4.2116 - -4.2116 psi site6 : 0.1746 0.0000 -4.2116 - -4.2116 psi site7 : 0.1746 0.0000 -4.2116 - -4.2116 psi site8 : 0.1746 0.0000 -4.2116 - -4.2116 psi site9 : 0.1746 0.0000 -4.2116 - -4.2116 psi 10 site10 : 0.1746 0.0000 -4.2116 - -4.2116 psi 11 site11 : 0.1746 0.0000 -4.2116 - -4.2116 psi 12 site12 : 0.1746 0.0000 -4.2116 - -4.2116 psi 13 site13 : 0.1746 0.0000 -4.2116 - -4.2116 psi 14 site14 : 0.1746 0.0000 -4.2116 - -4.2116 psi 15 site15 : 0.1746 0.0000 -4.2116 - -4.2116 psi 16 site16 : 0.7972 0.0000 0.2664 - 0.2664 psi 17 site17 : 0.7982 0.0000 0.2664 - 0.2664 psi 18 site18 : 0.7972 0.0000 0.2664 - 0.2664 psi 19 site19 : 0.7972 0.0000 0.2664 - 0.2664 psi 20 site20 : 0.7972 0.0000 0.2664 - 0.2664 psi 21 site21 : 0.7972 0.0000 0.2664 - 0.2664 psi 22 site22 : 0.7972 0.0000 0.2664 - 0.2664 psi 23 site23 : 0.7972 0.0000 0.2664 - 0.2664 psi 24 site24 : 0.7972 0.0000 0.2664 - 0.2664 psi 25 site25 : 0.7972 0.0000 0.2664 - 0.2664 psi 26 site26 : 0.7972 0.0000 0.2664 - 0.2664 psi 27 site27 : 0.7972 0.0000 0.2664 - 0.2664 psi 28 site28 : 0.7972 0.0000 0.2664 - 0.2664 psi 29 site29 : 0.7972 0.0000 0.2664 - 0.2664 psi 30 site30 : 0.7972 0.0000 0.2664 - 0.2664 psi 31 site31 : 0.7972 0.0000 0.2664 - 0.2664 psi 32 site32 : 0.7972 0.0000 0.2664 - 0.2664 psi 33 site33 : 0.7972 0.0000 0.2664 - 0.2664 psi 34 site34 : 0.7972 0.0000 0.2664 - 0.2664 psi 35 site35 : 0.7972 0.0000 0.2664 - 0.2664 psi 36 site36 : 0.7972 0.0000 0.2664 - 0.2664 psi 37 site37 : 0.7972 0.0000 0.2664 - 0.2664 psi 38 site38 : 0.7972 0.0000 0.2664 - 0.2664 psi 39 site39 : 0.7972 0.0000 0.2664 - 0.2664 psi 40 site40 : 0.7972 0.0000 0.2664 - 0.2664 psi 41 site41 : 0.1746 0.0000 -4.2116 - -4.2116 psi 42 site42 : 0.1746 0.0000 -4.2116 - -4.2116 psi 43 site43 : 0.1746 0.0000 -4.2116 - -4.2116 psi 44 site44 : 0.1746 0.0000 -4.2116 - -4.2116 psi 45 site45 : 0.1746 0.0000 -4.2116 - -4.2116 psi 46 site46 : 0.1746 0.0000 -4.2116 - -4.2116 psi 47 site47 : 0.1746 0.0000 -4.2116 - -4.2116 psi 48 site48 : 0.1746 0.0000 -4.2116 - -4.2116 psi 49 site49 : 0.1746 0.0000 -4.2116 - -4.2116 psi 50 site50 : 0.1746 0.0000 -4.2116 - -4.2116 =========================================================== DERIVED parameter - Psi-conditional = [Pr(occ | detection history)] Site psi-cond Std.err 95% conf interval P 20 psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond psi-cond 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 site1 site2 site3 site4 site5 site6 site7 site8 site9 site10 site11 site12 site13 site14 site15 site16 site17 site18 site19 site20 site21 site22 site23 site24 site25 site26 site27 site28 site29 site30 site31 site32 site33 site34 site35 site36 site37 site38 site39 site40 site41 site42 site43 site44 site45 site46 site47 site48 site49 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 0.1017 0.1013 0.1066 0.1052 0.1100 0.1042 1.0000 0.1006 0.1048 0.0919 1.0000 0.0958 0.1041 0.1025 0.0983 0.6075 0.5021 1.0000 1.0000 1.0000 0.5605 0.5819 0.5818 0.5711 0.0471 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.0719 0.0390 0.0413 1.0000 1.0000 1.0000 1.0000 0.0093 0.0092 1.0000 0.0055 0.0048 0.0027 0.0035 0.0009 0.0009 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.1017 - 0.1017 0.1013 - 0.1013 0.1066 - 0.1066 0.1052 - 0.1052 0.1100 - 0.1100 0.1042 - 0.1042 1.0000 - 1.0000 0.1006 - 0.1006 0.1048 - 0.1048 0.0919 - 0.0919 1.0000 - 1.0000 0.0958 - 0.0958 0.1041 - 0.1041 0.1025 - 0.1025 0.0983 - 0.0983 0.6075 - 0.6075 0.5021 - 0.5021 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 0.5605 - 0.5605 0.5819 - 0.5819 0.5818 - 0.5818 0.5711 - 0.5711 0.0471 - 0.0471 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 0.0719 - 0.0719 0.0390 - 0.0390 0.0413 - 0.0413 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 1.0000 - 1.0000 0.0093 - 0.0093 0.0092 - 0.0092 1.0000 - 1.0000 0.0055 - 0.0055 0.0048 - 0.0048 0.0027 - 0.0027 0.0035 - 0.0035 0.0009 - 0.0009 0.0009 - 0.0009 P 21 psi-cond 50 site50 : 0.0012 0.0000 0.0012 - 0.0012 =========================================================== Assessing Model Fit for Single-season model: History(cohort) Observed Expected Chi-square 000000000( 0) 32.0000 31.388191614 0.01 000010000( 6) 1.0000 0.624555034 0.23 001000000( 10) 1.0000 0.452436778 0.66 100100000( 17) 1.0000 0.086485445 9.65 001001000( 18) 1.0000 0.120543260 6.42 001000010( 19) 1.0000 0.139522521 5.31 010001100( 25) 1.0000 0.015898634 60.91 000000100( 26) 1.0000 0.886595821 0.01 010000000( 27) 1.0000 0.230635305 2.57 000001101( 28) 1.0000 0.053453076 16.76 000000011( 29) 1.0000 0.157891793 4.49 001101000( 30) 1.0000 0.125211370 6.11 001100000( 31) 1.0000 0.178197021 3.79 001111101( 32) 1.0000 0.013508682 72.04 100110001( 36) 1.0000 0.021354479 44.85 100110011( 37) 1.0000 0.015776359 61.40 001111110( 38) 1.0000 0.020632203 46.49 000000010( 39) 1.0000 0.823822342 0.04 000010110( 42) 1.0000 0.122597206 6.28 Test Statistic = 362.54 min(expect)=1.350868e-002 -Test Statistic (data) = 362.54 From 100 parametric bootstraps Probability of test statistic >= observed = 0.6238 Lowest simulated Test Stat = 115.1269 Average simulated Test Stat = 480.9720 Median simulated Test Stat = 399.8550 Highest simulated Test Stat = 3571.1543 Estimate of c-hat = 0.7538 (=TestStat/AvgTestStat) Distribution of simulated test statistics -TestStat Prop 267.8366 0.0100:* 446.3943 0.3600:***************************** 624.9520 0.4000:********************************* 803.5097 0.1300:*********** 982.0674 0.0400:**** 1160.6252 0.0400:**** 1339.1829 0.0100:* -avg param values for 100 bootstrap iterations -2.182017 2.669925 0.345372 0.055376 P 22 -16.052454 26.820643 -15.353365 0.467266 Bootstrap Variance-Covariance Matrix intercept Psi-Bar 14.3896 -14.3430 -0.0195 -0.0158 1.6321 -3.1590 1.8352 0.0046 -14.3430 14.9923 0.0133 0.0137 -1.3195 3.0363 -1.8641 0.0370 -0.0195 0.0133 0.0068 0.0019 -0.3464 0.5923 -0.3371 -0.0006 -0.0158 0.0137 0.0019 0.0007 -0.1082 0.1855 -0.1052 -0.0004 1.6321 -1.3195 -0.3464 -0.1082 18.4807 -31.9525 18.2288 0.0378 -3.1590 3.0363 0.5923 0.1855 -31.9525 70.8923 -43.9858 0.0014 1.8352 -1.8641 -0.3371 -0.1052 18.2288 -43.9858 27.9232 -0.0163 0.0046 0.0370 -0.0006 -0.0004 0.0378 0.0014 -0.0163 0.0078 Bootstrap estimate of SE for overall proportion of sites occupied = 0.0883 CPU time= 18 seconds (0.30 min) P 23 PHỤ LỤC 5: HÌNH ẢNH NGHIÊN CỨU TẠI THỰC ĐỊA Giám sát rừng tràm Giám sát rừng phi lao Giám sát bụi Giám sát bãi cát Giám sát rừng phi lao Giám sát lùm ven biển P 24 Giám sát bãi đất trống Tác giả thu mẫu Đánh dấu vị trí giám sát Đánh dấu vị trí giám sát Hang Nhơng cát sọc Dụng cụ thực địa P 25 PHỤ LỤC 6: HÌNH ẢNH NHÔNG CÁT SỌC TRÊN THỰC ĐỊA Trên cỏ Trong bụi dứa dại Trên bờ tường mộ Dưới gốc phi lao Dưới bụi Trên cát P 26 PHỤ LỤC 7: HÌNH ẢNH VỀ CÁC TÁC ĐỘNG CỦA CON NGƯỜI VÀ MÔI TRƯỜNG ĐẾN QUẦN THỂ NHÔNG CÁT SỌC Xây dựng khu nghỉ dưỡng sinh thái Khai thác gỗ rừng phi lao Rừng phi lao sau bão lớn Đốt rừng tràm sau khai thác gỗ Đào hang bắt Nhông cát sọc Đặt bẫy bắt Nhông cát sọc P 27 ... HỌC HUẾ TRƯỜNG ĐẠI HỌC SƯ PHẠM CAO THỊ THANH NGUYÊN NGHIÊN CỨU MẬT ĐỘ, BIẾN ĐỘNG SỐ LƯỢNG VÀ XÁC SUẤT PHÁT HIỆN LỒI NHƠNG CÁT SỌC (Leiolepis guentherpetersi) Ở VÙNG CÁT VEN BIỂN HUYỆN PHÚ LỘC, TỈNH... ? ?Nghiên cứu mật độ, biến động số lượng xác suất phát lồi Nhơng cát sọc (Leiolepis guentherpetersi) vùng cát ven biển huyện Phú Lộc, tỉnh Thừa Thiên Huế? ?? Mục đích nghiên cứu Đánh giá mật độ cá thể thông... mùa môi trường sống khác xác suất phát lồi Nhơng cát sọc (Leiolepis guentherpetersi) vùng cát ven biển huyện Phú Lộc, tỉnh Thừa Thiên Huế tạo sở khoa học để đề xuất môi trường sống cần ưu tiên