Tổng quan luận án

Đề tài luận án tiến sĩ lâm nghiệp: "Nghiên cứu cấu trúc và xây dựng mô hình tăng trưởng đường kính rừng tự nhiên lá rộng thường xanh một số khu rừng đặc dụng miền Bắc Việt Nam" do nghiên cứu sinh Nguyễn Thị Thu Hiền thực hiện tại Trường Đại học Nông lâm - Đại học Thái Nguyên, dưới sự hướng dẫn khoa học của PGS.TS. Trần Thị Thu Hà và PGS.TS. Trần Văn Con; bảo vệ năm 2015 thuộc chuyên ngành Lâm sinh (Mã số: 62.05).

  • Tính cấp thiết và khoảng trống nghiên cứu: Rừng tự nhiên nhiệt đới là hệ sinh thái luôn vận động thông qua các quá trình sinh trưởng, tái sinh và diễn thế phức tạp. Tại Việt Nam, độ che phủ rừng từng giảm từ 43% năm 1943 xuống còn 27% năm 1986, và tăng trở lại đạt 41% vào đầu năm 2014 nhờ các chương trình phục hồi rừng. Tuy nhiên, chất lượng rừng tự nhiên vẫn tiếp tục suy giảm cả về trữ lượng gỗ lẫn tính đa dạng sinh học. Rừng tự nhiên chiếm tới 74,51% tổng diện tích rừng cả nước. Để quản lý rừng tự nhiên bền vững, lâm nghiệp cần giải quyết ba vấn đề cốt lõi: (i) xác định tăng trưởng đường kính làm cơ sở định lượng đường kính khai thác tối thiểu và luân kỳ khai thác; (ii) lượng hóa tăng trưởng trữ lượng để xác định sản lượng khai thác cho phép hàng năm; và (iii) mô hình hóa động thái cấu trúc lâm phần nhằm dự báo trạng thái rừng tương lai. Hầu hết các nghiên cứu trước đây tại Việt Nam chủ yếu dựa vào số liệu thu thập một lần từ các ô tiêu chuẩn tạm thời và phương pháp giải tích thân cây hoặc đẽo vát cây tiêu chuẩn, dẫn đến độ chính xác bị hạn chế và số liệu động thái còn tản mạn.

  • Mục tiêu nghiên cứu: Góp phần hoàn thiện phương pháp luận về phân nhóm loài, tái sinh bổ sung, chuyển cấp và xác định mô hình mô phỏng các quá trình động thái bao gồm:

    1. Quá trình tăng trưởng đường kính;
    2. Quá trình chết;
    3. Quá trình tái sinh bổ sung của rừng tự nhiên lá rộng thường xanh ở một số khu rừng đặc dụng miền Bắc Việt Nam.

  • Đối tượng và phạm vi nghiên cứu:

    • Đối tượng nghiên cứu: Tầng cây cao có đường kính ngang ngực $D_{1.3} \ge 10\text{ cm}$ trong các lâm phần rừng tự nhiên lá rộng thường xanh.
    • Phạm vi không gian: 04 khu rừng đặc dụng miền Bắc Việt Nam với tổng số 21 ô tiêu chuẩn định vị (OTCĐV, mỗi ô diện tích 1 ha), bao gồm: Vườn Quốc gia (VQG) Ba Bể, tỉnh Bắc Kạn (06 OTCĐV); VQG Vũ Quang, tỉnh Hà Tĩnh (06 OTCĐV); VQG Xuân Sơn, tỉnh Phú Thọ (03 OTCĐV); và Khu Bảo tồn thiên nhiên (KBT) Hang Kia - Pà Cò, tỉnh Hòa Bình (06 OTCĐV).
    • Phạm vi thời gian: Chu kỳ theo dõi và thu thập số liệu lặp 5 năm (giai đoạn 2007 - 2012).

Tổng quan tài liệu và vị trí của luận án

  • Các hướng nghiên cứu trên thế giới:
    • Về cấu trúc tổ thành và đa dạng loài: Richards (1974) chỉ ra rằng trong rừng nhiệt đới, rất ít họ chiếm quá 10% tổng số loài và 10 họ lớn nhất chỉ chiếm 40 - 50% số loài. Laura Klappenbach (2001) khẳng định thành phần loài biến đổi liên tục qua chuỗi diễn thế sinh thái.
    • Về quy luật phân bố số cây theo đường kính ($N/D$): Meyer (1952) ứng dụng hàm giảm liên tục mô tả phân bố $N/D$. Prodan và Patatscase (1964) sử dụng phân bố Beta; Nelson (1964) sử dụng hàm Logarit chuẩn $\text{Ln}(a) = a_0 + b_0\text{Ln}(d)$. Các tác giả Suzuki (1971), Preussner (1974), Bock và Diener (1972) tiếp cận phân bố đường kính như quá trình ngẫu nhiên và áp dụng lý thuyết chuỗi Markov.
    • Về quy luật tương quan chiều cao - đường kính ($H/D$): Ivanov (1930), Krauter (1958), Naslund (1929), Assmann (1936), Hohenadl (1936), Michailov (1934, 1952), Prodan (1944, 1965), Krenn (1946) và Meyer (1952) đã thiết lập các dạng phương trình hồi quy như $h = a + b_1 d + b_2 d^2$.
    • Về phân nhóm loài: Swaine và Whitmore (1988) phân chia loài theo đặc tính ưa sáng - chịu bóng hoặc tiên phong - cực đỉnh; Alder (1995) phân 265 loài tại Tapajos thành 16 nhóm trong mô hình CAFOGRM; Ong và Kleine (1995) phân 436 loài tại Sabah thành 22 nhóm cho mô hình DIPSIM; Kohler và Huth (1998) dùng mô hình FORMIND; Picard và Franc (2003) phân nhóm tối ưu hóa sai số tích lũy; Phillips et al. (2002) phân 575 loài tại Kalimantan thành 10 nhóm.
    • Về mô hình sinh trưởng, tỷ lệ chết và tái sinh: Vanclay (1994) xây dựng hệ thống phương trình cân bằng động thái qua các đại lượng chuyển vào (ingrowth/recruitment), chuyển ra (outgrowth), chết (mortality) và khai thác (harvesting). Stage và Renner (1988) nhấn mạnh 80% sai lệch dự báo thể tích bắt nguồn từ sai số ước lượng tỷ lệ chết. Các lý thuyết tự tỉa thưa và mật độ giới hạn được phát triển bởi Reineke (1933) qua chỉ số mật độ lâm phần $\text{Ln}(N) + 1{,}6\text{Ln}(d_g) = \text{const}$, Yoda et al. (1963) qua phương trình $\text{Ln}(w) + 1{,}5\text{Ln}(N) = K$, Wilson (1951) và Mitchell (1975) qua $\text{Ln}(N) + 2\text{Ln}(h_T) = \text{const}$. Tỷ lệ chết trong rừng ẩm nhiệt đới được tổng hợp dao động từ 1% đến 5%/năm:
Tác giả Kiểu rừng Cây nhỏ ($D < 10\text{ cm}$) Cây lớn ($D \ge 10\text{ cm}$) Đặc điểm ghi nhận
Alder (1993) Rừng ẩm nhiệt đới nửa rụng lá (Bobiri, Ghana) - Loài chịu bóng: 0,3 - 0,7%
- Loài ưa sáng: 1,2 - 1,5%
- Loài tiên phong: 3 - 4,5%		 Toàn bộ các cỡ Đo đếm 15 năm; rừng bị khai thác mạnh sau 10 - 15 năm
Bryan (1980) Rừng Dầu (Sarawak) 4 - 5% 2 - 3% OTCĐV thiết lập sau khai thác dưới 10 năm; định kỳ 5 năm
Nicholson (1979) Rừng Dầu (Sabah) - 1,06% Rừng không khai thác
Mervart (1972) Rừng ẩm hỗn giao (Nigeria) 5,5% 3,5% 6.224 cây quan sát trong 5 năm
Singh (1981) Rừng ẩm thường xanh (Ấn Độ) - 1,11% Khác nhau giữa các loài, không tương quan tốc độ sinh trưởng
Swaine (1992) Rừng ẩm nhiệt đới nửa rụng lá (Kade, Ghana) 3,4 - 6,5% 1,5% Rừng không khai thác đo 19 năm; tỷ lệ chết cao ở cây bị chèn ép
Tang (1976) Rừng Dầu (Trengganu, W. Malaysia) 2,8 - 10,2% - Tỷ lệ chết tương quan với % tiết diện ngang khai thác ($R^2 = 0{,}62$)
Silva (1989) Rừng lục địa (Tapajos, Brazilian Amazon) 3 - 4% 1,5 - 2,5% Rừng khai thác có tỷ lệ chết cao hơn 1 - 2% sau khai thác

Các mô hình tái sinh và dự báo động thái cấu trúc khác được phát triển bởi Botkin et al. (1972) với JABOWA, Shifley et al. (1982, 1993), De-Miguel et al. (2013) sử dụng mô hình hiệu ứng cố định (Fixed-effects) và mô hình hỗn hợp, Liang et al. (2005) dùng mô hình ma trận cho loài Pseudotsuga menziesii và Tsuga heterophylla, Ek (1974) và Pukkala et al. (2009) mô hình hóa lâm phần không đồng tuổi tại Phần Lan.

  • Các nghiên cứu tại Việt Nam:

    • Cấu trúc tổ thành và đa dạng: Thái Văn Trừng (1978, 1999) phân định ưu hợp và phức hợp; Nguyễn Hồng Quân (1982) nghiên cứu rừng Kon Hà Nừng; Nguyễn Văn Trương (1983); Nguyễn Ngọc Lung (1991); Phạm Hoàng Hộ (1999); Nguyễn Hải Tuất (1991) dùng bảng chéo $2\times 2$ và phương pháp 6 cây; Hà Thị Mừng (2000); Bùi Đoàn (2001); Nguyễn Văn Thêm (2004).
    • Phân bố $N/D$ và tương quan $H/D$: Đồng Sỹ Hiền (1974) dùng hàm Meyer, Pearson và hàm $\text{Log}(h) = a + b\text{Log}(d_{1.3})$; Nguyễn Hải Tuất (1982) dùng phân bố giảm và phân bố khoảng cách; Lê Minh Trung (1991) so sánh phân bố Poisson, Weibull, Hyperbol, Meyer; Đào Công Khanh (1996) dùng hàm Schumacher; Vũ Đình Phương (1998); Trần Cẩm Tú (1999); Nguyễn Thành Mến (2005).
    • Phân nhóm loài và mô hình sinh trưởng: Phân loại truyền thống theo 3 nhóm sinh thái (ưa sáng, chịu bóng, trung tính), 3 nhóm tốc độ sinh trưởng hoặc 8 nhóm gỗ thương mại theo Quyết định của Bộ Lâm nghiệp (1977); phân nhóm theo hình số thân cây (Đồng Sỹ Hiền, 1974; Vũ Tiến Hinh, 2012). Đỗ Đình Sâm và cộng sự (2007) lập phương trình tăng trưởng đường kính theo vùng sinh thái. Trần Văn Con và cộng sự (2006, 2010) phân 20 loài cây tại Tây Nguyên thành 3 nhóm tăng trưởng dựa trên OTCĐV. Vũ Tiến Hinh và cộng sự (1987, 2010, 2011) mô phỏng động thái $N/D$ và tăng trưởng trữ lượng $Z_M$. Trần Văn Con (1991) xây dựng hàm tăng trưởng chiều cao rừng Khộp $h = e^a (S_i / e^a)(t_g / t)^c$ và suất tăng trưởng đường kính $P_d = -0{,}01044 + 0{,}1055 d^{-0{,}5}$ ($R^2 = 0{,}996$). Bảo Huy (1993) xây dựng quan hệ $M = 14{,}88416 G^{1{,}07126} S^{0{,}28900} IV^{0{,}10763}$ và hàm $Z_M$. Trịnh Khắc Mười (1993) nghiên cứu tương quan suất tăng trưởng thể tích và đường kính.

  • Khoảng trống nghiên cứu được lựa chọn: Hầu hết các nghiên cứu tăng trưởng rừng tự nhiên ở Việt Nam trước đây áp dụng phương pháp tạm thời (giải tích thớt gỗ, đẽo vát). Luận án lựa chọn giải quyết khoảng trống về thiếu hụt dữ liệu đo lặp định kỳ trên hệ sinh thái ô tiêu chuẩn định vị dài hạn, xây dựng mô hình toán học liên hoàn tích hợp 3 quá trình động thái: tăng trưởng đường kính (kết hợp vị thế tán cây và biến giả vùng), tỷ lệ chết và tái sinh bổ sung nhằm mô phỏng động thái cấu trúc lâm phần rừng đặc dụng miền Bắc.

Cơ sở lý thuyết và phương pháp nghiên cứu

  • Khung lý thuyết và khái niệm sử dụng:

    • Lý thuyết động thái cấu trúc lâm phần (Vanclay, 1994): Cấu trúc đường kính lâm phần biến đổi qua 4 quá trình: $$N_{k, t+1} = N_{k, t} + I_k - O_k - M_k - H_k$$ Hoặc chi tiết theo cây tái sinh bổ sung vào tầng cây cao ($R_k$) và không có khai thác ($H_k = 0$ đối với rừng đặc dụng): $$N_{k, t+1} = N_{k, t} + R_k + O_{k-1} - O_k - M_k$$ Trong đó: $N_{k,t}$ là số cây ở cỡ kính $k$ tại thời điểm $t$; $I_k$ là số cây chuyển vào; $O_k$ là số cây chuyển ra; $M_k$ là số cây chết; $R_k$ là số cây tái sinh bổ sung từ tầng cây tái sinh vào cấp kính nhỏ nhất của tầng cây cao ($D_{1.3} \ge 10\text{ cm}$). Dạng tỷ lệ phần trăm: $$N_{k, t+1} = R_k + N_{k-1, t} \cdot o_{k-1} + N_{k, t}(1 - o_k - m_k - h_k)$$
    • Phân cấp vị thế tán cây (Dawkins, 1958): Vị thế tầng tán ($PC$) được phân thành các cấp phản ánh mức độ nhận ánh sáng và cạnh tranh không gian dinh dưỡng.
    • Chỉ số tổ thành và đa dạng loài: Chỉ số quan trọng $IV%$ (Daniel Marmillod), hệ số hỗn loài $HL$, chỉ số Shannon - Wiener ($H'$), chỉ số Simpson/Gini ($D$), và họ chỉ số đa dạng Rényi ($H_\alpha$).

  • Phương pháp thu thập dữ liệu thực địa:

    • Thiết lập và điều tra trên hệ thống 21 OTCĐV cố định có diện tích 1 ha ($100\text{ m} \times 100\text{ m}$).
    • Mỗi ô 1 ha được chia thành 25 ô vuông phụ kích thước $20\text{ m} \times 20\text{ m}$.
    • Đo đếm toàn bộ cây thân gỗ tầng cao ($D_{1.3} \ge 10\text{ cm}$), đánh số cố định, xác định tọa độ tương đối $(x, y)$ trong ô $20\text{ m} \times 20\text{ m}$, đo đường kính ngang ngực tại vị trí 1,3 m ($D_{1.3}$), đo chiều cao vút ngọn ($H_{vn}$), xác định tên loài và phân cấp vị thế tầng tán theo Dawkins ($PC$).
    • Thu thập dữ liệu lặp lại giữa 2 kỳ đo đếm năm 2007 và năm 2012 để xác định lượng tăng trưởng thực tế ($z_d$), số cây chết ($M$) và số cây tái sinh bổ sung ($R$).

  • Phương pháp xử lý và phân tích số liệu:

    • Phân nhóm loài: Sử dụng phân tích cụm (Cluster Analysis) theo thuật toán K-Means kết hợp phân tích biệt thức (Discriminant Analysis) dựa trên các chỉ tiêu: lượng tăng trưởng đường kính bình quân năm ($z_d$), đường kính tối đa loài đạt được ($D_{max}$) cho 25 loài ưu thế có dung lượng mẫu quan sát $n \ge 50$, kết hợp với 3 dạng sống: gỗ lớn (GOL), gỗ trung bình (GOT), gỗ nhỏ (GON).
    • Mô hình hóa cấu trúc tĩnh: Kiểm nghiệm phân bố $N/D_{1.3}$ bằng hàm Khoảng cách và hàm Weibull; kiểm nghiệm tương quan $H_{vn}/D_{1.3}$ bằng các dạng hàm toán học.
    • Mô hình hóa tăng trưởng đường kính: Tuyến tính hóa phương trình tăng trưởng phi tuyến qua hàm logarit: $$\text{Ln}\left[\frac{z_d}{D_{max} - D}\right] = a + b \cdot \text{Ln}(D)$$ Mở rộng tích hợp vị thế tán cây ($PC$): $$\text{Ln}\left[\frac{z_d}{D_{max} - D}\right] = a + b \cdot \text{Ln}(D) + c \cdot PC$$ Mô hình hóa $D_{max}$ theo vị thế tán: $D_{max} = k + m \cdot PC$, chuyển đổi về phương trình phi tuyến gốc: $$z_d = \exp(a) \cdot \exp(c \cdot PC) \cdot [(k + m \cdot PC) - D] \cdot D^b$$ Đồng thời đưa biến giả vùng ($S_i$) vào phương trình: $$\text{Ln}\left[\frac{z_d}{D_{max} - D}\right] = a + b \cdot \text{Ln}(D) + c \cdot PC + d_i \cdot S_i$$
    • Đánh giá mô hình: Sử dụng tổng bình phương sai lệch ($RSS$), sai tiêu chuẩn hồi quy ($SE$), hệ số xác định ($R^2$), và kiểm định xác suất tiêu chuẩn ($Sig. f$).
    • Mô hình hóa quá trình chết và tái sinh bổ sung: Sử dụng hàm toán học (như hàm chữ S) mô phỏng số cây chết theo cỡ kính và hàm tái sinh bổ sung.

Nội dung chính theo từng chương

Chương 1: Tổng quan các vấn đề nghiên cứu

Chương này tổng thuật toàn diện các cơ sở lý luận và thực nghiệm trong và ngoài nước về cấu trúc lâm phần, phân nhóm sinh thái - sinh trưởng loài cây và các mô hình động thái rừng nhiệt đới. Tác giả phân tích các hạn chế của các mô hình mật độ giới hạn (Reineke, Yoda, Wilson) khi áp dụng cho rừng tự nhiên hỗn loài khác tuổi; tổng hợp các nguyên nhân gây chết tự nhiên (thiếu sáng, sâu bệnh, phơi sáng sau tạo lỗ trống, đổ gãy cơ học do gió bão) và biên độ tỷ lệ chết dao động 1 - 5%/năm ở vùng nhiệt đới. Chương 1 khẳng định việc thiếu vắng dữ liệu theo dõi định vị dài hạn tại Việt Nam là rào cản chính cho việc xây dựng mô hình dự báo động thái, từ đó xác lập các nội dung nghiên cứu trọng tâm của luận án.

Chương 2: Giới hạn, nội dung và phương pháp nghiên cứu

Chương 2 trình bày chi tiết thiết kế nghiên cứu thực nghiệm tại 4 khu rừng đặc dụng miền Bắc (VQG Ba Bể, VQG Vũ Quang, VQG Xuân Sơn, KBT Hang Kia - Pà Cò) với hệ thống 21 OTCĐV cố định 1 ha. Tác giả chuẩn hóa quy trình điều tra cây đứng tầng cao ($D_{1.3} \ge 10\text{ cm}$), chia ô phụ $20\text{ m} \times 20\text{ m}$, kỹ thuật đo tọa độ tương đối, xác định điểm đo $D_{1.3}$, phân cấp tán Dawkins ($PC$), và phương pháp xử lý dữ liệu đo lặp 5 năm (2007 - 2012). Chương này thiết lập hệ thống công thức toán học phân tích đa dạng sinh học ($IV%$, $HL$, $H'$, $D$, Rényi), kỹ thuật phân cụm K-Means, phân tích biệt thức và cấu trúc các phương trình hồi quy tăng trưởng có biến giả vị trí địa lý ($S_i$).

Chương 3: Kết quả và thảo luận

Chương 3 trình bày các kết quả định lượng cụ thể qua 4 nội dung chính:

  • Đặc điểm cấu trúc lâm phần và đa dạng loài:
    • Kết quả chỉ số quan trọng $IV%$ xác định các loài ưu thế tại các khu vực điều tra năm 2012.
    • Tính đa dạng loài thực vật được lượng hóa chi tiết qua hệ số hỗn loài $HL$, chỉ số $H'$, chỉ số $D$ và họ đường cong chỉ số Rényi ($H_\alpha$).
    • Mô phỏng cấu trúc $N/D_{1.3}$: Hàm Weibull (với tham số hình dạng $c \le 1$, phản ánh phân bố giảm liên tục đặc trưng của rừng nhiệt đới tự nhiên) và hàm Khoảng cách được kiểm nghiệm và mô phỏng phù hợp với phân bố thực nghiệm tại các OTCĐV.
    • Tương quan $H_{vn}/D_{1.3}$: Lựa chọn dạng phương trình tương quan tối ưu cho từng khu vực nghiên cứu (Ba Bể, Vũ Quang, Xuân Sơn, Hang Kia - Pà Cò).
  • Phân nhóm loài cây:
    • Áp dụng thuật toán K-Means trên tập dữ liệu 25 loài ưu thế ($n \ge 50$) với 2 biến $D_{max}$ và $z_d$, kiểm định qua phân tích biệt thức, sau đó tích hợp với 3 dạng sống (GOL, GOT, GON) để xác lập hệ thống 9 nhóm loài chức năng phục vụ mô hình hóa tăng trưởng cho toàn bộ lâm phần.
  • Xây dựng mô hình tăng trưởng đường kính, quá trình chết và tái sinh bổ sung:
    • Mô hình tăng trưởng đường kính: Ước lượng các tham số $a, b, c, d_i$ cho phương trình tăng trưởng đường kính bình quân năm ($z_d$). Việc bổ sung biến vị thế tầng tán ($PC$) và biến giả vùng ($S_i$) giúp giảm đáng kể tổng sai lệch bình phương ($RSS$) và sai tiêu chuẩn ($SE$) so với dạng phương trình chỉ chứa biến $D$.
    • Mô hình quá trình chết: Số cây chết ($M_{bb}, M_{hk}, M_{vq}, M_{xs}$) theo các cỡ đường kính được mô phỏng phù hợp nhất qua dạng hàm chữ S.
    • Mô hình tái sinh bổ sung: Xác lập phương trình hàm mô phỏng số lượng cây tái sinh bổ sung ($RTSBS$) vào cấp kính $10\text{ cm}$ cho các khu rừng đặc dụng.
  • Mô phỏng động thái cấu trúc lâm phần:
    • Tích hợp các hàm thành phần (tăng trưởng đường kính $z_d$, cây chuyển cấp $O$, cây tái sinh bổ sung $R$, cây chết $M$) để dự đoán động thái cấu trúc $N/D$, cấu trúc tiết diện ngang $G/D$ và tăng trưởng trữ lượng lâm phần $Z_M$ ($\text{m}^3/\text{ha}$) theo định kỳ thời gian tại VQG Ba Bể.

Kết quả và những đóng góp mới

  • Những đóng góp mới của luận án:

    1. Là công trình nghiên cứu có hệ thống đầu tiên tại Việt Nam về đồng thời 3 quá trình động thái: tăng trưởng đường kính, quá trình chết và quá trình tái sinh bổ sung trong rừng tự nhiên lá rộng thường xanh dựa trên mạng lưới ô tiêu chuẩn định vị đo lặp.
    2. Cung cấp hệ thống dữ liệu khoa học định lượng và các phương trình thực nghiệm phục vụ dự đoán tăng trưởng, diễn biến trữ lượng và động thái cấu trúc lâm phần, làm cơ sở khoa học cho quản lý rừng bền vững.
    3. Góp phần hoàn thiện phương pháp luận về phân nhóm loài thực vật (dựa trên $z_d$, $D_{max}$ và dạng sống), xác định tỷ lệ tái sinh bổ sung và chuyển cấp; áp dụng phương pháp nghiên cứu tăng trưởng đường kính không thông qua giải tích thân cây (chặt hạ cây mẫu) trong rừng tự nhiên lá rộng thường xanh.

  • Khuyến nghị và giải pháp đề xuất:

    • Vận dụng các mô hình động thái cấu trúc và phương trình tăng trưởng để dự báo xu thế phát triển tài nguyên rừng tại 4 khu rừng đặc dụng nghiên cứu.
    • Ứng dụng phương pháp luận và các mô hình sinh trưởng vào công tác quy hoạch, điều chế rừng và quản lý kinh doanh tại các trạng thái rừng sản xuất tự nhiên lá rộng thường xanh có điều kiện lập địa và cấu trúc tương đồng.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế:

    • Chu kỳ theo dõi số liệu thực nghiệm trên hệ thống OTCĐV mới dừng lại ở 1 định kỳ 5 năm (2007 - 2012).
    • Phạm vi không gian nghiên cứu tập trung tại 4 khu rừng đặc dụng thuộc miền Bắc và Bắc Trung Bộ (Hà Tĩnh), chưa bao quát toàn bộ các tiểu vùng sinh thái rừng tự nhiên khác trên toàn quốc.
    • Tầng cây khảo sát giới hạn ở cây gỗ có $D_{1.3} \ge 10\text{ cm}$, các thông tin về tầng cây tái sinh dưới $10\text{ cm}$ chủ yếu được tích hợp dưới dạng số lượng cá thể bổ sung vào tầng cây cao.

  • Hướng nghiên cứu tiếp:

    • Tiếp tục duy trì bảo vệ và tổ chức đo đếm định kỳ lặp lại ở các chu kỳ tiếp theo (10 năm, 15 năm, 20 năm) trên hệ thống 21 OTCĐV để hoàn thiện và hiệu chỉnh độ tin cậy của các tham số mô hình.
    • Mở rộng ứng dụng thử nghiệm mô hình hóa động thái trên các kiểu trạng thái rừng tự nhiên phục hồi sau khai thác và rừng sản xuất.

Giá trị tham khảo

  • Đối tượng tham khảo:
    • Nghiên cứu sinh, học viên cao học và sinh viên chuyên ngành Lâm sinh, Quản lý tài nguyên rừng, Điều tra quy hoạch rừng: Tham khảo phương pháp thiết lập OTCĐV 1 ha, kỹ thuật phân cụm loài K-Means kết hợp phân tích biệt thức và phương pháp mô hình hóa phi tuyến không giải tích.
    • Giảng viên các trường đại học khối nông lâm nghiệp: Bổ sung tư liệu thực tiễn và phương pháp luận cho giáo trình Điều tra rừng, Sản lượng rừng và Lâm sinh học.
    • Cán bộ quản lý tại các Vườn Quốc gia, Khu Bảo tồn thiên nhiên và các Ban Quản lý rừng phòng hộ/sản xuất: Ứng dụng các bảng tra tương quan $H/D$, hàm dự báo cấu trúc $N/D$, động thái $G/D$ và $Z_M$ để xây dựng phương án quản lý rừng bền vững.

Câu hỏi thường gặp

  • Câu hỏi 1: Đối tượng khảo sát chính của luận án là gì và được thu thập tại các địa bàn nào?
    Trả lời: Đối tượng nghiên cứu là tầng cây cao thân gỗ có đường kính ngang ngực $D_{1.3} \ge 10\text{ cm}$ trong rừng tự nhiên lá rộng thường xanh, được thu thập trên 21 ô tiêu chuẩn định vị cố định (mỗi ô 1 ha) tại 4 khu rừng đặc dụng: VQG Ba Bể (6 ô), VQG Vũ Quang (6 ô), VQG Xuân Sơn (3 ô) và KBT Hang Kia - Pà Cò (6 ô).

  • Câu hỏi 2: Phân nhóm loài cây trong luận án được thực hiện dựa trên các chỉ tiêu định lượng nào?
    Trả lời: Tác giả sử dụng phân tích cụm K-Means và phân tích biệt thức dựa trên 2 chỉ tiêu sinh trưởng chính là lượng tăng trưởng đường kính bình quân năm ($z_d$) và đường kính tối đa loài đạt được ($D_{max}$) trên 25 loài ưu thế ($n \ge 50$), sau đó kết hợp với 3 dạng sống (gỗ lớn - GOL, gỗ trung bình - GOT, gỗ nhỏ - GON) để xác lập hệ thống 9 nhóm loài cho lâm phần.

  • Câu hỏi 3: Biến số vị thế tầng tán cây ($PC$) được phân cấp theo cơ sở nào và đóng vai trò gì trong mô hình tăng trưởng đường kính?
    Trả lời: Vị thế tầng tán ($PC$) được phân cấp theo phương pháp của Dawkins (1958). Việc tích hợp biến $PC$ vào phương trình tăng trưởng $\text{Ln}[z_d/(D_{max} - D)] = a + b\text{Ln}(D) + cPC$ giúp mô hình phản ánh mức độ cạnh tranh không gian dinh dưỡng và ánh sáng của từng cá thể, làm giảm rõ rệt tổng sai lệch bình phương ($RSS$) và phương sai hồi quy ($SE$).

  • Câu hỏi 4: Quá trình chết và phân bố số cây theo đường kính trong lâm phần được mô phỏng bằng các dạng hàm toán học nào?
    Trả lời: Cấu trúc phân bố số cây theo cỡ kính ($N/D_{1.3}$) được mô phỏng bằng hàm Weibull (với tham số $c \le 1$) và hàm Khoảng cách; trong khi đó, số cây chết ($M$) theo cỡ đường kính trong chu kỳ 5 năm được mô phỏng phù hợp nhất qua dạng hàm chữ S.

  • Câu hỏi 5: Phương trình cân bằng cơ bản nào được áp dụng để mô phỏng động thái cấu trúc lâm phần qua các chu kỳ?
    Trả lời: Động thái thay đổi số cây theo cỡ kính được mô phỏng qua phương trình cân bằng: $N_{k, t+1} = N_{k, t} + R_k + O_{k-1} - O_k - M_k - H_k$, trong đó phản ánh số cây hiện tại ($N_{k,t}$), số cây tái sinh bổ sung ($R_k$), cây chuyển vào từ cấp trước ($O_{k-1}$), cây chuyển ra cấp sau ($O_k$), cây chết ($M_k$) và cây khai thác ($H_k$).

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Thị Thu Hiền đã giải quyết một vấn đề khoa học then chốt trong lâm sinh học thông qua việc định lượng hóa các quy luật cấu trúc và quá trình động thái của rừng tự nhiên lá rộng thường xanh tại 4 khu rừng đặc dụng miền Bắc Việt Nam. Dựa trên nguồn dữ liệu đo lặp 5 năm (2007 - 2012) từ mạng lưới 21 ô tiêu chuẩn định vị 1 ha, nghiên cứu đã xây dựng thành công hệ thống mô hình toán học tích hợp tăng trưởng đường kính cá thể, quy luật chết tự nhiên và tái sinh bổ sung. Kết quả nghiên cứu cung cấp cơ sở phương pháp luận hoàn chỉnh về phân nhóm loài và mô phỏng sinh trưởng không qua giải tích, đóng góp trực tiếp vào công tác dự báo tài nguyên và quản lý rừng bền vững tại Việt Nam.