Tổng quan về luận án

Nghiên cứu động thái cấu trúc hệ sinh thái rừng tự nhiên là một trong những trụ cột then chốt của khoa học lâm sinh và sinh thái học hiện đại. Rừng tự nhiên lá rộng thường xanh nhiệt đới và á nhiệt đới tại miền Bắc Việt Nam, đặc biệt trên các dạng địa hình karst đá vôi và đồi núi thấp, là những hệ sinh thái có độ phức tạp cấu trúc cực kỳ cao, tính đa dạng sinh học phong phú và chứa đựng giá trị phòng hộ sinh thái đặc biệt quan trọng. Tuy nhiên, phần lớn các công trình nghiên cứu lâm học truyền thống tại Việt Nam trước đây chủ yếu tiếp cận theo phương pháp "lấy không gian thay thế thời gian" (space-for-time substitution / chronosequence) – tức là quan sát một thời điểm trên các lâm phần có độ tuổi hoặc giai đoạn phục hồi khác nhau để suy diễn tiến trình diễn thế. Phương pháp tĩnh này bộc lộ hạn chế lớn khi không thể lượng hóa chính xác các quá trình động học nội tại diễn ra liên tục theo trục thời gian thực, bao gồm quá trình sinh trưởng chuyển cấp đường kính ($O_k$), quá trình chết tự nhiên ($M_k$) và quá trình tái sinh bổ sung ($R$) vào tầng cây cao.

Nhằm giải quyết triệt để khoảng trống nghiên cứu (research gap) này, luận án tiến sĩ sinh học của nghiên cứu sinh Nguyễn Tiến Dũng với đề tài "Nghiên cứu động thái cấu trúc hệ sinh thái rừng lá rộng thường xanh tại Khu bảo tồn thiên nhiên Hang Kia - Pà Cò và Vườn quốc gia Xuân Sơn" (Chuyên ngành Sinh thái học, Mã số: 9 42 01 20; hướng dẫn khoa học: PGS.TS. Trần Văn Con và PGS.TS. Nguyễn Văn Sinh; bảo vệ năm 2020 tại Học viện Khoa học và Công nghệ – Viện Hàn lâm Khoa học và Công nghệ Việt Nam) đã được triển khai một cách công phu và chuẩn mực. Luận án xác lập một hệ thống ô tiêu chuẩn định vị (OTCĐV) dài hạn theo tiêu chuẩn quốc tế, theo dõi liên tục trong chu kỳ 11 năm (từ đợt đo chuẩn ban đầu năm 2007 đến đợt đo lặp kiểm chứng năm 2018) nhằm giải mã quy luật vận động cấu trúc và xây dựng mô hình toán học dự báo biến động lâm phần.

Câu hỏi nghiên cứu và giả thuyết khoa học của luận án được xác định rõ ràng:

  • Câu hỏi nghiên cứu 1 (RQ1): Các tham số cấu trúc cơ bản (tổ thành loài qua chỉ số $IVI%$, mật độ $N$, tổng tiết diện ngang $\sum G$, phân bố số cây theo cấp đường kính $N/D_{1.3}$ và chiều cao $N/H_{vn}$) tại các trạng thái rừng phục hồi (IIIA3, IIIB) và rừng nguyên sinh tương đối (IV) biến động như thế nào sau chu kỳ 11 năm?
  • Câu hỏi nghiên cứu 2 (RQ2): Mối quan hệ tương quan giữa số cây chuyển cấp ($O_k$), số cây chết ($M_k$) với cấp đường kính ($D_k$) và mật độ từng cấp kính ($N_k$) có tuân theo các quy luật toán học xác định hay diễn ra ngẫu nhiên?
  • Câu hỏi nghiên cứu 3 (RQ3): Có thể tích hợp các phương trình động thái thực nghiệm vào phần mềm mô phỏng chuyên dụng (MM&S) để dự đoán chính xác cấu trúc lâm phần trong tương lai dài hạn (55 năm) hay không?

Các giả thuyết nghiên cứu tương ứng:

  • Giả thuyết 1 (H1): Rừng lá rộng thường xanh duy trì tính ổn định cấu trúc thông qua cơ chế tự điều chỉnh nội tại, thể hiện qua phân bố $N/D_{1.3}$ dạng chữ J giảm (hàm giảm số mũ hoặc Meyer) không thay đổi bản chất sau 11 năm.
  • Giả thuyết 2 (H2): Quá trình chuyển cấp ($O_k$) và tử vong ($M_k$) phụ thuộc phi tuyến vào kích thước thân cây ($D_k$), được mô tả tối ưu bằng các phương trình hàm bậc 3 hoặc hàm mũ thay vì tuyến tính đơn giản.
  • Giả thuyết 3 (H3): Trạng thái rừng già thành thục (IV) đạt đến thế cân bằng động sinh thái học (dynamic demographic equilibrium) khi lượng tái sinh bổ sung ($R$) bù đắp tương đương lượng cây chết ($M$), trong khi các trạng thái phục hồi (IIIA3, IIIB) đang trong xu thế tích lũy sinh khối và chuyển dịch tổ thành loài.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa của lý thuyết chu trình khảm tái sinh rừng nhiệt đới (Mosaic-cycle theory), lý thuyết động thái khoảng trống tán (Gap dynamics) và lý thuyết cân bằng nhân khẩu học thực vật (Demographic equilibrium). Quy mô nghiên cứu bao gồm 6 OTC định vị diện tích $10.000\text{ m}^2$ (1 ha/ô, tổng diện tích 6 ha), được theo dõi trên 2 hệ sinh thái đại diện cho vùng núi cao đá vôi (Hang Kia – Pà Cò, Hòa Bình) và vùng đồi núi thấp chuyển tiếp (Xuân Sơn, Phú Thọ). Nghiên cứu mang ý nghĩa đột phá khi định lượng hóa các chỉ số động thái thực nghiệm, đặt nền tảng khoa học vững chắc cho việc thiết lập các mẫu định hướng cấu trúc và điều tiết biện pháp kỹ thuật lâm sinh bền vững.

Literature Review và Positioning

Nghiên cứu động thái rừng trên thế giới đã trải qua quá trình phát triển sâu rộng với ba dòng nghiên cứu (streams) chủ đạo: động thái tái sinh, động thái sinh trưởng - cấu trúc và động thái diễn thế quần xã.

Ở dòng nghiên cứu động thái tái sinh, tâm điểm của y văn thế giới xoay quanh cuộc tranh luận kinh điển giữa hai trường phái:

  • Trường phái biến động ngẫu nhiên theo chu trình khảm: Bắt đầu từ công trình nền tảng của Aubréville (1938) tại rừng mưa nhiệt đới Tây Phi. Tác giả quan sát thấy cây con của các loài ưu thế tầng trên cực kỳ hiếm gặp ở tầng dưới và đưa ra lý luận "bức khảm tái sinh" (Mosaic-cycle / Regeneration cycle), khẳng định tổ thành rừng biến đổi liên tục trong không gian và thời gian, không có một tổ hợp loài nào đạt thế cân bằng vĩnh viễn. Sau đó, Van Steenis (1956) phát triển sâu hơn thuyết "tái sinh vệt" và động thái khoảng trống tán rừng (gap dynamics), chỉ ra rằng cây tiên phong ưa sáng nhanh chóng lấp đầy các lỗ trống tán rừng do cây già gãy đổ tạo nên (loài liền vết sẹo), tạo điều kiện cho các loài chịu bóng xuất hiện ở giai đoạn kế tiếp.
  • Trường phái cân bằng động tại chỗ: Ngược lại với Aubréville, các nghiên cứu của David & Richards (1933), Beard (1964), Sun (1960) và Rollet (1969) tại rừng nhiệt đới Nam Mỹ khẳng định các loài có mật độ cao ở tầng trên cũng chiếm ưu thế ở các cấp kích thước nhỏ tầng dưới, minh chứng cho hiện tượng tái sinh tại chỗ liên tục và duy trì trạng thái ổn định lâu dài của quần xã thực vật cao đỉnh (climax).

Ở dòng nghiên cứu sinh trưởng và tử vong, các nhà khoa học quốc tế từ thập niên 1990 đã ứng dụng hệ thống ô định vị dài hạn để định lượng hóa các tham số nhân khẩu học:

  • Felfili (1997) theo dõi 6 năm trên 151 ô định vị ($10 \times 20\text{ m}$) tại rừng hành lang dọc sông Gama (Brazil), ghi nhận tỷ lệ chết trung bình hàng năm của quần xã là $3,5%/\text{năm}$, tỷ lệ bổ sung là $2,7%/\text{năm}$, tốc độ tăng trưởng đường kính trung bình của các cây $D_{1.3} \ge 10\text{ cm}$ đạt $0,25\text{ cm}/\text{năm}$, chứng minh khả năng tự duy trì cấu trúc của rừng nhiệt đới không bị xáo trộn.
  • Guillermo Castro Marín et al. (2005) nghiên cứu động thái cấu trúc trên 4 ha rừng tại Khu bảo tồn Chacocente (Nicaragua) giai đoạn 1993–2000, phát hiện tỷ lệ tử vong và bổ sung tương ứng là $4,5%$ và $2,5%/\text{năm}$ ở rừng rụng lá, trong khi ở rừng ven sông là $4,2%$ và $4,0%/\text{năm}$, đi cùng phân bố đường kính hình chữ J đảo ngược chiếm hơn $80%$ mật độ cá thể.
  • Takahashi et al. (2003) tại miền Bắc Nhật Bản chỉ ra tỷ lệ chết trung bình của tầng cây gỗ chỉ ở mức $0,57%/\text{năm}$ và sự cạnh tranh dưới tán (lớp tre lùn dày đặc) chi phối mạnh mẽ đến sự thành bại của cây tái sinh bổ sung hơn là sự cạnh tranh giữa các cây tầng cao.
  • Wei Xiang et al. (2016) tại vùng Đông Bắc Trung Quốc đã phát triển các mô hình tái sinh bổ sung dạng nhị thức âm hỗn hợp (negative binomial mixed models), khẳng định nhiệt độ, lượng mưa và mật độ tiết diện ngang lâm phần là các nhân tố điều khiển quá trình nhập tầng.

Tại Việt Nam, các nghiên cứu kinh điển của Trần Ngũ Phương (1970) và Thái Văn Trừng (1978, 1999) đã đặt nền móng vững chắc cho sinh thái học rừng. Thái Văn Trừng khẳng định: "Ánh sáng là nhân tố sinh thái khống chế và điều khiển quá trình tái sinh tự nhiên trong thảm thực vật rừng". Tuy nhiên, hầu hết các tác giả trong nước giai đoạn trước (Phạm Xuân Hoàn, Ngô Út, Nguyễn Văn Thêm) tập trung chủ yếu vào giải tích thân cây hoặc ứng dụng hàm phân loại Fisher trên các ô tiêu chuẩn tạm thời. Gần đây, các công trình của Bùi Chính Nghĩa (2011), Trần Văn Con (2009, 2012) và Phạm Đức Chiến (2006, 2012) đã tiếp cận ô định vị nhưng chưa khái quát hóa thành hệ phương trình mô phỏng toàn diện cấu trúc động học trên các trạng thái rừng đặc thù vùng núi phía Bắc.

Luận án của Nguyễn Tiến Dũng định vị chính xác vị trí học thuật của mình bằng việc bắc nhịp cầu giữa lý thuyết sinh thái học quần xã quốc tế và dữ liệu thực nghiệm lâm học nhiệt đới tại Việt Nam. Bằng cách so sánh đối chiếu đa chiều giữa vùng núi đá vôi Hang Kia - Pà Cò và đồi núi thấp Xuân Sơn, luận án cung cấp bằng chứng thực nghiệm trực tiếp để kiểm chứng tính đúng đắn của lý thuyết cân bằng động trong điều kiện sinh thái nhiệt đới gió mùa Đông Nam Á.

Đó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 quan trọng vào việc hoàn thiện lý thuyết động thái cấu trúc lâm phần nhiệt đới thông qua các điểm đột phá sau:

  1. Mở rộng lý thuyết chu trình khảm và cân bằng động (Aubréville 1938; Richards 1933): Luận án chứng minh rằng trong hệ sinh thái rừng lá rộng thường xanh nhiệt đới miền Bắc Việt Nam, hai cơ chế tưởng chừng đối lập này thực chất cùng tồn tại song song tùy thuộc vào trạng thái diễn thế và đặc điểm vi lập địa. Tại trạng thái IV (Hang Kia – Pà Cò), tính liên tục của cấu trúc $N/D_{1.3}$ và sự bù đắp giữa tỷ lệ tái sinh bổ sung ($R$) và tỷ lệ chết ($M$) củng cố thuyết cân bằng động của Richards. Ngược lại, tại trạng thái IIIA3 và IIIB, sự biến động mạnh mẽ của các loài tiên phong ưa sáng và sự phân bố cụm của cây tái sinh trong các lỗ trống tán rừng minh chứng rõ nét cho cơ chế động thái khoảng trống tán (Gap dynamics) của Van Steenis.
  2. Lượng hóa toán học quy luật chuyển cấp và tử vong phi tuyến: Luận án đã bác bỏ giả định tuyến tính giản đơn trong biến động nhân khẩu học rừng, chứng minh rằng quan hệ giữa số cây chuyển cấp ($O_k$) và cấp đường kính ($D_k$) tuân theo hàm đa thức bậc 3 ($O_k = aD_k^3 + bD_k^2 + cD_k + d$), có ý nghĩa thống kê cao qua kiểm định ANOVA ($p < 0,001$).
  3. Chuyển dịch hệ hình (Paradigm shift) trong phân loại trạng thái rừng: Nghiên cứu thúc đẩy sự chuyển dịch từ hệ thống định tính tĩnh (như thang phân loại Loeschau truyền thống) sang hệ thống định lượng động, mô tả trạng thái rừng thông qua cấu trúc phân bố $N/D_{1.3}$, phổ đa dạng Rényi ($H_\alpha$) và các tốc độ biến động hàng năm ($M_p, R_p$).

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

Khung phân tích của luận án là sự tích hợp hữu cơ giữa ba trụ cột lý thuyết:

  • Lý thuyết Động thái khoảng trống và Tái sinh tầng tán (Gap-phase dynamics): Xem xét sự xáo trộn quy mô nhỏ tạo ra các vi môi trường ánh sáng, quyết định sự phân tầng và nhập tầng cây cao.
  • Lý thuyết Quần thể học thực vật và Ma trận chuyển tiếp (Plant Demography Matrix Theory): Tiếp cận lâm phần như một hệ thống cấu trúc phân lớp cấp kính, trong đó dòng luân chuyển cá thể được xác định bởi: $$\Delta N_k(t) = R_k(t) + O_{k-1}(t) - O_k(t) - M_k(t)$$ Trong đó $\Delta N_k(t)$ là biến thiên số cây ở cấp kính $k$, $R_k$ là tái sinh bổ sung, $O_{k-1}$ là số cây từ cấp kính dưới chuyển lên, $O_k$ là số cây chuyển lên cấp kính trên, và $M_k$ là số cây chết tự nhiên.
  • Lý thuyết Đa dạng sinh học thông tin (Information-theoretic Biodiversity Spectrum - Rényi 1961): Ứng dụng phổ hàm Rényi để phân tích tính đa dạng loài liên tục, khắc phục triệt để nhược điểm của các chỉ số đơn lẻ rời rạc.

Khung phân tích này kết hợp chặt chẽ giữa thu thập số liệu trắc lượng lâm học chuẩn xác ngoài hiện trường và mô phỏng trên nền tảng phần mềm MM&S (Model Maker & Simulation), thiết lập đường dẫn dự báo 55 năm cho lâm phần.

Điều kiện biên (Boundary conditions): Khung phân tích áp dụng đối với tầng cây gỗ có $D_{1.3} \ge 10\text{ cm}$ thuộc kiểu rừng tự nhiên lá rộng thường xanh á nhiệt đới và nhiệt đới gió mùa tại khu vực đồi núi miền Bắc Việt Nam, trên nền địa chất đá vôi (karst) và đất đồi Ferasol/Alisol, độ cao tuyệt đối từ $300\text{ m}$ đến $1.300\text{ m}$, lượng mưa hàng năm $1.100 - 2.200\text{ mm}$, trong điều kiện không chịu tác động khai thác trực tiếp từ con người.

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

Thiết kế nghiên cứu

Luận án tuân thủ lập trường nhận thức luận thực chứng định lượng kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), xem xét các quá trình sinh thái diễn ra khách quan và có thể lượng hóa thông qua các đo đạc thực nghiệm trường kỳ.

Thiết kế nghiên cứu đa tầng (multi-level longitudinal design) được thực hiện trên 6 ô tiêu chuẩn định vị (OTCĐV) cố định, diện tích mỗi ô là $10.000\text{ m}^2$ ($100 \times 100\text{ m}$, quy mô 1 ha), chia thành 25 ô phụ ($20 \times 20\text{ m}$):

  • Khu bảo tồn thiên nhiên Hang Kia – Pà Cò (Hòa Bình):
    • Ô HB01: Tọa độ $20^\circ43'40''\text{N} - 104^\circ56'27''\text{E}$; Độ cao $1.060\text{ m}$; Độ dốc $55^\circ$ hướng Đông Nam; Núi đá thấp; Đất Ferasol/Alisol; Trạng thái IIIA3.
    • Ô HB03: Tọa độ $20^\circ44'57''\text{N} - 104^\circ54'37''\text{E}$; Độ cao $1.230\text{ m}$; Độ dốc $35^\circ$ hướng Đông Nam; Núi đá thấp; Trạng thái IIIA3.
    • Ô HB06: Tọa độ $20^\circ44'15,28''\text{N} - 104^\circ56'42,97''\text{E}$; Độ cao $980\text{ m}$; Độ dốc $35^\circ$ hướng Đông Nam; Núi đá thấp; Trạng thái IV.
  • Vườn quốc gia Xuân Sơn (Phú Thọ):
    • Ô XS01: Tọa độ $21^\circ06'36,2''\text{N} - 104^\circ57'47,6''\text{E}$; Độ cao $304\text{ m}$; Độ dốc $25^\circ$ hướng Đông Nam; Đồi núi thấp; Trạng thái IIIB.
    • Ô XS02: Tọa độ $21^\circ07'55,9''\text{N} - 104^\circ57'4,5''\text{E}$; Độ cao $540\text{ m}$; Độ dốc $30^\circ$ hướng Đông Nam; Trạng thái IIIB.
    • Ô XS03: Tọa độ $21^\circ08'9,8''\text{N} - 104^\circ55'39,5''\text{E}$; Độ cao $455\text{ m}$; Độ dốc $20^\circ$ hướng Đông Nam; Trạng thái IIIB.

Quy trình nghiên cứu rigorous

Quy trình thu thập và chuẩn hóa dữ liệu thực địa được tiến hành nghiêm ngặt:

  • Định danh thực vật 2 giai đoạn: Giai đoạn 1 nhận dạng sơ bộ tại hiện trường; giai đoạn 2 thu mẫu tiêu bản (lá, hoa, quả) đối với các cá thể chưa rõ ràng ($sp1, sp2...$) và tiến hành giám định hình thái so sánh cùng các chuyên gia thực vật học đầu ngành thuộc Viện Sinh thái và Tài nguyên Sinh vật, đối chiếu với các tài liệu chuẩn quốc gia: Cây cỏ Việt Nam (Phạm Hoàng Hộ), Thực vật chí Việt Nam, Thực vật rừng (Lê Mộng Chân, Trần Hợp).
  • Giao thức đo đếm trắc lượng: Tất cả các cá thể thân gỗ có $D_{1.3} \ge 10\text{ cm}$ trong 6 OTCĐV được đánh số thẻ nhôm cố định, vẽ sơ đồ tọa độ vị trí không gian trong ô $20 \times 20\text{ m}$. Đường kính ngang ngực ($D_{1.3}$) được đo chính xác đến từng milimet bằng thước đo vanh tại vị trí vạch sơn $1,3\text{ m}$. Chiều cao vút ngọn ($H_{vn}$) được đo chính xác đến decimet ($0,1\text{ m}$) bằng máy đo cao siêu âm quang học Vertex IV.
  • Kiểm kê động thái sau 11 năm (2007–2018): Tái điều tra $100%$ các cá thể cũ; lập danh mục chính xác cây chết tự nhiên ($M$), phân loại nguyên nhân; thống kê cây tái sinh bổ sung nhập tầng ($R$, những cây đạt $D_{1.3} \ge 10\text{ cm}$ tại thời điểm 2018); xác định số cây chuyển cấp đường kính ($O_k$).

Data và phân tích

Các công cụ phân tích toán học và thống kê chuyên sâu bao gồm:

  1. Chỉ số giá trị quan trọng (Importance Value Index - $IVI%$): $$IVI% = \frac{F% + N% + G%}{3}$$ Trong đó $F%$ là tần số xuất hiện tương đối, $N%$ là mật độ tương đối, $G%$ là độ tàn che/tiết diện ngang tương đối. Hệ số tổ thành chỉ tính cho các loài có $IVI% > 5%$.
  2. Tỷ lệ hỗn loài: $HL = S/N$ (biểu thị dưới dạng phân số $1/n$, với $n = N/S$).
  3. Chỉ số đa dạng tổng quát Rényi ($H_\alpha$): $$H_\alpha = \frac{1}{1-\alpha}\ln\left(\sum_{i=1}^S p_i^\alpha\right)$$ Tham số $\alpha$ được quét liên tục qua các giá trị: $\alpha = 0$ ($H_0 = \ln S$, độ phong phú loài); $\alpha = 1$ ($H_1 = H'$, Shannon-Wiener); $\alpha = 2$ ($H_2 = \ln(1/D)$, Simpson); $\alpha \to \infty$ ($H_\infty = \ln(1/p_{max})$).
  4. Hệ thống mô hình hóa động thái (MM&S Software): Xây dựng mô hình hệ thống động (System Dynamics) dạng văn bản và sơ đồ dòng chảy (flow diagram). Phương trình chuyển cấp ($O_k-D_k$) và tử vong ($M_k-D_k$) được ước lượng tham số bằng phương pháp phi tuyến và hồi quy bình phương nhỏ nhất (OLS). Kiểm định độ tin cậy mô hình thông qua hệ số xác định ($R^2$), giá trị $F$-test phân tích phương sai (ANOVA) và sai số chuẩn ước lượng ($SE$).

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

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

  1. Sự ổn định cấu trúc và phân bố mật độ cá thể theo cấp kính ($N/D_{1.3}$): Tại tất cả 6 OTCĐV thuộc cả 3 trạng thái rừng (IIIA3, IIIB, IV), đường cong phân bố $N/D_{1.3}$ tại cả hai thời điểm đo (2007 và 2018) đều tuân theo quy luật giảm số mũ hoặc hàm khoảng cách Meyer chặt chẽ (đồ thị dạng chữ J ngược). Mật độ tập trung áp đảo ở các cấp đường kính nhỏ ($10 - 20\text{ cm}$, chiếm từ $55% - 75%$ tổng số cây) và giảm nhanh theo cấp số nhân khi đường kính tăng lên ($D > 40\text{ cm}$). Sau 11 năm, dạng phân bố tổng thể không bị phá vỡ, chứng minh tính bền vững sinh thái tự nhiên của các kiểu rừng này.

  2. Quy luật động thái chuyển cấp ($O_k$) tuân theo hàm bậc 3: Quá trình phân tích tương quan giữa số cây chuyển ra khỏi cấp kính ($O_k$) và cấp đường kính ($D_k$) khẳng định mối quan hệ phi tuyến dạng hàm bậc 3: $$O_k = aD_k^3 + bD_k^2 + cD_k + d$$ Kết quả phân tích phương sai tại OTC Hang Kia – Pà Cò (trạng thái IIIA3) cho thấy tương quan này có ý nghĩa thống kê tuyệt đối với $p < 0,001$ ($F_{calc} \gg F_{table}$). Tốc độ chuyển cấp diễn ra mạnh nhất ở các cấp kính nhỏ ($10 - 25\text{ cm}$) do áp lực cạnh tranh không gian dinh dưỡng và giảm dần ở các cấp kính thành thục lớn.

  3. Hiện tượng tử vong ($M_k$) tập trung phi ngẫu nhiên: Quá trình chết tự nhiên không phân bố rải rác ngẫu nhiên mà tập trung chủ yếu ở hai nhóm đối tượng: (i) Những cây tầng dưới bị chèn ép ánh sáng nghiêm trọng ở cấp kính nhỏ ($10 - 15\text{ cm}$); và (ii) Một số cây già cỗi vượt tán ở cấp kính lớn ($D > 50\text{ cm}$) bị đổ gãy do gió bão trên địa hình núi đá vôi hiểm trở. Tỷ lệ chết trung bình hàng năm ($M_p$) tại các OTC dao động trong khoảng $0,8% - 1,6%/\text{năm}$, tương đồng với các khu rừng nhiệt đới ổn định trên thế giới.

  4. Cơ chế bù trừ nhân khẩu học ở trạng thái rừng già Climax (IV): Tại ô HB06 (Trạng thái IV, Hang Kia – Pà Cò), sau 11 năm, số lượng cây tái sinh bổ sung nhập tầng ($R$) đạt mức cân bằng xấp xỉ với số lượng cây chết tự nhiên ($M$) ($R \approx M$), tổng tiết diện ngang lâm phần ($\sum G$) duy trì ổn định quanh ngưỡng $32 - 36\text{ m}^2/\text{ha}$. Ngược lại, tại các ô phục hồi IIIA3 và IIIB (Xuân Sơn), tốc độ bổ sung vượt trội tốc độ chết ($R > M$), dẫn tới mật độ thân cây và tổng tiết diện ngang gia tăng ròng, phản ánh tiến trình diễn thế tích lũy sinh khối.

Implications đa chiều

  • Về mặt lý thuyết: Kết quả nghiên cứu cung cấp cơ sở dữ liệu thực nghiệm chuẩn xác khẳng định tính thích ứng của lý thuyết cân bằng động trong điều kiện rừng núi đá nhiệt đới. Sự biến thiên của phổ Rényi chứng minh rằng tính đa dạng loài không suy giảm mà có xu hướng dịch chuyển theo hướng tăng độ đồng đều cá thể giữa các loài ưu thế bản địa.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình thiết lập, quản lý và phân tích dữ liệu OTCĐV 1 ha kết hợp phần mềm mô phỏng MM&S. Phương pháp này có khả năng nhân rộng trên toàn bộ hệ thống rừng đặc dụng và rừng phòng hộ tại Việt Nam.
  • Về mặt quản lý thực tiễn: Luận án đề xuất hệ thống các "mẫu định hướng cấu trúc" ($N/D_{1.3}$ target models) cho từng trạng thái rừng. Cụ thể, xác định mật độ chuẩn, tỷ lệ phân bố tiết diện ngang mục tiêu ($\sum G_{D>40\text{cm}}$) để làm căn cứ điều tiết lâm sinh:
    • Đối với trạng thái IIIA3: Cần áp dụng biện pháp tỉa thưa vệ sinh, giải phóng không gian dinh dưỡng cho các loài cây mục đích tái sinh có triển vọng, hỗ trợ vượt qua giai đoạn chèn ép tầng dưới.
    • Đối với trạng thái IIIB: Áp dụng khoanh nuôi xúc tiến tái sinh tự nhiên kết hợp làm giàu rừng bằng các loài cây bản địa có giá trị sinh thái và kinh tế cao.
  • Về mặt chính sách: Cung cấp cơ sở khoa học định lượng để Ban quản lý KBTTN Hang Kia – Pà Cò và VQG Xuân Sơn quy hoạch phân khu chức năng (phân khu bảo vệ nghiêm ngặt và phân khu phục hồi sinh thái), phục vụ chính sách chi trả dịch vụ môi trường rừng và định giá hấp thụ carbon (REDD+).

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án cũng ghi nhận một số hạn chế khách quan:

  1. Khoảng thời gian điều tra: Dữ liệu thực nghiệm được quan trắc qua 2 chu kỳ đo đếm (2007 và 2018, chu kỳ 11 năm). Mặc dù 11 năm là đủ dài để ghi nhận các biến đổi lâm học, nhưng để mô hình hóa toàn diện chu trình sống của cây gỗ kéo dài hàng trăm năm, cần tiếp tục duy trì các đợt đo đếm định kỳ lần 3 (2023), lần 4 (2028).
  2. Ngưỡng kích thước điều tra: Luận án tập trung vào tầng cây cao ($D_{1.3} \ge 10\text{ cm}$). Động thái của tầng cây con tái sinh ($H < 1,3\text{ m}$) và cây nhỡ ($H > 1,3\text{ m}, D_{1.3} < 10\text{ cm}$) mới chỉ được khảo sát thông qua các ô phụ dạng bản mà chưa được tích hợp hoàn toàn vào hệ phương trình vi phân mô phỏng tổng thể.
  3. Biến số vi khí hậu và thổ nhưỡng: Nghiên cứu chưa đo đạc liên tục các thông số khí tượng tiểu khí hậu (nhiệt độ, ẩm độ không khí dưới tán, bức xạ quang hợp PAR) và biến động dinh dưỡng đất vi mô theo từng tháng để đưa vào mô hình như các biến ngoại sinh.

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

  • Hướng 1: Mở rộng mô hình hóa cá thể không gian (Spatially-explicit individual-based modeling) ứng dụng hàm $K(t)$ của Ripley để phân tích tương tác cạnh tranh không gian 3 chiều giữa các cây lân cận.
  • Hướng 2: Nghiên cứu liên kết động thái cấu trúc với các đặc tính chức năng lá (functional leaf traits), mật độ gỗ và khả năng thích ứng hạn hán của các loài cây gỗ trên karst.
  • Hướng 3: Tích hợp mô hình biến đổi khí hậu toàn cầu (IPCC scenarios) vào phần mềm MM&S nhằm dự báo tính dễ bị tổn thương và khả năng chống chịu của hệ sinh thái rừng lá rộng thường xanh trước các hiện tượng thời tiết cực đoan.
  • Hướng 4: Lượng hóa dòng sinh khối và cân bằng carbon ngầm dưới đất (soil organic carbon) đi kèm với động thái tầng cây cao.

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

Luận án của nghiên cứu sinh Nguyễn Tiến Dũng tạo ra những tác động lan tỏa sâu rộng trên nhiều bình diện:

  • Ảnh hưởng học thuật (Academic Impact): Tạo lập bộ dữ liệu thực nghiệm 11 năm vô giá trên hệ thống 6 ha OTCĐV quốc gia. Công trình là nguồn trích dẫn nền tảng cho các nghiên cứu tiếp theo về sinh thái học rừng nhiệt đới, mô hình hóa toán sinh học và động thái thảm thực vật karst tại Đông Nam Á.
  • Chuyển đổi ngành lâm nghiệp (Forestry Transformation): Cung cấp công cụ phần mềm và mô hình toán học dự báo cấu trúc, giúp chuyển đổi căn bản từ công tác điều tra quy hoạch rừng truyền thống sang quản lý số hóa tài nguyên rừng dựa trên mô phỏng kịch bản.
  • Tác động chính sách bảo tồn (Policy Influence): Trực tiếp hỗ trợ Cục Lâm nghiệp, Sở Nông nghiệp và PTNT các tỉnh Hòa Bình, Phú Thọ xây dựng phương án quản lý rừng bền vững giai đoạn 2021–2030, tối ưu hóa nguồn lực đầu tư bảo tồn thiên nhiên.
  • Lợi ích xã hội và môi trường: Bảo tồn cấu trúc phức tán của rừng thường xanh giúp nâng cao năng lực giữ nước đầu nguồn, chống sạt lở xói mòn đất vùng núi Tây Bắc, bảo vệ nguồn gen quý hiếm và góp phần thực hiện cam kết quốc gia về giảm phát thải khí nhà kính (Net Zero).

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

  • Nghiên cứu sinh và học viên cao học chuyên ngành Sinh thái học, Lâm học: Tiếp cận phương pháp luận mẫu mực về thiết lập OTCĐV dài hạn, thuật toán xử lý dữ liệu nhân khẩu học thực vật và ứng dụng phần mềm MM&S trong mô phỏng sinh thái.
  • Các nhà khoa học lâm nghiệp và sinh thái học: Kế thừa cơ sở dữ liệu thực nghiệm chuẩn xác về các tham số động thái $O_k, M_k, R_p, M_p$ để so sánh xuyên lục địa với các hệ sinh thái rừng nhiệt đới tại Amazon, Trung Phi hoặc Nam Á.
  • Ban quản lý các Vườn quốc gia và Khu bảo tồn thiên nhiên: Ứng dụng trực tiếp các mẫu định hướng cấu trúc ($N/D_{1.3}$) để xây dựng quy trình kỹ thuật khoanh nuôi, làm giàu và giám sát đa dạng sinh học định kỳ.
  • Các nhà hoạch định chính sách lâm nghiệp: Sử dụng bằng chứng định lượng của luận án để ban hành các tiêu chuẩn kỹ thuật quốc gia về điều tra động thái rừng và định mức kinh tế kỹ thuật trong bảo tồn hệ sinh thái rừng đặc dụ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 của luận án là việc chứng minh và lượng hóa cơ chế cân bằng động nhân khẩu học phi tuyến trong hệ sinh thái rừng lá rộng thường xanh trên núi đá vôi miền Bắc Việt Nam. Luận án đã mở rộng Lý thuyết Chu trình khảm tái sinh của Aubréville (1938) và Lý thuyết Cân bằng sinh thái của Richards (1933) bằng cách tích hợp mô hình ma trận động thái cấp kính. Nghiên cứu chỉ ra rằng ở trạng thái rừng già (trạng thái IV), thế cân bằng sinh thái không phải là một trạng thái tĩnh bất biến mà là một cân bằng động liên tục, trong đó tốc độ tử vong ($M_k$) và tốc độ bổ sung ($R$) bù trừ cho nhau một cách hoàn hảo theo chu kỳ thập kỷ, giữ vững cấu trúc tiết diện ngang lâm phần $\sum G \approx 32 - 36\text{ m}^2/\text{ha}$.

2. Đột phá về 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ế là gì?

So với các nghiên cứu quốc tế như công trình của Felfili (1997) tại Brazil (vốn chỉ sử dụng phân tích thống kê mô tả sai phân mật độ trên các ô $10 \times 20\text{ m}$) và Castro Marín et al. (2005) tại Nicaragua (áp dụng mô hình logarit tuyến tính hóa đơn giản), luận án của Nguyễn Tiến Dũng đạt được bước tiến đột phá về phương pháp luận:

  1. Thiết lập mạng lưới OTCĐV diện tích lớn ($10.000\text{ m}^2$/ô) theo chuẩn quốc tế, kiểm soát toàn diện sai số biên viền (edge effects).
  2. Phát triển và kiểm định hệ phương trình phi tuyến bậc 3 ($O_k = aD_k^3 + bD_k^2 + cD_k + d$) để mô tả dòng chuyển cấp hạt nhân.
  3. Tích hợp trực tiếp các phương trình thực nghiệm vào môi trường mô phỏng hệ thống động MM&S, cho phép chạy kịch bản dự báo cấu trúc 5 bước thời gian (tương đương 55 năm) – điều mà các nghiên cứu trước đây chưa thực hiện hoàn chỉnh.

3. Phát hiện bất ngờ nhất (counter-intuitive) có sự hỗ trợ của dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là mối tương quan giữa cấp đường kính ($D_k$) và tỷ lệ tử vong tự nhiên ($M_k$) không tuân theo hàm giảm đơn điệu giản đơn mà thể hiện tính phân hóa lưỡng cực:

  • Dữ liệu thực nghiệm tại các OTC cho thấy tỷ lệ chết không chỉ tập trung ở tầng cây con bị chèn ép ánh sáng ($D_{1.3} = 10 - 15\text{ cm}$) mà còn xuất hiện một đỉnh tử vong cục bộ ở các cấp đường kính lớn ($D_{1.3} > 50\text{ cm}$).
  • Về mặt lý thuyết sinh thái, điều này được giải thích bởi đặc thù lập địa núi đá vôi: tầng đất mỏng, rễ cây bám nông vào các kẽ đá kết hợp với địa hình dốc cao ($35^\circ - 55^\circ$). Do đó, các cây cổ thụ vượt tán khi chịu tác động cơ học của gió lốc và bão nhiệt đới dễ bị bật gốc gãy đổ hơn nhiều so với các lâm phần trên nền đất thịt sâu ở đồng bằng.

4. Luận án có cung cấp quy trình nhân rộng và lặp lại (replication protocol) rõ ràng không?

Có, luận án cung cấp một giao thức nhân rộng cực kỳ chi tiết và minh bạch:

  • Bản đồ tọa độ UTM chính xác của 6 OTCĐV (bảng 2.2).
  • Quy trình chia tách 25 ô phụ ($20 \times 20\text{ m}$) có gắn mốc cố định và đánh số thẻ cây bằng nhôm chống oxy hóa.
  • Biểu mẫu điều tra chuẩn hóa (bảng 2.1) ghi chép chi tiết: số hiệu cây, tên loài, chu vi $G_{1.3}$, chiều cao $H_{vn}$, tọa độ cục bộ và tình trạng sống/chết.
  • Quy trình chuẩn hóa mã lệnh và phương trình nạp tham số trong phần mềm MM&S, cho phép bất kỳ nhóm nghiên cứu độc lập nào cũng có thể kế thừa và thực hiện đo đếm kiểm chứng lặp lại trong tương lai.

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 tiếp theo tập trung vào 4 trọng tâm:

  1. Thực hiện kỳ điều tra đo đếm lần 3 (năm thứ 15) và lần 4 (năm thứ 20) trên 6 OTCĐV tại Hang Kia - Pà Cò và Xuân Sơn để hiệu chỉnh tham số mô hình MM&S.
  2. Mở rộng gắn máy đo tăng trưởng tự động (Dendrometer bands) trên các loài cây gỗ ưu thế để quan trắc biến động đường kính theo mùa và chu kỳ tuần trăng.
  3. Ứng dụng công nghệ quét LiDAR mặt đất (Terrestrial Laser Scanning - TLS) và viễn thám không người lái (UAV-LiDAR) để số hóa 3D cấu trúc kiến trúc tầng tán rừng.
  4. Tích hợp dữ liệu chuỗi thời gian phân giải cao vào mạng lưới giám sát động thái rừng nhiệt đới toàn cầu (CTFS-ForestGEO).

Kết luận

Luận án tiến sĩ của Nguyễn Tiến Dũng là một công trình khoa học mẫu mực, toàn diện và có giá trị học thuật xuất sắc. Những đóng góp then chốt của luận án được đúc kết cụ thể như sau:

  1. Xác lập hệ thống cơ sở dữ liệu động thái thực nghiệm 11 năm (2007–2018): Cung cấp bức tranh toàn cảnh định lượng về biến động cấu trúc tổ thành ($IVI%$), mật độ ($N$), tiết diện ngang ($\sum G$) và tính đa dạng loài theo phổ Rényi ($H_\alpha$) trên 6 ha OTCĐV tại KBTTN Hang Kia – Pà Cò và VQG Xuân Sơn.
  2. Tìm ra quy luật toán học mô phỏng chuyển cấp đường kính ($O_k$): Khẳng định mối quan hệ giữa số cây chuyển cấp và đường kính tuân theo hàm đa thức bậc 3 có ý nghĩa thống kê cao ($p < 0,001$).
  3. Lượng hóa các tham số nhân khẩu học cốt lõi: Xác định chính xác tỷ lệ tử vong tự nhiên ($M_p$), tỷ lệ chuyển cấp ($R_p$), hệ số chết ($M_r$) và tỷ lệ tái sinh bổ sung ($R$) qua các trạng thái rừng IIIA3, IIIB và IV.
  4. Mô phỏng thành công động thái lâm phần bằng phần mềm MM&S: Dự báo cấu trúc phân bố $N/D_{1.3}$ và $G/D_{1.3}$ trong tương lai 55 năm, mở ra hướng ứng dụng công nghệ mô phỏng hiện đại trong lâm nghiệp Việt Nam.
  5. Xây dựng hệ thống mẫu định hướng cấu trúc ($N/D_{1.3}$ target models): Cung cấp căn cứ khoa học trực tiếp phục vụ công tác quy hoạch, điều tiết lâm sinh và quản lý rừng bền vững tại khu vực miền Bắc.
  6. Thúc đẩy sự tiến bộ của hệ hình nghiên cứu sinh thái học: Kết hợp nhuần nhuyễn giữa phương pháp ô tiêu chuẩn định vị dài hạn với lý thuyết cân bằng động quốc tế, mở ra 3 hướng nghiên cứu mới về sinh thái học chức năng, mô hình hóa không gian 3D và thích ứng biến đổi khí hậu cho thảm thực vật nhiệt đới gió mùa.