Khóa luận Nghiên cứu khả năng hấp thụ carbon của rừng trồng sa mộc (Cunnighamia lanceolata Hook ) ở các tuổi 5,7 tại huyện Mường Khương, tỉnh Lào Cai

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC NÔNG LÂM CAO VĂN VƯN NGHIÊN CỨU KHẢ NĂNG HẤP THỤ CARBON CỦA RỪNG TRỒNG SA MỘC (CUNNIGHAMIA LANCEOLATA HOOK) Ở CÁC TUỔI 5, TẠI HUYỆN MƯỜNG KHƯƠNG, TỈNH LÀO CAI KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC Hệ đào tạo : Chính quy Chuyên ngành : Quản lý tài nguyên rừng Khoa : Lâm Nghiệp Khóa học : 2015 - 2019 Thái Nguyên, năm 2019 ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC NÔNG LÂM CAO VĂN VƯN NGHIÊN CỨU KHẢ NĂNG HẤP THỤ CARBON CỦA RỪNG TRỒNG SA MỘC (CUNNIGHAMIA LANCEOLATA HOOK) Ở CÁC TUỔI 5, TẠI HUYỆN MƯỜNG KHƯƠNG, TỈNH LÀO CAI KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC Hệ đào tạo : Chính quy Chuyên ngành : QLTNR Lớp : K47 - QLTNR Khoa : Lâm Nghiệp Khóa học : 2015 - 2019 Giảng viên hướng dẫn : TS Nguyễn Thị Thu Hiền Thái Nguyên, năm 2019 i LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu khoa học thân Các số liệu kết nghiên cứu điều tra thực địa hồn tồn trung thực, chưa cơng bố tài liệu, có sai sót tơi xin chịu hoàn toàn trách nhiệm ii LỜI CẢM ƠN Được phân công Nhà trường, Khoa Lâm Nghiệp, Trường Đại học Nông Lâm Thái Nguyên, sau tháng thực tập em hồn thành khóa luận tốt nghiệp “Nghiên cứu khả hấp thụ carbon rừng trồng Sa mộc (Cunnighamia lanceolata Hook) tuổi 5,7 huyện Mường Khương, tỉnh Lào Cai” Để hoàn thành nhiệm vụ giao, nỗ lực học hỏi thân cịn có hướng dẫn tận tình thầy cô, bạn bè, cô chú, anh chị địa bàn huyện Mường Khương, tỉnh Lào Cai Xin cảm ơn cô giáo – TS Nguyễn Thị Thu Hiền, người hướng dẫn cho em suốt thời gian thực tập dẫn, định hướng cho em để em hoàn thành tốt nhiệm vụ Xin cảm ơn giúp đỡ cán Hạt Kiểm lâm huyện Mường Khương đồng chí kiểm lâm địa bàn thuộc xã huyện Mường Khương Đồng thời xin cảm ơn thầy cô giáo Khoa Lâm Nghiệp, phận liên quan thuộc Trường người thân gia đình bạn bè thân thiết giúp đỡ tơi q trình hồn thành khóa luận tốt nghiệp Mặc dù nỗ lực cố gắng, thuy nhiên hạn chế mặt kinh nghiệm thời gian trình độ nghiên cứu nên khóa luận khơng thể tránh khỏi thiếu sót Tơi mong nhận ý kiến đóng góp, bảo thầy cô bạn đọc khác để khóa luận hồn thiện Tơi xin chân thành cảm ơn! Thái Nguyên, ngày 30 tháng năm 2019 Sinh viên Cao Văn Vưn iii DANH MỤC CÁC BẢNG Bảng 4.1: Các tiêu sinh trưởng loài Sa mộc tuổi 27 Bảng 4.2 Thông tin sinh trưởng mẫu 28 Bảng 4.3 Kết cấu sinh khối tươi cá thể Sa mộc 29 Bảng 4.4 Cấu trúc sinh khối khô cá lẻ rừng trồng Sa mộc 32 Bảng 4.5 Bảng tính tốn sinh khối tươi, khô lâm phần Sa mộc tuổi tuổi 34 Bảng 4.6 Trữ lượng carbon cá thể Sa mộc tuổi 35 Bảng 4.7 Tổng trữ lượng Carbon hấp thụ lâm phần Sa mộc tuổi 36 Bảng 4.8 Lượng hóa giá trị thương mại từ tiêu CO2 tính cho rừng trồng Sa mộc tuổi 37 Bảng 4.9 Bảng tổng hợp kết lựa chọn dạng tương quan D1.3/SKT Sa mộc tuổi khu vực nghiên cứu 39 Bảng 4.10 Tổng hợp kết lựa chọn dạng tương quan D1.3/SKK Sa mộc tuổi khu vực nghiên cứu…………………………………… 40 Bảng 4.11 Tổng hợp kết lựa chọn dạng tương quan D1.3/carbon Sa mộc tuổi khu vực nghiên cứu 41 iv DANH MỤC CÁC HÌNH Hình 4.1: Tỷ lệ sinh khối tươi phận Sa mộc tuổi 30 Hình 4.2: Tỷ lệ sinh khối tươi phận Sa mộc tuổi 30 Hình 4.3: Cân thân tiêu chuẩn 31 Hình 4.4: Cân cành tiêu chuẩn 31 Hình 4.5: Tỷ lệ sinh khối khơ phận Sa mộc tuổi 32 Hình 4.6: Tỷ lệ sinh khối khơ phận Sa mộc tuổi 32 v MỤC LỤC LỜI CẢM ƠN ii DANH MỤC CÁC BẢNG iii DANH MỤC CÁC HÌNH iv MỤC LỤC v PHẦN MỞ ĐẦU 1.1 Đặt vấn đề 1.2 Mục tiêu yêu cầu đề tài PHẦN PHẦN TỔNG QUAN NGHIÊN CỨU 2.1 Cơ sở khoa học vấn đề nghiên cứu 2.2 Tình hình nghiên cứu Thế giới 2.2.1 Nghiên cứu sinh khối rừng trồng 2.2.2 Nghiên cứu khả tích lũy bon rừng trồng 2.3 Tình hình nghiên cứu nước 2.3.1 Nghiên cứu sinh khối rừng trồng 2.3.2 Nghiên cứu khả tích lũy bon rừng trồng 2.4 Tổng quan khu vực nghiên cứu 11 2.4.2 Điều kiện kinh tế - xã hội khu vực nghiên cứu 15 PHẦN ĐỐI TƯỢNG, PHẠM VI, NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 18 3.1 Đối tượng phạm vi nghiên cứu 18 3.2 Giới hạn đề tài 18 3.3 Nội dung nghiên cứu 18 3.4 Phương pháp nghiên cứu 19 3.4.1 Cơ sở phương phương pháp luận 19 3.4.2 Phương pháp thu thập số liệu 19 vi PHẦN KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 27 4.1 Kết đo đếm tiêu sinh trưởng lựa chọn mẫu 27 4.2 Nghiên cứu sinh khối cá thể lâm phần 28 4.2.1 Nghiên cứu kết cấu sinh khối tươi cá thể 29 4.2.2 Sinh khối khô cá thể 31 4.2.3 Nghiên cứu tổng sinh khối toàn lâm phần 33 4.3 Nghiên cứu trữ lượng bon cá thể lâm phần 34 4.3.1 Nghiên cứu trữ lượng Carbon hấp thụ cá thể Sa mộc 34 4.3.2 Nghiên cứu trữ lượng Carbon hấp thụ lâm phần Sa mộc 35 4.4 Lượng hóa giá trị thương mại từ hấp thụ CO2 rừng Sa mộc tuổi tuổi 36 4.5 Xây dựng mối tương quan sinh khối, trữ lượng carbon với tiêu sinh trưởng D1.3 38 4.5.1 Mối tương quan D1.3 tổng sinh khối tươi cá thể 38 4.5.2 Mối tương quan D1.3 với tổng sinh khối khô cá thể 40 4.5.3 Mối tương quan trữ lượng Carbon tích lũy cá thể với nhân tố điều tra lâm phần D1.3 41 PHẦN KẾT LUẬN VÀ KIẾN NGHỊ 43 5.1 Kết luận 43 5.2 Tồn 44 5.3 Kiến nghị 45 TÀI LIỆU THAM KHẢO 46 PHẦN MỞ ĐẦU 1.1 Đặt vấn đề Trên trái đất vấn đề quan tâm Biến đổi khí hậu Biến đổi khí hậu tồn cầu gia tăng nồng độ CO2 Trong khí gây nên tác động tiêu cực môi trường sống nước biển dâng, hạn hán, bão, lũ lụt Ở Việt Nam năm gần xuất nhiều dạng thời tiết cực đoan như: Hiện tượng thời tiết el nino, lũ quét, lũ ống, bão xuất với mật độ ngày dày Gây ảnh hưởng lớn đến đời sống người dân,thiệt hại lớn cho nhiều vùng nhiều quốc gia Biến đổi khí hậu ảnh hưởng chung đến tất quốc gia giới nên biến đổi khí hậu trở thành mối quan tâm hàng đầu với hàng loạt kiện như: phê chuẩn công ước biến đổi khí hậu (năm 1994) ký nghị định thư Kyoto (năm 1997) thành lập ban tư vấn điều hành quốc gia CDM (năm 2003) Những nghiên cứu nước nước khẳng định biến đổi khí hậu ảnh hưởng đến vùng biển nước ta Mực nước biển dâng làm cân sinh thái bị tác động mạnh, Việt Nam nước đứng thứ 10 nước chịu ảnh hưởng nhiều nước biển dâng Việt Nam có hoạt động tích cực để góp phần ngăn chặn nóng lên trái đất phạm vi tồn cầu, quan tâm đến việc trồng rừng, phủ xanh đất trống đồi núi trọc như: Chương trình 327, chương trình trồng triệu rừng 661, chương trình PAM nhiều chương trình bảo tồn khác nhà nước, tổ chức, cá nhân đầu tư xây dựng phát triển Vào đầu tháng năm 2018, Bộ NN&PTNT Việt Nam cho công bố trạng rừng tồn quốc tính đến ngày 31/12/2017 tổng diện tích rừng 14.415.381 ha, rừng tự nhiên 10.236.451 ha, rừng trồng 4.178.966 độ che phủ tương ứng 41,45% Có chế phát triển CDM cho phép nước phát triển đạt mục tiêu giảm phát thải bắt buộc họ thông qua đầu tư thương mại dự án trồng rừng nước phát triển nhằm hấp thụ khí CO2 từ khí quyển.Tuy nhiên việc trồng rừng nhằm hấp thụ CO2 theo chế phát triển CDM việc nghiên cứu định hướng giá trị lợi ích rừng mơi trường, đặc biệt khả hấp thụ CO2 vấn đề Việt Nam Để xác định lượng carbon hấp thụ địi hỏi phải có nghiên cứu định lượng khả hấp thụ CO2 làm sở xác định giá trị thương mại CO2 việc làm cần thiết Để góp phần giải vấn đề tiến hành nghiên cứu đề tài: “Nghiên cứu khả hấp thụ carbon rừng trồng sa mộc (Cunnighamia lanceolata Hook ) tuổi 5,7 huyện Mường Khương, tỉnh Lào Cai” 1.2 Mục tiêu yêu cầu đề tài - Xác định sinh khối lượng carbon tích lũy rừng trồng Sa mộc tuổi 5,7 huyện Mường Khương, tỉnh Lào Cai - Một số ứng dụng xác định sinh khối carbon tích lũy rừng trồng Sa mộc tuổi 5,7 huyện Mường Khương, tỉnh Lào Cai Sinh khối khô tuổi a The model requires all non-missing values to be positive Case Processing Summary N Total Cases Excluded Casesa Forecasted Cases Newly Created Cases 0 a Cases with a missing value in any variable are excluded from the analysis Variable Processing Summary Variables Dependent Independent SKK D1.3 Number of Positive Values Number of Zeros Number of Negative Values Number of Missing Values User-Missing 0 0 System-Missing 0 SKK Linear Model Summary R R Square 848 Adjusted R Square 719 Std Error of the Estimate 649 372 The independent variable is D1.3 ANOVA Sum of Squares Regression Residual Total df Mean Square 1.418 1.418 554 138 1.972 F Sig 10.243 033 The independent variable is D1.3 Coefficients Unstandardized Coefficients B D1.3 (Constant) Std Error 701 219 1.684 1.439 Standardized Coefficients t Sig Beta 848 3.201 033 1.170 307 Logarithmic Model Summary R R Square 850 Adjusted R Square 723 Std Error of the Estimate 654 369 The independent variable is D1.3 ANOVA Sum of Squares Regression Residual Total df Mean Square 1.426 1.426 546 137 1.972 F Sig 10.443 032 The independent variable is D1.3 Coefficients Unstandardized Coefficients B Std Error ln(D1.3) (Constant) Standardized Coefficients t Sig Beta 4.573 1.415 -2.294 2.653 850 3.231 032 -.865 436 Compound Model Summary R R Square 833 Adjusted R Square 694 Std Error of the Estimate 618 064 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 038 038 Residual 017 004 Total 054 F Sig 9.090 039 The independent variable is D1.3 Coefficients Unstandardized Coefficients B Std Error D1.3 1.121 042 (Constant) 2.957 736 The dependent variable is ln(SKK) Standardized Coefficients t Sig Beta 2.301 26.392 000 4.016 016 Power Model Summary R R Square 836 Adjusted R Square 699 Std Error of the Estimate 623 064 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 038 038 Residual 016 004 Total 054 F Sig 9.267 038 The independent variable is D1.3 Coefficients Unstandardized Coefficients B ln(D1.3) (Constant) Std Error 746 245 1.545 709 The dependent variable is ln(SKK) Standardized Coefficients t Sig Beta 836 3.044 038 2.178 095 Carbon a The model requires all non-missing values to be positive Case Processing Summary N Total Cases Excluded Casesa Forecasted Cases Newly Created Cases 0 a Cases with a missing value in any variable are excluded from the analysis Variable Processing Summary Variables Dependent Independent CARBON D1.3 Number of Positive Values Number of Zeros Number of Negative Values Number of Missing Values User-Missing 0 0 System-Missing 0 CARBON Linear Model Summary R R Square 908 Adjusted R Square 825 Std Error of the Estimate 781 112 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 237 237 Residual 050 013 Total 287 F Sig 18.875 012 The independent variable is D1.3 Coefficients Unstandardized Coefficients B D1.3 (Constant) Std Error 286 066 1.328 433 Standardized Coefficients t Sig Beta 908 4.345 012 3.065 037 Logarithmic Model Summary R R Square 907 Adjusted R Square 822 Std Error of the Estimate 777 113 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 236 236 Residual 051 013 Total 287 F Sig 18.444 013 The independent variable is D1.3 Coefficients Unstandardized Coefficients B Standardized Coefficients Std Error t Sig Beta ln(D1.3) 1.860 433 (Constant) -.281 812 907 4.295 013 -.347 746 Compound Model Summary R R Square 912 Adjusted R Square 831 Std Error of the Estimate 789 034 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 023 023 Residual 005 001 Total 028 F Sig 19.677 011 The independent variable is D1.3 Coefficients Unstandardized Coefficients B Std Error D1.3 1.093 022 (Constant) 1.780 236 The dependent variable is ln(CARBON) Standardized Coefficients t Sig Beta 2.488 49.638 000 7.553 002 Power Model Summary R R Square 910 Adjusted R Square 828 Std Error of the Estimate 785 035 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 023 023 Residual 005 001 Total 028 F Sig 19.234 012 The independent variable is D1.3 Coefficients Unstandardized Coefficients B ln(D1.3) (Constant) Std Error 580 132 1.078 267 The dependent variable is ln(CARBON) Standardized Coefficients t Sig Beta 910 4.386 012 4.031 016 Sinh khối khô tuổi a The model requires all non-missing values to be positive Case Processing Summary N Total Cases Excluded Casesa Forecasted Cases Newly Created Cases 0 a Cases with a missing value in any variable are excluded from the analysis Variable Processing Summary Variables Dependent Independent SKK D1.3 Number of Positive Values Number of Zeros Number of Negative Values Number of Missing Values User-Missing 0 0 System-Missing 0 SKK Linear Model Summary R R Square 979 Adjusted R Square 959 Std Error of the Estimate 949 102 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 989 989 Residual 042 010 1.031 Total F Sig 94.394 001 The independent variable is D1.3 Coefficients Unstandardized Coefficients B Std Error D1.3 (Constant) Standardized Coefficients 249 -6.796 2.167 979 Model Summary R Square 978 957 Adjusted R Square 946 Sig Beta 2.424 Logarithmic R t Std Error of the Estimate 105 9.716 001 -3.137 035 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 987 987 Residual 044 011 1.031 Total F Sig 89.435 001 The independent variable is D1.3 Coefficients Unstandardized Coefficients B ln(D1.3) (Constant) Standardized Coefficients Std Error t Sig Beta 20.941 2.214 -31.008 4.786 978 9.457 001 -6.479 003 Compound Model Summary R R Square 981 Adjusted R Square 961 Std Error of the Estimate 952 007 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 005 005 Residual 000 000 Total 005 F Sig 99.753 001 The independent variable is D1.3 Coefficients Unstandardized Coefficients B Standardized Coefficients Std Error 1.186 020 (Constant) 3.239 480 2.666 The dependent variable is ln(SKK) Power Model Summary R Square 980 Adjusted R Square 959 The independent variable is D1.3 .949 Sig Beta D1.3 R t Std Error of the Estimate 007 58.553 000 6.742 003 ANOVA Sum of Squares df Mean Square Regression 005 005 Residual 000 000 Total 005 F Sig 94.590 001 The independent variable is D1.3 Coefficients Unstandardized Coefficients B ln(D1.3) (Constant) Std Error 1.474 152 589 193 The dependent variable is ln(SKK) Standardized Coefficients t Sig Beta 980 9.726 001 3.053 038 Sinh khối tươi tuổi a The model requires all non-missing values to be positive Case Processing Summary N Total Cases Excluded Casesa Forecasted Cases Newly Created Cases 0 a Cases with a missing value in any variable are excluded from the analysis Variable Processing Summary Variables Dependent Independent SKT D1.3 Number of Positive Values Number of Zeros Number of Negative Values Number of Missing Values User-Missing 0 0 System-Missing 0 SKT Linear Model Summary R R Square 910 Adjusted R Square 829 Std Error of the Estimate 786 494 The independent variable is D1.3 ANOVA Sum of Squares Regression Residual Total df Mean Square 4.732 4.732 977 244 5.709 F Sig 19.379 012 The independent variable is D1.3 Coefficients Unstandardized Coefficients B D1.3 (Constant) Std Error 5.302 1.204 -15.964 10.460 Standardized Coefficients t Sig Beta 910 4.402 012 -1.526 202 Logarithmic Model Summary R R Square 911 Adjusted R Square 831 Std Error of the Estimate 789 491 The independent variable is D1.3 ANOVA Sum of Squares Regression Residual Total df Mean Square 4.743 4.743 966 241 5.709 F Sig 19.642 011 The independent variable is D1.3 Coefficients Unstandardized Coefficients B ln(D1.3) (Constant) Standardized Coefficients Std Error t Sig Beta 45.913 10.360 -69.153 22.390 911 4.432 011 -3.089 037 Compound Model Summary R R Square 914 Adjusted R Square 836 Std Error of the Estimate 795 016 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 005 005 Residual 001 000 Total 006 F Sig 20.384 011 The independent variable is D1.3 Coefficients Unstandardized Coefficients B Std Error D1.3 1.194 047 (Constant) 6.427 2.196 The dependent variable is ln(SKT) Standardized Coefficients t Sig Beta 2.495 25.414 000 2.926 043 Power Model Summary R R Square 915 Adjusted R Square 838 Std Error of the Estimate 797 016 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 005 005 Residual 001 000 Total 006 F Sig 20.691 010 The independent variable is D1.3 Coefficients Unstandardized Coefficients B Std Error ln(D1.3) 1.539 338 (Constant) 1.081 790 The dependent variable is ln(SKT) Standardized Coefficients t Sig Beta 915 4.549 010 1.368 243 Carbon a The model requires all non-missing values to be positive Case Processing Summary N Total Cases Excluded Casesa Forecasted Cases Newly Created Cases 0 a Cases with a missing value in any variable are excluded from the analysis Variable Processing Summary Variables Dependent Independent CARBON D1.3 Number of Positive Values Number of Zeros Number of Negative Values Number of Missing Values User-Missing 0 0 System-Missing 0 CARBON Linear Model Summary R R Square 953 Adjusted R Square 907 Std Error of the Estimate 884 107 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 450 450 Residual 046 011 Total 496 F Sig 39.215 003 The independent variable is D1.3 Coefficients Unstandardized Coefficients B Std Error D1.3 (Constant) Standardized Coefficients 261 -7.103 2.267 953 Model Summary R Square 950 903 Adjusted R Square 879 Sig Beta 1.635 Logarithmic R t Std Error of the Estimate 110 6.262 003 -3.133 035 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 448 448 Residual 048 012 Total 496 F Sig 37.200 004 The independent variable is D1.3 Coefficients Unstandardized Coefficients B ln(D1.3) (Constant) Standardized Coefficients Std Error t Sig Beta 14.104 2.312 -23.390 4.998 950 6.099 004 -4.680 009 Compound Model Summary R R Square 954 Adjusted R Square 910 Std Error of the Estimate 888 015 The independent variable is D1.3 ANOVA Sum of Squares df Mean Square Regression 009 009 Residual 001 000 Total 010 F Sig 40.672 003 The independent variable is D1.3 Coefficients Unstandardized Coefficients B D1.3 Standardized Coefficients Std Error 045 958 301 2.597 The dependent variable is ln(CARBON) Power Model Summary R R Square 952 Adjusted R Square 906 The independent variable is D1.3 .883 Sig Beta 1.259 (Constant) t Std Error of the Estimate 015 27.679 000 3.187 033 ANOVA Sum of Squares df Mean Square Regression 009 009 Residual 001 000 Total 010 F Sig 38.585 003 The independent variable is D1.3 Coefficients Unstandardized Coefficients B ln(D1.3) (Constant) Std Error 1.988 320 096 067 The dependent variable is ln(CARBON) Standardized Coefficients t Sig Beta 952 6.212 003 1.446 222 ... HỌC NÔNG LÂM CAO VĂN VƯN NGHIÊN CỨU KHẢ NĂNG HẤP THỤ CARBON CỦA RỪNG TRỒNG SA MỘC (CUNNIGHAMIA LANCEOLATA HOOK) Ở CÁC TUỔI 5, TẠI HUYỆN MƯỜNG KHƯƠNG, TỈNH LÀO CAI KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC... lượng carbon tích lũy rừng trồng Sa mộc tuổi 5,7 huyện Mường Khương, tỉnh Lào Cai - Một số ứng dụng xác định sinh khối carbon tích lũy rừng trồng Sa mộc tuổi 5,7 huyện Mường Khương, tỉnh Lào Cai. .. phần giải vấn đề tiến hành nghiên cứu đề tài: ? ?Nghiên cứu khả hấp thụ carbon rừng trồng sa mộc (Cunnighamia lanceolata Hook ) tuổi 5,7 huyện Mường Khương, tỉnh Lào Cai? ?? 1.2 Mục tiêu yêu cầu đề