Nghiên cứu này sử dụng bộ dữ liệu thu thập tại Vườn Quốc gia Bến En trong thời gian từ ngày 1/8/2010 đến ngày 3/11/2010, đó là bộ dữ liệu “Rừng nguyên sinh”. Dưới đây là những mô tả sơ bộ về dữ liệu này.
Bộ dữ liệu “Rừng nguyên sinh” được thu thập tại các khu vực rừng tự nhiên. Trong bộ dữ liệu này, chất lượng sinh trưởng của cây rừng được đánh giá theo phương pháp chuyên gia và được chia ở 3 mức khác nhau là: sinh trưởng tốt, sinh trưởng trung bình và sinh trưởng kém. Các tiêu chí điều tra với cây rừng bao gồm: đường kính 1m3; đường kính tán; chiều cao vút ngọn; độ cao ô tiêu chuẩn so với mực nước biển; mức độ quý hiếm của cây; kiểu rừng; ví trí của ơ tiêu chuẩn và tên cây. Mục tiêu của nghiên cứu là xây dựng mơ hình dự báo với biến phụ thuộc là chất lượng sinh trưởng của cây được mã hóa là biến tính trạng có thứ tự. Các biến độc lập bao gồm các số đo đường kính 1m3 (Dk1m3); đường kính tán (Dktan); chiều cao vút ngọn (Hvn); độ cao ô tiêu chuẩn so với mực nước biển(Docao); mức độ quý hiếm của cây (EN, CR, LR, V, VU và các cây không thuộc loại quý hiếm); kiểu rừng (IIb, IIIa1, IIIa2, nửa gỗ và một số cây không phân loại); và tên cây (Cúc, Hoa Hồng, Ngọc Lan, Sổ, Hoa Môi, Sau Sau, Thù Du và nhóm cây chưa được phân lớp) sau khi đã được mã hóa lại một cách thích hợp.
Trong số các biến độc lập kể trên, 4 biến : Dk1m3, Dktan, Hvn, Docao là các biến định lượng sẽ được đưa vào mơ hình một cách trực tiếp.
Đối với biến định tính là “mức độ quý hiếm” của cây, ta lập thành 3 biến giả “Nguycap1” ứng với giá trị EN; “Nguycap2” ứng với giá trị LR và CR; “Nguycap3” ứng với giá trị V và VU. Nhóm cây khơng thuộc loại q hiếm được lấy làm nhóm chứng của biến này.
Đối với biến định tính “kiểu rừng”, ta thành lập các biến giả “KieuIIb”; “KieuIIIa1”; “KieuIIIa2”; “KieuNuago”. Nhóm chứng cho biến này là nhóm cây khơng được phân loại.
Biến định tính “Tên cây” được lập ra các biến giả “LopCuc”; “LopHHong”; “LopHMoi”; “LopSSau”; “LopSo”; “LopThuDu”; “LopNgLan” và “LopKhac”. Lớp “Sổ” được lấy làm nhóm chứng cho biến này.
Dữ liệu được mơ tả đối với “rừng ngun sinh” có thể tóm tắt qua các bảng từ bảng 4.1 tới bảng 4.4:
Bảng 4.1: Phân loại chất lượng sinh trưởng.
Chất lượng sinh trưởng (theo đánh giá của chuyên gia) Phân loại Tần số Tỉ lệ %
Tốt 2496 65.6
Trung bình 1015 26.7
Kém 287 7.5
Không phân loại 4 0.1
Tổng cộng 3802 100
Dữ liệu thu được có tổng cộng 3802 cây rừng được đo đạc kích thước và đánh giá chất lượng sinh trưởng. Trong đó, có 2496 cây được các chuyên gia đánh giá có chất lượng sinh trưởng tốt, 1015 cây có chất lượng sinh trưởng trung bình, 287 cây được đánh giá phát triển kém và có 4 cây chưa được đánh giá về chất lượng sinh trưởng. Vì 4 cây này chiếm tỷ lệ rất nhỏ nên ta điều chỉnh số liệu này bằng cách gán 4 cây đó cho nhóm cây có chất lượng sinh trưởng trung bình. Bảng 4.2 tóm tắt lại việc phân lớp cây rừng trên. Như vậy, cây rừng thuộc phân lớp Hoa Hồng và phân lớp Sổ chiếm tỷ lệ lớn tương ứng 37.5% và 34.6%, tức là, các cây thuộc 2 lớp trên
là khá phổ biến tại khu vực điều tra. Những cây thuộc phân lớp Cúc và phân lớp Thù Du được bắt gặp với số lượng ít hơn lần lượt là 44 cây và 106 cây. Số cây chưa xác định được phân lớp là 115 cây .
Bảng 4.1: Phân loại chất lượng sinh trưởng.
Phân lớp cây rừng Phân loại Tần số Tỉ lệ % LopKhac 115 3.0 LopCuc 44 1.2 LopHHong 1426 37.5 LopHMoi 202 5.3 LopNgLan 456 12 LopSSau 136 3.6 LopSo 1317 34.6 LopThuDu 106 34.6
Các tiêu chí được đo đạc đối với mỗi cây là đường kính 1m3, chiều cao vút ngọn và đường kính tán.. Tuy nhiên, trong số liệu này có khoảng hơn 200 cây khơng được xác định đường kính 1m3. Các chỉ tiêu về chiều cao vút ngọn và đường kín tán thì đều được biết. Với biến độ cao, những ô tiêu chuẩn chưa được xác định độ cao thì độ cao của ơ đó được phục hồi bằng cách gán nó với độ cao trung bình của các ơ đã biết. Mô tả chi tiết các biến này được thể hiện trong bảng 4.3:
Bảng 4.3: Mô tả biến định lượng
Tên Biến Số quan sát
Minimum Maximum Trung bình Độ lệch chuẩn Dk1m3 (cm) 3696 5.0 180.0 25.052 22.2359 Hvn (cm) 3802 2.5 43.0 10.865 4.1774 Dktan (m) 3802 50 15.0 3.6443 1.92210 Docao (m) 3802 31.00 219.00 85.9324 36.45929
Trong số liệu trên, các cây rừng cũng được đánh giá mức độ bị đe dọa theo tiêu chuẩn của Liên minh Bảo tồn Thiên nhiên Quốc tế và tiêu chuẩn Việt Nam 2007.
Để xem xét ảnh hưởng của lớp cây này tới việc đánh giá chất lượng sinh trưởng của các chuyên gia, ta đưa vào 3 biến Nguycap1; Nguycap2 và Nguycap3. Số lượng các cây được cho chuyên gia cho là nguy cấp không nhiều. Phần lớn các cây điều tra được đều không bị đe dọa. Thông tin về các biến này được thể hiện ở bảng 4.4 sau:
Bảng 4.4: Mức độ nguy cấp của cây rừng
Tên Biến Tần số Tần suất
Nguycap1 579 15.2% Nguycap2 139 3.7% Nguycap3 49 1.3% Khongnguycap 3035 79.8% Tổng số 3802 100%