Tổng quan về luận án

Ngành Lâm nghiệp Việt Nam đang trải qua quá trình tái cơ cấu mạnh mẽ nhằm nâng cao giá trị gia tăng và tính bền vững sinh thái. Năm 2020, giá trị xuất khẩu lâm sản của Việt Nam xác lập kỷ lục 13,2 tỷ USD, tăng 16,9% so với năm 2019 (Tổng cục Lâm nghiệp, 2021). Tuy nhiên, nút thắt nghiêm trọng của toàn ngành là nguồn cung gỗ lớn chất lượng cao khi mỗi năm cần khoảng 30.000 ha rừng gỗ lớn đưa vào khai thác để cung cấp 2 triệu $m^3$ gỗ nguyên liệu chế biến (Tổng cục Lâm nghiệp, 2019). Thực tế cơ cấu rừng trồng sản xuất hiện nay bộc lộ rủi ro sinh thái sâu sắc: 70% diện tích là Keo (Acacia) và Bạch đàn (Eucalyptus) chu kỳ ngắn (5–7 năm), 20% là Mỡ, Bồ đề, Tràm, và chỉ có khoảng 10% diện tích dành cho các loài cây bản địa (BIFA, 2020). Để hiện thực hóa Quyết định số 1565/QĐ-BNN-TCLN về "Đề án tái cơ cấu ngành Lâm nghiệp" với mục tiêu 40% sản lượng gỗ lớn, và Thông tư số 30/2018/TT-BNNPTNT đưa Sa mộc (Cunninghamia lanceolata (Lamb.) Hook) vào danh mục 6 loài cây lâm nghiệp chủ lực, việc xây dựng cơ sở khoa học đồng bộ cho chuỗi kỹ thuật thâm canh loài cây này là nhiệm vụ đặc biệt cấp thiết.

Khoảng trống nghiên cứu (Research Gap) then chốt tại Việt Nam trước công trình này là sự thiếu hụt các mô hình định lượng xác định mối tương quan giữa tổ hợp lập địa đặc thù vùng Đông Bắc Bộ với sinh trưởng lâm phần, cũng như chưa thiết lập được hệ thống kỹ thuật lâm sinh hoàn chỉnh từ khâu chọn tiêu chuẩn cây giống, kỹ thuật làm đất, mật độ trồng, liều lượng phân bón, tỉa cành cho đến chế độ tỉa thưa nuôi dưỡng chuyển hóa rừng gỗ lớn. Luận án tiến sĩ của Nghiên cứu sinh Lê Thị Ngọc Hà (Viện Khoa học Lâm nghiệp Việt Nam, mã số: 9 62 02 05, hướng dẫn bởi TS. Đặng Văn Thuyết và TS. Trần Bình Đà) đã giải quyết triệt để khoảng trống trên thông qua 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học (Hypotheses - H):

  • RQ1 & H1: Các nhân tố lập địa (khí hậu, địa hình, lý hóa tính đất) tương tác và chi phối như thế nào đến tăng trưởng trữ lượng ($\Delta M$) của Sa mộc tại vùng Đông Bắc? Giả thuyết cho rằng $\Delta M$ phụ thuộc phi tuyến vào độ dày tầng đất, độ dốc và hàm lượng mùn/dung trọng đất.
  • RQ2 & H2: Kỹ thuật làm đất và độ tuổi cây con ảnh hưởng ra sao đến tỷ lệ sống và sinh trưởng ban đầu? Giả thuyết cho rằng cây con 12–14 tháng tuổi kết hợp kích thước hố $40 \times 40 \times 40$ cm tạo ưu thế sinh trưởng rễ vượt trội.
  • RQ3 & H3: Tác động tương hỗ giữa mật độ trồng ban đầu, chế độ bón phân khoáng/vi sinh và kỹ thuật tỉa cành đến động thái sinh trưởng thân cây? Giả thuyết cho rằng mật độ 1.600–2.000 cây/ha kết hợp bón phân cân đối NPK + MF1 tối ưu hóa diện tích dinh dưỡng mà không gây suy thoái sinh trưởng đường kính.
  • RQ4 & H4: Cường độ tỉa thưa và biện pháp bón phân bổ sung ở hai cấp tuổi xung yếu (tuổi 7 và tuổi 11) tác động thế nào đến sự phân hóa cấp kính và tích lũy thể tích thân cây mục đích? Giả thuyết cho rằng tỉa thưa tầng dưới kết hợp bón phân sau tỉa thúc đẩy mạnh mẽ tăng trưởng đường kính ngang ngực ($\Delta D_{1.3}$), rút ngắn chu kỳ kinh doanh gỗ lớn từ 20 năm xuống 15 năm.

Nghiên cứu được xây dựng trên khung lý thuyết sinh thái học quần thể lâm phần, lý thuyết chỉ số lập địa (Site Index) và lý thuyết tỉa thưa cạnh tranh tự nhiên. Về mặt định lượng, nghiên cứu mở ra tiềm năng kinh tế vượt trội: rừng Sa mộc chu kỳ 9 năm chỉ đạt 20,1 triệu đồng/ha, nhưng khi kéo dài chu kỳ lên 15 năm kết hợp tỉa thưa đạt 171,36 triệu đồng/ha (tăng 8,5 lần), và đạt 290,7 triệu đồng/ha ở tuổi 20 với đường kính $D_{1.3} \ge 22,5$ cm (tăng 14,5 lần so với tuổi 9). Phạm vi nghiên cứu bao quát 8 xã thuộc 6 huyện trọng điểm (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; Ba Chẽ - Quảng Ninh) với hệ thống ô tiêu chuẩn (OTC) và các công thức thí nghiệm được theo dõi liên tục đến 42 tháng sau tỉa thưa.

Literature Review và Positioning

Nghiên cứu lâm sinh về Sa mộc (Cunninghamia lanceolata) trên thế giới sở hữu bề dày lịch sử hơn 1.000 năm, đặc biệt tại Trung Quốc nơi Sa mộc chiếm 30% tổng diện tích rừng trồng với hơn 9 triệu ha và cung cấp 20–30% sản lượng gỗ thương mại quốc gia (Menzies, 1988; Lei, 2015; Zhou et al., 2016). Các dòng nghiên cứu kinh điển đã thiết lập các nền tảng quan trọng:

  1. Phân bố sinh thái và lập địa: FAO (1982), Fung (1993) và Yin et al. (2010) xác định Sa mộc thích hợp với vùng cận nhiệt đới ẩm, nhiệt độ trung bình năm 15–23°C, lượng mưa >1.200 mm, đất đỏ vàng phát triển trên đá mẹ granit, phiến thạch có pH 4,5–6,5. Gần đây, Yuhao Lu et al. (2015) ứng dụng mô hình sinh lý 3-PG chứng minh vùng năng suất cao (>450 $m^3$/ha ở tuổi 20) tập trung tại vùng Đông Nam Trung Quốc và có xu hướng dịch chuyển lên phía Bắc do biến đổi khí hậu.
  2. Chọn giống và cải thiện di truyền: Các chương trình khảo nghiệm xuất xứ quy mô quốc gia tại Trung Quốc từ thập niên 1970–1980 đã chọn lọc được 9 xuất xứ ưu việt (Dung Thủy, Lạc Xương, Kiến Âu, Cẩm Bính...) vượt trội 210–254% về trữ lượng so với nhóm trung bình (Hong et al., 1985; Fan & Yu, 1987). Vườn giống vô tính thế hệ thứ 3 giúp nâng cao sản lượng gỗ từ 33–72% so với quần thể ban đầu.
  3. Kỹ thuật thâm canh và tỉa thưa nuôi dưỡng: Tranh luận học thuật lớn tồn tại xung quanh phương thức trồng và mật độ: các tác giả như Wu et al. (2020) và Hong et al. (1985) ủng hộ trồng hỗn giao Sa mộc với cây lá rộng (Giổi bắc, Sát mộc) để cải thiện chu trình khoáng N, P, K và hạn chế thoái hóa đất; trong khi đó CCG of Chinese Fir (1981a) và Fang (1987) tập trung vào tối ưu hóa thâm canh thuần loài với mật độ 1.500–2.500 cây/ha kết hợp bón phân NPK. Về tỉa cành, Li et al. (2020) cho rằng tỉa cành ít ảnh hưởng đến quang hợp và sinh trưởng thân cây, trong khi Fung (1993) nhấn mạnh vai trò nâng cao phẩm cấp gỗ không mắt tật. Về tỉa thưa, Li R. (2020) chỉ ra tỉa thưa mạnh ở tuổi 6–10 thúc đẩy tăng trưởng sinh khối cây cá thể gấp 1,63–2,0 lần nhưng làm giảm tổng sinh khối tức thời của lâm phần.

Tại Việt Nam, các công trình trước đây của Lâm Phúc Cố (1996), Hà Văn Huy (1995), Lê Đình Khả et al. (2006), Đoàn Tiến Vinh et al. (2009), Nguyễn Bá Triệu (2016) và Phạm Minh Toại & Bùi Thế Đồi (2016) chủ yếu dừng lại ở việc khảo sát sinh trưởng sơ kỳ, thử nghiệm nông lâm kết hợp hoặc thăm dò xuất xứ vùng Tây Bắc. Nghiên cứu của Nguyễn Hữu Thiện (2011) đã đề xuất hàm sinh trưởng Weibull và mô hình chặt chuyển hóa tầng dưới cho rừng Sa mộc, song chưa kiểm chứng thực nghiệm đa nhân tố lập địa và các tác động thâm canh đồng bộ tại vùng Đông Bắc Bộ. Luận án của Lê Thị Ngọc Hà định vị chính xác vào khoảng trống này, kết nối chặt chẽ giữa các thành tựu lâm sinh quốc tế (Trung Quốc, Brazil, New Zealand) với điều kiện lập địa suy thoái vùng đồi núi Đông Bắc Việt Nam, đưa ra các luận cứ định lượng có giá trị dự báo cao.

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

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

Luận án đóng góp quan trọng vào việc kiểm chứng và mở rộng Lý thuyết Động thái lâm phần (Stand Dynamics Theory) và Quy luật tự tỉa thưa (Self-thinning Rule - định luật $y = k \cdot x^{-3/2}$ của Reineke/Yoda) trên đối tượng cây lá kim cận nhiệt đới di thực vào vùng nhiệt đới gió mùa. Nghiên cứu chứng minh rằng:

  • Phân hóa xã hội học cây rừng trong lâm phần Sa mộc diễn ra khốc liệt từ giai đoạn 6–7 năm tuổi khi độ tàn che đạt $>0,7$. Sự cạnh tranh ánh sáng và không gian dinh dưỡng khiến các cây yếu thế (cấp IV, V theo phân loại Kraft) bị suy giảm sinh trưởng nghiêm trọng, xác nhận mô hình phân hóa cấu trúc thân cây chịu sự chi phối mạnh mẽ của mật độ ban đầu.
  • Tác giả mở rộng Lý thuyết Chỉ số lập địa (Site Index Theory) bằng việc chứng minh sinh trưởng thể tích thân cây Sa mộc không chỉ phụ thuộc vào các biến vĩ mô (nhiệt độ, lượng mưa) mà bị kiểm soát trực tiếp bởi các chỉ số lý hóa đất tầng mặt ($0-30$ cm), đặc biệt là dung trọng ($D_{trong}$), dung tích hấp thu cation (CEC), hàm lượng mùn (OM) và tỷ lệ cấp hạt sét (Clay%).

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết lâm sinh học hiện đại:

  1. Lý thuyết Ổ sinh thái lập địa (Ecological Niche & Site-matching Framework): Phân tích tương quan đa biến giữa 11 nhân tố môi trường (địa hình, khí hậu, tính chất đất) với năng suất tích lũy sinh khối.
  2. Lý thuyết Tác động lâm sinh can thiệp (Silvicultural Disturbance & Stand Response): Đánh giá phản ứng sinh học của cây cá thể và toàn lâm phần đối với các biện pháp kỹ thuật: làm đất (kích thước hố), dinh dưỡng bổ sung (phân vô cơ NPK và vi sinh MF1), và giải tỏa áp lực cạnh tranh cơ giới (tỉa cành, tỉa thưa).
  3. Mô hình Allometric Động thái cấu trúc rừng: Ứng dụng các hàm toán học mô phỏng phân bố tần số đường kính $N/D_{1.3}$, tương quan $H_{vn} - D_{1.3}$, và tương quan tăng trưởng trữ lượng $\Delta M$ theo thời gian hồi phục $20, 32, 42$ tháng sau can thiệp.
       [ĐIỀU KIỆN LẬP ĐỊA ĐÔNG BẮC BỘ]
   (Địa hình, Khí hậu, Thổ nhưỡng, Độ phì đất)
  [CAN THIỆP KỸ THUẬT GÂY TRỒNG BAN ĐẦU]
Làm đất          Tuổi cây con         Mật độ trồng
(40x40x40 cm)   (12-14 tháng)    (1.600-2.000 cây/ha)
   [ĐIỀU TIẾT DINH DƯỠNG & KHÔNG GIAN SỐNG]
Bón lót & Thúc                    Tỉa cành định hình
(NPK + MF1 vi sinh)             (1/3 chiều cao dưới cành)
  [CHUYỂN HÓA RỪNG GỖ LỚN BẰNG TỈA THƯA]
Tỉa thưa Tuổi 7                  Tỉa thưa Tuổi 11
(Mật độ 1.000-1.200 cây/ha)     (Mật độ 800-1.000 cây/ha)
  [TỐI ƯU HÓA SINH TRƯỞNG & GIÁ TRỊ KINH TẾ]
 (Tăng trưởng D1.3, Hvn, Trữ lượng M, Tỷ lệ gỗ lớn >80%)

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

Thiết kế nghiên cứu

Luận án áp dụng thế giới quan thực chứng (Positivism) với thiết kế thực nghiệm định lượng nghiêm ngặt kết hợp giữa phương pháp điều tra hồi cứu lâm phần hiện có và bố trí thí nghiệm nhân tố đầy đủ (Factorial Randomized Complete Block Design - RCBD) trên hiện trường:

  • Hệ thống điều tra lập địa: Bố trí 24 ô tiêu chuẩn định hình ($500-1.000\ m^2$/OTC) tại 8 xã thuộc 6 huyện của 3 tỉnh Cao Bằng, Lạng Sơn, Quảng Ninh nhằm đánh giá thực trạng cấu trúc, trữ lượng và các chỉ số lập địa.
  • Hệ thống thí nghiệm kỹ thuật lâm sinh: Thiết lập tại huyện Ba Chẽ, tỉnh Quảng Ninh trên diện tích đồng nhất với 5 khối thí nghiệm chính:
    1. Thí nghiệm làm đất: 4 công thức kích thước hố ($30\times30\times40$ cm, $40\times40\times40$ cm, $50\times50\times40$ cm, $60\times60\times40$ cm), lặp lại 3 lần.
    2. Thí nghiệm tuổi cây con: 4 cấp tuổi ($9, 12, 15, 18$ tháng tuổi).
    3. Thí nghiệm mật độ trồng ban đầu: 5 cấp mật độ ($1.100, 1.600, 2.000, 2.500, 3.300$ cây/ha).
    4. Thí nghiệm bón phân: 6 công thức phân bón kết hợp phân khoáng NPK và phân hữu cơ vi sinh MF1.
    5. Thí nghiệm tỉa cành: 4 cường độ tỉa (đối chứng không tỉa, tỉa 1/3, tỉa 1/2 và tỉa 2/3 chiều cao tán).
  • Hệ thống thí nghiệm tỉa thưa và bón phân nuôi dưỡng: Thực hiện trên rừng 7 tuổi và 11 tuổi tại Ba Chẽ, theo dõi định kỳ sau $20, 32$ và $42$ tháng:
    • Rừng 7 tuổi: Các công thức mật độ để lại ($800, 1.000, 1.200$ cây/ha và đối chứng $1.600-2.000$ cây/ha) kết hợp thí nghiệm bón bổ sung phân NPK + vi sinh.
    • Rừng 11 tuổi: Các công thức mật độ để lại ($600, 800, 1.000$ cây/ha và đối chứng) kết hợp thí nghiệm bón phân sau tỉa.

Quy trình nghiên cứu và chuẩn mực thu thập số liệu

Quy trình điều tra trắc lượng rừng và phân tích mẫu đất tuân thủ tiêu chuẩn quốc gia và quốc tế:

  • Đo đếm lâm học: Đường kính ngang ngực ($D_{1.3}$) đo bằng thước kẹp kính chuyên dụng chính xác đến $0,1$ cm; đường kính gốc ($D_0$) đo tại vị trí sát cổ rễ; chiều cao vút ngọn ($H_{vn}$) và chiều cao dưới cành ($H_{dc}$) đo bằng máy đo cao Vertex IV chính xác đến $0,1$ m; đường kính tán ($D_t$) đo theo 2 hướng vuông góc. Thể tích thân cây ($V$) và trữ lượng lâm phần ($M$) tính toán thông qua hình số thân cây ($f$) thực nghiệm: $V = \frac{\pi}{4} \cdot D_{1.3}^2 \cdot H_{vn} \cdot f$.
  • Phân tích thổ nhưỡng: Lấy mẫu đất tại 2 tầng sâu ($0-30$ cm và $30-60$ cm) theo phương pháp đường chéo góc tại các OTC. Phân tích các chỉ tiêu: dung trọng ($D_{trong}$ - phương pháp ống dung tích chuẩn), thành phần cơ giới (tỷ lệ Sand, Silt, Clay - phương pháp ống hút Robinson), độ chua ($pH_{KCl}$ - đo điện thế), hàm lượng mùn tổng số (OM% - phương pháp Walkley-Black), đạm tổng số ($N%$ - Kjeldahl), lân dễ tiêu ($P_2O_5$ - phương pháp Oniani), kali trao đổi ($K_2O$ - quang kế ngọn lửa) và dung tích trao đổi cation (CEC - phương pháp Amoni axetat).

Data và phân tích thống kê

Toàn bộ cơ sở dữ liệu thực nghiệm được xử lý bằng các phần mềm thống kê chuyên sâu (SAS 9.4, SPSS 22.0 và R):

  • Phân tích phương sai một nhân tố và nhiều nhân tố (ANOVA) để kiểm định sự sai khác sinh trưởng giữa các công thức thí nghiệm.
  • Kiểm nghiệm sau phân tích phương sai bằng tiêu chuẩn Duncan (Duncan's Multiple Range Test) ở mức ý nghĩa $\alpha = 0,05$ và $\alpha = 0,01$.
  • Phân tích tương quan tuyến tính và phi tuyến từng bước (Stepwise Multiple Regression) để xây dựng mô hình liên hệ giữa tăng trưởng thể tích ($\Delta M$) với các nhân tố lập địa, đánh giá độ tin cậy thông qua hệ số tương quan ($R$), hệ số xác định ($R^2$), chỉ số kiểm định Fisher ($F$) và giá trị $p < 0,001$.

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: Kết quả phân tích hồi quy đa biến chứng minh tăng trưởng trữ lượng ($\Delta M$) của Sa mộc tại Đông Bắc Bộ chịu ảnh hưởng quyết định bởi 4 nhân tố cốt lõi: độ dày tầng đất ($H_{dat}$), dung trọng ($D_{trong}$), hàm lượng mùn tầng mặt ($OM%$) và độ dốc ($\alpha^\circ$). Phương trình tương quan đa biến đạt mức ý nghĩa thống kê rất cao ($R^2 > 0,85; p < 0,001$). Sa mộc sinh trưởng tối ưu tại lập địa có độ cao $300-800$ m, độ dốc $<20^\circ$, tầng đất dày $>60$ cm, thành phần cơ giới thịt nhẹ đến sét pha cát, dung trọng đất $<1,2\ g/cm^3$, hàm lượng mùn $>2,5%$.
  2. Quy chuẩn kích thước làm đất và tiêu chuẩn cây giống: Kỹ thuật đào hố $40 \times 40 \times 40$ cm và $50 \times 50 \times 40$ cm mang lại sinh trưởng đường kính ($D_0$) và chiều cao ($H_{vn}$) vượt trội có ý nghĩa thống kê ($p < 0,05$) so với hố $30 \times 30 \times 30$ cm truyền thống; việc tăng kích thước hố lên $60 \times 60 \times 40$ cm không làm tăng sinh trưởng rõ rệt nhưng làm tăng 35% chi phí nhân công. Cây con 12–15 tháng tuổi đạt tỷ lệ sống $>95%$, chiều cao xuất vườn $>35$ cm, đường kính cổ rễ $>0,45$ cm có sức đề kháng và tốc độ bứt ngọn vượt trội so với cây con 9 tháng tuổi ($p < 0,01$).
  3. Hiệu ứng mật độ trồng và dinh dưỡng phân bón: Mật độ ban đầu thích hợp nhất cho mục tiêu kinh doanh gỗ lớn kết hợp tỉa thưa là $1.600-2.000$ cây/ha. Mật độ $3.300$ cây/ha làm giảm $28,4%$ tăng trưởng đường kính cá thể sau 3 năm do cạnh tranh rễ sớm. Bón lót $200\ g\ NPK\ (5:10:3) + 200\ g$ phân vi sinh MF1/hố kết hợp bón thúc năm thứ 2 giúp tăng trưởng thể tích thân cây cao hơn $22,6-31,2%$ so với công thức đối chứng không bón ($p < 0,01$).
  4. Phát hiện nghịch lý về tỉa cành (Counter-intuitive Finding): Tỉa cành ở cường độ mạnh (1/2 đến 2/3 chiều cao tán) làm giảm nhẹ sinh trưởng đường kính ngang ngực trong 12 tháng đầu do mất diện tích quang hợp hữu hiệu. Tỉa cành hợp lý nhất là tỉa các cành chết, cành sâu bệnh và cành sát gốc không quá 1/3 chiều cao thân cây ở giai đoạn rừng bắt đầu giao tán (tuổi 4–5), đảm bảo thân cây tròn đều, giảm mắt tật mà không làm suy giảm sinh trưởng tích lũy.
  5. Đột phá về chế độ tỉa thưa chuyển hóa rừng 7 tuổi và 11 tuổi:
    • Đối với rừng 7 tuổi: Tỉa thưa tầng dưới đưa mật độ từ $1.800$ cây/ha về $1.000-1.200$ cây/ha giúp tăng trưởng đường kính ngang ngực ($\Delta D_{1.3}$) sau 42 tháng đạt $1,42-1,58$ cm/năm, cao hơn $45,8%$ so với đối chứng không tỉa ($0,97$ cm/năm). Bón thúc $0,5\ kg\ NPK + 0,5\ kg$ phân vi sinh/cây sau tỉa thưa làm tăng thêm $18,4%$ trữ lượng gỗ lớn so với không bón.
    • Đối với rừng 11 tuổi: Tỉa thưa đưa mật độ về $800-1.000$ cây/ha kết hợp bón phân thúc đẩy đường kính lâm phần đạt bình quân $22,8$ cm ở tuổi 15, nâng tỷ lệ gỗ lớn thương phẩm ($D_{1.3} \ge 20$ cm) lên $>82%$.

Implications đa chiều

  • Về mặt học thuật: Cung cấp bộ dữ liệu thực nghiệm chuẩn xác về phản ứng động thái học của Sa mộc di thực, đóng góp cơ sở dữ liệu quan trọng cho các mô hình sinh thái học rừng nhiệt đới và cận nhiệt đới.
  • Về phương pháp luận: Chuẩn hóa quy trình đánh giá thực nghiệm lâm sinh kết hợp chuỗi thời gian ($20, 32, 42$ tháng) với phân tích lý hóa đất đa biến, có thể nhân rộng nghiên cứu cho các loài cây gỗ lớn bản địa khác (Lim xanh, Giổi, Re, Vối thuốc).
  • Về thực tiễn sản xuất: Cung cấp "gói kỹ thuật hoàn chỉnh" cho các công ty lâm nghiệp và chủ rừng tại Cao Bằng, Lạng Sơn, Quảng Ninh: chọn đúng lập địa $\rightarrow$ cây giống $12-14$ tháng $\rightarrow$ hố $40 \times 40 \times 40$ cm $\rightarrow$ mật độ $1.600-2.000$ cây/ha $\rightarrow$ bón lót NPK + MF1 $\rightarrow$ tỉa thưa lần 1 ở tuổi 7 (còn $1.000-1.200$ cây/ha) $\rightarrow$ tỉa thưa lần 2 ở tuổi 11 (còn $800$ cây/ha) $\rightarrow$ thu hoạch gỗ lớn ở tuổi 15–20.
  • Về chính sách: Cung cấp luận cứ kỹ thuật trực tiếp để Bộ Nông nghiệp và Phát triển Nông thôn cập nhật Tiêu chuẩn Quốc gia (TCVN) về trồng rừng thâm canh Sa mộc gỗ lớn, thay thế quy trình kỹ thuật cũ ban hành từ năm 2002.

Limitations và Future Research

  1. Giới hạn không gian thực nghiệm tập trung: Các thí nghiệm chuyên sâu về kỹ thuật lâm sinh (làm đất, phân bón, tỉa cành, tỉa thưa) được bố trí chủ yếu tại huyện Ba Chẽ, tỉnh Quảng Ninh. Mặc dù Quảng Ninh mang tính đại diện cao cho vùng Đông Bắc, các biến dị tiểu khí hậu tại vùng núi cao Cao Bằng và Lạng Sơn có thể tạo ra những sai lệch nhỏ về biên độ phản ứng sinh trưởng.
  2. Giới hạn thời gian quan trắc sau can thiệp: Thời gian theo dõi thực nghiệm sau tỉa thưa đạt mốc 42 tháng (3,5 năm). Để đánh giá trọn vẹn toàn bộ luân kỳ kinh doanh gỗ lớn (15–20 năm), cần tiếp tục duy trì quan trắc đến kỳ khai thác chính.
  3. Chưa đánh giá cơ lý tính gỗ sau tỉa thưa: Nghiên cứu tập trung vào năng suất sinh khối và tăng trưởng kích thước hình học ($D_{1.3}, H_{vn}, V$), chưa tiến hành kiểm tra thực nghiệm môđun đàn hồi (MOE), môđun đứt gãy (MOR) và mật độ khối lượng thể tích gỗ dưới tác động của các chế độ bón phân và tỉa thưa khác nhau.

Chương trình nghiên cứu tiếp nối (Future Research Agenda):

  • Mở rộng mạng lưới thử nghiệm silviculture đa địa điểm (Multi-site validation) tại các đai cao $>800$ m thuộc Cao Bằng, Hà Giang và Lào Cai.
  • Nghiên cứu công nghệ chọn tạo giống phân tử kết hợp vi nhân giống vô tính (Tissue culture / Macro-cutting propagation) từ các cây mẹ ưu trội (Plus trees) đã xác định.
  • Ứng dụng mô hình sinh thái quá trình (Process-based models như 3-PG hoặc CABALA) tích hợp dữ liệu viễn thám LiDAR và vi khí hậu để dự báo năng suất rừng Sa mộc dưới các kịch bản biến đổi khí hậu (IPCC RCP 4.5 & 8.5).
  • Nghiên cứu động thái chu trình dinh dưỡng đất, hoạt tính enzyme và đa dạng vi sinh vật vùng rễ (Rhizosphere microbiome) trong rừng trồng Sa mộc luân kỳ 2 và luân kỳ 3 nhằm ngăn ngừa hiện tượng suy thoái đất lâm nghiệp.

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

  • Tác động học thuật (Academic Impact): Công trình tạo lập hệ thống dữ liệu lâm sinh học chuẩn xác về Sa mộc tại Việt Nam, đóng vai trò là tài liệu tham khảo cốt lõi cho các nghiên cứu sau đại học, dự kiến thu hút trích dẫn khoa học cao trong các tạp chí chuyên ngành Lâm nghiệp và Nông nghiệp nhiệt đới.
  • Chuyển đổi kinh tế - ngành gỗ (Industry Transformation): Việc chuyển đổi từ kinh doanh gỗ nhỏ sang rừng gỗ lớn thâm canh nâng cao giá trị gia tăng trên 1 đơn vị diện tích đất lâm nghiệp từ 20,1 triệu đồng/ha lên 290,7 triệu đồng/ha, tạo bước ngoặt cung ứng nguyên liệu gỗ xẻ, gỗ xây dựng cao cấp nội địa, giảm thiểu áp lực nhập khẩu gỗ tròn từ nước ngoài.
  • Tác động chính sách lâm nghiệp (Policy Influence): Cung cấp bằng chứng thực nghiệm vững chắc phục vụ công tác quy hoạch vùng trồng rừng sản xuất gỗ lớn tại 14 tỉnh trung du miền núi phía Bắc, góp phần hoàn thành chỉ tiêu 360.000 ha rừng gỗ lớn quốc gia theo Quyết định 1565/QĐ-BNN-TCLN.
  • Lợi ích xã hội và môi trường (Societal & Environmental Benefits): Giúp đồng bào các dân tộc thiểu số tại vùng cao Đông Bắc chuyển dịch phương thức canh tác từ quảng canh, khai thác non sang quản lý rừng bền vững chu kỳ dài. Rừng Sa mộc thâm canh hấp thụ từ $5,58-6,69$ tấn $CO_2$/ha/năm ở giai đoạn $7-11$ tuổi (Chen et al., 2017), tăng cường khả năng giữ nước, chống xói mòn đất vùng đầu nguồn xung yếu.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận nghiên cứu thực nghiệm lâm sinh chuẩn mực, hệ thống mô hình toán hồi quy lập địa - sinh trưởng và các khoảng trống lý thuyết cần tiếp tục khai phá.
  • Các Viện nghiên cứu và Trường Đại học Lâm nghiệp: Bộ giáo trình và tài liệu tham khảo chuyên sâu cho các chuyên ngành Lâm sinh, Quản lý tài nguyên rừng, Đo đếm và Điều tra quy hoạch rừng.
  • Các Công ty Lâm nghiệp, Ban Quản lý rừng và Doanh nghiệp chế biến gỗ: Sở hữu quy trình kỹ thuật định lượng chuẩn xác để áp dụng vào thực tế sản xuất thâm canh rừng gỗ lớn, kiểm soát chi phí đầu tư và dự báo chính xác sản lượng thu hoạch.
  • Người trồng rừng và Chủ trang trại vùng Đông Bắc: Nâng cao thu nhập bền vững trên một chu kỳ kinh doanh (tăng gấp 8,5 đến 14,5 lần giá trị kinh tế), tối ưu hóa việc sử dụng cây giống và phân bón.
  • Cơ quan quản lý Nhà nước (Cục Lâm nghiệp, Sở NN&PTNT các tỉnh): Công cụ đắc lực để ban hành định mức kinh tế kỹ thuật, phê duyệt các dự án trồng rừng gỗ lớn và xây dựng chính sách chi trả dịch vụ môi trường rừng (hấp thụ Carbon).

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: Luận án đã mở rộng Lý thuyết Ổ sinh thái lập địa (Site Index Theory) và Quy luật Cạnh tranh lâm phần (Stand Competition Theory) trong lâm sinh học nhiệt đới bằng việc xác lập mô hình tương quan phi tuyến đa biến giữa tổ hợp nhân tố lý hóa đất tầng mặt ($0-30$ cm: dung trọng, mùn, sét) kết hợp độ dốc địa hình với tốc độ tăng trưởng thể tích thân cây ($\Delta M$), chứng minh vai trò điều hòa quyết định của độ xốp và lý tính đất đối với sự phát triển hệ rễ cọc nông của Sa mộc.

2. Phương pháp nghiên cứu có tính đột phá gì khi so sánh với các công trình quốc tế? Trả lời: So với các nghiên cứu tại Trung Quốc (Fang, 1987; Li et al., 2020) vốn thường chỉ đánh giá đơn biến hoặc tập trung vào khảo nghiệm giống vô tính, luận án đã thiết kế hệ thống thực nghiệm nhân tố đầy đủ (Factorial RCBD) kết hợp theo dõi động thái đa thời gian ($20, 32, 42$ tháng) đồng thời trên 5 khâu kỹ thuật lâm sinh (làm đất, tuổi cây, mật độ, phân bón, tỉa cành) và 2 cấp tuổi tỉa thưa (tuổi 7 và tuổi 11) trên cùng một lập địa đồng nhất, tạo ra tính so sánh và kiểm soát sai số thực nghiệm cao nhất.

3. Phát hiện bất ngờ nhất từ số liệu thực nghiệm là gì? Trả lời: Phát hiện về phản ứng của cây Sa mộc đối với kỹ thuật tỉa cành: việc tỉa cành nhân tạo ở cường độ cao ($>1/3$ chiều cao tán) không làm tăng trưởng đường kính như kỳ vọng lâm sinh truyền thống mà ngược lại làm sụt giảm sinh trưởng thể tích thân cây trong ngắn hạn do đặc tính tự tỉa cành tự nhiên của loài rất tốt. Điều này giúp tối ưu hóa chi phí sản xuất bằng cách loại bỏ việc can thiệp tỉa cành sâu không cần thiết.

4. Quy trình nhân rộng và tái lặp nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Trả lời: Luận án cung cấp chi tiết toàn bộ thông số kỹ thuật có thể tái lặp 100%: tọa độ và đặc điểm lập địa 24 OTC điều tra, kích thước hố đào ($40\times40\times40$ cm), tiêu chuẩn cây con xuất vườn ($12-14$ tháng, cao $>35$ cm, $D_0 > 0,45$ cm), công thức phối trộn phân bón NPK và vi sinh MF1, phương pháp phân loại cây tỉa thưa theo 5 cấp Kraft và thuật toán xử lý phương sai ANOVA/Duncan trên SAS/R.

5. Luận án đã vạch ra chương trình nghiên cứu 10 năm tiếp theo như thế nào? Trả lời: Chương trình 10 năm tập trung vào: (i) Tiếp tục theo dõi các ô thí nghiệm tỉa thưa tại Ba Chẽ đến tuổi thành thục công nghệ 18–20 năm; (ii) Triển khai mô hình giải trình tự gen kết hợp nhân giống mô tế bào từ các cây Sa mộc Vua; (iii) Đánh giá biến tính cơ lý gỗ lớn phục vụ công nghiệp chế biến đồ gỗ xuất khẩu cao cấp; (iv) Xây dựng bản đồ lập địa số hóa ứng dụng GIS và thuật toán Machine Learning cho toàn bộ vùng Tây Bắc và Đông Bắc Bộ.

Kết luận

  1. Luận án đã xây dựng thành công cơ sở khoa học hoàn chỉnh và có hệ thống về lập địa và các biện pháp kỹ thuật trồng rừng thâm canh Sa mộc (Cunninghamia lanceolata) tại vùng Đông Bắc Bộ.
  2. Thiết lập mô hình hồi quy đa biến có giá trị dự báo định lượng cao về quan hệ giữa năng suất rừng với các nhân tố lập địa, xác định lập địa tối ưu cho Sa mộc: độ cao $300-800$ m, độ dốc $<20^\circ$, đất sâu $>60$ cm, tơi xốp, giàu mùn và thoát nước tốt.
  3. Chuẩn hóa trọn gói hệ thống kỹ thuật lâm sinh trồng thâm canh gỗ lớn: làm đất hố $40 \times 40 \times 40$ cm; cây con đem trồng đạt $12-14$ tháng tuổi; mật độ ban đầu $1.600-2.000$ cây/ha; bón lót $200\ g\ NPK + 200\ g$ vi sinh MF1; tỉa cành kỹ thuật dưới 1/3 thân cây ở tuổi $4-5$.
  4. Xác lập chế độ tỉa thưa nuôi dưỡng tối ưu hai giai đoạn: tỉa thưa lần 1 ở tuổi 7 (giữ lại $1.000-1.200$ cây/ha), tỉa thưa lần 2 ở tuổi 11 (giữ lại $800-1.000$ cây/ha) kết hợp bón phân bổ sung, giúp nâng tỷ lệ gỗ lớn lên $>80%$ và rút ngắn luân kỳ khai thác.
  5. Chứng minh hiệu quả kinh tế đột phá với giá trị thương phẩm đạt 290,7 triệu đồng/ha ở tuổi 20, cao gấp 14,5 lần so với rừng gỗ nhỏ chu kỳ 9 năm.
  6. Mở ra 3 hướng nghiên cứu lâm sinh học hiện đại mới: chọn giống phân tử kháng biến đổi khí hậu, sinh thái học chu trình dinh dưỡng đất rừng thâm canh luân kỳ dài, và ứng dụng viễn thám lượng hóa hấp thụ Carbon rừng trồng bản địa.