ĐẶT VẤN ĐỀ Rừng có vai trò như lá phổi xanh của con người, sự tồn tại của con người không tách khỏi môi trường sống mà rừng là một phần của môi trường sống đó. Trước tình hình như vậy, không chỉ Việt Nam mà trên toàn Thế Giới đã kí công ước bảo tồn đa dạng sinh học thông qua đại hội thượng đỉnh ở Ri O.de.Janero (1992). Theo công ước này bảo tồn đa dạng sinh học là cái mốc đánh dấu cam kết của các quốc gia trên toàn thế giới về bảo vệ đa dạng sinh học, bảo đảm sử dụng một cách bền vững nguồn tài nguyên sinh vật và nguồn lợi thu được phải chia công bằng. Chính vì vậy mà Liên Hợp Quốc đã kêu gọi các nước có rừng hãy bảo vệ rừng khi còn chưa muộn, đặc biệt là các nước đang phát triển đó chính là chiến lược toàn cầu. Việt Nam được coi là trung tâm đa dạng sinh học của vùng Đông Nam Châu Á, điều này thể hiện ở số lượng loài với 12000 loài thực vật, 273 loài thú, 180 loài bò sát, 773 loài chim, 80 loài ếch nhái và hàng ngàn loài động vật không xương sống và còn nhiều loài chưa được phát hiện. Trong sự nghiệp bảo vệ đa dạng sinh học ở các quốc gia, các Vườn quốc gia (VQG), các Khu bảo tồn (KBT) giữ một vai trò quan trọng. Muốn thực hiện được nhiệm vụ chính của mình điều đầu tiên là phải đánh giá được tính đa dạng sinh học một cách đầy đủ, đó là cơ sở khoa học để thực hiện nhiệm vụ như: Bảo tồn các loài quý hiếm, các loài nguy cấp, các nguồn gen hay các hệ sinh thái khác. Tuy nhiên, cũng như các VQG và KBTTN khác ở Việt Nam, VQG Pù Mát cũng đã và đang phải đối mặt với sự suy giảm nghiêm trọng về ĐDSH bởi các nguyên nhân trực tiếp và gián tiếp. Mặc dù vậy, cho tới nay VQG Pù Mát chưa có các biện pháp thật sự hữu hiệu để quản lý, bảo vệ một cách có hiệu quả hệ thực vật thuộc VQG. Xuất phát từ nhận thức đó, để góp phần vào sự nghiệp bảo vệ đa dạng sinh học và cấu trúc rừng, làm cơ sở cho việc quản lý rừng hiệu quả hơn, việc thực hiện đề tài: “Đặc trưng cấu trúc và đa dạng loài thực vật tại Phân khu phục hồi sinh thái – Vườn Quốc Gia Pù Mát – tỉnh Nghệ An” là hết sức cần thiết góp phần bổ sung thêm về lý thuyết sinh thái học rừng tự nhiên trên nhằm đề xuất một số giải pháp chăm sóc và nuôi dưỡng rừng tự nhiên phục vụ yêu cầu bảo tồn sinh học có hiệu quả.
ĐẶT VẤN ĐỀ Rừng có vai trị phổi xanh người, tồn người không tách khỏi môi trường sống mà rừng phần mơi trường sống Trước tình vậy, khơng Việt Nam mà tồn Thế Giới kí cơng ước bảo tồn đa dạng sinh học thông qua đại hội thượng đỉnh Ri O.de.Janero (1992) Theo công ước bảo tồn đa dạng sinh học mốc đánh dấu cam kết quốc gia toàn giới bảo vệ đa dạng sinh học, bảo đảm sử dụng cách bền vững nguồn tài nguyên sinh vật nguồn lợi thu phải chia cơng Chính mà Liên Hợp Quốc kêu gọi nước có rừng bảo vệ rừng chưa muộn, đặc biệt nước phát triển chiến lược tồn cầu Việt Nam coi trung tâm đa dạng sinh học vùng Đông Nam Châu Á, điều thể số lượng loài với 12000 loài thực vật, 273 lồi thú, 180 lồi bị sát, 773 lồi chim, 80 loài ếch nhái hàng ngàn loài động vật khơng xương sống cịn nhiều lồi chưa phát Trong nghiệp bảo vệ đa dạng sinh học quốc gia, Vườn quốc gia (VQG), Khu bảo tồn (KBT) giữ vai trò quan trọng Muốn thực nhiệm vụ điều phải đánh giá tính đa dạng sinh học cách đầy đủ, sở khoa học để thực nhiệm vụ như: Bảo tồn loài quý hiếm, loài nguy cấp, nguồn gen hay hệ sinh thái khác Tuy nhiên, VQG KBTTN khác Việt Nam, VQG Pù Mát phải đối mặt với suy giảm nghiêm trọng ĐDSH nguyên nhân trực tiếp gián tiếp Mặc dù vậy, VQG Pù Mát chưa có biện pháp thật hữu hiệu để quản lý, bảo vệ cách có hiệu hệ thực vật thuộc VQG Xuất phát từ nhận thức đó, để góp phần vào nghiệp bảo vệ đa dạng sinh học cấu trúc rừng, làm sở cho việc quản lý rừng hiệu hơn, việc thực đề tài: “Đặc trưng cấu trúc đa dạng loài thực vật Phân khu phục hồi sinh thái – Vườn Quốc Gia Pù Mát – tỉnh Nghệ An” cần thiết góp phần bổ sung thêm lý thuyết sinh thái học rừng tự nhiên nhằm đề xuất số giải pháp chăm sóc ni dưỡng rừng tự nhiên phục vụ yêu cầu bảo tồn sinh học có hiệu CHƯƠNG TỔNG QUAN VỀ VẤN ĐỀ NGHIÊN CỨU Trên giới Rừng nhiệt đới ẩm nơi mà tính đa dạng sinh học cao nhất, nhờ đa dạng phong phú mà hút nhiều nhà khoa học 1.1 Nghiên cứu cấu trúc rừng 1.1.1 Cấu trúc tổ thành Theo Richard P.W (1952), rừng mưa nhiệt đới, hecta ln có 40 lồi gỗ, có trường hợp cịn 100 lồi Nhiều lồi gỗ lớn sinh trưởng hỗn giao với theo tỷ lệ đồng đều, có có hai lồi chiếm ưu Trên giới có nhiều tác giả nghiên cứu sở sinh thái cấu trúc rừng tiêu biểu Baur G N (1964) [1] E.P Odum (1971) [13] Hai tác giả tập trung vào vấn đề sinh thái nói chung sở sinh thái kinh doanh rừng mưa nhiệt đới nói riêng Qua làm sáng tỏ khái niệm hệ sinh thái rừng, sở để nghiên cứu cấu trúc rừng đứng quan điểm sinh thái học 1.1.2 Nghiên cứu định lượng cấu trúc 1.1.2.1 Nghiên cứu quy luật phân bố số theo cỡ đường kính (N/D1.3) Là quy luật cấu trúc lâm phần nên nhiều nhà khoa học lâm học điều tra rừng nghiên cứu Các cơng trình tiêu biểu phải kể đến là: + Meyer (1934) (dẫn theo Nguyễn Hải Tuất, 1986) [15] mô tả quy luật phân bố N/D1.3 phương trình tốn học có dạng đường cong giảm liên tục gọi phương trình Meyer hay hàm Meyer + Ballell (1973) sử dụng hàm Weibull, Schiffel, Naslund (1936, 1937) xác lập phân bố Charlier cho phân bố N/D 1.3 lâm phần loài tuổi sau khép tán (dẫn theo Phạm Ngọc Giao, 1995) [6] 1.1.2.2 Nghiên cứu quy luật phân bố số theo chiều cao (N/Hvn) Phần lớn tác giả nghiên cứu cấu trúc lâm phần theo chiều thẳng đứng dựa vào phân bố số theo chiều cao Phương pháp áp dụng để nghiên cứu cấu trúc đứng rừng tự nhiên vẽ phẫu đồ đứng với kích thước khác tùy theo mục đích nghiên cứu Các phẫu đồ mang lại hình ảnh khái quát cấu trúc tầng tán, phân bố số theo chiều thắng đứng Từ rút nhận xét đề xuất ứng dụng thực tế Với phương pháp nhiều nhà nghiên cứu ứng dụng như: Richards P.W (1952) [21], Rolllet (1979) 1.1.2.3 Nghiên cứu quy luật tương quan chiều cao đường kính thân (Hvn /D1.3) Qua nghiên cứu nhiều tác giả cho thấy, chiều cao tương ứng với cỡ đường kính cho trước ln tăng theo tuổi, kết tự nhiên sinh trưởng Trong cỡ xác định, tuổi khác nhau, rừng thuộc cấp sinh trưởng khác nhau, cấp sinh trưởng giảm tuổi lâm phần tăng lên dẫn đến tỷ lệ Hvn/D1.3 tăng theo tuổi Từ đường cong quan hệ Hvn D1.3 thay đổi ln dịch chuyển phía tuổi lâm phần tăng Krauter.G (1958) Tiourin.A.V (1932) (dẫn theo Phạm Ngọc Giao, 1995) [6] nghiên cứu tương quan chiều cao đường kính ngang ngực dựa sở cấp đất cấp tuổi Naslund M (1929), Hohenadl W (1936), Michailov F (1934, 1952), Prodan.M (1944), Meyer.H.A (1952) (dẫn theo Phạm Ngọc Giao, 1995) [6], dùng phương pháp giải tích tốn học đề nghị sử dụng dạng phương trình để mơ tả quan hệ H/D h = a + b1.d + b2.d2 h = a + b1.d + b2.d2 + b3.d3 h - 1.3 = d2/(a + b.d)2 h = a + b.logd h = a + b1.d +b2.logd h = k.db Như vậy, để biểu thị tương quan chiều cao đường kính thân ta sử dụng nhiều dạng phương trình Song việc lựa chọn phương trình để biểu thị mối tương quan H - D1.3 tùy thuộc vào lồi trồng cụ thể 1.2 Nghiên cứu đa dạng sinh học Thuật ngữ đa dạng sinh học thuật ngữ mẻ dùng để tính phong phú sống trái đất hàng triệu loài thực vật, động vật, vi sinh vật, gen chứa đựng lồi hệ sinh thái vơ phức tạp tồn môi trường sống Ý nghĩa việc nghiên cứu đa dạng sinh học: qua việc nghiên cứu đa dạng sinh học loài đa dạng sinh học gen cho thấy loài, nguồn gen quý làm sở cho cơng tác nghiên cứu cải thiện giống lồi, có khả chống chịu với hoàn cảnh bất lợi mở rộng nơi sống loài ngày nâng cao suất, chất lượng chúng Các cơng trình nghiên cứu đa dạng sinh học giới có từ lâu: thực vật chí Đơng Dương (1905 – 1952) quyển, H.humber (1938 -1950) Ở Nga từ 1928 – 1932 xem thời kỳ mở đầu cho nghiên cứu hệ thực vật Cụ thể Talmachay AI cho “Chỉ cần điều tra diện tích đủ lớn để bao trùm phong phú nơi sống khơng có phân hóa mặt địa lý Ơng gọi hệ thực vật cụ thể vùng nhiệt đới ẩm thường 1500 – 2000 loài Ở Việt Nam 2.1 Nghiên cứu phân loại trạng thái rừng Mục đích chủ yếu phân loại trạng thái rừng nhằm xác định đối tượng rừng với đặc trưng cấu trúc cụ thể, từ lựa chọn, đề xuất biện pháp lâm sinh thích hợp để điều khiển, dẫn đắt rừng đạt trạng thái chuẩn Về phân loại rừng trước hết phải kể đến Loetschau (1966) [12] đưa hệ thống phân chia kiểu trạng thái cho kinh doanh rừng hỗn giao thường xanh rộng nhiệt đới Viện Điều tra Quy hoạch rừng dựa hệ thống phân loại Loeschau cải tiến cho phù hợp với đặc điểm rừng tự nhiên Việt Nam áp dụng hệ thống phân loại (QPN – 84) Tiếp theo Thái Văn Trừng (1978) [16] đứng quan điểm sinh thái chia rừng Việt Nam thành 14 kiểu thảm thực vật Đây cơng trình tổng qt, đáp ứng yêu cầu quy luật sinh thái Xuất phát từ tính đa dạng, phong phú rừng nhiệt đới, Thái Văn Trừng đưa kết luận: Không thể dùng quần hợp thực vật làm đơn vị phân loại tác giả kinh điển sử dụng vùng ơn đới Ơng đề xuất dùng kiểu thảm thực vật làm đơn vị phân loại lấy hình thái, cấu trúc quần thể làm tiêu chuẩn phân loại Trần Ngũ Phương (1985 – 1988) [14] đưa phương pháp phân chia rừng nhằm phục vụ công tác điều chế với đơn vị phân chia lô dựa nhân tố nhóm nhân tố sinh thái tự nhiên, giai đoạn phát triển suy thoái rừng, khả tái tạo rừng tái sinh tự nhiên, đặc điểm địa hình thổ nhưỡng Như vậy, tác giả cho rằng: Việc phân loại trạng thái rừng Việt Nam cần thiết công tác nghiên cứu sản xuất kinh doanh Tùy mục tiêu cụ thể mà lựa chọn phương pháp phân loại khác nhau, nhằm làm rõ đặc điểm đối tượng nghiên cứu 2.2 Nghiên cứu cấu trúc rừng 2.2.1 Cấu trúc tổ thành Đây nhân tố có ảnh hưởng định đến cấu trúc sinh thái hình thái rừng Tổ thành rừng tiêu quan trọng dùng để đánh giá mức độ đa dạng sinh học, tính ổn định, tính bền vững hệ sinh thái rừng Cấu trúc tổ thành nhiều nhà khoa học Việt Nam đề cập cơng trình nghiên cứu Bảo Huy (1993) [3], Đào Cơng Khanh (1995) nghiên cứu tổ thành loài rừng tự nhiên Đăklăc Hương Sơn – Hà Tĩnh xác định: Tỷ lệ tổ thành nhóm lồi mục đích, nhóm lồi hỗ trợ nhóm lồi phi mục đích cụ thể, từ đề xuất biện pháp khai thác thích hợp cho đối tượng theo hướng điều chỉnh tổ thành hợp lý 2.2.2 Quy luật phân bố số theo cỡ đường kính (N/D1.3) Thống kê cơng trình nghiên cứu rừng tự nhiên Việt Nam cho thấy: Phân bố N/D1.3 tầng cao (D ≥ 6cm) có hai dạng chính: - Dạng giảm liên tục có nhiều đỉnh phụ hình cưa - Dạng đỉnh chữ J Với dạng cụ thể, tác giả chọn mơ hình tốn học thích hợp để mô Đồng Sỹ Hiền (1974), lập biểu thể tích đứng rừng tự nhiên miền Bắc Việt Nam đưa kết luận: Dạng tổng quát phân bố N/D1.3 phân bố giảm, trình khai thác chọn thơ khơng theo quy tắc nên đường thực nghiệm có dạng hình cưa Với kiểu phân bố thực nghiệm vậy, tác giả dùng hàm Meyer họ đường cong Pearson để mô tả Nguyễn Hải Tuất (1986) [15] sử dụng hàm khoảng cách để mô tả phân bố thực nghiệm dạng đỉnh sát cỡ kính bắt đầu đo 2.2.3 Quy luật phân bố số theo chiều cao (N/Hvn) Những nghiên cứu Đồng Sỹ Hiền (1974) cho thấy: Phân bố số theo chiều cao (N/Hvn) lâm phần tự nhiên hay loài thường có nhiều đỉnh, phản ánh kết cấu phức tạp rừng chặt chọn Thái Văn Trừng (1978), cơng trình nghiên cứu đưa kết nghiên cứu cấu trúc tầng gỗ rừng loại IV 2.2.4 Tương quan chiều cao với đường kính (Hvn - D1.3) Trong điều tra kinh doanh rừng, việc nghiên cứu mối quan hệ có ý nghĩa đặc biệt quan trọng Thông qua tương quan H - D1.3, dựa vào giá trị cỡ kính để suy diễn giá trị chiều cao tương ứng mà không cần thiết đo cao tồn bộ, từ làm sở xác định trữ lượng chung lâm phần, xác định kiểu trạng thái rừng, kết cấu rừng…Từ đó, đề xuất biện pháp kỹ thuật lâm sinh để điều chỉnh kết cấu rừng tiến tới kết cấu rừng ổn định 2.3 Nghiên cứu đa dạng sinh học Nghiên cứu thực vật Việt Nam trước hết phải kể đến cơng trình: thực vật chí Nam Bộ Leureiro, thực vật chí rừng Nam Bộ Pierrel Một cơng trình lớn quy mô giá trị công trình nghiên cứu hệ thực vật Đơng Dương tác giả Pháp, kết “ Thực vật chí đại cương Đơng Dương” bao gồm tập sách có ý nghĩa lớn nhà thực vật học Việt Nam Tiếp theo bổ sung Humbert, đến thực vật chí Lào, Campuchia, Việt Nam xuất từ năm (1960) nước ta có đến 26 tập, sau Poct (1965) dựa thực vật chí đại cương Đơng Dương thống kê 5190 lồi Các tác giả Việt Nam đưa số công trình thảm thực vật tiêu biểu cơng trình lớn: - Thảm thực vật rừng Việt Nam Thái Văn Trừng (1963 – 1978) [16] Tác giả tổng kết cơng bố cơng trình nghiên cứu với 7004 lồi thực vật bậc cao có mạch thuộc 1850 họ Việt Nam ơng nhấn mạnh ưu ngành thực vật hạt kín hệ thực vật Việt Nam với 6336 loài chiếm 99,9%, 17227 chi chiếm 93,5% 239 họ chiếm 62,7% tổng số taxon bậc - Bước đầu nghiên cứu rừng Miền Bắc Việt Nam Trần Ngũ Phương (1970) [14] tiến hành phân loại rừng miền Bắc Việt Nam chia thành đai kiểu Ta thấy cơng trình đánh giá tổng quát cho toàn hệ thực vật Việt Nam đặc biệt “ gỗ rừng Việt Nam” Viện Điều tra Quy hoạch rừng xuất (1971 - 1988) [22], nhóm quan trọng định tồn hệ sinh thái rừng có ý nghĩa quan trọng đời sống người Nghiên cứu đa dạng thực vật VQG Pù Mát Năm 1993 Viện Điều Tra Quy hoạch rừng điều tra nghiên cứu hệ thực vật VQG Pù Mát để xây dựng Dự án thành lập Khu bảo tồn, bước đầu thống kê 986 loài thực vật thuộc 522 chi 153 họ; Đồng thời phân tích tính đa dạng đánh giá nguồn tài nguyên gồm 291 loài gỗ, 220 loài thuốc, 60 loài cảnh, 37 loài cho dầu béo, 96 loài ăn được, 34 làm rau, 30 lồi có chất độc 44 lồi có nguy bị đe dọa Năm 1998 - 2001 Nguyễn Nghĩa Thìn, Ngơ Trực Nhã Nguyễn Thị Hạnh [18] giới thiệu kết nghiên cứu thuốc số xã: Châu Khê, Lục Dạ, Môn Sơn thuộc vùng đệm VQG Pù Mát gồm 512 loài, 325 chi, 115 họ Các kết tác giả tổng kết, cơng bố tạp chí Quảng Tây Côn Minh – Trung Quốc năm 1999 Tạp chí di truyền ứng dụng tháng 4/1999 Năm 1998 Nguyễn Thị Quý điều tra thành phần loài Dương xỉ Khu BTTN Pù Mát lần xác định có 90 lồi, 42 chi, 23 họ phân bố sinh cảnh khác nhau, có 66,7% kinh tế Từ năm 2000 Nguyễn Nghĩa Thìn, Phân viện Điều tra Quy hoạch rừng Bắc Trung Bộ [17] số nhà khoa học khác tiếp tục điều tra hệ thực vật vùng lõi VQG Pù Mát hệ thực vật rừng Pù Mát gồm 1114 loài, 545 chi, 159 họ thuộc ngành mô tả ô tiêu chuẩn Năm 2001 Nguyễn Nghĩa Thìn cộng tiếp tục điều tra thống kê thêm 315 loài thực vật VQG Pù Mát [17] Năm 2004, Nguyễn Nghĩa Thìn, Nguyễn Thanh Nhàn công bố nghiên cứu đa dạng thực vật VQG Pù Mát với 2494 loài [19] CHƯƠNG MỤC TIÊU, ĐỐI TƯỢNG, PHẠM VI , NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Mục tiêu nghiên cứu - Góp phần định lượng số đặc trưng cấu trúc rừng thơng qua việc mơ hình hóa quy luật sinh học - Thử nghiệm số tiêu để đánh giá tính đa dạng lồi thuộc Phân khu phục hồi sinh thái - Đề xuất số giải pháp quản lý bảo vệ rừng để nâng cao khả bảo tồn đa dạng sinh học VQG Pù Mát 2.2 Đối tượng phạm vi nghiên cứu - Đối tượng nghiên cứu: Trạng thái IIIA1, IIB thuộc phân khu phục hồi sinh thái – VQG Pù Mát - Phạm vi nghiên cứu: Phân khu phục hồi – VQG Pù Mát 2.3 Nội dung nghiên cứu 2.3.1 Phân lập trạng thái rừng rừng 2.3.2 Đặc trưng cấu trúc tầng cao + Cấu trúc tổ thành + Phân bố số theo đường kính (N - D1.3) + Phân bố số theo chiều cao (N - Hvn) + Quy luật tương quan chiều cao đường kính (Hvn – D1.3) 2.3.3 Đặc trưng cấu trúc tầng tái sinh + Cấu trúc tổ thành + Phân bố số theo cấp chiều cao chất lượng + Mật độ tái sinh có triển vọng 2.3.4 Đặc trưng tính đa dạng lồi 2.3.4.1 Mức độ phong phú loài 2.3.4.2 Mức độ đa dạng loài - Hàm số liên kết Shannon – Wienr - Chỉ số Simpson - Phương pháp tính đa dạng lý thuyết thơng tin - Chỉ số hợp lý 2.3.5 Đề xuất số giải pháp cho công tác quản lý khu vực nghiên cứu 2.4 Phương pháp nghiên cứu 2.4.1 Phương pháp luận Rừng hệ sinh thái hoàn chỉnh, có mức ổn định tương đối, thành phần cấu thành hệ sinh thái ln có mối quan hệ chặt chẽ tách rời, phản ánh đặc điểm cấu trúc quần thể quần xã tương ứng Trong trình nghiên cứu đề tài tiến hành điều tra, thu thập số liệu đặc trưng rừng hệ thống tiêu chuẩn điển hình tạm thời Sử dụng cơng cụ tốn học để mơ hình hóa quy luật phân bố, tính tốn tiêu đa dạng sinh học nhằm hạn chế tính áp đặt chủ quan người nghiên cứu góp phần phản ánh quy luật chung lâm phần Từ đặc điểm cấu trúc rừng kết hợp với sách có liên quan hoạt động quản lý, sử dụng, xây dựng rừng địa phương, đề xuất số biện pháp quản lý rừng bền vững 2.4.2 Phương pháp điều tra ngoại nghiệp - Phương pháp điều tra sơ thám - Phương pháp điều tra thực địa Dựa vào tình hình thực tế rừng phân loại rừng theo phương pháp Loetschau để chọn đối tượng nghiên cứu + Điều tra tầng cao Mỗi trạng thái đặt OTC hình chữ nhật có diện tích 1000m (25x40m) đo đếm nhân tố điều tra: ◦ Xác định tên loài cho có D1.3 ≥ 6cm đánh giá phẩm chất ◦ Đo đường kính ngang ngực D1.3, Dt theo chiều ĐT, NB lấy trung bình ◦ Dùng thước Blumleiss đo chiều cao vút Hvn 10 OTC Xi fi Xd Xt Xi 12 12 12 20 12 16 24 16 20 28 20 24 α Fi.Xi^a Pi 5.65685 28.28427 0.052771 88.181 617.271 0.21134 6 316.22 3794.73 0.30637 733.36 5133.554 0.253082 1374.6 12371.5 0.128142 2270.16 0.03989 22 2270.16 2.5 λ Xi ^ a fll Kiểm tra 2.1636 8.6651 0.1267 12.5615 0.0251 10.3763 1.0986 5.2538 1.404404 1.6356 0.0017 4788.21 24215.54 0.99161 X052 3.84 Xi ^ a Fi.Xi^a Pi 0.052409 40.6560 2.6548 H+ OTC Xi fi Xd Xt Xi 4 5.65685 22.6 12 88.1816 617.3 12 10 316.228 20 12 16 14 733.365 24 16 20 18 1374.62 28 20 24 22 2270.16 16 40 α 4788.21 2.5 λ X Kiểm tra 8.40454 0.0237234 12.2161 4400.2 10.1391 5.17084 12371.5 0.129271 0.04064 1.62594 2270.2 23792.75 0.991322 39.65289 4111.0 0.0017 05 0.21011 0.30540 0.253477 fll 2.09634 3.84 90 0.0503002 1.6897101 1.5096212 3.2733547 H+ Phụ biểu 05: Mơ hình hóa phân bố N-D1.3 hàm khoảng cách Trạng thái IIIA1 – OTC D1.3 ft 12 16 20 24 28 32 36 Tổng Xi 10 16 2 47 ft*Xi 28 gama anpha 20 48 20 10 12 124 0.0851 0.6387 Pi fll k.tra 0.0851 0.3306 15.54 3.7296 0.2111 9.923 0.0006 0.1348 6.338 14.731 0.0861 4.048 0.0468 0.055 2.585 0.0351 1.651 0.0224 1.055 0.9603 45.14 18.508 X052= 3.84 H- OTC D1.3 12 16 20 24 28 32 36 ft Xi ft*Xi Pi fll k.tra 0 0.1373 10 10 0.3487 17.78 3.4052 17 34 0.2077 10.59 3.8772 10 30 0.1238 6.316 2.1485 12 0.0738 3.764 0.1619 10 0.044 2.246 6 0.0263 1.339 7 0.0157 0.799 51 28 109 0.9772 49.84 9.5928 gama 0.1373 X05 = 5.99 anpha 0.5963 H- OTC D1.3 ft Xi ft*Xi Pi fll k.tra 0 0.12 12 6 0.2933 14.67 5.1212 16 10 20 0.1956 9.778 0.0051 20 15 45 0.1304 6.519 11.036 24 28 0.0869 4.346 0.6158 28 20 0.0579 2.897 32 6 0.0386 1.931 36 7 0.0258 1.288 50 28 132 0.9485 47.42 gama anpha X052= 0.12 0.6667 16.778 5.99 H- 91 Trạng thái IIB – OTC D1.3 ft Xi ft*Xi Pi fll k.tra 0 0.0833 12 6 0.3399 12.24 2.5523 16 14 0.2139 7.699 0.0635 20 13 39 0.1346 4.844 10.584 24 20 0.0847 3.048 28 10 0.0533 1.918 36 15 89 0.9096 32.75 gama 0.0833 anpha 0.6292 ft Xi 13.199 H- OTC D1.3 ft*Xi Pi fll k.tra 0 0.1463 12 9 0.3782 15.51 2.73 16 14 28 0.2106 8.636 3.331 20 21 0.1173 4.81 1.0145 24 4 16 0.0653 2.679 28 5 0.0364 1.492 41 15 79 gama 0.1463 anpha X052= 7.0755 5.99 0.557 H- OTC D1.3 ft Xi ft*Xi Pi fll k.tra 0 0.1282 12 8 0.3752 14.63 3.0066 16 14 28 0.2137 8.335 3.8499 20 18 0.1217 4.748 1.0055 24 20 0.0693 2.705 28 5 0.0395 1.541 39 15 79 0.9477 36.96 gama 0.1282 anpha X052= 0.5696 7.862 5.99 H- Phụ biểu 06: Mơ hình hóa phân bố N-Hvn hàm khoảng cách 92 Trạng thái IIIA1 OTC Xi fi Xi ft*Xi 11 14 13 10 15 17 47 Pi fll k.tra 0.10638 16.18 0.34329 4.13814 9.947 28 0.211652 1.6508 0.13003 6.114 30 2.46888 0.07996 28 3.7585 2.54651 2.310 15 0.049155 109 0.1063 0.6146 gama anhpha 10.8043 X052= 5.99 H- OTC Xi fi Xi ft*Xi Pi fll 0.11764 k.tra 11 11 0.354517 11 16 13 15 32 0.212077 0.12686 27 6.4702 3.870 20 0.075894 17 19 51 18.08 10.81 2.7727 0.1176 anhpha 0.98909 0.26965 10 0.045401 2.3154 1.385 12 0.027159 112 gama 2.48472 6.51617 0.5982 X052= 5.99 H- OTC Xi fi Xi ft*Xi Pi fll 0.18367 7 11 13 17 34 18 93 k.tra 0.31702 15.534 0.19390 9.5014 0.11860 5.8115 4.68836 5.91804 0.00611 15 17 24 0.072543 3.5546 0.04437 20 2.1742 49 103 gama 0.1836 anhpha 13.797 0.6117 X052= 3.18446 5.99 H- Trạng thái IIB OTC Xi fi Xi ft*Xi Pi fll 0 12 12 11 16 13 15 15 4 16 17 10 0.305289 1.2704 36 0.11111 0.40634 14.62 k.tra 0.22059 7.9412 0.11974 4.311 0.06500 2.3402 69 gama 0.1388 anhpha 0.35201 0.00044 1.1963 1.54874 0.5429 X052= 3.84 H+ OTC Xi fi Xi ft*Xi Pi 0 7 11 15 30 13 18 15 anhpha 0.1291 89 0.1 5.164 3.08056 4.61831 1.53329 0.07688 3.0752 1.831 10 0.045762 40 gama 0.1 0.36404 14.562 0.21679 8.671 k.tra 24 17 fll 05 X = 9.23216 3.84 H- 0.5955 94 OTC Xi fi Xi ft*Xi Pi 0 0.128205 11 14 0.352869 0.21004 28 13 15 4 17 fll k.tra 13.76 8.191 3.32245 4.11853 18 0.125025 4.876 0.07441 16 2.9024 1.727 15 0.044297 84 39 gama 0.12821 anhpha 1.28433 8.72531 X052= 3.84 H- 0.5952 Phụ biểu 07: Mơ hình hóa phân bố N-Hvn hàm Weibull Trạng thái IIIA1 OTC Xi fi Xd Xt Xi Xi^a Fi.Xi^a Pi 1 27 11 14 125 13 10 343 15 8 729 17 10 12 11 1331 47 30 42 36 2556 3993 0.035 0.99 14497 X052= fll Kiểm tra 3.84 α 0.0032 1.66038 46.8272 3.5597651 λ 0.026 1.2040888 0.16 7.607760 216 1.99057087 0.31 14.86556 1750 0.0503985 0.30 14.39479 3430 1.34175058 0.15 7.094617 5103 0.1770452 H+ OTC Xi fi Xd Xt Xi Xi ^ a Fi.Xi^a 1 11 15.58846 11 16 13 55.9017 129.641 15 10 243 Pi fll 0.05136 2.619438 10.53348 171.473 0.206539 15.42861 894.42719 0.302522 1166.776 0.254556 12.982344 6.760270 1215 0.13255 95 Kiểm tra 1.1252262 0.0211611 1.2215872 0.0154988 401.311 609.338 17 10 12 11 19 12 14 13 51 42 α 56 49 1455.782 λ 802.6232 1218.676 5474.976 0.04291 0.00849 2.188015 0.433026 0.99892 50.945189 2.3834732 2.5 X052= 0.009 H+ 3.84 OTC Xi fi Xd Fi.Xi^a Pi fll Kiểm tra 0.10539 5.164255 1 63 0.254092 12.45051 0.1480282 17 25 13.40121 0.9664289 6 49 9.7367 1.4340513 8 81 5.220383 0.9484755 10 12 11 121 425 0.273494 0.19870 294 0.10653 486 0.04362 484 49 α 30 42 36 286 1761 0.981855 48.11088 15 λ Xt Xi Xi ^ a 2.137822 3.4969839 X052= 0.028 H+ 3.84 Trạng thái IIB OTC Xi fi Xd Xt Xi Xi^a Fi.Xi^a Pi fll Kiểm tra 0.169979 6.119253 0.2047192 0.284113 6.47 77.7 10.22807 0.3069723 12 11 15.43 13 27.33 15 8 41.90 17 10 12 11 36 α 30 42 36 58.93 151.0 λ 123.4 0.2465006 8.8740208 0.159481 5.741327 136.7 0.083447 3.004124 167.6 1.317599 117.9 0.0366 628.21 0.98 0.086084 0.087237 35.28 0.69 1.7 0.0573 X052= 96 3.84 H+ OTC Xi fi Xd Xt Xi Xi^a Fi.Xi^a Pi fll Kiểm tra 0.042394 1.696 19.4190 135.9 0.2029457 8.118 0.1434298 11 15 6 15 8 1156.9 0.3234729 0.270631 1147.9 12.939 13 77.12923 191.321 9 377.0985 5.006 17 2262.6 0.1251602 0.031168 1296.6 40 α λ 10 12 11 648.2765 1314.24 30 42 36 0.995773 X052= 2.7 0.006 6003.9 0.328317 10.825 2.1508096 0.487991 1.247 39.8309 3.84 3.110547 H+ OTC Xi fi Xd Xt Xi Xi ^ a 11 14 13 6 15 8 77.12923 191.321 17 39 α λ 19.4190 377.0985 Fi.Xi^a 135.93316 1079.8092 1147.9314 Pi fll Kiểm tra 0.0426 1.6622 0.2039 7.9514 0.5923327 0.3243 12.6475 0.1446412 10.539 0.2702 1.9553617 1508.3939 0.1242 1944.8293 0.0307 10 12 11 648.2765 1314.24 30 42 36 5821.89715 0.9959 4.8438 0.152874 1.1954 38.840 2.8452 2.7 X052= 0.0067 3.84 H+ Phụ biểu 08: Chỉ tiêu thống kê tương quan, hàm LogHvn = a+b.logD1.3 Trạng thái IIIA1 OTC1 SUMMARY OUTPUT Regression Statistics Multiple R R Square Adjusted R Square 0.987571 0.975296 Standard Error Observations 0.020986 0.971179 97 ANOVA F 236.8783 Significanc eF 4.76E-06 t Stat 8.5708385 15.3908517 P-value 0.000139 4.76E-06 Lower 95% 0.284886 0.458314 SS 0.140203834 MS 0.14020383 F 287.8144 Significanc eF 2.68E-06 0.002922797 0.143126631 Standard Error 0.048925689 0.037238445 0.00048713 P-value 0.000838 2.68E-06 Lower 95% 0.18176 0.540635 Regression df SS 0.104323494 MS 0.10432349 Residual 0.002642458 0.00044041 Total Coefficient s 0.398717 0.544953 0.106965951 Standard Error 0.046520212 0.035407582 Intercept LogD (X) OTC SUMMARY OUTPUT Regression Statistics Multiple R R Square Adjusted R Square Standard Error Observations ANOVA 0.989737 0.979579 0.976175 0.022071 df Regression Residual Total Intercept log D Coefficient s 0.301477 0.631754 98 t Stat 6.16193001 16.9650938 OTC SUMMARY OUTPUT Regression Statistics Multiple R 0.962792268 R Square 0.926968951 Adjusted R Square 0.914797109 Standard Error 0.04074408 Observations ANOVA df Regression Residual Total SS Significance F MS F 0.126426445 0.00166008 76.15683688 Coefficients 0.126426445 0.00996048 0.136386925 Standard Error 0.00012521 t Stat P-value Lower 95% Intercept 0.307812591 0.090318723 3.408070687 0.014352657 0.086810639 log D 0.599910994 0.068743616 8.726788463 0.00012521 0.431701426 Trạng thái IIB OTC SUMMARY OUTPUT Regression Statistics Multiple R R Square Adjusted R Square Standard Error Observations ANOVA 0.948753 0.900132 0.875165 0.047622 Significanc eF df Regression Residual Total Intercept log D SS MS F Coefficient s 0.265921 0.081763908 0.009071554 0.090835462 Standard Error 0.129897381 0.08176391 0.00226789 36.05288 0.003872 t Stat 2.04716186 P-value 0.110044 Lower 95% -0.09473 0.632558 0.105348924 6.00440506 0.003872 0.340062 99 OTC SUMMARY OUTPUT Regression Statistics Multiple R R Square Adjusted R Square Standard Error Observations ANOVA 0.941063884 0.885601234 0.857001542 0.046881672 df MS 0.068058605 0.002197891 Coefficients 0.008791565 0.07685017 Standard Error t Stat P-value Lower 95% Intercept 0.370522544 0.127877061 2.897490298 0.044229133 0.015478904 log D 0.577112905 0.103710411 5.564657401 0.005107842 0.289166643 Regression Residual Total F 30.965412 Significance F 0.005107842 SS 0.068058605 OTC SUMMARY OUTPUT Regression Statistics Multiple R R Square Adjusted R Square 0.994388479 0.988808447 Standard Error Observations ANOVA 0.016204125 0.986010558 SS 0.09279679 0.001050295 MS 0.09279679 0.000262574 F 353.4124066 Significance F 4.71454E-05 Coefficients 0.093847084 Standard Error t Stat P-value Lower 95% Intercept 0.250345057 0.044199274 5.664008337 0.00479066 0.127628199 log D 0.673884924 0.035846342 18.79926612 4.71454E-05 0.574359523 df Regression Residual Total 100 Phụ biểu 09: Nắn phân bố N-D1.3 theo hàm Weibul trạng thái IIIA1 trạng thái IIB 101 Phụ biểu 10: Nắn phân bố N-Hvn theo hàm Weibul trạng thái IIIA1 trạng thái IIB 102 Phụ biểu 11: Biểu đồ tương quan chiều cao đường kính thân hai trạng thái IIIA1 IIB 103 Lời cảm ơn Để hồn thành chương trình đào tạo Đại học khóa 2009 – 2011, đồng ý khoa Lâm học – Trường Đại học Lâm nghiệp, thực đề tài: “Đặc trưng cấu trúc đa dạng loài thực vật phân khu phục hồi sinh thái – Vườn Quốc Gia Pù Mát – tỉnh Nghệ An” Sau thời gian tiến hành, đến đề tài hoàn thành Nhân dịp cho phép gửi lời cảm ơn thầy giáo, cô giáo khoa Lâm học – trường Đại học Lâm nghiệp, BQL Vườn Quốc Gia Pù Mát tạo điều kiện thuận lợi cho tơi suốt q trình thu thập xử lý số liệu Đặc biệt cho bày tỏ lòng biết ơn sâu sắc đến Th.S Hà Quang Anh trực tiếp hướng dẫn, giúp đỡ, bảo, động viên tơi suốt q trình thực đề tài Cuối cùng, gửi lời cảm ơn chân thành đến bạn bè, gia đình động viên, giúp đỡ tơi hồn thành khóa luận Tuy nhiên, khn khổ thời gian kinh nghiệm cịn hạn chế, đề tài nghiên cứu số đặc điểm cấu trúc, đa dạng loài Phân khu phục hồi sinh thái – VQG Pù Mát – tỉnh Nghệ An Do vậy, đề tài không tránh khỏi thiếu sót Rất mong nhận ý kiến dóng góp quý báu thầy giáo, cô giáo bạn bè đồng nghiệp để khóa luận hồn thiện Tôi xin chân thành cảm ơn Hà nội, ngày 14 tháng 05 năm 2011 Sinh viên Nguyễn Thị Hương 104 ... đa dạng sinh học VQG Pù Mát 2.2 Đối tượng phạm vi nghiên cứu - Đối tượng nghiên cứu: Trạng thái IIIA1, IIB thuộc phân khu phục hồi sinh thái – VQG Pù Mát - Phạm vi nghiên cứu: Phân khu phục hồi. .. tồn hệ sinh thái rừng, giá trị khoa học đa dạng sinh học, địa chất, cảnh quan, loài động – thực vật hoang dã tồn sinh sống hệ sinh thái rừng Vì vậy, tính đa dạng sinh học quan trọng (đặc điểm khác... 4.4.3 Đa dạng dạng sống Việc nghiên cứu phổ dạng sống để từ thành lập phổ dạng sống hệ thực vật vùng cho ta hình ảnh chất sinh thái hệ thực vật vùng nghiên cứu Hệ thực vật Pù Mát giống hệ thực vật