Xay dựng mô hình phù hợp mô tả sinh trưởng thông 3 lá ở Lâm đồng

102 946 9
Xay dựng mô hình phù hợp mô tả sinh trưởng thông 3 lá ở Lâm đồng

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

Đề tài “Xác định những mô hình phù hợp nhất để mô tả quá trình sinh trưởng của rừng trồng Thông ba lá (Pinus keysia Royle ex Gordon)” được thực hiện tại Đức Trọng tỉnh Lâm Đồng. Thời gian nghiên cứu từ tháng 4 đến tháng 9 năm 2011. Mục tiêu của đề tài là xây dựng những mô hình thống kê phù hợp nhất để mô tả và phân tích quy luật sinh trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá trên 5 cấp đất khác nhau ở khu vực Đức Trọng tỉnh Lâm Đồng.Kết quả nghiên cứu đã chỉ ra rằng:(1) Nếu sử dụng 5 hàm (Korf, Schumacher, Gompertz, KorsunStrand và DrakinVuevski) để làm phù hợp với số liệu thực nghiệm về quá trình sinh trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá trên 5 cấp đất ở khu vực Đức Trọng tỉnh Lâm Đồng, thì hàm Korf là hàm phù hợp nhất.(2) Nếu sử dụng những mô hình khác nhau, thì đặc trưng sinh trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá trên 5 cấp đất cũng sẽ được báo cáo khác nhau.(3) Cùng một mô hình sinh trưởng, nếu sử dụng những phương pháp và tiêu chuẩn khác nhau để kiểm định mô hình phù hợp, thì đặc trưng sinh trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá trên 5 cấp đất cũng sẽ được báo cáo khác nhau.(4) Sinh trưởng đường kính thân cây Thông ba lá trong 30 năm đầu thay đổi tùy theo tuổi và cấp đất. Thời điểm xuất hiện ZDmax trên cả 5 cấp đất tại cấp tuổi 8. Thời điểm xuất hiện Dmax trên cấp đất IIV tại tuổi 16, còn cấp đất V tại cấp tuổi 18; trung bình tại tuổi 16. Tốc độ sinh trưởng đường kính trong giai đoạn 30 tuổi trên cấp đất I lớn hơn cấp đất II, III, IV và V tương ứng 1,07; 1,14; 1,22 và 1,32 lần.(5) Sinh trưởng chiều cao thân cây Thông ba lá trong 30 năm đầu thay đổi tùy theo tuổi và cấp đất. Thời điểm xuất hiện ZHmax trên cấp đất I, II và III tại cấp tuổi 4, còn cấp đất IV và V tại cấp tuổi 6; trung bình là cấp tuổi 4. Thời điểm xuất hiện Hmax trên cấp đất I, II và III tại cấp tuổi 8, còn cấp đất IV và V tại cấp tuổi 10; trung bình tại cấp tuổi 8. Tốc độ sinh trưởng chiều cao trong giai đoạn 30 tuổi trên cấp đất I lớn hơn cấp đất II, III, IV và V tương ứng 1,05; 1,10; 1,17 và 1,23 lần.(6) Thể tích thân cây Thông ba lá trong 30 năm đầu thay đổi tùy theo tuổi và cấp đất. Tốc độ sinh trưởng thể tích thân cây trong giai đoạn 30 tuổi trên cấp đất I lớn hơn cấp đất II, III, IV và V tương ứng 1,14; 1,28; 1,44 và 1,61 lần.Từ kết quả nghiên cứu, đề tài đã đề xuất những mô hình phù hợp để dự đoán sinh trưởng của Thông ba lá trên 5 cấp đất ở khu vực Đức Trọng tỉnh Lâm Đồng.Những thông tin về sinh trưởng và tăng trưởng của cây cá thể và quần thụ Thông ba lá (Pinus keysia Royle ex Gordon) đóng vai trò rất quan trọng đối với nghiên cứu khoa học và xây dựng những phương thức lâm sinh. Những thông tin này có thể thu thập và phân tích dựa trên những số liệu thực nghiệm. Tuy vậy, vì những số liệu thực nghiệm thường rất nhiều và sắp xếp lộn xộn, nên việc khai thác những thông tin hữu ích từ chúng thường gặp rất nhiều khó khăn, đặc biệt là khó xác định quy luật biến đổi của chúng. Mặt khác, từ những số liệu thực nghiệm, nhà lâm học cũng rất khó đưa ra những dự đoán chiều hướng biến đổi của rừng trong tương lai. Ngoài ra, nếu hàng năm cần phải thu thập những dữ liệu của nhiều quần thể Thông ba lá phân bố trên không gian rộng lớn, thì những tốn kém về người, vật tư và tài chính sẽ tăng lên rất lớn. Để khắc phục những khó khăn nêu trên, đồng thời vẫn nhận được những thông tin có giá trị, người ta có thể áp dụng toán học để giải quyết những vấn đề đặt ra. Bằng cách phân tích những đặc trưng của dãy số liệu thực nghiệm, người ta có thể xây dựng những mô hình thống kê để biểu thị quy luật biến đổi của các nhân tố điều tra trên cây cá thể và lâm phần. Sau đó, phân tích mô hình để tìm ra những thông tin có giá trị. Trước đây đã có một số công trình nghiên cứu về sinh trưởng của rừng tự nhiên và rừng trồng Thông ba lá; trong đó đáng kể là những nghiên cứu của Nguyễn Ngọc Lung (1988, 1999), Phó Đức Đỉnh (1995) và Lê Hồng Phúc (1995). Những nghiên cứu này tập trung chủ yếu vào rừng Thông ba lá tự nhiên, còn rừng trồng vẫn chưa được nghiên cứu đầy đủ. Hiện nay diện tích rừng trồng Thông ba lá ở Lâm Đồng đã được phát triển trên 70.000 ha; trong đó nhiều lâm phần đã có tuổi trên 30 năm. Tuy vậy, hiện nay ngành lâm nghiệp Lâm Đồng vẫn còn thiếu rất nhiều thông tin về đặc tính lâm học của rừng trồng Thông ba lá; trong đó bao gồm cả những thông tin về quy luật biến đổi đường kính, chiều cao và thể tích thân cây, đặc biệt là rừng trồng Thông ba lá có tuổi lớn hơn 30 năm. Vì thế, những nghiên cứu về quy luật biến đổi đường kính, chiều cao và thể tích thân cây Thông ba lá vẫn cần được đặt ra. Mặt khác, khi mô tả và phân tích quy luật biến đổi đường kính, chiều cao và thể tích thân cây Thông ba lá bằng các hàm phi tuyến (Korf, Gompertz và Schumacher), nhiều tác giả đã xác định các tham số những mô hình này bằng đại số ma trận. Về mặt thống kê, những tham số của hàm phi tuyến tính không chỉ được xác định theo phương pháp bình phương nhỏ nhất, mà còn cả phương pháp hồi quy tương quan phi tuyến. Mặt khác, cả hai phương pháp này cũng có thể giải theo những cách khác nhau. Nhận thấy rằng, chính vì chọn lựa những mô hình thống kê khác nhau và cách xác định những tham số của những mô hình cũng khác nhau, nên kết quả báo cáo về quy luật sinh trưởng của cây cá thể và lâm phần cũng khác nhau.Với mong muốn góp phần cung cấp thêm những thông tin hữu ích cho việc phân tích đặc điểm lâm học của rừng trồng Thông ba lá ở Lâm Đồng, đề tài “Xác định những mô hình phù hợp nhất để mô tả quá trình sinh trưởng của rừng trồng Thông ba lá (Pinus keysia Royle ex Gordon) ở khu vực Đức Trọng tỉnh Lâm Đồng” đã được đặt ra. 2. Mục tiêu nghiên cứu2.1. Mục tiêu chungXây dựng những mô hình thống kê phù hợp nhất để mô tả và phân tích quy luật sinh trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá trên 5 cấp đất khác nhau ở khu vực Đức Trọng tỉnh Lâm Đồng.2.2. Mục tiêu cụ thểĐể đạt được mục tiêu chung, đề tài đặt ra 4 mục tiêu cụ thể sau đây:(1)Xác định mô hình thống kê phù hợp nhất để mô tả quy luật sinh trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá trên 5 cấp đất khác nhau.(2)Định lượng những đặc trưng sinh trưởng và tăng trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá trên 5 cấp đất khác nhau.(3)Phân tích so sánh sự khác biệt về những đặc trưng sinh trưởng của Thông ba lá do ảnh hưởng của việc chọn mô hình thống kê khác nhau.(4)Phân tích so sánh sự khác biệt về những đặc trưng sinh trưởng của Thông ba lá do ảnh hưởng của việc chọn phương pháp ước lượng các tham số của mô hình sinh trưởng.3. Phạm vi nghiên cứuPhạm vi nghiên cứu của đề tài là rừng trồng Thông ba lá trong giai đoạn 30 tuổi thuộc 5 cấp khác nhau. Địa điểm nghiên cứu tại Đức Trọng tỉnh Lâm Đồng. Nội dung nghiên cứu chỉ tập trung vào chọn lựa những mô hình phù hợp nhất để mô tả quy luật biến đổi đường kính, chiều cao và thể tích thân cây Thông ba lá. Từ kết quả nghiên cứu, đề xuất một số mô hình thống kê phù hợp nhất để mô tả quy luật sinh trưởng của rừng Thông ba lá.4. Ý nghĩa của đề tàiNhững kết quả nghiên cứu của đề tài đưa lại những ý nghĩa sau đây:(1) Về lý luận, đề tài cung cấp những mô hình phù hợp để phân tích quy luật sinh trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá trên 5 cấp đất khác nhau. (2) Về thực tiễn, những kết quả nghiên cứu của đề tài không chỉ là căn cứ khoa học để mô tả, phân tích và dự đoán quá trình sinh trưởng của Thông ba lá, mà còn cung cấp những thông tin để xây dựng những chỉ tiêu kỹ thuật nuôi rừng, khai thác rừng và phân tích hiệu quả kinh doanh rừng Thông ba lá.

MỞ ĐẦU Đặt vấn đề Những thông tin sinh trưởng tăng trưởng cá thể quần thụ Thông ba (Pinus keysia Royle ex Gordon) đóng vai trò quan trọng nghiên cứu khoa học xây dựng phương thức lâm sinh Những thông tin thu thập phân tích dựa số liệu thực nghiệm Tuy vậy, số liệu thực nghiệm thường nhiều xếp lộn xộn, nên việc khai thác thông tin hữu ích từ chúng thường gặp nhiều khó khăn, đặc biệt khó xác định quy luật biến đổi chúng Mặt khác, từ số liệu thực nghiệm, nhà lâm học khó đưa dự đoán chiều hướng biến đổi rừng tương lai Ngoài ra, hàng năm cần phải thu thập liệu nhiều quần thể Thông ba phân bố không gian rộng lớn, tốn người, vật tư tài tăng lên lớn Để khắc phục khó khăn nêu trên, đồng thời nhận thông tin có giá trị, người ta áp dụng toán học để giải vấn đề đặt Bằng cách phân tích đặc trưng dãy số liệu thực nghiệm, người ta xây dựng mô hình thống kê để biểu thị quy luật biến đổi nhân tố điều tra cá thể lâm phần Sau đó, phân tích mô hình để tìm thông tin có giá trị Trước có số công trình nghiên cứu sinh trưởng rừng tự nhiên rừng trồng Thông ba lá; đáng kể nghiên cứu Nguyễn Ngọc Lung (1988, 1999), Phó Đức Đỉnh (1995) Lê Hồng Phúc (1995) Những nghiên cứu tập trung chủ yếu vào rừng Thông ba tự nhiên, rừng trồng chưa nghiên cứu đầy đủ Hiện diện tích rừng trồng Thông ba Lâm Đồng phát triển 70.000 ha; nhiều lâm phần có tuổi 30 năm Tuy vậy, ngành lâm nghiệp Lâm Đồng thiếu nhiều thông tin đặc tính lâm học rừng trồng Thông ba lá; bao gồm thông tin quy luật biến đổi đường kính, chiều cao thể tích thân cây, đặc biệt rừng trồng Thông ba có tuổi lớn 30 năm Vì thế, nghiên cứu quy luật biến đổi đường kính, chiều cao thể tích thân Thông ba cần đặt Mặt khác, mô tả phân tích quy luật biến đổi đường kính, chiều cao thể tích thân Thông ba hàm phi tuyến (Korf, Gompertz Schumacher), nhiều tác giả xác định tham số mô hình đại số ma trận Về mặt thống kê, tham số hàm phi tuyến tính không xác định theo phương pháp bình phương nhỏ nhất, mà phương pháp hồi quy tương quan phi tuyến Mặt khác, hai phương pháp giải theo cách khác Nhận thấy rằng, chọn lựa mô hình thống kê khác cách xác định tham số mô hình khác nhau, nên kết báo cáo quy luật sinh trưởng cá thể lâm phần khác Với mong muốn góp phần cung cấp thêm thông tin hữu ích cho việc phân tích đặc điểm lâm học rừng trồng Thông ba Lâm Đồng, đề tài “Xác định mô hình phù hợp để mô tả trình sinh trưởng rừng trồng Thông ba (Pinus keysia Royle ex Gordon) khu vực Đức Trọng tỉnh Lâm Đồng” đặt Mục tiêu nghiên cứu 2.1 Mục tiêu chung Xây dựng mô hình thống kê phù hợp để mô tả phân tích quy luật sinh trưởng đường kính, chiều cao thể tích thân Thông ba cấp đất khác khu vực Đức Trọng tỉnh Lâm Đồng 2.2 Mục tiêu cụ thể Để đạt mục tiêu chung, đề tài đặt mục tiêu cụ thể sau đây: (1) Xác định mô hình thống kê phù hợp để mô tả quy luật sinh trưởng đường kính, chiều cao thể tích thân Thông ba cấp đất khác (2) Định lượng đặc trưng sinh trưởng tăng trưởng đường kính, chiều cao thể tích thân Thông ba cấp đất khác (3) Phân tích so sánh khác biệt đặc trưng sinh trưởng Thông ba ảnh hưởng việc chọn mô hình thống kê khác (4) Phân tích so sánh khác biệt đặc trưng sinh trưởng Thông ba ảnh hưởng việc chọn phương pháp ước lượng tham số mô hình sinh trưởng Phạm vi nghiên cứu Phạm vi nghiên cứu đề tài rừng trồng Thông ba giai đoạn 30 tuổi thuộc cấp khác Địa điểm nghiên cứu Đức Trọng tỉnh Lâm Đồng Nội dung nghiên cứu tập trung vào chọn lựa mô hình phù hợp để mô tả quy luật biến đổi đường kính, chiều cao thể tích thân Thông ba Từ kết nghiên cứu, đề xuất số mô hình thống kê phù hợp để mô tả quy luật sinh trưởng rừng Thông ba Ý nghĩa đề tài Những kết nghiên cứu đề tài đưa lại ý nghĩa sau đây: (1) Về lý luận, đề tài cung cấp mô hình phù hợp để phân tích quy luật sinh trưởng đường kính, chiều cao thể tích thân Thông ba cấp đất khác (2) Về thực tiễn, kết nghiên cứu đề tài không khoa học để mô tả, phân tích dự đoán trình sinh trưởng Thông ba lá, mà cung cấp thông tin để xây dựng tiêu kỹ thuật nuôi rừng, khai thác rừng phân tích hiệu kinh doanh rừng Thông ba Chương TỔNG QUAN 1.1 Khái quát rừng Thông ba Theo Lê Mộng Chân Vũ Văn Dũng (1992), Thông ba phân bố tự nhiên Lâm Đồng, Đắc Lắc, Gia Lai Kontum Tổng diện tích rừng Thông ba khoảng 150.000ha Do phân bố vùng cao, nên rừng Thông ba có tác dụng phòng hộ đầu nguồn cho sông lớn, bảo vệ hồ thuỷ lợi, thuỷ điện; đồng thời tạo danh lam thắng cảnh phục vụ nghỉ dưỡng Năm 1988, dựa sở điều kiện tự nhiên sinh trưởng Thông ba lá, Nguyễn Ngọc Lung (1988), phân chia rừng Thông ba thành vùng – phạm vi thích hợp với độ cao từ 1.000 – 1.800 m; phạm vi mở rộng với độ cao 1.000 m 1.800 m phạm vi giới hạn với độ cao 600 m Phạm vi thích hợp trung tâm phân bố rừng Thông ba loài; Thông ba chiếm ưu tầng ưu sinh thái Phạm vi mở rộng vùng ngoại vi độ cao 1.000 m 1.800 m Khi phân bố độ cao 1.000 m, Thông ba mọc hỗn giao với Thông (Pinus merkusii) dầu trà beng (Dipterocapus obtusiforlius) Nếu phân bố độ cao 1.800 m, Thông ba mọc chung với loài thuộc họ Fagaceae thường hỗn giao theo đám Ở phạm vi giới hạn (dưới 600 m), Thông ba sinh trưởng Từ trước đến có số nghiên cứu thông Lâm đồng; đáng kể nghiên cứu ảnh hưởng khí hậu điều kiện lập địa đến sinh trưởng Thông ba Khi nghiên cứu ảnh hưởng điều kiện ngoại cảnh tới sinh trưởng Thông ba Lâm đồng, Nguyễn Ngọc Lung (1988, 1999) nhận thấy rằng, khác biệt đáng kể sinh trưởng chiều cao hai vùng sinh thái khác Đà lạt Bảo lộc, lượng tăng trưởng đường kính có khác mùa mưa mùa khô 1.2 Một số phương pháp mô hình hóa trình sinh trưởng rừng Việt Nam 1.2.1 Tình hình chung Trong lâm học, người ta phân biệt hai khái niệm sinh trưởng tăng trưởng quần thể gỗ hay lâm phần Sinh trưởng tăng lên kích thước gỗ theo tuổi Tăng trưởng cá thể chiều cao (H, m), đường kính thân (D1.3, cm), thể tích thân (V, m 3), tăng trưởng quần thụ trữ lượng (M, m3) theo tuổi biểu thị biến đổi lượng chúng sau năm, nghĩa chúng có đạo hàm bậc theo thời gian Vì sinh trưởng rừng lâm phần thay đổi theo thời gian, nên người ta gọi trình sinh trưởng gỗ lâm phần (Dẫn theo Nguyễn Văn Thêm, 2002) Lý thuyết lâm học rằng, sinh trưởng gỗ lâm phần xác định tập hợp nhiều yếu tố khí hậu (ánh sáng, nhiệt độ, mưa, gió…) phi khí hậu (địa hình, đất, sinh vật người) Theo Nguyễn Văn Thêm (2002) cách tổng quát xem sinh trưởng rừng lâm phần hàm số khí hậu, địa hình, đất, sinh vật, người thời gian, nghĩa sinh trưởng = f(khí hậu, địa hình, đất, sinh vật, người thời gian) Chính thế, nhiều nhà lâm học xây dựng mô hình toán để biểu thị trình sinh trưởng cá thể lâm phần theo thời gian Sau mô hình vận dụng để mô tả dự đoán nhân tố điều tra thân thời gian Để xây dựng mô hình biểu thị trình biến đổi theo tuổi nhân tố điều tra thân gỗ (đường kính, chiều cao, thể tích), người ta phân tích hình thái đường cong biến đổi đường kính, chiều cao thể tích thân cá thể toàn trình đời sống (Vũ Tiến Hinh ctv, 1992; Vũ Tiến Hinh, 2005; Chế Đình Lý, 1997; Nguyễn Hải Tuất, 1998; Nguyến Hải Tuất Nguyễn Trọng Bình, 2005) Nói chung, đường cong sinh trưởng đường kính, chiều cao thể tích thân toàn trình đời sống có dạng sigmoid hay hình chữ S Tuy vậy, theo Belov (1983)(Dẫn theo Nguyễn Văn Thêm, 2002), tùy theo giai đoạn sống lâm phần, đường cong sinh trưởng đường kính, chiều cao thể tích thân dạng đường thẳng, dạng hàm số mũ dạng bậc Nếu xét trình đời sống cây, đường cong sinh trưởng có dạng sigmoid Chính thế, nhiều nhà lâm học giới xây dựng mô hình toán khác để biểu thị trình biến đổi theo tuổi nhân tố điều tra thân gỗ (đường kính, chiều cao, thể tích) Năm 1999, Nguyễn Ngọc Lung áp dụng nhiều mô hình khác để mô tả trình sinh trưởng loài gỗ mọc nhanh Việt Nam, đặc biệt Thông ba Sau phân tích mô hình, Nguyễn Ngọc Lung (1999) cho mô hình Gompertz mô hình Schumacher hai mô hình phù hợp để mô tả trình sinh trưởng loài gỗ mọc nhanh Việt Nam Những mô hình nhiều tác giả sử dụng để phân tích quy luật sinh trưởng nhiều loài gỗ khác (Vũ Tiến Hinh, 2005; Bảo Huy, 1995; Lê Thị Hà, 2003; Chế Đình Lý, 1997; Vũ Nhâm, 1988; Lê Hồng Phúc, 1995; Phạm Trọng Thịnh, 2006) Theo Vũ Nhâm (1988), mô hình Korf vận dụng để lập biểu cấp đất cho rừng thông đuôi ngựa kinh doanh gỗ mỏ khu Đông Bắc Việt Nam 1.2.2 Những phương pháp áp dụng cho rừng Thông ba Từ trước đến có số công trình nghiên cứu sinh trưởng rừng Thông ba cách xây dựng mô hình thống kê Nguyễn Ngọc Lung (1988; 1999) công bố nghiên cứu sinh trưởng, suất sản lượng rừng Thông ba Việt Nam Sau phần kết nghiên cứu sử dụng để lập biểu dùng kinh doanh điều tra rừng Thông ba lá; chẳng hạn biểu cấp chiều cao, biểu thể tích, biểu cấp đất, biểu trình sinh trưởng… Để xây dựng loại biểu này, tác giả vận dụng nhiều mô hình toán khác để mô tả trình sinh trưởng, xây dựng biểu cấp đất biểu trình sinh trưởng rừng Thông ba Những mô hình áp dụng nhiều mô hình Gompertz (1925), Schumacher (1939), Drakin – Vuevski (1940), Kosun-Strand (1964) Korf (1973) Cũng theo cách thức trên, Viên Ngọc Hùng (1989), vận dụng mô hình Schumacher để lập biểu cấp đất Thông ba Lâm Đồng Lê Hồng Phúc (1995), sử dụng mô hình mũ biến đổi để mô tả trình biến đổi sinh khối phận mặt đất rừng Thông ba 1.2.3 Thảo luận chung Sau tổng quan số tài liệu có liên quan đến việc vận dụng mô hình toán để mô tả phân tích quy luật biến đổi đường kính, chiều cao thể tích thân theo tuổi, tác giả nhận thấy cần thảo luận thêm số vấn đề sau đây: (1) Cho đến nhiều nhà lâm học nước ta vận dụng mô hình toán để mô tả phân tích quy luật biến đổi đường kính, chiều cao thể tích thân theo tuổi Thông qua nhiều nghiên cứu, người ta nhận thấy rằng, tùy theo giai đoạn sống gỗ lâm phần, hình thái đường cong sinh trưởng khác Vì thế, việc vận dụng mô hình để mô tả trình biến đổi theo tuổi nhân tố điều tra thân gỗ (đường kính, chiều cao, thể tích) tùy theo loài giai đoạn tuổi Nói chung, người ta cho rằng, gỗ mọc nhanh Việt Nam, trình biến đổi theo tuổi nhân tố điều tra thân gỗ (đường kính, chiều cao, thể tích) mô tả hàm phi tuyến; mô hình Gompertz (1925), Schumacher (1939), Drakin – Vuevski (1940), Kosun-Strand (1964) Korf (1973) vận dụng nhiều Vì thế, xác định mô hình phù hợp để mô tả trình sinh trưởng Thông ba Lâm Đồng, tác giả vận dụng mô hình nói (2) Khi xây dựng mô hình phi tuyến, hầu hết tác giả xác định tham số mô hình cách chuyển chúng dạng tuyến tính ước lượng bình phương nhỏ xác định đại số ma trận Giải pháp bình phương nhỏ có đặc tính ước lượng phương sai không trệch nhỏ Tuy vậy, mô hình phi tuyến tuyến tính hóa phương pháp đại số sai số lệch chuẩn với phương sai số Để khắc phục nhược điểm hồi quy tuyến tính, người ta xác định tham số mô hình phi tuyến cách lặp lại nhiều lần xác định tham số mô hình cho “Tổng bình phương sai lệch không đổi = SSR” Theo đó, trước hết chương trình bắt đầu với ước lượng ban đầu tham số mô hình Từ tham số ban đầu này, chương trình tính toán SSR (tổng sai lệch bình phương từ giá trị quan sát so với đường cong) Sau chương trình lại điều chỉnh tham số mô hình để tạo đường cong bám sát điểm thực nghiệm Thủ tục lặp lại nhiều lần xác định tham số mô hình cho tổng sai lệch bình phương không đổi Cuối chương trình tóm tắt kết phù hợp Độ xác tham số phụ thuộc phần vào giá trị ước lượng ban đầu tiêu chuẩn dừng (thường SSR min) Điều có nghĩa là, với chuỗi số liệu nhau, giả định tham số ban đầu khác nhau, nhận tham số khác mô hình Phân tích cho thấy, mô hình thống kê dùng để mô tả dãy số liệu, giải phương pháp khác nhau, kết nhận mô hình với tham số khác Vì thế, để nhận mô hình phù hợp nhất, nhà nghiên cứu cần phải giả định tham số ban đầu hợp lý dùng tiêu chuẩn SSRmin để kiểm tra (3) Tác giả nhận thấy rằng, số trường hợp, nhiều lý khác nhau, nhà lâm học chọn hàm phi tuyến tính định để làm phù hợp với số liệu D-A, H-A V-A thực nghiệm Ngoài ra, tham số hàm phi tuyến tính xác định theo hai phương pháp – phương pháp bình phương nhỏ phương pháp hồi quy tương quan phi tuyến Mặt khác, hai phương pháp giải theo cách khác Chẳng hạn, phương pháp bình phương nhỏ nhất, người ta cố định tham số, tham số khác thay đổi Nhận thấy rằng, chọn lựa mô hình thống kê khác cách xác định tham số mô hình khác nhau, nên kết báo cáo quy luật sinh trưởng cá thể lâm phần khác Vì thế, đề tài hướng vào phân tích khác biệt kết báo cáo trình sinh trưởng đường kính, chiều cao thể tích thân Thông ba ảnh hưởng việc chọn hàm số phương pháp xác định tham số chúng (4) Tác giả nhận thấy rằng, có nhiều thông tin rừng Thông ba Tây Nguyên nói chung Lâm Đồng nói riêng, áp dụng cho khu vực cụ thể nhiều có sai lệch lớn Vì thế, để có thông tin hữu ích rừng Thông ba khu vực định, nhận thấy nghiên cứu cần phải đặt Chương NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Đối tượng nghiên cứu 2.1.1 Điều kiện tự nhiên khu vực nghiên cứu 2.1.1.1 Vị trí địa lý Khu vực nghiên cứu đề tài nằm khu vực huyện Đức Trọng có toạ độ địa lý sau: +Từ 1080 15'30" đến 108025’ độ kinh đông +Từ 11035’đến 110 37’ 30” độ vĩ bắc Ranh giới: + Phía đông giáp, huyện Ninh Sơn, tỉnh Ninh Thuận + Phía tây giáp, huyện Lâm Hà, tỉnh Lâm Đồng + Phía nam giáp, huyện Bắc Bình tỉnh Bình Thuận + Phía Bắc giáp, thành phố Đà Lạt, tỉnh Lâm Đồng 2.1.1.2 Địa hình Huyện Đức Trọng nằm vùng cao nguyên Di Linh, địa hình chia cắt mạnh, độ dốc cục lớn Độ cao tuyệt đối 850 m đến 1050 m, giảm dần từ nam xuống bắc Độ dốc binh quân 100 - 300 2.1.1.3 Đất Trong khu vực nghiên cứu có nhóm đất đai sau: + Nhóm đất phù sa: Được hình thành mẫu chất bồi đắp sông Đa Nhim, Đạ Queyon suối xung quanh khu vực, thành phần giới thịt trung bình có màu xám đen, nhóm đất thích hợp nhiều loại trồng như: bắp, rau, đậu đỗ, mía, dâu, ăn trái… 10 Hàm (1) Phương trình (2) Korf H(II) = 54,22566*exp(-5,84275*A^-0,59363) Gompertz H(II) = 25,72672*exp(-2,71205*exp(-0,12462*A)) Kosun-Strand H(II) = A^2/(1,68887 + 0,41114*A + 0,02496*A^2) Drakin-Vuevxki H(II) = 28,32621*(1-exp(-0,07772*A))^1,39019 Công thức (3) (3.36) (3.38) (3.39) (3.40) Bảng 3.72 Đặc trưng sinh trưởng H Thông ba cấp đất I Hàm Lượng tăng trưởng hàng năm ZHmax (m) A(năm) H(m) 1,80 5,1 1,90 3,8 Korf Gompert z Korsun Drakin Thực tế 1,64 1,54 1,63 4 10 Lượng tăng trưởng bình quân năm ∆Hmax (m) A(năm) H(m) 1,45 11,6 1,90 3,8 5,2 5,5 12,3 1,43 1,43 1,32 8 14 11,4 11,4 18,5 Bảng 3.73 Đặc trưng sinh trưởng H Thông ba cấp đất II HI (cm) Hàm Korf Gompert z Korsun Drakin Thực tế HI (cm) ZH ∆H (m/năm) DI (cm) Lượng tăng trưởng hàng năm ZHmax (m) A(năm) H(m) 1,52 4,2 1,55 3,1 Hàm Korf 1,42 1,33 1,55 6 A(năm) 7,1 7,2 7,3 ZH ∆H (m/năm) ZH ∆H (m/năm) Lượng tăng trưởng bình quân năm ∆Hmax (m) A(năm) H(m) 1,24 9,9 1,55 3,1 1,22Hàm Korsun-Strand 1,21 1,24 HI (cm) A(năm) ZH ∆H (m/năm) Hàm Gompertz Hàm Drakin-Vuevski A(năm) A(năm) 88 sinh trưởng chiều cao thân Hình 3.52 Đồ thị so sánh trình Thông ba cấp đất I làm phù hợp với hàm (Korf, KorsunStrand, Drakin-Vuevski Gompertz) 9,7 9,7 9,9 ZH ∆H (m/năm) HII (m) HII (m) Hàm Korf ZH ∆H (m/năm) Hàm Korsun-Strand A(năm) HII (cm) ZH ∆H (m/năm) A(năm) HII (m) ZH ∆H (m/năm) ) Hàm Gompertz Hàm Drakin-Vuevski A(năm) A(năm) Hình 3.53 Đồ thị so sánh trình 89 sinh trưởng chiều cao thân Thông ba cấp đất II làm phù hợp với hàm (Korf, Korsun-Strand, Drakin-Vuevski Gompertz) Từ so sánh đặc trưng sinh trưởng D H Thông ba làm phù hợp với hàm (Korf, Korsun-Strand, Drakin-Vuevski Gompertz), rút nhận định sau đây: (a) Khi áp dụng hàm (Korf, Korsun-Strand, Drakin-Vuevski Gompertz) để mô tả trình sinh trưởng Thông ba lá, kết nhận khác (b) Để chọn mô hình phù hợp với số liệu thực tế, cần phải xác định rõ tiêu chuẩn đánh giá phù hợp mô hình hay tiêu chuẩn dừng Việc chọn tiêu chuẩn dừng khác dẫn đến kết luận khác Vì thế, để chọn mô hình phù hợp với số liệu thực tế, bên cạnh tiêu chuẩn kiểm định tính phù hợp mô hình, cần phải lưu ý đến hai đặc trưng quan trọng - thời điểm đạt lượng tăng trưởng thường xuyên hàng năm lớn lượng tăng trưởng bình quân năm lớn 3.5.2 Sự khác hai phương pháp ước lượng hệ số mô hình sinh trưởng Về lý thuyết, hệ số mô hình sinh trưởng xác định theo hai phương pháp – phương pháp bình phương nhỏ phương pháp hồi quy phi tuyến tính Mặt khác, mô hình phù hợp chọn dựa tiêu chuẩn kiểm định khác hay tiêu chuẩn dừng khác Vì thế, phần so sánh khác biệt kết dự đoán trình sinh trưởng Thông ba ảnh hưởng việc chọn phương pháp ước lượng tham số mô hình sinh trưởng tiêu chuẩn dừng Để làm rõ vấn đề đặt đây, bảng 3.74 3.75 dẫn kết so sánh trình sinh trưởng D Thông ba cấp đất I mô tả hàm Korf; ba tham số hàm Korf ước 90 lượng theo phương pháp bình phương nhỏ hồi quy phi tuyến tính Ngoài ra, mô hình Korf phù hợp chọn dựa theo hai tiêu chuẩn kiểm định – R2max) SSRmin Bảng 3.74 Phân tích hồi quy tương quan D-A cấp đất I mô hình Korf m B c R (2) (3) (4) (5) 40 26,99064 1,42473 -0,9362 60 11,45252 0,87839 -0,9798 80 9,45706 0,71198 -0,9856 100 8,73825 0,62397 -0,9875 120 8,39579 0,56778 -0,9883 140 8,21112 0,52807 -0,9889 160 8,10540 0,49815 -0,989 180 8,04395 0,47457 -0,9891 200 8,00911 0,45538 -0,9892 Phi tuyến 121,63213 8,36199 0,56370 Bảng 3.75 Những thống kê sai lệch mô hình Korf Bình phương nhỏ Phương pháp (1) Bình phương nhỏ Phương pháp (1) Phi tuyến m ±Se SSR(*) MAE (cm) MAPE (%) (2) (3) (4) (5) (6) 40 0,3159 3,6929 2,25 16,18 60 0,1058 0,4142 1,60 9,13 80 0,0722 0,1928 1,57 8,15 100 0,0588 0,1278 1,59 8,15 120 0,0516 0,0985 1,59 8,14 140 0,0471 0,0820 1,59 8,14 160 0,0439 0,0715 1,60 8,18 180 0,0416 0,0642 1,60 8,25 200 0,0399 0,0588 1,60 8,31 1,50 7,77 (*) (*) Ghi chú: SSR = [ln(ln(m/Y))tn - ln(ln(m/Y))lt]^2 R^2 (6) 87,65 96,01 97,14 97,52 97,69 97,77 97,82 97,84 97,85 96,34 SSR (7) 283,66 167,49 149,97 145,14 144,07 144,42 145,36 146,55 147,83 143,61 Kết phân tích hồi quy tương quan D với A Thông ba cấp đất I theo hai phương pháp bình phương nhỏ hồi quy phi tuyến tính ghi lại Bảng 3.74 3.75 91 Đối với phương pháp bình phương nhỏ nhất, để ước lượng ba tham số m, b c hàm Korf, trước hết cố định tham số m, sau tìm tham số b c Phân tích số liệu Bảng 3.74 cho thấy, giả định m tăng dần từ 40 đến 200, R tăng dần từ 87,65% đến 97,85% Theo khuynh hướng biến đổi R nhận thấy, m tăng lên 200, R tiếp tục tăng lên tiệm cận giá trị 98% Để xác định tham số m, b c tương ứng với R 98%, cần phải tiếp tục thay đổi tham số m Những tính toán cho thấy, dù m tăng lên lớn R đạt tới giá trị 98% Điều cho thấy, sử dụng R 2max tiêu chuẩn dừng, xác định xác tham số mô hình Korf Tuy vậy, xem giá trị 97,85% xấp xỉ 98%, tham số m mô hình Korf giả định 200 Theo đó, sử dụng R 2max tiêu chuẩn dừng, trình sinh trưởng đường kính thân Thông ba cấp đất I mô tả hàm Korf có dạng: D(I) = 200*exp(-8,00911*A^-0,45538) (3.84) R2 = 97,85%; Se = ±0,0399; MAE = 1,60 cm; MAPE = 8,31% Ngược lại, sử dụng SSRmin tiêu chuẩn dừng, trình sinh trưởng đường kính thân Thông ba cấp đất I mô tả hàm Korf có dạng: D(I) = 120*exp(-8,39579*A^-0,56778) (3.85) R2 = 97,69%; Se = ±0,0516; MAE = 1,59 cm; MAPE = 8,14% Nếu sử dụng phương pháp hồi quy tương quan phi tuyến để ước lượng ba tham số m, b c hàm Korf, trình sinh trưởng đường kính thân Thông ba cấp đất I mô tả hàm Korf có dạng: D(I) = 121,63213*exp(-8,36199*A^-0,56370) (3.86) R2 = 96,34%; MAE = 1,50 cm; MAPE = 7,77% Như vậy, sử dụng R2max tiêu chuẩn dừng, mô hình 3.84 mô hình phù hợp nhất, mô hình 3.86 mô hình phù hợp Ngược lại, sử dụng SSRmin tiêu chuẩn dừng, mô hình 3.86 mô hình phù hợp nhất, mô hình 3.84 mô hình phù hợp 92 Theo tác giả, để thấy rõ mô hình phù hợp nhất, cần phải so sánh đặc trưng sinh trưởng đường kính thân Thông ba khảo sát từ ba hàm 3.84, 3.85 3.86 với số liệu thực nghiệm Điều đặt vì, mục đích chọn hàm phân tích hồi quy tương quan xác định mô hình phù hợp với số liệu thực nghiệm Theo đó, khảo sát ba hàm 3.84, 3.85 3.86 cho thấy (Bảng 3.76 - 3.79; Hình 3.54-3.56), hai hàm 3.85 3.86 đặc trưng ZD max tuổi đạt ZDmax, ∆Dmax tuổi đạt ∆Dmax phù hợp với số liệu thực tế Trái lại, sử dụng hàm 3.84 để mô tả trình sinh trưởng D Thông ba cấp đất I, đặc trưng ZDmax tuổi đạt ZDmax, ∆Dmax tuổi đạt ∆Dmax có sai khác lớn so với số liệu thực tế Bảng 3.76 Những đặc trưng tăng trưởng đường kính thân Thông ba khảo sát hàm khác Bình phương nhỏ Phương pháp (1) Phi tuyến Thực tế m (2) 40 60 80 100 120 140 160 180 200 ZDmax (cm) (3) 2,51 1,95 1,80 1,72 1,67 1,64 1,61 1,59 1,58 1,68 1,63 A(năm) (4) 8 8 8 8 8 ∆Dmax (5) 1,52 1,39 1,35 1,33 1,32 1,31 1,30 1,30 1,30 1,34 1,32 A(năm) (6) 12 14 14 14 16 16 16 16 16 16 16 Bảng 3.77 Quá trình sinh trưởng đường kính thân Thông ba cấp đất I khảo sát hàm Korf với tiêu chuẩn dừng R2max D(I) = 200*exp(-8,00911*A^-0,45538) 93 (3.84) A(năm) (1) 10 12 14 16 18 20 22 24 26 28 30 D(cm) (2) 0,6 2,8 5,8 8,9 12,1 15,1 18,0 20,7 23,4 25,8 28,2 30,4 32,5 34,5 36,5 DI(cm) ZD(cm/năm) ΔD(cm/năm) (3) (4) 0,29 0,29 1,12 0,71 1,48 0,97 1,58 1,12 1,57 1,21 1,51 1,26 1,45 1,29 1,37 1,30 1,30 1,30 1,24 1,29 1,17 1,28 1,11 1,27 1,06 1,25 1,01 1,23 0,96 1,22 Pd% (5) 50,0 39,7 25,6 17,6 13,0 10,0 8,0 6,6 5,6 4,8 4,2 3,7 3,3 2,9 2,6 Kd (6) 0,206 0,488 0,647 0,741 0,800 0,839 0,868 0,888 0,904 0,917 0,927 0,935 0,942 0,947 ZD ΔD (cm/năm) Hình 3.54 Đồ thị mô tả trình sinh trưởng đường kính thân Thông ba cấp đất I Mô hình dự đoán hàm Korf; hệ số xác định theo hồi quy tương quan tuyến tính Tiêu chuẩn dừng R2max A (năm) Bảng 3.78 Quá trình sinh trưởng đường kính thân Thông ba cấp đất I khảo sát hàm Korf với tiêu chuẩn dừng SSRmin D(I) = 120*exp(-8,39579*A^-0,56778) A(năm) (1) D(cm) (2) 0,4 2,6 5,8 9,1 ZD(cm/năm) ΔD(cm/năm) (3) (4) 0,21 0,21 1,11 0,66 1,57 0,96 1,67 1,14 94 (3.85) Pd% (5) 50,0 42,1 27,2 18,4 Kd (6) 0,158 0,456 0,633 12,4 15,5 18,4 21,1 23,6 25,9 28,1 30,1 32,0 33,8 35,5 10 12 14 16 18 20 22 24 26 28 30 DI(cm) 1,64 1,55 1,45 1,35 1,25 1,17 1,09 1,02 0,95 0,90 0,84 1,24 1,29 1,31 1,32 1,31 1,30 1,28 1,26 1,23 1,21 1,18 13,2 10,0 7,9 6,4 5,3 4,5 3,9 3,4 3,0 2,6 2,4 0,736 0,800 0,842 0,872 0,894 0,910 0,922 0,932 0,940 0,947 0,953 ZD ΔD (cm/năm) Hình 3.55 Đồ thị mô tả trình sinh trưởng đường kính thân Thông ba cấp đất I Mô hình dự đoán hàm Korf; hệ số xác định theo hồi quy tương quan tuyến tính Tiêu chuẩn dừng SSRmin A (năm) Bảng 3.79 Quá trình sinh trưởng đường kính thân Thông ba cấp đất I mô tả hàm Korf phù hợp D(I) = 121,63213*exp(-8,36199*A^-0,56370) A(năm) D(cm) ZD(cm/năm) ΔD(cm/năm) (1) (2) (3) 0,4 (3.86) Pd% Kd (4) (5) (6) 0,21 0,21 50,0 - 2,6 1,11 0,66 42,0 0,160 5,8 1,57 0,96 27,1 0,457 9,1 1,67 1,14 18,3 0,634 10 12,4 1,63 1,24 13,2 0,736 12 15,5 1,55 1,29 10,0 0,800 95 14 18,4 1,45 1,31 7,9 0,842 16 21,1 1,35 1,32 6,4 0,872 18 23,6 1,26 1,31 5,3 0,894 20 25,9 1,17 1,30 4,5 0,910 22 28,1 1,09 1,28 3,9 0,922 24 30,2 1,02 1,26 3,4 0,932 26 32,1 0,96 1,23 3,0 0,940 28 33,9 0,90 1,21 2,7 0,947 30 35,6 0,85 1,19 2,4 0,952 DI(cm) ZD ΔD (cm/năm) Hình 3.56 Đồ thị mô tả trình sinh trưởng đường kính thân Thông ba cấp đất I mô hình Korf Các hệ số hàm Korf xác định theo hồi quy tương quan phi tuyến tính Tiêu chuẩn dừng SSRmin A (năm) Từ kết nghiên cứu mục 3.5.1 3.5.2, đến nhận định chung sau: (1) Khi áp dụng hàm phi tuyến tính (Korf, Korsun-Strand, Drakin-Vuevski Gompertz) để mô tả trình sinh trưởng đường kính thân, chiều cao thể tích thân Thông ba giai đoạn 30 năm đầu, kết nhận khác Vì thế, để chọn hàm số phù hợp với số liệu thực nghiệm, cần phải kiểm định tính phù hợp hàm theo tiêu chuẩn SSRmin (2) Nếu sử dụng phương pháp bình phương sai lệch nhỏ để ước lượng tham số hàm phi tuyến tính, trình tính toán phải thực qua nhiều bước Việc tính toán dừng lại thỏa mãn điều kiện SSRmin 96 (4) So với phương pháp bình phương sai lệch nhỏ nhất, phương pháp hồi quy phi tuyến tính cho phép ước lượng tham số hàm phi tuyến tính với sai lệch nhỏ (5) Nếu chọn phương pháp ước lượng tham số hàm phi tuyến tính tiêu chuẩn dừng khác nhau, mô hình ước lượng phù hợp khác Vì thế, mô tả dự đoán trình sinh trưởng cá thể lâm phần hàm (Korf, Korsun-Strand, Drakin-Vuevski Gompertz), cần phải rõ phương pháp ước lượng tham số tiêu chuẩn dừng KẾT LUẬN VÀ ĐỀ NGHỊ 1.Kết luận Từ kết nghiên cứu, đến kết luận sau đây: (1) Trong số hàm (Korf, Schumacher, Gompertz, Korsun-Strand Drakin-Vuevski), hàm Korf hàm phù hợp để mô tả trình sinh trưởng đường kính, chiều cao thể tích thân Thông ba cấp đất khu vực Đức Trọng tỉnh Lâm Đồng (2) Sinh trưởng đường kính thân Thông ba 30 năm đầu thay đổi tùy theo tuổi cấp đất Thời điểm xuất ZD max cấp đất cấp tuổi Thời điểm xuất ∆Dmax cấp đất I-IV tuổi 16, cấp đất V cấp tuổi 18; trung bình tuổi 16 Suất tăng trưởng đường kính Pd (%) cấp đất suy giảm nhanh theo tuổi Tốc độ sinh trưởng đường kính giai đoạn 30 tuổi cấp đất I lớn cấp đất II, III, IV V tương ứng 1,07; 1,14; 1,22 1,32 lần (3) Sinh trưởng chiều cao thân Thông ba 30 năm đầu thay đổi tùy theo tuổi cấp đất Thời điểm xuất ZH max cấp đất I, II III cấp tuổi 4, cấp đất IV V cấp tuổi 6; trung bình cấp tuổi Thời điểm xuất 97 ∆Hmax cấp đất I, II III cấp tuổi 8, cấp đất IV V cấp tuổi 10; trung bình cấp tuổi Suất tăng trưởng chiều cao Ph (%) cấp đất suy giảm nhanh theo tuổi Tốc độ sinh trưởng chiều cao giai đoạn 30 tuổi cấp đất I lớn cấp đất II, III, IV V tương ứng 1,05; 1,10; 1,17 1,23 lần (4) Thể tích thân Thông ba 30 năm đầu thay đổi tùy theo tuổi cấp đất Suất tăng trưởng thể tích thân cấp đất suy giảm nhanh theo tuổi Tốc độ sinh trưởng thể tích thân giai đoạn 30 tuổi cấp đất I lớn cấp đất II, III, IV V tương ứng 1,14; 1,28; 1,44 1,61 lần (5) Nếu sử dụng hàm (Korf, Schumacher, Gompertz, KorsunStrand Drakin-Vuevski) để mô tả trình sinh trưởng đường kính, chiều cao thể tích thân Thông ba cấp đất, kết nhận khác Vì thế, để biết hàm phù hợp nhất, cần phải thực kiểm định mặt thống kê; tiêu chuẩn kiểm định thích hợp tổng bình phương sai lệch nhỏ Đề nghị Đề tài áp dụng hàm (Korf, Schumacher, Gompertz, Korsun-Strand Drakin-Vuevski) để làm phù hợp với trình sinh trưởng đường kính, chiều cao thể tích thân Thông ba giai đoạn 30 năm đầu Vì thế, tác giả kiến nghị quan tâm đến việc xây dựng mô hình sinh trưởng đường kính, chiều cao thể tích thân Thông ba lá, cần tiếp tục thử nghiệm tính phù hợp hàm lâm phần Thông ba có tuổi lớn 30 năm Ngoài ra, bên cạnh việc áp dụng hàm (Korf, Schumacher, Gompertz, Korsun-Strand Drakin-Vuevski), tìm kiếm khả ứng dụng hàm số mũ khác hàm Sigmoid (S-curve), Assmann-Franz, Terasaki 98 TÀI LIỆU THAM KHẢO Bảo Huy, 1995 Dự đoán sản lượng rừng tếch Đắc Lắc Tạp chí lâm nghiệp, Hà Nội Bùi Việt Hải, 2002 Thống kê lâm nghiệp Tủ sách Trường đại học nông lâm Tp Hồ Chí Minh 150 trang Chế Đình Lý, 1997 Mô hình hóa Lâm nghiệp Tài liệu dành cho học viên cao học ngành Lâm nghiệp, Đại học Quốc gia Tp Hồ Chí Minh 147 trang Kleinbaun D.G; Kupper L.L 1998 Applied Regression Analysis and Other Multivariable Methods Brooks/Cole Publishing Company 800 Page Lê Hồng Phúc, 1995 Nghiên cứu sinh khối rừng thông ba (Pinus keysia Royle ex Gordon) Đà Lạt Lâm Đồng Tạp chí lâm nghiệp Số 9/1995 Lê Mộng Chân, Vũ Văn Dũng, 1992 Thực vật và thực vật đặc sản rừng, Trường Đại học Lâm nghiệp 180 trang Lê Thị Hà, 2003 Đánh giá khả ứng dụng phương trình sinh trưởng vào mô tả dự đoán sinh trưởng cho số loài trồng nước ta Luận Văn thạc sỹ, Đại học lâm nghiệp Việt Nam 80 trang Nguyễn Ngọc Lung, 1988 Nghiên cứu sở khoa học kỹ thuật để kinh doanh tổng hợp rừng thông ba Tây Nguyên Viện khoa học lâm nghiệp Việt Nam 70 trang 99 Nguyễn Ngọc Lung, 1999 Nghiên cứu tăng trưởng sản lượng rừng trồng áp dụng cho rừng thông ba (Pinus keysia Royle ex Gordon) Việt Nam Nhà xuất Nông Nghiệp, Hà Nội 250 trang 10 Nguyễn Hải Tuất, 1982 Thống kê toán học lâm nghiệp Nhà xuất Nông Nghiệp, Hà Nội 324 trang 11 Nguyễn Hải Tuất Nguyễn Trọng Bình, 2005 Khai thác sử dụng SPSS để xử lý số liệu nghiên cứu lâm nghiệp, Nhà xuất Nông Nghiệp, Hà Nội 12 Nguyễn Văn Thêm, 2002 Sinh thái rừng Nhà xuất Nông Nghiệp Chi nhánh Tp Hồ Chí Minh 374 trang 13 Nguyễn Văn Thêm, 2003 Phản ứng Thông ba (Pinus keysia Royle ex Gordon) yếu tố khí hậu Lạc Dương – Lâm Đồng Tạp chí nông nghiệp phát triển nông thôn, Hà Nội Số 2/2003 145 trang 14 Nguyễn Văn Thêm, 2000 Hướng dẫn sử dụng Statgraphics 3.0 & 5.1 để xử lý thông tin lâm học Nhà xuất Nông Nghiệp Chi nhánh Tp Hồ Chí Minh 179 trang 15 Nguyễn Văn Thêm, 2004 Lâm sinh học Nhà xuất Nông Nghiệp, Hà Nội 179 trang 16 Ngô Đình Quế, 1983 Đất rừng thông ba Lâm Đồng, Trong sách “Một số kết nghiên cứu khoa học kỹ thuật lâm nghiệp 1976-1985, Viện khoa học lâm nghiệp Nhà xuất Nông nghiệp, Hà Nội 17 Phạm Trong Thịnh, 2006 Nghiên cứu số yếu tố làm sở cho việc tỉa thưa rừng đước (Rhizophora apiculata Blume) trồng vùng ven biển Nam Bộ, tóm tắt luận án tiến sỹ nông nghiệp 24 trang 18 Phạm Trọng Nhân, 2001 Ảnh hưởng nhân tố khí hậu đến sinh trưởng thông ba (Pinus keysia Royle ex Gordon) Đà Lạt tỉnh Lâm Đồng Luận Văn thạc sỹ, Trường đại học nông lâm nghiệp Tp Hồ Chí Minh 75 trang 100 19 Phó Đức Đỉnh (1995), Nghiên cứu tái sinh tự nhiên Thông ba Lâm Đồng, tóm tắt luận án tiến sĩ nông nghiệp 24 trang 20 Thái Văn Trừng (1998), Những hệ sinh thái rừng nhiệt đới Việt Nam, Nhà xuất Khoa học Kỹ thuật, Hà Nội 412 trang 21 Tiêu chuẩn ngành 04-TCN-66-2003 (2003), Biểu điều tra kinh doanh rừng trồng 14 loài chủ yếu Nhà xuất Nông Nghiệp, Hà Nội 182 trang 22 Vancley J.K., 1989 A growth model for north Queensland rainforests For Ecol Manage 27:245-271 23 Vancley J.K., 1991 Aggregating tree species to developing diamater incement equations for tropical Ecol Manage 42:143-168 24 Vancley J.K., 1994 Modelling forest growth and yield CAB International, Wallingford UK 304 Page 25 Viện điều tra quy hoạch rừng, 1995 Số tay điều tra rừng Nhà xuất Nông Nghiệp, Hà Nội 252 trang 26 Viên Ngọc Hùng, 1989 Nghiên cứu xây dựng biểu cấp đất thông ba Lâm Đồng, Trong sách “Một số kết nghiên cứu khoa học kỹ thuật lâm nghiệp 1976-1985”, Viện khoa học lâm nghiệp Nhà xuất Nông nghiệp, Hà Nội 27 Võ Văn Hồng, Trần Văn Hùng (2006), Tăng trưởng rừng, Cẩm nang lâm nghiệp, Cục phát triển Lâm Nghiệp Việt Nam 250 trang 28 Vũ Tiến Hinh tác giả khác 1992 Điều tra rừng, Trường đại học lâm nghiệp, Hà Nội 203 trang 29 Vũ Tiến Hinh, 2005 Sản lượng rừng Nhà xuất Nông Nghiệp, Hà Nội 212 trang 30 Vũ Nhâm, 1988 Lập biểu cấp đất cho rừng thông đuôi ngựa kinh doanh gỗ mỏ khu Đông Bắc, Tạp chí lâm nghiệp, Hà Nội Số 3/1988 175 trang 101 102 [...]... tăng trưởng D và H của Thông ba lá để chỉ ra sự khác nhau giữa 5 hàm sinh trưởng 2.3.3.3 So sánh sự khác nhau giữa hai phương pháp ước lượng các hệ số của mô hình sinh trưởng Thực chất của phần này là trả lời câu hỏi: Các hệ số của hàm sinh trưởng được xác định bằng phương pháp bình phương nhỏ nhất và phương pháp hồi quy tương quan phi tuyến có dẫn đến những báo cáo khác nhau về những đặc trưng sinh trưởng. .. Drakin-Vuevski) với hai dãy số liệu D-A và H-A của Thông ba lá Cách thức xác định các hệ số của 5 hàm sinh trưởng được thực hiện như bước 2 của mục 2.3.3.1 Tiếp theo, khảo sát 5 mô hình sinh trưởng (Korf, Schumacher, Gompertz, Korsun-Strand và Drakin-Vuevski) để xác định những đặc trưng sinh trưởng và tăng trưởng D và H thân cây Những chỉ tiêu dùng để so sánh 5 hàm sinh trưởng là 18 ZYmax và tuổi đạt ZYmax,... nghiên cứu bắt đầu từ tháng 05 năm 2011 đến tháng 09 năm 2011 2.2 Nội dung nghiên cứu (1) Mô hình sinh trưởng đường kính thân cây Thông ba lá (2) Mô hình sinh trưởng chiều cao thân cây Thông ba lá (3) Mô hình sinh trưởng thể tích thân cây Thông ba lá (4) Đặc điểm sinh trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá (5) Sự khác nhau giữa những mô hình sinh trưởng 2.3 Phương pháp nghiên cứu 2.3.1... (2.13) Sự khác biệt về khuynh hướng sinh trưởng D, H và V thân cây Thông ba lá trên 5 cấp đất khác nhau được kiểm định thông qua so sánh điểm chặn và độ dốc của các mô hình bằng thống kê F Khi điểm chặn và độ dốc của các mô hình khác nhau (P < 0,05 hoặc 0,01), thì cấp đất có ảnh hưởng đến sinh trưởng D, H và V thân cây Thông ba lá Ngược lại, khi điểm chặn và độ dốc của các mô hình không có khác biệt... (11) suất tăng trưởng chiều cao 13 thân cây (Ph%), (12) suất tăng trưởng thể tích thân cây (PV%), (13) nhịp điệu sinh trưởng đường kính thân cây (Kd), (14) nhịp điệu sinh trưởng chiều cao thân cây (Kh), (15) nhịp điệu sinh trưởng thể tích thân cây (Kv) + Đối với cấp đất, đề tài xác định cấp đất của rừng Thông ba lá theo “Tiêu chuẩn ngành 04-TCN-66-2003”, bao gồm 5 cấp từ cấp I (tốt nhất) đến cấp V (xấu... tích và thuyết minh Bước 3 Xác định những đặc trưng sinh trưởng đường kính, chiều cao và thể tích thân cây Thông ba lá Nội dung tính toán ở đây nhằm trả lời hai câu hỏi chính sau đây: (a) Sinh trưởng D, H và V thân cây Thông ba lá trong toàn bộ khu vực nghiên cứu và trên 5 cấp đất có những đặc điểm gì? (b) Cấp đất khác nhau có ảnh hưởng đến sinh trưởng D, H và V thân cây Thông ba lá như thế nào? Để... chuyển từ giai đoạn sinh trưởng nhanh sang giai đoạn sinh trưởng chậm Để trả lời câu hỏi (b), trình tự xử lý số liệu như sau: + Trước hết, tập hợp dãy số liệu D - A, H - A và V – A của những cây giải tích đại diện cho những lâm phần mọc trên 5 cấp đất khác nhau Kế đến, mã hóa dãy số liệu D - A, H - A và V – A từ cấp đất I đến cấp đất V bằng những trị số nguyên tương ứng từ 1 đến 5 + Kế đến, chuyển dãy... làm rõ vấn đề đặt ra, đề tài chỉ sử dụng mô hình sinh trưởng đường kính thân cây Theo đó, trình tự xử lý số liệu như sau: Trước hết, sử dụng hai phương pháp bình phương nhỏ nhất và hồi quy tương quan phi tuyến để ước lượng các tham số của hàm sinh trưởng D phù hợp nhất Ở đây hàm sinh trưởng D phù hợp nhất đã được xác định ở bước 2 của mục 2.3.3.1 Kế đến, sử dụng hai tiêu chuẩn R2max và SSRmin để chọn... phương pháp ước lượng những tham số của hàm sinh trưởng D Tiếp theo, khảo sát hai mô hình sinh trưởng D phù hợp nhất để xác định những đặc trưng ZDmax và tuổi đạt ZDmax, ∆Dmax và tuổi đạt ∆Dmax Cuối cùng, phân tích so sánh những đặc trưng ZD max và tuổi đạt ZDmax, ∆Dmax và tuổi đạt ∆Dmax để chỉ ra sự khác nhau giữa hai phương pháp ước lượng các hệ số của mô hình sinh trưởng 2.4.3 Công cụ xử lý số liệu... Trọng Bình, 2005; Kleinbaun và Kupper, 1998; Vancley, 1989; 1991 và Vancley, 1994) Chương 3 20 KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 3.1 Mô hình sinh trưởng đường kính của Thông ba lá 3.1.1 Mô hình sinh trưởng đường kính của Thông ba lá trên cấp đất I Để làm rõ quá trình sinh trưởng đường kính thân cây Thông ba lá trên cấp đất từ I, đã làm phù hợp 5 hàm (Korf, Schumacher, Gompertz, Kosun-Strand và Drakin-Vuevski)

Ngày đăng: 09/08/2016, 12:56

Từ khóa liên quan

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan