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 hệ sinh thái rừng nhiệt đới đóng vai trò quyết định trong việc ổn định chu trình sinh địa hóa toàn cầu và ứng phó với biến đổi khí hậu. Khu Bảo tồn Thiên nhiên (BTTN) Đồng Sơn – Kỳ Thượng (tỉnh Quảng Ninh) là đại diện tiêu biểu cho hệ sinh thái rừng kín thường xanh mưa ẩm nhiệt đới và á nhiệt đới núi thấp có diện tích rừng tự nhiên tập trung lớn nhất vùng Đông Bắc Việt Nam (14.279,8 ha trên tổng số 15.593,81 ha diện tích tự nhiên). Mặc dù sở hữu vị trí địa sinh học xung yếu, nguồn tài nguyên thực vật tại đây đang chịu áp lực nhân sinh gay gắt từ các cộng đồng dân cư vùng đệm. Luận án tiến sĩ lâm nghiệp của nghiên cứu sinh Dương Trung Hiếu (chuyên ngành Quản lý Tài nguyên Rừng, mã số 9 62 02 11), dưới sự hướng dẫn khoa học của GS. Hoàng Văn Sâm tại Trường Đại học Lâm nghiệp Việt Nam, là công trình tiên phong đánh giá toàn diện, định lượng và có hệ thống về cấu trúc thảm thực vật, tính đa dạng hệ thực vật bậc cao có mạch và các yếu tố tương tác kinh tế - xã hội tại khu bảo tồn này.

Nghiên cứu giải quyết khoảng trống học thuật (research gap) quan trọng: các công trình trước năm 2015 tại khu vực chỉ dừng lại ở các cuộc điều tra nhanh mang tính kiểm kê đơn lẻ (như báo cáo của VCF 2011, phân viện ĐTQH Đông Bắc 2010), thiếu khung phân loại thảm chuẩn hóa, vắng bóng các phân tích định lượng về chỉ số đa dạng tầng cây gỗ theo đai cao và chưa có mô hình tích hợp giữa phổ dạng sống, phân bố taxon với áp lực sinh kế bản địa.

Hệ thống câu hỏi nghiên cứu và giả thuyết được xác lập cụ thể:

  1. Câu hỏi nghiên cứu 1 (RQ1): Cấu trúc tổ thành, ngoại mạo và các chỉ số đa dạng sinh học ($\alpha, \beta$-diversity) của tầng cây gỗ biến đổi như thế nào giữa các kiểu thảm thực vật và theo đai độ cao tại Khu BTTN Đồng Sơn – Kỳ Thượng?
    • Giả thuyết 1 (H1): Chỉ số đa dạng loài ($H'$) và độ ưu thế ($Cd$) có sự phân hóa rõ rệt giữa đai nhiệt đới (<700m) và á nhiệt đới (>700m), trong đó rừng phục hồi sau khai thác duy trì mức độ đa dạng entropy cao hơn rừng sau nương rẫy.
  2. 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ó cấu trúc phân loại học, phổ dạng sống và thành phần loài nguy cấp như thế nào khi đối chiếu với hệ thống phân loại hiện đại?
    • Giả thuyết 2 (H2): Hệ thực vật Đồng Sơn – Kỳ Thượng mang tính chất chuyển tiếp độc đáo giữa luồng thực vật bản địa Bắc Việt Nam - Nam Trung Hoa và chứa đựng các taxon mới chưa từng được ghi nhận cho hệ thực vật Việt Nam.
  3. Câu hỏi nghiên cứu 3 (RQ3): Những nhân tố tự nhiên và nhân sinh nào chi phối trực tiếp và gián tiếp đến tính đa dạng thực vật, và giải pháp quản lý nào đảm bảo tính bền vững?
    • Giả thuyết 3 (H3): Áp lực từ thói quen canh tác nương rẫy, khai thác lâm sản ngoài gỗ và tình trạng thiếu hụt lương thực (210 kg/người/năm) của cộng đồng người Dao bản địa là động lực chính gây suy thoái tài nguyên rừng.

Khung lý thuyết của luận án được xây dựng dựa trên:

  • Học thuyết Sinh vật địa lý quần thể và Sinh thái phát sinh quần thể của Sukachev V. và Thái Văn Trừng (1978, 1999).
  • Lý thuyết Cực đỉnh khí hậu (Climax Theory) của Clements F. E. và Braun E. L. (1947).
  • Khung phân loại thảm thực vật quốc tế của UNESCO (1973) và tiếp cận đo lường đa dạng loài của Whittaker R. H. (1975).
  • Hệ thống phân loại dạng sống thích nghi sinh thái của Raunkiaer C. (1934).

Đóng góp đột phá của công trình được định lượng cụ thể: Xây dựng danh lục chuẩn hóa gồm 1.246 loài, 688 chi, 180 họ thuộc 5 ngành thực vật bậc cao có mạch; bổ sung chính thức 01 loài thực vật mới cho Hệ thực vật Việt Nam là Dó vân nam (Aquilaria yunnanensis S. Huang, thuộc họ Trầm Thymelaeaceae); đồng thời bổ sung 218 loài, 71 chi, 12 họ mới cho danh lục Khu BTTN Đồng Sơn – Kỳ Thượng. Nghiên cứu thực hiện trong khung thời gian thực địa từ tháng 8/2015 đến tháng 12/2019, bao phủ diện tích 15.593,81 ha trải dài trên 5 xã: Đồng Sơn, Kỳ Thượng, Đồng Lâm, Vũ Oai và Hòa Bình.

Literature Review và Positioning

Nghiên cứu thảm thực vật và hệ thực vật trên thế giới đã trải qua quá trình tiến hóa lý thuyết sâu sắc. Ở quy mô quốc tế, các trường phái phân loại thảm thực vật ban đầu tập trung vào ngoại mạo sinh thái như Rubel (1912, 1930), Ellenberg (1967) và Shimwell (1972). Phân loại thảm thực vật nhiệt đới được định hình bởi Schimper (1898), Champion (1936) với việc chia đai nhiệt độ, và Aubréville A. với hệ thống phân chia theo độ tàn che. Đến năm 1973, UNESCO đã ban hành khung phân loại toàn cầu dựa trên cấu trúc, ngoại mạo từ bậc Lớp quần hệ (Formation Class) đến Phân quần hệ (Subformation). Về mặt định lượng, Whittaker (1975) đã đặt nền móng phân tách đa dạng alpha ($\alpha$), beta ($\beta$) và gamma ($\gamma$), mở đường cho việc ứng dụng chỉ số Shannon - Wiener ($H'$), Simpson ($Cd$), và gần đây là Entropy Rényi ($H_\alpha$) trong phân tích diễn thế phục hồi rừng (Bautista-Cruz et al., 2007).

Tại Việt Nam, các công trình kinh điển của Chevalier A., Maurand P. (1953), Loschau (1960), Trần Ngũ Phương (1970) và đặc biệt là công trình "Các hệ sinh thái rừng nhiệt đới ở Việt Nam" của Thái Văn Trừng (1978, 1999) đã xác lập hệ thống phân loại sinh thái phát sinh quần thể gồm 14 kiểu thảm thực vật chính. Các nghiên cứu kiểm kê taxon thực vật bậc cao có mạch gắn liền với các công trình của Lecomte H. (1907–1952), Takhtajan A. (1978, 1997, 2009), Phạm Hoàng Hộ (1999–2003) với 11.611 loài, Nguyễn Nghĩa Thìn (2004, 2007) với 11.428 loài, và Trung tâm Nghiên cứu Tài nguyên & Môi trường (2001, 2003, 2005) với Danh lục các loài thực vật Việt Nam (gần 20.000 loài).

Trong y văn tồn tại các cuộc tranh luận học thuật sâu sắc:

  1. Tranh luận về tiêu chí phân loại thảm thực vật: Trường phái cấu trúc - ngoại mạo (UNESCO, 1973; Phan Kế Lộc, 1985) nhấn mạnh tính trực quan và khả năng thể hiện trên bản đồ viễn thám, trong khi trường phái sinh thái phát sinh (Thái Văn Trừng, 1978, 1999) cho rằng nếu tách rời các nhân tố đá mẹ, thổ nhưỡng và lịch sử di cư của khu hệ thực vật thì không thể phản ánh đúng bản chất diễn thế và tổ thành rừng nhiệt đới Việt Nam.
  2. Tranh luận về tác động nhân sinh đối với đa dạng loài: Hoang Van Sam et al. (2009, 2011) tại VQG Bến En lập luận rằng các nhân tố môi trường và con người tác động độc lập; sự can thiệp của con người làm suy giảm nghiêm trọng mật độ các loài cây gỗ kinh tế và cây Sách Đỏ nhưng không làm thay đổi đáng kể chỉ số đa dạng loài tổng thể ($H'$). Ngược lại, các nghiên cứu của MacKinnon et al. (1992) và Trần Thế Liên (2002) cho rằng áp lực nhân sinh trực tiếp phá vỡ cân bằng sinh thái, làm suy thoái cấu trúc tầng tán và đẩy nhanh tốc độ tuyệt chủng cục bộ.

Định vị nghiên cứu trong dòng chảy học thuật quốc tế được thể hiện qua bảng đối sánh thực nghiệm:

  • So sánh với nghiên cứu Vườn quốc gia Doi Suthep-Pui (miền Bắc Thái Lan; Maxwell et al.): Với diện tích 261 $\text{km}^2$, Doi Suthep-Pui ghi nhận 2.084 loài thực vật có mạch, trong đó kiểu rừng thường xanh chiếm ưu thế với 930 loài. Tại Đồng Sơn – Kỳ Thượng (155,9 $\text{km}^2$), mật độ ghi nhận đạt 1.246 loài, thể hiện mật độ đa dạng loài trên đơn vị diện tích tương đương và đặc trưng cho đai rừng kín thường xanh núi thấp Đông Bắc Á.
  • So sánh với nghiên cứu của Gagnepain (1926) và Takhtajan (1978) về khu hệ thực vật Đông Dương: Luận án tái khẳng định tính chất giao thoa giữa luồng thực vật bản địa Bắc Việt Nam - Nam Trung Hoa (Dunn, 1912; Hu Shiu-ying, 2008) với sự hiện diện áp đảo của các họ Euphorbiaceae, Lauraceae, Rubiaceae, Orchidaceae.

Công trình này thúc đẩy chuyên ngành Quản lý Tài nguyên Rừng bằng việc thiết lập cơ sở dữ liệu thực vật học chuẩn hóa đầu tiên cho một khu bảo tồn vùng Đông Bắc, cung cấp bằng chứng định lượng về diễn thế phục hồi sau nương rẫy thông qua chỉ số Entropy Rényi đa chiều.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đóng góp vào việc kiểm chứng và mở rộng các lý thuyết lâm học và sinh thái học kinh điển:

  • Mở rộng lý thuyết Sinh thái phát sinh quần thể của Thái Văn Trừng (1978, 1999): Nghiên cứu đã phân định và mô tả chi tiết 2 kiểu thảm thực vật chính tại vùng Đông Bắc: Kiểu rừng kín thường xanh mưa ẩm nhiệt đới (ở đai cao <700m) và Kiểu rừng kín thường xanh mưa ẩm á nhiệt đới núi thấp (ở đai cao >700m), bổ sung các phân loại chi tiết cho các kiểu phụ thứ sinh nhân tác (phục hồi sau khai thác kiệt và phục hồi sau nương rẫy, cháy rừng).
  • Kiểm chứng Lý thuyết Cực đỉnh sinh thái (Clements, 1916): Các dẫn liệu về tổ thành tầng cây gỗ và chỉ số tương đồng ($SI$) chứng minh rằng rừng thứ sinh phục hồi sau nương rẫy có quỹ đạo phục hồi chậm hơn và cấu trúc đơn giản hơn nhiều so với rừng phục hồi sau khai thác chọn, do sự suy thoái nghiêm trọng của tầng đất mặt Feralit giàu mùn.
  • Kiểm chứng Mô hình Phổ sinh học thích nghi của Raunkiaer (1934): Dữ liệu thực tế tại khu vực thiết lập một phổ dạng sống sinh thái đặc thù cho vùng núi Đông Bắc Việt Nam, khẳng định sự thống trị tuyệt đối của nhóm chồi trên ($Ph$) trong điều kiện khí hậu nhiệt đới ẩm gió mùa.

Mô hình lý thuyết được lượng hóa thông qua hệ thống mệnh đề khoa học:

  • Mệnh đề 1 ($P1$): Tính đa dạng cấu trúc loài tầng cây gỗ tỷ lệ thuận với tính ổn định của tầng đất và thời gian phục hồi không có tác động xâm lấn của con người.
  • Mệnh đề 2 ($P2$): Tỷ lệ các loài thực vật có giá trị dược liệu và lâm sản ngoài gỗ phân bố tập trung ở đai nhiệt đới núi thấp (<700m) chịu rủi ro suy thoái cao gấp 3,2 lần so với đai á nhiệt đới (>700m) do khả năng tiếp cận địa hình của cộng đồng địa phương.

Khung phân tích độc đáo

Khung phân tích của luận án là sự tích hợp liên ngành giữa 3 trụ cột lý thuyết: (1) Lý thuyết cấu trúc quần xã thực vật định lượng (Whittaker, Curtis & McIntosh); (2) Phân loại học hình thái kết hợp địa lý sinh vật (Takhtajan, APG); (3) Sinh thái học nhân văn và kinh tế tài nguyên rừng (Social-Ecological Systems Framework).

+-------------------------------------------------------------------------+
|                  KHUNG PHÂN TÍCH LIÊN NGÀNH CỦA LUẬN ÁN                 |
+-------------------------------------------------------------------------+
                                     |
       +-----------------------------+-----------------------------+
       |                                                           |
       v                                                           v
+------------------------------------+             +------------------------------------+
|  QUẦN XÃ THỰC VẬT ĐỊNH LƯỢNG       |             |  PHÂN LOẠI & ĐỊA LÝ SINH VẬT       |
|  - Thảm thực vật (Thái Văn Trừng)  |             |  - 5 ngành thực vật bậc cao        |
|  - Chỉ số IVI, H', Cd, SI          |             |  - Phổ dạng sống Raunkiaer (SB)    |
|  - Entropy Rényi (H_alpha)         |             |  - Đánh giá Sách Đỏ, IUCN, NĐ 06   |
+------------------------------------+             +------------------------------------+
       |                                                           |
       +-----------------------------+-----------------------------+
                                     |
                                     v
                 +---------------------------------------+
                 |  TƯƠNG TÁC NHÂN SINH & KINH TẾ XÃ HỘI |
                 |  - Khảo sát PRA tại 5 xã vùng đệm     |
                 |  - Sức ép sinh kế người Dao (97%)     |
                 |  - Mức thiếu hụt lương thực (210kg/ng)|
                 +---------------------------------------+
                                     |
                                     v
                 +---------------------------------------+
                 |  GIẢI PHÁP BẢO TỒN TÀI NGUYÊN BỀN VỮNG|
                 +---------------------------------------+

Khung phân tích xác lập rõ các điều kiện biên (Boundary conditions): Đối tượng nghiên cứu giới hạn ở thực vật bậc cao có mạch tự nhiên và cây gỗ rừng trồng trong phạm vi ranh giới hành chính 15.593,81 ha của Khu BTTN Đồng Sơn – Kỳ Thượng, với giới hạn đai cao tối đa là 1.096m (đỉnh Thiên Sơn).

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu vận dụng thế giới quan thực chứng (Positivism) kết hợp với chủ nghĩa thực tại sinh thái (Ecological Realism), áp dụng phương pháp hỗn hợp (Mixed Methods Research Design) đan xen giữa định lượng lâm học chuẩn xác và điều tra xã hội học nông thôn có sự tham gia (PRA).

Thiết kế đa cấp độ (Multi-level design) được cấu trúc chặt chẽ:

  • Cấp độ cảnh quan/quần xã (Macro-level): Phân tích các kiểu thảm thực vật rừng dựa trên phân tích giải đoán ảnh vệ tinh LANDSAT TM5 (độ phân giải không gian 25m), đối soát bản đồ địa hình và mạng lưới thủy văn lưu vực sông Ba Chẽ và sông Man.
  • Cấp độ lâm phần/quần thể (Meso-level): Thiết lập hệ thống ô tiêu chuẩn (OTC) định lượng diện tích $1.000\text{ m}^2$ ($20\text{m} \times 50\text{m}$ hoặc $25\text{m} \times 40\text{m}$) trên các tuyến điều tra đại diện theo các đai độ cao và trạng thái rừng.
  • Cấp độ cá thể/loài (Micro-level): Thu thập mẫu tiêu bản thực vật, đo đếm đường kính thân cây ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), chiều cao dưới cành ($H_{dc}$), đường kính tán ($D_t$) và thiết lập các ô dạng bản (ODB) $4\text{ m}^2$ ($2\text{m} \times 2\text{m}$) để đánh giá tầng cây tái sinh và thảm tươi.

Quy trình nghiên cứu rigorous

Chiến lược lấy mẫu thực địa tuân thủ nghiêm ngặt các tiêu chuẩn sinh thái học rừng:

  • Tuyến điều tra thực địa: Thiết lập mạng lưới các tuyến điều tra cắt ngang các dạng địa hình, dông núi chính (dông Mo 852,5m; Am Váp 1.051m; Thiên Sơn 1.096m; đèo Dài 597,3m; khe Ru 566,9m) đi qua các sinh cảnh rừng giàu, rừng trung bình, rừng nghèo, rừng phục hồi sau khai thác và nương rẫy.
  • Thu thập và xử lý mẫu vật: Mẫu tiêu bản thực vật được thu thập theo đúng quy chuẩn thực vật học (mẫu có hoa, quả, chồi non, ghi chép màu sắc, mùi vị, nhựa mủ). Tiêu bản được xử lý cồn, sấy khô và định danh tại Phòng Tiêu bản Thực vật (HNU và VNF) thuộc Trường Đại học Lâm nghiệp và Viện Sinh thái và Tài nguyên Sinh vật.
  • Phương pháp chuyên gia và đối chiếu danh pháp: Giám định tên khoa học thông qua bộ Cây cỏ Việt Nam (Phạm Hoàng Hộ), Thực vật chí Việt Nam (11 tập), Cẩm nang nhận biết các họ hạt kín (Nguyễn Tiến Bân), Flora of China (Wu Zheng-yi & Raven) và phân loại hiện đại APG IV (2016).
  • Kiểm soát độ tin cậy và giá trị: Dữ liệu điều tra xã hội học được thực hiện thông qua phỏng vấn bán cấu trúc 34 hộ dân sống trong vùng lõi và hàng trăm lượt hộ dân tại 22 thôn bản thuộc 5 xã lân cận, đảm bảo phép tam giác đạc dữ liệu (Data Triangulation) giữa số liệu báo cáo của Hạt Kiểm lâm, biên bản vi phạm lâm luật giai đoạn 2015–2018 và lời khai thực tế của người dân.

Data và phân tích

Đặc điểm mẫu nghiên cứu: Toàn bộ dữ liệu lâm học thu thập trên diện tích điều tra đại diện cho 14.279,8 ha rừng tự nhiên. Dữ liệu xã hội học phản ánh thực trạng của 1.504 người dân địa phương với cơ cấu 79,5% – 97% là dân tộc Dao.

Các công thức toán học và chỉ số định lượng áp dụng:

  1. Chỉ số giá trị quan trọng (Importance Value Index - IVI) của Curtis & McIntosh (1950): $$IVI (%) = \frac{N% + G% + F%}{3}$$ Trong đó: $N%$ là mật độ tương đối; $G%$ là diện tích tiết diện thân tương đối; $F%$ là tần số xuất hiện tương đối.
  2. Chỉ số đa dạng loài Shannon - Wiener ($H'$): $$H' = -\sum_{i=1}^{S} P_i \ln(P_i) \quad \text{với } P_i = \frac{n_i}{N}$$
  3. Chỉ số mức độ chiếm ưu thế Simpson ($Cd$): $$Cd = \sum_{i=1}^{S} P_i^2$$
  4. Chỉ số tương đồng Sorensen ($SI$): $$SI = \frac{2C}{A + B} \times 100%$$ Trong đó: $A, B$ là số loài của hai quần xã so sánh; $C$ là số loài chung.
  5. Chỉ số Entropy Rényi đa chiều ($H_\alpha$): $$H_\alpha = \frac{1}{1-\alpha} \ln \left(\sum_{i=1}^{S} P_i^\alpha \right)$$ Được phân tích với các tham số $\alpha = 0, 0.5, 1, 2, \infty$ nhằm đánh giá toàn diện độ phong phú loài và mức độ đồng đều mà không bị thiên lệch bởi các loài cực kỳ hiếm hay cực kỳ phổ biến.
  6. Phần mềm phân tích: Dữ liệu được tính toán và mô hình hóa bằng phần mềm Microsoft Excel chuyên dụng cho lâm nghiệp kết hợp các gói phân tích sinh thái học và phần mềm MapInfo/ArcGIS để số hóa bản đồ tuyến điều tra và phân bố thảm thực vật.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Phát hiện và công bố loài thực vật bổ sung cho hệ thực vật Việt Nam: Luận án đã phát hiện, thu thập mẫu chuẩn và giám định thành công loài Dó vân nam (Aquilaria yunnanensis S. Huang), thuộc họ Trầm (Thymelaeaceae). Đây là ghi nhận đầu tiên xác lập sự hiện diện của loài này tại Việt Nam, mở rộng vùng phân bố địa lý sinh vật của chi Aquilaria từ vùng Tây Nam Trung Quốc xuống vùng núi Đông Bắc Việt Nam. Đồng thời, nghiên cứu bổ sung tới 218 loài, 71 chi và 12 họ thực vật bậc cao có mạch cho danh lục Khu BTTN Đồng Sơn – Kỳ Thượng.

  2. Thống kê chuẩn hóa cấu trúc hệ thực vật: Xác lập danh lục toàn diện với 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:

    • Ngành Mộc lan (Magnoliophyta): Chiếm ưu thế tuyệt đối với 1.157 loài (92,86%), 629 chi (91,42%), 151 họ (83,89%).
    • Ngành Dương xỉ (Polypodiophyta): 65 loài (5,22%), 43 chi, 20 họ.
    • Ngành Thông (Pinophyta): 14 loài (1,12%), 10 chi, 6 họ.
    • Ngành Thông đất (Lycopodiophyta): 9 loài (0,72%), 5 chi, 2 họ.
    • Ngành Khuyết lá thông (Psilotophyta): 1 loài (0,08%), 1 chi, 1 họ.

  3. Tính toán phổ dạng sống (Biological Spectrum - SB): Thiết lập phổ dạng sống thực nghiệm cho hệ thực vật Đồng Sơn – Kỳ Thượng: $$SB = 78,41% Ph + 4,09% Ch + 6,58% Hm + 5,38% Cr + 5,54% Th$$ Trong nhóm cây chồi trên ($Ph$), các loài cây chồi trên nhỏ ($Mi$) chiếm 32,42%, cây chồi trên vừa ($Me$) chiếm 22,55%, cây gỗ lớn ($Mg$) chiếm 8,27%, dây leo thân gỗ ($Lp$) chiếm 16,84% và phụ sinh ($Ep$) chiếm 8,30%. Kết quả này phản ánh đặc trưng sinh thái của rừng mưa nhiệt đới gió mùa điển hình với tầng tán rậm rạp và độ ẩm không khí cao quanh năm.

  4. Đa dạng công dụng và nguồn gen nguy cấp:

    • Giá trị sử dụng: Trong số 1.246 loài, có tới 892 loài cây có ích (chiếm 71,59% tổng số loài), phân thành các nhóm: Cây làm thuốc có 643 loài (51,60%); Cây cho gỗ có 271 loài (21,75%); Cây làm cảnh có 168 loài (13,48%); Cây cho lương thực, thực phẩm có 138 loài (11,08%).
    • Tình trạng bảo tồn: Ghi nhận 69 loài thực vật quý hiếm nguy cấp, trong đó có 43 loài trong Sách Đỏ Việt Nam (2007), 31 loài trong Danh mục Đỏ IUCN (2019), 23 loài thuộc Nghị định 06/2019/NĐ-CP và 16 loài thuộc các phụ lục CITES (như Đỉnh tùng Cephalotaxus mannii, Kim giao Nageia fleuryi, Ba kích Morinda officinalis, Trầm hương Aquilaria crassna).

  5. Phân hóa đa dạng tầng cây gỗ theo đai cao và trạng thái diễn thế: Chỉ số đa dạng Shannon-Wiener ($H'$) ở đai nhiệt đới (<700m) đạt giá trị trung bình từ 2,85 đến 3,42; trong khi ở đai á nhiệt đới (>700m) dao động từ 2,45 đến 2,98. Rừng phục hồi sau khai thác có chỉ số tương đồng $SI$ so với rừng nguyên sinh đạt 62,4%, trong khi rừng phục hồi sau nương rẫy chỉ đạt $SI = 38,7%$, với ưu thế thuộc về các loài cây tiên phong ưa sáng kém giá trị như Hu đay (Trema orientalis), Màng tang (Litsea cubeba).

Implications đa chiều

  • Về lý thuyết: Cung cấp cứ liệu thực chứng củng cố lý thuyết về mối quan hệ giữa gradient độ cao và tính đa dạng loài ở vùng nhiệt đới gió mùa ven biển; khẳng định vai trò của dải rừng Đông Triều như một hành lang sinh học quan trọng lưu giữ các nguồn gen cổ sinh học hạt trần và hạt kín nguyên thủy.
  • Về phương pháp luận: Chứng minh tính ưu việt của việc kết hợp chỉ số Entropy Rényi với các chỉ số truyền thống ($H', Cd$) trong việc loại bỏ độ nhiễu thống kê khi đánh giá diễn thế rừng sau nương rẫy.
  • Về ứng dụng thực tiễn: Cung cấp danh mục 643 loài cây thuốc bản địa và định vị các điểm phân bố tập trung của 10 loài dược liệu quý giá trị cao (như Ba kích, Sa nhân tím, Khôi nhung), tạo tiền đề cho các dự án phát triển dược liệu dưới tán rừng.
  • Về hoạch định chính sách: Cung cấp cơ sở khoa học để Ban Quản lý Khu bảo tồn và UBND thành phố Hạ Long điều chỉnh quy hoạch phân khu chức năng (phân khu bảo vệ nghiêm ngặt, phân khu phục hồi sinh thái) và xây dựng chính sách giao khoán bảo vệ rừng gắn liền với hỗ trợ lương thực cho 34 hộ vùng lõi và 22 thôn bản vùng đệm.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn ghi nhận một số giới hạn nghiên cứu (Limitations):

  1. Giới hạn về phân tích phân tử sinh học: Việc xác định danh tính loài Dó vân nam và các taxon mới chủ yếu dựa trên phương pháp hình thái học so sánh đối chiếu mẫu chuẩn (Morphological taxonomy); chưa thực hiện giải trình tự gen DNA Barcoding (các đoạn gen rbcL, matK, ITS) để xác lập cây phát sinh chủng loại phân tử chi tiết.
  2. Độ bao phủ theo mùa: Do địa hình chia cắt hiểm trở và thời tiết mưa lũ cực đoan tập trung vào tháng 7–8 (chiếm 80% lượng mưa năm), việc thu thập mẫu quả/hoa của một số loài thực vật biểu sinh và thân thảo ngắn ngày trong mùa đông khô hạn chưa bao phủ trọn vẹn 100% chu kỳ sinh học.
  3. Định lượng sinh khối và trữ lượng carbon: Luận án tập trung chuyên sâu vào đa dạng thành phần loài và cấu trúc quần xã tầng cây gỗ, chưa tích hợp mô hình đo đếm sinh khối trên mặt đất (AGB) và khả năng hấp thụ carbon của các trạng thái rừng.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  • Hướng 1: Ứng dụng kỹ thuật chỉ thị phân tử DNA Barcoding và hệ gen học để 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 tại khu vực.
  • Hướng 2: Thiết lập các ô định vị sinh thái học lâu dài (Permanent Forest Dynamics Plots - 1 ha hoặc 5 ha) theo chuẩn mạng lưới ForestGEO/CTFS để giám sát biến động động thái rừng, tỷ lệ tử vong và sinh trưởng của cây gỗ dưới tác động của biến đổi khí hậu.
  • Hướng 3: Đánh giá chi tiết giá trị dịch vụ hệ sinh thái (Ecosystem Services Valuation), đặc biệt là chức năng phòng hộ đầu nguồn cung cấp nước ngọt cho hồ Cao Vân (dung tích 5.000 $\text{m}^3$) phục vụ sinh hoạt cho toàn bộ đô thị Hạ Long.
  • Hướng 4: Nghiên cứu thực nghiệm các mô hình nhân giống vô tính và thuần dưỡng các loài cây thuốc nguy cấp có giá trị thương phẩm cao phục vụ sinh kế cộng đồng bản địa.

Tác động và ảnh hưởng

Công trình tạo ra những tác động lan tỏa trên nhiều bình diện:

  • Ảnh hưởng học thuật: Dự báo sẽ đóng góp số lượng trích dẫn quan trọng cho các công trình nghiên cứu về thực vật chí Việt Nam, phân loại học họ Trầm (Thymelaeaceae) và quản lý rừng đặc dụng vùng Đông Bắc. Cung cấp dữ liệu chuẩn mực cho cơ sở dữ liệu đa dạng sinh học quốc gia.
  • Tác động đến chính sách quản lý: Kết quả của luận án đã trực tiếp phục vụ cho công tác rà soát, điều chỉnh quy hoạch rừng đặc dụng giai đoạn 2020–2030 của tỉnh Quảng Ninh theo Quyết định 2041/QĐ-UB; cung cấp luận cứ thực tiễn để chuyển đổi diện tích rừng nghèo kiệt sang trạng thái khoanh nuôi xúc tiến tái sinh tự nhiên có làm giàu rừng.
  • Lợi ích kinh tế - xã hội: Tạo cơ sở khoa học để triển khai các chương trình chi trả dịch vụ môi trường rừng (PES) và xây dựng chuỗi giá trị lâm sản ngoài gỗ (LSNG), giúp giải quyết bài toán thiếu đói lúc giáp hạt (thiếu lương thực 3–4 tháng/năm) cho đồng bào dân tộc Dao, giảm thiểu các vụ vi phạm khai thác gỗ và săn bắt trái phép (giảm tải số vụ vi phạm lâm luật ghi nhận từ 2015–2018).

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giới học thuật lâm nghiệp: Thừa hưởng hệ thống dữ liệu sơ cấp chuẩn mực về 1.246 loài thực vật, phương pháp luận tích hợp giữa chỉ số $H'$ truyền thống và Entropy Rényi, cùng hệ thống tọa độ các tuyến điều tra thực địa để tiếp tục triển khai các đề tài chuyên sâu.
  • Ban Quản lý Khu BTTN Đồng Sơn – Kỳ Thượng và Hạt Kiểm lâm: Sở hữu bộ cẩm nang tra cứu và danh lục các loài nguy cấp, bản đồ hiện trạng phân bố các loài quý hiếm để tổ chức tuần tra rừng có trọng tâm, bảo vệ nghiêm ngặt các điểm nóng đa dạng sinh học.
  • Chính quyền địa phương (TP Hạ Long và tỉnh Quảng Ninh): Có cơ sở dữ liệu khoa học chính xác để ban hành các quyết định đầu tư cơ sở hạ tầng, phát triển du lịch sinh thái có kiểm soát tại đỉnh Thiên Sơn - Am Váp gắn liền với các chương trình mục tiêu quốc gia về phát triển kinh tế vùng đồng bào dân tộc thiểu số.
  • Cộng đồng dân cư vùng đệm (đặc biệt là 22 thôn bản thuộc 5 xã): Được tiếp cận các hướng dẫn kỹ thuật bảo tồn tri thức bản địa về 643 loài cây thuốc, mở ra cơ hội tham gia các hợp tác xã trồng và sơ chế dược liệu dưới tán rừng tự nhiên.

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? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc chuẩn hóa và mở rộng Lý thuyết Sinh thái phát sinh quần thể của Thái Văn Trừng (1978, 1999) vào cấu trúc đai cao đặc thù của vùng núi đất Đông Bắc chịu ảnh hưởng khí hậu biển. Luận án đã mô hình hóa sự thay thế của các ưu hợp thực vật từ đai nhiệt đới (<700m) lên đai á nhiệt đới (>700m) và lượng hóa chính xác phổ dạng sống sinh thái $SB = 78,41% Ph + 4,09% Ch + 6,58% Hm + 5,38% Cr + 5,54% Th$, chứng minh tính quy luật của sự thích nghi hình thái thực vật trước gió mùa lạnh Đông Bắc.

2. Đổi mới phương pháp luận của luận án khi so sánh với ít nhất 2 công trình trước đây? Trả lời: So với các nghiên cứu trước đây tại khu vực chỉ dừng ở phương pháp kiểm kê tuyến định tính đơn thuần (như Đỗ Xuân Trường, 2011 chỉ khảo sát cây gỗ; VCF, 2011 khảo sát nhanh), luận án đã: (1) Ứng dụng hệ thống ô tiêu chuẩn định lượng đa cấp ($1.000\text{ m}^2$ và $4\text{ m}^2$) phân tích đầy đủ các chỉ số $IVI, H', Cd, SI$; (2) Tiên phong ứng dụng chỉ số Entropy Rényi ($H_\alpha$) đa chiều theo phương pháp của Bautista-Cruz et al. (2007) để đánh giá độ dốc phân hóa loài giữa các trạng thái diễn thế rừng phục hồi sau nương rẫy và sau khai thác kiệt.

3. Phát hiện thực nghiệm đáng ngạc nhiên nhất là gì và dữ liệu nào chứng minh? Trả lời: Phát hiện bất ngờ nhất là việc tìm thấy loài Dó vân nam (Aquilaria yunnanensis S. Huang) – một loài trước đây chỉ được ghi nhận phân bố hẹp tại các vùng núi đá vôi và nhiệt đới thuộc tỉnh Vân Nam (Trung Quốc), nay xuất hiện tự nhiên tại rừng kín thường xanh mưa ẩm núi thấp Đồng Sơn – Kỳ Thượng. Dữ liệu hình thái hoa, quả và cấu trúc gỗ đã khẳng định loài mới bổ sung cho Hệ thực vật Việt Nam, nâng tổng số loài bổ sung cho khu bảo tồn lên 218 loài (tăng 21,2% so với các danh lục tổng hợp trước đó).

4. Luận án có cung cấp quy trình tái lặp nghiên cứu (Replication Protocol) không? Trả lời: Có. Luận án mô tả chi tiết quy trình tái lặp thực nghiệm bao gồm: Hệ thống tọa độ địa lý các đỉnh núi chuẩn (Thiên Sơn 1.096m, Am Váp 1.051m, Đèo Mo 852,5m), danh sách chi tiết các tuyến điều tra cắt qua 5 xã, kích thước và phương pháp bố trí ô tiêu chuẩn $1.000\text{ m}^2$, tiêu chuẩn đo đếm cây gỗ ($D_{1.3} \ge 6\text{ cm}$), quy trình xử lý cồn/sấy tiêu bản và thuật toán tính toán các chỉ số sinh thái học thông qua công thức toán học chuẩn hóa.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào? Trả lời: Lộ trình 10 năm tập trung vào 3 trọng tâm: (1) Giai đoạn 2021–2024: Số hóa cơ sở dữ liệu mẫu vật (Herbarium database) và thực hiện giải trình tự gen DNA Barcoding toàn bộ 69 loài thực vật nguy cấp; (2) Giai đoạn 2025–2028: Xây dựng ô định vị sinh thái vĩnh cửu 1 ha để theo dõi động thái quần xã và chu trình dinh dưỡng đất rừng Feralit có mùn; (3) Giai đoạn 2029–2030: Hoàn thiện mô hình đồng quản lý rừng và thương mại hóa bền vững 5 loài cây thuốc bản địa mũi nhọn nhằm xóa bỏ hoàn toàn áp lực xâm hại vùng lõi khu bảo tồn.

Kết luận

  1. Luận án đã xây dựng thành công bộ dữ liệu khoa học toàn diện, chuẩn hóa và cập nhật nhất về tính đa dạng thực vật tại Khu BTTN Đồng Sơn – Kỳ Thượng với 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.
  2. Công trình đạt được phát hiện mang tầm quốc gia khi bổ sung 01 loài mới cho Hệ thực vật Việt Nam (Aquilaria yunnanensis S. Huang), đồng thời bổ sung 218 loài, 71 chi, 12 họ mới cho khu hệ thực vật địa phương.
  3. Nghiên cứu đã xác lập và mô tả chi tiết cấu trúc 2 kiểu thảm thực vật tự nhiên chính (Rừng kín thường xanh mưa ẩm nhiệt đới và Rừng kín thường xanh mưa ẩm á nhiệt đới núi thấp) cùng các trạng thái thứ sinh diễn thế, lượng hóa thành công các chỉ số đa dạng sinh học ($H', Cd, IVI, SI$) và Entropy Rényi tầng cây gỗ.
  4. Thiết lập phổ dạng sống thực nghiệm $SB = 78,41% Ph + 4,09% Ch + 6,58% Hm + 5,38% Cr + 5,54% Th$, phản ánh bản chất sinh thái của hệ sinh thái rừng mưa nhiệt đới gió mùa điển hình vùng Đông Bắc.
  5. Kiểm kê, phân loại giá trị sử dụng của 892 loài cây có ích (trong đó có 643 loài cây thuốc) và đánh giá hiện trạng bảo tồn của 69 loài thực vật nguy cấp ghi danh trong Sách Đỏ Việt Nam, Danh lục Đỏ IUCN và Nghị định 06/2019/NĐ-CP.
  6. Xác định rõ nguyên nhân trực tiếp và gián tiếp gây suy giảm tài nguyên rừng xuất phát từ tình trạng kinh tế - xã hội khó khăn của cộng đồng người Dao bản địa (chiếm 79,5% – 97% dân số, thiếu lương thực 3–4 tháng/năm), từ đó đề xuất hệ thống 4 nhóm giải pháp đồng bộ về quản lý bảo vệ, phát triển sinh kế, nâng cao nhận thức và ứng dụng khoa học kỹ thuật lâm nghiệp nhằm bảo tồn và phát triển bền vững tài nguyên thực vật khu vực.