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

Trong bối cảnh nhu cầu sử dụng gỗ phục vụ xây dựng và gia dụng ngày càng gia tăng, việc thực thi chính sách đóng cửa rừng tự nhiên cùng với nguồn cung ứng gỗ nhập khẩu dần bị thu hẹp đặt ra yêu cầu cấp thiết về việc phát triển rừng trồng sản xuất. Cây Xoan ta (Melia azedarach L.), thuộc họ Xoan (Meliaceae), là loài cây gỗ mọc nhanh bản địa, có khả năng thích ứng rộng, gỗ nhẹ mềm ít bị mối mọt, phục vụ hiệu quả cho nhu cầu làm nhà và đóng đồ gia dụng. Theo Quyết định số 16/2005/QĐ-BNN ngày 15/3/2005 của Bộ trưởng Bộ Nông nghiệp và Phát triển Nông thôn, Xoan ta được xếp vào danh mục các loài cây trồng rừng sản xuất chủ yếu cho 6/9 vùng sinh thái lâm nghiệp tại Việt Nam.

Nhằm cung cấp cơ sở khoa học cho công tác quản lý, nuôi dưỡng và dự báo sản lượng rừng giống, đề tài "Nghiên cứu sinh trưởng của rừng trồng Xoan ta (Melia azedarach L.) trồng trong mô hình rừng giống cây đầu dòng tại Viện Nghiên cứu và Phát triển lâm nghiệp - Trường Đại học Nông Lâm Thái Nguyên" được thực hiện với các mục tiêu cụ thể:

  1. Đánh giá sinh trưởng, tăng trưởng của rừng giống Xoan ta thuộc Viện Nghiên cứu và Phát triển lâm nghiệp - Trường Đại học Nông Lâm Thái Nguyên.
  2. Phân tích các quy luật kết cấu lâm phần rừng giống Xoan ta tại khu vực nghiên cứu.
  3. Lập mô hình biểu diễn mối quan hệ giữa các chỉ tiêu sản lượng rừng, mật độ lâm phần loài Xoan ta làm cơ sở xây dựng mô hình sản lượng và dự đoán trữ lượng gỗ lâm phần.

Đối tượng nghiên cứu là lâm phần Xoan ta 5 tuổi có nguồn gốc từ hạt giống thu hái từ cây trội tại tỉnh Thanh Hóa, gieo ươm và trồng tại Viện Nghiên cứu và Phát triển lâm nghiệp. Phạm vi nghiên cứu tập trung vào các đặc điểm sinh trưởng đường kính thân cây ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), đường kính tán ($D_t$) và cấu trúc lâm phần giai đoạn 5 tuổi trong khoảng thời gian thực tập từ 30/11/2015 đến 16/4/2016. Địa điểm nghiên cứu đặt tại xã Quyết Thắng, thành phố Thái Nguyên, tỉnh Thái Nguyên.

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

Khung lý thuyết của nghiên cứu dựa trên các nguyên lý điều tra rừng, phân loại học thực vật và lý thuyết mô hình hóa cấu trúc - sinh trưởng lâm phần. Tác giả kế thừa hệ thống lý thuyết về mô hình sinh trưởng thực nghiệm (empirical model), mô hình động thái (process model) và mô hình hỗn hợp (hybrid model) từ các nghiên cứu quốc tế (Vanclay, Landsberg & Gower) và các nghiên cứu lâm sinh tại Việt Nam (Trịnh Đức Huy, Nguyễn Ngọc Lung, Bảo Huy, Đào Công Khanh, Vũ Văn Thông). Về mặt pháp lý, đề tài bám sát định hướng cơ cấu cây trồng theo Quyết định 16/2005/QĐ-BNN của Bộ Nông nghiệp và Phát triển Nông thôn.

Phương pháp nghiên cứu kết hợp giữa điều tra thực địa toàn diện (ngoại nghiệp) và phân tích thống kê toán học (nội nghiệp):

  • Kế thừa dữ liệu và điều tra ngoại nghiệp: Thu thập hồ sơ lý lịch rừng trồng thiết lập tháng 7/2011 với quy mô 3,5 ha, tổng số cây ban đầu 4.260 cây, mật độ 1.100 cây/ha (cự ly $3\text{ m} \times 3\text{ m}$), hố đào $50 \times 50 \times 50\text{ cm}$, bón lót phân chuồng, NPK và vôi bột theo sơ đồ khối vuông la-tinh. Tiến hành điều tra tỉ mỉ toàn diện 100% số cây hiện còn trong lâm phần ($n = 2.158\text{ cây}$). Đo chu vi $C_{1.3}$ quy đổi sang đường kính ngang ngực $D_{1.3} = C_{1.3} / 3{,}14$; đo chiều cao vút ngọn $H_{vn}$ bằng thước sào kết hợp thước Blume-Leiss; đo đường kính tán $D_t$ theo hai hướng trực giao Đông - Tây và Nam - Bắc; phân loại phẩm chất cây theo 3 cấp: Tốt, Trung bình, Xấu.
  • Phương pháp nội nghiệp và mô hình hóa:
    • Tính toán các tham số thống kê mô tả: giá trị trung bình ($\bar{X}$), phương sai ($S^2$), độ lệch chuẩn ($S$), hệ số biến động ($S%$).
    • Mô phỏng quy luật phân bố số cây theo cỡ kính ($N/D_{1.3}$) bằng hàm mật độ xác suất Weibull với tham số dạng $\alpha = 3$ (dạng đối xứng), kiểm định độ tương thích giữa tần số thực nghiệm và lý thuyết bằng tiêu chuẩn $\chi^2$ (Chi-bình phương) của Pearson ở mức ý nghĩa 5%.
    • Mô phỏng động thái sinh trưởng đường kính và chiều cao theo tuổi ($A$) bằng hàm phi tuyến Schumacher: $Y = b_0 \cdot \exp(-b_1 / A^{b_2})$ thông qua phân tích hồi quy phi tuyến trên phần mềm SPSS 20.
    • Xác định các chỉ tiêu tăng trưởng: lượng tăng trưởng thường xuyên hàng năm ($Z_t$), tăng trưởng định kỳ ($Z_{n.t}$), tăng trưởng bình quân chung ($\Delta_t$), tăng trưởng bình quân định kỳ ($\Delta_{nt}$) và suất tăng trưởng ($p%$).
    • Phân tích tương quan hai chiều giữa $H_{vn} - D_{1.3}$ và $D_t - D_{1.3}$ bằng công cụ hồi quy tuyến tính trong Microsoft Excel.

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

Phần 1: Mở đầu

Trình bày tính cấp thiết của việc phát triển nguồn nguyên liệu gỗ từ rừng trồng thay thế gỗ tự nhiên, phân tích giá trị kinh tế và sinh thái của loài Xoan ta, xác định 3 mục tiêu nghiên cứu cụ thể, đồng thời nêu rõ ý nghĩa học thuật trong việc rèn luyện kỹ năng điều tra lâm học và ý nghĩa thực tiễn đối với kỹ thuật chăm sóc, quản lý rừng giống.

Phần 2: Tổng quan nghiên cứu

Tổng hợp các cơ sở khoa học về phân loại thực vật học, hình thái, sinh thái học của loài Xoan ta. Khái quát bức tranh nghiên cứu trong nước và quốc tế về mô hình sinh trưởng (hàm Gompertz, Schumacher), cấu trúc phân bố số cây theo cấp kính và chiều cao (Weibull, Beta, Gamma, Meyer), quy luật tương quan hình thái ($H/D$, $D_t/D$), và các nghiên cứu sinh khối rừng.

Chương này cũng mô tả chi tiết điều kiện khu vực nghiên cứu tại xã Quyết Thắng: địa hình đồi bát úp lượn sóng độ dốc $10 - 15^\circ$, khí hậu nhiệt đới gió mùa lượng mưa trung bình $2.000 - 2.500\text{ mm/năm}$, nền đất feralit phát triển trên đá phiến thạch sét ($Fs$) có tính chất chua mạnh ($\text{pH} < 4{,}0$) và nghèo dinh dưỡng:

Độ sâu tầng đất (cm) Mùn (%) N tổng số (%) $\text{P}_2\text{O}_5$ (%) $\text{K}_2\text{O}$ (%) N dễ tiêu (mg/100g) $\text{P}_2\text{O}_5$ dễ tiêu (mg/100g) $\text{K}_2\text{O}$ dễ tiêu (mg/100g) pH
1 - 10 1,766 0,024 0,241 0,035 3,64 4,56 0,90 3,5
10 - 30 0,670 0,058 0,211 0,060 3,06 0,12 0,12 3,9
30 - 60 0,711 0,034 0,131 0,107 0,107 3,04 3,04 3,7

Phần 3: Đối tượng, phạm vi, nội dung và phương pháp nghiên cứu

Xác lập khung phương pháp luận của đề tài, định rõ quy mô đối tượng nghiên cứu là $2.158$ cây Xoan ta 5 tuổi trên diện tích 3,5 ha. Trình bày chi tiết các bước thu thập số liệu ngoại nghiệp, công thức tính toán toán học thống kê, thuật toán nắn phân bố Weibull, phương pháp ước lượng tham số hàm Schumacher trên SPSS, các công thức xác định chỉ tiêu tăng trưởng lâm phần ($Z_t, Z_{n.t}, \Delta_t, p%$) và phương pháp kiểm định thống kê Pearson.

Phần 4: Kết quả nghiên cứu và thảo luận

Đây là chương nội dung trọng tâm, cung cấp toàn bộ dữ liệu định lượng thu được từ đợt điều tra thực tế:

  • Đặc trưng sinh trưởng và phẩm chất: Lâm phần 5 tuổi có mật độ bình quân $617\text{ cây/ha}$, đường kính bình quân $\bar{D}{1.3} = 10{,}09\text{ cm}$ ($S% = 22{,}22%$), chiều cao vút ngọn bình quân $\bar{H}{vn} = 9{,}13\text{ m}$ ($S% = 21{,}41%$). Phân cấp chất lượng ghi nhận $939\text{ cây}$ loại Tốt (chiếm $43{,}5%$), $855\text{ cây}$ loại Trung bình (chiếm $39{,}6%$) và $364\text{ cây}$ loại Xấu (chiếm $16{,}9%$). Tỷ lệ cây phẩm chất trung bình và xấu còn cao do ảnh hưởng của phân hóa tự nhiên và mật độ khuyết cây.
  • Quy luật sinh trưởng và tăng trưởng đường kính, chiều cao:
    • Mô hình sinh trưởng đường kính Schumacher: $D_{1.3} = 44{,}314 \cdot \exp(-3{,}384 / A^{0{,}532})$ với hệ số xác định $R^2 = 0{,}921$. Lượng tăng trưởng đường kính đạt đỉnh ở giai đoạn 2 - 3 tuổi ($Z_d = 5{,}0\text{ cm}$, suất tăng trưởng $p_d = 64{,}52%$), sau đó giảm dần ở tuổi 4 ($Z_d = 1{,}4\text{ cm}$) và tuổi 5 ($Z_d = 0{,}94\text{ cm}$). Lượng tăng trưởng bình quân chung $\Delta_d = 2{,}018\text{ cm/năm}$.
    • Mô hình sinh trưởng chiều cao Schumacher: $H_{vn} = 22{,}185 \cdot \exp(-3{,}070 / A^{0{,}809})$ với $R^2 = 0{,}913$. Tăng trưởng chiều cao mạnh nhất ở tuổi 3 ($Z_h = 4{,}85\text{ m}$, suất tăng trưởng $p_h = 66{,}89%$), sau đó giảm xuống $Z_h = 1{,}35\text{ m}$ ở tuổi 4 và $Z_h = 0{,}5\text{ m}$ ở tuổi 5. Tăng trưởng bình quân chung $\Delta_h = 1{,}82\text{ m/năm}$.
Tuổi ($A$) $Z_d$ (cm) $Z_{n.d}$ (cm) $\Delta_{nd}$ (cm/năm) $\Delta_d$ (cm/năm) $p_d$ (%) $Z_h$ (m) $Z_{n.h}$ (m) $\Delta_{nh}$ (m/năm) $\Delta_h$ (m/năm) $p_h$ (%)
1 2,25 - - 2,018 100,00 1,85 - - 1,82 100,00
2 0,50 - - 2,018 18,19 0,55 - - 1,82 20,83
3 5,00 5,50 (1-3) 1,83 2,018 64,52 4,85 5,40 (1-3) 1,80 1,82 66,89
4 1,40 - - 2,018 15,30 1,35 - - 1,82 15,70
5 0,94 2,34 (3-5) 1,17 2,018 9,32 0,50 1,85 (3-5) 0,925 1,82 5,49
  • Cấu trúc phân bố số cây theo cỡ kính ($N/D_{1.3}$): Tổng diện tích ngang của lâm phần đạt $G = 18{,}082\text{ m}^2$. Phân bố tập trung chủ yếu ở các cỡ kính $7 - 13\text{ cm}$, cao nhất tại cỡ kính $9\text{ cm}$ ($705\text{ cây}$) và $11\text{ cm}$ ($673\text{ cây}$). Ứng dụng hàm Weibull với tham số scale $\lambda = 0{,}0032$, kết quả kiểm định cho thấy giá trị thực nghiệm $\chi^2_{tn} = 9{,}1502 < \chi^2_{0{,}05; k=4} = 9{,}4877$, khẳng định hàm Weibull hoàn toàn phù hợp để mô phỏng cấu trúc đường kính của lâm phần.
Cỡ kính $D_{1.3}$ (cm) Khoảng biến thiên Số cây thực nghiệm ($f_t$) Tổng diện ngang $g_i$ ($\text{m}^2$) Tần số lý thuyết Weibull ($f_{ll}$) $\chi^2$ kiểm định Pearson
5 4 - 6 44 0,102 54,209 1,9228
7 6 - 8 320 1,249 343,190 1,5675
9 8 - 10 705 4,449 674,850 1,3473
11 10 - 12 673 6,421 662,190 0,1766
13 12 - 14 344 4,583 333,830 0,3095
15 14 - 16 67 1,169 80,868 2,3783
17 16 - 18 5 0,109 8,511 1,4482
Tổng cộng 2.158 18,082 2.157,648 $\chi^2_{tn} = 9{,}1502$
  • Đề xuất giải pháp kỹ thuật: Đề xuất 2 nhóm biện pháp gồm kỹ thuật lâm sinh (tỉa cành tạo thân thẳng, tỉa thưa vệ sinh loại bỏ $16{,}9%$ cây phẩm chất xấu, bón vôi và phân hữu cơ cải tạo độ chua cho đất feralit) và quản lý bảo vệ (ngăn chặn chăn thả gia súc, phòng trừ sâu bệnh thân lá).

Phần 5: Kết luận và đề nghị

Khái quát các kết quả định lượng về chỉ tiêu sinh trưởng, kết cấu lâm phần và tính tương thích của mô hình Schumacher, Weibull; chỉ ra các tồn tại về thời gian theo dõi và kiến nghị tiếp tục thu thập dữ liệu trong các năm tiếp theo.

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

Nghiên cứu mang lại hệ thống dữ liệu thực nghiệm chi tiết về lâm phần rừng giống Xoan ta 5 tuổi tại Thái Nguyên:

  • Xác định được các chỉ tiêu sinh trưởng thực tế: mật độ còn lại $617\text{ cây/ha}$, đường kính trung bình $\bar{D}{1.3} = 10{,}09\text{ cm}$, chiều cao vút ngọn trung bình $\bar{H}{vn} = 9{,}13\text{ m}$, tổng diện tích ngang $G = 18{,}082\text{ m}^2$ trên diện tích 3,5 ha ($2.158\text{ cây}$).
  • Làm rõ quy luật sinh trưởng: Cây Xoan ta có giai đoạn tăng trưởng bùng phát mạnh nhất về cả đường kính và chiều cao ở năm tuổi thứ 3 ($Z_d = 5{,}0\text{ cm}$, $Z_h = 4{,}85\text{ m}$), sau đó tốc độ tăng trưởng giảm rõ rệt ở tuổi 4 và 5.
  • Thiết lập thành công hai phương trình sinh trưởng phi tuyến theo hàm Schumacher có độ tin cậy cao:
    • Sinh trưởng đường kính: $D_{1.3} = 44{,}314 \cdot \exp(-3{,}384 / A^{0{,}532})$ ($R^2 = 0{,}921$).
    • Sinh trưởng chiều cao: $H_{vn} = 22{,}185 \cdot \exp(-3{,}070 / A^{0{,}809})$ ($R^2 = 0{,}913$).
  • Chứng minh tính quy luật của cấu trúc lâm phần: Phân bố số cây theo cấp đường kính ($N/D_{1.3}$) tuân theo quy luật phân bố Weibull với $\alpha = 3$, được kiểm định chặt chẽ qua tiêu chuẩn Pearson ($\chi^2_{tn} = 9{,}1502 < \chi^2_{0{,}05} = 9{,}4877$).
  • Về giải pháp thực tiễn: Đề xuất phương án tỉa thưa loại bỏ $16{,}9%$ cây phẩm chất xấu, tỉa cành nuôi dưỡng thân thẳng và xử lý hóa tính đất chua (nâng độ pH từ mức $3{,}5 - 3{,}9$) nhằm nâng cao chất lượng rừng giống lấy hạt đầu dòng.

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

  • Hạn chế: Nghiên cứu mới chỉ thực hiện trên một mô hình rừng giống thuần loài 5 tuổi tại một vị trí lập địa cụ thể thuộc Viện Nghiên cứu và Phát triển lâm nghiệp (Trường Đại học Nông Lâm Thái Nguyên). Các dự báo sinh trưởng từ tuổi 6 trở đi chỉ là giá trị nội suy lý thuyết từ mô hình toán, chưa được kiểm chứng thực tế qua toàn bộ chu kỳ kinh doanh.
  • Hướng nghiên cứu tiếp: Cần duy trì hệ thống đo đếm định kỳ hàng năm trên lâm phần ở các giai đoạn tuổi tiếp theo để hiệu chỉnh các tham số mô hình Schumacher, hoàn thiện biểu sản lượng và đánh giá năng suất tạo hạt giống của cây đầu dòng khi rừng bước vào giai đoạn thành thục sinh sản.

Giá trị tham khảo

Bản khóa luận là tài liệu tham khảo hữu ích cho sinh viên, học viên cao học và cán bộ nghiên cứu chuyên ngành Lâm sinh, Quản lý tài nguyên rừng và Điều tra quy hoạch rừng. Các nội dung có giá trị tham khảo nổi bật bao gồm:

  • Quy trình ứng dụng hàm Schumacher trong phần mềm SPSS để mô hình hóa sinh trưởng cây rừng.
  • Phương pháp nắn phân bố thực nghiệm bằng hàm phân bố xác suất Weibull và các bước kiểm định độ phù hợp $\chi^2$ Pearson.
  • Cơ sở dữ liệu định lượng về đặc tính tăng trưởng theo tuổi của loài Xoan ta phục vụ cho việc lập kế hoạch kỹ thuật tỉa thưa, bón phân và chăm sóc rừng trồng gỗ nhanh.

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

1. Mật độ và tỷ lệ phân cấp chất lượng của rừng giống Xoan ta ở tuổi 5 đạt mức nào?
Mật độ lâm phần ở tuổi 5 đạt $617\text{ cây/ha}$ (tương ứng $2.158\text{ cây}$ trên diện tích 3,5 ha). Tỷ lệ phân cấp chất lượng gồm có $43{,}5%$ cây tốt ($939\text{ cây}$), $39{,}6%$ cây trung bình ($855\text{ cây}$) và $16{,}9%$ cây xấu ($364\text{ cây}$).

2. Phương trình mô tả sinh trưởng đường kính và chiều cao của lâm phần có dạng như thế nào?
Hai phương trình sinh trưởng được xác lập theo hàm Schumacher:

  • Đường kính: $D_{1.3} = 44{,}314 \cdot \exp(-3{,}384 / A^{0{,}532})$ với hệ số xác định $R^2 = 0{,}921$.
  • Chiều cao: $H_{vn} = 22{,}185 \cdot \exp(-3{,}070 / A^{0{,}809})$ với hệ số xác định $R^2 = 0{,}913$.

3. Ở giai đoạn tuổi nào cây Xoan ta trong mô hình có tốc độ tăng trưởng nhanh nhất?
Cây Xoan ta tăng trưởng nhanh nhất ở giai đoạn từ 2 đến 3 tuổi, với lượng tăng trưởng thường xuyên hàng năm về đường kính đạt $Z_d = 5{,}0\text{ cm}$ (suất tăng trưởng $64{,}52%$) và chiều cao đạt $Z_h = 4{,}85\text{ m}$ (suất tăng trưởng $66{,}89%$).

4. Kết quả nắn phân bố số cây theo cỡ đường kính ($N/D_{1.3}$) theo hàm Weibull được kiểm định ra sao?
Phân bố $N/D_{1.3}$ được nắn theo hàm Weibull với tham số dạng $\alpha = 3$ và $\lambda = 0{,}0032$. Giá trị kiểm định Chi-bình phương tính toán được là $\chi^2_{tn} = 9{,}1502$, nhỏ hơn giá trị tra bảng ở mức ý nghĩa $5%$ với bậc tự do $k = 4$ ($\chi^2_{0{,}05} = 9{,}4877$), chứng minh hàm Weibull phù hợp với phân bố thực nghiệm.

5. Điều kiện đất đai tại khu vực nghiên cứu có đặc điểm gì cần lưu ý đối với lâm sinh?
Khu vực nghiên cứu có nền đất feralit phát triển trên đá phiến thạch sét có độ chua cao ($\text{pH}$ tầng mặt $1 - 10\text{ cm}$ là $3{,}5$; tầng $10 - 30\text{ cm}$ là $3{,}9$), hàm lượng mùn và các chất dinh dưỡng dễ tiêu ($N, \text{P}_2\text{O}_5, \text{K}_2\text{O}$) ở mức thấp, đòi hỏi phải bón vôi khử chua và bổ sung phân bón cải tạo đất khi chăm sóc rừng.

Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Xuân Diện đã đánh giá toàn diện các đặc trưng sinh trưởng, tăng trưởng và quy luật kết cấu của rừng giống Xoan ta 5 tuổi tại Viện Nghiên cứu và Phát triển lâm nghiệp - Trường Đại học Nông Lâm Thái Nguyên. Bằng phương pháp điều tra tỉ mỉ $2.158$ cây và áp dụng các mô hình toán học thực nghiệm, đề tài đã mô phỏng thành công động thái sinh trưởng qua hàm Schumacher ($R^2 > 0{,}91$) và quy luật cấu trúc đường kính qua hàm Weibull ($\chi^2_{tn} = 9{,}1502$). Các số liệu và giải pháp kỹ thuật lâm sinh được đề xuất cung cấp căn cứ thực tiễn quan trọng cho việc nuôi dưỡng, tỉa thưa và quản lý mô hình rừng giống cây đầu dòng đạt hiệu quả cao.