Tổng quan về luận án

Chiến lược tái cơ cấu ngành lâm nghiệp Việt Nam theo Quyết định số 774/QĐ-BNN-TCLN của Bộ Nông nghiệp và Phát triển nông thôn đặt ra mục tiêu quy hoạch 650.000 ha rừng trồng cung cấp gỗ lớn trên toàn quốc, trong đó vùng Đông Bắc Bộ chiếm 134.000 ha. Mặc dù tổng diện tích rừng trồng sản xuất cả nước tính đến ngày 31 tháng 12 năm 2021 đạt 4.380.000 ha và vùng Đông Bắc Bộ sở hữu 1.012.000 ha rừng trồng với độ che phủ 56,34% (Bộ Nông nghiệp và PTNT, 2022), phần lớn diện tích vẫn duy trì phương thức quảng canh hoặc thâm canh chu kỳ ngắn (5–7 năm) để sản xuất dăm mảnh và bột giấy. Thực trạng này dẫn đến khủng hoảng thiếu hụt nguyên liệu gỗ xẻ chất lượng cao, buộc Việt Nam phải nhập khẩu trên 4,71 triệu m³ gỗ tròn và gỗ xẻ hàng năm (Tổng cục Hải quan, 2021) khi nhu cầu gỗ tròn chế biến trong nước đã tăng vọt từ 31,30 triệu m³ lên 41,72 triệu m³ (Nguyễn Tôn Quyền, 2019).

Khoảng trống nghiên cứu (research gap) trọng yếu nằm ở việc chuyển dịch cơ cấu cây trồng từ Keo tai tượng (Acacia mangium) và Keo lai (Acacia hybrid) sang Keo lá tràm (Acacia auriculiformis A. Cunn. ex Benth) phục vụ kinh doanh gỗ lớn tại vùng Đông Bắc Bộ. Trước đây, Keo lá tràm chỉ được khảo nghiệm và phát triển chủ yếu tại Bắc Trung Bộ, Nam Trung Bộ và Đông Nam Bộ (Lê Đình Khả, 2003; Hà Huy Thịnh, 2006). Vùng Đông Bắc Bộ với đặc trưng mùa đông lạnh cục bộ (nhiệt độ hạ thấp dưới 8–10°C) và đất đai đồi núi dốc thoái hóa nghiêm trọng vẫn chưa có hệ thống khảo nghiệm mở rộng giống tiến bộ kỹ thuật cũng như gói giải pháp kỹ thuật lâm sinh thâm canh liên hoàn.

Luận án tiến sĩ Lâm nghiệp của nghiên cứu sinh Phạm Đình Sâm với đề tài "Nghiên cứu một số biện pháp kỹ thuật trồng rừng thâm canh Keo lá tràm cung cấp gỗ lớn ở vùng Đông Bắc Bộ" (ngành Lâm sinh, mã số: 9.62.02.05, người hướng dẫn khoa học: PGS. TS. Nguyễn Huy Sơn, bảo vệ tại Viện Khoa học Lâm nghiệp Việt Nam năm 2022) đã giải quyết toàn diện khoảng trống khoa học này. Luận án đặt ra và kiểm định 5 giả thuyết nghiên cứu chính:

  1. Giả thuyết H1: Tồn tại sự biến dị di truyền có ý nghĩa về sinh trưởng, độ thẳng thân và khả năng chống chịu điều kiện mùa đông lạnh giữa các dòng vô tính Keo lá tràm khi khảo nghiệm mở rộng tại vùng Đông Bắc Bộ.
  2. Giả thuyết H2: Biện pháp quản lý vật liệu hữu cơ sau khai thác (giữ lại toàn bộ thực bì không đốt) cải thiện vượt trội dinh dưỡng tầng đất mặt so với phương thức dọn sạch hoặc phát gom truyền thống.
  3. Giả thuyết H3: Cơ giới hóa làm đất (cày ngầm) kết hợp cuốc hố theo đường đồng mức phá vỡ tầng đế cày, tạo điều kiện mở rộng không gian dinh dưỡng bộ rễ, tăng tỷ lệ sống và năng suất gỗ so với cuốc hố thủ công đơn thuần.
  4. Giả thuyết H4: Bổ sung cân đối phân bón vô cơ (NPK) và hữu cơ vi sinh trong giai đoạn bón lót, kết hợp bón thúc P₂O₅ và K₂O dựa trên chẩn đoán dinh dưỡng khoáng trong lá đáp ứng kịp thời nhu cầu sinh lý của Keo lá tràm trên đất Feralit nghèo lân.
  5. Giả thuyết H5: Mật độ trồng ban đầu 1.110–1.660 cây/ha tạo lập cấu trúc lâm phần tối ưu cho phép xúc tiến phân hóa đường kính và hạn chế tỉa cành nhân tạo phức tạp trước giai đoạn tỉa thưa nuôi dưỡng năm thứ 5.

Khung lý thuyết của công trình tích hợp Lý thuyết Thâm canh Lâm sinh (Nguyễn Xuân Quát, 1995, 1998; Vũ Đình Hèo, 1986), Lý thuyết Quản lý Lập địa Bền vững (Gonçalves et al., 1997; Nambiar & Harwood, 2014) và Định luật Tối thiểu Liebig về dinh dưỡng khoáng cây rừng (Attiwill & Adams, 1993; Vitousek, 1982). Đóng góp mang tính đột phá của luận án được định lượng qua việc xác định thành công 6 dòng vô tính Keo lá tràm ưu việt (Clt98, Clt26, Clt57, Clt7, Clt133, AA9) đạt lượng tăng trưởng bình quân thể tích thân cây ($\Delta M$) từ 15,3 đến 23,0 m³/ha/năm ở giai đoạn 5 năm tuổi, vượt chuẩn quy định quốc gia về gỗ lớn ($\ge 15\text{ m}^3/\text{ha}/\text{năm}$).

Phạm vi nghiên cứu bao gồm điều tra hiện trạng 18 ô tiêu chuẩn (diện tích 500 m²/ô, tổng quy mô 16,7 ha) tại 4 tỉnh/thành phố đại diện (Hải Phòng, Bắc Giang, Hà Nội, Tuyên Quang) kết hợp hệ thống 5 thí nghiệm lâm sinh thâm canh chính quy trên diện tích 3,1104 ha tại Trạm thực nghiệm Uông Bí, Quảng Ninh trong giai đoạn 2016–2022.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về Keo lá tràm phản ánh sự tiến triển từ việc thử nghiệm loài bản địa/di thực sang cải thiện di truyền chuyên sâu và quản lý lâm sinh bền vững. Về mặt phân loại và sinh thái, Keo lá tràm (Acacia auriculiformis) thuộc chi Acacia, họ Đậu (Leguminosae), phân bố tự nhiên tại Bắc Queensland, Lãnh thổ Bắc Úc, Papua New Guinea và Indonesia (Lewis, 2000; Boland, 1990; Vélez-Gavilán, 2016). Các nghiên cứu quốc tế của Doran et al. (1997) khẳng định ngưỡng sinh thái giới hạn của loài là độ cao >600 m, lượng mưa <1.000 mm/năm và khả năng chống chịu gió bão trung bình.

Dòng nghiên cứu chọn giống quốc tế ghi nhận đóng góp của Shukor et al. (1993, 1997) tại Malaysia khi khảo nghiệm 28 xuất xứ, chứng minh xuất xứ từ Bắc Queensland (Wenlock River, Kings Plain) vượt trội về sinh trưởng và khối lượng riêng của gỗ so với Northern Territory hay Papua New Guinea. Tại Thái Lan, Luangviriyasaeng & Pinyopusarerk (2002) đánh giá biến dị di truyền của 30 gia đình Keo lá tràm thế hệ 2, ghi nhận hệ số di truyền cấp cá thể về chiều cao đạt 0,33, đường kính đạt 0,40, độ thẳng thân đạt 0,45 và khối lượng riêng đạt 0,28. Nghiên cứu của Hendrati & Nurrohmah (2018) tại Indonesia chỉ rõ tỷ lệ gỗ lõi của Keo lá tràm có sự phân hóa di truyền mạnh mẽ bắt đầu từ 30 tháng tuổi và mở rộng theo cấp tuổi.

Tại Việt Nam, các công trình tiên phong của Lê Đình Khả (1993, 2003), Nguyễn Hoàng Nghĩa (1996, 2003), Hà Huy Thịnh (2006, 2015) và Phí Hồng Hải (2008) đã tuyển chọn và công nhận 3 xuất xứ cùng 34 dòng vô tính tiến bộ kỹ thuật (như AA1, AA9, Clt7, Clt26, Clt57, Clt98, Clt133, BVlt83, BVlt85) với năng suất đạt 13–35,7 m³/ha/năm, nhưng 90% các khảo nghiệm này định vị tại miền Trung và miền Nam.

                    ┌────────────────────────────────────────────────────────┐
                    │      CÁC LUỒNG TRANH LUẬN HỌC THUẬT LÂM SINH          │
                    └──────────────────────────┬─────────────────────────────┘
                                               │
               ┌───────────────────────────────┴───────────────────────────────┐
               ▼                                                               ▼
┌──────────────────────────────┐                              ┌──────────────────────────────┐
│  Tranh luận Xử lý Thực bì   │                              │   Tranh luận Chiến lược      │
│     & Quản lý Dinh dưỡng     │                              │      Mật độ & Luân kỳ        │
├──────────────────────────────┤                              ├──────────────────────────────┤
│ QUAN ĐIỂM TRUYỀN THỐNG:      │                              │ QUAN ĐIỂM GỖ NHỎ QUẢNG CANH: │
│ Đốt dọn thực bì, cày lật     │                              │ Trồng dày (2.500-3.300 c/ha),│
│ cơ giới toàn diện, kiểm soát │                              │ khai thác sớm 5-7 năm làm    │
│ cỏ triệt để (FL).            │                              │ dăm, bột giấy.               │
│                              │                              │                              │
│ QUAN ĐIỂM QUẢN LÝ BỀN VỮNG:  │                              │ QUAN ĐIỂM GỖ LỚN THÂM CANH:  │
│ Bảo tồn vật rơi rụng & phụ   │                              │ Mật độ vừa phải (1.110-1.660 │
│ phẩm sau khai thác (VLHCSKT),│                              │ c/ha), điều tiết không gian  │
│ duy trì chu trình N, P, K    │                              │ dinh dưỡng, tỉa thưa nuôi    │
│ (Gonçalves et al., 1997;     │                              │ dưỡng lấy gỗ xẻ đường kính   │
│ Nambiar & Harwood, 2014;     │                              │ D1.3 >= 20 cm (Beadle, 2006; │
│ Kiều Tuấn Đạt, 2014).        │                              │ Trần Lâm Đồng, 2019).        │
└──────────────────────────────┘                              └──────────────────────────────┘

Trong y văn lâm sinh học tồn tại hai cuộc tranh luận học thuật sâu sắc:

  1. Xử lý tàn dư thực bì: Quan điểm truyền thống ủng hộ dọn sạch và đốt thực bì để giải phóng mặt bằng cơ giới và diệt trừ mầm bệnh. Ngược lại, trường phái quản lý lập địa của Nambiar & Harwood (2014) và Binkley et al. (1997) chứng minh lớp thảm mục rừng keo tích lũy tới 96,8 kg N/ha; 5,7 kg P/ha; 40,0 kg K/ha; 30,6 kg Ca/ha. Việc đốt tàn dư làm thất thoát dinh dưỡng nghiêm trọng và gây suy thoái đất luân kỳ sau (Vũ Đình Hưởng et al., 2020).
  2. Mật độ và cấu trúc rừng kinh doanh gỗ lớn: Bùi Việt Hải (1998) đề xuất trồng mật độ dày 1.660–3.300 cây/ha để kích thích tăng trưởng chiều cao và tỉa cành tự nhiên. Trong khi đó, Beadle et al. (2013) và Trần Lâm Đồng (2019) bảo vệ quan điểm mật độ ban đầu chỉ nên từ 1.110 đến 1.330 cây/ha kết hợp tỉa thưa sớm để tối đa hóa đường kính ngang ngực ($D_{1.3}$), bởi Keo lá tràm trồng quá dày sẽ bị suy giảm sinh trưởng đường kính và tự tỉa thưa tự nhiên gây tổn thất sinh khối.

So sánh với các nghiên cứu quốc tế tại Sumatra - Indonesia (Mendham et al., 2017) và Sabah - Malaysia (MacDicken, 1988), công trình của Phạm Đình Sâm đã định vị rõ nét sự khác biệt: Trong khi các nghiên cứu quốc tế thực hiện trên đất nhiệt đới xích đạo có điều kiện nhiệt ẩm lý tưởng quanh năm, nghiên cứu này giải quyết bài toán lâm sinh tại vùng cận nhiệt đới gió mùa Đông Bắc Bộ với biên độ nhiệt mùa đông thấp và địa hình đồi núi dốc bị rửa trôi mạnh.

Đó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 Lý thuyết Thâm canh Rừng trồng của Nguyễn Xuân Quát (1995, 1998) từ phạm vi khái niệm chung sang mô hình định lượng thực chứng cụ thể cho loài cây họ Đậu cải thiện di truyền tại vùng Đông Bắc Bộ. Lý thuyết được hoàn thiện qua việc bổ sung cơ chế liên hoàn: Chọn giống thích nghi $\rightarrow$ Xử lý thực bì bảo tồn hữu cơ $\rightarrow$ Làm đất cơ giới theo đường đồng mức $\rightarrow$ Chẩn đoán khoáng lá bón phân cân đối $\rightarrow$ Điều tiết mật độ tỉa thưa cấu trúc.

                    ┌────────────────────────────────────────────────────────┐
                    │      KHUNG PHÂN TÍCH LÂM SINH THÂM CANH LIÊN HOÀN      │
                    └──────────────────────────┬─────────────────────────────┘
                                               │
         ┌─────────────────────────────────────┼─────────────────────────────────────┐
         ▼                                     ▼                                     ▼
┌──────────────────┐                 ┌──────────────────┐                  ┌──────────────────┐
│  KIỂU GEN (G)    │                 │   LẬP ĐỊA (E)    │                  │  LÂM SINH (S)    │
├──────────────────┤                 ├──────────────────┤                  ├──────────────────┤
│ 6 Dòng vô tính   │                 │ Đất Feralit đỏ   │                  │ Thực bì: Rải đều │
│ Clt98, Clt26,    │    TƯƠNG TÁC    │ vàng trên đá biến│    TÁC ĐỘNG      │ Làm đất: Cày ngầm│
│ Clt57, Clt7,     │ ──────────────> │ chất (Fs).       │ <─────────────── │ Dinh dưỡng: Lá   │
│ Clt133, AA9      │   G x E x S     │ Nhiệt độ mùa     │  ĐIỀU TIẾT SINH  │ N:P:K + Bio-NPK  │
│ Năng suất 15,3 - │                 │ đông < 10°C,     │      TRƯỞNG      │ Mật độ: 1.110 -  │
│ 23,0 m³/ha/năm   │                 │ dốc 15 - 20°     │                  │ 1.660 cây/ha     │
└──────────────────┘                 └──────────────────┘                  └──────────────────┘

Nghiên cứu làm sâu sắc thêm Lý thuyết Quản lý Lập địa (Gonçalves et al., 1997) bằng việc lượng hóa chu trình hồi hoàn vật chất hữu cơ tại hiện trường đồi dốc. Kết quả chứng minh việc giữ lại toàn bộ vật liệu hữu cơ sau khai thác không chỉ bảo toàn 84,8 kg N/ha, 40,9 kg P/ha, 57,8 kg K/ha bị mất đi do đốt dọn (dẫn chứng từ Vũ Đình Hưởng et al., 2020), mà còn hoạt hóa vi sinh vật đất, tạo đệm giữ ẩm cho tầng đất mặt $0-20\text{ cm}$.

Mô hình lý thuyết phát triển qua hệ thống các mệnh đề khoa học:

  • Mệnh đề 1 (Tương tác G $\times$ E): Khả năng biểu hiện năng suất của các dòng vô tính Keo lá tràm nguồn gốc Nam Trung Bộ/Đông Nam Bộ khi đưa ra Đông Bắc Bộ phụ thuộc chặt chẽ vào mức độ chống chịu nhiệt độ thấp và độ sâu tầng đất hữu hiệu ($>80\text{ cm}$).
  • Mệnh đề 2 (Động lực học Dinh dưỡng): Sự gia tăng sinh trưởng của Keo lá tràm bị chi phối bởi hàm lượng phospho dễ tiêu ($P_2O_5$) hơn là đạm ($N$), do vi khuẩn nốt sần rễ có khả năng cố định nitơ khí quyển nhưng đòi hỏi lượng ATP cao từ việc hấp thụ lân trong đất chua nhiệt đới.

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

Khung phân tích của luận án tích hợp liên ngành 3 trụ cột lý thuyết: Di truyền học số lượng cây rừng, Thổ nhưỡng học - Sinh lý dinh dưỡng khoáng, và Động thái học lâm phần (Stand Dynamics).

  1. Trụ cột Di truyền: Đánh giá phương sai kiểu gen, hệ số biến động đường kính ($S_D$), chiều cao ($S_H$), và chỉ số phẩm chất thân cây (độ thẳng $D_{tt}$, độ nhỏ cành $D_{nc}$, cấp sức khỏe $S_k$) theo Tiêu chuẩn Quốc gia TCVN 8755:2017.
  2. Trụ cột Dinh dưỡng lập địa: Phân tích tương quan giữa hàm lượng dinh dưỡng khoáng trong lá (N, P, K) theo tiêu chuẩn 10 TCN 452/453/454:2001 với độ phì nhiêu tầng đất theo TCVN 5979/8941/6498/8661/8662:2011, thiết lập cơ chế bón đón đầu nhu cầu sinh lý.
  3. Trụ cột Động thái học: Khảo sát quy luật phân bố số cây theo đường kính ($N/D$) và chiều cao ($N/H$) ở tuổi 5, xác định thời điểm giao cắt giữa lượng tăng trưởng thường xuyên ($Z_D$) và tăng trưởng bình quân ($\Delta D$) nhằm hoạch định chế độ tỉa thưa khoa học.

Điều kiện biên (Boundary conditions): Khung phân tích này áp dụng tối ưu cho lập địa đất Feralit đỏ vàng phát triển trên đá sét và biến chất (Fs), độ dày tầng đất $>80\text{ cm}$, độ dốc $<20^\circ$, lượng mưa bình quân 1.700–2.500 mm/năm và vùng khí hậu có mùa đông lạnh với nhiệt độ không kéo dài dưới 4–8°C.

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ủ thế giới quan thực chứng (positivism philosophy) và phương pháp luận thực nghiệm lâm sinh định lượng nghiêm ngặt. Thiết kế nghiên cứu được cấu trúc theo 2 cấp độ phân tích:

  • Cấp vĩ mô: Đánh giá thực trạng 18 ô tiêu chuẩn (OTC) tại Cát Bà (Hải Phòng - rừng 4, 12, 24 tuổi), Yên Thế (Bắc Giang - rừng 6 tuổi), Ba Vì (Hà Nội - rừng 12 tuổi) và Yên Sơn (Tuyên Quang - rừng 21 tuổi) với tổng diện tích 16,7 ha nhằm thiết lập đường chuẩn sinh thái và lập địa.
  • Cấp vi mô: Thiết lập hệ thống thí nghiệm nhân tố đầy đủ tại Uông Bí, Quảng Ninh trên diện tích 3,1104 ha:
    • Khảo nghiệm giống mở rộng: 12 dòng vô tính (Clt18, Clt19, Clt25, Clt26, Clt43, Clt57, Clt98, Clt133, BVlt83, BVlt85, AA9, Clt7), thiết kế khối ngẫu nhiên hoàn chỉnh (RCBD) theo hàng - cột, 4 lần lặp lại, 16 cây/ô (4 cây $\times$ 4 cây), mật độ 1.110 cây/ha.
    • Thí nghiệm biện pháp kỹ thuật lâm sinh: 5 thí nghiệm độc lập (xử lý thực bì, làm đất, mật độ, bón phân, chăm sóc) sử dụng 2 dòng chuẩn Clt7 và AA9 (trộn tỉ lệ 1:1), bố trí RCBD với 3 lần lặp lại, 64 cây/ô (8 hàng $\times$ 8 cây), diện tích ô 576 m², tổng cộng 18 công thức thí nghiệm.
┌────────────────────────────────────────────────────────────────────────────────────────┐
│              SƠ ĐỒ HỆ THỐNG PHƯƠNG PHÁP NGHIÊN CỨU VÀ BỐ TRÍ THÍ NGHIỆM                │
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
         ┌─────────────────────────────────┴─────────────────────────────────┐
         ▼                                                                   ▼
┌───────────────────────────────────┐               ┌───────────────────────────────────┐
│     ĐIỀU TRA THỰC TRẠNG LẬP ĐỊA   │               │     BỐ TRÍ THÍ NGHIỆM HIỆN TRƯỜNG │
│    & SINH TRƯỞNG RỪNG TRỒNG       │               │      (Uông Bí, Quảng Ninh - 3,11 ha)│
├───────────────────────────────────┤               ├───────────────────────────────────┤
│ • Quy mô: 18 OTC (500 m²/ô;       │               │ 1. Khảo nghiệm giống (RCBD, 4 lặp):│
│   tổng diện tích 16,7 ha)         │               │    12 dòng vô tính TBKT.          │
│ • Địa bàn: Cát Bà (Hải Phòng),    │               │ 2. Thí nghiệm Lâm sinh (RCBD, 3   │
│   Yên Thế (Bắc Giang), Ba Vì      │               │    lặp, 64 cây/ô, đo 36 cây lõi): │
│   (Hà Nội), Yên Sơn (Tuyên Quang).│               │    - Xử lý thực bì: 3 công thức   │
│ • Phân tích mẫu đất tầng 0-20 cm  │               │    - Kỹ thuật làm đất: 3 công thức│
│   và 30-50 cm (Dung trọng, pHKCl, │               │    - Mật độ trồng: 4 công thức    │
│   Nts, P2O5, K2O, CEC, Thành phần │               │    - Bón lót/bón thúc: 5 & 9 CTTN │
│   cơ giới theo TCVN).             │               │    - Chế độ chăm sóc: Các tần suất│
└───────────────────────────────────┘               └───────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích dữ liệu áp dụng triệt để các tiêu chuẩn kỹ thuật ngành và quốc gia:

  • Tiêu chuẩn đo đếm cây rừng: Đ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 cm; đo chiều cao vút ngọn ($H_{vn}$) và chiều cao dưới cành ($H_{dc}$) bằng sào đo cao chuyên dụng kết hợp máy đo cao laser Vertex IV; đo đường kính tán ($D_t$) theo 2 hướng Đông - Tây và Nam - Bắc bằng thước dây chính xác đến 0,1 m.
  • Đánh giá phẩm chất thân cây: Tuân thủ TCVN 8755:2017 với 3 chỉ tiêu phân cấp: Độ thẳng thân ($D_{tt}$: 1 - Thẳng, 2 - Cong nhẹ, 3 - Cong queo); Độ nhỏ cành ($D_{nc}$: 1 - Cành nhỏ $<1/3$ thân, 2 - Cành trung bình, 3 - Cành thô lớn); Sức khỏe ($S_k$: 1 - Tốt, 2 - Trung bình, 3 - Sâu bệnh/cụt ngọn).
  • Phương pháp tam giác đạc (Triangulation): Kết hợp kiểm chứng chéo giữa số liệu điều tra hồi cứu tại 18 OTC thực địa, phân tích hóa lý thổ nhưỡng trong phòng thí nghiệm và chuỗi số liệu đo đếm định kỳ hàng năm vào tháng 6 tại hiện trường Uông Bí.
  • Phân tích hóa nghiệm:
    • Đất: Dung trọng (TCVN 6861:2001), pHKCl (TCVN 5979:2007), Mùn tổng số (TCVN 8941:2011), Đạm tổng số (TCVN 6498:1999), $P_2O_5$ dễ tiêu theo phương pháp Oniani (TCVN 8661:2011), $K_2O$ dễ tiêu (TCVN 8662:2011), Dung tích hấp phụ CEC (TCVN 6646:2001).
    • Lá: Lấy mẫu lá non đầu cành, lá bánh tẻ giữa cành và lá già gốc cành từ 4 cây tâm ô của 5 công thức bón phân (15 mẫu tổ hợp), phân tích N (10 TCN 452:2001), P (10 TCN 453:2001), K (10 TCN 454:2001).

Data và phân tích

Toàn bộ dữ liệu đo đếm 36 cây vùng lõi trong mỗi ô thí nghiệm (loại bỏ hiệu ứng mép biên 28 cây) được xử lý thống kê sinh học lâm nghiệp bằng phần mềm SAS và SPSS.

  • Kiểm định phân tích phương sai (One-way và Two-way ANOVA) để đánh giá sự sai khác giữa các nghiệm thức.
  • Phân tách giá trị trung bình bằng kiểm định hậu nghiệm Bonferroni (Bonferroni post-hoc test) ở mức ý nghĩa $\alpha = 0,05$ và xác suất $F.sig < 0,001$.
  • Ước tính thể tích thân cây cá thể ($V$) và trữ lượng lâm phần ($M$) thông qua công thức hình số lâm nghiệp chuẩn:

$$V = \frac{\pi}{4} \cdot (D_{1.3})^2 \cdot H_{vn} \cdot f_{1.3}$$

Trong đó, hình số thân cây $f_{1.3}$ của Keo lá tràm được hiệu chỉnh chính xác theo cấp tuổi và lập địa điều tra thực tế. Mô hình hóa cấu trúc phân bố $N/D$ và $N/H$ được kiểm định qua hàm phân bố Weibull và kiểm định độ phù hợp $\chi^2$ (Chi-square test).

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

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

Kết quả nghiên cứu thực nghiệm 5 năm tại vùng Đông Bắc Bộ đã mang lại 5 nhóm phát hiện khoa học mang tính bước ngoặt:

┌────────────────────────────────────────────────────────────────────────────────────────┐
│          BẢNG TỔNG HỢP NĂNG SUẤT VÀ TĂNG TRƯỞNG CỦA 6 DÒNG VÔ TÍNH ƯU VIỆT             │
│                      TẠI UÔNG BÍ, QUẢNG NINH (GIAI ĐOẠN 5 TUỔI)                        │
├─────────┬──────────────┬──────────────┬───────────────┬────────────────┬───────────────┤
│ Tên     │ Tỷ lệ sống   │ Đường kính   │ Chiều cao     │ Năng suất gỗ   │ Độ vượt NS so │
│ Dòng    │ TLS (%)      │ D1.3 (cm)    │ Hvn (m)       │ (m³/ha/năm)    │ với TBKN (%)  │
├─────────┼──────────────┼──────────────┼───────────────┼────────────────┼───────────────┤
│ Clt98   │ 93,8         │ 14,2 ± 0,38  │ 12,8 ± 0,27   │ 23,0           │ +39,4         │
│ Clt26   │ 90,6         │ 13,8 ± 0,41  │ 12,4 ± 0,31   │ 20,8           │ +26,1         │
│ Clt57   │ 92,2         │ 13,5 ± 0,36  │ 12,1 ± 0,29   │ 19,4           │ +17,6         │
│ Clt7    │ 93,8         │ 13,0 ± 0,44  │ 11,8 ± 0,33   │ 17,2           │ +4,2          │
│ Clt133  │ 87,5         │ 12,9 ± 0,40  │ 11,6 ± 0,30   │ 16,5           │ +0,0          │
│ AA9     │ 89,1         │ 12,6 ± 0,42  │ 11,5 ± 0,35   │ 15,3           │ -7,3          │
│ TBKN    │ 89,2         │ 12,7 ± 0,40  │ 11,6 ± 0,32   │ 16,5           │ 0,0           │
└─────────┴──────────────┴──────────────┴───────────────┴────────────────┴───────────────┘
  1. Phân hóa năng suất và khả năng thích ứng của các dòng vô tính: Xác định được 6 dòng vô tính Keo lá tràm thích ứng xuất sắc tại Đông Bắc Bộ, đạt năng suất từ 15,3 đến 23,0 m³/ha/năm ở tuổi 5. Nổi bật nhất là dòng Clt98 ($D_{1.3} = 14,2\text{ cm}$; $H_{vn} = 12,8\text{ m}$; năng suất 23,0 m³/ha/năm, vượt 39,4% so với trung bình khảo nghiệm TBKN) và dòng Clt26 ($D_{1.3} = 13,8\text{ cm}$; $H_{vn} = 12,4\text{ m}$; năng suất 20,8 m³/ha/năm, vượt 26,1%). Về chất lượng thân cây, các dòng Clt98, Clt26, Clt57 đạt tỷ lệ cây thân thẳng (loại 1) $>80%$, cành nhánh nhỏ mịn, không bị cụt ngọn trong các đợt rét hại.
  2. Hiệu lực vượt trội của xử lý thực bì không đốt (CT1): Công thức phát trắng, chặt băm nhỏ thực bì $<0,5\text{ m}$ rải đều toàn diện đạt đường kính $D_{1.3} = 13,2\text{ cm}$, chiều cao $H_{vn} = 12,1\text{ m}$, tỷ lệ sống 93,5% và năng suất 18,2 m³/ha/năm; cao hơn rõ rệt ($p < 0,01$) so với công thức phát gom theo băng (16,4 m³/ha/năm) và công thức phát theo rạch 2 m (14,1 m³/ha/năm). Việc giữ ẩm và chống rửa trôi đất dốc của lớp thảm tàn dư băm nhỏ đóng vai trò quyết định.
  3. Đột phá về kỹ thuật làm đất cơ giới kết hợp bậc thang: Công thức làm đất bằng cày ngầm sâu 40 cm, rộng 1,5 m kết hợp cuốc hố $30 \times 30 \times 30\text{ cm}$ (CT1) và làm bậc thang rộng 1,0 m theo đường đồng mức giữ lại tầng đất mặt kết hợp hố $40 \times 40 \times 40\text{ cm}$ (CT2) giúp cây sinh trưởng vượt trội ($D_{1.3} = 13,4\text{ cm}$; năng suất đạt 19,1–19,3 m³/ha/năm) so với cuốc hố thủ công $40 \times 40 \times 40\text{ cm}$ đơn thuần (15,8 m³/ha/năm).
  4. Cơ chế dinh dưỡng khoáng lá và hiệu quả bón phân chuyên biệt: Phân tích lá 12 tháng tuổi chỉ ra hàm lượng phospho trong lá ở công thức đối chứng chỉ đạt 0,11% (rất nghèo), trong khi đạm đạt 2,35% và kali đạt 0,86%. Công thức bón lót tổng hợp 1,0 kg hữu cơ vi sinh Sông Gianh + 0,3 kg NPK (16:16:8)/hố kết hợp bón thúc vào năm 2 và năm 3 với liều lượng 0,1 kg $K_2O$ + 0,2 kg $P_2O_5$/cây hoặc 0,5 kg NPK (16:16:8)/cây làm tăng trữ lượng gỗ lên 24,6% so với đối chứng không bón ($p < 0,001$).
  5. Quy luật phân hóa mật độ và cấu trúc lâm phần tuổi 5: Tại mật độ 1.660 cây/ha, mặc dù trữ lượng tổng số đạt $102,5\text{ m}^3/\text{ha}$, đường kính bình quân chỉ đạt 11,8 cm và bắt đầu xuất hiện hiện tượng tự tỉa thưa, cạnh tranh ánh sáng gay gắt. Trong khi đó, mật độ 1.110 cây/ha (cự ly $3 \times 3\text{ m}$) đạt $D_{1.3} = 13,2\text{ cm}$, tỷ lệ cây có $D_{1.3} \ge 15\text{ cm}$ đạt trên 35%, là mật độ tối ưu nhất để nuôi dưỡng rừng gỗ lớn. Phân tích phân bố $N/D$ chứng minh rừng 5 tuổi mật độ 1.660 cây/ha đòi hỏi phải tỉa thưa 40–50% số cây (chặt cây kiệt bì, cong queo, sâu bệnh) để đưa mật độ về 800–900 cây/ha.

Implications đa chiều

  • Về mặt học thuật: Cung cấp bằng chứng thực nghiệm vững chắc xác nhận tính đúng đắn của mô hình quản lý dinh dưỡng lập địa đa luân kỳ tại vùng á nhiệt đới, làm phong phú thêm cơ sở dữ liệu về tính thích ứng sinh thái của loài Acacia auriculiformis.
  • Về phương pháp luận: Chuẩn hóa quy trình tích hợp giữa chẩn đoán sinh hóa khoáng lá cây rừng với thiết kế thí nghiệm lâm sinh thâm canh nhân tố, mở ra hướng tiếp cận định lượng có độ chính xác cao cho các nghiên cứu lâm nghiệp nhiệt đới.
  • Về ứng dụng thực tiễn: Cung cấp gói kỹ thuật thâm canh hoàn chỉnh (Standard Operating Procedure - SOP) từ khâu chọn dòng vô tính (Clt98, Clt26, Clt57), xử lý mặt bằng, làm đất cơ giới, chế độ phân bón đến kỹ thuật nuôi dưỡng tỉa thưa cho các công ty lâm nghiệp, ban quản lý rừng phòng hộ/đặc dụng và hàng trăm ngàn hộ gia đình trồng rừng tại Quảng Ninh, Bắc Giang, Lạng Sơn, Thái Nguyên, Tuyên Quang.
  • Về chính sách: Là luận cứ khoa học trực tiếp hỗ trợ Bộ Nông nghiệp và Phát triển nông thôn cùng các địa phương điều chỉnh quy hoạch vùng trồng rừng gỗ lớn, ban hành quy trình kỹ thuật quốc gia về trồng thâm canh Keo lá tràm cung cấp gỗ lớn theo Quyết định số 774/QĐ-BNN-TCLN.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Thời gian theo dõi: Thí nghiệm lâm sinh mới dừng lại ở giai đoạn 5 năm tuổi (thời điểm giao thời chuyển sang nuôi dưỡng rừng gỗ lớn), chưa theo dõi trọn vẹn toàn bộ luân kỳ kinh doanh gỗ lớn (10–12 năm).
  2. Phạm vi địa lý thực nghiệm: Hệ thống thí nghiệm chuyên sâu mới được bố trí tập trung tại Uông Bí - Quảng Ninh, cần tiếp tục kiểm chứng mở rộng trên các dạng lập địa đặc thù khác của vùng Đông Bắc Bộ (như đất Feralit trên đá vôi thoái hóa hoặc vùng núi cao nhiệt độ mùa đông $<5^\circ\text{C}$).
  3. Chỉ tiêu cơ lý gỗ: Chưa thực hiện đánh giá các chỉ tiêu công nghệ gỗ xẻ thành phẩm ở tuổi thành thục khai thác chính (như tỷ lệ nứt nẻ khi sấy, mô đun đàn hồi MOE, ứng suất uốn tĩnh MOR và độ co rút thể tích).
  4. Biến động khí hậu cực đoan: Chưa mô phỏng tác động của các hiện tượng thời tiết bất thường (bão mạnh đổ bộ, hạn hán kéo dài hoặc sương muối lịch sử) đối với cấu trúc lâm phần sau tỉa thưa.

Chương trình nghiên cứu tiếp nối (Future Research Agenda) cần tập trung vào:

  • Hướng 1: Tiếp tục đo đếm, theo dõi động thái tăng trưởng của lâm phần thí nghiệm từ tuổi 6 đến tuổi 12; đánh giá hiệu lực của các cường độ tỉa thưa nuôi dưỡng (giữ lại 600, 800, 1.000 cây/ha).
  • Hướng 2: Phân tích đánh giá chất lượng gỗ cơ lý, tỷ lệ lợi dụng gỗ xẻ và hiệu quả kinh tế - tài chính (NPV, IRR, B/C) của chuỗi giá trị gỗ xẻ đóng đồ mộc cao cấp từ Keo lá tràm thâm canh so với Keo tai tượng và Keo lai.
  • Hướng 3: Ứng dụng công nghệ GIS, viễn thám và học máy (Machine Learning) để xây dựng bản đồ phân vùng lập địa thích ứng tối ưu cho Keo lá tràm trên toàn vùng Đông Bắc Bộ.
  • Hướng 4: Nghiên cứu sâu về cân bằng carbon (carbon footprint), khả năng hấp thụ $CO_2$ và cải thiện độ phì sinh học đất của mô hình thâm canh Keo lá tràm phục vụ thị trường tín chỉ carbon rừng.

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

Nghiên cứu của NCS. Phạm Đình Sâm tạo ra những tác động lan tỏa mạnh mẽ:

  • Tác động học thuật: Các kết quả của luận án đã được công bố trên các tạp chí khoa học chuyên ngành uy tín (Tạp chí Khoa học Lâm nghiệp), ước tính tạo lập nguồn trích dẫn quan trọng cho các nghiên cứu sinh, học viên cao học và chuyên gia lâm học trong và ngoài nước nghiên cứu về chi Acacia.
  • Tái cấu trúc ngành lâm nghiệp: Tạo bước ngoặt thay đổi tập quán sản xuất từ kinh doanh gỗ nhỏ làm dăm xuất khẩu giá trị thấp (50–60 USD/tấn) sang sản xuất gỗ lớn đóng đồ gỗ gia dụng cao cấp với giá trị thương phẩm cao gấp 2,5–3,5 lần trên một đơn vị diện tích chu kỳ 10–12 năm.
  • Ảnh hưởng chính sách công: Cung cấp luận cứ kỹ thuật phục vụ việc sửa đổi, bổ sung Tiêu chuẩn Quốc gia TCVN 11567-1:2016 về rừng trồng keo gỗ lớn và hoàn thiện Đề án Phát triển kinh tế lâm nghiệp bền vững tại các tỉnh miền núi phía Bắc.
  • Lợi ích xã hội và môi trường: Nâng cao thu nhập bền vững cho đồng bào dân tộc thiểu số và người dân làm nghề rừng tại các tỉnh Đông Bắc; bảo vệ môi trường sinh thái thông qua việc giữ lại thảm phủ, hạn chế xói mòn rửa trôi hàng triệu tấn đất mặt và tăng trữ lượng hấp thụ carbon rừng trồng.
  • Tầm vóc quốc tế: Đóng góp một điển hình thành công (case-study) về chuyển dịch giống và lâm sinh thích ứng cận nhiệt đới vào mạng lưới nghiên cứu keo nhiệt đới quốc tế (IUFRO, ACIAR, CIFOR).

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

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                        CƠ CẤU ĐỐI TƯỢNG HƯỞNG LỢI TỪ LUẬN ÁN                           │
├──────────────────────────┬─────────────────────────────────────────────────────────────┤
│ Nhóm đối tượng           │ Giá trị và Lợi ích cụ thể tiếp nhận                         │
├──────────────────────────┼─────────────────────────────────────────────────────────────┤
│ Nghiên cứu sinh &        │ Tiếp cận khung phân tích lâm sinh định lượng liên hoàn;     │
│ Học viên Cao học         │ cơ sở dữ liệu thực nghiệm chuẩn mực về biến dị di truyền    │
│                          │ và động thái lập địa của chi Acacia.                        │
├──────────────────────────┼─────────────────────────────────────────────────────────────┤
│ Các nhà Khoa học         │ Kế thừa luận cứ lý thuyết về tương tác G x E x S; phương     │
│ Lâm nghiệp & Thổ nhưỡng  │ pháp chẩn đoán khoáng lá điều tiết phân bón chuyên biệt     │
│                          │ cho cây gỗ họ Đậu trên đất đồi dốc.                         │
├──────────────────────────┼─────────────────────────────────────────────────────────────┤
│ Doanh nghiệp &           │ Tiếp nhận 6 dòng vô tính vượt trội (Clt98, Clt26, Clt57...) │
│ Công ty Lâm nghiệp       │ và quy trình công nghệ thâm canh, nâng năng suất rừng gỗ    │
│                          │ lớn từ 10-12 m³ lên 18-23 m³/ha/năm, gia tăng lợi nhuận.    │
├──────────────────────────┼─────────────────────────────────────────────────────────────┤
│ Cơ quan Quản lý          │ Bộ tài liệu căn cứ khoa học chuẩn xác phục vụ thẩm định,    │
│ & Hoạch định Chính sách  │ phê duyệt quy hoạch rừng gỗ lớn và xây dựng chính sách hỗ   │
│                          │ trợ tín dụng lâm nghiệp dài hạn.                            │
└──────────────────────────┴─────────────────────────────────────────────────────────────┘

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? Trả lời: Đóng góp lý thuyết độc đáo nhất là việc mở rộng và kiểm chứng thực nghiệm thành công Lý thuyết Thâm canh Lâm sinh (Nguyễn Xuân Quát, 1995, 1998) kết hợp Lý thuyết Quản lý Lập địa Bền vững (Gonçalves et al., 1997; Nambiar & Harwood, 2014) trên đối tượng Keo lá tràm tại vùng biên giới sinh thái cận nhiệt đới Đông Bắc Bộ. Luận án đã mô hình hóa toán học mối quan hệ tương tác phức hợp giữa Kiểu gen thích ứng lạnh ($G$) $\times$ Lập địa đồi dốc nghèo lân ($E$) $\times$ Tổ hợp biện pháp kỹ thuật không đốt, cày ngầm và bón đón đầu nhu cầu sinh lý ($S$).

2. Điểm mới đột phá về phương pháp luận nghiên cứu so với các công trình trước đây là gì? Trả lời: So với các nghiên cứu trước đây của Bùi Việt Hải (1998) hay Nguyễn Huy Sơn (2006) vốn chỉ tiếp cận đơn lẻ từng biện pháp (hoặc chỉ khảo nghiệm giống, hoặc chỉ thử nghiệm phân bón), luận án đã:

  • Tích hợp phương pháp chẩn đoán hóa sinh dinh dưỡng khoáng trong lá (N, P, K) theo mùa sinh trưởng để định hướng chính xác công thức phân bón thúc (phát hiện đất thiếu lân nghiêm trọng đối với Keo lá tràm).
  • Thiết kế hệ thống thí nghiệm nhân tố liên hoàn chuẩn mực (RCBD với diện tích đủ lớn, đo đếm vùng lõi 36 cây khử sai số biên) kết hợp khảo sát hồi quy 18 OTC thực địa đa cấp tuổi (4 đến 24 tuổi).

3. Phát hiện bất ngờ nhất (counter-intuitive result) từ số liệu thực nghiệm là gì? Trả lời: Phát hiện bất ngờ nhất là mặc dù Keo lá tràm là loài cây có vi khuẩn nốt sần rễ (Rhizobium) cố định đạm rất mạnh, việc chỉ bón phân NPK thông thường không đem lại hiệu quả tăng trưởng đột phá bằng việc bón kết hợp phân hữu cơ vi sinh với phân lân nung chảy ($P_2O_5$) và kali ($K_2O$). Hàm lượng phospho trong lá ở rừng không bón chỉ đạt 0,11%, chứng minh phospho chính là yếu tố giới hạn sinh thái (ecological limiting factor) quyết định tốc độ sinh trưởng của Keo lá tràm trên đất Feralit Đông Bắc Bộ.

4. Luận án có cung cấp quy trình nhân bản (replication protocol) đầy đủ không? Trả lời: Có, luận án cung cấp quy trình công nghệ lâm sinh cực kỳ chi tiết, chuẩn hóa ở mức độ sẵn sàng chuyển giao (SOP), bao gồm:

  • Tiêu chuẩn cây mô giống: 3 tháng tuổi, đường kính cổ rễ $D_{00} = 0,2-0,3\text{ cm}$, chiều cao $H_{vn} = 25-30\text{ cm}$.
  • Quy cách xử lý thực bì: Chặt băm $<0,5\text{ m}$ rải đều toàn diện, tuyệt đối không đốt.
  • Quy cách làm đất: Cày ngầm sâu 40 cm hoặc làm bậc thang 1,0 m theo đường đồng mức kết hợp hố $40 \times 40 \times 40\text{ cm}$.
  • Định mức phân bón: Bón lót 1,0 kg phân hữu cơ vi sinh + 0,3 kg NPK (16:16:8)/hố; bón thúc 0,1 kg $K_2O$ + 0,2 kg $P_2O_5$/cây hoặc 0,5 kg NPK (16:16:8)/cây vào đầu mùa mưa năm 2 và năm 3.
  • Mật độ trồng ban đầu: 1.110–1.660 cây/ha và quy chế tỉa thưa ở tuổi 5.

5. Luận án đã phác thảo chương trình nghiên cứu 10 năm tiếp theo như thế nào? Trả lời: Luận án đã vạch ra lộ trình 10 năm (2022–2032) với 3 trọng tâm: (i) Tiếp tục quan trắc các ô thí nghiệm đến tuổi 10–12 để xác định tuổi thành thục công nghệ gỗ xẻ và hoàn thiện chế độ tỉa thưa lần 2; (ii) Kiểm định chất lượng gỗ cơ lý thành phẩm (độ bền uốn, tỷ trọng, độ co rút); (iii) Xây dựng mô hình chuỗi cung ứng gỗ xẻ Keo lá tràm gắn chứng chỉ quản lý rừng bền vững (FSC/PEFC) và đánh giá dịch vụ hệ sinh thái hấp thụ carbon.

Kết luận

  1. Luận án đã giải quyết căn cơ bài toán khoa học và thực tiễn về phát triển rừng trồng thâm canh Keo lá tràm (Acacia auriculiformis) cung cấp gỗ lớn tại vùng Đông Bắc Bộ, lấp đầy khoảng trống nghiên cứu sinh thái - lâm sinh tại vùng cận nhiệt đới miền Bắc Việt Nam.
  2. Tuyển chọn thành công 6 dòng vô tính ưu việt (Clt98, Clt26, Clt57, Clt7, Clt133, AA9) sinh trưởng vượt trội, đạt năng suất thể tích thân cây từ 15,3 đến 23,0 m³/ha/năm ở tuổi 5, trong đó dòng Clt98 và Clt26 có phẩm chất thân cây thẳng đẹp, cành nhỏ, kháng chịu rét hại tốt, đáp ứng hoàn hảo tiêu chuẩn quốc gia về rừng gỗ lớn.
  3. Xác lập hệ thống biện pháp kỹ thuật lâm sinh thâm canh đồng bộ và liên hoàn: (i) Phát trắng thực bì băm nhỏ rải đều bảo tồn dinh dưỡng hữu cơ; (ii) Cày ngầm sâu 40 cm hoặc làm bậc thang 1,0 m kết hợp cuốc hố theo đường đồng mức; (iii) Mật độ trồng ban đầu tối ưu 1.110 cây/ha ($3 \times 3\text{ m}$) hoặc 1.660 cây/ha có tỉa thưa ở tuổi 5; (iv) Bón lót 1,0 kg hữu cơ vi sinh + 0,3 kg NPK (16:16:8)/hố kết hợp bón thúc P và K theo chẩn đoán khoáng lá.
  4. Nâng tầm mô hình quản lý rừng trồng từ quảng canh gỗ nhỏ sang thâm canh gỗ lớn giá trị cao, tạo nền tảng mở ra 3 hướng nghiên cứu mới: Công nghệ chế biến gỗ xẻ Keo lá tràm cao cấp; Bản đồ số hóa lập địa thích ứng biến đổi khí hậu; Kinh tế học carbon và dịch vụ chi trả môi trường rừng.
  5. Khẳng định giá trị thực tiễn và tính khả thi chuyển giao công nghệ cao, đóng góp thiết thực vào chiến lược tái cơ cấu ngành lâm nghiệp, nâng cao thu nhập cho người làm rừng và nâng cao vị thế của khoa học lâm nghiệp Việt Nam trên trường quốc tế.