a)Trạng thái rừng IIIA1
Kết quả tính toán một số đặc trưng mẫu của phân bố N/D1.3 trạng thái rừng IIIA1 được tổng hợp trong bảng 4.5.
Bảng 4.5. Một số đặc trưng mẫu của phân bố N/D1.3 trạng thái rừng IIIA1
Các đặc trưng mẫu của D1.3 OTC 1 OTC 2
Trung bình mẫu
D 1.3 (cm) 25.03 19.99 Sai số của trung bình mẫu Sx (cm) 0.96 0.69
Trung vị mẫu Me 21.96 18.14
Mode Mo 20.05 17.83
Sai tiêu chuẩn mẫu S 9.88 7.20
Phương sai mẫu S2 97.53 51.78
Độ nhọn của phân bố Ex 0.65 3.90
Độ lệch của phân bố Sk 1.03 1.64
Phạm vi biến động R (cm) 43.61 38.74
Giá trị quan sát nhỏ nhất Dmin (cm) 12.10 9.96 Giá trị quan sát lớn nhất Dmax (cm) 55.70 48.70
Dung lượng quan sát N (cây) 106 110
Qua bảng 4.5 cho thấy đường kính trung bình ở 2 OTC của trạng thái rừng IIIA1 tại khu vực nghiên cứu dao động từ 19,99 – 25,03 cm, sai tiêu chuẩn S = 7,20 – 9,88 cm; phương sai S2 = 51,75 – 97,93 cm2; phạm vi biến động đường kính từ 38,74cm đến 43,61cm. Cả 2 OTC đều có giá trị độ nhọn phân bố Ex > 0, chứng tỏ phân bố N/D1.3 nhọn hơn so với phân bố chuẩn. 2 OTC cũng có giá trị độ lệch phân bố SK > 0, chứng tỏ phân bố N/D1.3 của trạng thái rừng IIIA1 có dạng 1 đỉnh lệch trái so với phân bố chuẩn.
b)Trạng thái rừng IIIA2
Kết quả tính toán một số đặc trưng mẫu của phân bố N/D1.3 trạng thái rừng IIIA2 được tổng hợp trong bảng 4.6.
Bảng 4.6. Một số đặc trưng mẫu của phân bố N/D1.3 trạng thái rừng IIIA2
Các đặc trưng mẫu của D1.3 OTC 3 OTC 4
Trung bình mẫu
D 1.3 (cm) 28.05 26.71 Sai số của trung bình mẫu Sx (cm) 1.21 1.02
Trung vị mẫu Me 24.19 25.15
Mode Mo 21.01 18.46
Sai tiêu chuẩn mẫu S 13.05 11.18
Phương sai mẫu S2 170.37 124.99
Độ nhọn của phân bố Ex 0.59 0.29
Độ lệch của phân bố Sk 1.15 0.91
Phạm vi biến động R (cm) 54.75 48.54
Giá trị quan sát nhỏ nhất Dmin (cm) 11.46 10.98 Giá trị quan sát lớn nhất Dmax (cm) 66.21 59.52
Dung lượng quan sát N (cây) 117 121
Qua bảng 4.6 cho thấy đường kính trung bình ở 2 OTC của trạng thái rừng IIIA2 tại khu vực nghiên cứu dao động từ 26,71 – 28,05 cm, sai tiêu chuẩn S = 11,18 – 13,05 cm; phương sai S2 = 124,99 – 170,37 cm2; phạm vi biến động đường kính từ 48,54 cm đến 54,75cm. Cả 2 OTC đều có giá trị độ nhọn phân bố Ex > 0, chứng tỏ phân bố N/D1.3 nhọn hơn so với phân bố chuẩn. 2 OTC cũng có giá trị độ lệch phân bố SK > 0, chứng tỏ phân bố N/D1.3 của trạng thái rừng IIIA2 có dạng 1 đỉnh lệch trái so với phân bố chuẩn.
c) Trạng thái rừng IIIB
Kết quả tính toán một số đặc trưng mẫu của phân bố N/D1.3 trạng thái rừng IIIB được tổng hợp trong bảng 4.7.
Bảng 4.7. Một số đặc trưng mẫu của phân bố N/D1.3 trạng thái rừng IIIB
Các đặc trưng mẫu của D1.3 OTC 5 OTC 6
Trung bình mẫu
D 1.3 (cm) 30.63 30.77 Sai số của trung bình mẫu Sx (cm) 1.78 1.72
Trung vị mẫu Me 22.92 24.83
Mode Mo 22.28 22.28
Sai tiêu chuẩn mẫu S 20.22 19.49
Phương sai mẫu S2 408.86 379.86
Độ nhọn của phân bố Ex 2.59 4.81
Độ lệch của phân bố Sk 1.73 1.96
Phạm vi biến động R (cm) 90.40 104.09
Giá trị quan sát nhỏ nhất Dmin (cm) 10.82 6.37 Giá trị quan sát lớn nhất Dmax (cm) 101.22 110.45
Dung lượng quan sát N (cây) 129 129
Qua bảng 4.7 cho thấy đường kính trung bình ở 2 OTC của trạng thái rừng IIIB tại khu vực nghiên cứu dao động từ 30,63 – 30,77 cm, sai tiêu chuẩn S = 19,49 – 20,22 cm; phương sai S2 = 379,86 – 408,86 cm2; phạm vi biến động đường kính từ 90,40 cm đến 104,09 cm. Cả 2 OTC đều có giá trị độ nhọn phân bố Ex > 0, chứng tỏ phân bố N/D1.3 nhọn hơn so với phân bố chuẩn. 2 OTC cũng có giá trị độ lệch phân bố SK > 0, chứng tỏ phân bố N/D1.3 của trạng thái rừng IIIB có dạng 1 đỉnh lệch trái so với phân bố chuẩn.
4.3.2. Quy luật phân bố số cây theo cỡ đường kính (N/D1.3 )
Phân bố N/D1.3 thể hiện quy luật sắp xếp các thành phần cấu tạo nên quần thể cây rừng trong không gian và theo thời gian. Đây là cơ sở quan trọng cho việc thống kê, dự đoán trữ lượng, sản lượng rừng nên nó là quy luật quan trọng trong kết cấu lâm phần. Từ quy luật cấu trúc này, có thể đánh giá được kết cấu của rừng, hướng đề xuất các biện pháp kỹ thuật lâm sinh thích hợp để xây dựng quần xã thực vật có năng suất và tính ổn định cao. Thông qua mật độ của từng cấp kính có thể biết được rừng đang ở trạng thái nào, dự báo xu hướng phát triển. Phân bố N/D1.3 được mô phỏng bằng phân bố giảm dạng hàm Meyer, phân bố khoảng cách và phân bố Weibull với hai và ba tham số, kết quả đã lựa chọn được dạng phân bố phù hợp đó là phân bố Weibull ba tham số. Kết quả nắn phân bố N/D1.3 theo phân bố Weibull ba tham số được tổng hợp tại bảng 4.8 dưới đây.
Bảng 4.8. Kết quả mô phỏng phân bố N/D1.3 theo hàm Weibull ba tham số của trạng thái rừng IIIA1, IIIA2 và IIIB tại khu vực nghiên cứu Trạng thái OTC β α µ D p- value α luận Kết IIIA1 1 1.8 18.2 6.8 0.082 0.08 0.05 H0+ 2 1.4 12.7 9.7 0.07 0.374 0.05 H0+ IIIA2 3 1.3 19.9 9.8 0.043 0.904 0.05 H0+ 4 1.3 18.4 10.1 0.04 0.989 0.05 H0+ IIIB 5 1.4 19.5 10.5 0.059 0.348 0.05 H0+ 6 1.5 29.5 4.3 0.108 0.004 0.05 H0– Kết quả mô phỏng phân bố N/D1.3 bằng hàm Weibull ba tham số ở Bảng 4.8 cho thấy, có 5/6 ô đo đếm có giá trị p-value > 0,05, nghĩa là phân bố Weibull ba tham số mô phỏng tốt cho phân bố N/D cho 5/6 OTC.
Kết quả nghiên cứu này cũng phù hợp với kết quả của nhiều tác giả nghiên cứu về phân bố N/D lâm phần rừng tự nhiên như nghiên cứu phân bố N/D của rừng tự nhiên lá rộng thường xanh ở Kon Hà Nừng (Lê Sáu, 1996), ở Hương Sơn (Đào Công Khanh, 1996) hay rừng tự nhiên trạng thái IIIA tại huyện An Lão, tỉnh Bình Định (Phạm Quý Vân, 2018) đều chỉ ra rằng hàm Weibull mô phỏng tốt nhất phân bố N/D. Cao Thị Thu Hiền và Nguyễn Hồng Hải (2018), Nguyễn Quang Phúc (2019) [53] cũng kết luận phân bố Weibull mô phỏng tốt cho phân bố N/D.
Sự phù hợp giữa phân bố lý thuyết và phân bố thực nghiệm theo hàm Weibull ba tham số được minh họa trong Hình 4.1 sau đây.
Kết quả Hình 4.1 cho thấy, hình dạng phân bố N/D1.3 ở các OTC khá giống nhau và số lượng cây đạt cực đại tại cỡ kính tại cỡ đường kính 12 cm hoặc cỡ đường kính 16 cm) và giảm dần khi cỡ đường kính tăng lên, điều này cho thấy phân bố N/D1.3 tại khu vực nghiên cứu theo quy luật phân bố giảm đặc trưng cho rừng tự nhiên hỗn loài khác tuổi. Cây có đường kính trên 70 cm có ở OTC 3, OTC 4, OTC 5, OTC 6.
Các cây gỗ lớn (đường kính ngang ngực ≥ 70 cm) đóng vai trò quan trọng trong lưu trữ carbon trong các khu rừng nhiệt đới và có mối liên quan nhiều với điều kiện thời tiết và khí hậu (Clark và Clark 1996). Tuy nhiên, có rất ít cây gỗ lớn trong các khu rừng nhiệt đới ở Đông Nam Á. Trong nghiên cứu này cho thấy mật độ cây gỗ lớn dao động từ 2 cây/ha đến 14 cây/ha chiếm tỷ lệ phần trăm rất thấp (< 3%). Tỷ lệ cây lớn này thấp hơn so với các khu rừng nhiệt đới khác như ở rừng nhiệt đới Tanzania, cây lớn chiếm 4,5% tổng số thân cây (Huang và cộng sự, 2003).
0 5 10 15 20 25 12 20 28 36 44 52 60 68 76 84 92 ft, flt D1.3
OTC 4 Trạng thái IIIA2
ft 0 5 10 15 20 25 30 35 40 12 16 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80 84 88 92 96 1 0 0 ft, flt D1.3
OTC 5 Trạng thái IIIB
ft flt 0 5 10 15 20 25 12 20 28 36 44 52 60 68 76 84 92 ft, flt D1.3
OTC 3 Trạng thái IIIA2
ft flt 0 10 20 30 40 50 60 8 12 16 20 24 28 32 36 40 44 48 52 56 60 64 68 ft, flt D1.3
OTC 2 Trạng thái IIIA1
ft flt 0 10 20 30 40 50 60 70 8 12 16 20 24 28 32 36 40 44 48 52 56 60 ft, flt D1.3
OTC 1 Trạng thái IIIA1
ft flt
Hình 4.1. Phân bố N/D1.3 của đối tượng nghiên cứu theo hàm Weibull ba tham số