Tổng quan về luận án

Nghiên cứu động thái cấu trúc hệ sinh thái rừng lá rộng thường xanh đóng vai trò nền tảng trong khoa học lâm nghiệp và sinh thái học nhiệt đới, đặc biệt trong bối cảnh diện tích rừng tự nhiên Việt Nam dù đạt 14,6 triệu ha năm 2019 với độ che phủ 41,89% nhưng chất lượng trữ lượng và tính đa dạng sinh học vẫn suy giảm nghiêm trọng. Luận án tiến sĩ sinh học chuyên ngành Sinh thái học (Mã số: 9 42 01 20) của nghiên cứu sinh Nguyễn Tiến Dũng, thực hiện tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam dưới sự hướng dẫn của TS. Trần Văn Con và TS. Nguyễn Văn Sinh (2020), mang tính tiên phong khi chuyển đổi phương thức nghiên cứu từ tĩnh sang động thông qua việc theo dõi thực nghiệm dài hạn trên hệ thống ô tiêu chuẩn định vị (OTCĐV).

Research gap cụ thể mà luận án giải quyết xuất phát từ thực trạng phần lớn các công trình sinh thái rừng tại Việt Nam trước đây chủ yếu áp dụng phương pháp "tiếp cận chuỗi thời gian" (lấy không gian thay thế thời gian) trên các ô tiêu chuẩn tạm thời, vốn bộc lộ nhiều sai số do không phản ánh được tính kế tục cơ chế sinh thái nội tại. Thiếu hụt nghiêm trọng các bộ dữ liệu quan trắc định vị liên tục về tốc độ chuyển cấp kính ($O_k$), tỷ lệ tử vong tự nhiên ($M_k$) và lượng tái sinh bổ sung ($R$) vào tầng cây cao ($D_{1.3} \ge 10\text{ cm}$) của các trạng thái rừng tự nhiên phục hồi sau khai thác và rừng nguyên sinh trên núi đá vôi vùng Tây Bắc.

Luận án đặt ra ba câu hỏi nghiên cứu và giả thuyết tương ứng:

  1. Câu hỏi 1 (RQ1): Các nhân tố cấu trúc hình thái ($N/D_{1.3}$, $N/H_{vn}$), cấu trúc sinh thái ($IVI%$, $HL_1$) và tính đa dạng loài biến đổi theo quy luật nào qua chu kỳ theo dõi 11 năm? Giả thuyết 1 (H1): Cấu trúc tổ thành và phân bố số cây theo cấp đường kính ($N/D_{1.3}$) duy trì dạng giảm liên tục (chữ J ngược hoặc hàm số mũ), phản ánh trạng thái cân bằng động của hệ sinh thái rừng lá rộng thường xanh.
  2. Câu hỏi 2 (RQ2): Có tồn tại mối tương quan phi tuyến định lượng giữa kích thước cây ($D_{1.3}$) với số cây chuyển cấp ($O_k$) và số cây chết ($M_k$) hay không? Giả thuyết 2 (H2): Mối quan hệ chuyển cấp $O_k - D_k$ và chết $M_k - D_k$ tuân theo phương trình hàm bậc 3 với độ tin cậy thống kê cao ($p < 0,05$).
  3. Câu hỏi 3 (RQ3): Có thể tích hợp các tham số động thái thực nghiệm vào phần mềm chuyên dụng để dự đoán cấu trúc lâm phần trong trung và dài hạn? Giả thuyết 3 (H3): Mô hình hóa động lực học hệ thống (System Dynamics) thông qua phần mềm MM&S cho phép dự báo chính xác động thái $N/D_{1.3}$ và $G/D_{1.3}$ ở các bước thời gian 11 năm đến 55 năm.

Khung lý thuyết của luận án tích hợp: (i) Lý thuyết cấu trúc phân tầng thẳng đứng của Thái Văn Trừng (1978, 1999); (ii) Lý thuyết động thái lỗ trống tàn che (Gap-phase dynamics) của Van Steenis (1956); (iii) Lý thuyết entropy tổng quát Rényi (1961); và (iv) Lý thuyết mô hình hóa hệ thống động lực rừng của Hartmut Bossel & Holger Krieger (1991). Luận án được thực hiện trên 6 OTCĐV quy chuẩn (diện tích mỗi ô từ 1 ha chia thành các ô phụ $20\text{ m} \times 20\text{ m}$) trong khung thời gian 11 năm (2007–2018) tại Khu bảo tồn thiên nhiên Hang Kia – Pà Cò (Hòa Bình) và Vườn quốc gia Xuân Sơn (Phú Thọ).

Literature Review và Positioning

Các công trình nghiên cứu động thái rừng nhiệt đới trên thế giới đã hình thành nên những trường phái lý thuyết kinh điển song hành với những tranh luận khoa học sâu sắc. Điển hình là cuộc tranh luận giữa hai quan điểm tái sinh rừng: Lý luận bức khảm tái sinh (Mosaic / Cyclic Regeneration Theory) của Aubréville (1938) tại rừng mưa nhiệt đới châu Phi cho rằng "tổ hợp loài cây của rừng mưa đều không cố định trong không gian và theo thời gian; không có một tổ hợp của loài cây nào có thể đạt được thế cân bằng sinh thái với hoàn cảnh một cách vĩnh viễn và ổn định", cây mẹ tầng ưu thế thường vắng bóng cây con ở tầng dưới; đối lập với quan sát của Davis & Richards (1933), Beard (1964), Sun (1960) và Rollet (1969) tại rừng nhiệt đới Nam Mỹ khẳng định hiện tượng tái sinh tại chỗ liên tục, nơi các loài chiếm ưu thế ở tầng cao luôn hiện diện với mật độ tương ứng ở các cấp thể tích nhỏ dưới tán.

Van Steenis (1956) làm sâu sắc thêm sinh thái học phục hồi khi chỉ ra quy luật tái sinh vệt và vai trò của các loài cây tiên phong ưa sáng – được mệnh danh là "loài liền vết sẹo" (scar-healing species) – xuất hiện tại các lỗ trống tàn che trước khi bị thay thế bởi các loài chịu bóng trường thọ. Trong nước, Thái Văn Trừng (1978, 1999) xác định ánh sáng là nhân tố sinh thái khống chế cấu trúc 5 tầng của rừng mưa mùa, trong khi Trần Ngũ Phương (1970) chứng minh diễn thế thứ sinh rừng Lim xanh (Erythrophleum fordii) luôn gắn liền với điều kiện lập địa địa phương. Trần Văn Con (2009) và Nguyễn Đắc Triển (2015) bước đầu ghi nhận hiện tượng tích tụ số loài theo cấp kính tăng dần tại miền Bắc Việt Nam.

Luận án định vị chính xác khoảng trống học thuật bằng việc đối chiếu trực tiếp với hai nghiên cứu quốc tế tiêu biểu:

  1. So sánh với nghiên cứu của Jeanine Maria Felfili (1997) tại rừng hành lang ven sông Gama (Brazil) với 151 ô định vị ($10\text{ m} \times 20\text{ m}$), theo dõi cây có chu vi $G_{bh} \ge 31\text{ cm}$ ($D_{1.3} > 10\text{ cm}$) ghi nhận tỷ lệ chết chung 3,5%/năm và tỷ lệ bổ sung 2,7%/năm. Luận án xác định tại Hang Kia – Pà Cò và Xuân Sơn, tỷ lệ chết tự nhiên hàng năm dao động thấp hơn đáng kể ($0,5% - 1,8%/\text{năm}$), tương đồng với cơ chế tự tỉa thưa ở các trạng thái IIIA3, IIIB và IV trên núi đá vôi Ferasol/Alisol.
  2. So sánh với nghiên cứu của Koichi Takahashi et al. (2003) tại Nhật Bản trên ô định vị 1 ha, nơi loài ưu thế tầng tán Quercus crispula có tỷ lệ bổ sung gần như bằng 0 do lớp tre lùn dưới tán khống chế, trong khi các loài tầng phụ (Acer mono, Abies sachalinensis) tái sinh mạnh. Ngược lại, dữ liệu của luận án cho thấy tại rừng lá rộng thường xanh phía Bắc Việt Nam, các loài ưu thế tầng cao như Re gừng (Cinnamomum sp.), Dẻ gai (Castanopsis sp.), Kháo gai (Machilus sp.) vẫn duy trì khả năng tái sinh bổ sung ổn định vào tầng cây cao khi xuất hiện lỗ trống vi khí hậu thuận lợi.
+-----------------------------------------------------------------------------------+
| KHUNG ĐỊNH VỊ HỌC THUẬT (THEORETICAL & METHODOLOGICAL POSITIONING)                |
+-----------------------------------------------------------------------------------+
|  Trường phái Cổ điển (Aubréville 1938, Richards 1933, Van Steenis 1956)           |
|  - Tranh luận: Bức khảm tuần hoàn vs Tái sinh liên tục tại chỗ                    |
|  - Cơ chế: Động thái lỗ trống tàn che (Gap-phase dynamics)                        |
|  Sinh thái học Cấu trúc Rừng VN (Thái Văn Trừng 1978, Trần Ngũ Phương 1970)       |
|  - Cấu trúc phân tầng A1-A2-A3; Hệ thống phân loại trạng thái Loeschau            |
|  - Hạn chế lịch sử: Phương pháp giải tích cây, tĩnh học, thay thế không gian      |
|  ĐÓNG GÓP TIÊN PHONG CỦA LUẬN ÁN NGUYỄN TIẾN DŨNG (2020)                          |
|  - Quan trắc động thái thực nghiệm 11 năm (2007-2018) trên 6 OTCĐV 1 ha           |
|  - Định lượng hóa: Tái sinh bổ sung (R), Chuyển cấp (Ok), Tự tỉa thưa (Mk)        |
|  - Mô hình hóa toán học bậc 3 & Mô phỏng System Dynamics bằng phần mềm MM&S       |
+-----------------------------------------------------------------------------------+

Đó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à kiểm chứng thực nghiệm học thuyết cân bằng động trong quần xã thực vật nhiệt đới qua các luận điểm cốt lõi:

  • Mở rộng lý thuyết phân bố $N/D_{1.3}$: Chứng minh rằng phân bố số cây theo cấp đường kính của quần xã cây gỗ rừng tự nhiên duy trì sự ổn định cấu trúc dạng hàm số mũ giảm ($N = N_0 \cdot e^{-k \cdot D}$) hoặc dạng chữ J ngược qua hơn một thập kỷ, khẳng định khả năng tự điều chỉnh (homeostasis) mạnh mẽ của hệ sinh thái rừng nhiệt đới khi không có tác động nhân vi.
  • Phát triển lý thuyết tích tụ loài: Củng cố luận điểm của Trần Văn Con (2009) rằng số lượng loài có xu hướng tăng dần từ lớp cây tái sinh lên tầng cây cao nhờ sự hội tụ của nhiều đợt phát tán hạt giống và tái sinh lỗ trống qua thời gian dài, bác bỏ giả thuyết cho rằng tầng cây cao chỉ phản ánh tổ thành đơn điệu của một thế hệ duy nhất.
  • Chuyển dịch Paradigm (Paradigm Shift): Chuyển hóa phương pháp luận nghiên cứu cấu trúc rừng tại Việt Nam từ mô tả định tính, tĩnh học sang định lượng động thái học (Quantitative System Dynamics), xác lập các phương trình phi tuyến mô tả tốc độ đào thải và thay thế cá thể theo cấp kính.

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

Khung phân tích của luận án kết hợp đa chiều ba hệ thống lý thuyết:

  1. Hệ thống chỉ số sinh thái quần xã: Sử dụng chỉ số mức độ quan trọng $IVI%$ tổng hợp từ 3 cấu phần: tần số xuất hiện tương đối ($F%$), mật độ tương đối ($N%$) và tiết diện ngang tương đối ($G%$): $$IVI% = \frac{F% + N% + G%}{3}$$ Cùng hệ số hỗn loài $HL_1 = S/N$ để xác định cấu trúc ưu thế sinh thái.
  2. Hệ entropy tổng quát Rényi ($H_\alpha$): $$H_\alpha = \frac{\ln\left(\sum_{i=1}^{s} p_i^\alpha\right)}{1 - \alpha}$$ Khảo sát liên tục với $\alpha \in {0,1; 0,25; 0,5; 1; 2; 3; 4; 5; 6; 7; 8; \infty}$, cho phép đánh giá toàn diện từ độ giàu loài ($\alpha = 0$, $H_0 = \ln S$), chỉ số Shannon-Wiener ($\alpha = 1$, $H_1 = -\sum p_i \ln p_i$), chỉ số ưu thế Simpson ($\alpha = 2$, $H_2 = \ln(1/D)$) đến độ ưu thế của loài chiếm ưu thế nhất ($\alpha = \infty$, $H_\infty = \ln(1/p_{max})$).
  3. Khung động lực học chuyển cấp và mô phỏng số MM&S: Tiếp cận hệ thống toàn lâm phần (Stand-level system dynamics), mô hình hóa mối quan hệ giữa các biến trạng thái (số cây ở cấp kính $k$, $N_k$) với các biến dòng chảy: dòng vào là tái sinh bổ sung $R$ và chuyển cấp vào $I_k$, dòng ra là chuyển cấp ra $O_k$ và cây chết tự nhiên $M_k$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế trên nền tảng triết học hiện thực phản biện (Critical Realism) kết hợp thực chứng định lượng (Positivism). Luận án áp dụng phương pháp nghiên cứu trên hệ thống ô tiêu chuẩn định vị dài hạn (Permanent Sample Plot - PSP) theo tiêu chuẩn của Liên đoàn các Tổ chức Nghiên cứu Lâm nghiệp Quốc tế (IUFRO).

Thiết kế mẫu gồm 6 OTCĐV kích thước lớn ($100\text{ m} \times 100\text{ m} = 1\text{ ha}$ mỗi ô), được phân chia chi tiết thành 25 ô phụ ($20\text{ m} \times 20\text{ m}$) để định vị không gian tọa độ từng cá thể cây gỗ.

  • Khu bảo tồn thiên nhiên Hang Kia – Pà Cò (Hòa Bình):
    • HB01: $20^\circ 43'40''\text{N}, 104^\circ 56'27''\text{E}$; Độ cao: 1060 m; Độ dốc: $55^\circ$; Hướng: Đông Nam; Đất: Ferasol/Alisol trên núi đá vôi; Trạng thái rừng: IIIA3.
    • HB03: $20^\circ 44'57''\text{N}, 104^\circ 54'37''\text{E}$; Độ cao: 1230 m; Độ dốc: $35^\circ$; Hướng: Đông Nam; Đất: Ferasol/Alisol; Trạng thái rừng: IIIA3.
    • HB06: $20^\circ 44'15,28''\text{N}, 104^\circ 56'42,97''\text{E}$; Độ cao: 980 m; Độ dốc: $35^\circ$; Hướng: Đông Nam; Đất: Ferasol/Alisol; Trạng thái rừng: IV (rừng nguyên sinh ít tác động).
  • Vườn Quốc gia Xuân Sơn (Phú Thọ):
    • XS01: $21^\circ 06'36,2''\text{N}, 104^\circ 57'47,6''\text{E}$; Độ cao: 304 m; Độ dốc: $25^\circ$; Hướng: Đông Nam; Đất: Ferasol/Alisol trên đồi núi thấp; Trạng thái rừng: IIIB.
    • XS02: $21^\circ 07'55,9''\text{N}, 104^\circ 57'04,5''\text{E}$; Độ cao: 540 m; Độ dốc: $30^\circ$; Hướng: Đông Nam; Đất: Ferasol/Alisol; Trạng thái rừng: IIIB.
    • XS03: $21^\circ 08'09,8''\text{N}, 104^\circ 55'39,5''\text{E}$; Độ cao: 455 m; Độ dốc: $20^\circ$; Hướng: Đông Nam; Đất: Ferasol/Alisol; Trạng thái rừng: IIIB.

Quy trình nghiên cứu rigorous

Quy trình thu thập số liệu tuân thủ nghiêm ngặt hai giai đoạn đo đếm chuẩn hóa:

  1. Đo đếm nền tảng (Năm 2007): Thu thập toàn bộ các cá thể cây gỗ tầng cao có đường kính ngang ngực $D_{1.3} \ge 10\text{ cm}$ (đo chu vi $G_{bh}$ bằng thước vanh chuyên dụng chia vạch đến milimet, chuyển đổi sang $D_{1.3}$). Đo chiều cao vút ngọn ($H_{vn}$) bằng thước đo cao siêu âm Vertex IV (độ chính xác decimet). Xác định tên loài thực vật qua 2 bước: giám định tại chỗ và thu thập tiêu bản (lá, hoa, quả) đối với các mẫu chưa rõ danh pháp, giám định đối chiếu với các tài liệu chuẩn quốc gia (Cây cỏ Việt Nam của Phạm Hoàng Hộ, Thực vật chí Việt Nam, Tên cây rừng Việt Nam).
  2. Đo đếm lặp lại và quan trắc động thái (Năm 2018): Sau 11 năm, tiến hành kiểm kê toàn bộ cây cũ: xác định chính xác số hiệu cây còn sống ($N_s$), đo đếm lại $D_{1.3}$ và $H_{vn}$ để xác định sinh trưởng và chuyển cấp ($O_k$); kiểm kê cây chết tự nhiên ($M_k$, phân tích nguyên nhân đổ gãy, tàn lụi); thống kê và đánh số mới toàn bộ cây tái sinh bổ sung vào tầng cây cao ($R$, những cây đạt $D_{1.3} \ge 10\text{ cm}$ xuất hiện mới). Khảo sát lớp cây bụi thảm tươi bằng 5 ô dạng bản $5\text{ m} \times 5\text{ m}$ tại mỗi OTCĐV.

Data và phân tích

Phân tích dữ liệu ứng dụng các kỹ thuật toán học thống kê và mô hình hóa sinh thái:

  • Xử lý đa dạng loài và cấu trúc: Thiết lập các bảng tần số phân bố $N/D_{1.3}$ (cự ly cấp kính 4–5 cm) và phân bố $N/H_{vn}$ (cự ly cấp chiều cao 2 m). Tính toán các chỉ số Rényi, Shannon-Wiener và Simpson mẫu: $$D = 1 - \sum_{i=1}^s \frac{n_i(n_i - 1)}{N(N - 1)}$$
  • Mô hình hóa phi tuyến chuyển cấp và cây chết: Sử dụng phương pháp hồi quy đa thức và phân tích phương sai (ANOVA) để kiểm định sự tồn tại của phương trình tương quan bậc 3 dạng: $$O_k = a \cdot D_k^3 + b \cdot D_k^2 + c \cdot D_k + d$$ $$M_k = a' \cdot D_k^3 + b' \cdot D_k^2 + c' \cdot D_k + d'$$ Kiểm tra hệ số xác định $R^2$, giá trị kiểm định Fisher ($F$) và mức ý nghĩa thống kê $p < 0,05$.
  • Phần mềm và công cụ ứng dụng: Xử lý thống kê trên R và Excel; trực quan hóa không gian địa lý trên Google Earth và GIS; mô phỏng động thái cấu trúc lâm phần trong tương lai bằng phần mềm mô hình hóa hệ thống MM&S (Model Maker & Solver), thiết lập mô hình chạy 5 bước thời gian (mỗi bước $t = 1 = 11\text{ năm}$, tổng cộng 55 năm dự báo).

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

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

  1. Sự ổn định vượt bậc của phân bố cấu trúc $N/D_{1.3}$: Qua 11 năm quan sát thực chứng, phân bố số cây theo cấp đường kính tại cả 6 OTCĐV thuộc 3 trạng thái rừng (IIIA3, IIIB, IV) đều duy trì cấu trúc chữ J ngược ổn định. Quá trình tái sinh bổ sung ($R$) bù đắp thỏa đáng cho số cây chuyển cấp ($O_k$) và số cây chết ($M_k$) ở các cấp kính nhỏ ($10 - 20\text{ cm}$), giúp cấu trúc lâm phần không bị phá vỡ.
  2. Quy luật chuyển cấp phi tuyến tuân theo hàm bậc 3: Số lượng cây chuyển ra khỏi cấp kính ($O_k$) đạt cực đại ở các cấp kính nhỏ ($12 - 20\text{ cm}$), sau đó giảm dần và tiệm cận mức ổn định ở các cấp kính lớn ($D_{1.3} > 40\text{ cm}$). Phân tích phương sai chứng minh phương trình bậc 3 mô phỏng quan hệ $O_k - D_k$ và $M_k - D_k$ có ý nghĩa thống kê rất cao ($p < 0,01$), thể hiện rõ ràng cơ chế cạnh tranh ánh sáng gay gắt ở tầng tán dưới.
  3. Quy luật tử vong tự nhiên ($M_k$) mang tính đặc thù lập địa: Tỷ lệ chết tập trung chủ yếu ở các cá thể cây nhỏ bị chèn ép tán ($D_{1.3} < 16\text{ cm}$) do hiện tượng tự tỉa thưa sinh học, trong khi ở cấp kính lớn ($D_{1.3} > 50\text{ cm}$), số lượng cây chết rất thấp nhưng khi xảy ra thường do các biến cố đổ gãy cơ học trên sườn núi đá dốc ($35^\circ - 55^\circ$).
  4. Hiện tượng tích tụ tính đa dạng loài theo thời gian: Phân tích entropy Rényi chứng minh chỉ số đa dạng loài ở trạng thái rừng nguyên sinh IV (OTC HB06) duy trì giá trị cao nhất và ổn định nhất ($H_\alpha$ cao ở mọi bậc $\alpha$), trong khi các trạng thái rừng phục hồi IIIA3 và IIIB có xu hướng tiệm cận dần về độ giàu loài của rừng nguyên sinh nhờ sự xuất hiện của các loài cây chịu bóng trung gian.
+------------------------------------------------------------------------------------+
| SƠ ĐỒ DÒNG ĐỘNG HỌC LÂM PHẦN TẦNG CÂY CAO (SYSTEM DYNAMICS FLOW DIAGRAM)           |
+------------------------------------------------------------------------------------+
|                                                                                    |
|                                                     (Chuyển cấp O1)                |
|                                                     (Chuyển cấp O2)                |
|                                                           ...                      |
|                                                                                    |
+------------------------------------------------------------------------------------+

Implications đa chiều

  • Về lý luận sinh thái: Cung cấp bằng chứng thực nghiệm vững chắc chứng minh tính bền vững nội tại của hệ sinh thái rừng lá rộng thường xanh nhiệt đới gió mùa trên núi đá vôi, củng cố học thuyết diễn thế hồi phục của sinh thái học Việt Nam.
  • Về phương pháp luận: Chuẩn hóa quy trình quan trắc định vị 11 năm và phương pháp ứng dụng phần mềm MM&S trong mô phỏng động thái cấu trúc lâm phần, có thể chuyển giao áp dụng cho các hệ sinh thái rừng khộp, rừng ngập mặn và rừng rụng lá tại Đông Nam Á.
  • Về quản lý và thực tiễn lâm sinh: Kết quả mô phỏng cấu trúc $N/D_{1.3}$ và $G/D_{1.3}$ cung cấp căn cứ định lượng để xây dựng các "mẫu cấu trúc định hướng" (Target Stand Models). Cho phép các Ban quản lý rừng xác định chính xác cường độ tỉa thưa nuôi dưỡng rừng (chỉ tác động vào các cấp kính có mật độ dư thừa so với đường cong lý thuyết), hạn chế tối đa việc can thiệp vào các loài cây ưu thế sinh thái.
  • Về chính sách bảo tồn: Cung cấp cơ sở khoa học để đề xuất chính sách phân vùng bảo vệ nghiêm ngặt đối với các diện tích rừng trạng thái IV và IIIB tại Hang Kia – Pà Cò và Xuân Sơn, định hướng chi trả dịch vụ môi trường rừng gắn liền với năng lực hấp thụ carbon dựa trên sinh khối tích lũy qua tiết diện ngang ($\sum G$).

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  1. Quy mô không gian của mạng lưới OTCĐV: Mặc dù 6 ô tiêu chuẩn định vị 1 ha được thiết lập công phu, cỡ mẫu vẫn còn giới hạn trong phạm vi 2 khu vực địa lý đại diện (Hòa Bình và Phú Thọ), chưa bao phủ toàn bộ các kiểu phụ thổ nhưỡng của đai rừng nhiệt đới gió mùa toàn miền Bắc.
  2. Chu kỳ quan sát thực nghiệm: Chu kỳ 11 năm (2007–2018) phản ánh rõ nét các biến động chuyển cấp và chết tự nhiên ở tầng cây cao, nhưng đối với một số loài cây gỗ siêu đại thụ có tuổi thọ hàng trăm năm, chu kỳ này cần được kéo dài hơn để ghi nhận đầy đủ chu kỳ sống (life-cycle dynamics).
  3. Biến số ngoại cảnh vi khí hậu: Nghiên cứu chưa tích hợp các trạm cảm biến vi khí hậu đo đạc liên tục cường độ bức xạ quang hợp (PAR), độ ẩm đất theo thời gian thực tại các tầng tán.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Duy trì đo đếm chu kỳ 3 (2029) trên hệ thống 6 OTCĐV hiện hữu để thiết lập chuỗi dữ liệu 22 năm.
  • Ứng dụng công nghệ viễn thám LiDAR quét mặt đất (Terrestrial Laser Scanning - TLS) và ảnh chụp Drone phân giải cao để giải mã cấu trúc không gian 3D tầng tán.
  • Mở rộng mô hình toán học tích hợp biến đổi khí hậu (nhiệt độ cực trị, lượng mưa mùa khô) vào phần mềm MM&S để mô phỏng khả năng chống chịu và dự trữ carbon của rừng.

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

Luận án tạo ra tác động học thuật sâu sắc với việc công bố các bộ dữ liệu gốc chuẩn xác, đóng góp vào mạng lưới quan sát sinh thái học rừng nhiệt đới quốc tế (như ForestGEO / CTFS). Kết quả nghiên cứu được chuyển giao trực tiếp cho Ban quản lý Khu BTTN Hang Kia – Pà Cò và Vườn quốc gia Xuân Sơn để hoàn thiện Đề án Quản lý rừng bền vững giai đoạn 2021–2030.

Về mặt xã hội và môi trường, các giải pháp kỹ thuật nuôi dưỡng và bảo vệ rừng đề xuất từ luận án giúp duy trì bể chứa carbon ổn định, bảo vệ nguồn nước ngầm đầu nguồn cho các lưu vực sông Đà và sông Hồng, đồng thời tạo tiền đề khoa học vững chắc cho việc triển khai thị trường tín chỉ carbon rừng tại các cộng đồng dân tộc bản địa Mường, Dao và H'Mông sinh sống tại vùng đệm.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận thiết lập ô tiêu chuẩn định vị IUFRO chuẩn hóa, công thức toán học phân tích entropy Rényi và kỹ thuật mô phỏng phần mềm MM&S.
  • Các nhà sinh thái học và lâm học: Kế thừa bộ số liệu thực chứng 11 năm về tốc độ sinh trưởng đường kính, tỷ lệ chuyển cấp và tỷ lệ chết tự nhiên của rừng nhiệt đới gió mùa.
  • Ban quản lý các Vườn quốc gia và Khu bảo tồn thiên nhiên: Sử dụng các mẫu cấu trúc định hướng để giám sát, đánh giá diễn thế phục hồi và tổ chức tuần tra bảo tồn nghiêm ngặt.
  • Nhà hoạch định chính sách Lâm nghiệp (Cục Lâm nghiệp, Cục Kiểm lâm): Ứng dụng các dẫn liệu định lượng để xây dựng quy trình kỹ thuật quốc gia về phục hồi rừng tự nhiên nghèo kiệt và hướng dẫn quản lý rừng bền 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ì? Đó là việc kiểm chứng và định lượng hóa thành công phương trình chuyển cấp phi tuyến bậc 3 kết hợp mô hình động lực học hệ thống MM&S trên số liệu quan sát thực chứng 11 năm, mở rộng lý thuyết cấu trúc rừng nhiệt đới của Thái Văn Trừng từ dạng mô tả tĩnh sang mô hình toán học dự báo biến đổi cấu trúc theo thời gian.

  2. Điểm đổi mới phương pháp luận so với các nghiên cứu trước đây tại Việt Nam? Khác với phương pháp truyền thống lấy không gian thay thế thời gian của các tác giả trước thập niên 2000, luận án thiết lập phương pháp quan trắc định vị lặp lại trên 6 OTCĐV quy chuẩn 1 ha, áp dụng hệ thống phân tích entropy đa quy mô Rényi ($H_\alpha$) thay vì chỉ dùng các chỉ số đơn lẻ Shannon hay Simpson.

  3. Phát hiện gây ngạc nhiên nhất từ dữ liệu thực nghiệm là gì? Đó là sự bảo toàn gần như tuyệt đối của đường cong phân bố số cây theo cấp kính $N/D_{1.3}$ (dạng hàm mũ / chữ J ngược) tại các trạng thái rừng qua 11 năm. Bất chấp tỷ lệ đào thải do cạnh tranh sinh học ở các cấp kính nhỏ, dòng tái sinh bổ sung ($R$) tự nhiên đã tự động cân bằng chính xác với tỷ lệ tử vong và chuyển cấp, phản ánh năng lực tự phục hồi cực kỳ bền vững của hệ sinh thái rừng lá rộng thường xanh bản địa.

  4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) rõ ràng không? Có. Luận án mô tả chi tiết tọa độ GPS, độ cao, hướng dốc, quy cách chia ô phụ $20\text{ m} \times 20\text{ m}$, quy trình đánh mã số cây, phương pháp đo vanh chu vi milimet, đo chiều cao siêu âm Vertex và thuật toán nạp biến trạng thái vào phần mềm MM&S.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Tiếp tục đo đếm chu kỳ 3 (2029) để hoàn thiện chuỗi dữ liệu 22 năm, tích hợp mô hình phân tích không gian điểm Ripley's K(t) và đo đạc sinh khối tán lá bằng viễn thám LiDAR nhằm tối ưu hóa mô hình dự báo hấp thụ carbon.

Kết luận

Luận án của nghiên cứu sinh Nguyễn Tiến Dũng đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp cốt lõi:

  1. Xác lập bộ cơ sở dữ liệu thực chứng 11 năm (2007–2018) về cấu trúc và động thái của 6 quần xã thực vật rừng lá rộng thường xanh tại Hang Kia – Pà Cò và Xuân Sơn.
  2. Định lượng hóa thành công các quy luật động thái tái sinh bổ sung ($R$), chuyển cấp ($O_k$) và chết tự nhiên ($M_k$) của tầng cây cao qua hệ phương trình hồi quy bậc 3 có ý nghĩa thống kê cao.
  3. Ứng dụng thành công phần mềm chuyên dụng MM&S trong mô phỏng dự báo động thái cấu trúc lâm phần $N/D_{1.3}$ và $G/D_{1.3}$ ở các bước thời gian 11–55 năm.
  4. Xây dựng hệ thống mẫu cấu trúc định hướng làm cơ sở khoa học tin cậy cho các biện pháp kỹ thuật lâm sinh và phục hồi rừng tự nhiên bền vững.
  5. Nâng tầm phương pháp luận nghiên cứu sinh thái học rừng nhiệt đới tại Việt Nam tiệm cận với các tiêu chuẩn quan trắc định vị dài hạn của mạng lưới lâm nghiệp quốc tế.