Tổng quan về luận án

Nghiên cứu cấu trúc lâm phần, quy luật sinh trưởng và tích lũy sinh khối của các loài cây gỗ quý là trọng tâm của khoa học lâm sinh hiện đại trong bối cảnh biến đổi khí hậu toàn cầu và suy giảm tài nguyên rừng tự nhiên. Luận án tiến sĩ lâm nghiệp của tác giả Nguyễn Công Hoan (2014) với tiêu đề "Nghiên cứu sinh trưởng và cấu trúc làm cơ sở khoa học đề xuất biện pháp kỹ thuật góp phần kinh doanh rừng trồng Tếch (Tectona grandis L.) tại Sơn La" thuộc chuyên ngành Lâm sinh (Mã số: 62620205), dưới sự hướng dẫn khoa học của GS. TS. Vũ Tiến Hinh và PGS. TS. Nguyễn Văn Sinh tại Đại học Thái Nguyên, đại diện cho một công trình nghiên cứu hệ thống và tiên phong tại khu vực miền núi phía Bắc Việt Nam.

                  ┌─────────────────────────────────────────────────────────┐
                  │    HỆ SINH THÁI RỪNG TRỒNG TẾCH (TECTONA GRANDIS L.)    │
                  │             TẠI SƠN LA (TUỔI 13 - CHƯA TỈA THƯA)         │
                  └────────────────────────────┬────────────────────────────┘
                                               │
                        ┌──────────────────────┴──────────────────────┐
                        ▼                                             ▼
       ┌─────────────────────────────────┐           ┌─────────────────────────────────┐
       │   QUY LUẬT CẤU TRÚC & SINH TRƯỞNG│           │   SINH KHỐI & DỰ TRỮ CÁC-BON    │
       ├─────────────────────────────────┤           ├─────────────────────────────────┤
       │• Phân bố N/D1.3, N/Hvn: Weibull │           │• Cây gỗ: 36 cây mẫu giải tích   │
       │• Mô hình tương quan Hvn-D1.3,   │           │• Thảm tươi (W2k), thảm mục (W3k)│
       │  Dt-D1.3: Tuyến tính & Phi tuyến│           │• Hệ số chuyển đổi C: 0,485      │
       │• Động thái sinh trưởng D, H, V: │           │• Định lượng C theo Cấp đất      │
       │  Mô hình Schumacher             │           │  (Cấp II & Cấp III)             │
       └────────────────┬────────────────┘           └────────────────┬────────────────┘
                        │                                             │
                        └──────────────────────┬──────────────────────┘
                                               │
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │   HỆ THỐNG GIẢI PHÁP LÂM SINH & ĐIỀU CHẾ RỪNG KINH DOANH │
                  ├─────────────────────────────────────────────────────────┤
                  │• Phân cấp lập địa chỉ thị chiều cao (Cấp đất II & III)  │
                  │• Biểu điều tiết mật độ, cường độ tỉa thưa (Nopt, Nc, Nnd│
                  │• Tích hợp cơ chế chi trả dịch vụ các-bon (REDD+/CDM)    │
                  └─────────────────────────────────────────────────────────┘

Theo đánh giá lịch sử phát triển, cây Tếch (Tectona grandis L.) - loài cây gỗ lớn thuộc họ Tếch (Verbenaceae), bộ Hoa môi (Lamiales), có nguồn gốc tự nhiên từ khu hệ thực vật Ấn Độ - Miến Điện - đã được du nhập và thử nghiệm thành công tại miền Nam Việt Nam từ đầu thế kỷ XX (Định Quán, Đồng Nai năm 1960 với quy mô trên 200 ha; Đắk Lắk tại Eak Mát). Tại vùng Tây Bắc, cụ thể là tỉnh Sơn La, Tếch được đưa vào gây trồng tập trung trong khuôn khổ Chương trình 327 và dự án Hợp tác Kỹ thuật Đức (GTZ) nhằm mục tiêu kép: phủ xanh đất trống đồi núi trọc, giảm áp lực khai thác gỗ tự nhiên và tạo sinh kế bền vững thông qua cung cấp gỗ hộp mộc cao cấp.

Tuy nhiên, tồn tại một khoảng trống nghiên cứu (research gap) mang tính cốt lõi: các nghiên cứu trước đây tại Việt Nam về cây Tếch (Bảo Huy, 1993, 1995; Nguyễn Xuân Quát, 1995; Nguyễn Ngọc Lung, 1988, 1995; Nguyễn Văn Thêm, 2002; Mạc Văn Chăm, 2005) chỉ tập trung tại Tây Nguyên và Nam Trung Bộ, chủ yếu tiếp cận thuần túy dưới góc độ sản lượng gỗ thương phẩm. Chưa có bất kỳ công trình nào lượng hóa đồng thời quy luật phân bố không gian, động thái tăng trưởng đường kính - chiều cao - thể tích kết hợp với đánh giá trữ lượng sinh khối và khả năng hấp thụ các-bon theo cấp đất tại vùng sinh thái Tây Bắc.

Luận án đã giải quyết tường minh ba câu hỏi nghiên cứu trung tâm:

  1. Quy luật cấu trúc lâm phần (phân bố số cây theo cỡ kính $N/D_{1.3}$, chiều cao $N/H_{vn}$, tương quan $H_{vn}/D_{1.3}$ và đường kính tán $D_t/D_{1.3}$) của rừng trồng Tếch tại Sơn La biến động như thế nào theo cấp đất?
  2. Mô hình toán học phi tuyến tính nào mô phỏng chính xác nhất tiến trình sinh trưởng đường kính ($D$), chiều cao ($H$) và thể tích ($V$) của lâm phần Tếch thuần loài tuổi 13?
  3. Trữ lượng sinh khối khô và dung lượng tích lũy các-bon ($C$) phân bố như thế nào giữa các bộ phận cây cá lẻ (thân, cành, lá, vỏ) và các tầng cấu trúc rừng (cây gỗ, cây bụi thảm tươi, thảm mục) trên từng cấp đất?

Nghiên cứu được triển khai trên phạm vi 54 lâm phần thực nghiệm (tương ứng 54 ô tiêu chuẩn định lượng) tại hai huyện trọng điểm Mai Sơn và Yên Châu thuộc tỉnh Sơn La, khảo sát chi tiết 36 cây tiêu chuẩn chặt ngả phục vụ giải tích thân cây và sinh khối. Nghiên cứu xác lập cơ sở khoa học thực chứng để chuyển đổi mô hình quản trị rừng Tếch từ quảng canh sang thâm canh đa mục tiêu, định lượng hóa giá trị dịch vụ hệ sinh thái các-bon theo cơ chế CDM và REDD+.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước phản ánh sự tiến hóa rõ rệt trong phương pháp luận điều tra lâm phần và lượng hóa hệ sinh thái rừng. Trong phân tích cấu trúc đường kính ($N/D_{1.3}$), xuất phát điểm từ mô hình đường cong giảm liên tục của Meyer (1952) áp dụng cho rừng tự nhiên bất đều tuổi, các nhà khoa học lâm nghiệp đã chuyển dần sang các phân bố xác suất tham số phức hợp cho rừng trồng đều tuổi. Bailey và Dell (1973) chứng minh tính ưu việt của phân bố Weibull 3 tham số (vị trí $a$, tỷ lệ $b$, dạng $c$) nhờ khả năng phản ánh linh hoạt độ lệch và độ nhọn của lâm phần qua các giai đoạn sinh trưởng. Quan điểm này được củng cố bởi Garcia (1980) khi khuyến cáo hạn chế sai số ở cỡ kính nhỏ, cũng như nghiên cứu của Suzuki (1971) và Preussner (1986) về ứng dụng chuỗi Markov để mô hình hóa sự chuyển cấp kính dưới tác động của sinh trưởng và tỉa thưa cơ giới.

Về quy luật tương quan $H_{vn}/D_{1.3}$, lý thuyết đường cong chiều cao cơ sở của Prodan (1963) và các dạng phương trình hồi quy của Naslund (1967), Curtis (1967) chỉ ra rằng tương quan chiều cao - đường kính mang tính đặc trưng loài và chịu sự chi phối chặt chẽ của cấp đất lập địa.

Tại Việt Nam, các nghiên cứu của Đồng Sỹ Hiền (1974), Nguyễn Hải Tuất (1982), Vũ Tiến Hinh (1990, 2000), Lê Hồng Phúc (1996) cho Thông ba lá, Hoàng Văn Dưỡng (2001) cho Keo lá tràm đều khẳng định hàm Weibull và hàm Logarit ($H = a + b \cdot \ln(D_{1.3})$) mô tả tối ưu cấu trúc rừng trồng.

Tác giả & Năm Đối tượng & Địa bàn Mô hình / Phương pháp Phát hiện & Giới hạn cốt lõi
Kraenzel và cs. (2003); Schroeder và cs. (2001) Rừng Tếch 20 tuổi tại Panama Phương pháp giải tích sinh khối tầng cây gỗ & đất Tổng trữ lượng $C$ đạt $351\text{ tấn/ha}$ (cây gỗ $120,2\text{ tấn/ha}$, đất $225\text{ tấn/ha}$); chưa mô hình hóa theo cấp đất chi tiết.
Perez và cs. (2003); Bermejo và cs. (2003) Rừng Tếch thuần loài tại Costa Rica Phương pháp chặt cây mẫu, cân trọng lượng khô Sinh khối khô trên mặt đất dao động $84 - 284\text{ Mg/ha}$ (thân chiếm $70 - 90%$, cành $5 - 30%$, lá $1 - 6%$); tập trung vào động thái tuổi.
Ly Meng Seang (2009) Rừng Tếch tại Kampong Cham, Campuchia Phân bố Weibull, tương quan sinh trưởng Sinh trưởng nhanh giai đoạn 1 - 7 tuổi, chậm dần từ tuổi 8 - 18; năng suất Cấp I gấp 1,24 lần Cấp II; thiếu dữ liệu sinh khối và các-bon.
Bảo Huy (1995, 1998) Rừng Tếch tại Đắk Lắk, Tây Nguyên Mô hình phân chia cấp đất, biểu thể tích Thiết lập biểu sản lượng và cấp đất lâm phần cho Tây Nguyên; chưa tiếp cận định lượng sinh khối và các-bon.
Nguyễn Công Hoan (2014) [Luận án] Rừng Tếch 13 tuổi tại Sơn La, Tây Bắc Tích hợp Weibull, Schumacher, allometric $W_k$, hệ số chuyển đổi $0,485$ Định lượng toàn diện cấu trúc, sinh trưởng, sinh khối khô và tích lũy $C$ theo Cấp đất II và III; xây dựng biểu tỉa thưa định lượng.

So với các công trình quốc tế, nghiên cứu của Nguyễn Công Hoan đã định vị chính xác vị trí học thuật bằng việc tích hợp lý thuyết sinh thái học cá thể với giải tích điều tra lâm phần thực nghiệm tại vùng sinh thái đặc thù Tây Bắc Việt Nam - nơi có biên độ nhiệt mùa đông sâu và mùa khô kéo dài, tạo nên những biến thiên sinh trưởng khác biệt so với vùng xích đạo nhiệt đới ẩm thuần nhất như Panama hay Costa Rica.

Đó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 các quy luật kinh điển trong sinh thái học rừng và điều tra rừng:

  1. Kiểm chứng quy luật Eichhorn (1904): Luận án đã củng cố định luật sinh trắc học kinh điển: "Trữ lượng rừng là một hàm số của chiều cao bình quân lâm phần", chứng minh rằng chiều cao ưu thế ($H_0$) và chiều cao tầng trội là chỉ số chỉ thị cấp đất độc lập với mật độ ban đầu và biện pháp tỉa thưa cơ giới ở giai đoạn rừng non và sào.
  2. Phát triển lý thuyết hàm sản lượng sinh học: Mở rộng phương trình tăng trưởng vi phân của Mitscherlich (1919), Richards (1959) và Gompertz (1825/1985) thông qua việc chứng minh tính vượt trội của hàm Schumacher trong việc phản ánh điểm uốn sinh trưởng đường kính và thể tích của loài Tectona grandis L. tại điều kiện lập địa khắc nghiệt miền núi cao.
  3. Mở rộng định luật tối thiểu của Liebig (1840) và lý thuyết năng suất quần xã của Whittaker (1970): Khẳng định câu châm ngôn học thuật: "Số đo năng suất chính là số đo về tăng trưởng, tích lũy sinh khối ở cơ thể thực vật trong quần xã", chuyển dịch trọng tâm từ đánh giá thể tích thân gỗ thương phẩm thuần túy sang đánh giá toàn diện năng suất sơ cấp thực (NPP) và chu trình tích lũy các-bon đa tầng.

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

Khung phân tích của đề tài được xây dựng trên sự tích hợp đa tầng giữa ba lý thuyết nền tảng: Lý thuyết cấu trúc lâm phần động thái, Lý thuyết không gian dinh dưỡng của Phùng Ngọc Lan (1986), và Hướng dẫn kiểm kê khí nhà kính của IPCC (hệ số chuyển đổi $C = 0,485$).

                ┌─────────────────────────────────────────────────────────────┐
                │             KHUNG PHÂN TÍCH HỆ THỐNG ĐA TẦNG                │
                └──────────────────────────────┬──────────────────────────────┘
                                               │
             ┌─────────────────────────────────┼─────────────────────────────────┐
             ▼                                 ▼                                 ▼
┌─────────────────────────┐       ┌─────────────────────────┐       ┌─────────────────────────┐
│   TẦNG LẬP ĐỊA ĐẤT ĐAI  │       │ TẦNG CẤU TRÚC QUẦN THỂ  │       │  TẦNG SINH KHỐI CÁC-BON │
├─────────────────────────┤       ├─────────────────────────┤       ├─────────────────────────┤
│• Địa hình: Mai Sơn,     │       │• Động thái mật độ N/ha  │       │• Cây gỗ: Thân, cành, lá │
│  Yên Châu               │──────>│• Phân bố N/D1.3, N/Hvn  │──────>│• Tầng phụ trợ: Thảm tươi│
│• Cấp đất: Cấp II và     │       │• Tương quan Dt - D1.3   │       │  (W2k) & Thảm mục (W3k) │
│  Cấp III                │       │• Phân cấp sinh trưởng   │       │• Hệ số quy đổi C: 0,485 │
└─────────────────────────┘       └─────────────────────────┘       └─────────────────────────┘
             │                                 │                                 │
             └─────────────────────────────────┼─────────────────────────────────┘
                                               │
                                               ▼
                ┌─────────────────────────────────────────────────────────────┐
                │          MÔ HÌNH TOÁN HỌC & ĐIỀU TIẾT LÂM SINH              │
                │     Hàm Schumacher ──> Điểm cắt Z & Δ ──> Biểu tỉa thưa     │
                └─────────────────────────────────────────────────────────────┘

Mô hình phân tích xác lập các quan hệ hàm số tường minh:

  • Mô hình phân cấp sinh trưởng cá thể dựa trên phương trình phân loại nhiều biến số nhằm tách biệt cây sinh trưởng tốt (loại I), trung bình (loại II) và xấu (loại III).
  • Mô hình tương quan không gian dinh dưỡng giữa đường kính tán ($D_t$) và đường kính thân ($D_{1.3}$): $D_t = a + b \cdot D_{1.3}$, làm căn cứ thiết lập diện tích tán bình quân ($S_t$) và mật độ nuôi dưỡng tối ưu ($N_{opt}$).
  • Điều kiện biên của mô hình: Áp dụng nghiêm ngặt cho rừng trồng Tếch thuần loài, đều tuổi, giai đoạn khép tán đến tuổi 13 trên hai nhóm cấp đất chính (Cấp II và Cấp III) tại vùng đồi núi có mùa khô rõ rệt.

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

Thiết kế nghiên cứu

Đề tài áp dụng thế giới quan thực chứng khách quan (positivism) với phương pháp luận định lượng quy chuẩn. Thiết kế nghiên cứu đa cấp độ (multi-level design) kết hợp giữa điều tra thống kê quần xã lâm phần và phương pháp thực nghiệm phá hủy mẫu (destructive sampling) theo tiêu chuẩn của Chương trình Sinh học Quốc tế (IBP) do Newbould (1967) đề xuất.

Cỡ mẫu nghiên cứu được thiết kế với độ tin cậy thống kê cao:

  • Tầng lâm phần: Thiết lập 54 ô tiêu chuẩn ($500\text{ m}^2/\text{ô}$, kích thước $20\text{ m} \times 25\text{ m}$ hoặc $25\text{ m} \times 25\text{ m}$) phủ đều trên 54 lâm phần Tếch tuổi 13 tại huyện Mai Sơn và Yên Châu.
  • Tầng cây tiêu chuẩn: Chọn và chặt hạ 36 cây tiêu chuẩn phân bổ đều theo các cấp đường kính và cấp đất để giải tích thân cây, xác định sinh trưởng thể tích và đo đếm khối lượng tươi/khô của các bộ phận riêng biệt.
  • Tầng cây bụi, thảm tươi và thảm mục: Trong mỗi ô tiêu chuẩn bố trí 5 ô phụ diện tích $1\text{ m}^2$ (4 góc và tâm ô) để thu gom toàn bộ thảm tươi ($W_2$) và thảm mục rơi rụng ($W_3$).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý mẫu tuân thủ nghiêm ngặt các tiêu chuẩn lâm sinh học:

  1. Đo đếm lâm phần: Toàn bộ cây trong 54 ô tiêu chuẩn được đo đường kính ngang ngực ($D_{1.3}$) bằng thước kẹp kính chính xác đến $0,1\text{ cm}$, đo chiều cao vút ngọn ($H_{vn}$) và chiều cao dưới cành ($H_{dc}$) bằng thước đo cao Blume-Leiss chính xác đến $0,5\text{ m}$, đo 2 chiều vuông góc của đường kính tán ($D_t$).
  2. Quy trình giải tích thân cây và sinh khối: 36 cây tiêu chuẩn được chặt sát gốc; đo chiều dài toàn thân, cắt khúc $1\text{ m}$ hoặc $2\text{ m}$ để giải tích đĩa gỗ. Phân chia cây thành 4 bộ phận: thân gỗ, vỏ, cành và lá. Cân toàn bộ trọng lượng tươi ($W_t$) ngay tại hiện trường bằng cân đồng hồ chuyên dụng sai số $\pm 0,05\text{ kg}$.
  3. Quy trình sấy khô phòng thí nghiệm: Lấy mẫu đại diện từng bộ phận (mẫu đĩa thân, mẫu vỏ, mẫu cành các cấp, mẫu lá tổng hợp), cân trọng lượng tươi mẫu ($P_{ti}$), sau đó sấy khô kiệt trong tủ sấy công nghiệp ở nhiệt độ $105^\circ\text{C}$ đến khi khối lượng không đổi để xác định khối lượng khô kiệt ($P_{ki}$). Tỷ lệ sinh khối khô được tính theo công thức: $$P_{ki}% = \frac{P_{ki}}{P_{ti}} \times 100$$
  4. Định lượng các-bon: Áp dụng hệ số chuyển đổi các-bon chuẩn hóa quốc tế $C_{factor} = 0,485$ (tương đương $48,5%$ sinh khối khô kiệt) cho tất cả các phân đoạn sinh khối.
                              ┌────────────────────────────────────────┐
                              │     KHẢO SÁT HIỆN TRƯỜNG LÂM PHẦN      │
                              │ 54 OTC (Mai Sơn & Yên Châu, Sơn La)    │
                              └──────────────────┬─────────────────────┘
                                                 │
                        ┌────────────────────────┴────────────────────────┐
                        ▼                                                 ▼
       ┌─────────────────────────────────┐               ┌─────────────────────────────────┐
       │     ĐO ĐẾM KHÔNG PHÁ HỦY        │               │   GIẢI TÍCH PHÁ HỦY MẪU (IBP)   │
       ├─────────────────────────────────┤               ├─────────────────────────────────┤
       │• Đo D1.3, Hvn, Hdc, Dt toàn bộ  │               │• 36 cây tiêu chuẩn chặt hạ      │
       │  cây trong 54 OTC               │               │• Cắt khúc giải tích đĩa thân    │
       │• Thu mẫu 5 ô thảm tươi/mục (1m2)│               │• Phân đoạn: Thân, Vỏ, Cành, Lá  │
       │• Phân định Cấp đất II và III    │               │• Cân trọng lượng tươi ngay (Wt) │
       └────────────────┬────────────────┘               └────────────────┬────────────────┘
                        │                                                 │
                        └────────────────────────┬────────────────────────┘
                                                 │
                                                 ▼
                              ┌────────────────────────────────────────┐
                              │     PHÂN TÍCH PHÒNG THÍ NGHIỆM         │
                              │ Sấy 105°C đến khối lượng không đổi     │
                              │ Xác định Pki% -> Sinh khối khô Wk      │
                              │ Quy đổi Các-bon (C = Wk x 0,485)       │
                              └──────────────────┬─────────────────────┘
                                                 │
                                                 ▼
                              ┌────────────────────────────────────────┐
                              │   MÔ HÌNH HÓA TOÁN HỌC & TỐI ƯU HÓA    │
                              │ SPSS / Nonlinear Regression / R2       │
                              │ Weibull -> Schumacher -> Biểu lâm phần │
                              └────────────────────────────────────────┘

Data và phân tích

Dữ liệu thực nghiệm được xử lý bằng phần mềm chuyên dụng SPSS, Statistica và các module hồi quy phi tuyến trên nền tảng Microsoft Excel Solver:

  • Kiểm định phân bố thực nghiệm: Sử dụng phân bố Weibull 3 tham số để nắn dòng phân bố số cây theo đường kính $f(d)$ và chiều cao $f(h)$: $$f(x) = \frac{c}{b} \left( \frac{x-a}{b} \right)^{c-1} \exp\left[ -\left( \frac{x-a}{b} \right)^c \right]$$ với kiểm định mức độ phù hợp $\chi^2$ (Chi-square test) và tiêu chuẩn Kolmogorov-Smirnov ở mức ý nghĩa $p < 0,05$.
  • Mô hình hóa sinh trưởng: Khảo sát và so sánh 4 hàm phi tuyến kinh điển (Gompertz, Richards, Korf và Schumacher). Hàm Schumacher dạng tuyến tính hóa: $$\ln Y = a - \frac{b}{t}$$ cho hệ số xác định ($R^2$) cao nhất, vượt trội trong việc mô tả tiến trình biến đổi của $D_{1.3}$, $H_{vn}$ và $V$ theo tuổi $t$.
  • Đánh giá tăng trưởng: Tính toán chính xác lượng tăng trưởng thường xuyên hàng năm ($Z_d, Z_h, Z_v$) và tăng trưởng bình quân chung ($\Delta_d, \Delta_h, \Delta_v$) để xác định chính xác thời điểm giao nhau của hai đường cong - chỉ thị tuổi thành thục số lượng và thời điểm can thiệp tỉa thưa kỹ thuật.

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

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

Kết quả phân tích từ 54 lâm phần thực nghiệm và 36 cây tiêu chuẩn giải tích tại Sơn La mang lại 5 hệ thống phát hiện định lượng mang tính đột phá:

    ┌──────────────────────────────────────────────────────────────────────────────────┐
    │              TỔNG HỢP CÁC KẾT QUẢ ĐỊNH LƯỢNG ĐỘT PHÁ CỦA LUẬN ÁN                 │
    └────────────────────────────────────────┬─────────────────────────────────────────┘
                                             │
      ┌──────────────────────────────────────┼──────────────────────────────────────┐
      ▼                                      ▼                                      ▼
┌─────────────────────────┐    ┌─────────────────────────┐    ┌─────────────────────────┐
│ CẤU TRÚC PHÂN BỐ & MÔ   │    │  ĐỘNG THÁI SINH TRƯỞNG  │    │ SINH KHỐI & DỰ TRỮ C    │
│ HÌNH TƯƠNG QUAN LÂM PHẦN│    │   MÔ HÌNH SCHUMACHER    │    │ CẤP ĐẤT II vs CẤP III   │
├─────────────────────────┤    ├─────────────────────────┤    ├─────────────────────────┤
│• N/D1.3: 100% phù hợp   │    │• D1.3 tăng mạnh tuổi 1-6│    │• Sinh khối khô cây gỗ:  │
│  hàm Weibull (lệch trái │    │• Chiều cao giao điểm:   │    │  Cấp II: 104,18 tấn/ha  │
│  chiếm ưu thế).         │    │  Tuổi 5 (Cấp II: Zh=1,8m│    │  Cấp III: 77,55 tấn/ha  │
│• Tương quan H-D: Mô hình│    │  Cấp III: Zh=1,5m/năm). │• Tích lũy Các-bon rừng: │
│  Logarit đạt R2=0,78-   │    │• Thể tích giao điểm:    │    │  Cấp II: 54,23 tấn C/ha │
│  0,89.                  │    │  Tuổi 11-12 (Cấp II: Zv=│    │  Cấp III: 40,84 tấn C/ha│
│• Tương quan Dt-D1.3: Dt │    │  12,5 m3/ha/năm).       │• Thân chiếm 68,5-71,2%  │
│  = 0,215 + 0,186.D1.3   │    │• Cấp II vượt Cấp III từ │  tổng sinh khối cây.    │
│  (R2 > 0,82).           │    │  28% đến 35% thể tích.  │• Thảm mục/tươi: 5-8%.   │
└─────────────────────────┘    └─────────────────────────┘    └─────────────────────────┘
  1. Quy luật phân bố cấu trúc lâm phần: Phân bố số cây theo đường kính ngang ngực ($N/D_{1.3}$) và chiều cao vút ngọn ($N/H_{vn}$) của 54 lâm phần Tếch tuổi 13 đều tương thích cao với hàm phân bố Weibull ($p > 0,05$). Đa số các lâm phần có dạng phân bố lệch trái (hệ số bất đối xứng $c > 0$), chứng tỏ mật độ cây cá thể tập trung chủ yếu ở các cỡ kính nhỏ và trung bình, phản ánh tình trạng cạnh tranh không gian gay gắt do rừng chưa được tỉa thưa kịp thời.
  2. Xác lập tương quan sinh trắc học tin cậy:
    • Tương quan giữa chiều cao vút ngọn ($H_{vn}$) và đường kính ($D_{1.3}$) tuân theo quy luật hàm Logarit tự nhiên: $H_{vn} = a + b \cdot \ln(D_{1.3})$ với hệ số xác định $R^2$ dao động từ $0,78$ đến $0,89$.
    • Tương quan giữa đường kính tán ($D_t$) và đường kính thân ($D_{1.3}$) thể hiện mối quan hệ tuyến tính chặt chẽ: $D_t = 0,215 + 0,186 \cdot D_{1.3}$ ($R^2 > 0,82$), cung cấp công thức toán học chuẩn xác để tính toán diện tích chiếm đất của tán cây phục vụ xác định mật độ nuôi dưỡng.
  3. Quy luật động thái sinh trưởng theo cấp đất:
    • Tiến trình sinh trưởng đường kính và chiều cao diễn ra rất nhanh trong 6 năm đầu, sau đó giảm tốc rõ rệt từ tuổi thứ 8.
    • Điểm cắt giữa tăng trưởng thường xuyên ($Z_h$) và tăng trưởng bình quân ($\Delta_h$) của chiều cao xuất hiện sớm ở giai đoạn 4 - 5 tuổi (đạt đỉnh $1,8\text{ m/năm}$ trên Cấp đất II và $1,5\text{ m/năm}$ trên Cấp đất III).
    • Điểm cắt tăng trưởng thể tích ($Z_v = \Delta_v$) xảy ra ở tuổi 11 - 12 đối với Cấp đất II ($Z_v\text{ cực đại} = 12,5\text{ m}^3/\text{ha/năm}$) và tuổi 12 - 13 đối với Cấp đất III ($Z_v\text{ cực đại} = 9,1\text{ m}^3/\text{ha/năm}$).
  4. Cấu trúc phân đoạn sinh khối cây cá lẻ:
    • Phân tích 36 cây mẫu cho thấy tỷ lệ chất khô kiệt bình quân đạt $49,5 - 52,3%$ trọng lượng tươi.
    • Cấu trúc sinh khối khô của cây gỗ Tếch tập trung chủ yếu ở thân ($68,5 - 71,2%$), tiếp theo là cành ($15,4 - 17,8%$), rễ/vỏ ($8,2 - 10,1%$) và thấp nhất là lá ($4,1 - 5,6%$). Phương trình tương quan sinh khối khô toàn cây dạng allometric $W_k = a \cdot (D_{1.3}^2 H_{vn})^b$ đạt độ chính xác $R^2 > 0,94$.
  5. Định lượng trữ lượng tích lũy các-bon lâm phần:
    • Trên Cấp đất II: Tổng sinh khối khô đạt $108,42\text{ tấn/ha}$ (trong đó tầng cây gỗ chiếm $104,18\text{ tấn/ha}$, thảm tươi và thảm mục chiếm $4,24\text{ tấn/ha}$), tích lũy tương ứng $54,23\text{ tấn C/ha}$ (tương đương hấp thụ $198,82\text{ tấn CO}_2\text{ tương đương/ha}$).
    • Trên Cấp đất III: Tổng sinh khối khô đạt $81,62\text{ tấn/ha}$ (tầng cây gỗ chiếm $77,55\text{ tấn/ha}$, thảm tươi và thảm mục chiếm $4,07\text{ tấn/ha}$), tích lũy tương ứng $40,84\text{ tấn C/ha}$ (tương đương $149,74\text{ tấn CO}_2\text{ tương đương/ha}$).

Implications đa chiều

  • Về mặt học thuật: Cung cấp bộ cơ sở dữ liệu thực nghiệm chuẩn hóa đầu tiên về sinh trắc học và sinh khối rừng trồng Tếch vùng núi cao phía Bắc, làm cầu nối so sánh với các hệ sinh thái Tếch tại Tây Nguyên, Đông Nam Á và Trung - Nam Mỹ.
  • Về mặt kỹ thuật lâm sinh: Chứng minh rằng rừng Tếch 13 tuổi tại Sơn La đang rơi vào trạng thái suy giảm sinh trưởng nghiêm trọng do mật độ ban đầu quá dày ($1.660 - 2.500\text{ cây/ha}$) mà không được tỉa thưa. Luận án đề xuất quy trình tỉa thưa 2 lần: Lần 1 ở tuổi 5 - 6 (cường độ chặt $40 - 50%$ số cây, đưa mật độ về $1.000 - 1.100\text{ cây/ha}$), Lần 2 ở tuổi 10 - 11 (cường độ chặt $30 - 35%$, giữ lại $600 - 700\text{ cây/ha}$ nuôi dưỡng đến chu kỳ khai thác chính 25 - 30 năm).
  • Về mặt chính sách và kinh tế môi trường: Tạo luận cứ khoa học để tỉnh Sơn La và các địa phương Tây Bắc xây dựng đề án tín chỉ các-bon rừng trồng, thu hút nguồn tài chính chi trả dịch vụ môi trường rừng (PFES), gia tăng giá trị kinh tế trên mỗi hecta rừng Tếch từ $20 - 30%$ ngoài doanh thu gỗ tròn.

Limitations và Future Research

Mặc dù đạt được những kết luận mang tính hệ thống cao, luận án ghi nhận một số giới hạn thực tiễn và phương pháp luận:

  1. Giới hạn chuỗi thời gian điều tra: Nghiên cứu thực hiện trên các lâm phần Tếch đồng loạt ở tuổi 13 (tuổi cao nhất hiện có tại thời điểm nghiên cứu ở Sơn La). Việc ngoại suy tiến trình sinh trưởng cho giai đoạn thành thục công nghệ ($25 - 35\text{ tuổi}$) dựa trên mô hình Schumacher có thể tiềm ẩn sai số nhất định do chưa có số liệu quan trắc dài hạn thực tế của rừng già.
  2. Phạm vi lập địa: Mẫu nghiên cứu mới chỉ phân lập được 2 cấp đất (Cấp II và Cấp III), chưa xuất hiện Cấp đất I (rất tốt) do đặc thù điều kiện địa hình dốc và thổ nhưỡng phân tán tại vùng khảo sát Mai Sơn, Yên Châu.
  3. Cơ chế đo đếm bể các-bon đất: Luận án tập trung định lượng sinh khối trên mặt đất (cây gỗ, cây bụi, thảm tươi, thảm mục), chưa tiến hành khoan lấy mẫu đất sâu ($0 - 100\text{ cm}$) để phân tích tổng hàm lượng các-bon hữu cơ trong đất (SOC - Soil Organic Carbon).

Định hướng nghiên cứu tiếp theo (Future Research Agenda):

  • Thiết lập hệ thống định vị ô tiêu chuẩn định vị vĩnh viễn (Permanent Sampling Plots - PSP) để theo dõi động thái tăng trưởng liên tục sau khi thực hiện các công thức tỉa thưa thực nghiệm.
  • Nghiên cứu bổ sung giải tích sinh khối hệ rễ cọc và rễ chùm dưới mặt đất thông qua phương pháp đào hố chuyên sâu nhằm hoàn thiện phương trình allometric toàn phần ($Biomass_{total}$).
  • Nghiên cứu tương tác giữa đặc điểm cấu tạo thớ gỗ, tỷ trọng gỗ Tếch theo các cấp đất và chế độ tỉa thưa nhằm đánh giá chất lượng cơ lý của gỗ phục vụ công nghiệp chế biến xuất khẩu.

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

  • Tác động học thuật: Công trình là tài liệu tham khảo chuẩn mực trong giảng dạy đại học và sau đại học chuyên ngành Lâm sinh, Quản lý tài nguyên rừng tại các trường đại học chuyên ngành (Đại học Lâm nghiệp, Đại học Nông Lâm Thái Nguyên). Cung cấp hệ thống phương trình hồi quy sinh thái phục vụ các mô hình dự báo sinh quyển quốc gia.
  • Chuyển đổi kỹ thuật ngành lâm nghiệp địa phương: Định hình lại tập quán trồng rừng Tếch tự phát tại Sơn La và Tây Bắc. Thay thế phương thức trồng thuần dày đặc không chăm sóc bằng quy trình kỹ thuật lâm sinh chuẩn hóa từ khâu chọn đất (ưa đất Feralit phát triển trên đá vôi, bazan, tầng dày $> 80\text{ cm}$, thoát nước tốt) đến kỹ thuật tỉa cành, tỉa thưa định kỳ.
  • Giá trị kinh tế - xã hội: Góp phần chuyển dịch cơ cấu cây trồng lâm nghiệp giá trị cao, nâng cao thu nhập cho đồng bào các dân tộc thiểu số (Thái, Kinh, Mông, Mường) tại Sơn La thông qua mô hình Nông lâm kết hợp (Taungya) trồng xen hoa màu, đậu đỗ trong 3 năm đầu khi rừng Tếch chưa khép tán.

Đố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 mẫu mực về kết hợp giữa giải tích điều tra lâm học truyền thống với kỹ thuật lượng hóa sinh khối hiện đại theo chuẩn quốc tế (IBP, IPCC).
  • Chủ rừng, công ty lâm nghiệp và ban quản lý rừng phòng hộ: Nắm vững biểu tra mật độ, biểu cấp đất và lịch thời gian tỉa thưa chuẩn xác nhằm tối ưu hóa đường kính gỗ thương phẩm xuất khẩu (đạt quy cách đường kính bạnh gốc $> 40\text{ cm}$).
  • Nhà hoạch định chính sách lâm nghiệp cấp tỉnh và Trung ương: Có số liệu thực chứng để xây dựng hạn ngạch phát thải, định giá dịch vụ hấp thụ các-bon rừng trong quy hoạch phát triển lâm nghiệp bền vững vùng Tây Bắc tầm nhìn 2030 - 2050.

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ì?

Luận án đã chứng minh sự phù hợp tuyệt đối của mô hình sinh trưởng Schumacher trong việc mô phỏng tiến trình tích lũy thể tích ($V$) và đường kính ($D_{1.3}$) của loài Tectona grandis L. tại vùng sinh thái núi cao phía Bắc Việt Nam, mở rộng lý thuyết sản lượng sinh thái của Eichhorn và Whittaker sang một loài cây gỗ du nhập có chu kỳ rụng lá mùa đông.

2. Điểm mới về phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu của Bảo Huy (1995) tại Đắk Lắk và Ly Meng Seang (2009) tại Campuchia (chỉ dừng lại ở cấu trúc và biểu thể tích gỗ tròn), luận án đã thực hiện phương pháp tiếp cận toàn diện theo tiêu chuẩn IBP: giải tích đồng thời 36 cây mẫu phá hủy toàn phần, sấy khô kiệt $105^\circ\text{C}$, phân đoạn chi tiết 4 bộ phận cây gỗ kết hợp với 270 ô tiêu chuẩn thảm tươi/thảm mục để xây dựng hệ thống phương trình allometric sinh khối khô và bản đồ các-bon theo cấp đất.

3. Phát hiện bất ngờ nhất từ số liệu thực nghiệm là gì?

Tốc độ tăng trưởng chiều cao của cây Tếch tại Sơn La đạt đỉnh rất sớm ở tuổi 4 - 5 ($Z_h = 1,8\text{ m/năm}$ trên Cấp đất II), tương đương với vùng phân bố tự nhiên tại Ấn Độ và Myanmar. Tuy nhiên, nếu sau tuổi thứ 6 không thực hiện tỉa thưa lần 1, rừng sẽ rơi vào trạng thái "đình trệ sinh trưởng" do tán lá che phủ cạnh tranh cực độ, làm giảm $35 - 40%$ tiềm năng tăng trưởng đường kính ở giai đoạn sau.

4. Giao thức nhân bản (Replication Protocol) có được cung cấp rõ ràng không?

Luận án cung cấp hệ thống giao thức chi tiết từng bước: từ kỹ thuật chọn ô tiêu chuẩn $500\text{ m}^2$, phương pháp chọn cây tiêu chuẩn theo công thức phân cấp Draudt, quy trình cưa khúc lấy đĩa gỗ giải tích $1\text{ m}$, kỹ thuật cân đo tươi - khô, đến thuật toán lập trình hàm Weibull và Schumacher trên phần mềm thống kê, đảm bảo khả năng lặp lại hoàn toàn tại các vùng sinh thái tương đồng.

5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Trọng tâm là xây dựng mô hình quản trị số hóa rừng Tếch: tiếp tục giải tích thân cây ở tuổi 20, 25 và 30; tích hợp dữ liệu viễn thám LiDAR và ảnh vệ tinh đa phổ để dự báo trữ lượng sinh khối tán lá tự động; và thiết lập mô hình kinh doanh tín chỉ các-bon tự nguyện (VCS) kết hợp thương mại hóa gỗ Tếch có chứng chỉ bền vững FSC.

Kết luận

  1. Công trình luận án tiến sĩ của Nguyễn Công Hoan (2014) là nghiên cứu toàn diện, tiên phong và có độ tin cậy khoa học cao nhất về sinh thái học lâm vật và cấu trúc rừng trồng Tếch (Tectona grandis L.) tại Sơn La và miền Bắc Việt Nam.
  2. Xác lập vững chắc quy luật phân bố số cây theo đường kính và chiều cao tuân theo phân bố Weibull 3 tham số; quan hệ $H_{vn}/D_{1.3}$ tuân theo hàm Logarit tự nhiên; quan hệ $D_t/D_{1.3}$ mang tính chất tuyến tính thuận.
  3. Chứng minh mô hình Schumacher là công cụ toán học tối ưu để dự báo sinh trưởng đường kính, chiều cao và thể tích lâm phần Tếch theo cấp đất.
  4. Lượng hóa chính xác trữ lượng sinh khối khô và tích lũy các-bon: Cấp đất II đạt $108,42\text{ tấn sinh khối/ha}$ ($54,23\text{ tấn C/ha}$), Cấp đất III đạt $81,62\text{ tấn sinh khối/ha}$ ($40,84\text{ tấn C/ha}$); trong đó thân cây đóng vai trò bể chứa trung tâm ($68,5 - 71,2%$).
  5. Xây dựng hệ thống giải pháp kỹ thuật lâm sinh hoàn chỉnh: thiết lập biểu điều tiết mật độ, quy định cường độ tỉa thưa 2 giai đoạn (tuổi 5 - 6 và tuổi 10 - 11), làm cơ sở khoa học cốt lõi cho kinh doanh rừng trồng Tếch bền vững và khai thác dịch vụ môi trường rừng đa mục tiêu.