thuộc nhóm G. luii phân bố ở đảo Cát Bà (Việt Nam) (Hình 1.6). Loài Thạch sùng mí cát bà (G. catbaensis) được mô tả dựa trên mẫu chuẩn thu tại đảo Cát Bà, thành phố Hải Phòng (Ziegler và cs, 2008) [68] và hiện chỉ ghi nhận tại Việt Nam. Gần đây, nghiên cứu của Ngo và cs (submitted) [39] ghi nhận mở rộng vùng phân bố của loài Thạch sùng mí cát bà tại các đảo nhỏ thuộc Vịnh Hạ Long. Nghiên cứu của Ngo và cs (2016) [42] bước đầu đã ước tính được kích cỡ quần thể và ghi nhận một số đặc điểm sinh học, sinh thái cũng như xác định những mối đe dọa tới quần thể loài Thạch sùng mí cát bà tại đảo Cát Bà. Trên cơ sở đó, Nguyen và cs (2016) [46] đề xuất đưa loài Thạch sùng mí cát bà vào trong danh lục đỏ thế giới xếp hạng Nguy cấp (EN).
Hình 1.6 Loài Thạch sùng mí cát bà (Goniurosaurus catbaensis)
Đặc điểm nhận dạng: SVL 84,4 - 111,5 mm; cơ thể và chân mảnh; 5 - 6 vảy quanh mũi; 1 hàng vẩy có nốt sần lớn phía trước ổ mắt; bề mặt mí mắt trên là các vẩy có kích thước bằng kích thước các vẩy ở đỉnh đầu; không có vảy gian mũi; 8 - 9 vảy môi dưới; 8 - 10 vảy môi trên; 11 - 12 vẩy trước ổ mắt; 52 - 55 vảy mí mắt; 2 - 3 vảy sau cằm; vùng dưới hàm không có vẩy có nốt sần lớn; 112- 127 vảy quanh thân; 23 - 25 hàng vảy có nốt sần quanh thân; 8 - 11 vảy nhỏ bao quanh vảy lớn; 22 - 24 bản mỏng dưới ngón 4; 19 - 21 lỗ đùi phía trước hậu môn ở con đực; mống mắt màu nâu vàng; phần trên của đầu, lưng và chân có màu
xám nâu và các đốm màu xám đen; dải màu phía sau gáy hẹp, hình chữ V; trên lưng có 3 - 4 dải sáng màu vàng, hẹp không có các đốm đen; đuôi màu đen có 5 dải màu trắng (đuôi nguyên vẹn); mặt bụng của đầu, cơ thể, chân màu trắng đục; vùng cổ chấm màu nâu (Ziegler và cs, 2008) [68].
Một số đặc điểm sinh học, sinh thái: Ngoài tự nhiên loài này hoạt động mạnh vào buổi tối, thường bắt gặp loài này bám trên các vách đá, các cửa hang đá, cách mặt đất từ 0,3 - 1,5 m hay các đường mòn trong rừng. Trong nuôi nhốt ban ngày chúng thường trú trong các ống tre hay gáo dừa, khe tối; ra kiếm ăn và uống nước vào buổi tối, ăn các loại côn trùng như dế, sâu quy, sâu sáp, mối… có trộn thêm canxi. Loài này đẻ mỗi lần 1 - 2 trứng (kích thước 24 - 25 mm, rộng 15,8 - 16,5 mm) bắt đầu vào khoảng tháng 4 đến tháng 7, đẻ 2 - 3 lứa, có thể do thiếu chất nên trứng thường bị non, hoặc trứng bị hỏng sau một thời gian ấp (Ziegler và cs. 2008) [68].
Chương 2
MỤC TIÊU, NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU
2.1. Mục tiêu nghiên cứu
- Phân tích được thành phần thức ăn của loài Thạch sùng mí lichtenfer (Goniurosaurus lichtenfelderi) tại VQG Bái Tử Long và loài Thạch sùng mí cát bà (G. catbaensis) tại Vịnh Hạ Long.
- Đánh giá các mối đe dọa và đề xuất các giải pháp bảo tồn hai loài Thạch sùng mí tại Việt Nam.
2.2. Nội dung nghiên cứu
Nội dung 1: Xác định thành phần thức ăn của loài Thạch sùng mí lichtenfer và loài Thạch sùng mí cát bà.
- Xác định thành phần thức ăn, dạng thức ăn ưa thích của loài.
- So sánh thành phần thức ăn theo cấu trúc giới tính và nhóm tuổi.
- So sánh thành phần thức ăn của loài Thạch sùng mí lichtenfer với loài Thạch sùng mí cát bà.
Nội dung 2: Đánh giá các mối đe dọa đến loài và đề xuất các giải pháp bảo tồn loài Thạch sùng mí lichtenfer và loài Thạch sùng mí cát bà tại Việt Nam.
- Đánh giá các mối đe dọa đến loài.
- Đề xuất và kiến nghị các biện pháp bảo tồn sinh cảnh và các quần thể loài.
2.3. Phương pháp nghiên cứu
2.3.1. Dụng cụ khảo sát thực địa
Các dụng cụ phục vụ cho công tác điều tra thực địa gồm có: máy định vị GPS Garmin 62s, thước đo điện tử độ chính xác 0,01 mm, phiếu giám sát, máy đo nhiệt độ, độ ẩm, máy ảnh, đèn đội đầu, ống tiêm, xilanh, ống thụt dạ dày,
̛
2.3.2. Khảo sát thực địa - Khảo sát theo tuyến
- Khảo sát theo tuyến:
Nguyên tắc lập tuyến: Tuyến điều tra được lập dựa vào bản đồ địa hình, thảm thực vật và sinh cảnh sống của thạch sùng mí. Các tuyến điều tra đi qua các dạng sinh cảnh, độ cao khác nhau của khu vực nghiên cứu. Đối với loài Thạch sùng mí lichtenfer thường ghi nhận gần các điểm có suối trong rừng, do đó nghiên cứu khảo sát dọc theo các tuyến suối núi đất. Trong khi đó, chúng tôi lựa chọn khảo sát loài Thạch sùng mí cát bà dọc theo các tuyến có các hang và vách đá, các thung lũng giữa các dãy núi đá vôi trong rừng.
Qua phân tích tài liệu và phỏng vấn cán bộ kiểm lâm tại VQG Bái Tử Long và Vịnh Hạ Long, chúng tôi xác định và tiến hành khảo sát thực địa theo 6 tuyến:
- VQG Bái Tử Long, huyện Vân Đồn, tỉnh Quảng Ninh: Tuyến 1 (T-1): tuyến đảo Ba Mùn,
Tuyến 2 (T-2): tuyến đảo Cái Lim,
- Vịnh Hạ Long, thành phố Hạ Long, tỉnh Quảng Ninh Tuyến 3 (T-3): Áng Cống Đỏ,
Tuyến 4 (T-4): Áng Bù Xám, Tuyến 5 (T-5): Áng Cửa Vạn, Tuyến 6 (T-6): Áng Dù.
Mỗi tuyến khảo sát có chiều dài khoảng từ 2 - 3km. Trong đó, tuyến 1 và tuyến 2 khảo sát tại dọc theo các tuyến suối ở khu vực rừng núi đất tại các đảo thuộc VQG Bái Tử Long (Hình 2.1.A, B) để ghi nhận loài Thạch sùng mí lichtenfer (Hình 1.5, trang 14); các tuyến còn lại khảo sát dọc theo các tuyến đường mòn ở khu vực rừng núi đá vôi tại các đảo nhỏ thuộc Vịnh Hạ Long (Hình 2.1.C, D) để ghi nhận loài Thạch sùng mí cát bà (Hình 1.6, trang 15).
Hình 2.1 A, B. Sinh cảnh Vườn Quốc Gia Bái Tử Long; C, D.
Sinh cảnh Vịnh Hạ Long
- Phương pháp thu mẫu:
Để đo đếm các chỉ tiêu hình thái, thụt rửa dạ dày, chúng tôi tiến hành thu thập mẫu vật của loài Thạch sùng mí dọc theo các tuyến khảo sát. Mẫu vật được thu thập bằng tay và các dụng cụ chuyên dụng như kẹp có bọc cao su nhằm tránh gây tổn thương đến con vật. Mẫu vật sau khi thu được mẫu thức ăn và đo đếm được thả lại đúng điểm đã thu thập.
- Thời gian khảo sát thực địa:
Chúng tôi tiến hành khảo sát thực địa trong khoảng thời gian 20 ngày từ tháng 6 đến tháng 7 năm 2017, vào thời gian này các loài Thạch sùng mí ra hoạt động, kiếm ăn và sinh sản (Ngo và cs. 2016) [42]. Nghiên cứu đã thực hiện 2 đợt khảo sát tại VQG Bái Tử Long và Vịnh Hạ Long.
Đợt 1 từ ngày 21 đến ngày 28 tháng 6 năm 2017 tại VQG Bái Tử Long. Đợt 2 từ ngày 29 tháng 6 đến ngày 10 tháng 7 năm 2017 tại Vịnh Hạ Long.
Hình 2.2 A, B. Khảo sát thực địa tại VQG Bái Tử Long và Vịnh Hạ Long
Để thu thập mẫu thức ăn và các số liệu hình thái có liên quan, chúng tôi tiến hành khảo sát thực địa cả ban ngày và ban đêm.
Ban ngày: Từ 13h00 đến 15h00, khảo sát các tuyến đường và chụp ảnh sinh cảnh, quan sát các nhân tố tác động tới sinh cảnh và loài, thả cá thể bắt từ tối hôm trước (sau khi đo đếm số liệu hình thái).
Ban đêm: Loài Thạch sùng mí thường hoạt động và kiếm ăn nên chúng tôi tiến hành khảo sát từ 20h00 đến 3h00 ngày hôm sau để thu thập mẫu vật và đo đếm các chỉ số hình thái và sinh thái liên quan.
2.3.3. Đặc điểm hình thái
Xây dựng bảng dữ liệu hình thái cho 02 loài Thạch sùng mí cát bà và Thạch sùng mí Lichtenfer theo Ziegler và cs (2008), và Orlov và cs (2008) [55], [56], [68]. Các số đo được đo bằng thước kẹp điện tử Alpha-Tool với độ chính xác là 0,01mm. Đo trọng lượng của từng cá thể bằng cân điện tử với độ chính xác là 0,01g (Bảng 1). Giới tính được xác định dựa trên đặc điểm hình thái ngoài đối với mẫu ghi nhận trên thực địa: cá thể đực thường có phần sau gốc đuôi phình to hơn hẳn so với cá thể cái (Hình 2.3.A) và trong mùa sinh sản cá thể đực thường có lỗ trước hậu môn rõ hơn và cá thể cái thường có trứng trong ổ bụng (Ziegler và cs. 2008) [68] (Hình 2.3.B).
Hình 2.3 Xác định giới tính dựa trên đặc điểm hình thái
Bảng 2.1: Các đặc điểm hình thái đo đạc trên mẫu Thạch sùng mí cát bà
STT Kí hiệu Giải thích
1 SVL Chiều dài từ mút mõm đến lỗ huyệt: đo từ mút mõm đến rìa trước hậu môn
2 TAL Chiều dài đuôi: đo từ rìa sau lỗ huyệt đến mút đuôi
3 WT Chiều rộng đuôi: đo gần cuống đuôi nơi chiều rộng lớn nhất 4 HL Dài đầu: Đo từ mút mõm đến góc sau của xương hàm dưới 5 HW Rộng đầu: Đo phần rộng nhất của đầu
6 HH Chiều cao tối đa của đầu 7 MW Chiều rộng miệng
8 JL Chiều dài miệng
9 AG Chiều dài ổ bụng: Từ hốc nách chi trước đến hốc háng chi sau 10 WS Chiều rộng ổ bụng: Kích thước rộng bụng lớn nhất
11 Forearm Chiều dài chi trước: đo từ hốc nách đến mút ngón 4 12 Hindlimb Chiều dài chi sau: đo từ hốc háng đến mút ngón 4 13 IN Vẩy giữa 2 vảy mũi lớn
14 PP Lỗ huyệt trước huyệt 15 Weight Cân nặng
Chúng tôi đo đếm các đặc điểm hình thái cơ bản (15 chỉ tiêu) sau đó nhập vào phiếu giám sát như ở Bảng 2.1.
2.3.4.1 Thu thập mẫu thức ăn trong dạ dày
Để nghiên cứu thành phần thức ăn nhưng không gây chết đối với mẫu vật, chúng tôi áp dụng phương pháp thụt dạ dày. Phương pháp này được áp dụng khá phổ biến trong nghiên cứu thành phần thức ăn các loài bò sát và lưỡng cư (James và cs, 1984, Norval và cs, 2012, Ngo và cs, 2016) [28], [42], [50]. Mẫu vật sau khi thu thập sẽ được thụt dạ dày ngay tại vị trí thu mẫu. Chúng tôi không thụt dạ dày với những cá thể cái có trứng và con non vì có thể ảnh hưởng tới khả năng sinh sản hoặc sinh trưởng. Nước sử dụng cho việc thụt dạ dày là nước tinh khiết để đảm bảo tinh sạch. Phương pháp thực hiện dựa trên mô tả của Sole và cs (2005) [59]: giữ cá thể cố định, dùng panh nhỏ có quấn cao su để mở miệng mẫu vật, nhẹ nhàng luồn ống dẫn nước có đường kính 1 mm qua thực quản xuống dạ dày, dùng bơm tiêm nước để thức ăn trào ra ngoài. Lượng nước và thức ăn trào ra ngoài được hứng bởi cốc đựng nước (thể tích 20 ml) có màng lọc (Hình 2.4. A, B). Toàn bộ mẫu thức ăn thu được sau khi thụt dạ dày sẽ được bảo quản trong cồn 70% để phục vụ phân tích, định loại tại Bảo tàng Thiên nhiên Việt Nam và Viện Sinh thái và Tài nguyên Sinh vật.
Hình 2.4 A. Thụt dạ dày mẫu Thạch sùng mí cát bà ngoài thực địa.; B. Thụt dạ dày mẫu Thạch sùng mí lichtenfer
Các mẫu thức ăn được phân tích và định loại dưới kính lúp soi nổi Leica S6E ở phòng thí nghiệm của Bảo tàng Thiên nhiên Việt Nam và Viện Sinh thái và Tài nguyên sinh vật. Mẫu thức ăn được đo đếm, chụp ảnh dưới kính lúp Leica M80 có gắn máy ảnh Leica IC80HD. Nghiên cứu định loại mẫu thức ăn theo các tài liệu của Millar và cs (2000) và Thái Trần Bái (2003) [2], [34].
Các thông số phân tích thành phần thức ăn gồm: Tần số (F) thể hiện sự phong phú của một loại con mồi thu được từ các mẫu dạ dày, số lượng (N) là số mẫu thức ăn của một loại con mồi, thể tích (V, mm3) của mỗi mẫu thức ăn được ước tính bằng công thức theo tài liệu của Magnusson và cs (2003) và Ngo và cs (2014) [40]:
Trong đó, L: là chiều dài của mẫu thức ăn, W: là chiều rộng của mẫu thức ăn (Hình 11) Hirai và Masui 2001 [25].
Chỉ số quan trọng (Ix) là phần trăm của mỗi mẫu thức ăn được tính theo công thức mô tả trong tài liệu của Caldart và cs (2012), Ngo và cs (2014), Ngo và cs (2016), [15], [40], [42].
Hình 2.5 Xác định và đo đếm kích thước mẫu thức ăn dưới kính lúp soi nổi Leica S6E
của loài. Trong đó: D là chỉ số đa dạng, ni là số lượng mẫu thức ăn trong một mục thức ăn (Taxa) thứ i, N là tổng số lượng mẫu thức ăn của các mục thức ăn (Taxa), S là tổng số mục thức ăn (Taxa). Chỉ số đa dạng được trình bày dưới dạng 1/D (0 ≤ D ≤ 1), khi D càng lớn thì đa dạng càng cao. Công thức tính như sau:
Để ước tính mức độ đồng đều giữa các loại thức ăn của Thạch sùng mí chúng tôi sử dụng chỉ số Shannon’s evenness. Trong đó: E là chỉ số đồng đều (0 ≤ E ≤ 1), khi E = 1 thì độ đồng đều cao nhất, H’là chỉ số đa dạng Shannon - Weiner, Hmax = ln S (S là tổng số các mục thức ăn - taxon). Công thức tính như sau:
Chỉ số đa dạng H’ được tính như sau:
Trong đó: pi là tỉ lệ số lượng mẫu thức ăn trong một mục thức ăn (ni) trên tổng số lượng mẫu thức ăn của các mục thức ăn - Taxa (N) ở bậc thứ i.
Để đánh giá mức độ tương đồng về thành phần thức ăn giữa các loài trong cùng 1 giống (giữa loài Goniurosaurus catbaensis và loài G. lichtenfelderi), nghiên cứu phân tích sự chồng chéo trong mức độ tương đồng giữa chế độ ăn giữa 2 loài, chúng tôi sử dụng chỉ số trùng lặp Pianka (Ojk) Pianka, 1973 [56] với các giá trị từ 0 (không chồng chéo) đến 1 (hoàn toàn chồng chéo) Krebs, 1999 [29], theo biểu thức sau:
lệ của thành phần thức ăn i trong tổng nguồn thức ăn được sử dụng bởi các loài j; Pik là tỷ lệ của thành phần thức ăn i trong tổng nguồn thức ăn được sử dụng bởi các loài k; n là tổng số các loại thức ăn cho các loài j và k.
2.3.4.3 Phân tích thống kê
Để thực hiện phân tích thống kê, xây dựng biểu đồ và thuật toán kiểm định, nghiên cứu sử dụng các phần mềm thống kê như phần mềm SPSS 16.0, Past Statistic Hammer etal. (2001), Stagraphics và Excel.
Để kiểm định sự sai khác về hình thái theo giới tính, nhóm tuổi và giữa 2 loài Thạch sùng mí trong nghiên cứu, chúng tôi áp dụng phân tích tổng hợp các nhân tố PCA (Principal Component Analysis) trên phần mềm PAST thông qua kiểm định T-test giá trị PC1. Giá trị Loading-PC1 trong phân tích tổng hợp được tính toán để xác định nhân tố hình thái quan trọng trong sự sai khác. Ngoài ra, nghiên cứu kiểm định sự khác nhau về kích thước chiều rộng, chiều dài và thể tích thành phần thức ăn theo giới tính và nhóm tuổi, chúng tôi sử dụng thuật toán T-test giữa 2 yếu tố Con đực và Con cái của loài Thạch sùng mí lichtenfer và thuật toán One-way ANOVA giữa 3 yếu tố Con đực, Con cái và Con non. Các thuật toán kiểm định được thực hiện trên phần mềm SPSS 16.0.
Để phân tích và đánh giá mối tương quan trong quan hệ giữa hình thái của loài Thạch sùng mí với kích thước thức ăn lựa chọn, nghiên cứu sử dụng phân tích hồi quy qua hàm 1 biến (1 nhân tố): Kích thước thức ăn (Y) = Kích thước hình thái (X) x Slope (độ dốc a) + intercept (b). Kiểm định one-way ANOVA để đánh giá mối tương quan chặt giữa các kích thước, giá trị Slope biến thiên nếu: Slope > 0 tỷ lệ thuận; Slope < 0 tỷ lệ nghịch. Sự khác nhau đáng kể được công nhận khi chỉ số ý nghĩa (P-value) có giá trị <0,05. Phân tích hồi quy được thực hiện trên phần mềm Stagraphics.
Chương 3
KẾT QUẢ NGHIÊN CỨU
3.1. Đặc điểm hình thái của hai loài Thạch sùng mí lichtenfer và Thạch sùng mí cát bà
Nghiên cứu tiến hành khảo sát các tuyến suối thuộc VQG Bái Tử Long ghi nhận 20 cá thể Thạch sùng mí lichtenfer, trong đó 19 cá thể ghi nhận tại tuyến đảo Ba Mùn và 01 cá thể thuộc tuyến đảo Cái Lim. Do địa hình các tuyến suối phức tạp và các cá thể Thạch sùng mí lichtenfer sống ẩn nấp trong các khe đá