Tổng quan về luận án

Nghiên cứu cấu trúc rừng tự nhiên và tính đa dạng sinh học là nền tảng cốt lõi trong khoa học lâm sinh nhằm giải mã quy luật vận động nội tại của quần xã thực vật, từ đó xây dựng các giải pháp phục hồi và quản trị rừng bền vững. Luận án tiến sĩ lâm nghiệp của nghiên cứu sinh Quốc Minh Dũng (2025) với đề tài "Nghiên cứu cấu trúc rừng lá rộng thường xanh và tính đa dạng thực vật thân gỗ tại Vườn quốc gia Phia Oắc - Phia Đén tỉnh Cao Bằng" (Mã số: 9620205; Người hướng dẫn: PGS. Trần Quốc Hưng và PGS. Nguyễn Hồng Hải) đã tiếp cận toàn diện một trong những hệ sinh thái rừng núi cao á nhiệt đới và nhiệt đới gió mùa đặc thù nhất vùng Đông Bắc Việt Nam. Vườn quốc gia (VQG) Phia Oắc - Phia Đén sở hữu tổng diện tích tự nhiên 10.593,5 ha với độ che phủ rừng lên đến 81%, trong đó diện tích rừng gỗ tự nhiên đạt 6.146,5 ha nhưng phần lớn diện tích (hơn 4.000 ha) lại thuộc trạng thái rừng nghèo và nghèo kiệt do lịch sử khai thác và các tác động nhân vi kéo dài.

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án hướng tới giải quyết là sự thiếu vắng các công trình định lượng dài hạn ở quy mô vi cấu trúc không gian và lịch sử tiến hóa quần xã. Phần lớn các công trình tiền trạm tại khu vực chỉ dừng lại ở khảo sát tĩnh, phạm vi hẹp hoặc kiểm kê danh lục mà chưa tích hợp phân tích mô hình điểm không gian (spatial point pattern analysis), cấu trúc phát sinh chủng loại (phylogenetic community structure) và phân bố độ phong phú loài (Species Abundance Distribution - SAD) để định lượng động thái lâm phần giai đoạn 2015–2020.

Luận án thiết lập 4 câu hỏi nghiên cứu ($RQ$) và 3 giả thuyết khoa học ($H$):

  • $RQ_1$: Cấu trúc không gian, cấu trúc phân bố đường kính ($N/D_{1.3}$) và chiều cao ($N/H_{vn}$) của lâm phần rừng lá rộng thường xanh (LRTX) nghèo và nghèo kiệt mang đặc trưng gì?
  • $RQ_2$: Mức độ đa dạng loài và cấu trúc phát sinh chủng loại của thực vật thân gỗ chịu sự chi phối của những cơ chế sinh thái nào?
  • $RQ_3$: Biến động về mật độ, sinh khối trên mặt đất (AGB), tỷ lệ chết và bổ sung diễn ra theo chiều hướng nào trong chu kỳ 5 năm (2015–2020)?
  • $RQ_4$: Những giải pháp lâm sinh nào có cơ sở khoa học vững chắc để xúc tiến quá trình phục hồi các trạng thái rừng này?
  • $H_1$: Phân bố không gian của các loài cây gỗ chủ yếu thể hiện dạng phân bố cụm ở cự ly nhỏ do phát tán hạn chế và chuyển dần sang phân bố ngẫu nhiên ở cấp kích thước lớn do cạnh tranh tự tỉa thưa.
  • $H_2$: Quần xã thực vật rừng nghèo và nghèo kiệt thể hiện sự co cụm phát sinh loài (phylogenetic clustering) dưới tác động của cơ chế lọc môi trường khắc nghiệt vùng núi cao.
  • $H_3$: Quá trình tích lũy sinh khối và hồi phục cấu trúc lâm phần diễn ra theo hướng gia tăng tổng tiết diện ngang ($\sum G/\text{ha}$) dù mật độ cây tầng cao có thể biến thiên do tỷ lệ chết tự nhiên.

Khung lý thuyết của nghiên cứu tích hợp đa chiều: Thuyết trung tính thống nhất về đa dạng sinh học và địa sinh học (Unified Neutral Theory of Biodiversity - Hubbell, 2001), Thuyết phân hóa ổ sinh thái (Niche Differentiation Theory - Tilman, 1982), Thuyết chết phụ thuộc mật độ Janzen-Connell (Janzen, 1970; Connell, 1971), và Thuyết bảo tồn ổ sinh thái phát sinh chủng loại (Phylogenetic Niche Conservatism - Wiens et al., 2010). Đối tượng và phạm vi nghiên cứu được thực hiện trên 02 Ô định vị (ODV) theo dõi sinh thái quốc gia số 12 và 13 tại VQG Phia Oắc - Phia Đén, thu thập dữ liệu định lượng 2 giai đoạn (2015 và 2020), tạo lập cơ sở dữ liệu thực nghiệm chuẩn xác phục vụ bảo tồn và phát triển lâm nghiệp sinh thái.

Literature Review và Positioning

Nghiên cứu cấu trúc rừng tự nhiên trên bình diện quốc tế đã trải qua sự chuyển dịch mạnh mẽ từ các mô hình tĩnh mô tả tầng tán sang phân tích động lực học không gian đa biến và sinh thái học tiến hóa. Trường phái cấu trúc truyền thống tập trung mô hình hóa phân bố thực nghiệm bằng các hàm toán học như Weibull, Meyer hay phân bố Giảm số mũ (Cao Thị Thu Hiền và cs., 2013; Phùng Đình Trung và cs., 2019). Tuy nhiên, các phương pháp này không thể hiện được vị trí tương đối và tương tác cạnh tranh hay tương hỗ giữa các cá thể. Để giải quyết hạn chế này, trường phái phân tích mô hình điểm không gian (spatial point pattern analysis) với việc áp dụng hàm $K$-Ripley và hàm tương quan cặp $g(r)$ đã trở thành công cụ chủ đạo (Getzin và cs., 2008; Lan và cs., 2012; Wiegand & Moloney, 2014).

Về mặt tranh luận học thuật, tồn tại mâu thuẫn sâu sắc giữa hai trường phái: Thuyết ổ sinh thái (Niche Theory) cho rằng chọn lọc môi trường (environmental filtering) và sự phân hóa không gian sống quyết định sự chung sống của các loài (Tilman, 1982); trong khi Thuyết trung tính (Neutral Theory) khẳng định sự đa dạng được duy trì bởi các quá trình ngẫu nhiên sinh - tử, phát tán hạn chế (dispersal limitation) và trôi dạt sinh thái (Hubbell, 2001). Gần đây, việc tích hợp phân tích cấu trúc phát sinh chủng loại (Webb và cs., 2002; Cavender-Bares và cs., 2009) và phân bố độ phong phú loài (SAD) đã mở ra cầu nối để kiểm định sự tương tác giữa hai cơ chế này.

Đặt trong tương quan với các nghiên cứu quốc tế:

  • So sánh với nghiên cứu của Fibich và cs. (2016) trên rừng mưa nhiệt đới nguyên sinh và thứ sinh tại New Guinea, mô hình phân bố không gian của các loài chiếm ưu thế trong rừng thứ sinh thường biểu hiện dạng cụm mạnh mẽ do hạn chế phát tán sau xáo trộn, trái ngược với cấu trúc ngẫu nhiên trong rừng nguyên sinh già cỗi.
  • So sánh với công trình của Lan và cs. (2012) trên ô định vị 20 ha tại Tây Song Bản Nạp (Trung Quốc), tương tác tương hỗ chiếm ưu thế ở cự ly hẹp ($< 5\text{ m}$) giữa các cây con, trong khi áp lực cạnh tranh loại trừ gia tăng rõ rệt ở tầng cây sào và cây thành thục.
  • So sánh với Culmsee và cs. (2012) tại vùng núi Malesia, tính chất co cụm phát sinh loài (phylogenetic clustering) thường chiếm ưu thế ở đai núi thấp và trung bình, phản ánh tác động mạnh mẽ của yếu tố địa hình và khí hậu đến sự quy tụ các nhóm loài có họ hàng gần.

Tại Việt Nam, các nghiên cứu của Nguyễn Hồng Hải và cs. (2014, 2018) tại A Lưới hay Trần Văn Đô và cs. (2017, 2020) trên 6 vùng sinh thái đã tiếp cận mô hình điểm không gian và sinh khối AGB, nhưng việc kết hợp đồng thời phân tích phát sinh chủng loại (dựa trên cây phân loại APG IV), mô hình hóa SAD qua tiêu chuẩn AIC và theo dõi biến động 5 năm trên các ô định vị sinh thái quốc gia tại vùng núi cao Phia Oắc - Phia Đén là một đóng góp tiên phong, định vị rõ nét khoảng trống tri thức mà luận án đã lấp đầy.

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

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

Luận án mang lại những bước tiến quan trọng cho lý thuyết sinh thái học quần xã rừng nhiệt đới gió mùa:

  • Mở rộng Thuyết trung tính và Thuyết ổ sinh thái: Luận án cung cấp bằng chứng thực nghiệm chứng minh rằng cấu trúc rừng phục hồi sau khai thác tại VQG Phia Oắc - Phia Đén không vận hành đơn lẻ theo một cơ chế thuần nhất, mà là sự giao thoa liên tục (niche-neutral continuum). Trong đó, các quá trình ngẫu nhiên chi phối sự phân bố của lớp cây non thông qua hiện tượng phát tán hạn chế, còn cơ chế chọn lọc ổ sinh thái định hình cấu trúc tầng cây cao thành thục.
  • Chứng minh cơ chế bảo vệ nhóm loài (Herd protection mechanism) và tử vong phụ thuộc mật độ (Density-dependent mortality): Phân tích không gian đa biến xác nhận sự chuyển hóa hình thái phân bố từ dạng cụm ở cây non sang dạng ngẫu nhiên ở cây thành thục, cung cấp bằng chứng củng cố giả thuyết Janzen-Connell trên hệ sinh thái rừng núi cao miền Bắc Việt Nam.
  • Đóng góp cho lý thuyết sinh thái học tiến hóa: Việc phát hiện chỉ số họ hàng ròng ($NRI$) và chỉ số taxon gần nhất ($NTI$) mang giá trị dương khẳng định hiện tượng co cụm phát sinh chủng loại (phylogenetic clustering), củng cố giả thuyết về tính bảo tồn ổ sinh thái phát sinh loài (Phylogenetic Niche Conservatism) trong điều kiện lập địa khắc nghiệt vùng núi cao.

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

Khung phân tích của luận án được xây dựng dựa trên sự liên kết 4 chiều cấu trúc:

  • Khung cấu trúc kích thước truyền thống: Tích hợp các hàm phân bố $N/D_{1.3}$ và $N/H_{vn}$ (Weibull 3 tham số, Meyer) để đánh giá tính liên tục của quá trình tái sinh và phân tầng tán.
  • Khung cấu trúc không gian điểm (Spatial Point Pattern Framework): Ứng dụng hàm tương quan cặp $g(r)$ đơn biến và hai biến, kết hợp mô hình vô hiệu hóa điểm Poisson đồng nhất (CSR - Complete Spatial Randomness) và Poisson không đồng nhất (Heterogeneous Poisson) để bóc tách ảnh hưởng của môi trường sống không đồng nhất và tương tác sinh học giữa các cá thể.
  • Khung cấu trúc phát sinh chủng loại (Phylogenetic Community Framework): Xây dựng cây phát sinh chủng loại bằng Phylomatic dựa trên bộ khung siêu cây Megatree APG IV, định lượng đa dạng phát sinh loài ($PD$), khoảng cách cặp trung bình ($MPD$), khoảng cách taxon gần nhất trung bình ($MNTD$), từ đó tính toán $NRI$ và $NTI$.
  • Khung phân bố độ phong phú loài (Species Abundance Distribution - SAD): Mô hình hóa đồng thời 6 hàm phân bố (Broken stick, Lognormal, Zipf, Mandelbrot, Niche preemption, Neutral) và kiểm định bằng tiêu chí thông tin Akaike (AIC) cùng kiểm định Kolmogorov-Smirnov.

Điều kiện biên (boundary conditions): Khung phân tích được xác lập tối ưu cho quần xã thực vật rừng tự nhiên thứ sinh lá rộng thường xanh nghèo và nghèo kiệt tại đai độ cao á nhiệt đới núi trung bình và cao (từ 700 m đến 1.931 m so với mực nước biển), nơi chịu tác động mạnh của địa hình chia cắt, độ dốc $25^\circ - 30^\circ$ và khí hậu nhiệt đới gió mùa có mùa đông lạnh.

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

Thiết kế nghiên cứu

  • Bản thể luận và Triết lý nghiên cứu: Luận án tiếp cận theo trường phái Thực chứng luận (Positivism) kết hợp Chủ nghĩa hiện thực phê phán (Critical Realism), xem xét các cấu trúc và quy luật sinh thái là khách quan, có thể đo lường và lượng hóa chính xác thông qua hệ thống thực nghiệm và thống kê mô hình.
  • Thiết kế đa cấp (Multi-level design):
    • Cấp vĩ mô: Toàn bộ hệ sinh thái rừng LRTX tại VQG Phia Oắc - Phia Đén (10.593,5 ha).
    • Cấp trung mô: 02 Ô định vị sinh thái quốc gia (ODV 12 và ODV 13) đại diện cho trạng thái rừng nghèo và nghèo kiệt, phân chia thành các ô nghiên cứu (ONC) và ô đo đếm (ODD) kích thước $20\text{ m} \times 20\text{ m}$.
    • Cấp vi mô: Từng cá thể cây gỗ ($D_{1.3} \ge 6.0\text{ cm}$) và cây tái sinh được định vị tọa độ cực $(x, y)$ chính xác trong không gian 2D.

Quy trình nghiên cứu rigorous

  • Chiến lược lấy mẫu và thu thập dữ liệu: Kế thừa và mở rộng số liệu điều tra chuẩn hóa trên 02 ODV sinh thái số 12 và 13 trong hai chu kỳ 2015 và 2020 do Phân viện Điều tra - Quy hoạch Rừng Đông Bắc Bộ và NCS thực hiện.
  • Giao thức điều tra tầng cây cao: Đo đếm toàn diện tất cả các cá thể cây có $D_{1.3} \ge 6.0\text{ cm}$; ghi nhận tên loài, đường kính ngang ngực $D_{1.3}$ (thước kẹp kính chuyên dụng), chiều cao vút ngọn $H_{vn}$ (máy đo cao Vertex IV), đường kính tán ($D_t$), tọa độ không gian tương đối $(x, y)$ và chất lượng cây (tốt, trung bình, xấu).
  • Giao thức điều tra tầng tái sinh: Thiết lập các ô dạng bản (ODB) kích thước $5\text{ m} \times 5\text{ m}$ hoặc $2\text{ m} \times 2\text{ m}$ trong các ODD để đếm số lượng, xác định loài, nguồn gốc (chồi/hạt), chất lượng và phân cấp chiều cao cây non ($H < 0.5\text{ m}$; $0.5 - 1.0\text{ m}$; $1.0 - 2.0\text{ m}$; $> 2.0\text{ m}$).
  • Tam giác đạc phương pháp (Triangulation) & Độ tin cậy: Kết hợp điều tra thực địa với giám sát viễn thám GIS hiện trạng rừng; đối chiếu phân loại thực vật với tài liệu chuẩn (Danh lục thực vật Việt Nam, Khóa định loại thực vật, APG IV); kiểm định độ tin cậy của các mô hình thống kê qua phép thử $\chi^2$ (Chi-square test), kiểm định Kolmogorov-Smirnov ($p > 0.05$) và 999 mô phỏng Monte Carlo.

Data và phân tích

  • Kỹ thuật phân tích cấu trúc lâm phần: Mô hình hóa cấu trúc đường kính ($N/D_{1.3}$) và chiều cao ($N/H_{vn}$) bằng các hàm Weibull 3 tham số, hàm khoảng cách và hàm Meyer trên phần mềm chuyên dụng và R.
  • Kỹ thuật phân tích mô hình điểm không gian: Sử dụng phần mềm Programita và gói lệnh spatstat trong môi trường R để tính toán hàm tương quan cặp $g(r)$, thiết lập dải bao tin cậy 99% (Monte Carlo envelope) để xác định kiểu phân bố (cụm, ngẫu nhiên hay đều) và mối quan hệ giữa các cấp tuổi (cây non - cây sào - cây thành thục).
  • Phân tích đa dạng loài và phát sinh chủng loại: Sử dụng các gói lệnh vegan và picante trong R để tính toán chỉ số Shannon-Wiener ($H'$), Simpson ($D$), Margalef ($d$), Pielou ($J'$), Renyi diversity profiles, đường cong Robbins, chỉ số $MPD, MNTD, NRI, NTI$.
  • Phân tích phân bố độ phong phú loài (SAD): Ứng dụng gói lệnh sads và fitdistrplus trong R để kiểm tra mức độ phù hợp của 6 mô hình lý thuyết, so sánh sự tối ưu thông qua giá trị chênh lệch Akaike ($\Delta AIC$).
  • Định lượng sinh khối trên mặt đất (AGB): Áp dụng phương trình tương quan sinh trắc allometric chuẩn cho rừng nhiệt đới ẩm để ước tính sinh khối và trữ lượng carbon giai đoạn 2015–2020.

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

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

  • Quy luật phân bố $N/D_{1.3}$ và $N/H_{vn}$ tuân thủ chặt chẽ hàm Weibull 3 tham số: Cấu trúc mật độ theo cỡ đường kính biểu hiện đường cong giảm sâu dạng chữ "J ngược", với mật độ cá thể tập trung cực lớn ở các cấp kính nhỏ ($6 - 12\text{ cm}$) và giảm nhanh ở các cấp kính lớn ($> 30\text{ cm}$). Phân bố chiều cao $N/H_{vn}$ có dạng một đỉnh lệch trái, khẳng định lâm phần đang trong giai đoạn hồi sinh mãnh liệt sau xáo trộn với trữ lượng cây non dồi dào.
  • Phát hiện quy luật chuyển dịch không gian từ phân bố cụm sang ngẫu nhiên: Kết quả phân tích hàm $g(r)$ chỉ ra rằng $80.4%$ các loài ưu thế và toàn bộ tầng cây non thể hiện phân bố dạng cụm (aggregated pattern) ở cự ly hẹp ($0 - 10\text{ m}$) do cơ chế phát tán hạn chế hạt giống quanh cây mẹ. Ngược lại, khi chuyển sang giai đoạn cây sào và cây thành thục, phân bố chuyển dịch rõ rệt sang dạng ngẫu nhiên (random pattern) do áp lực tự tỉa thưa và cạnh tranh không gian dinh dưỡng. Tương tác giữa cây non và cây thành thục chủ yếu là tương hỗ hoặc độc lập ở cự ly $< 5\text{ m}$, minh chứng cho cơ chế bảo vệ nhóm loài.
  • Xác lập hiện tượng co cụm phát sinh loài ($NRI > 0$, $NTI > 0$): Tại cả 2 ô định vị 12 và 13, cấu trúc phát sinh chủng loại cho thấy các loài cùng tồn tại có quan hệ họ hàng tiến hóa gần gũi hơn so với kỳ vọng ngẫu nhiên từ nhóm loài khu vực ($p < 0.05$). Điều này chứng minh quá trình chọn lọc môi trường (environmental filtering) - đặc biệt là điều kiện địa hình dốc và khí hậu núi cao - đóng vai trò chủ đạo định hình thành phần loài hơn là cạnh tranh loại trừ.
  • Mô hình trung lập và Log-normal chiếm ưu thế trong phân bố độ phong phú loài (SAD): Kiểm định qua 6 mô hình cho thấy mô hình trung tính (Neutral model với $\Delta AIC$ trung bình $= 1.3$) và Log-normal vượt trội hoàn toàn so với mô hình Niche preemption hay Broken stick ($\Delta AIC > 20$). Điều này chứng tỏ sự phong phú loài trong lâm phần rừng nghèo kiệt chịu tác động lớn của các biến cố phát tán ngẫu nhiên kết hợp với sự tích lũy loài theo bậc tiến hóa.
  • Biến động tích cực về cấu trúc và sinh khối giai đoạn 2015–2020: Trong chu kỳ 5 năm, tỷ lệ cây tái sinh bổ sung (recruitment) đạt mức cao, bù đắp vượt bậc so với tỷ lệ cây chết tự nhiên (tập trung chủ yếu ở cấp kính nhỏ do tự tỉa thưa). Tổng tiết diện ngang ($\sum G/\text{ha}$) và sinh khối trên mặt đất (AGB) ghi nhận mức tăng trưởng dương rõ rệt, chứng minh năng lực phục hồi nội tại (ecological resilience) rất lớn của kiểu rừng LRTX khi được bảo vệ nghiêm ngặt.

Implications đa chiều

  • Đóng góp học thuật: Bổ sung bộ dữ liệu định lượng chuẩn xác và dài hạn về động thái rừng núi cao Việt Nam vào cơ sở dữ liệu sinh thái mạng lưới ô định vị rừng toàn cầu (ForestGEO/CTFS network).
  • Đổi mới phương pháp luận: Thiết lập quy trình mẫu mực tích hợp phân tích điểm không gian 2D, cây phát sinh chủng loại APG IV và mô hình SAD trong nghiên cứu lâm sinh học tại Việt Nam, mở rộng tiềm năng áp dụng cho các hệ sinh thái rừng nhiệt đới khác.
  • Ứng dụng thực tiễn lâm sinh:
    • Xác định cự ly điều chỉnh mật độ trong các biện pháp tỉa thưa nuôi dưỡng rừng: dựa vào bán kính phân bố cụm $5 - 10\text{ m}$ để loại bỏ cây phẩm chất xấu, giảm áp lực cạnh tranh cho cây mục đích.
    • Lựa chọn loài phục hồi rừng: Ưu tiên các loài bản địa ưu thế thuộc họ Thầu dầu (Euphorbiaceae), họ Dẻ (Fagaceae), họ Re (Lauraceae) có khả năng tái sinh hạt tốt và giá trị bảo tồn cao.
  • Khuyến nghị chính sách: Cung cấp cơ sở định lượng để Ban quản lý VQG Phia Oắc - Phia Đén điều chỉnh phân khu chức năng bảo vệ nghiêm ngặt, xây dựng đề án chi trả dịch vụ môi trường rừng gắn với hấp thụ và lưu giữ carbon (C-PFES).

Limitations và Future Research

  • Giới hạn về chu kỳ thời gian và cỡ mẫu: Luận án mới theo dõi biến động cấu trúc trong chu kỳ 5 năm (2015–2020) trên 02 ô định vị sinh thái quốc gia. Mặc dù số liệu chi tiết nhưng chu kỳ 5 năm là tương đối ngắn đối với vòng đời kéo dài hàng trăm năm của rừng tự nhiên nhiệt đới.
  • Giới hạn về dữ liệu vi khí hậu và thổ nhưỡng tầng sâu: Nghiên cứu chưa tích hợp chuỗi số liệu đo đạc liên tục về vi khí hậu dưới tán rừng (nhiệt độ, độ ẩm vi mô, bức xạ tán che) và các chỉ tiêu hóa lý đất tầng sâu theo mùa.
  • Chương trình nghiên cứu tương lai (Future Research Agenda):
    • Kéo dài chu kỳ đo đếm định kỳ 5 năm một lần tại các ODV đến các năm 2025, 2030 để thiết lập chuỗi dữ liệu động thái sinh thái 15–20 năm.
    • Ứng dụng công nghệ LiDAR quét laser 3D từ thiết bị bay không người lái (UAV) để lượng hóa cấu trúc vòm tán đa tầng và sinh khối chính xác ở quy mô cảnh quan.
    • Tích hợp kỹ thuật giải trình tự gen môi trường (eDNA metabarcoding) để đánh giá tương tác giữa hệ vi sinh vật đất và mạng lưới rễ cây rừng trong quá trình tái sinh.
    • Nghiên cứu sâu hơn về chức năng sinh thái (functional traits) như diện tích lá chuyên biệt (SLA), tỷ trọng gỗ, khả năng dẫn truyền nước để liên kết cấu trúc không gian với khả năng thích ứng biến đổi khí hậu.

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

  • Tác động học thuật: Luận án tạo ra bước ngoặt trong nghiên cứu lâm học định lượng tại Việt Nam, mở ra hướng nghiên cứu giao thoa giữa lâm sinh học, toán sinh học và sinh thái học phân tử, dự kiến mang lại các công bố quốc tế uy tín trên các tạp chí thuộc hệ thống ISI/Scopus (Q1/Q2).
  • Chuyển đổi quản trị rừng: Cung cấp phương pháp luận chuẩn xác cho lực lượng kiểm lâm và cán bộ kỹ thuật lâm nghiệp trong việc đánh giá thực trạng tài nguyên rừng mà không chỉ dựa trên các chỉ tiêu trữ lượng giản đơn.
  • Tác động chính sách: Đóng góp luận cứ khoa học trực tiếp cho Bộ Nông nghiệp và PTNT, UBND tỉnh Cao Bằng trong việc hoàn thiện Phương án Quản lý rừng bền vững VQG Phia Oắc - Phia Đén đến năm 2030, hỗ trợ thực thi Chiến lược phát triển lâm nghiệp Việt Nam giai đoạn 2021–2030 và cam kết đạt phát thải ròng bằng "0" (Net Zero) vào năm 2050.

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

  • Nghiên cứu sinh và Giảng viên Lâm nghiệp: Thụ hưởng khung phương pháp luận tiên tiến, bộ số liệu thực nghiệm chuẩn tắc về cấu trúc không gian, phát sinh loài và mô hình SAD để phục vụ giảng dạy cao học và phát triển đề tài nghiên cứu chuyên sâu.
  • Ban quản lý Vườn quốc gia và Khu bảo tồn thiên nhiên: Sở hữu tài liệu khoa học làm căn cứ định vị các khu vực cần xúc tiến tái sinh tự nhiên, lựa chọn tập đoàn cây bản địa trong phục hồi rừng và tối ưu hóa chi phí tuần tra bảo vệ rừng.
  • Các nhà hoạch định chính sách lâm nghiệp: Có cơ sở dữ liệu xác thực để xây dựng hạn mức hấp thụ carbon rừng, định giá dịch vụ hệ sinh thái và ban hành các quy trình kỹ thuật quốc gia về phục hồi rừng tự nhiên nghèo kiệt.

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ì?
    Đóng góp độc đáo nhất là việc mở rộng Thuyết trung tính (Hubbell) và Thuyết ổ sinh thái (Tilman) thông qua việc chứng minh sự tồn tại của hiện tượng co cụm phát sinh loài ($NRI > 0, NTI > 0$) kết hợp với cơ chế phân bố cụm do phát tán hạn chế ở rừng LRTX nghèo kiệt vùng núi cao Phia Oắc - Phia Đén. Nghiên cứu xác lập rằng trong điều kiện lập địa khắc nghiệt sau suy thoái, chọn lọc môi trường đóng vai trò chi phối tuyệt đối cấu trúc quần xã.

  2. Phương pháp nghiên cứu có điểm gì cải tiến so với các công trình trước đây?
    So với các nghiên cứu lâm học truyền thống tại Việt Nam vốn chỉ dừng ở hàm phân bố tĩnh (như Cao Thị Thu Hiền, 2013; Phùng Đình Trung, 2019), luận án đã nâng cấp toàn diện bằng cách tích hợp phân tích điểm không gian 2D đa biến (hàm $g(r)$, kiểm định CSR qua mô phỏng Monte Carlo), phân tích cấu trúc phát sinh chủng loại trên nền siêu cây Megatree APG IV, và kiểm định 6 mô hình phân bố độ phong phú loài (SAD) bằng tiêu chí thông tin Akaike (AIC).

  3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?
    Phát hiện bất ngờ là mật độ tương tác tương hỗ giữa tầng cây non và cây thành thục chiếm ưu thế ở cự ly $< 5\text{ m}$ thay vì quan hệ cạnh tranh gay gắt. Điều này khẳng định cơ chế bảo vệ nhóm loài (herd protection) và hiệu ứng che chở vi khí hậu của cây tầng trên giữ vai trò quyết định giúp cây con vượt qua các điều kiện khắc nghiệt của mùa đông vùng núi cao Phia Oắc.

  4. Luận án có cung cấp giao thức tái lặp nghiên cứu (replication protocol) không?
    Có. Toàn bộ thiết kế hệ thống ô định vị sinh thái (ODV), kích thước ô nghiên cứu ($20\text{ m} \times 20\text{ m}$), ô dạng bản ($5\text{ m} \times 5\text{ m}$), thuật toán xử lý không gian trên Programita/R, mã lệnh phân tích phát sinh loài picante và các tiêu chuẩn kiểm định AIC đều được mô tả chi tiết, minh bạch trong chương phương pháp và phụ lục, cho phép các nhà khoa học độc lập tái lập hoàn toàn trên các hệ sinh thái khác.

  5. Chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?
    Chương trình nghiên cứu 10 năm (2025–2035) tập trung vào 3 trọng tâm: (i) Tiếp tục quan trắc định kỳ 5 năm trên 02 ODV để đánh giá động thái tự tỉa thưa và tăng trưởng AGB dài hạn; (ii) Tích hợp công nghệ cảm biến vi khí hậu tự động và quét 3D LiDAR vòm tán; (iii) Nghiên cứu cơ chế trao đổi chất và đa dạng chức năng (functional diversity) của các loài cây bản địa quý hiếm để ứng phó với biến đổi khí hậu toàn cầu.

Kết luận

  1. Luận án đã giải mã toàn diện các đặc trưng cấu trúc phân bố $N/D_{1.3}$ và $N/H_{vn}$ của rừng LRTX nghèo và nghèo kiệt tại VQG Phia Oắc - Phia Đén, chứng minh sự phù hợp tối ưu của hàm phân bố Weibull 3 tham số và quy luật giảm số mũ của rừng tự nhiên đang hồi phục.
  2. Xác lập bằng chứng khoa học vững chắc về quy luật phân bố không gian: chuyển dịch từ phân bố dạng cụm ở cây non sang phân bố ngẫu nhiên ở cây thành thục do áp lực tự tỉa thưa và hạn chế phát tán.
  3. Đột phá trong việc chứng minh hiện tượng co cụm phát sinh loài ($NRI > 0, NTI > 0$), khẳng định vai trò quyết định của chọn lọc môi trường trong cấu trúc quần xã thực vật núi cao.
  4. Chứng minh mô hình trung lập và Log-normal mô tả chính xác phân bố độ phong phú loài (SAD) với giá trị $\Delta AIC < 2$, chỉ ra tính đa cơ chế trong duy trì đa dạng sinh thái rừng thứ sinh.
  5. Lượng hóa cụ thể động thái tích cực giai đoạn 2015–2020: sinh khối AGB và tổng tiết diện ngang thân cây tăng trưởng ổn định, tỷ lệ tái sinh bổ sung vượt trội tỷ lệ cây chết.
  6. Đề xuất hệ thống giải pháp kỹ thuật lâm sinh có luận cứ khoa học chuẩn xác, từ nuôi dưỡng rừng, điều chỉnh cự ly không gian đến bảo tồn nguồn gen các loài cây gỗ bản địa quý hiếm, đóng góp thiết thực cho sự nghiệp quản lý rừng bền vững và bảo tồn đa dạng sinh học quốc gia.