Đa dạng về thành phần loài bò sát theo họ

Một phần của tài liệu (LUẬN văn THẠC sĩ) đa dạng các loài bò sát (reptilia) và ếch nhái (amphibia) tại vườn quốc gia cúc phương, tỉnh ninh bình​ (Trang 45)

Qua hình 4.1 cho thấy, trong 19 họ bò sát ở VQG Cúc Phƣơng thì họ Rắn nƣớc Colubridae đa dạng nhất với 22 loài (chiếm 27,5% tổng số loài bò sát ghi nhận

ở VQG Cúc Phƣơng). Có 6 họ chỉ có 1 loài (chiếm 1,25% tổng số loài ghi nhận gồm: họ Thằn lằn giun (Dibamidae), họ Thằn lằn chính thức (Lacertidae), họ Kỳ đà (Varanidae), họ Trăn (Pythonidae), họ Rắn mống (Xenopeltidae), họ Rắn hổ đất nâu (Lamprophiidae), họ Rùa núi (Testunidae).

Hình 4.1. Đa dạng thành phần loài bò sát theo các họ

4.1.3. Mô tả các loài bò sát ghi nhận mới tại VQG Cúc Phương

So sánh với công bố trƣớc đây đề tài đã ghi nhận đƣợc 5 loài bò sát mới cho VQG Cúc Phƣơng (hình 4.2), bao gồm: họ Tắc kè có 2 loài (Cyrtodactylus bobrovi, Cyrtodactylus soni); họ Rắn nƣớc có 1 loài (Lycodon rufozonatus); Họ Rắn ri có 1 loài (Myrrophis chinensis) và họ Rắn lục có 1 loài (Protobothrops cornutus). Riêng loài Cyrtodactylus soni có 2 mẫu lƣu trữ tại Bảo tàng VQG Cúc Phƣơng nhƣng một số thông tin về điểm thu mẫu, sinh cảnh chƣa rõ ràng nên đề tài không tiến hành mô tả chi tiết ở đây.

Hình 4.2. Các loài bò sát ghi nhận mới cho VQG Cúc Phƣơng Họ Tắc kè Gekkonidae

1. Thằn lằn chân ngón bobrovi - Cyrtodactylus bobrovi Nguyen, Le,

Pham, Ngo, Hoang, Pham & Ziegler, 2016 (Hình 4.3).

Mẫu phân tích: Một con cái trƣởng thành NHQ.17.122. thu ở tọa độ 20021'957’’ vĩ Bắc 105030'064’’ kinh Đông, ngày 30/9/2017 tại xã Thạch Lâm, huyện Thạch Thành, tỉnh Thanh Hóa.

Đặc điểm hình thái: Con cái trƣởng thành, chiều dài thân (SVL) 86.5 mm, chiều dài đuôi (TaL) 64.7 mm (đuôi mới), đầu thon dài, đầu dài hơn rộng (HW/HL 0.69), khác biệt rõ ràng với cổ, vùng trƣớc mắt lõm, mõm dài (SE/HL 0.46), dài hơn đƣờng kính ổ mắt (OD/SE 0.51); tấm trên mõm nhỏ; mắt to (OD/HL 0.23), đồng t thẳng đứng; tai nhỏ, hình trái xoan (ED/HL 0.06); vảy mút mõm rộng hơn dài; lỗ mũi tròn; tấm cằm hình tam giác; có 2 vảy sau cằm, mở rộng; vảy môi trên 11-10;

vảy môi dƣới 11-10. Vảy lƣng nhỏ hình hạt, có các nốt sần hình tròn xuất hiện ở đầu và lƣng, mỗi nốt sần đƣợc bao quanh bởi 9-10 vảy nhỏ; vảy bụng mịn, tròn, trung bình lớn hơn 2 đến 3 lần vảy lƣng, 41 hàng vảy quanh thân; phần bên thân không có nốt sần và nếp gấp giữa thân và bụng không rõ ràng; 198 vảy bụng từ cằm đến lỗ huyệt; xuất hiện vảy đùi lớn dần; có 7 lỗ trƣớc huyệt xếp theo hình chữ V; không xuất hiện lỗ đùi; có 2 nốt sần lớn ở mỗi bên phía sau lỗ huyệt, giữa các ngón không có màng; có 21 tấm dƣới ngón thứ 4 của chi trƣớc và 24 tấm dƣới ngon thứ 4 của chi sau.

Màu sắc mẫu sống: lƣng màu xám hơi nâu, mặt trên đầu có các chấm đen nhỏ, có 2 sọc vằn nhỏ màu nâu xuất hiện ở vùng mũi; xuất hiện 1 sọc rộng màu tối kéo dài từ góc phía sau của mắt đi qua màng nhĩ và tới hai đốm đen lớn ở trên cổ, có viền màu trắng sữa; trên thân có năm vằn ngang rộng màu đen phân bố không đều giữa chi trƣớc và chi sau; mặt trên của chi trƣớc và chi sau có các đai nhỏ màu đen; mặt trên của đuôi mới tái tạo hơi đen, không xuất hiện khoanh; cằm, họng và bụng màu kem.

Ghi chú: Mẫu cái thu đƣợc ở Cúc Phƣơng có xuất hiện lỗ trƣớc huyệt, trong khi mẫu cái từ mô tả gốc của Nguyen et al. (2015) [46], không ghi nhận đặc điểm này.

Đặc điểm sinh thái: Mẫu đƣợc thu vào lúc 6 giờ 16 tối trên mỏm đá vôi cao xấp xỉ 1m so với mặt đất, tại độ cao 205 so với mực nƣớc biển. Môi trƣờng sống là rừng nguyên sinh trên núi đá vôi. Nhiệt độ lúc thu mẫu vào khoảng 28.7 độ C với độ ẩm tƣơng đối khoảng 92%.

Phân bố: Thằn lằn chân ngón bobrovi đƣợc mô tả tại Việt Nam (Hòa Bình) (Nguyen et al. 2015). Đây là loài lần đầu tiên đƣợc ghi nhận tại tỉnh Thanh Hóa cũng nhƣ Vƣờn quốc gia Cúc Phƣơng.

Hình 4.3. Thằn lằn chân ngón bobrovi (Cyrtodactylus bobrovi)

Họ Rắn nƣớc Colubridae

2.Rắn lệch đầu hoa - Lycodon rufozonatus (Cantor, 1842) (Hình 4.4)

Mẫu phân tích: Một con cái NHQ17.18, thu ở tọa độ 20o17.385’ vĩ Bắc 105o35.078’ kinh Đông, ngày 26/8/2017 tại xã Thành Yên, huyện Thạch Thành, tỉnh Thanh Hóa, thuộc VQG Cúc Phƣơng.

Đặc điểm hình thái: Đặc điểm hình thái của mẫu vật phù hợp với mô tả của Nguyễn Văn Sáng (2007) [17]. Dài thân (SVL 980 mm); dài đuôi (TaL 220 mm); Đƣờng tiếp giáp của 2 tấm gian mũi bằng một n a đƣờng tiếp giáp của hai tấm trƣớc trán; 1 tấm trán; 2 tấm đỉnh lớn. Tấm má tiếp giáp mắt và ở dƣới tấm trƣớc mắt; 1 tấm trƣớc mắt; 2 tấm sau mắt; 2+4 tấm thái dƣơng. Môi trên 8 tấm, tấm thứ 3, 4, 5 tiếp giáp mắt. Môi dƣới 10 tấm, có 5 tấm tiếp giáp tấm sau cằm trƣớc; 2 đôi tấm sau cằm mà tấm trƣớc dài hơn so với tấm sau. Vảy thân: 19-17-15 hàng, những vảy trên lƣng hơi rõ gờ; 225 tấm bụng; tấm hậu môn nguyên; 84 tấm dƣới đuôi.

Màu sắc mẫu khi còn sống: Đầu đen, có 2 vệt trắng nhỏ ở trên tấm trán. Mặt dƣới đầu trắng nhạt. Lƣng xám đen có 52 vòng nâu đỏ rất mảnh chạy ngang thân. Bụng màu trắng đục rải rác có các chấm đen nhỏ.

Đặc điểm sinh thái: Mẫu Rắn lệch đầu hoa thu đƣợc vào khoảng 20h30- 21h00 tại trên mặt đất, chân núi đá, rừng thứ sinh, độ cao 102 m so với mực nƣớc biển. Nhiệt độ không khí lúc thu mẫu khoảng 250C, độ ẩm 92%.

Phân bố:

Việt Nam: Lai Châu, Lào Cai, Hà Giang, Cao Bằng, Lạng Sơn, Ninh Bình (Tràng An), từ Nghệ An đến Thừa Thiên-Huế (Luu et al. 2015) [38]. Ghi nhận phân bố mới của loài cho tỉnh Thanh Hóa và VQG Cúc Phƣơng.

Thế giới: Trung Quốc (Luu et al. 2015) [38].

Ghi chú: 01 mẫu cái (NHQ17.18) từ VQG Cúc Phƣơng có số hàng vẩy thân là 23-21- 17 trong khi mẫu cái theo mô tả của Luu et al. 2015 là 21-21- 15.

Việt Nam: Tuyên Quang, Vĩnh Phúc, Hà Tĩnh, Quảng Bình, Quảng Trị

Hình 4.4. Rắn lệch đầu hoa (Lycodon rufozonatus)

Họ Rắn bồng Homalopsidae

4.Rắn bồng trung quốc Myrrophis chinensis (Gray, 1842) (Hình 4.5.)

Mẫu phân tích: Một con đực NHQ17.130, thu ở tọa độ 20o15.460’vĩ Bắc, 105o42.560’kinh Đông, ngày 20/11/2017 tại tiểu khu 20, xóm Mạc, Cúc Phƣơng, Nho Quan, Ninh Bình.

Đặc điểm hình thái: Đặc điểm hình thái của mẫu vật phù hợp với mô tả của Ziegler et al. (2007) [73]. Dài thân (SVL 269 mm); dài đuôi (TaL 45 mm); 1 tấm gian mũi; 2 tấm trƣớc trán; 1 tấm trán, dài gấp 2 lần rộng; 2 tấm đỉnh; tấm mũi chia một phần; 1 tấm má, không tiếp giáp với tâm gian mũi và mắt; 1 tấm trƣớc mắt dài; 2 tấm sau mắt; 1+2+3 tấm thái dƣơng. Môi trên 7 tấm, tấm thứ 4 tiếp giáp mắt. Môi

dƣới 10 tấm, 4 tấm tiếp giáp tấm sau cằm trƣớc; 2 đôi tấm sau cằm mà tấm trƣớc dài hơn tấm sau, 2 tấm sau cách nhau bởi 6 tấm nhỏ. Vảy thân: 23-23-19, nhẵn; 147 vảy bụng; tấm hậu môn chia hai; 48 tấm dƣới đuôi, kép.

Màu sắc khi mẫu còn sống: Trên đầu có vệt xám to từ đỉnh đầu tới cổ. Lƣng xám có các hàng chấm xám đen xếp xiên và nhiều về phía hai bên sƣờn. Hai bên sƣờn có hai đƣờng màu trắng vàng chạy từ đầu tới mút đuôi. Bụng màu xám nhạt có những vệt xám đen chạy ngang qua nơi tiếp giáp 2 tấm bụng.

Đặc điểm sinh thái: Mẫu Rắn bồng trung quốc thu đƣợc vào khoảng 22h00- 22h30 ở vũng nƣớc, dƣới đầm lầy, rừng thứ sinh, nơi có nƣớc quanh năm, độ cao 155m so với mực nƣớc biển. Nhiệt độ không khí lúc thu mẫu khoảng 180C, độ ẩm 70%.

Hình 4.5. Rắn bồng Trung Quốc (Myrrophis chinensis)

Phân bố:

Việt Nam: Thái Nguyên, Vĩnh Phúc, Hải Dƣơng, Bắc Giang, Sơn La, Hà Nội, Hòa Bình, Bắc Trung Bộ, Đắc Lắk, Đắc Nông (Nguyen et al. 2009) [42]. Đây là ghi nhận phân bố mới của loài này ở tỉnh Ninh Bình và VQG Cúc Phƣơng.

Thế giới: Trung Quốc, Đài Loan (Nguyen et al. 2009) [42]

Họ Rắn lục Viperidae

5. Rắn lục sừng -Protobothrops cornutus (Smith, 1930) (Hình 4.6)

Mẫu phân tích: Hai con cái, mẫu NHQ17.18, NHQ17.59 thu ở tọa độ 20o17.531’ vĩ Bắc, 105o34.696’ kinh Đông, ngày 26/8/2017 tại khu vực Thành Yên, Thanh Hóa và mẫu NHQ17.59, thu ở tọa độ 20o17.738’ vĩ Bắc, 105o40.026’ kinh Đông, ngày 01/9/2017 tại động Ngƣời xƣa, thuộc VQG Cúc Phƣơng.

Đặc điểm hình thái: Đặc điểm hình thái của mẫu vật phù hợp với mô tả của Luu et al. (2015). Dài thân (SVL 550-580 mm); dài đuôi (TaL 135-142 mm); dài đầu (HL 22,2-26); rộng đầu (HW 19-19,4). Những vẩy trên đầu nhỏ, nhẵn ở đỉnh đầu và có gờ ở phần sau đầu. Có 13-14 vảy nhỏ giữa 2 tấm trên mắt. Tấm trên mắt chia ra thành 3 hay 4 tấm nhỏ và kéo dài làm thành phần phụ giống cái sừng trên mắt; 4 tấm sau mắt; 1 tấm dƣới mắt dài. Môi trên 10-11 tấm, tấm thứ 1 và tấm mũi tách biệt, tấm thứ 3 lớn hơn, tấm thứ 4 cách tấm dƣới mắt bởi 3 hàng vảy nhỏ. Môi dƣới 10-12 tấm, có1 hoặc 3 tấm tiếp giáp tấm sau cằm trƣớc. Vảy thân: 21-21-15 ở mẫu NHQ17.59 và 23-21-17 ở mẫu NHQ17.18, có gờ; tấm bụng 194-195; tấm hậu môn đơn; 69-77 tấm dƣới đuôi kép.

Màu sắc mẫu khi còn sống: Đầu có những vết nâu nhỏ. Lƣng nâu xám có 2 hàng đốm sẫm hình tứ giác viền đen chạy dọc lƣng, nối lại thành những vạch ngang chính giữa lƣng. Một dãy vết hơi trắng ở mỗi bên sƣờn. Bụng hơi trắng pha nâu.

Đặc điểm sinh thái: Mẫu Rắn lục sừng thu đƣợc vào khoảng 19h00-20h00 chân đồi đất của động Ngƣời xƣa, độ cao 298 m so với mực nƣớc biển. Nhiệt độ không khí lúc thu mẫu khoảng 260C, độ ẩm 95%.

Hình 4.6. Rắn lục sừng (Protobothrops cornutus)

4.1.4. Khả năng phát hiện loài các loài bò sát

Qua hình 4.7 cho thấy, số loài phát hiện tăng nhanh từ ngày thứ nhất cho đến ngày thứ 11, sau đó tăng chậm đến ngày 23; từ ngày 24 đến ngày 40 số loài tăng ít và có dấu hiệu ổn định theo thời gian. Điều này cho thấy số lƣợng loài có thể dao động trong khoảng này. Danh lục bò sát cập nhật so với báo cáo cũ là 80 loài,

trong đó có 44 loài trong báo cáo cũ không ghi nhận ở nghiên cứu này. Điều này cho thấy, môi trƣờng sống của một số loài đã bị tác động, sinh cảnh đang dần thay đổi hoặc cũng có thể nhiều loài là đối tƣợng bị săn bắt quá mức nên số lƣợng bị suy giảm qua nhiều năm. Tuy nhiên do thời gian điều tra ngắn và chỉ vào một số thời điểm nhất định trong năm nên chƣa thể đánh giá hết đƣợc mức độ đa dạng thành phần loài. Vì vậy, nếu tiếp tục nỗ lực điều tra thì có khả năng ghi nhận thêm loài.

Hình 4.7. Đƣờng cong phát hiện loài bò sát tại VQG Cúc Phƣơng

4.1.5. Đa dạng về sinh cảnh sống và phân bố của bò sát tại VQG Cúc Phương

4.1.5.1. Đa dạng các loài bò sát theo sinh cảnh

Thống kê số lƣợng cá thể của mỗi loài trong tổng số 36 loài ghi nhận ngoài thực địa theo 5 dạng sinh cảnh (bảng 4.2).

Bảng 4.2. Phân bố của bò sát theo các sinh cảnh tại VQG Cúc Phƣơng

TT Tên loài Sinh cảnh

1 2 3 4 5

1 Ba ba trơn (Pelodiscus sinensis) x

2 Liu điu chỉ (Takydromus sexlineatus) x x x x

3 Nhông em-ma (Calotes emma) x x x

4 Nhông xanh (Calotes versicolor) x x x

5 Ô rô vẩy (Acanthosaura lepidogaster) x x x

6 Rắn bông chì (Hypsiscopus plumbea) x

7 Rắn bồng trung quốc (Myrrophis chinensis) x

8 Rắn cạp nong (Bungarus fasciatus) x x

9 Rắn cƣờm (Chrysopelea ornate) x x

10 Rắn hổ mây ham tơn (Pareas hamptoni) x x x

11 Rắn hoa cỏ nhỏ (Rhabdophis subminiatus) x x x

12 Rắn khiếm trung quốc (Oligodon chinensis) x x

13 Rắn khuyết cạp nong (Lycodon fasciatus) x x

14 Rắn lệch đầu kinh tuyến (Lycodon meriodinalis) x x

15 Rắn lệch đầuhoa (Lycodon rufozonatus) x

16 Rắn lục cƣờm (Protobothrops mucrosquamatus) x x

17 Rắn lục mép trắng (Trimeresurus albolabris) x x

18 Rắn lục sừng (Protobothrops cornutus) x x x

19 Rắn lục xanh (Trimeresurus stejnegeri) x x x

20 Rắn mống (Xenopeltis unicolor) x

21 Rắn nhiều đai (Cyclophiops multicinctus) x x x

22 Rắn rào đốm (Boiga multomaculata) x x

23 Rắn ráo thƣờng (Ptyas korros) x x x x

24 Rắn roi thƣờng (Ahaetulla prasina) x x x

25 Rắn sọc khoanh (Elaphe moellendorffii) x x x

26 Tắc kè (Gekko reevesii) x

TT Tên loài Sinh cảnh

1 2 3 4 5

28 Thạch sùng đuôi sần (Hemidactylus frenatus) x

29 Thằn lằn bay đốm (Draco maculatus) x x

30 Thằn lằn phê nô ba vạch (Sphenomorphus tritaeniatus) x x

31 Thằn lằn bóng đuôi dài (Eutropis longicaudata) x x

32 Thằn lằn chân ngón bobrovi (Cyrtodactylus bobrovi) x

33

Thằn lằn chân ngón cúc phƣơng (Cyrtodactylus cucphuongensis)

x

x

34 Thằn lằn tai ba vì (Tropidophorus baviensis) x

35 Rắn bồng chì (Hypsiscopus plumbea) x x

36 Rắn sãi thƣờng (Amphiesma stolatum) x x

Tổng số 16 20 25 3 11

Ghi chú: 1: Rừng trên núi đá vôi (SC1); 2: Rừng trên núi đất (SC2); 3: Rừng trên núi đá lẫn núi đất (SC3); 4: Sông suối, ao hồ, đồng ruộng (SC4); 5: Nương rẫy, trảng cỏ, cây bụi (SC5).

Hình 4.8. Đa dạng thành phần loài bò sát theo sinh cảnh

Qua bảng 4.2. và hình 4.8. cho thấy các loài bò sát ở VQG Cúc Phƣơng phân bố chủ yếu ở 3 dạng sinh cảnh chính đó là SC1; SC2; SC3 chiếm tới 81%. SC4 chỉ ghi nhận 3 loài chiếm 4%.

4.1.5.2. Đa dạng các loài bò sát theo KVNC

Qua hình 4.9 cho thấy thành phần loài giữa ba tỉnh chênh lệch khá lớn, cụ thể tỉnh Thanh Hóa chiếm diện tích nhỏ nhất (chiến 22,5% tổng diện tích Vƣờn) nhƣng thành phần loài đa dạng nhất với 74 loài; trong khi đó tỉnh Ninh Bình chiếm diện tích hơn 50 % nhƣng chỉ ghi nhận 68 loài; tỉnh Hòa Bình (chiếm 26,4 % tổng diện tích Vƣờn) ghi nhận đƣợc 56 loài. Tỉnh Thanh Hóa chiếm diện tích nhỏ nhƣng thành phần loài đa dạng có thể do, khu vực này giáp ranh với KBT Ngọc Sơn- Ngổ Luông (đặc trƣng rừng núi đá) nên có sự giao thoa giữa các loài động vật, ngoài ra khu vực này dân cƣ thƣa thớt và sống xa rừng hơn so với hai tỉnh còn lại.

Hình 4.9. So sánh thành phần loài bò sát giữa 3 tỉnh

4.1.6. So sánh sự tương đồng về thành phần loài bò sát của VQG Cúc Phương và các khu bảo tồn khác có sinh cảnh tương đồng

Để góp phần đánh giá tính đa dạng khu hệ bò sát tại VQG Cúc Phƣơng đề tài tiến hành so sánh thành phần loài bò sát của khu vực với một số Khu bảo tồn và Vƣờn Quốc gia khác. Đây là các khu vực tƣơng đối đồng nhất về địa hình, sinh thái với khu vực nghiên cứu. Kết quả thể hiện trong bảng 4.3 và hình 4.11.

Bảng 4.3. So sánh số lƣợng taxon bò sát tại VQG Cúc Phƣơng với các khu bảo tồn khác có cùng hệ sinh thái núi đá vôi

TT Khu vực Số bộ Số họ Số loài Nguồn

1 VQG Cúc

Phƣơng 2 19 80

Tổng hợp kết quả từ nghiên cứu này 2 VQG Ba Bể 2 16 33 Nguyen et al. (2009) 3 KBT Kim Hỷ 2 8 27 Phạm Thị Kim Dung

(2015)

4 VQG Cát Bà 2 17 45 Nguyen et al. (2011) 5 VQG Phong

Nha- Kẻ Bàng 2 18 101 Luu et al. (2013) 6 KBT Du Già 2 9 17 Bain & Nguyen (2004) 7 Khu danh thắng Tràng An 2 12 34 Trần Thị Hồng Ngọc và cs (2017) 8 KBT ĐNN Vân Long 2 8 13 Lê Trung Dũng và cs (2016)

Hình 4.10. Mức độ đa dạng các họ bò sát giữa các khu bảo tồn khác có cùng hệ sinh thái núi đá vôi

Qua hình trên cho thấy VQG Cúc Phƣơng có sự đa dạng về họ cao nhất so với các VQG và KBT còn lại nhƣ VQG Cát Bà (17 họ); VQG Ba Bể (16 họ); KBT Vân Long và KBT Kim Hỷ ít đa dạng nhất với 8 họ.

Hình 4.11. Mức độ đa dạng thành phần loài bò sát giữa các khu bảo tồn có cùng hệ sinh thái núi đá vôi

Qua hình 4.11 cho thấy VQG Cúc Phƣơng có sự đa dạng về thành phần loài rất cao chỉ đứng sau VQG Phong Nha – Kẻ Bàng

So sánh chi tiết về thành phần loài để tính toán hệ số tƣơng tự giữa 08 khu vực (bảng 4.4. và hình 4.12). Kết quả phân tích thống kê s dụng phần mềm PAST Statistic (Hammer et al. 2001) cho thấy mức độ tƣơng đồng về thành phần loài cao nhất giữa VQG Cúc Phƣơng với VQG Phong Nha – Kẻ Bàng (djk = 0,47); giữa VQG Cúc Phƣơng và Quần thể danh thắng Tràng An (djk=0,46) và thấp nhất giữa VQG Cúc Phƣơng và KBT ĐNN Vân Long (djk= 0,20).

Bảng 4.4. Hệ số tƣơng tự về thành phần loài bò sát giữa các khu bảo tồn có cùng hệ sinh thái núi đá vôi

0 CP VL TA BB CP DG KH PN-KB CP 1 VL 0.2 1 TA 0.46018 0.22222 1 BB 0.39286 0.090909 0.50746 1 CB 0.45161 0.25 0.4557 0.51282 1 DG 0.20833 0.071429 0.31373 0.36 0.35484 1 KH 0.28302 0.21053 0.45902 0.46667 0.55556 0.40909 1 PN- KB 0.47514 0.15929 0.39706 0.35556 0.39456 0.18487 0.34109 1

Ghi chú: CP: Cúc Phương; VL: Vân Long; TA: Tràng An; CB: Cát Bà; DG: Du Già; KH: Kim Hỷ; PN- KB:

Một phần của tài liệu (LUẬN văn THẠC sĩ) đa dạng các loài bò sát (reptilia) và ếch nhái (amphibia) tại vườn quốc gia cúc phương, tỉnh ninh bình​ (Trang 45)

Tải bản đầy đủ (PDF)

(90 trang)