Tổng quan về luận án
Bảo tồn đa dạng sinh học (ĐDSH) và duy trì tính toàn vẹn của các dịch vụ hệ sinh thái rừng đang là trọng tâm cấp bách trong bối cảnh suy thoái tài nguyên toàn cầu và biến đổi khí hậu diễn biến phức tạp. Khu Bảo tồn Thiên nhiên (BTTN) Đồng Sơn – Kỳ Thượng (tỉnh Quảng Ninh) với diện tích 15.593,81 ha thuộc địa bàn 5 xã (Đồng Sơn, Kỳ Thượng, Đồng Lâm, Vũ Oai, Hòa Bình) đại diện cho hệ sinh thái rừng kín thường xanh núi thấp tập trung lớn nhất vùng Đông Bắc Việt Nam. Mặc dù sở hữu vị trí địa sinh học xung yếu với độ cao biến thiên từ chân đồi lên đỉnh Thiên Sơn (1.096 m) và Am Váp (1.051 m), tài nguyên thực vật tại đây phải đối mặt với áp lực nhân sinh gay gắt từ 1.504 nhân khẩu vùng đệm và vùng lõi (trong đó đồng bào dân tộc Dao chiếm tới 79,5% toàn vùng và 97% trong vùng lõi), dẫn đến nguy cơ chia cắt sinh cảnh và suy giảm nguồn gen nguy cấp.
Khoảng trống nghiên cứu (research gap) lớn trước khi luận án được thực hiện là sự thiếu hụt các đánh giá định lượng có hệ thống về cấu trúc thảm thực vật, biến động chỉ số đa dạng sinh thái theo đai cao và trạng thái diễn thế, cũng như chưa có một danh lục thực vật bậc cao có mạch cập nhật toàn diện theo các hệ thống phân loại hiện đại. Các công trình trước đây của Phân viện Điều tra Quy hoạch Rừng Đông Bắc (2010, 2011) hay Quỹ Bảo tồn Việt Nam - VCF (2011) chỉ dừng lại ở các cuộc điều tra nhanh, rời rạc giữa thực vật thân gỗ (546 loài) và thân thảo (617 loài), chưa tích hợp đa tầng giữa cấu trúc sinh thái phát sinh với phân tích rủi ro nhân tạo.
Để giải quyết triệt để khoảng trống này, luận án xác định hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:
- Câu hỏi nghiên cứu 1 (RQ1): Cấu trúc phân tầng, tổ thành loài và các chỉ số đa dạng sinh học tầng cây gỗ biến đổi như thế nào giữa các kiểu thảm thực vật nhiệt đới và á nhiệt đới núi thấp tại Khu BTTN Đồng Sơn – Kỳ Thượng?
- Câu hỏi nghiên cứu 2 (RQ2): Hệ thực vật bậc cao có mạch tại khu vực có quy mô taxon, dạng sống sinh học (life-form) và phổ công dụng như thế nào khi đối chiếu với hệ thực vật chung của Việt Nam và các vùng lân cận?
- Câu hỏi nghiên cứu 3 (RQ3): Những nhân tố tự nhiên (đai cao, đất Feralit) và nhân sinh (khai thác gỗ, lâm sản ngoài gỗ, chăn thả) tác động như thế nào đến tính đa dạng và nguy cơ tuyệt chủng của các loài thực vật quý hiếm?
- Giả thuyết nghiên cứu 1 (H1): Chỉ số đa dạng loài ($H'$) và độ giàu loài giảm dần theo đai cao và chịu tác động tiêu cực rõ rệt tại các trạng thái rừng thứ sinh phục hồi sau nương rẫy so với rừng ít bị tác động.
- Giả thuyết nghiên cứu 2 (H2): Hệ thực vật Khu BTTN Đồng Sơn – Kỳ Thượng mang đặc trưng điển hình của hệ thực vật nhiệt đới gió mùa Đông Bắc, có tỷ lệ chồi trên (Phanerophytes) chiếm ưu thế áp đảo và sở hữu các taxon mới chưa từng được ghi nhận trong hệ thực vật Việt Nam.
Khung lý thuyết của công trình được xây dựng dựa trên sự tích hợp giữa học thuyết Sinh vật địa lý quần thể (Sukachev, 1944), hệ sinh thái của Tansley (1935), quan điểm sinh thái phát sinh quần thể thảm thực vật rừng nhiệt đới của Thái Văn Trừng (1978, 1999), kết hợp phân loại ngoại mạo UNESCO (1973) và hệ thống phổ dạng sống Raunkiaer (1934). Đóng góp đột phá được định lượng của luận án bao gồm việc hoàn thiện danh lục 1.246 loài thực vật bậc cao có mạch thuộc 688 chi, 180 họ của 5 ngành thực vật; phát hiện và bổ sung 01 loài mới cho hệ thực vật Việt Nam là Dó vân nam (Aquilaria yunnanensis S. Huang, họ Thymelaeaceae); đồng thời bổ sung mới 218 loài, 71 chi và 12 họ cho danh lục Khu BTTN Đồng Sơn – Kỳ Thượng. Công trình được tiến hành công phu trong khung thời gian từ tháng 8/2015 đến tháng 12/2019 trên toàn bộ phạm vi 15.593,81 ha của khu bảo tồn.
Literature Review và Positioning
Nghiên cứu thảm thực vật trên thế giới khởi đầu từ các mô tả ngoại mạo của Patsoxki (1915), Rubel (1930) và Ellenberg (1967). Trường phái châu Âu lấy thành phần loài làm tiêu chuẩn hạt nhân (Braun-Blanquet, 1928), trong khi trường phái Bắc Mỹ định hình dựa trên học thuyết Cực đỉnh (Climax Theory) của Clements gắn liền với tính ổn định khí hậu lâu dài. Tại khu vực nhiệt đới, Schimper (1898), Champion (1936), Aubréville (1956) và UNESCO (1973) đã chuẩn hóa khung phân loại dựa trên ngoại mạo, cấu trúc và độ tàn che của tầng ưu thế sinh thái.
Trong tài liệu học thuật tồn tại cuộc tranh luận kéo dài giữa hai luồng quan điểm:
- Quan điểm tiền định cực đỉnh (Clementsian Climax Paradigm) cho rằng thảm thực vật là một siêu cơ thể (super-organism) tiến tới trạng thái cân bằng xác định bởi khí hậu.
- Quan điểm phân bố cá thể/gradient sinh thái (Individualistic Hypothesis của Gleason và Gradient Analysis của Whittaker, 1975) khẳng định sự phân bố loài là liên tục dọc theo các gradient môi trường và chịu chi phối bởi các biến động ngẫu nhiên.
Tại Việt Nam, các nghiên cứu kinh điển của Chevalier (1918), Maurand (1953), Loschau (1960), Trần Ngũ Phương (1970) và đặc biệt là công trình "Những hệ sinh thái rừng nhiệt đới ở Việt Nam" của Thái Văn Trừng (1978, 1999) đã đặt nền móng sinh thái phát sinh quần thể vững chắc. Thái Văn Trừng đã phân lập 14 kiểu thảm thực vật dựa trên 5 nhóm nhân tố sinh thái: địa lý - địa hình, khí hậu - thủy văn, đá mẹ - thổ nhưỡng, khu hệ thực vật và nhân sinh. Về phân loại học và kiểm kê hệ thực vật, các bộ Thực vật chí Đông Dương (Lecomte, 1907–1951), Thực vật chí Campuchia, Lào và Việt Nam (1960–2018), Cây cỏ Việt Nam (Phạm Hoàng Hộ, 1999–2003) và Danh lục các loài thực vật Việt Nam (2001–2005) đã ghi nhận khoảng 11.373 đến 13.766 loài có mạch.
Định vị nghiên cứu của luận án nằm tại điểm giao thoa giữa sinh thái học thực vật định lượng và quản lý tài nguyên rừng bền vững. Trong khi các nghiên cứu quốc tế như công trình của Maxwell & Elliott tại Vườn Quốc gia Doi Suthep-Pui (Thái Lan, diện tích 261 km², ghi nhận 2.084 loài) và nghiên cứu của Zhu Hua et al. tại rừng nhiệt đới gió mùa Tây Song Bản Nạp (Vân Nam, Trung Quốc) tập trung vào phân tích mối quan hệ giữa gradient độ cao và tính đặc hữu, luận án này áp dụng phương pháp tiếp cận đa chiều để so sánh cấu trúc thảm thực vật và hệ thực vật vùng núi Đông Bắc Việt Nam với các hệ sinh thái tương đồng trong khu vực Đông Dương và Nam Trung Quốc. Luận án giải quyết triệt để sự thiếu nhất quán giữa các bảng phân loại trạng thái lâm nghiệp thuần túy (như phân loại Loschau I–IV) với phân loại cấu trúc sinh thái chuyên sâu, đưa ra luận cứ xác đáng để giải thích cơ chế diễn thế thứ sinh phục hồi sau khai thác và nương rẫy.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án mở rộng và làm giàu thêm học thuyết Sinh vật địa lý quần thể (Sukachev) và lý thuyết Sinh thái phát sinh quần thể (Thái Văn Trừng) thông qua việc lượng hóa cấu trúc tổ thành và mức độ đa dạng sinh học theo lát cắt đai cao tại vùng Đông Bắc Việt Nam. Kết quả nghiên cứu xác nhận sự tồn tại của ranh giới chuyển tiếp sinh thái (ecotone) rõ nét ở độ cao 700 m trên nền đất Feralit mùn phát triển trên sa thạch khối, phân tách hai kiểu thảm thực vật chính: Rừng kín thường xanh mưa ẩm nhiệt đới (dưới 700 m) và Rừng kín thường xanh mưa ẩm á nhiệt đới núi thấp (trên 700 m).
Về mặt lý thuyết tiến hóa và thích nghi, luận án kiểm chứng mô hình phổ dạng sống chuẩn ($SN$) của Raunkiaer (1934):
$$SN = 46\text{ Ph} + 9\text{ Ch} + 26\text{ Hm} + 6\text{ Cr} + 13\text{ Th}$$
Kết quả phổ sinh thái ($SB$) thực tế tại Đồng Sơn – Kỳ Thượng cho thấy tỷ lệ nhóm chồi trên ($\text{Ph}$) vượt trội đạt xấp xỉ 80%, khẳng định điều kiện nhiệt ẩm nhiệt đới gió mùa tạo ưu thế tuyệt đối cho các dạng thân gỗ và dây leo hóa gỗ, trong khi nhóm chồi ẩn ($\text{Cr}$) và chồi một năm ($\text{Th}$) chiếm tỷ lệ rất thấp. Phát hiện loài Aquilaria yunnanensis S. Huang tại khu vực mở rộng biên giới phân bố địa lý sinh học của chi Aquilaria từ cao nguyên Vân Nam xuống vùng đồi núi thấp Đông Bắc Việt Nam, cung cấp bằng chứng thực nghiệm quan trọng cho lý thuyết di lưu thực vật (floristic migration) giữa khu hệ Hoa Nam và Đông Dương.
Khung phân tích độc đáo
Khung phân tích của luận án là sự hợp nhất của ba trụ cột lý thuyết:
- Lý thuyết Sinh thái phát sinh quần thể: Sử dụng hệ thống 5 nhóm nhân tố môi trường của Thái Văn Trừng để phân loại và định danh chính xác các kiểu thảm tự nhiên và kiểu phụ thứ sinh nhân tác.
- Lý thuyết Đa dạng sinh thái Whittaker & Phổ Entropy đa bậc Renyi: Kết hợp các chỉ số alpha truyền thống ($H'$, $Cd$, $SI$) với đồ thị phổ đa dạng Renyi ($H_\alpha$), cho phép so sánh độ giàu loài và tính đồng đều ở các bậc trọng số khác nhau mà không bị thiên lệch bởi kích thước mẫu.
- Lý thuyết Dạng sống và Thích nghi sinh thái học Raunkiaer: Phân tích cấu trúc thích nghi của toàn bộ hệ thực vật với biến động thời tiết cực đoan (mùa đông khô lạnh, biên độ nhiệt ngày đêm 5–8°C, nhiệt độ cực hạ xuống 0°C).
graph TD
A[Khung phân tích đa dạng thực vật] --> B[Cấu trúc thảm thực vật]
A --> C[Đặc điểm hệ thực vật Taxon]
A --> D[Động lực nhân sinh & Bảo tồn]
B --> B1[Phân loại sinh thái phát sinh Thái Văn Trừng 1978/1999]
B --> B2[Định lượng đa dạng: IVI, Shannon H', Simpson Cd, Renyi Profile]
B --> B3[Gradient đai cao: <700m Nhiệt đới vs >700m Á nhiệt đới]
C --> C1[Kiểm kê 1.246 loài, 688 chi, 180 họ thuộc 5 ngành]
C --> C2[Phổ dạng sống Raunkiaer: Xác lập công thức SB]
C --> C3[Phân loại giá trị sử dụng 428 cây thuốc & Danh lục đỏ SĐVN/IUCN]
D --> D1[Áp lực sinh kế: 1.504 dân cư, 79.5% dân tộc Dao]
D --> D2[Tác động: Khai thác gỗ, LSNG, Nương rẫy, Chăn thả]
D --> D3[Giải pháp bảo tồn: Quản lý, Sinh kế, Kỹ thuật lâm sinh]
Điều kiện biên (boundary conditions) của khung phân tích được xác lập chặt chẽ: Áp dụng cho hệ sinh thái rừng nhiệt đới và á nhiệt đới trên đất núi đất Feralit biến tính chịu tác động của khí hậu gió mùa ven biển Đông Bắc Việt Nam, có độ cao địa hình không vượt quá 1.100 m.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu tiếp cận theo hệ hình thực chứng sau hiện thực (post-positivism) kết hợp hiện thực phê phán (critical realism), phối hợp chặt chẽ giữa điều tra định lượng thực địa với phân tích kinh tế - xã hội định tính. Thiết kế nghiên cứu đa tầng (multi-level design) được tổ chức theo:
- Tầng vĩ mô (Landscape level): Khảo sát trên toàn bộ 15.593,81 ha qua 5 xã, tích hợp phân tích bản đồ địa hình, bản đồ thổ nhưỡng (22 dạng đất thuộc 4 nhóm chính) và hệ thống thủy văn lưu vực sông Ba Chẽ và sông Man.
- Tầng trung gian (Community/Formation level): Phân định các đai độ cao sinh thái (<700 m và >700 m) và các trạng thái rừng (rừng tự nhiên nghèo 3.788,29 ha, rừng trung bình 2.028,87 ha, rừng chưa có trữ lượng 5.670,15 ha, rừng hỗn giao gỗ - tre nứa 2.264,76 ha).
- Tầng vi mô (Plot/Sample level): Thiết lập mạng lưới ô tiêu chuẩn định lượng trên các tuyến điều tra thực địa đại diện.
Quy trình nghiên cứu thực nghiệm
- Thiết lập tuyến và ô mẫu: Thiết lập hệ thống tuyến điều tra thực địa cắt qua các sinh cảnh đặc trưng từ độ cao 80 m đến 1.096 m. Bố trí các ô tiêu chuẩn (OTC) định lượng diện tích $500\text{ m}^2$ ($20\text{ m} \times 25\text{ m}$) hoặc $1.000\text{ m}^2$ ($20\text{ m} \times 50\text{ m}$) để điều tra tầng cây gỗ lớn có đường kính ngang ngực ($D_{1.3}$) $\ge 6\text{ cm}$ và chiều cao vút ngọn ($H_{vn}$).
- Điều tra tầng tái sinh và thảm tươi: Trong mỗi OTC, bố trí 5 ô dạng bản (ODB) diện tích $4\text{ m}^2$ ($2\text{ m} \times 2\text{ m}$) tại 4 góc và tâm ô để kiểm kê cây tái sinh ($H_{vn} < 2\text{ m}$), thảm tươi và cây bụi.
- Thu thập và định danh mẫu tiêu bản: Thu thập mẫu vật thực vật có đầy đủ cơ quan sinh dưỡng và sinh sản (hoa, quả). Quy trình xử lý mẫu, sấy khô, tẩm cồn và bảo quản tuân thủ nghiêm ngặt tiêu chuẩn quốc tế. Giám định tên khoa học dựa trên phương pháp hình thái học so sánh đối chiếu với các bộ Thực vật chí chuẩn (Flora of China, Thực vật chí Việt Nam, Cây cỏ Việt Nam) và tra cứu chuyên gia tại Viện Sinh thái và Tài nguyên sinh vật, Trường Đại học Lâm nghiệp.
- Điều tra xã hội học sinh thái: Tiến hành phỏng vấn sâu và khảo sát bằng bảng hỏi bán cấu trúc tại 22 thôn bản thuộc 5 xã để thu thập dữ liệu về tri thức bản địa sử dụng cây thuốc, cây lấy gỗ và lượng hóa số vụ vi phạm lâm luật giai đoạn 2015–2018.
Dữ liệu và các chỉ số thống kê toán học
Dữ liệu điều tra được chuẩn hóa và tính toán thông qua các phần mềm chuyên dụng như Bio-DAP, R (chương trình vegan) và Excel:
- Chỉ số giá trị quan trọng (Importance Value Index - IVI):
$$IVI% = \frac{N% + G% + F%}{3}$$
Trong đó: $N%$ là mật độ tương đối; $G%$ là độ tàn che/tiết diện ngang tương đối; $F%$ là tần suất xuất hiện tương đối.
- Chỉ số đa dạng loài Shannon – Wiener ($H'$):
$$H' = -\sum_{i=1}^{S} P_i \ln(P_i)$$
Trong đó $P_i = n_i / N$ là tỷ lệ cá thể loài $i$ trên tổng số cá thể của quần xã.
- Chỉ số mức độ ưu thế Simpson ($Cd$):
$$Cd = \sum_{i=1}^{S} P_i^2$$
- Chỉ số tương đồng Sorensen ($SI$):
$$SI = \frac{2C}{A + B}$$
Trong đó $A, B$ là số loài của hai quần xã đem so sánh; $C$ là số loài chung xuất hiện ở cả hai quần xã.
- Phổ đa dạng Entropy Renyi ($H_\alpha$):
$$H_\alpha = \frac{1}{1-\alpha} \ln\left(\sum_{i=1}^{S} P_i^\alpha\right) \quad (\text{với } \alpha \ge 0, \alpha \neq 1)$$
Phát hiện đột phá và implications
Những phát hiện then chốt
1. Hoàn thiện danh lục thực vật bậc cao có mạch với quy mô vượt bậc:
Nghiên cứu đã ghi nhận chính thức 1.246 loài thuộc 688 chi, 180 họ của 5 ngành thực vật bậc cao có mạch tại Khu BTTN Đồng Sơn – Kỳ Thượng. Ngành Ngọc lan (Magnoliophyta) chiếm ưu thế áp đảo tuyệt đối với 1.157 loài (chiếm 92,86%), 629 chi (91,42%) và 150 họ (83,33%). Các ngành khác chiếm tỷ trọng nhỏ nhưng giữ vai trò quan trọng về mặt tiến hóa: Polypodiophyta (Dương xỉ) có 68 loài (5,46%); Pinophyta (Hạt trần) có 11 loài (0,88%); Lycopodiophyta (Thông đất) có 9 loài (0,72%) và Psilotophyta (Quyết lá thông) có 1 loài (0,08%).
pie title Cơ cấu các ngành thực vật bậc cao có mạch (%)
"Magnoliophyta (Ngọc lan)": 92.86
"Polypodiophyta (Dương xỉ)": 5.46
"Pinophyta (Hạt trần)": 0.88
"Lycopodiophyta (Thông đất)": 0.72
"Psilotophyta (Quyết lá thông)": 0.08
2. Đóng góp phát hiện mới cho hệ thực vật Việt Nam và khu vực:
Luận án đã ghi nhận và bổ sung 01 loài thực vật mới cho hệ thực vật Việt Nam là loài Dó vân nam (Aquilaria yunnanensis S. Huang) thuộc họ Trầm (Thymelaeaceae). Đồng thời, bổ sung mới 218 loài, 71 chi và 12 họ cho hệ thực vật Khu BTTN Đồng Sơn – Kỳ Thượng so với các danh lục công bố trước đây, làm thay đổi căn bản hiểu biết khoa học về mức độ giàu có tài nguyên thực vật vùng Đông Bắc.
3. Phân lập cấu trúc thảm thực vật theo quy luật đai cao và phân tích chỉ số đa dạng:
- Đai dưới 700 m: Kiểu rừng kín thường xanh mưa ẩm nhiệt đới (Rkx) có cấu trúc 3-5 tầng tán phức tạp. Chỉ số Shannon-Wiener ($H'$) của tầng cây gỗ dao động từ $2,85$ đến $3,42$; chỉ số tập trung ưu thế Simpson ($Cd$) thấp ($0,08 - 0,14$), phản ánh tính hỗn loài cao, không có loài ưu thế tuyệt đối. Các họ chiếm ưu thế bao gồm Euphorbiaceae, Lauraceae, Moraceae, Rubiaceae, Fagaceae.
- Đai trên 700 m: Kiểu rừng kín thường xanh mưa ẩm á nhiệt đới núi thấp (Rka) có cấu trúc tầng tán đơn giản hơn (2-3 tầng), ưu thế thuộc về các loài thuộc họ Giẻ (Fagaceae), Re (Lauraceae), Chè (Theaceae) và sự xuất hiện của các loài Thông tre (Podocarpus neriifolius), Dẻ tùng sọc trắng (Amentotaxus argotaenia). Chỉ số $H'$ tại đai này đạt $2,45 - 2,90$, chỉ số tương đồng $SI$ giữa hai đai độ cao chỉ đạt $0,42$, chứng minh sự phân hóa rõ rệt về thành phần loài dọc theo gradient nhiệt ẩm.
4. Đặc trưng phổ dạng sống sinh học (SB) thích nghi cao độ:
Toàn bộ 1.246 loài được phân loại theo hệ thống Raunkiaer xác lập công thức phổ dạng sống sinh học thực tế:
$$SB = 79,29% \text{ Ph} + 5,38% \text{ Ch} + 6,58% \text{ Hm} + 5,22% \text{ Cr} + 3,53% \text{ Th}$$
Trong nhóm cây chồi trên ($\text{Ph}$), cây gỗ vừa ($\text{Me}$) chiếm 28,45%, cây gỗ nhỏ ($\text{Mi}$) chiếm 24,12%, cây gỗ lớn ($\text{Mg}$) chiếm 8,35%, và cây dây leo hóa gỗ ($\text{Lp}$) chiếm 12,68%. Sự vượt trội của nhóm $\text{Ph}$ chứng minh sinh cảnh rừng nhiệt đới ẩm điển hình duy trì sinh trưởng liên tục quanh năm.
5. Kiểm kê giá trị tài nguyên và mức độ đe dọa nguồn gen quý hiếm:
Nghiên cứu xác định được 428 loài cây có giá trị làm thuốc thuộc 330 chi và 125 họ; 511 loài cây cho gỗ; 182 loài cây làm cảnh; 145 loài cây ăn được. Tổng hợp mức độ nguy cấp ghi nhận 74 loài quý hiếm, trong đó:
- 43 loài có tên trong Sách Đỏ Việt Nam (2007) (như Ba kích - Morinda officinalis, Đỉnh tùng - Cephalotaxus mannii, Kim giao - Nageia fleuryi).
- 31 loài trong Danh lục Đỏ thế giới (IUCN).
- 15 loài thuộc Nghị định 06/2019/NĐ-CP (nhóm IA và IIA).
Implications đa chiều
- Về mặt lý thuyết: Cung cấp cơ sở dữ liệu thực nghiệm chuẩn xác khẳng định vị trí trung chuyển và giao thoa của hệ thực vật Đông Bắc Việt Nam trong chuỗi hành lang sinh thái Đông Á - Đông Nam Á.
- Về phương pháp luận: Chuẩn hóa quy trình tích hợp giữa phân loại sinh thái phát sinh (Thái Văn Trừng) và đồ thị đa dạng Renyi trong đánh giá diễn thế rừng sau nương rẫy, mở ra phương pháp tiếp cận định lượng có thể tái lập cho các khu rừng đặc dụng khác tại Việt Nam.
- Về quản lý và thực tiễn: Định vị chính xác các điểm nóng bảo tồn (hotspots) tại các sườn núi Am Váp, đỉnh Thiên Sơn và các lưu vực phòng hộ đầu nguồn cấp nước cho Hồ Cao Vân (công trình cấp nước sinh hoạt trọng yếu cho thành phố Hạ Long với dung tích 5.000 $\text{m}^3$).
- Về chính sách lâm nghiệp: Cung cấp cơ sở khoa học để UBND tỉnh Quảng Ninh và Ban Quản lý Khu BTTN xây dựng phương án giao khoán bảo vệ rừng bền vững, phát triển mô hình trồng thâm canh dược liệu bản địa (Ba kích, Sa nhân tím, Trầm hương) tại vùng đệm nhằm giải tỏa sức ép khai thác tự nhiên từ 5.790 lao động nông nghiệp địa phương.
Limitations và Future Research
- Hạn chế về công nghệ phân tử sinh học: Việc định danh 1.246 loài thực vật bậc cao có mạch trong luận án chủ yếu dựa trên phương pháp hình thái học so sánh truyền thống và đối chiếu mẫu vật tiêu bản. Nghiên cứu chưa có điều kiện áp dụng kỹ thuật giải trình tự DNA (DNA Barcoding) trên quy mô lớn đối với các nhóm phân loại phức tạp (như Orchidaceae, Lauraceae) để làm sáng tỏ cây phát sinh chủng loại phân tử.
- Giới hạn thời gian quan trắc chu kỳ diễn thế sinh khối: Các số liệu về tầng cây gỗ, cây tái sinh và thảm mục được thu thập trong khung thời gian 2015–2019. Đây là lát cắt thời gian tĩnh, chưa phản ánh đầy đủ chu kỳ động lực diễn thế tự nhiên dài hạn (long-term permanent plot monitoring) kéo dài 10–20 năm dưới tác động của biến đổi khí hậu cực đoan.
- Phạm vi không gian vùng vi khí hậu: Nghiên cứu chưa đo lường định lượng liên tục các biến số vi khí hậu (nhiệt độ đất, độ ẩm không khí tầng dưới tán, cường độ bức xạ quang hợp) tại từng tiểu sinh cảnh khe suối sâu so với dông núi yên ngựa.
Chương trình nghiên cứu tương lai (Future Research Agenda):
- Hướng 1: Ứng dụng chỉ thị phân tử DNA đa đoạn (matK, rbcL, ITS) để giám định chính xác và nghiên cứu đa dạng di truyền các quần thể loài nguy cấp như Aquilaria yunnanensis, Cephalotaxus mannii.
- Hướng 2: Thiết lập ô mẫu định vị lâu dài (Permanent Forest Dynamics Plot) diện tích 1–5 ha theo tiêu chuẩn mạng lưới CTFS-ForestGEO để theo dõi tỷ lệ tử vong, sinh trưởng đường kính và khả năng hấp thụ carbon ($CO_2$) của các kiểu rừng kín thường xanh.
- Hướng 3: Ứng dụng viễn thám siêu phổ (Hyperspectral RS) kết hợp LiDAR quét từ thiết bị bay không người lái (UAV) để lập bản đồ phân bố 3D chi tiết sinh khối và phát hiện sớm các điểm suy thoái cấu trúc tán rừng.
- Hướng 4: Nghiên cứu chiết tách hoạt chất sinh học và đánh giá dược tính chuyên sâu của 428 loài cây thuốc bản địa, phục vụ thương mại hóa lâm sản ngoài gỗ gắn với bảo tồn nguồn gen tại chỗ (in-situ).
Tác động và ảnh hưởng
Tác động học thuật (Academic Impact)
Công trình đặt ra chuẩn mực nghiên cứu mới cho chuyên ngành Quản lý tài nguyên rừng tại Việt Nam bằng cách kết hợp nhuần nhuyễn giữa thực vật học cổ điển và sinh thái học định lượng. Bộ dữ liệu 1.246 loài và phát hiện loài Aquilaria yunnanensis S. Huang là nguồn trích dẫn nền tảng cho các công trình nghiên cứu về thực vật chí Đông Nam Á, sinh thái học rừng nhiệt đới và tiến hóa thực vật. Tiềm năng trích dẫn ước tính đạt từ 80–120 lượt trích dẫn khoa học trên các tạp chí Scopus/ISI chuyên ngành lâm nghiệp và đa dạng sinh học trong giai đoạn 5–10 năm tới.
Chuyển đổi ngành Lâm nghiệp và Thực tiễn quản lý
Luận án chuyển đổi phương thức quản lý từ bảo vệ thụ động, thuần túy hành chính sang quản lý sinh thái dựa trên bằng chứng khoa học (evidence-based conservation management). Kết quả điều tra chi tiết diện tích 14.279,8 ha rừng tự nhiên (trong đó rừng nghèo 3.788,29 ha, rừng trung bình 2.028,87 ha) giúp Ban Quản lý Khu BTTN quy hoạch chính xác các phân khu bảo vệ nghiêm ngặt, phân khu phục hồi sinh thái và hành lang đa dạng sinh học.
Tác động chính sách và Lợi ích xã hội
Cung cấp luận cứ vững chắc để cơ quan lập pháp địa phương (HĐND, UBND tỉnh Quảng Ninh) ban hành các chính sách chi trả dịch vụ môi trường rừng (PES), cơ chế đồng quản lý rừng đặc dụng với cộng đồng dân tộc Dao tại 5 xã vùng cao. Về mặt an ninh nguồn nước, việc bảo vệ nghiêm ngặt thảm thực vật rừng tại lưu vực hồ Cao Vân đảm bảo cung cấp nguồn nước ngọt an toàn, bền vững với dung tích 5.000 $\text{m}^3$ phục vụ sinh hoạt và phát triển công nghiệp du lịch của toàn bộ thành phố Hạ Long.
Đối tượng hưởng lợi
graph LR
User[Đối tượng hưởng lợi] --> A[Nghiên cứu sinh & Nhà khoa học]
User --> B[Ban Quản lý Rừng & Chi cục Kiểm lâm]
User --> C[Cộng đồng địa phương 5 xã]
User --> D[Nhà hoạch định chính sách Quảng Ninh]
A --> A1[Cơ sở dữ liệu 1.246 loài & Loài mới Aquilaria yunnanensis]
B --> B1[Bản đồ phân vùng sinh thái & Điểm nóng bảo tồn loài quý hiếm]
C --> C1[Sinh kế bền vững: Trồng dược liệu, chi trả DVMTR]
D --> D1[Quyết sách phân bổ ngân sách & Quy hoạch không gian bảo tồn]
- Nghiên cứu sinh và Giới học thuật: Tiếp cận cơ sở dữ liệu chuẩn hóa về 180 họ, 688 chi, 1.246 loài thực vật bậc cao có mạch; hệ thống công thức tính toán chỉ số sinh thái ($IVI, H', Cd, SI, H_\alpha$) và các khoảng trống tri thức về sinh học phân tử cần tiếp tục khai phá.
- Ban Quản lý Khu BTTN và Cán bộ Kiểm lâm: Sở hữu danh mục chi tiết 74 loài nguy cấp (Sách Đỏ Việt Nam, IUCN, NĐ 06/2019/NĐ-CP), vị trí phân bố theo đai cao và hiện trạng 22 dạng đất để xây dựng phương án tuần tra, phòng cháy chữa cháy rừng và ngăn chặn khai thác gỗ lậu có trọng điểm.
- Doanh nghiệp R&D Dược liệu và Lâm sản: Thừa hưởng danh lục 428 loài cây thuốc và 511 loài cây gỗ, trong đó có 39 loài cây thuốc có giá trị kinh tế cao, mở ra tiềm năng hợp tác phát triển vùng trồng dược liệu sạch dưới tán rừng.
- Cộng đồng dân tộc thiểu số tại địa phương: Hưởng lợi từ các dự án hỗ trợ chuyển giao kỹ thuật nhân giống cây bản địa, thay thế tập quán canh tác nương rẫy du canh bằng mô hình nông lâm kết hợp, ổn định sinh kế cho 1.475 hộ dân trong vùng.
Câu hỏi chuyên sâu
1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?
Đóng góp lý thuyết độc đáo nhất là việc mở rộng và kiểm chứng thực nghiệm học thuyết Sinh vật địa lý quần thể (Sukachev, 1944) và lý thuyết Sinh thái phát sinh thảm thực vật (Thái Văn Trừng, 1978, 1999) trên không gian sơn văn Đông Bắc Việt Nam. Luận án đã mô hình hóa sự phân cắt rõ nét của hai kiểu thảm thực vật theo ranh giới đai cao 700 m tương ứng với hai nhóm đất Feralit (đất Feralit mùn trên núi cao sa thạch và đất Feralit đỏ vàng vùng đồi). Đồng thời, việc phát hiện loài Aquilaria yunnanensis S. Huang thuộc họ Thymelaeaceae mở rộng lý thuyết về hành lang di lưu thực vật xuyên biên giới giữa khu hệ Cổ nhiệt đới (Paleotropical) và khu hệ Đông Á.
2. Tính đổi mới trong phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?
So với công trình của Maxwell & Elliott (2001) tại Vườn Quốc gia Doi Suthep-Pui (Thái Lan) chỉ thuần túy thống kê tần suất xuất hiện và độ phong phú loài theo phương pháp ô mẫu truyền thống, và nghiên cứu của Zhu Hua et al. (2006) tại Vân Nam (Trung Quốc) tập trung vào tỷ lệ đặc hữu chi thực vật, luận án này thể hiện sự đổi mới phương pháp luận vượt trội:
- Tích hợp đồng thời phương pháp phân loại cấu trúc ngoại mạo 5 tầng tán với phương pháp định lượng đa chỉ số ($IVI$, Shannon $H'$, Simpson $Cd$, Sorensen $SI$).
- Lần đầu tiên áp dụng công cụ phổ Entropy đa bậc Renyi ($H_\alpha$) để phân tích độ nhạy của tính đa dạng thực vật tầng cây gỗ trước các tác động nhân vi (nương rẫy, khai thác chọn), khắc phục hoàn toàn nhược điểm phụ thuộc vào kích thước mẫu của các chỉ số đơn lẻ.
3. Phát hiện bất ngờ nhất (counter-intuitive/surprising finding) kèm bằng chứng dữ liệu?
Phát hiện bất ngờ nhất là: Mặc dù rừng tự nhiên tại Khu BTTN Đồng Sơn – Kỳ Thượng bị tác động khai thác lâm sản và nương rẫy phân mảnh nghiêm trọng trong lịch sử (rừng nghèo và rừng chưa có trữ lượng chiếm tới 9.458,44 ha / 14.279,8 ha rừng tự nhiên, tương đương 66,2%), nhưng chỉ số đa dạng loài Shannon-Wiener ($H'$) của tầng cây gỗ tại các kiểu rừng thứ sinh phục hồi sau khai thác vẫn duy trì ở mức cao ($H' = 2,98 - 3,25$), không chênh lệch quá lớn so với các trạng thái rừng ít bị tác động ($H' = 3,42$).
Giải thích lý thuyết: Hiện tượng này phản ánh Giả thuyết Xáo trộn Trung gian (Intermediate Disturbance Hypothesis của Connell, 1978), khi các khoảng trống tán rừng (canopy gaps) sau khai thác chọn đã kích thích mạnh mẽ mật độ tái sinh của các loài cây gỗ tiên phong ưa sáng thuộc họ Euphorbiaceae, Moraceae xen lẫn các loài cây chịu bóng của rừng nguyên sinh, làm gia tăng cục bộ độ giàu loài tạm thời trước khi quần thể đạt trạng thái cạnh tranh loại trừ.
4. Luận án có cung cấp quy trình tái lập (Replication Protocol) minh bạch không?
Có. Luận án cung cấp quy trình tái lập hoàn chỉnh và chuẩn hóa:
- Mô tả tọa độ địa lý chính xác (Từ $107^\circ 00'30''$ đến $107^\circ 14'00''$ Đ; từ $21^\circ 04'00''$ đến $21^\circ 11'00''$ B) cùng sơ đồ bản đồ mạng lưới tuyến điều tra thực địa cắt qua 5 xã.
- Quy chuẩn hóa kích thước ô tiêu chuẩn OTC ($20\text{ m} \times 25\text{ m}$ hoặc $20\text{ m} \times 50\text{ m}$), tiêu chí đo đếm cây gỗ ($D_{1.3} \ge 6\text{ cm}$, đo $H_{vn}$, $H_{dc}$ bằng thước Blume-Leiss), ô dạng bản ($2\text{ m} \times 2\text{ m}$).
- Toàn bộ công thức toán học ($IVI, H', Cd, SI, H_\alpha$) và danh mục tài liệu khóa định loại thực vật chuyên khảo (Thực vật chí Việt Nam, Cây cỏ Việt Nam) được liệt kê chi tiết, cho phép các nhà nghiên cứu độc lập tái lập quy trình tại bất kỳ vùng đồi núi nào tại Việt Nam.
5. Chương trình nghiên cứu 10 năm (10-year Research Agenda) được vạch ra như thế nào?
Chương trình nghiên cứu 10 năm tiếp theo tập trung vào 3 giai đoạn mang tính hệ thống:
- Giai đoạn 1 (Năm 1–3): Xây dựng ngân hàng mã vạch DNA (DNA Barcode Library) cho toàn bộ 74 loài thực vật quý hiếm và 428 loài cây thuốc; thiết lập 02 ô định vị động thái rừng thường trực diện tích 1 ha tại đai nhiệt đới và á nhiệt đới.
- Giai đoạn 2 (Năm 4–6): Ứng dụng công nghệ ảnh vệ tinh đa thời gian Sentinel-2 kết hợp cảm biến LiDAR quét định kỳ nhằm mô hình hóa tốc độ phục hồi sinh khối và khả năng hấp thụ carbon của 5.670 ha rừng chưa có trữ lượng.
- Giai đoạn 3 (Năm 7–10): Triển khai chương trình bảo tồn chuyển chỗ (ex-situ) và nghiên cứu thuần hóa, chiết xuất hợp chất thứ cấp từ loài mới Aquilaria yunnanensis và các loài đặc hữu phục vụ công nghệ y dược phẩm quốc gia.
Kết luận
- Xác lập hệ thống danh mục thực vật bậc cao có mạch đồ sộ: Luận án đã điều tra, kiểm kê và chuẩn hóa danh lục gồm 1.246 loài thuộc 688 chi, 180 họ của 5 ngành thực vật bậc cao có mạch, khẳng định Khu BTTN Đồng Sơn – Kỳ Thượng là trung tâm đa dạng thực vật quan trọng nhất vùng Đông Bắc Việt Nam.
- Ghi nhận phát hiện mới cho hệ thực vật Việt Nam: Bổ sung chính thức 01 loài mới cho hệ thực vật Việt Nam là Dó vân nam (Aquilaria yunnanensis S. Huang, họ Thymelaeaceae), đồng thời bổ sung 218 loài, 71 chi, 12 họ mới cho danh lục khu bảo tồn.
- Lượng hóa chi tiết cấu trúc thảm thực vật theo quy luật đai cao: Phân định rõ 2 kiểu thảm thực vật chính theo ranh giới 700 m; xác định các chỉ số sinh thái tầng cây gỗ đạt mức đa dạng cao ($H' = 2,85 - 3,42$ ở đai nhiệt đới; $H' = 2,45 - 2,90$ ở đai á nhiệt đới núi thấp), phân tích thành công động lực tổ thành loài qua chỉ số giá trị quan trọng $IVI$ và phổ đa dạng Renyi.
- Chuẩn hóa đặc trưng sinh học dạng sống: Xác lập công thức phổ dạng sống sinh học thực tế $SB = 79,29% \text{ Ph} + 5,38% \text{ Ch} + 6,58% \text{ Hm} + 5,22% \text{ Cr} + 3,53% \text{ Th}$, chứng minh sự thích nghi tuyệt đối của hệ thực vật với điều kiện khí hậu nhiệt đới mưa mùa ẩm.
- Đánh giá toàn diện giá trị sử dụng và nguồn gen nguy cấp: Kiểm kê 428 loài cây thuốc, 511 loài cây lấy gỗ và nhận diện 74 loài thực vật quý hiếm nằm trong Sách Đỏ Việt Nam (43 loài), Danh lục Đỏ IUCN (31 loài) và Nghị định 06/2019/NĐ-CP (15 loài).
- Đề xuất hệ thống giải pháp bảo tồn đồng bộ: Xây dựng 4 nhóm giải pháp chiến lược khả thi gồm: Tăng cường quản lý bảo vệ rừng trên nền tảng tuần tra GIS; Nâng cao sinh kế và giao khoán bảo vệ rừng gắn với chi trả DVMTR cho 5 xã vùng đệm; Tuyên truyền giáo dục môi trường cho đồng bào dân tộc Dao; và Ứng dụng kỹ thuật lâm sinh khoanh nuôi xúc tiến tái sinh tự nhiên trên 5.670,15 ha rừng nghèo kiệt.
Công trình tạo bước tiến quan trọng trong hệ hình quản lý tài nguyên rừng nhiệt đới, mở ra ba nhánh nghiên cứu đột phá mới: Sinh thái học phân tử và mã vạch DNA thực vật vùng Đông Bắc; Động thái sinh khối rừng trong bối cảnh biến đổi khí hậu; và Thương mại hóa bền vững nguồn tài nguyên cây thuốc bản địa gắn liền với bảo tồn đa dạng sinh học cấp quốc gia và quốc tế.