Bối cảnh và vấn đề nghiên cứu

Trong chiến lược phát triển lâm nghiệp và chương trình trồng mới 5 triệu héc-ta rừng (Dự án 661) của Việt Nam, việc lựa chọn loài cây trồng phù hợp với điều kiện lập địa nhằm nâng cao độ che phủ đất trống đồi núi trọc và cung cấp nguyên liệu cho công nghiệp chế biến gỗ là yêu cầu cấp bách. Cây Thông đuôi ngựa (Pinus massoniana Lamb, thuộc họ Thông - Pinaceae) là loài cây gỗ lá kim bản địa có khả năng sinh trưởng trên các sườn dốc khô cằn, đất nghèo dinh dưỡng nhiều đá lẫn, đóng vai trò là cây trồng rừng kinh tế chủ lực tại các tỉnh vùng Đông Bắc như Lạng Sơn, Quảng Ninh, Bắc Giang. Gỗ cây được sử dụng làm trụ mỏ, cột điện, bột giấy, đồ gia dụng, trong khi nhựa thông là nguồn nguyên liệu xuất khẩu có giá trị thương phẩm cao.

Tại xã Gia Cát, huyện Cao Lộc, tỉnh Lạng Sơn, sau hơn một thập kỷ triển khai trồng rừng, diện tích Thông đuôi ngựa thuần loài đã mở rộng đáng kể. Tuy nhiên, thực tế sản xuất cho thấy trữ lượng và chất lượng rừng tại đây còn thấp, chưa phát huy hết tiềm năng kinh tế do công tác quản lý và các biện pháp kỹ thuật lâm sinh tác động (đặc biệt là mật độ trồng và chế độ chăm sóc, tỉa thưa) chưa được áp dụng hợp lý. Nhằm xác định cơ sở khoa học về sinh thái và kinh tế cho loài cây này tại địa phương, đề tài khóa luận tập trung giải quyết các mục tiêu và nhiệm vụ cụ thể:

  1. Đánh giá tình hình sinh trưởng về đường kính ngang ngực, chiều cao vút ngọn, chiều cao dưới cành, đường kính tán và trữ lượng của hai mô hình rừng trồng Thông đuôi ngựa thuần loài ở các mật độ khác nhau.
  2. Đánh giá chất lượng rừng và xác định quy luật kết cấu lâm phần (phân bố số cây theo đường kính, chiều cao và quy luật tương quan chiều cao - đường kính).
  3. Đánh giá hiệu quả kinh tế của các mô hình thông qua hệ thống chỉ tiêu tài chính (NPV, BCR, IRR) trong chu kỳ kinh doanh 25 năm.
  4. Đề xuất một số giải pháp kỹ thuật lâm sinh thích hợp nhằm nâng cao năng suất, chất lượng và hiệu quả kinh tế của rừng trồng Thông đuôi ngựa.

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

  • Đối tượng: Hai mô hình rừng trồng Thông đuôi ngựa (Pinus massoniana Lamb) thuần loài 14 tuổi (trồng tháng 4 năm 1998) do Ủy ban nhân dân xã Gia Cát quản lý, gồm Mô hình 1 có mật độ hiện tại 900 cây/ha và Mô hình 2 có mật độ hiện tại 1060 cây/ha.
  • Không gian nghiên cứu: Địa bàn xã Gia Cát, huyện Cao Lộc, tỉnh Lạng Sơn.
  • Thời gian nghiên cứu: Thu thập số liệu thực địa vào năm 2012 (thời điểm rừng đạt 14 tuổi); tính toán hiệu quả kinh tế xác định theo chu kỳ kinh doanh 25 năm (dự kiến khai thác chính vào năm 2023 - 2024).

Cơ sở lý thuyết và phương pháp

Khóa luận kế thừa hệ thống lý thuyết về điều tra quy luật sinh trưởng lâm phần và phương pháp luận đánh giá hiệu quả kinh tế đầu tư lâm nghiệp của FAO (Hans M. Gregersen và Arnoldo H. Contreras) cùng các giáo trình chuyên ngành trong nước.

Hệ thống phương pháp và công cụ nghiên cứu cụ thể:

  • Phương pháp ngoại nghiệp:

    • Bố trí ô tiêu chuẩn (OTC): Áp dụng phương pháp so sánh loại trừ trên các vị trí đại diện (sườn dốc, chân và đỉnh đồi). Thiết lập 6 OTC điển hình tạm thời (mỗi mô hình 3 OTC), diện tích mỗi OTC là 500 m² (kích thước $25 \times 20\text{ m}$), hình chữ nhật lập theo phương pháp Pitago, chiều dài song song và chiều rộng vuông góc với đường đồng mức; sai số khép góc $< 1/200$ chu vi.
    • Thu thập chỉ tiêu sinh trưởng cây đứng: Đo đường kính ngang ngực ($D_{1.3}$) bằng thước kẹp kính chính xác đến mm theo hai hướng Đông - Tây và Nam - Bắc; đo chiều cao vút ngọn ($H_{vn}$) và chiều cao dưới cành ($H_{dc}$) bằng thước Blumeleiss; đo đường kính tán ($D_t$) bằng thước dây qua hình chiếu tán thẳng đứng.
    • Phân cấp chất lượng cây đứng: Chia thành 3 cấp phẩm chất gồm cây tốt (cấp A), cây trung bình (cấp B) và cây xấu (cấp C) dựa trên độ thẳng thân cây, mức độ phát triển tán và tình trạng sâu bệnh.
    • Điều tra thực bì và thảm mục: Mỗi OTC bố trí 5 ô dạng bản (ODB) diện tích $4\text{ m}^2$ ($2 \times 2\text{ m}$) tại 4 góc và 1 ô ở tâm để xác định thành phần loài, độ che phủ và phẩm chất sinh trưởng.
    • Điều tra thổ nhưỡng: Đào 2 phẫu diện đất đại diện (rộng 70 - 90 cm, dài 1,2 - 1,5 m, sâu 1,2 m hoặc hết tầng B) để phân tích hình thái tầng phát sinh và các chỉ tiêu lý tính.

  • Phương pháp nội nghiệp và thống kê sinh trắc:

    • Tính toán các đại lượng mẫu trung bình, phương sai, độ lệch chuẩn và hệ số biến động ($S%$) bằng chức năng Data Analysis - Descriptive Statistics trong Excel theo tài liệu của Nguyễn Hải Tuất (1982) và Vũ Tiến Hinh cùng cộng sự (2006).
    • Kiểm tra sự sai khác trung bình giữa các OTC và giữa hai mô hình bằng kiểm định tiêu chuẩn $U$ của phân bố chuẩn ($|U| > 1,96$ biểu thị sai khác có ý nghĩa thống kê ở mức xác suất 95%).
    • Chia tổ thực nghiệm theo công thức kinh nghiệm của Brooks và Carruthers; mô hình hóa quy luật phân bố số cây theo cỡ kính ($N/D_{1.3}$) và theo chiều cao ($N/H_{vn}$) bằng hàm phân bố xác suất liên tục Weibull: $$F(x) = \alpha \lambda x^{\alpha-1} e^{-\lambda x^\alpha}$$ Kiểm tra độ tương thích giữa tần số lý thuyết và thực nghiệm bằng tiêu chuẩn Pearson ($\chi^2$).
    • Mô phỏng tương quan thực nghiệm chiều cao - đường kính ($H_{vn}/D_{1.3}$) theo dạng hàm logarit: $$H_{vn} = a + b \log D_{1.3}$$
    • Xác định thể tích cây đứng cá thể và trữ lượng lâm phần: $$V_i = \frac{1}{4} \pi D_{1.3}^2 H_{vn} F$$ (với hình số thân cây rừng trồng $F = 0,45$). Trữ lượng lâm phần tính bằng công thức: $$M/\text{ha} = M_{OTC} \times \frac{10000}{S_{OTC}}$$
    • Đánh giá chất lượng lâm phần bằng kiểm định $\chi^2$ độc lập (hoặc Kruskal-Wallis) cho nhiều mẫu.
    • Đánh giá hiệu quả tài chính: Cân đối dòng tiền chi phí và thu nhập cho 1 ha rừng trong chu kỳ kinh doanh 25 năm với tỷ lệ chiết khấu $r = 3%$ (lãi suất vay vốn lâm nghiệp năm 1998), tính toán các chỉ tiêu: Giá trị hiện tại thuần (NPV), Tỷ số lợi ích - chi phí (BCR) và Suất thu hồi vốn nội bộ (IRR).

  • Nguồn dữ liệu: Số liệu điều tra thực địa năm 2012 tại xã Gia Cát; hồ sơ thiết kế trồng rừng năm 1998; tài liệu thống kê kinh tế - xã hội địa phương; phỏng vấn hộ gia đình và đơn giá khai thác lâm sản thực tế trên địa bàn.

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

Phần I: Tổng quan vấn đề nghiên cứu

Chương này tổng kết các công trình khoa học trong và ngoài nước về sinh trưởng cây rừng, phương pháp đánh giá hiệu quả kinh tế và đặc điểm loài Thông đuôi ngựa. Trên thế giới, quá trình sinh trưởng được hệ thống hóa từ nhận định của Leonardo da Vinci thế kỷ XV đến mô hình Gompertz (1895), Odum (1974), Beketop (1968), Romell (1932); phương pháp kinh tế dựa trên khung hướng dẫn của FAO và các công trình của Hans M. Gregersen. Tại Việt Nam, tác giả tổng hợp kết quả nghiên cứu sinh trưởng của Phùng Ngọc Lan, Nguyễn Ngọc Lung, Vũ Tiến Hinh, Lê Mộng Chân; hiệu quả kinh tế của Đặng Trung Thuận, Vương Văn Quỳnh; các nghiên cứu chuyên sâu về Thông đuôi ngựa của Vũ Nhâm, Đồng Sỹ Hiền, Võ Đại Hải, Đặng Thiện Triều. Đồng thời, chương trình bày chi tiết đặc điểm hình thái học, đặc tính sinh thái, phân bố và giá trị kinh tế của Thông đuôi ngựa trong khai thác gỗ trụ mỏ, nhựa dầu và tiềm năng tích lũy carbon.

Phần II: Mục tiêu - Đối tượng - Giới hạn - Nội dung và phương pháp nghiên cứu

Chương trình bày mục tiêu tổng quát, nhiệm vụ cụ thể, đối tượng nghiên cứu (hai mô hình rừng thông thuần loài 14 tuổi ở mật độ 900 cây/ha và 1060 cây/ha). Nội dung quy định đầy đủ các bước thực nghiệm ngoại nghiệp (lập ô tiêu chuẩn 500 m², kỹ thuật đo đếm chỉ tiêu lâm học, điều tra thực bì qua 5 ô dạng bản, đào phẫu diện đất) và các thuật toán nội nghiệp (tính sai số, kiểm định $U$, hàm phân bố Weibull, tương quan logarit, hình số $F=0,45$, công thức tính NPV, BCR, IRR).

Phần III: Điều kiện cơ bản của khu vực nghiên cứu

Chương mô tả đặc điểm tự nhiên và kinh tế - xã hội xã Gia Cát, huyện Cao Lộc. Xã có diện tích tự nhiên 3675 ha, địa hình đồi núi dốc chia cắt mạnh, độ cao trung bình trên 300 m. Khí hậu á nhiệt đới mát mẻ, nhiệt độ trung bình năm 21°C, lượng mưa trung bình năm 1392 mm (tập trung 70% từ tháng 5 đến tháng 9), độ ẩm 82%, có sương muối vào mùa đông. Diện tích đất có rừng là 1613,38 ha (chiếm 43,90% diện tích tự nhiên). Toàn xã có 914 hộ, 4554 nhân khẩu, nhân khẩu nông nghiệp chiếm 95,74%, thu nhập bình quân 2,3 triệu đồng/người/năm. Rừng nghiên cứu được trồng vào tháng 4/1998 với mật độ ban đầu 1800 cây/ha ($1,8 \times 3,0\text{ m}$), sử dụng cây con rễ trần 6 tháng tuổi, cuốc hố thủ công $40 \times 40 \times 40\text{ cm}$, bón lót 0,3 kg NPK/hố; chăm sóc 3 năm đầu, bón thúc năm 1; phòng trừ sâu róm thông bằng thuốc Vophatox, Thiophat và bảo tồn thiên địch (ong mắt đỏ, ong trứng đen).

Phần IV: Kết quả và thảo luận

Chương tập trung phân tích thực nghiệm các chỉ tiêu sinh trưởng, cấu trúc, thổ nhưỡng và hiệu quả tài chính:

  • Sinh trưởng đường kính ngang ngực ($D_{1.3}$): Mô hình 1 đạt bình quân 17,64 cm, Mô hình 2 đạt 16,78 cm. Kiểm định $U = 2,23 > 1,96$, chứng minh mật độ 900 cây/ha cho sinh trưởng đường kính cá thể vượt trội (chênh lệch 0,86 cm) do không gian dinh dưỡng lớn hơn.
  • Sinh trưởng chiều cao ($H_{vn}, H_{dc}$): Chiều cao vút ngọn bình quân ở Mô hình 2 (13,01 m) cao hơn Mô hình 1 (12,30 m) với $U = 3,284 > 1,96$. Chiều cao dưới cành của Mô hình 2 đạt 9,45 m, cao hơn Mô hình 1 (8,69 m) với $U = 3,587 > 1,96$. Mật độ dày đã thúc đẩy quá trình tỉa cành tự nhiên và cạnh tranh ánh sáng vươn cao.
  • Đường kính tán ($D_t$): Mô hình 1 đạt 3,90 m, lớn hơn đáng kể so với Mô hình 2 (3,68 m) với kiểm định $U = 3,52 > 1,96$.
  • Trữ lượng lâm phần: Tại thời điểm 14 tuổi, thể tích một cây ở Mô hình 1 ($0,1444\text{ m}^3$) lớn hơn Mô hình 2 ($0,1419\text{ m}^3$), nhưng do số lượng cây đứng nhiều hơn nên trữ lượng Mô hình 2 đạt $151,41\text{ m}^3/\text{ha}$ (tăng trưởng bình quân $10,82\text{ m}^3/\text{ha/năm}$), vượt Mô hình 1 với $129,96\text{ m}^3/\text{ha}$ (tăng trưởng bình quân $9,28\text{ m}^3/\text{ha/năm}$). Dự tính hết chu kỳ 25 năm, trữ lượng Mô hình 2 đạt $270,75\text{ m}^3/\text{ha}$, Mô hình 1 đạt $232,00\text{ m}^3/\text{ha}$.
  • Chất lượng rừng: Mô hình 1 có tỷ lệ cây tốt 65,19%, trung bình 20,00%, xấu 14,81%. Mô hình 2 có tỷ lệ cây tốt 61,25%, trung bình 23,75%, xấu 15,00%. Kiểm định $\chi^2 = 1,57 < \chi^2_{0.05; 2} = 5,99$, cho thấy sự khác biệt về chất lượng cây giữa hai mô hình không có ý nghĩa thống kê.
  • Quy luật phân bố: Phân bố số cây theo đường kính $N/D_{1.3}$ ở cả hai mô hình đều phù hợp với hàm Weibull (Mô hình 1 có $\chi^2 = 1,97 < 11,07$; Mô hình 2 có $\chi^2 = 4,19 < 9,49$) và đều có dạng lệch trái ($\alpha < 3$). Ngược lại, phân bố số cây theo chiều cao $N/H_{vn}$ không tuân theo quy luật phân bố Weibull (Mô hình 1 có $\chi^2 = 42,76 > 11,07$; Mô hình 2 có $\chi^2 = 28,25 > 12,59$).
  • Tương quan $H_{vn}/D_{1.3}$: Phương trình thực nghiệm dạng logarit thể hiện độ tin cậy khá: Mô hình 1 có $H = -9,37 + 17,48 \log D_{1.3}$ ($R^2 = 0,498$); Mô hình 2 có $H = 0,69 + 0,04 \log D_{1.3}$ ($R^2 = 0,678$).
  • Thổ nhưỡng: Đất thuộc nhóm đất feralit vàng phát triển trên đá mẹ phiến thạch sét, tầng trung bình, thành phần cơ giới thịt nặng, kết cấu viên, hơi chặt, ít đá lẫn, tầng mặt mỏng do xói mòn; lớp thảm mục dày 2 - 3,5 cm.
  • Hiệu quả kinh tế (chu kỳ 25 năm, lãi suất 3%): Mô hình 1 đạt NPV = 277.235.644 đồng/ha, BCR = 6,25, IRR = 27%, thu nhập bình quân đạt 11,09 triệu đồng/ha/năm. Mô hình 2 đạt NPV = 288.695.240 đồng/ha, BCR = 6,33, IRR = 27%, thu nhập bình quân đạt 11,55 triệu đồng/ha/năm.

Phần V: Kết luận - Tồn tại - Kiến nghị

Chương tổng kết các phát hiện cốt lõi về tác động của mật độ tới sinh trưởng cá thể và sản lượng quần thể, đồng thời nêu rõ các hạn chế nghiên cứu và đề xuất 7 nhóm giải pháp kỹ thuật lâm sinh cho địa phương.

Kết quả và đóng góp

Số liệu thực nghiệm tổng hợp

Chỉ tiêu nghiên cứu Đơn vị tính Mô hình 1 (900 cây/ha) Mô hình 2 (1060 cây/ha) Giá trị kiểm định / Mức sai khác
Đường kính ngang ngực bình quân ($D_{1.3}$) cm 17,64 16,78 $|U| = 2,23 > 1,96$ (sai khác rõ rệt)
Chiều cao vút ngọn bình quân ($H_{vn}$) m 12,30 13,01 $|U| = 3,284 > 1,96$ (sai khác rõ rệt)
Chiều cao dưới cành bình quân ($H_{dc}$) m 8,69 9,45 $|U| = 3,587 > 1,96$ (sai khác rõ rệt)
Đường kính tán bình quân ($D_t$) m 3,90 3,68 $|U| = 3,52 > 1,96$ (sai khác rõ rệt)
Thể tích cây đứng trung bình ($V$) $\text{m}^3/\text{cây}$ 0,1444 0,1419 Mô hình 1 lớn hơn
Trữ lượng lâm phần 14 tuổi ($M$) $\text{m}^3/\text{ha}$ 129,96 151,41 Mô hình 2 cao hơn $21,45\text{ m}^3/\text{ha}$
Lượng tăng trưởng trữ lượng bình quân ($\Delta M$) $\text{m}^3/\text{ha/năm}$ 9,28 10,82 Mô hình 2 cao hơn $1,54\text{ m}^3/\text{ha/năm}$
Dự tính trữ lượng cuối chu kỳ 25 năm $\text{m}^3/\text{ha}$ 232,00 270,75 Mô hình 2 cao hơn $38,75\text{ m}^3/\text{ha}$
Tỷ lệ cây phẩm chất tốt (A) % 65,19 61,25 $\chi^2 = 1,57 < 5,99$ (không sai khác)
Tỷ lệ cây phẩm chất trung bình (B) % 20,00 23,75 $\chi^2 = 1,57 < 5,99$ (không sai khác)
Tỷ lệ cây phẩm chất xấu (C) % 14,81 15,00 $\chi^2 = 1,57 < 5,99$ (không sai khác)
Tham số phân bố $N/D_{1.3}$ (Weibull) - $\alpha=2,2; \lambda=0,015$ $\alpha=2,5; \lambda=0,006$ Cả hai đều phù hợp ($\chi^2 < \chi^2_{0.05}$)
Hệ số xác định tương quan $H_{vn}/D_{1.3}$ $R^2$ 0,498 0,678 Tương quan trung bình đến tương đối chặt
Giá trị hiện tại thuần (NPV) đồng/ha 277.235.644 288.695.240 Mô hình 2 cao hơn 11,46 triệu đồng/ha
Tỷ số lợi ích - chi phí (BCR) lần 6,25 6,33 Mô hình 2 cao hơn
Suất thu hồi vốn nội bộ (IRR) % 27% 27% Cả hai mô hình đều có tính khả thi cao
Thu nhập bình quân năm triệu đồng/ha/năm 11,09 11,55 Mô hình 2 cao hơn

Hệ thống giải pháp kỹ thuật lâm sinh đề xuất

Dựa trên kết quả điều tra, tác giả đề xuất các giải pháp kỹ thuật lâm sinh cụ thể:

  1. Tỉa thưa điều chỉnh mật độ: Tiến hành tỉa thưa ở Mô hình 2 nhằm giảm bớt mật độ cây, mở rộng không gian dinh dưỡng giúp cây tăng trưởng đường kính tốt hơn; tận dụng sản phẩm tỉa thưa để tạo nguồn thu phụ.
  2. Phát dọn thực bì: Thực bì và dây leo phát triển mạnh cạnh tranh trực tiếp dinh dưỡng với cây thông, do đó cần phát dọn định kỳ kết hợp tạo băng cản lửa nhằm phòng chống cháy rừng trong mùa hanh khô.
  3. Tập huấn kỹ thuật cạo nhựa: Người dân địa phương còn cạo mủ tự do, sai kỹ thuật làm suy kiệt thân cây. Cần mở các lớp tập huấn kỹ thuật trích nhựa đạt tiêu chuẩn, có sự hướng dẫn của cán bộ khuyến nông, khuyến lâm xã.
  4. Phòng trừ sâu róm thông: Giám sát chặt chẽ dịch sâu róm thông; khi có dịch cần bắt sâu thủ công hoặc phun thuốc trừ sâu sinh học/hóa học nồng độ thích hợp (Vophatox, Thiophat), bảo vệ các loài thiên địch ký sinh (ong mắt đỏ, ong trứng đen); dọn dẹp cây chết, cành lá nhiễm bệnh để triệt tiêu nguồn lây.
  5. Trồng bổ sung loài cây bản địa: Tại Mô hình 1 (mật độ thưa), đề xuất trồng dặm bổ sung các loài cây ưa bóng giai đoạn non và phù hợp đất đai địa phương như hồi, dẻ ăn quả, re hương, sồi phảng, màng tang để chuyển hóa dần thành rừng hỗn loài, khắc phục nhược điểm của rừng thuần loài.
  6. Xúc tiến tái sinh tự nhiên: Tận dụng khả năng tái sinh hạt tốt của Thông đuôi ngựa dưới tán rừng bằng cách phát luỗng thực bì, cuốc rãnh hoặc cày xới đất cục bộ giúp hạt dễ tiếp xúc với đất khoáng, sớm nảy mầm và tránh côn trùng phá hại, chuẩn bị thế hệ rừng kế thừa cho chu kỳ sau.
  7. Phòng cháy chữa cháy rừng: Xây dựng phương án phòng cháy nghiêm ngặt vào các tháng cao điểm mùa khô và mùa hè do thảm mục lá thông dày và trong thân, lá chứa nhiều tinh dầu dễ bắt lửa.

Đóng góp thực tiễn của công trình

Công trình đã cung cấp bộ dẫn liệu định lượng về sinh trưởng cá thể, trữ lượng và cơ cấu phẩm chất của Thông đuôi ngựa 14 tuổi tại vùng núi huyện Cao Lộc. Nghiên cứu chứng minh rằng dù mô hình mật độ thưa (900 cây/ha) cho đường kính thân cây và tán lá phát triển vượt trội, mô hình mật độ dày hơn (1060 cây/ha) lại đạt tổng sản lượng gỗ và giá trị tài chính tổng hợp (NPV, BCR) cao hơn trên đơn vị diện tích trong chu kỳ 25 năm, tạo căn cứ để các nhà quản lý rừng địa phương lựa chọn phương thức kinh doanh rừng phù hợp với mục tiêu lấy gỗ lớn hay tối ưu hóa tổng thu nhập.

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

Khóa luận tự xác định hai tồn tại kỹ thuật xuất phát từ giới hạn thời gian và điều kiện thực địa:

  • Dự tính kinh tế chưa phản ánh đầy đủ rủi ro: Các mô hình rừng nghiên cứu mới ở tuổi 14, chưa kết thúc chu kỳ kinh doanh khai thác chính (25 năm). Toàn bộ tính toán dòng tiền từ năm thứ 15 đến năm thứ 25 mang tính dự báo lý thuyết, chưa lượng hóa được các tổn thất kinh tế nếu xảy ra biến cố thiên tai, cháy rừng hay dịch hại sâu róm bùng phát trên diện rộng.
  • Thiếu hụt dữ liệu vi khí hậu: Nghiên cứu chưa đo đếm và phân tích chuyên sâu mối tương quan định lượng giữa các nhân tố vi khí hậu dưới tán (cường độ ánh sáng, nhiệt độ không khí, độ ẩm tầng đất) với tốc độ sinh trưởng của lâm phần.

Hướng mở rộng nghiên cứu được kiến nghị:

  • Tiếp tục khảo sát, theo dõi sinh trưởng và lượng hóa hiệu quả kinh tế trên nhiều loài cây trồng rừng khác tại địa phương để so sánh.
  • Mở rộng phạm vi đánh giá sang các chỉ tiêu giá trị môi trường sinh thái, bao gồm khả năng phòng hộ đầu nguồn, chống xói mòn và năng lực hấp thụ, tích lũy carbon của rừng trồng Thông đuôi ngựa trên các dạng lập địa khác nhau.

Giá trị tham khảo

Bản khóa luận có giá trị tham khảo thiết thực cho:

  • Sinh viên và học viên ngành Lâm học, Quản lý tài nguyên rừng, Kinh tế nông lâm nghiệp: Cung cấp tài liệu tham khảo về quy trình điều tra lâm học thực địa, kỹ thuật lập ô tiêu chuẩn, phương pháp xử lý số liệu thống kê sinh trắc bằng Excel Descriptive Statistics, ứng dụng kiểm định $U$, phân bố Weibull và tính toán các chỉ số tài chính (NPV, BCR, IRR) cho dự án lâm nghiệp.
  • Kỹ sư lâm sinh và cán bộ kỹ thuật tại các ban quản lý rừng, công ty lâm nghiệp, trạm khuyến nông - khuyến lâm miền núi: Tham khảo các định mức chi phí trồng, chăm sóc, bảo vệ và định mức thu hoạch lâm sản (nhựa, gỗ tỉa thưa, gỗ khai thác chính) của Thông đuôi ngựa chu kỳ 25 năm; vận dụng các khuyến nghị về mật độ trồng, kỹ thuật tỉa thưa và giải pháp xúc tiến tái sinh tự nhiên dưới tán.

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

1. Tại sao rừng trồng Thông đuôi ngựa ở mật độ 900 cây/ha có đường kính lớn hơn nhưng trữ lượng lại thấp hơn mô hình 1060 cây/ha?

Ở mật độ 900 cây/ha, cây rừng có không gian dinh dưỡng và tán lá rộng hơn ($D_t = 3,90\text{ m}$ so với $3,68\text{ m}$), giúp khả năng quang hợp tốt hơn và đường kính ngang ngực đạt cao hơn ($17,64\text{ cm}$ so với $16,78\text{ cm}$, chênh lệch 0,86 cm với kiểm định $U = 2,23 > 1,96$). Thể tích bình quân một cây ở mật độ 900 cây/ha lớn hơn ($0,1444\text{ m}^3$ so với $0,1419\text{ m}^3$). Tuy nhiên, do số lượng cây đứng trên 1 ha của mô hình 1060 cây/ha nhiều hơn, tổng trữ lượng thân cây của mô hình 1060 cây/ha đạt $151,41\text{ m}^3/\text{ha}$, cao hơn mức $129,96\text{ m}^3/\text{ha}$ của mô hình 900 cây/ha.

2. Quy luật phân bố số cây theo đường kính ($N/D_{1.3}$) và chiều cao ($N/H_{vn}$) của hai mô hình thể hiện như thế nào?

Phân bố số cây theo đường kính $N/D_{1.3}$ của cả hai mô hình đều tuân theo hàm phân bố Weibull với mức ý nghĩa 95% (Mô hình 1 có $\chi^2 = 1,97 < \chi^2_{0.05} = 11,07$; Mô hình 2 có $\chi^2 = 4,19 < \chi^2_{0.05} = 9,49$) và đều có dạng lệch trái do tham số hình dạng $\alpha < 3$ ($\alpha = 2,2$ và $\alpha = 2,5$). Ngược lại, phân bố số cây theo chiều cao $N/H_{vn}$ không phù hợp với phân bố lý thuyết Weibull (Mô hình 1 có $\chi^2 = 42,76$; Mô hình 2 có $\chi^2 = 28,25$, đều vượt ngưỡng tra bảng $\chi^2_{0.05}$). Thực nghiệm cho thấy đường cong chiều cao của Mô hình 1 có dạng gần đối xứng, trong khi Mô hình 2 có một đỉnh hơi lệch trái.

3. Chất lượng cây rừng giữa hai mật độ có sự khác biệt mang ý nghĩa thống kê hay không?

Không có sự khác biệt mang ý nghĩa thống kê. Tỷ lệ cây phẩm chất tốt và trung bình chiếm chủ đạo ở cả hai mô hình (Mô hình 1 đạt 85,19%; Mô hình 2 đạt 85,00%). Tỷ lệ cây xấu ở Mô hình 1 là 14,81% và Mô hình 2 là 15,00%. Khi kiểm định tính sai dị chất lượng bằng tiêu chuẩn $\chi^2$, kết quả đạt $\chi^2 = 1,57$, nhỏ hơn giá trị tra bảng $\chi^2_{0.05; 2} = 5,99$, chứng minh sự phân hóa phẩm chất cây đứng giữa hai cấp mật độ là không đáng kể.

4. Kết quả đánh giá hiệu quả kinh tế chu kỳ 25 năm nghiêng về mô hình nào?

Cả hai mô hình đều đạt hiệu quả tài chính cao với suất thu hồi vốn nội bộ IRR đạt 27% (vượt xa mức lãi suất vay vốn 3% thời điểm năm 1998). Mô hình 2 (1060 cây/ha) có các chỉ tiêu tài chính nhỉnh hơn: NPV đạt 288.695.240 đồng/ha (so với 277.235.644 đồng/ha ở Mô hình 1), BCR đạt 6,33 (so với 6,25 ở Mô hình 1) và thu nhập bình quân đạt 11,55 triệu đồng/ha/năm (so với 11,09 triệu đồng/ha/năm). Vì vậy, xét trên mục tiêu tối ưu hóa sản lượng và doanh thu trồng rừng sản xuất, Mô hình 2 được đánh giá là lựa chọn phù hợp hơn.

Kết luận

Khóa luận tốt nghiệp của tác giả Dương Thị Tuyền đã phản ánh toàn diện tình hình sinh trưởng, đặc điểm cấu trúc lâm phần và hiệu quả tài chính của hai mô hình rừng trồng Thông đuôi ngựa 14 tuổi tại xã Gia Cát, huyện Cao Lộc, tỉnh Lạng Sơn. Kết quả phân tích sinh trắc học và kiểm định thống kê khẳng định mật độ trồng có ảnh hưởng rõ rệt đến sự phân hóa kích thước cá thể và trữ lượng lâm phần: mật độ thưa (900 cây/ha) ưu thế về đường kính và kích thước tán lá, trong khi mật độ dày (1060 cây/ha) ưu thế về chiều cao vút ngọn, khả năng tỉa cành tự nhiên, tổng trữ lượng gỗ và giá trị hiện tại thuần (NPV). Hệ thống số liệu thực nghiệm và 7 giải pháp kỹ thuật lâm sinh được đề xuất trong công trình là tài liệu tham khảo có giá trị ứng dụng thực tiễn cho công tác quản lý, kinh doanh rừng trồng Thông đuôi ngựa bền vững tại vùng Đông Bắc Việt Nam.