Mô hình hóa vùng phân bố tiềm năng của loài rùa hồ gươm (rafetus swinhoei (gray, 1873))

LỜI CẢM ƠN Đƣợc đồng ý Ban Giám Hiệu trƣờng Đại Học Lâm Nghiệp, Khoa Quản Lý Tài Nguyên Rừng Môi Trƣờng thầy PGS.TS Vũ Tiến Thịnh, thầy Ths Trần Văn Dũng, tơi tiến hành khóa luận tốt nghiệp: ―Mơ hình hóa vùng phân bố tiêm loài Rùa hồ gƣơm (Rafetus swinhoei)‖ Nhân dịp này, cho phép tơi đƣợc bày tỏ lịng biết ơn sâu sắc tới hƣớng dẫn giúp đỡ nhiệt tình thầy PGS.TS Vũ Tiến Thịnh, thầy Ths Trần Văn Dũng quan tâm giúp đỡ thầy cô môn Động vật rừng thuộc Khoa Quản Lý Tài Nguyên Rừng Môi Trƣờng Đến khóa luận hồn thành, tơi xin chân thành cảm ơn tới giúp đỡ quý báu Mặc dù cố gắng trình thực nhƣng kiến thức, kinh nghiệm thân, điều kiện thời gian, nhƣ tài liệu tham khảo cịn hạn chế nên khóa luận khơng tránh khỏi thiếu sót định Kính mong nhận đƣợc ý kiến quý báu góp ý, bổ sung thầy giáo bạn, để khóa luận hoàn thiện Xin chân thành cảm ơn! Hà Nội, ngày 10 tháng 05 năm 2018 Sinh viên thực Nguyễn Hữu Khánh i MỤC LỤC LỜI CẢM ƠN i MỤC LỤC ii DANH MUC TỪ VIẾT TẮT iv DANH MỤC BẢNG v DANH LỤC HÌNH vi ĐẶT VẤN ĐỀ CHƢƠNG I TỔNG QUAN VỀ ĐỀ TÀI NGHIÊN CỨU 1.1 Hệ thống phân loại Rùa hồ gƣơm 1.1.1 Lớp bò sát (Reptilia) 1.1.2 Bộ Rùa (Testudines) 1.1.3 Họ Rùa Mai Mềm hay Ba Ba (Trionychidae) 1.2 Rùa hồ gƣơm (Rafetus swinhoei) 1.3 Mơ hình Entropy cực đại (MaxEnt) CHƢƠNG II MỤC TIÊU, ĐỐI TƢỢNG, NỘI DUNG VÀ PHƢƠNG PHÁP 12 NGHIÊN CỨU 12 2.1 Mục tiêu nghiên cứu 12 2.2 Đối tƣợng nghiên cứu 12 2.3 Nội dung nghiên cứu 12 2.4 Phƣơng pháp nghiên cứu 12 2.4.1 Phƣơng pháp kế thừa số liệu thứ cấp 12 2.4.2 Phƣơng pháp xử lí số liệu phần mềm MaxEnt 16 CHƢƠNG III ĐẶC ĐIỂM KHU VỰC NGHIÊN CỨU 19 3.1 Điều kiện tự nhiên Việt Nam 19 3.2 Khu vực Đông Bắc 22 3.3 Khu vực Tây Bắc 23 3.4 Khu vực đồng sông Hồng 25 3.5 Khu vực Bắc Trung Bộ 27 CHƢƠNG IV KẾT QUẢ VÀ THẢO LUẬN 29 4.1 Đánh giá tính xác mơ hình 29 ii 4.2 Vùng phân bố tiềm loài Rùa hồ gƣơm 30 4.2.1 Vị trí vùng phân bố tiềm lồi Rùa hồ gƣơm 30 4.2.2 Diện tích vùng phân bố tiềm lồi Rùa hồ gƣơm 32 4.3 Các yếu tố tự nhiên ảnh hƣởng mạnh đến vùng phân bố tiềm loài Rùa hồ gƣơm 33 4.4 Đề xuất nghiên cứu, bảo tồn loài Rùa hồ gƣơm 34 4.4.1 Đề xuất khu vực ƣu tiên tìm kiếm lồi Rùa hồ gƣơm Việt Nam 34 4.4.2 Đề xuất giải pháp bảo tồn loài Rùa hồ gƣơm 36 CHƢƠNG V KẾT LUẬN – TỒN T I – KIẾN NGH 38 5.1 Kết luận 38 5.2 Tồn 38 5.3 Kiến nghị 38 TÀI TIỆU THAM KHẢO PHỤ LỤC iii DANH MUC TỪ VIẾT TẮT Viết tắt Viết đầy đủ ATP Chƣơng trình Bảo tồn Rùa Châu Á AUC Vùng diện tích dƣới đƣờng cong CI Tổ the chức CurveBảo tồn thiên nhiên Quốc tế Công ƣớc Buôn bán Quốc tế Các loài Động, Thực vật Hoang dã CITES CR Nguy cấp Cực kì nguy cấp (Criticalt Endangered) eDNA Kỹ thuật gen môi trƣờng EN Nguy cấp (Endangered) ENMs Mô hình ổ sinh thái ICZN Luật Quốc tế Danh pháp Đông vật IMC Tổ chức bảo tồn Indo-Myanmar IUCN Tổ chức bảo tồn thiên nhiên giới MaxEnt Mô hình Entropy cực đại NDVI Chỉ số thực vật SĐVN Sách Đỏ Việt Nam TCC Trung tâm Bảo tồn Rùa Cúc Phƣơng TCF Quỹ bảo tồn rùa iv DANH MỤC BẢNG Bảng 2.1 Các biến liệu môi trƣờng đƣợc sử dụng mơ hình (các biến bơi đậm đƣợc sử dụng tính tốn cuối mơ hình) 15 Bảng 2.2 Thang phân chia mức độ thích hợp sau chạy 17 phần mềm Maxent 17 Bảng 4.1 Bảng mức độ tham gia biến khí hậu 33 v DANH LỤC HÌNH Hình 1.1 Rùa hồ gƣơm bị bắt Đồng Mô, Sơn Tây, Hà Nội Hình 1.2 Cá thể Rùa hồ gƣơm (giải Thƣợng Hải) Trung Quốc Hình 1.3 Cá thể Rùa hồ gƣơm hồ Đồng Mô Hình 2.1 Tọa độ điểm ghi nhận loài chuẩn bị cho phần mềm MaxEnt 17 Hình 2.2 Giao diện phần mềm MaxEnt nhập biến khí hậu 18 Hình 3.1 Bản đồ Việt Nam 21 Hình 4.1 Biểu đồ diện tích dƣới đƣờng cong (AUC) mơ hình MaxEnt cho loài Rùa hồ gƣơm 29 Hình 4.2 Bản đồ vùng phân bố thích hợp lồi Rùa hồ gƣơm 31 Hình 4.3 Biểu đồ diện tích phân bố thích hợp lồi Rùa hồ gƣơm 32 Hình 4.4 Biểu đồ ảnh hƣởng biến môi trƣờng đến vùng phân bố tiềm loài Rùa hồ gƣơm 34 Hình 4.5 Bản đồ khu vực ƣu tiên nghiên cứu, bảo tồn loài Rùa hồ gƣơm 35 vi ĐẶT VẤN ĐỀ Mơ hình hóa vùng phân bố (Species distribution Modeling) mơ hình ổ sinh thái cơng cụ sử dụng để giúp hiểu thêm vùng phân bố loài (Nazeri et al.,2012) Mơ hình ổ sinh thái (ENMs) phƣơng pháp sử dụng liệu ghi nhận vị trí xuất lồi kết hợp với liệu mơi trƣờng, từ tạo mơ hình thể đƣợc khu vực địa lý đáp ứng đƣợc yêu cầu sinh thái lồi ENMs dựa liệu phân bố loài tại, biến giá trị môi trƣờng ô lƣới sử dụng thuật tốn xử lý để xác định lƣới phù hợp hay khơng phù hợp với lồi (Raxworthy et al., 2007; Hirzel et al., 2002) MaxEnt phần mềm sử dụng phƣơng pháp dự đốn để mơ vùng phân bố tiềm loài từ thơng tin có (Phillips et al., 2006) Hiện nay, MaxEnt mơ hình phổ biến traong xây dựng mơ hình ổ sinh thái (phân bố) lồi có nhiều nghiên cứu sử dụng để mơ hình hóa vùng phân bố lồi thú Loài Rùa hồ gƣơm Rafetus swinhoei (Gray, 1873), loài rùa nƣớc lớn giới, Việt Nam đƣợc coi rùa huyền thoại, mang ý nghĩa to lớn Loài Rùa hồ gƣơm tên thƣờng gọi rùa Hồ Gƣơm, rùa Hoàn Kiếm hay giải Sin-hoe (Swinhoe’s softshell turtle), thuộc họ rùa Mai Mềm (Ba Ba) (Trionychidae), Rùa (Testudines) Những kết điều tra nghiên cứu cho thấy vùng phân bố loài hẹp với số lƣợng cá thể Năm 2011, Quỹ bảo tồn rùa (TCF) xếp Rùa hồ gƣơm 25 loài rùa bị đe dọa cao giới Tổ chức bảo tồn thiên nhiên quốc tế (IUCN) năm 2018 Sách Đỏ Việt Nam năm 2007 xếp loài Rùa hồ gƣơm vào mức nguy cấp (CR) Loài Rùa hồ gƣơm đứng đầu danh sách loài rùa đối mặt với nguy tuyệt chủng giới Việt Nam Hiện nay, lồi Rùa hồ gƣơm cịn lại cá thể ni Trung Quốc cá sống ngồi tự nhiên Việt Nam Hiện nay, có nghiên cứu lồi Rùa hồ gươm, nghiên cứu chủ yếu tập trung vào điều tra ghi nhận có mặt hay điều tra kích thƣớc quần thể, mật độ định danh loài Các nghiên cứu cung cấp thơng tin hữu ích đặc điểm sinh thái lồi nhƣng chƣa có nghiên cứu dự đoán vùng phân bố tiềm lồi Rùa hồ gƣơm nên thơng tin vùng phân bố chúng hạn chế Vậy nên, thực đề tài nghiên cứu khoa học: “Mơ hình hóa vùng phân bố tiềm lồi Rùa hồ gƣơm (Rafetus swinhoei (Gray, 1873))” Nghiên cứu có ý nghĩa to lớn lồi Rùa hồ gƣơm Qua nghiên cứu biết đƣợc nơi có điều kiện tự nhiên đáp ứng đƣợc thích nghi, tập tính chúng, sống đâu sinh trƣởng, phát triển trì nịi giống tốt Nghiên cứu hỗ trợ chuyên gia, nhà bảo tồn khoanh vùng tìm kiếm cá thể ngồi tự nhiên tìm giải pháp để bảo tồn loài Rùa hồ gƣơm quý CHƢƠNG I TỔNG QUAN VỀ ĐỀ TÀI NGHIÊN CỨU 1.1 Hệ thống phân loại Rùa hồ gƣơm 1.1.1 Lớp bò sát (Reptilia) - Lớp bò sát nhóm động vật ―biến nhiệt‖ nên chúng cần phải làm ấm cách sử dụng nhiệt bên ngồi (mặt trời … ) Lớp bị sát đa dạng hình dạng thể:  Dạng điển hình bị sát thấy thằn lằn cá sấu với phần đầu phần cổ rõ ràng, bốn chi dài khoẻ, nằm ngang nâng đƣợc thân khỏi mặt đất đuôi dài Một số lồi thằn lằn chun hố với đời sống có thêm màng da bên thân giúp việc nhảy chuyền từ cành sang cành khác (tắc kè bay) Một số thằn lằn sống chui luồng khe, hốc đất có chi tiêu giảm (liu điu)…  Rắn nhóm thằn lằn chuyên hóa đặc biệt với đời sống trƣờn đất có thân dài, thiếu chi  Nhóm rùa có dạng biến đổi thể đƣợc bảo vệ giáp xƣơng Cổ dài nhƣng thân đuôi tƣơng đối ngắn Một số lồi rùa nƣớc (vích, đồi mồi) có chi trƣớc biến thành bơi chèo, khác xa dạng chi năm ngón điển hình - Ngày nay, theo Benton (2004) Colin Tudge (2000) lồi thuộc lớp bị sát đƣợc ghi nhận bộ:  Bộ Thằn lằn Đầu mỏ (Rhynchocephalia) (các loài tuatara New Zealand): lồi  Bộ Có vảy (Squamata): gồm lồi thằn lằn, rắn amphisbaenia ("bị sát giống bọ"), có khoảng 7.900 loài  Bộ Cá sấu (Crocodylia): gồm loài cá sấu thực sự, cá sấu caiman cá sấu Mỹ, có 23 lồi  Bộ Rùa (Testudines): gồm lồi rùa, ba ba, vích, đồi mồi v.v., có khoảng 300 loài 1.1.2 Bộ Rùa (Testudines) Bộ rùa (Testudines) (Anquetin, 2009) lồi bị sát thuộc nhóm chỏm siêu Chelonia (hay Testudinata) Trong tiếng Việt, loài thuộc rùa đƣợc gọi nhiều tên khác nhƣ rùa, ba ba, giải, vích, đồi mồi Trong văn hóa, rùa đƣợc nhìn nhận dƣới nhiều gốc độ khác tùy theo văn hóa Rùa đƣợc xem vật mang nhiều ý nghĩa biểu trƣng văn hóa Bộ chia làm nhóm sinh thái nhóm sống mặt đất di chuyển chậm ăn thực vật nhóm sống dƣới nƣớc di chuyển nhanh, ăn thịt Tất loài rùa đẻ trứng thở phổi, xƣơng sống gắn liền với mai 1.1.3 Họ Rùa Mai Mềm hay Ba Ba (Trionychidae) Tên gọi tiếng Việt họ bị sát thuộc phân lớp Khơng cung (Anapsida), Rùa (Testudines) Thân dài đến m, chân có móng Phiến giáp bụng hở, khơng liền với mai lƣng Vỏ phủ lƣợt da mềm, hô hấp phụ da hỗ trợ cho phổi Đƣợc xác định chủ yếu ăn thịt, rùa nƣớc Họ đƣợc đặc trƣng vỏ thiếu biểu bì thay lớp phủ mềm da thuộc da Lớp vỏ thấp, trịn thƣờng thiếu xƣơng ngoại biên cận biên scutes (Carrol, 1988) 1.2 Rùa hồ gƣơm (Rafetus swinhoei) Năm 2011, Quỹ bảo tồn rùa (TCF) xếp Rùa hồ gƣơm 25 loài rùa bị đe dọa cao giới Tổ chức bảo tồn thiên nhiên quốc tế (IUCN) năm 2018 Sách Đỏ Việt Nam năm 2007 xếp loài Rùa hồ gƣơm vào mức nguy cấp (CR) Luật pháp bảo vệ: - Tất loài rùa cạn rùa nƣớc Việt Nam đƣợc bảo vệ Luật Bảo vệ Phát triển Rừng số 29/2004/QH11 Trong quy định đƣợc phép sở hữu bn bán lồi động vật hoang dã có nguồn gốc hợp pháp Tài liệu tiếng Anh: 10 Anquetin J., Barrett P M., Jones M E H., Moore-Fay S., Evans S E., (2009), A new stem turtle from the Middle Jurassic of Scotland: new insights into the evolution and palaeoecology of basal turtles, Proceedings of the Royal Society of London: Biological series: Vol 276, tr 879- 886, doi:10.1098/rspb.2008.1429 11 Benton, M J (2004), Vertebrate Paleontology, ấn lần thứ Blackwell Science Ltd 12 Bett N.N., Blair M.E., Sterling E.J (2012) Ecological Niche Conservatism in Doucs (Genus Pygathrix) - Int J Primatol 33: 972–988 13 Broxton P.D., Zeng X., Sulla-Menashe D., Troch P.A (2014) A Global Land Cover Climatology Using MODIS Data - J Appl Meteor Climatol 53: 1593 -1605 doi: http://dx.doi.org/10.1175/JAMC-D-13-0270.1 14 Broxton P.D., Zeng X., Sulla-Menashe D., Troch P.A (2014) A Global Land Cover Climatology Using MODIS Data - J Appl Meteor Climatol 53: 1593 -1605 doi: http://dx.doi.org/10.1175/JAMC-D-13-0270.1 15 Carrol, R.L, (1988), Vertebrate paleontology and evolution - 698 pages W.H Freeman & Company 16 Colin Tudge (2000) The Variety of Life Ấn Đại học Oxford ISBN 019-860426-2 17 Elith J ,(2000), Quantitative methods for modeling species habitat: comparative performance and an application to Australian plants In: Ferson S, Burgman M, editors Quantitative methods for conservation biology Springer, New York, pp 39–58 18 Farkas B and Webb R G (2003), ―Rafetus leloii Ha Dinh Duc, 2000 — a invalid species of softshell turtle from Hoan Kiem Lake, Hanoi, Vietnam (Reptilia, Testudines, Trionychidae),‖ Zool Abh., 53, 107 – 112 19 Farkas BL (1992), Wiederentdeckung eines Exemplars von Rafetus swinhoei (GRAY, 1873) in Naturhistorischen Museum Wien, Salamandra 28(2): 145-152 20 Fritz, U (2011) Order Testudines Batsch, 1788 In: Zhang, Z.-Q (Ed.), Animal biodiversity: An outline of higher-level classification and survey of taxonomic richness Zootaxa, 3148, 61–62 21 Gray JE (1873) XXII Notes on Chinese Mud-Tortaises (Trionychidae), with the Description of a new Species sent to the British Meseum by Mr Swinhoe, and Observations on the Male Organ of this Family The Annals and Magazine of Natural History (Fourth Series) 68: 156-160 22 Ha Dinh Duc ,(2000), A news pecies of the genus Rafetus (Family: Trionychidae) from Hoan Kiem lake of Vietnam Archaeology Magazine (Vietnam) 4:104 23 Hijmans R.J., Cameron S.E., Parra J.L., Jones P.G and Jarvis A (2005) Very high resolution interpolated climate surfaces for global land areas - Int J Climatol 25: 1965-1978 24 Hirzel, A.H., Hausser, J., Chessel, D., Perrin, N., (2002) Ecological-niche factor analysis: how to compute habitat-suitability maps without absence data? Ecology 83, 2027–2036 25 Kumar S., Stohlgren T.J (2009) Maxent modeling for predicting suitable habitat for threatened and endangered tree Canacomyrica monticola in New Caledonia - J Ecol Nat Environ 1: 094-098 26 Lauritoetal.-2005-A trionychidae (reptilia: testudines, cryptodira) from the Pliocene of costa rica, southern central America 27 Le D M and Pritchard P (2009), ―Genetic variability of the critically endangered softshell turtle, Rafetus swinhoei: a preliminary report,‖ in: Ngo D C., Ta H T., Le N N., Hoang X Q., Vo V P., Nguyen V S., Nguyen V T., Le T T., and Tran Q D (eds.), Proc of the 1st Natl Sci.Workshop on Amphibians and Reptiles in Vietnam, Hue Univ Press, Hue, pp 84 – 92 28 Le T B., Le Q H., Tran M L., Phan T H., Phan M T., Tran T T H, Pham T T., Nguyen D T., Nong V H., Phan V C., Dinh D K., Truong N H., and Ha D D (2010) ―Comparative morphological and DNA analysis of specimens of giant freshwater soft-shelled turtle in Vietnam related to Hoan Kiem turtle,‖ Tap chi Cong nghe Sinh hoc (J Biotechnol.), 8(3A), 949 – 954 29 Le T B., Le Q H., Tran M L., Phan T H., Phan M T., Tran T T H, Pham T T., Nguyen D T., Nong V H., Phan V C., Dinh D K., Truong N H., and Ha D D (2010) ―Comparative morphological and DNA analysis of specimens of giant freshwater soft-shelled turtle in Vietnam related to Hoan Kiem turtle,‖ Tap chi Cong nghe Sinh hoc (J Biotechnol.), 8(3A), 949 – 954 30 Merow C., Smith M.J., Silander J.A.Jr (2013) A practical guide to MaxEnt for modeling species’ distributions: what it does, and why inputs and settings matter, - Ecography 36: 1058–1069 doi: 10.1111/j.1600-0587.2013.07872.x 31 Meylan P A R G Webb (1988) (Rafetus swinhoei (Gray) 1873, a valid species of living soft-shelled turtle (family Trionychidae) from China Journal of Hepatology 22:118-119 32 Morrissey, J.F and J.L Sumich (2012) Introduction to the Biology of Marine Life Jones & Bartlett Learning 466pp 33 Nazeri M, Jusoff K, Madani N, Mahmud AR, Bahman AR et al., (2012) Predictive Modeling and Mapping of Malayan Sun Bear (Helarctos malayanus) Distribution Using Maximum Entropy PLoS ONE 7(10): e48104 doi:10.1371/journal.pone.0048104 34 Peterson A.T., Papes M., Eaton M (2007) Transferability and model evolution in ecological niche modeling: a comparison of GARP and MAXENT - Ecography 30: 550–560 35 Phillips S.J., Anderson R.P., and Schapire R.E (2006) Maximum entropy modeling of species geographic distributions - Ecol Model 190: 231-259 36 Phillips S.J., Dudík M., Schapire R.E (2017) [Internet] Maxent software for modeling species niches and distributions (Version 3.4.1) Available from url: http://biodiversityinformatics.amnh.org/open_source/maxent/ Accessed on 2017 - 10 - 37 Pritchard PCH (2001) Observation of body size, sympatry, and niche divergence in softshell turtles (Trionychidae) Chelonian Conservation and Biology 4(1): 5-27 38 Raxworthy CJ, Ingram C, Rabibisoa N, Pearson R (2007) Species delimitation applications for ecological niche modeling: a review and empirical evaluation using Phelsuma day Gecko Groups from Madagascar Systematic Biology 56:907–923 PHỤ LỤC Các biểu đồ ảnh hƣởng biến môi trƣờng đến vùng phân bố tiềm loài rùa Hồ gƣơm Ảnh thực địa Ảnh hồ Đồng Mô, Sơn Tây, Hà Nội ngày 17/03/2018 Ảnh hồ Hoàn Kiếm, Hoàn Kiếm, Hà Nội ngày 25/03/2018 Bảng tọa độ điểm ghi nhận có mặt loài Rùa hồ gƣơm Tọ Độ Điểm G i N ận Loài Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei X Y 21.074332 105.465900 21.077676 105.465857 21.048477 105.479770 Địa Điểm Nguồn Ghi Đang 21.068265 21.060737 105.465771 105.477280 Hồ Đồng Mô, Sơn Tây, Hà Nội cá thể ATP (Asian Turtle Program) sống tự nhiên 21.073391 105.466157 21.075398 105.473760 21.074512 105.468904 Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei 21.072710 105.460321 21.069867 105.477873 21.134771 105.436916 21.125885 105.434770 21.123043 105.435585 21.125404 105.438547 21.124844 105.442237 21.028054 105.852408 21.027693 105.852461 cá thể sống tự nhiên Hồ Xuân Khanh, Sơn Tây, Hà Nội Hồ Hoàn Kiếm, Hoàn Kiếm, Hà Nội vừa đƣợc ATP (Asian Turtle Program) Ha Dinh Duc (2000) phát Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei 21.027768 105.852048 21.031203 105.851415 21.031023 105.853014 21.026642 105.852418 21.026171 105.852847 21.159897 105.378148 Hồ Suối Hai, Ba Vì, Hà Nội 21.156375 105.372139 21.173316 105.335164 Đầm Long, Ba Vì, Hà Nội 21.175600 105.332458 Lê et al (2010); Farks Webb (2003) Lê et al (2010); Farks Webb (2003) Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei 21.930171 104.609361 Ngã Ba Văn Yên, Yên Bái 21.947274 104.568567 21.795418 104.789702 Trấn Yên, Yên Bái 21.766870 104.793076 21.567656 105.012809 Đầm Ao Châu, Phú Thọ 21.568914 105.012874 19.926332 105.400600 19.922224 105.415368 20.020693 105.465809 Lam Kinh, Thọ Xuân, Thanh Hóa Quảng Phú, Thọ Xuân, Thanh Hóa Lê et al (2010); Farks Webb (2003) Lê et al (2010); Farks Webb (2003) Lê et al (2010); Farks Webb (2003) Lê et al (2010); Farks Webb (2003) Lê et al (2010); Farks Webb (2003) Rafetus swinhoei Rafetus swinhoei Rafetus swinhoei 20.009154 105.471649 20.800846 105.347028 Quỳnh Lâm, Hịa Bình 20.803951 105.350053 Lê et al (2010); Farks Webb (2003)