Luận án tiến sĩ nghiên cứu bổ sung cơ sở khoa học về kỹ thuật trồng rừng phòng hộ trên các dạng lập địa chính vùng cát ven biển các tỉnh hà tĩnh, quảng bình và quảng trị

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP & PTNT VIỆN KHOA HỌC LÂM NGHIỆP VIỆT NAM LÊ ĐỨC THẮNG NGHIÊN CỨU BỔ SUNG CƠ SỞ KHOA HỌC VỀ KỸ THUẬT TRỒNG RỪNG PHÒNG HỘ TRÊN CÁC DẠNG LẬP ĐỊA CHÍNH VÙNG CÁT VEN BIỂN CÁC TỈNH HÀ TĨNH, QUẢNG BÌNH VÀ QUẢNG TRỊ LUẬN ÁN TIẾN SĨ LÂM NGHIỆP Hà Nội - 2022 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP & PTNT VIỆN KHOA HỌC LÂM NGHIỆP VIỆT NAM LÊ ĐỨC THẮNG NGHIÊN CỨU BỔ SUNG CƠ SỞ KHOA HỌC VỀ KỸ THUẬT TRỒNG RỪNG PHÒNG HỘ TRÊN CÁC DẠNG LẬP ĐỊA CHÍNH VÙNG CÁT VEN BIỂN CÁC TỈNH HÀ TĨNH, QUẢNG BÌNH VÀ QUẢNG TRỊ Ngành đào tạo: Lâm sinh Mã ngành: 62 02 05 LUẬN ÁN TIẾN SĨ LÂM NGHIỆP Người hướng dẫn khoa học: PGS TS NGƠ ĐÌNH QUẾ GS TS NGUYỄN XUẤT QUÁT Hà Nội, 2022 i LỜI CAM ĐOAN Luận án hoàn thành Viện Khoa học Lâm nghiệp Việt Nam theo Chương trình đào tạo nghiên cứu sinh khóa 27 Tơi xin cam đoan cơng trình nghiên cứu khoa học thân tôi, luận án thực thời gian từ năm 2015 đến năm 2021 hướng dẫn khoa học PGS TS Ngơ Đình Quế GS TS Nguyễn Xuân Quát Các liệu, kết phân tích, đánh giá hình, biểu đồ trình bày luận án hoàn toàn trung thực chưa cơng bố cơng trình Nội dung luận án có kế thừa sử dụng phần kết Đề tài độc lập cấp Nhà nước: “Nghiên cứu phát triển bền vững hệ thống rừng phòng hộ dải ven biển Bắc Trung bộ”, mã số ĐTĐL 2012.T/33 Viện Nghiên cứu Phát triển Vùng - Bộ Khoa học Cơng nghệ chủ trì, đó, tác giả thư ký khoa học Đề tài NCS trực tiếp thực công việc thiết kế, bố trí thí nghiệm; xây dựng thí nghiệm trồng rừng; điều tra, thu thập số liệu ngoại nghiệp thí nghiệm xây dựng tỉnh Hà Tĩnh, Quảng Bình Quảng Trị việc tổng hợp, xử lý, phân tích, đánh giá tham gia viết sản phẩm chuyên đề, báo báo cáo tổng kết đề tài Các thông tin, liệu kết liên quan, tác giả cho phép sử dụng đơn vị chủ trì, chủ nhiệm đề tài thành viên tham gia thực đề tài để công bố luận án Hà Nội, tháng 04 năm 2022 Tác giả luận án NCS Lê Đức Thắng ii LỜI CẢM ƠN Trong q trình thực luận án, tơi nhận nhiều giúp đỡ, động viên Thầy cô, quan, bạn bè đồng nghiệp Tôi xin chân thành cảm ơn sâu sắc giúp đỡ quý báu đó: Ban lãnh đạo Viện Khoa học Lâm nghiệp Việt Nam, Ban Khoa học, Đào tạo Hợp tác Quốc tế, Viện Nghiên cứu Lâm sinh, Viện Nghiên cứu Phát triển Vùng - nơi NCS công tác,…đã tạo điều kiện thuận lợi cho NCS suốt q trình học tập nghiên cứu để hồn thiện luận án Trong trình thực luận án, NCS nhận hướng dẫn, giúp đỡ, bảo tận tình, chu đáo quý Thầy giáo hướng dẫn khoa học PGS TS Ngơ Đình Quế GS TS Nguyễn Xuân Quát Qua đây, NCS xin bày tỏ lòng biết ơn sâu sắc kính trọng đến quý Thầy dành nhiều thời gian công sức để hướng dẫn tận tình giúp đỡ NCS hồn thành luận án Tôi xin trân trọng cảm ơn nhà khoa học: GS TS Võ Đại Hải, PGS TS Phạm Xuân Hoàn, PGS TS Đặng Thái Dương, PGS TS Nguyễn Huy Sơn, … số nhà khoa học khác có ý kiến góp ý quý báu cho luận án Xin gửi lời cảm ơn chân thành sâu sắc tới quý Thầy, gia đình, bạn bè, đồng nghiệp động viên, hỗ trợ vật chất tinh thần để tơi hồn thành tốt luận án Một lần xin gửi lời cảm ơn chân thành sâu sắc giúp đỡ quý báu đỡ quý báu Hà Nội, tháng 04 năm 2022 Tác giả luận án NCS Lê Đức Thắng iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC TỪ, KÝ HIỆU VIẾT TẮT vi DANH MỤC CÁC BẢNG, BIỂU viii DANH MỤC CÁC HÌNH, BIỂU ĐỒ x MỞ ĐẦU 1 Tính cấp thiết vấn đề nghiên cứu Mục tiêu nghiên cứu 2.1 Mục tiêu lý luận 2.2 Mục tiêu thực tiễn Đối tượng, phạm vi giới hạn nghiên cứu 3.1 Đối tượng nghiên cứu 3.2 Phạm vi nghiên cứu Ý nghĩa khoa học thực tiễn luận án 4.1 Ý nghĩa khoa học 4.2 Ý nghĩa thực tiễn Những đóng góp luận án Cấu trúc bố cục luận án CHƯƠNG TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Trên giới 1.1.1 Cơ sở khoa học tiêu chí phân chia lập địa vùng cát ven biển 1.1.2 Trồng rừng phòng hộ vùng cát ảnh hưởng biện pháp kỹ thuật đến sinh trưởng phát triển rừng phịng hộ chắn gió, chắn cát bay 1.1.3 Cơ sở khoa học việc đánh giá hiệu phòng hộ đai rừng chắn gió, chắn cát bay 11 1.2 Ở nước 14 1.2.1 Cơ sở khoa học tiêu chí phân chia lập địa vùng cát ven biển 14 iv 1.2.2 Trồng rừng phòng hộ vùng cát ảnh hưởng biện pháp kỹ thuật đến sinh trưởng phát triển rừng phịng hộ chắn gió, chắn cát bay 16 1.2.3 Cơ sở khoa học việc đánh giá hiệu phòng hộ đai rừng chắn gió, chắn cát bay 24 1.3 Thảo luận chung 26 CHƯƠNG NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 28 2.1 Nội dung nghiên cứu 28 2.2 Phương pháp nghiên cứu 29 2.2.1 Quan điểm, cách tiếp cận phương pháp luận nghiên cứu 29 2.2.2 Phương pháp nghiên cứu kỹ thuật sử dụng 30 2.3 Đặc điểm điều kiện tự nhiên khu vực nghiên cứu 44 CHƯƠNG KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 50 3.1 Cơ sở khoa học xây dựng tiêu chí phân chia nhóm dạng lập địa vùng cát ven biển 50 3.1.1 Cơ sở xây dựng tiêu chí phân chia nhóm dạng lập địa 50 3.1.2 Phân chia nhóm dạng lập địa trồng rừng phòng hộ vùng cát ven biển 69 3.2 Đánh giá biện pháp kỹ thuật trồng rừng phòng hộ sinh trưởng lâm phần rừng trồng số nhóm dạng lập địa vùng cát ven biển 74 3.2.1 Đánh giá biện pháp kỹ thuật trồng rừng vùng cát ven biển 74 3.2.2 Sinh trưởng loài trồng rừng phòng hộ số dạng lập địa vùng cát ven biển 82 3.2.3 Một số thành công, tồn học kinh nghiệm công tác trồng rừng, bảo vệ phát triển rừng vùng cát ven biển 86 3.3 Ảnh hưởng biện pháp kỹ thuật trồng rừng phịng hộ số nhóm dạng lập địa vùng cát ven biển 89 3.3.1 Ảnh hưởng phân bón lót kết hợp chất giữ ẩm đến sinh trưởng Phi lao nhóm dạng lập địa II 89 3.3.2 Ảnh hưởng phân bón lót kết hợp chất giữ ẩm đến sinh trưởng Keo liềm nhóm dạng lập địa III1 95 v 3.3.3 Sinh trưởng Keo liềm Keo tràm nhóm dạng lập địa III1 99 3.3.4 Ảnh hưởng trồng phù trợ vị trí trồng (trước sau đai rừng Phi lao có) đến sinh trưởng Phi lao trồng nhóm dạng lập địa III2 100 3.3.5 Ảnh hưởng phân bón lót kết hợp chất giữ ẩm đến sinh trưởng Keo liềm nhóm dạng lập địa III3 102 3.3.6 Ảnh hưởng dạng lập địa trồng rừng đến sinh trưởng Phi lao 109 3.3.7 Ảnh hưởng dạng lập địa trồng rừng đến sinh trưởng Keo liềm 112 3.4 Hiệu phòng hộ đai rừng vùng cát ven biển khu vực nghiên cứu 115 3.4.1 Hiệu chắn gió đai rừng 115 3.4.2 Hiệu chắn cát đai rừng 118 3.4.3 Hiệu cải thiện đất đai rừng 121 3.4.4 Hiệu tích lũy carbon hấp thụ CO2 đai rừng 123 3.5 Giải pháp kỹ thuật trồng rừng phòng hộ phù hợp, hiệu bền vững nhóm dạng lập địa vùng cát ven biển 132 KẾT LUẬN, TỒN TẠI, VÀ KHUYẾN NGHỊ 142 Kết luận 142 Tồn 143 Khuyến nghị 143 DANH MỤC CÁC CÔNG TRÌNH KHOA HỌC ĐÃ CƠNG BỐ LIÊN QUAN ĐẾN LUẬN ÁN 144 TÀI LIỆU THAM KHẢO 145 PHỤ LỤC 158 vi DANH MỤC CÁC TỪ, KÝ HIỆU VIẾT TẮT Từ viết tắt Từ viết đầy đủ BQL Ban quản lý BTB Bắc Trung BĐKH Biến đổi khí hậu Cc Cồn cát trắng vàng C Đất cát biển CTTN Cơng thức thí nghiệm D0 Đường kính gốc D1.3 Đường kính thân vị trí 1,3m ∆D1.3 Tăng trưởng bình qn chung đường kính thân DT Đường kính tán ĐC Đối chứng ĐDSH Đa dạng sinh học FAO Tổ chức Nông lương Liên Hiệp Quốc HVN Chiều cao ∆HVN Tăng trưởng bình quân chung chiều cao KH&CN Khoa học Công nghệ KTC Khoảng tin cậy M Trữ lượng ∆M Tăng trưởng bình quân chung trữ lượng NA Số liệu trống không NĐ Nghị định NLKH Nông lâm kết hợp NN&PTNT Nông nghiệp Phát triển Nơng thơn OTC Ơ tiêu chuẩn QĐ-TTg Quyết định Thủ tướng Chính phủ RĐD Rừng đặc dụng vii Từ viết tắt Từ viết đầy đủ BQL Ban quản lý RPH Rừng phòng hộ RSX Rừng sản xuất SD Sai tiêu chuẩn TT Thông tư TCVN Tiêu chuẩn Việt Nam UBND Ủy ban nhân dân UNDP Chương trình Phát triển Liên hợp quốc viii DANH MỤC CÁC BẢNG, BIỂU Bảng 3.1 Các dạng địa hình, địa mạo vùng cát ven biển 53 Bảng 3.2 Diện tích phân bố loại đất cát ven biển khu vực nghiên cứu 55 Bảng 3.3 Thành phần cấp hạt mẫu đất vùng cát ven biển khu vực nghiên cứu 56 Bảng 3.4 Một số tính chất hóa tính đất vùng cát ven biển khu vực nghiên cứu 62 Bảng 3.5 Diện tích phân bố loại đất cát ven biển theo độ cao 64 Bảng 3.6 Khả thoát nước, giữ nước đất cát vùng ven biển 65 Bảng 3.7 Trạng thái thực bì vùng đất cát ven biển 68 Bảng 3.8 Tổng hợp yếu tố cấu thành nhóm dạng lập địa vùng cát ven biển khu vực nghiên cứu 70 Bảng 3.9 Tổng hợp nhóm dạng lập địa vùng cát ven biển khu vực nghiên cứu 72 Bảng 3.10 Diện tích rừng đất, cát ven biển khu vực nghiên cứu 75 Bảng 3.11 Kết thực nhiệm vụ bảo vệ phát triển rừng chắn cát, chắn gió ven biển giai đoạn 2015 - 2020 77 Bảng 3.12 Tổng hợp biện pháp kỹ thuật áp dụng trồng rừng vùng cát ven biển khu vực nghiên cứu 78 Bảng 3.13 Một số tiêu sinh trưởng lồi trồng rừng số dạng lập địa vùng cát ven biển khu vực nghiên cứu 83 Bảng 3.14 Lượng tăng trưởng bình quân chung khoảng tin cậy 95% loài trồng rừng số dạng lập địa vùng cát ven biển khu vực nghiên cứu 85 Bảng 3.15 Sinh trưởng đường kính gốc, chiều cao đường kính tán Phi lao trồng nhóm dạng lập địa II Lệ Thủy 89 Bảng 3.16 Số cành dài 50cm, tỷ lệ chết tỷ lệ sống Phi lao trồng nhóm dạng lập địa II huyện Lệ Thủy 91 Bảng 3.17 Sinh trưởng đường kính gốc, chiều cao đường kính tán Phi lao trồng nhóm dạng lập địa II Triệu Phong 92 Bảng 3.18 Số cành dài 50cm, tỷ lệ chết tỷ lệ sống Phi lao trồng nhóm dạng lập địa II Triệu Phong 93 Bảng 3.19 Sinh trưởng đường kính gốc, chiều cao đường kính tán Keo liềm trồng nhóm dạng lập địa III1 Cẩm Xuyên 96 Bảng 3.20 Tỷ lệ sống, số thân chính, số cành 50cm lượng vật rơi rụng Keo liềm trồng nhóm dạng lập địa III1 Cẩm Xuyên 98 12 $means CT1 CT2 CT3 DC Hvn_m 1.0644444 1.2842308 1.3987097 0.8855556 std r LCL UCL 0.2867103 0.8672334 1.261655 0.2936007 26 1.1682020 1.400260 0.3229008 31 1.2924493 1.504970 0.1969207 0.6883446 1.082767 Min 0.60 0.43 0.57 0.70 Max 1.50 1.76 1.87 1.30 Q25 0.9000 1.1325 1.1750 0.7700 Q50 1.10 1.27 1.44 0.80 $groups CT3 CT2 CT1 DC Hvn_m groups 1.3987097 a 1.2842308 ab 1.0644444 bc 0.8855556 c attr(,"class") [1] "group" > dt=aov(Dtan_tb~CTTN,data=TP14) > summary(dt) Df Sum Sq Mean Sq F value Pr(>F) CTTN 0.942 0.31398 11.49 3.15e-06 *** Residuals 71 1.940 0.02732 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 22 observations deleted due to missingness > dt1=LSD.test(dt) > dt1 Linear mixed model fit by REML ['lmerMod'] Formula: DelDt ~ (1 | Dong) Data: DelDt1 REML criterion at convergence: 458.2522 Random effects: Groups Name Std.Dev Dong (Intercept) 0.0000 Residual 0.3259 Number of obs: 759, groups: Dong, 20 Fixed Effects: (Intercept) 0.4456 convergence code 0; optimizer warnings; lme4 warnings > dt1=LSD.test(dt,"CTTN") > dt1 $statistics MSerror Df Mean CV 0.02732181 71 0.658 25.12053 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means CT1 CT2 CT3 DC Dtan_tb 0.5788889 0.6792308 0.7416129 0.3877778 std r LCL UCL 0.2281143 0.4690272 0.6887506 0.1214882 26 0.6145938 0.7438677 0.1831411 31 0.6824177 0.8008081 0.1361780 0.2779161 0.4976395 $groups Dtan_tb groups CT3 0.7416129 a CT2 0.6792308 ab CT1 0.5788889 b DC 0.3877778 c attr(,"class") [1] "group" > c=aov(Bough_50cm~CTTN,data=TP14) > s=aov(Song_chet~CTTN,data=TP14) Min 0.23 0.45 0.35 0.25 Max 0.90 0.88 1.00 0.63 Q25 0.4000 0.5625 0.6000 0.3000 Q50 0.63 0.69 0.80 0.35 Q75 0.70 0.75 0.88 0.48 Q75 1.16 1.48 1.65 0.94 13 > summary(s) Df Sum Sq Mean Sq F value Pr(>F) CTTN 3.217 1.0724 7.23 0.000205 *** Residuals 93 13.793 0.1483 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ > s1=LSD.test(s,"CTTN") > s1 $statistics MSerror Df Mean CV 0.1483143 93 0.7731959 49.80832 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means Song_chet CT1 1.0000000 CT2 1.0000000 CT3 0.6078431 DC 0.8181818 std r LCL UCL Min Max Q25 Q50 Q75 0.0000000 0.7450787 1.2549213 1 1 0.0000000 26 0.8500175 1.1499825 1 1 0.4930895 51 0.5007548 0.7149315 1 0.4045199 11 0.5875968 1.0487668 1 1 $groups Song_chet groups CT1 1.0000000 a CT2 1.0000000 a DC 0.8181818 ab CT3 0.6078431 b attr(,"class") 2.2 Giai đoạn 24 tháng tuổi > summary(d) Df Sum Sq Mean Sq F value Pr(>F) CTTN 18.79 6.263 36.88 d1=LSD.test(d,"CTTN") > d1 $statistics MSerror Df Mean CV 0.1698261 113 1.802222 22.86619 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means CT1 CT2 CT3 DC Dgoc 2.031034 1.846190 2.011915 0.931500 std 0.4830953 0.4290044 0.4343543 0.1248694 $groups Dgoc groups CT1 2.031034 a CT3 2.011915 a CT2 1.846190 a r 29 21 47 20 LCL 1.8794247 1.6680279 1.8928244 0.7489377 UCL 2.182644 2.024353 2.131005 1.114062 Min 1.02 1.08 1.27 0.64 Max 3.09 2.67 3.50 1.15 Q25 1.7500 1.5000 1.7350 0.8825 Q50 2.01 1.85 1.97 0.95 Q75 2.26 2.01 2.29 1.02 14 DC 0.931500 b attr(,"class") [1] "group" > h=aov(Hvn~CTTN,data=PLTP) > h=aov(Hvn_m~CTTN,data=PLTP) > summary(h) Df Sum Sq Mean Sq F value Pr(>F) CTTN 7.295 2.4316 18.03 1.24e-09 *** Residuals 113 15.237 0.1348 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ > h1=LSD.test(h,"CTTN") > h1 $statistics MSerror Df Mean CV 0.1348376 113 1.25547 29.2482 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means CT1 CT2 CT3 DC Hvn_m 1.397241 1.345238 1.361277 0.707000 std 0.4743108 0.3272861 0.3605408 0.2071892 r 29 21 47 20 LCL 1.2621491 1.1864859 1.2551607 0.5443274 UCL 1.5323336 1.5039903 1.4673925 0.8696726 Min 0.40 0.85 0.52 0.40 Max 2.40 1.95 2.10 1.25 Q25 1.1000 1.1000 1.0750 0.5575 Q50 1.35 1.32 1.30 0.70 Q75 1.6500 1.5000 1.5800 0.8125 $groups Hvn_m groups CT1 1.397241 a CT3 1.361277 a CT2 1.345238 a DC 0.707000 b attr(,"class") [1] "group" > dt=aov(Dtan_tb~CTTN,data=PLTP) > summary(dt) Df Sum Sq Mean Sq F value Pr(>F) CTTN 7.166 2.3885 35.96 dt1=LSD.test(dt,"CTTN") > dt1 $statistics MSerror Df Mean CV 0.06642651 113 1.012051 25.46644 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means Dtan_tb CT1 1.153793 CT2 1.016667 CT3 1.149574 DC 0.478500 std 0.3228712 0.2238377 0.2697016 0.1122157 $groups Dtan_tb groups CT1 1.153793 a CT3 1.149574 a CT2 1.016667 a r 29 21 47 20 LCL 1.0589740 0.9052410 1.0750934 0.3643227 UCL 1.2486122 1.1280923 1.2240555 0.5926773 Min 0.65 0.50 0.65 0.30 Max 1.98 1.45 1.73 0.65 Q25 0.950 0.940 0.975 0.380 Q50 1.13 1.03 1.11 0.50 Q75 1.300 1.180 1.365 0.580 15 DC 0.478500 b attr(,"class") [1] "group" > c=aov(Bough_50cm~CTTN,data=PLTP) > summary(c) Df Sum Sq Mean Sq F value Pr(>F) CTTN 513.4 171.14 14.04 7.63e-08 *** Residuals 113 1377.1 12.19 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ > c1=LSD.test(c,"CTTN") > c1 $statistics MSerror Df Mean CV 12.18672 113 10.70085 32.62309 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means Bough_50cm CT1 11.20690 CT2 10.47619 CT3 12.34043 DC 6.35000 std 4.34560 2.78602 3.55868 2.41214 r LCL UCL Min Max Q25 Q50 Q75 29 9.922591 12.491203 20 10 14 21 8.966952 11.985429 16 10 12 47 11.331594 13.349257 23 12 14 20 4.803491 7.896509 12 $groups Bough_50cm groups CT3 12.34043 a CT1 11.20690 ab CT2 10.47619 b DC 6.35000 c attr(,"class") [1] "group" > s=aov(Song_chet~CTTN,data=PLTP) > summary(s) Df Sum Sq Mean Sq F value Pr(>F) CTTN 2.17 0.7229 4.282 0.00573 ** Residuals 244 41.19 0.1688 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ > s1=LSD.test(s,"CTTN") > s1 $statistics MSerror Df Mean CV 0.168796 244 0.7741935 53.06786 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means Song_chet CT1 0.7037037 CT2 0.9591837 CT3 0.7289720 DC 0.7631579 std r LCL 0.4609109 54 0.5935772 0.1999149 49 0.8435749 0.4465823 107 0.6507377 0.4308515 38 0.6318784 $groups Song_chet groups CT2 0.9591837 a DC 0.7631579 b UCL Min Max Q25 Q50 Q75 0.8138302 1 1.0747924 1 1 0.8072062 1 0.8944374 1 1 16 CT3 0.7289720 CT1 0.7037037 b b attr(,"class") [1] "group" 2.3 KEO LÁ LIỀM LẬP ĐỊA B VÀ C1 > library(agricolae) > d=aov(dg~Lapdia,data=KL) > summary(d) Df Sum Sq Mean Sq F value Pr(>F) Lapdia 29.49 29.486 24.47 1.86e-06 *** Residuals 164 197.63 1.205 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ observations deleted due to missingness > d=LSD.test(d,"Lapdia") > d $statistics MSerror Df Mean CV 1.205073 164 4.213855 26.05116 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none Lapdia 0.05 $means dg std r LCL UCL Min Max Q25 Q50 Q75 B 3.845000 1.033559 94 3.621433 4.068567 1.66 6.62 3.12 3.675 4.5725 C1 4.695417 1.176562 72 4.439967 4.950866 2.16 7.83 3.82 4.620 5.4100 $groups dg groups C1 4.695417 a B 3.845000 b attr(,"class") [1] "group" > h=aov(hvn~Lapdia,data=KL) > summary(h) Df Sum Sq Mean Sq F value Pr(>F) Lapdia 8.98 8.979 32.34 5.81e-08 *** Residuals 164 45.54 0.278 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ observations deleted due to missingness > h1=LSD.test(h,"Lapdia") > h1 $statistics MSerror Df Mean CV 0.2776672 164 1.651205 31.91254 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none Lapdia 0.05 $means hvn std r LCL UCL Min Max Q25 Q50 Q75 B 1.447660 0.4841341 94 1.340344 1.554975 0.75 3.0 1.0500 1.35 1.800 C1 1.916944 0.5782383 72 1.794325 2.039564 0.65 3.1 1.4375 1.85 2.325 $groups hvn groups 17 C1 1.916944 B 1.447660 a b attr(,"class") [1] "group" > t=aov(dt~Lapdia,data=KL) > summary(t) Df Sum Sq Mean Sq F value Pr(>F) Lapdia 1.17 1.1657 5.285 0.0228 * Residuals 164 36.18 0.2206 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ observations deleted due to missingness > t1=LSD.test(t,"Lapdia") > t1 $statistics MSerror Df Mean CV 0.220586 164 2.078554 22.5958 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none Lapdia 0.05 $means dt std r LCL UCL Min Max Q25 Q50 Q75 B 2.005213 0.5082581 94 1.909562 2.100864 0.90 3.65 1.7000 1.925 2.31 C1 2.174306 0.4137048 72 2.065014 2.283597 1.25 3.50 1.9225 2.165 2.40 $groups dt groups C1 2.174306 a B 2.005213 b attr(,"class") [1] "group" > fix(KL) > th=aov(Thanchinh~Lapdia,data=KL) > summary(th) Df Sum Sq Mean Sq F value Pr(>F) Lapdia 3.11 3.1125 4.157 0.0431 * Residuals 164 122.79 0.7487 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ observations deleted due to missingness > th1=LSD.test(th,"Lapdia") > th1 $statistics MSerror Df Mean CV 0.7487261 164 1.975904 43.7921 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none Lapdia 0.05 $means Thanchinh std r LCL UCL Min Max Q25 Q50 Q75 B 2.095745 0.8559614 94 1.919522 2.271968 2 C1 1.819444 0.8773582 72 1.618091 2.020798 2 $groups Thanchinh groups B 2.095745 a C1 1.819444 b attr(,"class") [1] "group" 18 > fix(KL) > c=aov(c50~Lapdia,data=KL) > summary(c) Df Sum Sq Mean Sq F value Pr(>F) Lapdia 133 133.4 4.616 0.0331 * Residuals 164 4739 28.9 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ observations deleted due to missingness > c1=LSD.test(c,"Lapdia") > c1 $statistics MSerror Df Mean CV 28.89721 164 13.29518 40.43279 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none Lapdia 0.05 $means c50 std r LCL UCL Min Max Q25 Q50 Q75 B 12.51064 5.134209 94 11.41585 13.60542 27 12 16 C1 14.31944 5.676307 72 13.06853 15.57036 29 10 14 18 $groups c50 groups C1 14.31944 a B 12.51064 b attr(,"class") [1] "group" > pt=aov(phth~Lapdia,data=KL) > summary(pt) Df Sum Sq Mean Sq F value Pr(>F) Lapdia 3.56 3.560 23.18 3.33e-06 *** Residuals 164 25.19 0.154 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ observations deleted due to missingness > pt1=LSD.test(pt,"Lapdia") > pt1 $statistics MSerror Df Mean CV 0.1536138 164 0.7771084 50.43518 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none Lapdia 0.05 $means phth std r LCL UCL Min Max Q25 Q50 Q75 B 0.6489362 0.4798621 94 0.5691154 0.728757 1 C1 0.9444444 0.2306689 72 0.8532405 1.035648 1 1 $groups phth groups C1 0.9444444 a B 0.6489362 b attr(,"class") [1] "group" > s=aov(sc~Lapdia,data=KL) > summary(s) Df Sum Sq Mean Sq F value Pr(>F) Lapdia 0.005 0.00471 0.106 0.745 19 Residuals 172 7.627 0.04435 > s1=LSD.test(s,"Lapdia") > s1 $statistics MSerror Df Mean CV 0.04434578 172 0.954023 22.0733 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none Lapdia 0.05 $means sc std r LCL UCL Min Max Q25 Q50 Q75 B 0.9494949 0.2200991 99 0.9077193 0.9912706 1 1 C1 0.9600000 0.1972788 75 0.9120034 1.0079966 1 1 $groups sc groups C1 0.9600000 a B 0.9494949 a attr(,"class") [1] "group" 2.4 KEO LÁ LIỀM LẬP ĐỊA C3 TẠI LỆ THỦY > d=aov(stump_diameter~CTTN,data=LT13) > summary(d) Df Sum Sq Mean Sq F value Pr(>F) CTTN 0.45 0.1484 0.339 0.797 Residuals 223 97.76 0.4384 16 observations deleted due to missingness > library(agricolae) > d1=LSD.test(d,"CTTN") > d1 $statistics MSerror Df Mean CV 0.4383911 223 2.642863 25.05279 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means stump_diameter CT 2.660294 CT 2.703971 CT 2.606264 DC 2.601176 std 0.6693347 0.5833105 0.6681844 0.7784552 r 34 68 91 34 LCL 2.436524 2.545741 2.469484 2.377406 UCL 2.884065 2.862200 2.743043 2.824947 $groups stump_diameter groups CT 2.703971 a CT 2.660294 a CT 2.606264 a DC 2.601176 a attr(,"class") [1] "group" > h=aov(tree_height~CTTN,data=LT13) > summary(h) Df Sum Sq Mean Sq F value Pr(>F) CTTN 0.15 0.05006 0.965 0.41 Min 1.72 1.59 1.37 1.46 Max 4.81 4.30 5.43 4.62 Q25 2.2300 2.2525 2.1000 1.9950 Q50 2.48 2.58 2.51 2.53 Q75 2.9675 3.1425 2.9300 2.9125 20 Residuals 223 11.56 0.05186 16 observations deleted due to missingness > h1=LSD.test(h,"CTTN") > h1 $statistics MSerror Df Mean CV 0.0518596 223 0.9486344 24.00577 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means tree_height CT 0.9008824 CT 0.9652941 CT 0.9401099 DC 0.9858824 std 0.2211042 0.2092023 0.2336783 0.2525400 r 34 68 91 34 LCL 0.8239185 0.9108725 0.8930658 0.9089185 UCL 0.9778462 1.0197158 0.9871540 1.0628462 Min 0.60 0.53 0.43 0.60 Max 1.55 1.33 1.92 1.47 Q25 0.725 0.800 0.795 0.805 Q50 0.875 0.940 0.940 0.940 Q75 1.0375 1.1400 1.1000 1.2300 $groups tree_height groups DC 0.9858824 a CT 0.9652941 a CT 0.9401099 a CT 0.9008824 a attr(,"class") [1] "group" > dt=aov(canopy_diameter~CTTN,data=LT13) > summary(dt) Df Sum Sq Mean Sq F value Pr(>F) CTTN 0.351 0.11708 1.73 0.162 Residuals 223 15.095 0.06769 16 observations deleted due to missingness > dt1=LSD.test(dt,"CTTN") > dt1 $statistics MSerror Df Mean CV 0.067692 223 1.100441 23.64297 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means canopy_diameter CT 1.071471 CT 1.158676 CT 1.084615 DC 1.055294 std 0.2543814 0.2402045 0.2810666 0.2452115 r 34 68 91 34 LCL 0.9835399 1.0965001 1.0308677 0.9673634 UCL 1.159401 1.220853 1.138363 1.143225 $groups canopy_diameter groups CT 1.158676 a CT 1.084615 a CT 1.071471 a DC 1.055294 a attr(,"class") [1] "group" > th=aov(main_trunk~CTTN,data=LT13) > summary(th) Df Sum Sq Mean Sq F value Pr(>F) CTTN 1.34 0.4463 0.403 0.751 Min 0.67 0.75 0.48 0.45 Max 1.68 1.83 2.03 1.65 Q25 0.915 0.950 0.845 0.935 Q50 1.035 1.145 1.090 1.050 Q75 1.3000 1.3000 1.3050 1.1875 21 Residuals 223 246.83 1.1069 16 observations deleted due to missingness > th1=LSD.test(th,"CTTN") > th1 $statistics MSerror Df Mean CV 1.106854 223 2.550661 41.247 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means main_trunk CT 2.411765 CT 2.588235 CT 2.527473 DC 2.676471 std 1.0478737 1.0683408 0.9468152 1.2725681 r 34 68 91 34 LCL 2.056201 2.336814 2.310134 2.320907 UCL Min Max Q25 Q50 Q75 2.767328 2.0 2.839657 3.0 2.744811 3.0 3.032034 2.5 $groups main_trunk groups DC 2.676471 a CT 2.588235 a CT 2.527473 a CT 2.411765 a attr(,"class") [1] "group" > c=aov(bough_50_cm~CTTN,data=LT13) > summary(c) Df Sum Sq Mean Sq F value Pr(>F) CTTN 2.78 0.9261 1.251 0.292 Residuals 223 165.09 0.7403 16 observations deleted due to missingness > c1=LSD.test(c,"CTTN") > c1 $statistics MSerror Df Mean CV 0.7403125 223 1.903084 45.21158 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means bough_50_cm CT 1.647059 CT 1.985294 CT 1.934066 DC 1.911765 std 0.7739060 0.9540588 0.8406637 0.7926804 r 34 68 91 34 LCL 1.356269 1.779674 1.756321 1.620975 UCL Min Max Q25 Q50 Q75 1.937849 2 2.190914 2 2.111811 2.202555 2 $groups bough_50_cm groups CT 1.985294 a CT 1.934066 a DC 1.911765 a CT 1.647059 a attr(,"class") [1] "group" > pt=aov(phan_than~CTTN,data=LT13) > summary(pt) Df Sum Sq Mean Sq F value Pr(>F) CTTN 0.47 0.1560 0.814 0.487 Residuals 223 42.71 0.1915 22 16 observations deleted due to missingness > pt1=LSD.test(pt,"CTTN") > pt1 $statistics MSerror Df Mean CV 0.1915366 223 0.2555066 171.2868 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means phan_than CT 0.2058824 CT 0.2647059 CT 0.2307692 DC 0.3529412 std 0.4104256 0.4444566 0.4236593 0.4850713 r 34 68 91 34 LCL 0.05797222 0.16011762 0.14035918 0.20503104 UCL Min Max Q25 Q50 Q75 0.3537925 0 0.3692941 0 0.3211793 0 0.5008513 0 $groups phan_than groups DC 0.3529412 a CT 0.2647059 a CT 0.2307692 a CT 0.2058824 a attr(,"class") [1] "group" > s=aov(song_chet~CTTN,data=LT13) > summary(s) Df Sum Sq Mean Sq F value Pr(>F) CTTN 0.423 0.14103 2.321 0.0759 Residuals 239 14.523 0.06077 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ > s1=LSD.test(s,"CTTN") > s1 $statistics MSerror Df Mean CV 0.06076742 239 0.9341564 26.38857 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means song_chet CT 1.0000000 CT 0.9189189 CT 0.9009901 DC 1.0000000 std r LCL 0.0000000 34 0.9167184 0.2748228 74 0.8624678 0.3001650 101 0.8526700 0.0000000 34 0.9167184 UCL Min Max Q25 Q50 Q75 1.0832816 1 1 0.9753700 1 1 0.9493102 1 1 1.0832816 1 1 $groups song_chet groups CT 1.0000000 a DC 1.0000000 a CT 0.9189189 ab CT 0.9009901 b attr(,"class") [1] "group" > sk=aov(sk_tong~CTTN,data=LT13) > summary(sk) Df Sum Sq Mean Sq F value Pr(>F) CTTN 5.3 1.778 0.339 0.797 Residuals 223 1169.7 5.245 16 observations deleted due to missingness 23 > sk1=LSD.test(sk,"CTTN") > sk1 $statistics MSerror Df Mean CV 5.245419 223 5.482555 41.7741 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means sk_tong CT 5.545000 CT 5.693529 CT 5.354835 DC 5.340000 std 2.316262 2.016537 2.311448 2.693313 r 34 68 91 34 LCL 4.770962 5.146202 4.881705 4.565962 UCL 6.319038 6.240857 5.827966 6.114038 Min 2.28 1.84 1.07 1.40 Max 12.97 11.21 15.12 12.31 Q25 4.0500 4.1325 3.6100 3.2525 Q50 4.930 5.260 5.040 5.095 Q75 6.6075 7.2125 6.4700 6.4150 $groups sk_tong groups CT 5.693529 a CT 5.545000 a CT 5.354835 a DC 5.340000 a attr(,"class") [1] "group" > fix(LT13) > v=aov(litter_fall~CTTN,data=LT13) > summary(v) Df Sum Sq Mean Sq F value Pr(>F) CTTN 0.00134 0.0004463 0.363 0.78 Residuals 223 0.27390 0.0012282 16 observations deleted due to missingness > v1=LSD.test(v,"CTTN") > v1 $statistics MSerror Df Mean CV 0.001228234 223 0.1254185 27.94339 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means litter_fall CT 0.1261765 CT 0.1288235 CT 0.1232967 DC 0.1235294 std 0.03481681 0.03098047 0.03537278 0.04155189 r 34 68 91 34 LCL 0.1143321 0.1204483 0.1160568 0.1116850 UCL 0.1380209 0.1371988 0.1305366 0.1353738 Min 0.08 0.07 0.06 0.06 $groups litter_fall groups CT 0.1288235 a CT 0.1261765 a DC 0.1235294 a CT 0.1232967 a attr(,"class") [1] "group" > v=aov(litter_fall~CTTN,data=LT23) > summary(v) Df Sum Sq Mean Sq F value Pr(>F) CTTN 0.0865 0.02884 9.273 7.43e-06 *** Residuals 268 0.8335 0.00311 Max 0.24 0.21 0.27 0.23 Q25 0.1000 0.1075 0.1000 0.0900 Q50 0.12 0.12 0.12 0.12 Q75 0.1400 0.1525 0.1400 0.1400 24 Signif codes: ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 27 observations deleted due to missingness > v1=LSD.test(v,"CTTN") > v1 $statistics MSerror Df Mean CV 0.00310996 268 0.2072794 26.90427 $parameters test p.ajusted name.t ntr alpha Fisher-LSD none CTTN 0.05 $means litter_fall CT 0.1867164 CT 0.1992500 CT 0.2338667 DC 0.2078000 std 0.05214746 0.05566116 0.06573685 0.04272790 r 67 80 75 50 LCL 0.1733026 0.1869743 0.2211884 0.1922723 UCL 0.2001303 0.2115257 0.2465450 0.2233277 Min 0.07 0.11 0.10 0.12 Max 0.33 0.34 0.40 0.31 Q25 0.15 0.16 0.19 0.18 Q50 0.18 0.19 0.23 0.20 $groups litter_fall groups CT 0.2338667 a DC 0.2078000 b CT 0.1992500 bc CT 0.1867164 c attr(,"class") [1] "group" III PHÂN TÍCH KHẢ NĂNG CHẮN GIĨ CỦA CÁC ĐAI RỪNG > md=escalc(n1i = N,n2i = N,m1i=Vt,m2i = Vs,sd1i=st,sd2i=ss,data=CG,measure = " SMD",append=TRUE) > summary(md) Models Year Dairung N Vt st Vs ss E se K sk Que 2008 Bach dan 3.10 0.00 2.70 0.00 0.13 0.00 1.10 0.00 Que 2008 Keo chiu han 3.35 0.21 2.90 0.14 0.14 0.01 1.15 0.07 Thuyet & Hung 2005 Keo chiu han 5.50 0.00 1.27 0.21 0.77 0.04 4.43 0.67 LA 2018 Keo la liem 116 2.62 0.53 1.24 0.77 0.53 0.27 4.00 4.96 Lieu 2017 Keo la liem 7.16 0.11 5.53 0.81 0.23 0.11 1.31 0.19 Que 2008 Keo la tram 3.40 0.17 2.65 0.20 0.22 0.04 1.28 0.04 Thuyet 2004 Keo la tram 14 2.19 1.53 0.78 0.22 0.49 0.28 3.02 2.55 Que 2008 Klt + Kll 3.70 0.00 2.50 0.00 0.32 0.00 1.50 0.00 Que 2008 Keo tai tuong 3.00 0.00 2.20 0.00 0.27 0.00 1.40 0.00 10 Que 2008 Phi lao 3.40 0.29 2.35 0.29 0.31 0.03 1.48 0.05 11 Thuyet 2004 Phi lao 41 4.84 1.28 2.18 0.91 0.54 0.17 2.51 1.03 12 Thuyet & Hung 2005 Phi lao 5.50 0.00 1.60 0.00 0.71 0.00 3.40 0.00 13 Que 2008 Xoan AD 3.40 0.14 2.75 0.07 0.19 0.01 1.25 0.07 yi vi sei zi pval ci.lb ci.ub NA NA NA NA NA NA NA 1.4226 1.2530 1.1194 1.2709 0.2038 -0.7713 3.6165 22.7288 43.7165 6.6118 3.4376 0.0006 9.7698 35.6878 2.0810 0.0266 0.1630 12.7655 summary(md) Models Year Que; Bach dan Que; Keo chiu han Thuyet & Hung; Keo chiu han LA; Keo la liem Lieu; Keo la liem Que; Keo la tram Thuyet; Keo la tram Que; Klt+Kll Que; Keo tai tuong 10 Que; Phi lao 11 Thuyet; Phi lao 12 Thuyet & Hung; Phi lao 13 Que; Xoan AD E1 se1 K sk K1 sk1 0 1.10 0.03 0 0 1.15 0.07 0 0 4.43 0.67 0 0 4.00 4.96 0 0 1.31 0.19 0 0 1.28 0.04 0 0 3.02 2.55 0 0 1.50 0.04 0 0 1.40 0.03 0 10 0 1.48 0.05 0 11 0 2.51 1.03 0 12 0 3.40 0.50 0 13 0 1.25 0.07 0 Dairung N Vt st 2008 Bach dan 3.10 0.10 2.70 0.10 2008 Keo chiu han 3.35 0.21 2.90 0.14 2005 Keo chiu han 5.50 0.10 1.27 0.21 2018 Keo la liem 2.62 0.53 1.24 0.77 2017 Keo la liem 7.16 0.11 5.53 0.81 2008 Keo la tram 3.40 0.17 2.65 0.20 2004 Keo la tram 14 2.19 1.53 0.78 0.22 2008 Klt + Kll 3.70 0.10 2.50 0.10 2008 Keo tai tuong 3.00 0.10 2.20 0.10 2008 Phi lao 3.40 0.29 2.35 0.29 2004 Phi lao 41 4.84 1.28 2.18 0.91 2005 Phi lao 5.50 0.10 1.60 0.10 2008 Xoan AD 3.40 0.14 2.75 0.07 yi vi sei zi pval -10.3725 14.4486 3.8011 -2.7288 0.0064 -11.1704 16.5972 4.0740 -2.7419 0.0061 -21.7213 39.9846 6.3233 -3.4351 0.0006 -2.5616 0.6067 0.7789 -3.2886 0.0010 -2.8158 0.4425 0.6652 -4.2331