Lâm học CHỌN TUỔI CƠ SỞ THÍCH HỢP ĐỂ ƯỚC LƯỢNG CHỈ SỐ LẬP ĐỊA ĐỐI VỚI RỪNG TRỒNG KEO LAI (Acacia auriculiformis*mangium) Ở TỈNH ĐỒNG NAI Trần Thị Ngoan1, Lê Bá Toàn2 Phân hiệu Trường Đại học Lâm nghiệp Hội Khoa học kỹ thuật Lâm nghiệp TP HCM TÓM TẮT Bài báo giới thiệu kết nghiên cứu phân chia cấp số lập địa rừng trồng Keo lai tỉnh Đồng Nai Mục tiêu nghiên cứu xác định tuổi sở thích hợp để phân chia số lập địa rừng trồng Keo lai Nghiên cứu tiến hành với rừng trồng Keo lai từ - 10 tuổi lập địa khác tỉnh Đồng Nai Cấp số lập địa (SI) rừng trồng Keo lai phân chia theo chiều cao ưu Số liệu thu thập bao gồm 111 mẫu mẫu điển hình với kích thước 0,1 ha; 108 sử dụng để xây dựng hàm SI, để kiểm tra khả ứng dụng hàm SI Chiều cao ưu xác định phương pháp giải tích thân Kết nghiên cứu tuổi sở thích hợp để phân chia số lập địa rừng trồng Keo lai năm Chỉ số lập địa rừng trồng Keo lai phân chia thành cấp từ I – III; số lập địa I, II III tương ứng với chiều cao 24, 20 16 m Để nâng cao hiệu kinh doanh, chủ rừng cần tập trung trồng rừng Keo lai cấp số lập địa I II Từ khóa: Cây ưu thế, số lập địa, đường cong số lập địa, rừng trồng Keo lai, tuổi sở I ĐẶT VẤN ĐỀ Chỉ số lập địa chiều cao trung bình ưu đồng ưu loài gỗ định tuổi sở (Monserud, 1984) Sở dĩ lâm học sử dụng chiều cao trung bình ưu đồng ưu để biểu thị cho số lập địa phụ thuộc vào chất lượng lập địa, mà không phụ thuộc vào mật độ quần thụ Tuy vậy, chiều cao trung bình ưu đồng ưu nhận giá trị thấp quần thụ có mật độ cao thấp Chỉ số lập địa sử dụng để đánh giá suất tương đối hay suất tiềm lập địa định Chọn tuổi sở thích hợp vấn đề quan trọng việc xác định số lập địa Nếu chọn tuổi sở nhỏ, kết nhận số lập địa bị trệch Trái lại, chọn tuổi sở cao, kết dự đốn số lập địa với biến động lớn Thông thường tuổi sở chọn thời điểm sinh trưởng chiều cao gỗ quần thụ khơng phụ thuộc vào mật độ Nói chung, tuổi sở chọn tương ứng với chu kỳ kinh doanh sau lượng tăng trưởng bình quân năm chiều cao đạt cao (Monserud, 1984; Larsen, 1999) Tuổi sở chọn tuổi mà chiều cao tuổi khác dự đoán từ hàm số lập địa với sai lệch nhỏ (Onyekwelu, 2003) Giá trị số lập địa thay đổi tùy theo loài Đối với lồi cây, số lập địa phân chia thành – cấp (Monserud, 1984; Larsen, 1999; Vũ Tiến Hinh, 2005) Cho đến có số cơng trình nghiên cứu chọn giống, kỹ thuật trồng nuôi dưỡng rừng trồng Keo lai (Lê Đình Khả, 2000) Theo Nguyễn Huy Sơn ctv (2006), tuổi thành thục công nghệ rừng trồng Keo lai Đông Nam Bộ năm Theo thống kê đến năm 2016, diện tích rừng trồng Keo lai tỉnh Đồng Nai tập trung chủ yếu khu vực Vĩnh Cửu, Xuân Lộc, Định Quán với tổng diện tích 23.557 Tuy vậy, chưa có nghiên cứu phân chia số lập địa rừng trồng Keo lai phạm vi toàn tỉnh Đồng Nai Bài báo giới thiệu kết nghiên cứu phân chia số lập địa rừng trồng Keo lai tỉnh Đồng Nai Kết nghiên cứu sở khoa học cho quản lý rừng phương thức lâm sinh II PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Vật liệu nghiên cứu TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 6-2017 51 Lâm học Nghiên cứu thực phạm vi toàn tỉnh Đồng Nai Tọa độ địa lý: 10030’03” - 11034’57” vĩ độ Bắc, 106045’30” 107035’00” kinh độ Đông Khu vực nghiên cứu nằm vùng khí hậu nhiệt đới gió mùa cận xích đạo Hàng năm khí hậu phân chia thành mùa mưa khô rõ rệt Mùa mưa từ tháng đến tháng 10, cịn mùa khơ từ tháng từ 11 năm trước đến tháng năm sau Tổng lượng xạ cao ổn định (390 - 556 cal/cm2/ngày) Nhiệt độ không khí trung bình 22,00C Lượng mưa trung bình năm 2.100 mm Độ ẩm khơng khí trung bình 80% Tổng tích ơn lớn (9.417 - 9.782oC/năm), số nắng nhiều (2.475,7 giờ/năm), bão Địa hình đồi thấp với độ cao từ 50 - 350 m so với mặt biển Đất có nhóm chính: đất hình thành đá bazan, đất hình thành đá phiến sét, đất hình thành phù sa cổ Đối tượng nghiên cứu rừng trồng Keo lai từ – 10 tuổi 2.2 Phương pháp thu thập số liệu (i) Số lượng mẫu ô mẫu: Trong nghiên cứu này, số lập địa rừng trồng Keo lai phân chia theo chiều cao ưu (H0, m) tuổi sở (A0, năm) Những mẫu thu thập từ ô mẫu 1.000 m2 thuộc rừng trồng Keo lai tuổi từ - 10 năm Tổng số mẫu thu thập 111 Những mẫu chọn từ quần thụ sinh trưởng tốt đến quần thụ sinh trưởng xấu Chỉ tiêu H0 bình quân tuổi xác định phương pháp giải tích thân (ii) Giải tích mẫu: Những giải tích sinh trưởng phát triển bình thường, khơng bị sâu bệnh hay cụt ngọn, thân thẳng tán trịn Những giải tích chặt hạ vị trí cách mặt đất 10 cm Mỗi giải tích thu thập đường kính ngang ngực (D, cm) chiều cao toàn thân (H0, m) Sau chặt hạ, tất mẫu đo đạc D chiều dài tồn thân thước dây với độ xác 0,10 cm Thân phân chia thành phân đoạn có chiều dài 1,0 m, 52 riêng đoạn cịn khoảng 0,5 – 1,0 m Sau cưa thớt giải tích vị trí 0,0 m; 1,0 m; 1,3 m; 2,0 m; 3,0 m; 4,0 m… Trên phân đoạn, thu thớt đường kính đầu lớn để xác định tuổi đạt đến chiều cao phân đoạn Những thớt giải tích tập hợp theo giải tích; sau ghi thứ tự vị trí thớt mặt thớt hướng phía 2.3 Phương pháp xử lý số liệu (i) Tính đặc trưng thống kê mơ tả H0 mẫu: Giá trị bình quân (H0bq, m), sai lệch chuẩn ước lượng (±S, m), hệ số biến động (CV%), giá trị nhỏ (H0min, m), giá trị lớn (H0max, m) biên độ H0max – H0min (ii) Xây dựng hàm ước lượng H0 = f(A) bình qn chung tồn rừng trồng Keo lai tỉnh Đồng Nai: Hàm H0 = f(A) bình quân chung ước lượng hàm Schumacher (1939) (Hàm 4) Số lượng mẫu để xây dựng hàm H0 = f(A) bình quân chung 108 Mức độ phù hợp hàm (3’) với số liệu thực nghiệm (H0TN, m) đánh giá theo hệ số xác định (r2), sai lệch chuẩn ước lượng (S), sai số tuyệt đối trung bình (MAE) sai số tuyệt đối trung bình theo phần trăm (MAPE) H0 = a*exp(-b/A^c) (1) Hay Ln(H0) = Ln(a) - b/A^c (2) Đặt Ln(a) = b0 b = b1, ta có: Ln(H0) = b0 + b1/A^c (3) Hay H0 = exp(b0 + b1/A^c) (3’) b0 = Ln(H0) – b1/A^c (4) Theo định nghĩa, số lập địa (SI) H0 A0 (năm) Theo đó, từ hàm (3’) hàm (4), hàm SI = f(A) xác định dạng hàm (5), tham số b0 xác định theo hàm (6) Ln(SI) = b0 + b1/A0^c (5) b0 = Ln(SI) – b1/A0^c (6) Thay b0 hàm (6) vào hàm (3’) biến đổi để nhận hàm (7) – (9) Hàm (9) hàm SI A0 hay H0 A0 Ln(H0) = Ln(SI) – b1/A0^c + b1/A^c (7) Ln(SI) = Ln(H0) – b1/A0^c + b1/A^c (8) TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ 6-2017 Lâm học SI = exp(Ln(H0) – b1(1/A^c -1/A0^c)) (9) (iii) Xác định tuổi sở (A0): Tuổi sở thích hợp chọn thời điểm mà hàm SI sử dụng để chuyển H0 A0 H0 tuổi A (năm) với tổng sai lệch bình phương nhỏ (SSRmin) Tuổi sở thích hợp kiểm định khoảng A = – 10 năm (iv) Xác định số lượng số SI: Số lượng số SI dựa theo biên độ biến động H0 A0 (Hmax – H0min) Độ dốc (tham số b1) hàm SI chọn Theo đó, trước hết xác định b1 từ hàm (3) Sau thay b1 H0 số SI A0 vào hàm (9) để nhận hàm SI Khả ứng dụng hàm SI kiểm định từ giải tích khơng tham gia xây dựng mơ hình Mức độ phù hợp hàm SI ba cấp số SI với hàm H0 = f(A) không tham gia xây dựng mơ hình kiểm định phương pháp so sánh ngang điểm chặn độ dốc mơ hình Sau xây dựng biểu số SI đường cong SI Biểu số SI đường cong SI xác định cách thay A H0 số SI A0 vào hàm (9) Từ lập bảng SI vẽ đường cong SI III KẾT QUẢ VÀ THẢO LUẬN 3.1 Những đặc trưng thống kê chiều cao ưu Những đặc trưng thống kê chiều cao ưu rừng trồng Keo lai từ – 10 tuổi khu vực nghiên cứu dẫn bảng Bảng Đặc trưng thống kê chiều cao ưu rừng trồng Keo lai từ – 10 tuổi tỉnh Đồng Nai A (năm) N (cây) H0 (m) ±S CV% H0min H0max H0max-H0min (1) (2) (3) (4) (5) (6) (7) (8) 108 2,6 0,4 15,4 2,1 3,3 1,2 108 6,7 1,6 23,9 4,0 9,1 5,1 108 11,0 1,9 17,3 7,1 14,2 7,1 108 14,0 2,4 17,1 9,4 17,9 8,5 108 15,9 2,9 18,2 10,5 20,4 9,9 108 18,0 2,7 15,0 12,2 22,7 10,5 108 18,5 3,7 20,0 12,3 24,4 12,1 108 20,3 3,8 18,7 13,6 26,8 13,2 108 21,9 3,3 15,1 14,8 27,6 12,8 10 108 22,4 3,3 14,7 15,7 29,0 13,3 Phân tích số liệu bảng cho thấy, giá trị H0 trung bình ưu biến đổi từ 2,6 m tuổi đến 22,4 m tuổi 10 Biên độ dao động H0 2,1 đến 3,3 m tuổi 15,7 đến 29,0 m tuổi 10 Hệ số biến động chiều cao (CV%) nhận giá trị cao (23,9%) tuổi 2, thấp (14,7%) tuổi 10 Nói chung, H0 rừng trồng Keo lai có biến động lớn theo tuổi điều kiện lập địa Vì thế, phân chia rừng trồng Keo lai tỉnh Đồng Nai thành cấp số lập địa khác việc làm cần thiết 3.2 Xây dựng hàm ước lượng H0 = f(A) rừng trồng Keo lai Những phân tích hồi quy tương quan cho thấy, hàm ước lượng H0 = f(A) rừng trồng Keo lai từ – 10 tuổi tỉnh Đồng Nai có dạng hàm (10) H0 = exp(3,65344 – 2,76734/A^0,707464) (10) r2 = 83,4%; S = ±2,8; MAE = 2,2; MAPE = 16,2% Bằng cách khảo sát hàm (10), xác định q trình sinh trưởng H0 rừng trồng Keo lai từ tuổi - 10 năm (Bảng 2; Hình 1) Từ cho thấy, lượng tăng trưởng thường xuyên hàng năm (ZH0) gia tăng dần từ TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 6-2017 53 Lâm học tuổi (2,4 m/năm) đạt cao tuổi (4,7 m/năm); sau giảm dần đến tuổi 10 (0,9 m/năm) Lượng tăng trưởng bình quân năm ( H0) gia tăng dần từ tuổi (2,4 m/năm) đạt cao tuổi (3,6 m/năm); sau giảm dần đến tuổi 10 (2,2 m/năm) Suất tăng trưởng chiều cao (Ph0%) giảm nhanh từ tuổi (100%) đến tuổi (14,1%) tuổi 10 (4,1%) Vì thế, tuổi thời kỳ H0 chuyển từ giai đoạn sinh trưởng nhanh sang giai đoạn sinh trưởng chậm Bảng Tăng trưởng H0 (m) rừng trồng Keo lai từ – 10 tuổi tỉnh Đồng Nai A (năm) (1) 10 H0 (m) (2) 2,4 7,1 10,8 13,7 15,9 17,7 19,2 20,4 21,5 22,4 ZH0 (m/năm) (3) 2,4 4,7 3,7 2,9 2,2 1,8 1,5 1,2 1,1 0,9 H0 (m/năm) (4) 2,4 3,5 3,6 3,4 3,2 3,0 2,7 2,6 2,4 2,2 Ph0% (5) 100,0 65,8 34,5 20,9 14,1 10,2 7,7 6,1 5,0 4,1 H0 (m) 5.0 4.5 4.0 3.5 3.0 ZH0 (m) 2.5 H0 (m) 2.0 1.5 1.0 0.5 10 A (năm) Hình Đồ thị biểu diễn tăng trưởng chiều cao ưu (H0, m) rừng trồng Keo lai tỉnh Đồng Nai 3.3 Xác định tuổi sở để xây dựng số lập địa Tuổi sở (A0, năm) thích hợp chọn thời điểm mà H0 tuổi A dự đoán theo hàm SI = f(A) với SSRmin Tại khu vực nghiên cứu, biên độ tuổi rừng trồng Keo lai dao động từ – 10 năm Vì thế, tuổi sở thích hợp kiểm định từ A = - 10 năm Kết qủa kiểm định sai lệch dự đoán H0 từ A0 = 54 – 10 năm hàm (10) ghi lại bảng bảng Từ cho thấy, SSR nhận giá trị cao A = năm (2,87), thấp A = năm (0,86) Vì thế, tuổi tuổi sở để xây dựng hàm SI đường cong SI rừng trồng Keo lai tỉnh Đồng Nai Tuổi phù hợp với tuổi thành thục công nghệ rừng trồng Keo lai miến Đông Nam Bộ (Nguyễn Huy Sơn ctv, 2006) TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 6-2017 Lâm học Bảng Dự đoán chiều cao ưu tuổi – 10 hàm SI = f(A) A0 = - 10 năm A (năm) H0 (m) (Thực tế) (1) 10 (2) 18,0 18,5 20,3 21,9 22,4 (3) 18,0 17,1 17,6 18,0 17,7 Giá trị dự đoán H0 (m) A0 (năm) (4) (5) (6) 19,5 20,8 21,9 18,5 19,7 20,7 19,1 20,3 21,4 19,5 20,8 21,9 19,2 20,4 21,5 10 (7) 22,8 21,6 22,3 22,8 22,4 Bảng Kiểm định sai lệch dự đoán H0 (m) tuổi – 10 A0 = - 10 năm A (năm) (1) 10 Tổng số (2) 0,00 0,87 0,17 0,00 0,10 1,14 Giá trị SSR tương ứng với A0 (năm) (3) (4) (5) 1,02 0,23 0,00 0,00 0,36 1,37 0,31 0,00 0,29 1,09 0,27 0,00 0,45 0,01 0,18 2,87 0,86 1,84 3.4 Xác định số lượng số lập địa tuổi sở Từ số liệu bảng cho thấy, biên độ dao động H0 tuổi 13,2 m (Hmax – H0min = 26,8 - 13,6 m) (làm tròn 13,0 m) Tại tuổi 8, H0 bình quân rừng trồng Keo lai 20,4 m (làm trịn 20,0 m) Thơng thường sai số đo chiều cao thân dao động từ ±0,5 – ±1,0 m Nếu phân chia biên độ dao động H0 tuổi (13,0 m) thành cấp, cấp 4,0 m Khoảng cách H0 hai cấp số lập địa kế cận (4,0 m) lớn - lần so với sai số đo chiều cao thân Vì thế, rừng trồng Keo lai 10 (6) 0,16 0,61 0,02 0,19 0,00 0,98 tỉnh Đồng Nai phân chia thành cấp số lập địa (I, II, III) dựa theo H0; khoảng cách hai cấp số lập địa kế cận A0 4,0 m Ba cấp số SI A0 (8 năm) nhận giá trị tương ứng 24,0 m, 20,0 m 16,0 m Các số SI ranh giới cấp số lập địa I II, II III tương ứng 22 m 18 m Tương tự, số SI biên độ cấp số lập địa III 14 m, biên độ cấp số lập địa I 26 m 3.5 Xây dựng hàm số lập địa rừng trồng Keo lai Bảng Các hàm số SI rừng trồng Keo lai Cấp SI Hàm số lập địa (1) (2) I(Trên) SI = exp((Ln(26) - 2,76734*(1/A^0,707464 – 0,22967))) (11) I SI = exp((Ln(24) - 2,76734*(1/A^0,707464 – 0,22967))) (12) II - I SI = exp((Ln(22) - 2,76734*(1/A^0,707464 – 0,22967))) (13) II SI = exp((Ln(20) - 2,76734*(1/A^0,707464 – 0,22967))) (14) II - III SI = exp((Ln(18) - 2,76734*(1/A^0,707464 – 0,22967))) (15) III SI = exp((Ln(16) - 2,76734*(1/A^0,707464 – 0,22967))) (16) III(Dưới) SI = exp((Ln(14) - 2,76734*(1/A^0,707464 – 0,22967))) (17) TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ 6-2017 Kí hiệu 55 Lâm học Trong nghiên cứu này, ba hàm SI chọn với độ dốc (b1) Những phân tích thống kê cho thấy, tham số b1 = 2,76734 1/8^0,707464 = 0,22967 Theo định nghĩa số SI, hàm SI (SIG) hàm SI giới hạn (SIGH) cấp số lập địa I – III ghi lại bảng 3.6 Kiểm định mức độ phù hợp hàm số lập địa So sánh ba hàm SI = f(A) ba số lập địa ba hàm H0 = f(A) mẫu không tham gia xây dựng mơ hình cho thấy, điểm chặn chúng khơng có khác biệt rõ rệt (P = 0,2390 cấp đất I; P = 0,2849 cấp đất II; P = 0,2608 cấp đất III) Tương tự, độ dốc chúng khơng có khác biệt rõ rệt (P = 0,5700 cấp đất I; P = 0,6113 cấp đất II; P = 0,3799 cấp đất III) Điều chứng tỏ hàm (11) – (17) sử dụng để xây dựng biểu đường cong số lập địa rừng trồng Keo lai tỉnh Đồng Nai 3.7 Xây dựng biểu đường cong SI rừng trồng Keo lai Bảng Biểu số lập địa rừng trồng Keo lai từ - 10 tuổi tỉnh Đồng Nai A (năm) (1) 10 I(Trên) (2) 3,1 9,0 13,7 17,4 20,2 22,5 24,4 26,0 27,3 28,5 H0 (m) ba cấp số lập địa I - II II II - III (4) (5) (6) 2,6 2,4 2,1 7,6 6,9 6,2 11,6 10,6 9,5 14,7 13,4 12,0 17,1 15,5 14,0 19,0 17,3 15,6 20,6 18,8 16,9 22,0 20,0 18,0 23,1 21,0 18,9 24,1 21,9 19,7 I (3) 2,8 8,3 12,7 16,0 18,7 20,8 22,5 24,0 25,2 26,3 H0 (m) III (7) 1,9 5,5 8,5 10,7 12,4 13,8 15,0 16,0 16,8 17,5 III(Dưới) (8) 1,7 4,9 7,4 9,4 10,9 12,1 13,1 14,0 14,7 15,3 SI (m) 30 24 25 20 20 16 15 10 5 10 A (năm) Hình Đường cong chiều cao ưu (H0, m) số lập địa (SI, m) rừng trồng Keo lai tỉnh Đồng Nai 56 TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ 6-2017 Lâm học Biểu dự đoán số lập địa đường cong số lập địa xây dựng cách thay tuổi rừng trồng Keo lai (A, năm) từ - 10 năm vào hàm (11) – (17) (Bảng 6, Hình 2) IV KẾT LUẬN Rừng trồng Keo lai tỉnh Đồng Nai phân chia thành ba cấp số lập địa Tuổi sở thích hợp để phân chia ba cấp số lập địa tuổi Các hàm biểu số lập địa rừng trồng Keo lai xây dựng dựa tuổi sở Biểu số lập địa rừng trồng Keo lai ứng dụng phạm vi từ tuổi – 10 năm Để nâng cao hiệu kinh doanh, chủ rừng cần tập trung trồng rừng Keo lai cấp số lập địa I II TÀI LIỆU THAM KHẢO Larsen, D.R (1999) Site index Natural Resource Biometrics Construction of site index equations for Pinus sylvestris L using permanent plot data in Sweden The School of Natural Resources, University of Missouri-Columbia Lê Đình Khả (2000) Nốt sần khả cải tạo đất Keo lai Tạp chí Lâm nghiệp, Số 6/2000 Monserud, R.A (1984) Height growth and site index curves for inland Doughlas-fir based on stem analysis data and forest habitat type Forest Science 30: 943–965 Nguyễn Huy Sơn, Nguyễn Văn Thịnh, Bùi Thanh Hằng, Nguyễn Thanh Minh, Phan Minh Sáng (2006) Nghiên cứu đặc điểm sinh trưởng Keo lai tuổi thành thục công nghệ rừng trồng Keo lai Đông Nam Bộ Tạp chí Khoa học Lâm nghiệp, Số 4/2006 Onyekwelu, J, C (2003) Choosing appropriate index age for estimating site index of Gmelina arborea timber plantations in Oluwa forest reserve Food, Agriculture & Environment Vol.1(3&4): 286-290 Vũ Tiến Hinh (2005) Sản lượng rừng Nxb Nông nghiệp, Hà Nội DETERMINATING AN APPROPRIATE INDEX AGE FOR ESTIMATING SITE INDEX OF ACACIA HYBRID PLANTATIONS AT DONG NAI PROVINCE Tran Thi Ngoan1, Le Ba Toan2 Vietnam National University of Forestry - Southern Campus Ho Chi Minh City Science and Technology Association of Forestry SUMMARY This article introduces the results of research on site index classes for Acacia hybrid plantations at Dong Nai province The objectives of research is identify appropriate site index age to divide up the site index for Acacia hybrid plantations at Dong Nai province In this study, the site index (SI) for Acacia hybrid plantations are divided according to the total height of dominant trees The collected data includes 111 sample trees in the typical plots with size of 0.1 ha; in which 108 trees were used to build the function SI, trees to examine possibilities of the functions SI The total height of dominant trees are determined by stem analysis methods Research results show that, appropriate index age to divide up the site index for Acacia hybrid plantations at Dong Nai province is years old The site index of Acacia hybrid plantations are divided into classes from IIII; in which classes of I, II and III respectively of 24, 20 and 16 m To improve the effectiveness of the business, the forest owner focus on forest planting of Acacia hybrid on the site classes of I or II Keywords: Acacia hybrid, dominant trees, site index, site index age, site index curve Ngày nhận Ngày phản biện Ngày định đăng : 27/10/2017 : 16/11/2017 : 25/11/2017 TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ 6-2017 57