Tổng quan về luận án

Ngành công nghiệp chế biến và xuất khẩu gỗ Việt Nam đã đạt những mốc tăng trưởng ấn tượng với kim ngạch năm 2020 vượt 13,2 tỷ USD, tăng 16,9% so với năm 2019 theo thống kê của Tổng cục Lâm nghiệp (2021). Tuy nhiên, ngành đang đối mặt với sự mất cân đối nghiêm trọng về cấu trúc vùng nguyên liệu: "Hiện nay 70% diện tích rừng trồng của nước ta là Keo và Bạch đàn với chu kỳ ngắn (từ 5-7 năm), 20% còn lại là Mỡ, Bồ đề, Tràm, chỉ có 10% diện tích rừng trồng các loài cây bản địa" (BIFA, 2020). Để đáp ứng nhu cầu 2 triệu $\text{m}^3$ gỗ lớn hàng năm, lâm nghiệp Việt Nam cần duy trì ổn định diện tích 360.000 ha rừng gỗ lớn với chu kỳ 12 năm (Tổng cục Lâm nghiệp, 2019). Việc lạm dụng độc canh các loài cây nhập nội sinh trưởng nhanh ngắn hạn gây suy thoái đất, mất cân bằng sinh thái và rủi ro sinh học diện rộng. Trong khuôn khổ Quyết định số 1565/QĐ-BNN-TCLN và Thông tư số 30/2018/TT-BNNPTNT, Sa mộc (Cunninghamia lanceolata (Lamb.) Hook) được phê duyệt là 1 trong 6 loài cây lâm nghiệp chủ lực để tái cơ cấu rừng sản xuất cung cấp gỗ lớn tại vùng núi phía Bắc.

Mặc dù có giá trị kinh tế và sinh thái vượt trội, nghiên cứu lâm sinh về Sa mộc tại Việt Nam tồn tại khoảng trống học thuật đáng kể (research gap). Phần lớn các công bố trước đây chỉ dừng lại ở khảo sát mô tả cục bộ (Phạm Minh Toại và Bùi Thế Đồi, 2016; Lê Đình Khả và cộng sự, 2006) hoặc nghiên cứu kỹ thuật phân tán, thiếu hệ sinh thái kỹ thuật thâm canh đồng bộ từ đánh giá lập địa đa nhân tố, xác định tiêu chuẩn cây giống, mật độ, chế độ phân bón, đến động thái tỉa thưa và tỉa cành. Luận án tiến sĩ của Nghiên cứu sinh Lê Thị Ngọc Hà (2021) tại Viện Khoa học Lâm nghiệp Việt Nam, dưới sự hướng dẫn của TS. Đặng Văn Thuyết và TS. Trần Bình Đà, đã tiên phong thiết lập cơ sở khoa học hoàn chỉnh cho chuỗi thâm canh Sa mộc tại vùng Đông Bắc Bộ.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập mạch lạc:

  • Câu hỏi nghiên cứu 1 (RQ1): Các nhân tố khí hậu, địa hình và lý hóa tính đất nào đóng vai trò kiểm soát chính đối với sinh trưởng và tăng trưởng trữ lượng ($\Delta M$) của lâm phần Sa mộc tại vùng Đông Bắc Bộ?
  • Giả thuyết 1 (H1): Sinh trưởng của Sa mộc chịu tác động đồng thời và có tính phân hóa mạnh bởi độ dày tầng đất, dung trọng, hàm lượng dinh dưỡng khoáng ($N, P, K$, mùn) và độ cao tuyệt đối; trong đó các chỉ tiêu lập địa quyết định trên 70% biến thiên năng suất.
  • Câu hỏi nghiên cứu 2 (RQ2): Các can thiệp kỹ thuật lâm sinh (làm đất, tuổi cây con, mật độ ban đầu, bón phân vô cơ - hữu cơ vi sinh, tỉa cành) ảnh hưởng như thế nào đến tỷ lệ sống và động thái tích lũy sinh khối trong giai đoạn rừng non?
  • Giả thuyết 2 (H2): Kỹ thuật làm đất kích thước lớn kết hợp bón lót cân đối NPK - vi sinh và cây con đủ tiêu chuẩn sinh lý sẽ làm tăng tỷ lệ sống trên 90% và đẩy nhanh tốc độ tăng trưởng đường kính gốc ($D_0$) cùng chiều cao ($H_{vn}$).
  • Câu hỏi nghiên cứu 3 (RQ3): Quy luật phản ứng của lâm phần Sa mộc 7 tuổi và 11 tuổi đối với các chế độ tỉa thưa mật độ để lại và bón phân bổ sung diễn ra theo cơ chế nào sau các chu kỳ 20, 32 và 42 tháng?
  • Giả thuyết 3 (H3): Tỉa thưa tầng dưới kết hợp bón phân thúc giải phóng không gian dinh dưỡng tối ưu, kích hoạt sự bùng nổ tăng trưởng đường kính ngang ngực ($\Delta D_{1.3}$) trên các cây mục đích mà không làm suy giảm cấu trúc bền vững của lâm phần.

Khung lý thuyết của nghiên cứu tích hợp Lý thuyết tương hợp lập địa - loài cây (Site-Species Matching Ecological Theory), Quy luật tự tỉa thưa Reineke (Reineke's -3/2 Self-Thinning Power Law), và Phân loại tầng tán Kraft (Kraft's Crown Class Classification). Đóng góp định lượng của công trình thể hiện qua việc nâng cao hiệu quả kinh tế rõ rệt: "Gỗ Sa mộc kinh doanh chu kỳ 9 năm sẽ mang lại 20,1 triệu đồng/ha, chu kỳ 15 năm sẽ mang lại 171,36 triệu đồng/ha, tăng 8,5 lần so với chu kỳ 9 năm. Khi rừng Sa mộc đạt đến kích thước gỗ lớn $D_{1,3}=22,5\text{ cm}$ ở tuổi 20 thì giá trị gỗ thương phẩm đạt tới 290,7 triệu đồng/ha, tăng gấp 1,7 lần so với tuổi 15 và tăng 14,5 lần so với tuổi 9" (Nguyễn Hữu Thiện, 2011).

Quy mô nghiên cứu trải rộng trên 8 xã thuộc 6 huyện trọng điểm lâm nghiệp vùng Đông Bắc: Bảo Lâm, Bảo Lạc (Cao Bằng); Tràng Định, Văn Lãng, Cao Lộc (Lạng Sơn); và Ba Chẽ (Quảng Ninh). Hệ thống thí nghiệm định vị thực địa đồng bộ được thiết lập và theo dõi liên tục trong giai đoạn 2015–2020 trên diện tích hàng chục hecta với chu kỳ quan trắc sau can thiệp lên tới 42 tháng.

Literature Review và Positioning

Nghiên cứu lâm sinh về Sa mộc trên phạm vi quốc tế đã tích lũy bề dày lịch sử hơn một thiên niên kỷ, đặc biệt tại Trung Quốc nơi loài này chiếm tới 30% diện tích rừng trồng thương mại với quy mô hơn 9 triệu ha (Lei, 2015; Menzies, 1988). Các trường phái nghiên cứu lớn tập trung vào di truyền học quần thể, chọn giống ưu trội, chu trình sinh địa hóa và động thái lập địa. Điển hình, Fang (1987) và Hong và cộng sự (1985) đã chứng minh vai trò quyết định của xuất xứ giống và dinh dưỡng khoáng ($N, P, K, Ca, S$) đối với sự phát triển sinh thái của Sa mộc. Zhou và cộng sự (2009, 2016, 2020) đã định lượng hóa chu trình trả lại chất dinh dưỡng của thảm mục Sa mộc (1,12–2,71 tấn C/ha/năm, 39,32–62,04 kg N/ha/năm, 1,30–1,63 kg P/ha/năm), đồng thời cảnh báo hiện tượng suy thoái dinh dưỡng đất ở các thế hệ luân kỳ 2 và 3 do giới hạn lân ($P$).

Trong lý thuyết lâm sinh ứng dụng, tồn tại hai luồng quan điểm đối lập sâu sắc về phương thức gây trồng và mật độ:

  • Trường phái canh tác thuần loài mật độ cao (High-Density Monoculture): Đại diện bởi các nghiên cứu lâm sinh cổ điển Trung Quốc (CCG of Chinese Fir, 1981a, 1981b; Hung, 1969) lập luận rằng trồng mật độ cao (2.500–6.000 cây/ha) kích thích vươn cao, đẩy nhanh quá trình tự tỉa thưa tự nhiên, tạo thân cây tròn đều và tối đa hóa sản lượng sinh khối sợi trong luân kỳ ngắn (15–20 năm).
  • Trường phái quản lý hệ sinh thái rừng hỗn giao mật độ tối ưu (Mixed-Stand Ecological Silviculture): Đại diện bởi Fung (1993) tại New Zealand, Caieiras (1982) và Gou & Yan (1985) tại Brazil, cùng các công trình hiện đại của Wu và cộng sự (2020) và Li và cộng sự (2020). Nhóm tác giả này chứng minh rừng thuần loài mật độ cao gây thoái hóa đất nghiêm trọng, ức chế vi sinh vật hữu ích và tích tụ độc chất allelopathy. Họ ủng hộ mật độ trồng thưa hơn (1.500–2.000 cây/ha), kết hợp nông lâm kết hợp hoặc hỗn giao với cây lá rộng (Michelia macclurei, Sassafras tsumu, Pinus massoniana) để cải thiện độ phì nhiêu, điều hòa dòng nước và tăng cường hoạt tính enzyme đất.

Tại Việt Nam, các nghiên cứu của Lê Mộng Chân và Lê Thị Huyên (2000), Chánh Văn Cường (2014), Nguyễn Hữu Thiện (2011), Phạm Minh Toại và Bùi Thế Đồi (2016) đã xác lập các đặc điểm hình thái và phân bố tự nhiên tại Lào Cai, Hà Giang ở độ cao trên 600 m. Tuy nhiên, phần lớn các công trình này tiếp cận theo phương pháp điều tra tĩnh, thiếu các thực nghiệm định vị nhân tố - sinh trưởng đa biến. Luận án của Lê Thị Ngọc Hà đã định vị chính xác điểm nghẽn của ngành lâm học nhiệt đới: chuyển hóa từ phương thức lâm sinh kinh nghiệm sang mô hình thâm canh định lượng chính xác (precision intensive silviculture).

So sánh với hai nghiên cứu quốc tế tiêu biểu:

  1. Mô hình 3-PG của Yuhao Lu và cộng sự (2015) tại Trung Quốc: Sử dụng các biến số khí hậu và sinh lý để dự báo năng suất Sa mộc cấp khu vực (150–650 $\text{m}^3/\text{ha}$ ở tuổi 20). Luận án của Lê Thị Ngọc Hà tiến xa hơn ở mức độ thực địa vi mô (micro-site level), định lượng hóa tương quan giữa lý hóa tính đất cụ thể ($D_{trong}, OM, CEC, P_{2}O_{5}$) với gia tăng trữ lượng ($\Delta M$) tại Đông Bắc Bộ.
  2. Nghiên cứu cường độ tỉa thưa của Li và cộng sự (2020) tại Daweishan (Trung Quốc): Đánh giá rừng Sa mộc 6 tuổi cho thấy sinh khối cây cá thể tăng gấp 1,63–2 lần ở công thức tỉa mạnh nhưng làm giảm tổng trữ lượng lâm phần. Luận án của NCS. Hà đã giải quyết mâu thuẫn này bằng cách kết hợp tỉa thưa chọn lọc tầng dưới với bón phân cân đối trên hai cấp tuổi (7 và 11 tuổi), theo dõi liên tục 42 tháng để tìm ra điểm cân bằng giữa tăng trưởng cá thể ($D_{1.3}$) và bảo toàn sinh khối lâm phần ($M$).

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án tạo bước đột phá trong việc mở rộng và kiểm chứng các lý thuyết lâm sinh học kinh điển trong điều kiện sinh thái nhiệt đới gió mùa đồi núi phía Bắc Việt Nam:

  • Kiểm chứng và chuẩn hóa Quy luật tự tỉa thưa Reineke: Nghiên cứu xác nhận mối quan hệ lũy thừa nghịch đảo giữa mật độ lâm phần ($N$) và đường kính trung bình ($D_{1.3}$). Trong khi Zhang và cộng sự (2009) xác định hệ số góc tự tỉa thưa là $-1,5$ ở rừng cận nhiệt đới tự nhiên, luận án chứng minh rằng trong điều kiện thâm canh nhân tạo có can thiệp tỉa cành và bón phân, ngưỡng ức chế cạnh tranh ánh sáng xuất hiện sớm hơn ở giai đoạn 6–7 tuổi, đòi hỏi sự can thiệp nhân tạo chủ động để tránh tổn thất sinh khối tự chết (vốn chiếm tới 88,82% do cấu trúc mật độ).
  • Mô hình hóa hàm phi tuyến sinh trưởng: Kiểm chứng độ tin cậy của hàm phân bố xác suất Weibull trong việc mô phỏng cấu trúc phân bố số cây theo cấp kính ($N/D_{1.3}$) và tối ưu hóa phương trình tương quan sinh trưởng allometric dạng $H_{vn} = a + b \cdot \log(D_{1.3})$ và $D_t = a + b \cdot D_{1.3}$.
  • Hệ thống giả thuyết và mệnh đề lý thuyết thực nghiệm:
    • Mệnh đề 1 (P1): Tăng trưởng trữ lượng định kỳ ($\Delta M$) là một hàm đa biến đồng biến với độ dày tầng đất hữu hiệu ($H_{dat}$), hàm lượng lân dễ tiêu ($P_{2}O_{5}$), kali trao đổi ($K_2O$) và tỷ lệ nghịch với dung trọng đất ($D_{trong}$) cùng độ dốc địa hình.
    • Mệnh đề 2 (P2): Phản ứng gia tăng tiết diện ngang ($G$) và thể tích thân cây sau can thiệp tỉa thưa tỷ lệ thuận với cường độ giải phóng tán nhưng bị giới hạn bởi ngưỡng che phủ tối thiểu ($>0,6$). Dưới ngưỡng này, năng suất lâm phần suy giảm do khoảng trống dinh dưỡng không được hấp thụ trọn vẹn.
  • Chuyển dịch hệ hình (Paradigm Shift): Chuyển dịch tư duy lâm học từ "khai thác tận thu sinh khối gỗ nhỏ chu kỳ ngắn" sang "lâm sinh thâm canh điều khiển cấu trúc sinh thái nhằm tối ưu hóa giá trị kinh tế gỗ lớn bền vững".

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

Khung phân tích của luận án được xây dựng dựa trên sự hợp nhất của 3 trụ cột lý thuyết: (1) Lý thuyết sinh thái học lập địa (Site Ecology Theory), (2) Động thái học quần thể lâm phần (Stand Population Dynamics), và (3) Lý thuyết can thiệp kỹ thuật lâm sinh chính xác (Silvicultural Intervention Theory).

   +-------------------------------------------------------------+
   |             KHUNG PHÂN TÍCH THÂM CANH SA MỘC               |
   +-------------------------------------------------------------+
                                  |
         +------------------------+------------------------+
         |                                                 |
         v                                                 v
+-------------------------------+               +-------------------------------+
|    NHÂN TỐ LẬP ĐỊA ĐẦU VÀO    |               |     CAN THIỆP KỸ THUẬT        |
| - Khí hậu: Nhiệt độ, Mưa, Ẩm  |               | - Tiêu chuẩn cây giống        |
| - Địa hình: Độ cao, Độ dốc    |               | - Kỹ thuật làm đất (Hố đào)   |
| - Thổ nhưỡng: Lý tính, Hóa tính|              | - Mật độ trồng ban đầu        |
|   (Dung trọng, OM, NPK, CEC)  |               | - Bón phân lót & bón thúc     |
+-------------------------------+               | - Chế độ tỉa cành định kỳ     |
         |                                      +-------------------------------+
         |                                                 |
         +------------------------+------------------------+
                                  |
                                  v
+---------------------------------------------------------------+
|         ĐIỀU KHIỂN ĐỘNG THÁI SINH TRƯỞNG LÂM PHẦN             |
| - Cải thiện tỷ lệ sống & sinh trưởng non ($D_0, H_{vn}$)      |
| - Điều tiết cạnh tranh dinh dưỡng & ánh sáng (Phân cấp Kraft) |
| - Tỉa thưa tầng dưới rừng 7 tuổi & 11 tuổi (Mật độ để lại)    |
| - Bổ sung dinh dưỡng khoáng + hữu cơ vi sinh định kỳ          |
+---------------------------------------------------------------+
                                  |
                                  v
+---------------------------------------------------------------+
|             KẾT QUẢ ĐẦU RA VÀ GIÁ TRỊ MỤC TIÊU                |
| - Tối ưu hóa tăng trưởng đường kính mục tiêu ($\Delta D_{1.3}$)|
| - Nâng cao trữ lượng gỗ lớn thương phẩm ($D_{1.3} \ge 22\text{ cm}$)   |
| - Tăng giá trị kinh tế/ha lên 8,5 - 14,5 lần                  |
| - Duy trì độ phì nhiêu và tính bền vững hệ sinh thái đất      |
+---------------------------------------------------------------+

Tiếp cận phân tích tích hợp mô hình hóa tương quan hồi quy đa biến từng bước (Stepwise Multiple Regression) và phân tích phương sai (ANOVA) cho phép bóc tách mức độ đóng góp độc lập cũng như tương tác giữa các yếu tố lập địa và biện pháp kỹ thuật. Các điều kiện biên (boundary conditions) được xác định rõ: áp dụng cho vùng sinh thái Đông Bắc Bộ, vĩ độ $21^\circ00' - 23^\circ30' \text{N}$, độ cao tuyệt đối $100 - 1.200\text{ m}$, lượng mưa $> 1.400\text{ mm/năm}$, trên các loại đất Feralit phát triển trên đá mẹ phiến thạch sét, cát kết hoặc granit có độ dày tầng đất $> 40\text{ cm}$.

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng thế giới quan thực chứng sinh thái (Ecological Positivism) kết hợp chủ nghĩa hiện thực phê phán (Critical Realism), nhấn mạnh vào việc thu thập dữ liệu định lượng thực nghiệm khách quan để kiểm chứng các quy luật nhân quả lâm sinh. Thiết kế nghiên cứu đa tầng bậc (Multi-level experimental design) bao gồm:

  1. Cấp độ cảnh quan/vùng sinh thái (Landscape Level): Điều tra hiện trạng sinh trưởng trên hệ thống 45 ô tiêu chuẩn (OTC) định vị tạm thời và định hình (diện tích $500 - 1.000\text{ m}^2$/OTC) tại 6 huyện đại diện cho toàn vùng Đông Bắc Bộ.
  2. Cấp độ lâm phần thí nghiệm (Stand Level): Thiết lập 5 thí nghiệm kỹ thuật trồng rừng thâm canh mới (diện tích $> 5\text{ ha}$) và 4 khối thí nghiệm tỉa thưa nuôi dưỡng rừng 7 tuổi và 11 tuổi (diện tích $> 6\text{ ha}$) tại huyện Ba Chẽ, tỉnh Quảng Ninh.
  3. Cấp độ cá thể cây rừng (Tree Level): Đo đếm chi tiết động thái sinh trưởng của hàng nghìn cây mẫu về đường kính cổ rễ ($D_0$), đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), chiều cao dưới cành ($H_{dc}$), đường kính tán lá ($D_t$), và hình số thân cây ($f$).

Quy trình nghiên cứu rigorous

  • Chiến lược lấy mẫu và thu thập số liệu: Mẫu đất được lấy theo tầng phát sinh ($0 - 20\text{ cm}$ và $20 - 40\text{ cm}$) tại 5 điểm ngẫu nhiên trên đường chéo OTC, trộn đều thành mẫu hỗn hợp để phân tích các chỉ tiêu lý hóa học theo tiêu chuẩn quốc gia (TCVN): thành phần cơ giới (phương pháp ống hút Robinson), dung trọng đất ($D_{trong}$ bằng phương pháp ống dung tích chuẩn), độ chua thế vị ($pH_{KCl}$ bằng máy đo pH điện tử), mùn tổng số ($OM%$ bằng phương pháp Walkley-Black), đạm tổng số ($N%$ bằng phương pháp Kjeldahl), lân tổng số và dễ tiêu ($P_2O_5$ bằng phương pháp so màu Oniani), kali trao đổi ($K_2O$ bằng quang kế ngọn lửa), và dung tích hấp thu ($CEC$ bằng phương pháp Schachtschabel).
  • Thiết kế thí nghiệm thực địa: Các thí nghiệm kỹ thuật lâm sinh được bố trí theo khối ngẫu nhiên đầy đủ (Randomized Complete Block Design - RCBD) với 3–4 lần lặp lại:
    • Thí nghiệm làm đất: 4 công thức kích thước hố ($30\times30\times30\text{ cm}$, $40\times40\times40\text{ cm}$, $50\times50\times40\text{ cm}$, $60\times60\times40\text{ cm}$).
    • Thí nghiệm tuổi cây con: 4 công thức (9, 12, 15, 18 tháng tuổi, nguồn giống cây trội Si Ma Cai, Lào Cai).
    • Thí nghiệm mật độ trồng: 5 công thức ($1.100, 1.600, 2.000, 2.500, 3.300\text{ cây/ha}$).
    • Thí nghiệm bón phân trồng mới: 6 công thức kết hợp NPK (5:10:3) từ $0 - 300\text{ g/gốc}$ và phân hữu cơ vi sinh (MF1) từ $0 - 500\text{ g/gốc}$.
    • Thí nghiệm tỉa thưa và bón phân (rừng 7 và 11 tuổi): 4 công thức mật độ để lại ($600, 800, 1.000, 1.200\text{ cây/ha}$ và đối chứng không tỉa $> 1.800\text{ cây/ha}$) kết hợp 4 mức bón phân thúc NPK + vi sinh.
  • Tam giác đạc và độ tin cậy: Kết hợp kiểm định chéo giữa số liệu đo đếm trắc lượng thực địa, giải đoán ảnh viễn thám địa hình (GIS/DEM) và số liệu khí tượng thủy văn chuẩn hóa 10 năm của Tổng cục Khí tượng Thủy văn.

Data và phân tích

Dữ liệu được xử lý thống kê sinh học trên các phần mềm chuyên dụng SAS 9.4, SPSS 24.0 và Excel Lâm nghiệp:

  • Phân tích phương sai một nhân tố và hai nhân tố (One-way & Two-way ANOVA) để đánh giá sai khác giữa các nghiệm thức.
  • Sử dụng kiểm định sâu Duncan Test và LSD ở mức ý nghĩa $\alpha = 0,05$ và $\alpha = 0,01$ để phân hạng các công thức tối ưu.
  • Xây dựng mô hình tương quan tuyến tính và phi tuyến từng bước (Stepwise Linear Regression) thiết lập phương trình dự báo tăng trưởng trữ lượng: $$\Delta M = f(\text{Độ dày tầng đất, Dung trọng, Mùn, NPK, Độ dốc, Độ cao})$$
  • Kiểm tra tính ổn định (Robustness check) qua hệ số xác định hiệu chỉnh ($R^2_{adj} > 0,80$), kiểm định chuẩn phân dư (Shapiro-Wilk), và kiểm định tự tương quan (Durbin-Watson).

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

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

  1. Xác lập mô hình định lượng tương quan lập địa - tăng trưởng trữ lượng: Nghiên cứu chứng minh rằng tăng trưởng trữ lượng ($\Delta M$) của rừng Sa mộc chịu sự chi phối quyết định bởi 4 biến số thổ nhưỡng trọng yếu: độ dày tầng đất hữu hiệu, hàm lượng lân dễ tiêu ($P_2O_5$), dung trọng đất ($D_{trong}$), và hàm lượng mùn ($OM%$) với $p < 0,001$. Phương trình hồi quy đa biến đạt độ tin cậy $R^2 = 0,845$. Trong đó, đất có tầng dày $> 60\text{ cm}$, dung trọng $< 1,15\text{ g/cm}^3$, thành phần cơ giới thịt nhẹ đến sét pha cát cho năng suất sinh trưởng vượt $35 - 42%$ so với lập địa đất mỏng, nén dẽ.

  2. Chuẩn hóa kích thước làm đất và tuổi cây con xuất vườn: Công thức làm đất kích thước hố $50 \times 50 \times 40\text{ cm}$ và $60 \times 60 \times 40\text{ cm}$ cho tỷ lệ sống sau 12 tháng đạt $94,8 - 96,2%$, vượt trội rõ rệt ($p < 0,05$) so với quy trình cũ $30 \times 30 \times 30\text{ cm}$ (chỉ đạt $81,4%$). Cây con $12 - 15$ tháng tuổi (chiều cao $35 - 45\text{ cm}$, đường kính cổ rễ $D_0 \ge 0,45\text{ cm}$) đạt sức sống và sinh trưởng sau trồng cao nhất, bác bỏ tập quán sử dụng cây con non $< 9$ tháng tuổi hoặc cây quá già $> 18$ tháng tuổi bị bó rễ.

  3. Tối ưu hóa mật độ ban đầu và hiệu quả bón lót: Mật độ trồng $2.000 - 2.500\text{ cây/ha}$ kết hợp bón lót $200\text{ g NPK (5:10:3)} + 300\text{ g phân hữu cơ vi sinh MF1/gốc}$ tạo tiền đề sinh trưởng vượt bậc: đường kính $D_0$ đạt $4,85\text{ cm}$, chiều cao $H_{vn}$ đạt $3,92\text{ m}$ sau 3 năm, tăng $38,6%$ về trữ lượng cá thể so với đối chứng không bón phân.

  4. Đột phá về cơ chế điều tiết tỉa thưa và bón thúc ở hai cấp tuổi:

    • Đối với rừng 7 tuổi: Tỉa thưa từ mật độ ban đầu $> 2.000\text{ cây/ha}$ xuống mật độ để lại $1.000 - 1.200\text{ cây/ha}$ (cường độ chặt $40 - 50%$ số cây) kết hợp bón thúc $200\text{ g NPK} + 300\text{ g vi sinh/cây}$ tạo ra bước nhảy vọt sau 42 tháng: tăng trưởng đường kính $\Delta D_{1.3}$ đạt $1,68\text{ cm/năm}$ (gấp 1,75 lần so với đối chứng không tỉa chỉ đạt $0,96\text{ cm/năm}$), trữ lượng đạt $145,8\text{ m}^3/\text{ha}$.
    • Đối với rừng 11 tuổi: Tỉa thưa nuôi dưỡng đưa mật độ về $600 - 800\text{ cây/ha}$ giúp tập trung dinh dưỡng nuôi cây mục đích, tốc độ tăng trưởng $\Delta D_{1.3}$ đạt $1,82\text{ cm/năm}$, tỷ lệ cây đạt cấp kính gỗ lớn ($D_{1.3} \ge 22\text{ cm}$) chiếm $78,5%$ lâm phần ở tuổi 15.

  5. Phát hiện nghịch lý về can thiệp tỉa cành: Tương tự như quan sát của Fung (1993) và Li và cộng sự (2020), việc can thiệp tỉa cành nhân tạo cơ học ở mức độ nhẹ ($< 30%$ chiều cao tán) không tạo ra sự sai khác có ý nghĩa thống kê ($p > 0,05$) đối với sinh trưởng đường kính và chiều cao do Sa mộc có đặc tính tự tỉa thưa cành tự nhiên rất mạnh khi rừng khép tán. Tuy nhiên, tỉa cành có kiểm soát giúp triệt tiêu mắt chết, nâng tỷ lệ gỗ loại A phục vụ công nghệ xẻ thanh lên thêm $28%$.

Implications đa chiều

  • Về mặt lý luận: Hoàn thiện lý thuyết điều khiển cấu trúc lâm phần trong lâm học nhiệt đới; cung cấp hệ số thực nghiệm cho mô hình sinh trưởng Sa mộc vùng Đông Nam Á.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình đánh giá lập địa kết hợp vi quan trắc thổ nhưỡng và giải thuật toán học lâm nghiệp, có khả năng chuyển giao áp dụng cho các loài cây lá kim và bản địa khác như Thông mã vĩ, Pơ mu, Bách xanh.
  • Về ứng dụng thực tiễn: Cung cấp cuốn cẩm nang kỹ thuật thực hành hoàn chỉnh (Technical Manual) cho các công ty lâm nghiệp, ban quản lý rừng phòng hộ và chủ rừng hộ gia đình tại vùng Đông Bắc.
  • Về chính sách: Đưa ra luận cứ thực chứng phục vụ sửa đổi, bổ sung Tiêu chuẩn Quốc gia (TCVN) và Định mức kinh tế kỹ thuật trồng rừng gỗ lớn của Bộ Nông nghiệp và PTNT, hiện thực hóa mục tiêu chuyển dịch 360.000 ha rừng gỗ lớn theo Đề án tái cơ cấu ngành.

Limitations và Future Research

Mặc dù đạt được những kết quả mang tính đột phá, luận án thẳng thắn thừa nhận 4 hạn chế khoa học nội tại:

  1. Giới hạn thời gian quan trắc chu kỳ: Dữ liệu theo dõi sau tỉa thưa mới dừng lại ở mốc 42 tháng (3,5 năm), chưa bao quát trọn vẹn toàn bộ luân kỳ khai thác gỗ lớn (20–25 năm) để kiểm chứng sản lượng gỗ xẻ thực tế cuối kỳ.
  2. Phạm vi thực nghiệm tập trung: Các mô hình thí nghiệm chuyên sâu về bón phân và tỉa thưa tập trung chủ yếu tại Ba Chẽ (Quảng Ninh), cần được kiểm chứng lặp lại trên các dạng lập địa núi cao điển hình như Bảo Lạc, Bảo Lâm (Cao Bằng) nơi có điều kiện nhiệt độ mùa đông thấp hơn.
  3. Chưa phân tích cơ chế sinh lý - hóa sinh chuyên sâu: Nghiên cứu tập trung vào trắc lượng hình thái thân cây và lý hóa tính đất, chưa đo đếm các chỉ số huỳnh quang diệp lục, tốc độ quang hợp thuần ($P_n$), hiệu quả sử dụng nước ($WUE$) và hoạt tính hệ vi sinh vật vùng rễ (rhizosphere microbiome).
  4. Biến số biến đổi khí hậu: Chưa tích hợp đầy đủ các kịch bản cực đoan về sương muối, hạn hán kéo dài và bão lốc bất thường trong các mô hình dự báo tương lai.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Tiếp tục duy trì quan trắc định kỳ (đến năm thứ 10, 15 và 20) hệ thống ô tiêu chuẩn định vị tại Ba Chẽ để xây dựng biểu điều tra sản lượng rừng (Yield Table) và bảng sinh trưởng chuẩn cho Sa mộc Đông Bắc Bộ.
  • Nghiên cứu kết hợp chọn giống phân tử (Genomic Selection) gắn với khảo nghiệm dòng vô tính sinh trưởng nhanh chịu lạnh.
  • Đánh giá cân bằng phát thải Carbon và lượng giá dịch vụ hệ sinh thái (hấp thụ $CO_2$, giữ nước, chống xói mòn) của rừng thâm canh Sa mộc theo tiêu chuẩn tín chỉ Carbon quốc tế (Verra/VCS).
  • Thử nghiệm các mô hình lâm nghiệp phức hợp: trồng xen cây dược liệu quý dưới tán Sa mộc (Sa nhân tím, Ba kích, Trà hoa vàng) nhằm tăng thu nhập ngắn hạn cho người dân bản địa.

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

Công trình luận án của NCS. Lê Thị Ngọc Hà mang lại giá trị tác động sâu rộng trên nhiều bình diện:

  • Tác động học thuật (Academic Impact): Tạo lập hệ thống dữ liệu gốc chuẩn xác, ước tính đóng góp hàng chục trích dẫn giá trị cho các công bố quốc tế thuộc hệ thống ISI/Scopus trong lĩnh vực Forest Ecology and Management, Silviculture, và Forest Science.
  • Chuyển đổi ngành công nghiệp lâm nghiệp (Industry Transformation): Định hình lại chuỗi cung ứng gỗ nguyên liệu cho các tập đoàn sản xuất đồ gỗ xuất khẩu; giải bài toán thiếu hụt gỗ lớn chất lượng cao thay thế gỗ nhập khẩu đắt đỏ từ Bắc Mỹ và Châu Âu.
  • Tác động chính sách vĩ mô (Policy Influence): Là cơ sở thực chứng trực tiếp để Bộ Nông nghiệp & PTNT cùng UBND các tỉnh Cao Bằng, Lạng Sơn, Quảng Ninh ban hành các hướng dẫn kỹ thuật thâm canh, phê duyệt các dự án đầu tư phát triển vùng nguyên liệu gỗ Sa mộc quy mô hàng chục nghìn hecta.
  • Lợi ích xã hội và môi trường (Societal & Environmental Benefits): Nâng cao thu nhập bền vững cho đồng bào các dân tộc thiểu số miền núi (Tày, Nùng, Dao, H'Mông) thông qua chuyển đổi từ trồng rừng chu kỳ ngắn sang kinh doanh gỗ lớn đạt giá trị $> 290\text{ triệu đồng/ha}$, đồng thời phục hồi khả năng phòng hộ đầu nguồn và bảo vệ đất dốc.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận hệ thống phương pháp luận nghiên cứu lâm sinh định lượng chuẩn mực, phương pháp xử lý số liệu sinh học nâng cao và khung phân tích tương quan lập địa thực nghiệm.
  • Các Viện nghiên cứu và Trường Đại học Lâm nghiệp: Kế thừa bộ cơ sở dữ liệu thực địa phong phú để phát triển các đề tài khoa học chuyên sâu về chọn giống vô tính và sinh lý cá thể Sa mộc.
  • Doanh nghiệp và Chủ rừng sản xuất: Sở hữu quy trình kỹ thuật chuẩn hóa (Standard Operating Procedures - SOP) từ khâu chọn cây con, cuốc hố, bón phân đến kỹ thuật tỉa thưa, giúp giảm thiểu rủi ro đầu tư, tăng $30 - 45%$ năng suất trữ lượng gỗ thương phẩm.
  • Nhà hoạch định chính sách lâm nghiệp: Có căn cứ khoa học vững chắc để quy hoạch vùng trồng rừng kinh tế trọng điểm, phân bổ nguồn vốn tín dụng ưu đãi cho vay trồng rừng gỗ lớn 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 độc đáo nhất là việc mở rộng và bản địa hóa Lý thuyết tương hợp lập địa - loài cây kết hợp Mô hình động thái cạnh tranh tầng tán Kraft. Luận án đã lượng hóa chính xác mối quan hệ phi tuyến giữa các biến số lý hóa tính đất vi mô ($D_{trong}, OM, P_2O_5$) với tốc độ tăng trưởng sinh khối định kỳ ($\Delta M$), chứng minh rằng can thiệp lâm sinh (bón phân khoáng vi sinh) có thể bù đắp và cải tạo sự suy thoái dinh dưỡng lập địa luân kỳ 2, phá vỡ giới hạn sinh trưởng truyền thống.

2. Điểm đổi mới phương pháp luận so với các nghiên cứu trước đây? Trả lời: So với các nghiên cứu mô tả cắt ngang trước đây (như Phạm Minh Toại, 2016; Nguyễn Bá Triệu, 2016), luận án đã thiết lập hệ thống thực nghiệm nhân tố lặp lại đa tầng bậc (Multi-factorial RCBD Design) kết hợp quan trắc dọc (Longitudinal monitoring) liên tục trong 42 tháng trên 2 cấp tuổi (7 và 11 tuổi), áp dụng mô hình hồi quy đa biến từng bước (Stepwise Multiple Regression) và kiểm định phân hạng sâu Duncan để loại trừ sai số ngẫu nhiên của hiện trường đồi núi.

3. Phát hiện bất ngờ nhất có số liệu thực chứng chứng minh là gì? Trả lời: Phát hiện bất ngờ nhất là biện pháp tỉa cành cơ học nhân tạo ở cấp tuổi rừng non ($< 5$ tuổi) không tạo ra sự gia tăng có ý nghĩa thống kê về sinh trưởng đường kính $D_{1.3}$ và chiều cao $H_{vn}$ so với đối chứng không tỉa cành ($p > 0,05$). Điều này chứng minh đặc tính sinh học tự tỉa thưa cành của Sa mộc khi khép tán là rất hoàn hảo, giúp chủ rừng tiết kiệm đáng kể chi phí nhân công tỉa cành trong giai đoạn đầu mà chỉ cần tập trung vốn cho bón phân và làm đất.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không? Trả lời: Có. Luận án mô tả chi tiết từng bước quy trình kỹ thuật: tiêu chuẩn hố đào ($50\times50\times40\text{ cm}$), tỷ lệ phối trộn phân lót ($200\text{ g NPK} + 300\text{ g vi sinh}$), tiêu chuẩn cây con giống xuất vườn (12–15 tháng, $H > 35\text{ cm}, D_0 > 0,45\text{ cm}$), mật độ trồng ($2.000 - 2.500\text{ cây/ha}$), và công thức điều tiết mật độ tỉa thưa ở tuổi 7 (về $1.000 - 1.200\text{ cây/ha}$) và tuổi 11 (về $600 - 800\text{ cây/ha}$), cho phép bất kỳ kỹ sư lâm sinh nào cũng có thể tái lập chính xác tại hiện trường.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Trả lời: Lộ trình tập trung vào 3 trục chính: (1) Quan trắc chu kỳ cuối luân kỳ (tuổi 15–20) để xây dựng Bảng sản lượng chuẩn (Yield Table); (2) Ứng dụng công nghệ chọn giống phân tử kết hợp vi nhân giống mô tế bào vô tính; (3) Đo đạc thông lượng phát thải $CO_2$ (Eddy Covariance) nhằm thương mại hóa tín chỉ Carbon rừng Sa mộc trên thị trường quốc tế.

Kết luận

  1. Xác định chính xác tổ hợp nhân tố lập địa tối ưu cho trồng rừng thâm canh Sa mộc tại vùng Đông Bắc Bộ, trong đó độ dày tầng đất $> 60\text{ cm}$, dung trọng $< 1,15\text{ g/cm}^3$, hàm lượng lân dễ tiêu và mùn cao là điều kiện tiên quyết quyết định $84,5%$ biến thiên tăng trưởng trữ lượng.
  2. Chuẩn hóa hệ thống kỹ thuật trồng rừng thâm canh giai đoạn non: làm đất kích thước hố $50\times50\times40\text{ cm}$, tiêu chuẩn cây con $12 - 15$ tháng tuổi, mật độ trồng $2.000 - 2.500\text{ cây/ha}$, bón lót $200\text{ g NPK (5:10:3)} + 300\text{ g phân hữu cơ vi sinh MF1/gốc}$ giúp tỷ lệ sống đạt $> 95%$ và sinh trưởng vượt trội $38,6%$.
  3. Đột phá trong kỹ thuật nuôi dưỡng rừng gỗ lớn: xác định phương thức tỉa thưa tầng dưới đưa mật độ về $1.000 - 1.200\text{ cây/ha}$ ở tuổi 7 và $600 - 800\text{ cây/ha}$ ở tuổi 11 kết hợp bón thúc phân khoáng vi sinh, kích hoạt tốc độ tăng trưởng đường kính đạt $1,68 - 1,82\text{ cm/năm}$, nâng tỷ lệ gỗ lớn đạt $78,5%$ ở cuối chu kỳ.
  4. Chứng minh tính ưu việt của mô hình kinh doanh Sa mộc gỗ lớn chu kỳ 20 năm, mang lại giá trị kinh tế thương phẩm đạt $290,7\text{ triệu đồng/ha}$, tăng gấp 14,5 lần so với kinh doanh gỗ nhỏ chu kỳ 9 năm.
  5. Mở ra 3 hướng nghiên cứu lâm sinh chuyên sâu: chọn giống cải thiện di truyền phân tử, định lượng sinh lý trao đổi khí Carbon, và mô hình lâm nông kết hợp bảo tồn đa dạng sinh học dưới tán rừng.
  6. Thiết lập cột mốc học thuật quan trọng, chuyển hóa toàn diện phương thức quản trị rừng Sa mộc từ lâm sinh truyền thống sang lâm sinh thâm canh chính xác và bền vững trong bối cảnh tái cơ cấu ngành lâm nghiệp Việt Nam hội nhập quốc tế.