Nghiên cứu sử dụng cỏ vetiver (vetiveria zizanioides (l ) nash) để cải tạo đất bị ô nhiễm pb, as sau khai thác khoáng sản ở tỉnh thái nguyên

BỘ GIÁO DỤC VÀ ĐÀO TẠO ĐẠI HỌC THÁI NGUYÊN LƢƠNG THỊ THUÝ VÂN NGHIÊN CỨU SỬ DỤNG CỎ VETIVER (VETIVERIA ZIZANIOIDES (L.) NASH) ĐỂ CẢI TẠO ĐẤT BỊ Ô NHIỄM Pb, As SAU KHAI THÁC KHOÁNG SẢN Ở TỈNH THÁI NGUYÊN Chuyên ngành: Trồng trọt Mã số: 62.62.01.01 LUẬN ÁN TIẾN SĨ NÔNG NGHIỆP Người hướng dẫn khoa học: PGS.TS Lƣơng Văn Hinh TS.NCVCC Trần Văn Tựa THÁI NGUYÊN - 2012 Tai ngay!!! Ban co the xoa dong chu nay!!! i LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Các số liệu, kết nêu luận án trung thực Những kết luận luận án chưa công bố cơng trình khác Tác giả luận án Lương Thị Thúy Vân ii LỜI CẢM ƠN * Tôi xin trân trọng cảm ơn tới: - Ban Giám hiệu, Phòng Quản lý đào tạo Sau đại học, khoa Tài nguyên Môi trường trường Đại học Nông lâm - Đại học Thái Nguyên giúp đỡ, tạo điều kiện tốt cho tơi học tập, nghiên cứu hồn thành luận án - Ban Giám hiệu, Ban chủ nhiệm khoa Sinh - KTNN trường Đại học Sư phạm - Đại học Thái Nguyên tạo điều kiện giúp đỡ tơi q trình học tập, nghiên cứu * Với lịng biết ơn chân thành, tơi xin gửi lời cảm ơn sâu sắc tới: - PGS.TS Lương Văn Hinh người thầy hướng nghiên cứu, hướng dẫn tận tình, động viên giúp đỡ bước tơi q trình nghiên cứu hồn thành luận án - TS NCVCC Trần Văn Tựa người thầy trực tiếp hướng dẫn, giúp đỡ nghiên cứu hồn thành luận án - Tơi xin chân thành cảm ơn Bộ Khoa học Công nghệ chủ nhiệm đề tài cấp Nhà nước KC.08.04/06-10 tài trợ kinh phí để thực nghiên cứu - Tôi xin chân thành cảm ơn đồng nghiệp, bạn bè gia đình tạo điều kiện giúp đỡ, động viên, khích lệ tơi q trình nghiên cứu luận án Xin chân thành cảm ơn ! Thái Nguyên, tháng năm 2012 Nghiên cứu sinh Lương Thị Thúy Vân iii Lương Thị Thúy Vân MỤC LỤC Lời cam đoan i Lời cảm ơn ii Mục lục iii Danh mục ký hiệu chữ viết tắt vii Danh mục bảng viii MỞ ĐẦU 1 Tính cấp thiết đề tài Mục tiêu nghiên cứu đề tài Đối tượng phạm vi nghiên cứu 3.1 Đối tượng nghiên cứu 3.2 Phạm vi nghiên cứu Ý nghĩa đề tài 4.1 Ý nghĩa khoa học đề tài 4.2 Ý nghĩa thực tiễn đề tài Những đóng góp đề tài Chƣơng TỔNG QUAN CÁC VẤN ĐỀ NGHIÊN CỨU 1.1 Ô nhiễm kim loại nặng đất 1.1.1 Khái niệm kim loại nặng 1.1.2 Đất ô nhiễm kim loại nặng hoạt động khai thác khoáng sản 1.1.3 Tiêu chuẩn đánh giá mức độ ô nhiễm đất kim loại nặng 1.1.4 Sự tồn chuyển hóa nguyên tố Pb, As đất 12 1.2 Sử dụng thực vật xử lý đất ô nhiễm kim loại nặng 17 1.2.1 Khái niệm chung 17 iv 1.2.2 Các biện pháp sử dụng thực vật xử lý kim loại nặng đất 18 1.2.3 Tiêu chuẩn loài thực vật sử dụng để xử lý kim loại nặng đất 22 1.2.4 Phương pháp xử lý thực vật sau tích lũy chất nhiễm 22 1.2.5 Ưu điểm hạn chế biện pháp sử dụng thực vật xử lý kim loại nặng đất 22 1.2.6 Một số kết nghiên cứu sử dụng thực vật để xử lý đất ô nhiễm kim loại nặng giới Việt Nam 25 1.3 Một số đặc điểm cỏ Vetiver tình hình nghiên cứu sử dụng cỏ Vetiver cải tạo đất ô nhiễm 27 1.3.1 Một số đặc điểm cỏ Vetiver 27 1.3.2 Tình hình nghiên cứu, sử dụng cỏ Veitver cải tạo đất ô nhiễm giới Việt Nam 30 Chƣơng ĐỐI TƢỢNG, NỘI DUNG VÀ PHƢƠNG PHÁP NGHIÊN CỨU 36 2.1 Đối tượng nghiên cứu 36 2.2 Nội dung nghiên cứu 36 2.3 Phương pháp nghiên cứu 36 2.3.1 Phương pháp đánh giá trạng ô nhiễm kim loại nặng đất 36 2.3.2 Phương pháp bố trí thí nghiệm 38 2.3.3 Phương pháp thu mẫu xác định tiêu sinh trưởng 42 2.3.4 Phương pháp nghiên cứu phịng thí nghiệm 43 Chƣơng KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 45 3.1 Tình trạng nhiễm kim loại nặng mơi trường đất q trình khai thác khoáng sản khu vực khai thác thiếc (xã Hà Thượng, huyện Đại Từ) khu vực khai thác chì - kẽm (xã Tân Long, huyện Đồng Hỷ) tỉnh Thái Nguyên 45 v 3.1.1 Khu vực khai thác thiếc xã Hà Thượng - Đại Từ - tỉnh Thái Nguyên 46 3.1.2 Khu vực khai thác quặng Pb - Zn làng Hích, Đồng Hỷ, Thái Nguyên 48 3.2 Nghiên cứu khả chống chịu tích lũy Pb, As cỏ Vetiver trồng đất nhiễm q trình khai thác khống sản 50 3.2.1 Nghiên cứu khả chống chịu tích lũy Pb cỏ Vetiver trồng đất ô nhiễm 51 3.2.2 Nghiên cứu khả sinh trưởng tích lũy As cỏ Vetiver trồng đất ô nhiễm 62 3.3 Nghiên cứu ảnh hưởng phân bón đến khả sinh trưởng khả tích lũy Pb, As cỏ Vetiver trồng đất ô nhiễm 70 3.3.1 Ảnh hưởng phân bón đến khả phân nhánh cỏ Vetiver trồng đất ô nhiễm Pb ô nhiễm As 71 3.3.2 Ảnh hưởng phân bón đến chiều cao thân chiều dài rễ cỏ Vetiver trồng đất ô nhiễm Pb ô nhiễm As 73 3.3.3 Ảnh hưởng phân bón đến sinh khối cỏ Vetiver trồng đất ô nhiễm Pb ô nhiễm As 77 3.3.4 Ảnh hưởng phân bón đến khả tích lũy Pb, As cỏ Vetiver trồng đất ô nhiễm 82 3.4 Nghiên cứu ảnh hưởng mật độ trồng đến sinh trưởng, phát triển tích lũy Pb, As cỏ Vetiver 88 3.4.1 Ảnh hưởng mật độ trồng đến khả đẻ nhánh cỏ Vetiver đất ô nhiễm Pb ô nhiễm As 88 3.4.2 Ảnh hưởng mật độ trồng đến chiều cao thân chiều dài rễ cỏ Vetiver đất ô nhiễm Pb ô nhiễm As 90 3.4.3 Ảnh hưởng mật độ trồng đến sinh khối cỏ Vetiver trồng đất ô nhiễm Pb ô nhiễm As 92 3.4.4 Ảnh hưởng mật độ đến khả tích lũy Pb, As cỏ Vetiver trồng đất ô nhiễm 98 3.5 Nghiên cứu ảnh hưởng chu kỳ thu hoạch đến sinh khối khả tích lũy Pb, As cỏ Vetiver trồng đất ô nhiễm 100 vi 3.5.1 Ảnh hưởng chu kỳ thu hoạch đến sinh khối cỏ Vetiver trồng đất ô nhiễm Pb ô nhiễm As 100 3.5.2 Ảnh hưởng chu kỳ thu hoạch đến khả tích lũy Pb, As cỏ Vetiver 102 3.6 So sánh tính chất đất nhiễm Pb, As trước sau trồng cỏ Vetiver 104 3.7 Đề xuất biện pháp kỹ thuật trồng cỏ Vetiver cải tạo đất ô nhiễm kim loại nặng trình khai thác khống sản 105 3.7.1 Thời vụ trồng cỏ 105 3.7.2 Làm đất 105 3.7.3 Mật độ trồng, phương pháp trồng 106 3.7.4 Chăm sóc 106 3.7.5 Bón phân 107 3.7.6 Các lứa cắt 107 3.7.7 Xử lý sinh khối 108 KẾT LUẬN VÀ ĐỀ NGHỊ 110 Kết luận 110 Đề nghị 111 TÀI LIỆU THAM KHẢO 112 DANH MỤC CÁC CƠNG TRÌNH ĐÃ ĐƢỢC CƠNG BỐ CỦA TÁC GIẢ CÓ LIÊN QUAN ĐẾN LUẬN ÁN 123 PHỤ LỤC 124 vii DANH MỤC KÝ HIỆU CÁC CHỮ VIẾT TẮT Bộ KH CNMT: Bộ Khoa học Công nghệ môi trường BOD: Biochemical Oxygen Demand (chỉ số nhu cầu oxy sinh hóa) CEC: Cation Exchange Capacity (Dung tích trao đổi cation) COD: Chemical Oxygen Demand (chỉ số nhu cầu oxy hóa học) cs: cộng DAP: Diamino phosphate DW: Dry weight (Khối lượng khô) EEA: European Environment Agency (Cục môi trường châu Âu) FW: Fresh weight (Khối lượng tươi) HT: Hà Thượng ICP-MS: Inductively-Coupled Plasma-Mass Spectrometry (Phương pháp khối phổ plasma cảm ứng) ppm: past per million (Nồng độ phần triệu) ppb: past per billion (Nồng độ phần tỷ) QCVN: Quy chuẩn Việt Nam SAS: Statistical Analysis System (Phần mềm phân tích thống kê) SKK: Sinh khối khơ TCE: Tricloroetylen TCVN: Tiêu chuẩn Việt Nam TKV: Tập đoàn Than Khoáng sản Việt Nam TL: Tân Long TNHH NN: Trách nhiệm hữu hạn Nhà nước viii DANH MỤC BẢNG Bảng 1.1 Hàm lượng kim loại nặng số loại đất khu mỏ hoang Songcheon Bảng 1.2 Hàm lượng kim loại nặng chất thải số mỏ vàng điển hình Úc Bảng 1.3 Đánh giá mức ô nhiễm kim loại đất Hà Lan Bảng 1.4 Hàm lượng tối đa cho phép kim loại nặng xem độc thực vật đất nông nghiệp Bảng 1.5 Đánh giá ô nhiễm đất mặt kim loại nặng Ba Lan 10 Bảng 1.6 Giới hạn tối đa cho phép hàm lượng tổng số As, Cd, Cu, Pb Zn đất (tầng đất mặt) 11 Bảng 1.7 Khả linh động số nguyên tố kim loại nặng đất 12 Bảng 1.8 Hàm lượng As thực phẩm trồng 16 Bảng 1.9 So sánh ngưỡng chịu kim loại nặng cỏ Vetiver loài cỏ khác 30 Bảng 2.1 Vị trí đặc điểm điểm lấy mẫu 37 Bảng 2.2 Một số tính chất đất thí nghiệm chậu 38 Bảng 2.3 Cơng thức thí nghiệm nghiên cứu khả chống chịu tích lũy Pb, As cỏ Vetiver trồng đất ô nhiễm 39 Bảng 2.4 Một số tính chất đất thí nghiệm ngồi đồng ruộng khu vực nghiên cứu 40 Bảng 2.5 Cơng thức thí nghiệm nghiên cứu ảnh hưởng phân bón đến sinh trưởng, phát triển khả tích lũy Pb, As cỏ Vetiver trồng đất ô nhiễm 40 Bảng 2.6 Cơng thức thí nghiệm nghiên cứu ảnh hưởng mật độ trồng đến sinh trưởng, phát triển khả tích lũy Pb, As cỏ Vetiver trồng đất ô nhiễm 41 ix Bảng 3.1 Hàm lượng kim loại nặng đất khu vực khai thác thiếc xã Hà Thượng - Đại Từ - tỉnh Thái Nguyên 47 Bảng 3.2 Hàm lượng kim loại nặng đất khu vực khai thác quặng Pb - Zn xã Tân Long, Đồng Hỷ, Thái Nguyên 50 Bảng 3.3 Ảnh hưởng hàm lượng Pb đất đến phân nhánh cỏ Vetiver 52 Bảng 3.4 Ảnh hưởng hàm lượng Pb đất đến chiều cao thân chiều dài rễ cỏ Vetiver 54 Bảng 3.5 Ảnh hưởng hàm lượng Pb đất đến sinh khối cỏ Vetiver 56 Bảng 3.6 Hàm lượng Pb thân cỏ Vetiver giai đoạn sinh trưởng khác 58 Bảng 3.7 Hàm lượng Pb rễ cỏ Vetiver giai đoạn sinh trưởng khác 59 Bảng 3.8 Ảnh hưởng hàm lượng As đất đến đến khả phân nhánh cỏ Vetiver 62 Bảng 3.9 Ảnh hưởng hàm lượng As đất đến chiều cao thân chiều dài rễ cỏ Vetiver 64 Bảng 3.10 Ảnh hưởng hàm lượng As đất đến tạo thành sinh khối cỏ Vetiver 65 Bảng 3.11 Hàm lượng As thân cỏ Vetiver giai đoạn sinh trưởng khác 66 Bảng 3.12 Hàm lượng As rễ cỏ Vetiver giai đoạn sinh trưởng khác 67 Bảng 3.13 Ảnh hưởng phân bón đến khả phân nhánh cỏ Vetiver trồng đất ô nhiễm Pb 71 Bảng 3.14 Ảnh hưởng phân bón đến phân nhánh cỏ Vetiver trồng đất ô nhiễm As 72 Bảng 3.15 Ảnh hưởng phân bón đến chiều cao thân cỏ Vetiver trồng đất ô nhiễm Pb 73 162 B B B B B B 5.3200 3 5.1667 4.9067 The SAS System Analysis of Variance Procedure Dependent Variable: HÀM LƯỢNG AS TRONG THÂN LÁ CỎ VETIVER (150 NGÀY) Source DF Sum of Squares Mean Square F Value Pr > F Model 11 46.08990000 4.18999091 27.05 0.0001 Error 18 2.78858000 0.15492111 Corrected Total 29 48.87848000 R-Square C.V Root MSE YIELD17 Mean 0.942949 4.523100 0.39360019 8.70200000 Source DF Anova SS Mean Square F Value Pr > F REP 0.51062000 0.25531000 1.65 0.0202 TRT 45.57928000 5.06436444 32.69 0.0001 T tests (LSD) for variable: HÀM LƯỢNG AS TRONG THÂN LÁ CỎ VETIVER (150 NGÀY) NOTE: This test controls the type I comparisonwise error rate not the experimentwise error rate Alpha= 0.05 df= 18 MSE= 0.154921 Critical Value of T= 2.10 Least Significant Difference= 0.6752 Means with the same letter are not significantly different T Grouping Mean N TRT A 10.4767 A A 10.4733 A A 10.4500 10 B 8.5900 3 B B 8.5700 B B 8.4467 C 7.5633 C C 7.5300 C C 7.4667 C C 7.4533 163 The SAS System Analysis of Variance Procedure Dependent Variable: HÀM LƯỢNG AS TRONG THÂN LÁ CỎ VETIVER (210 NGÀY) Source DF Sum of Squares Mean Square F Value Pr > F Model 11 47.72888000 4.33898909 4.16 0.0037 Error 18 18.78064000 1.04336889 Corrected Total 29 66.50952000 R-Square C.V Root MSE YIELD18 Mean 0.717625 9.596527 1.02145430 10.64400000 Source DF Anova SS Mean Square F Value Pr > F REP 0.14496000 0.07248000 0.07 0.0031 TRT 47.58392000 5.28710222 5.07 0.0017 T tests (LSD) for variable: HÀM LƯỢNG AS TRONG THÂN LÁ CỎ VETIVER (210 NGÀY) NOTE: This test controls the type I comparisonwise error rate not the experimentwise error rate Alpha= 0.05 df= 18 MSE= 1.043369 Critical Value of T= 2.10 Least Significant Difference= 1.7522 Means with the same letter are not significantly different T Grouping Mean N TRT A 12.5533 10 A A 12.4467 A A 12.3233 B 10.5133 B B 10.4400 3 B B 10.4367 B B 9.6567 B B 9.6067 B B 9.2767 B B 9.1867 The SAS System Analysis of Variance Procedure Dependent Variable: HÀM LƯỢNG As TRONG RỄ CỎ VETIVER (90 NGÀY) Source DF Sum of Squares Mean Square F Value Pr > F Model 11 605.46720333 55.04247303 2.72 0.0291 Error 18 364.80049333 20.26669407 Corrected Total 29 970.26769667 R-Square C.V Root MSE YIELD1 Mean 0.624021 10.76063 4.50185451 41.83633333 Source DF Anova SS Mean Square F Value Pr > F 164 REP TRT 282.35010667 323.11709667 141.17505333 35.90189963 6.97 1.77 0.0057 0.1443 T tests (LSD) for variable: HÀM LƯỢNG As TRONG RỄ CỎ VETIVER (90 NGÀY) NOTE: This test controls the type I comparisonwise error rate not the experimentwise error rate Alpha= 0.05 df= 18 MSE= 20.26669 Critical Value of T= 2.10 Least Significant Difference= 7.7225 Means with the same letter are not significantly different T Grouping Mean N TRT A 50.757 B 42.347 B B 42.050 10 B B 41.770 B B 41.620 B B 41.567 B B 41.413 3 B B 39.630 B B 39.547 B B 37.663 The SAS System Analysis of Variance Procedure Dependent Variable: HÀM LƯỢNG As TRONG RỄ CỎ VETIVER (150 NGÀY) Source DF Sum of Squares Mean Square F Value Pr > F Model 11 1545.25328333 140.47757121 4.23 0.0034 Error 18 597.57325333 33.19851407 Corrected Total 29 2142.82653667 R-Square C.V Root MSE YIELD5 Mean 0.721128 9.090283 5.76181517 63.38433333 Source DF Anova SS Mean Square F Value Pr > F REP 6.20514667 3.10257333 0.09 0.0412 TRT 1539.04813667 171.00534852 5.15 0.0015 T tests (LSD) for variable: HÀM LƯỢNG As TRONG RỄ CỎ VETIVER (150 NGÀY) NOTE: This test controls the type I comparisonwise error rate not the experimentwise error rate Alpha= 0.05 df= 18 MSE= 33.19851 Critical Value of T= 2.10 Least Significant Difference= 9.8838 Means with the same letter are not significantly different T Grouping Mean N TRT 165 B B B B B B B B B B B D D D D D D D A A A A A A A A A A A E E E E E C C C C C C C 72.520 72.397 68.387 10 67.793 65.277 63.540 3 62.610 56.597 54.663 50.060 The SAS System Analysis of Variance Procedure Dependent Variable: HÀM LƯỢNG As TRONG RỄ CỎ VETIVER (210 NGÀY) Source DF Sum of Squares Mean Square F Value Pr > F Model 11 1293.70199667 117.60927242 2.15 0.0424 Error 18 985.72075333 54.76226407 Corrected Total 29 2279.42275000 R-Square C.V Root MSE YIELD6 Mean 0.567557 5.934603 7.40015298 124.69500000 Source DF Anova SS Mean Square F Value Pr > F REP 33.07398000 16.53699000 0.30 0.0430 TRT 1260.62801667 140.06977963 2.56 0.0430 T tests (LSD) for variable: HÀM LƯỢNG As TRONG RỄ CỎ VETIVER (210 NGÀY) NOTE: This test controls the type I comparisonwise error rate not the experimentwise error rate Alpha= 0.05 df= 18 MSE= 54.76226 Critical Value of T= 2.10 Least Significant Difference= 12.694 Means with the same letter are not significantly different T Grouping Mean N TRT A 139.153 A B A 129.067 10 B A B A C 127.383 B A C B A C 127.127 B C B C 125.443 B C 166 B B B B B C C C C C C C C C 123.590 123.400 3 119.807 116.047 115.933 THÍ NGHIỆM VỀ MẬT ĐỘ TRÊN ĐẤT Ô NHIỄM As The SAS System Analysis of Variance Procedure Class Level Information Class Levels Values REP 3 TRT 4 Number of observations in data set = 12 The SAS System Analysis of Variance Procedure Dependent Variable: SỐ NHÁNH CỎ VETIVER GIAI ĐOẠN 90 NGÀY Source DF Sum of Squares Mean Square F Value Pr > F Model 4.37445000 0.87489000 0.29 0.0107 Error 7.94575000 2.99095833 Corrected Total 11 22.32020000 R-Square C.V Root MSE YIELD1 Mean 0.195986 14.49655 1.72943873 11.93000000 Source DF Anova SS Mean Square F Value Pr > F REP 1.39065000 0.69532500 0.23 0.0194 TRT 2.98380000 0.99460000 0.33 0.8027 Duncan's Multiple Range Test for variable: SỐ NHÁNH CỎ VETIVER GIAI ĐOẠN 90 NGÀY NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 2.990958 Number of Means Critical Range 3.455 3.581 3.643 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 12.777 3 A A 11.777 A A 11.667 A A 11.500 167 The SAS System Analysis of Variance Procedure Dependent Variable: SỐ NHÁNH CỎ VETIVER GIAI ĐOẠN 150 NGÀY Source DF Sum of Squares Mean Square F Value Pr > F Model 29.68880833 5.93776167 0.71 0.0371 Error 50.08108333 8.34684722 Corrected Total 11 79.76989167 R-Square C.V Root MSE YIELD2 Mean 0.372181 13.63154 2.88909107 21.19416667 Source DF Anova SS Mean Square F Value Pr > F REP 4.11231667 2.05615833 0.25 0.0392 TRT 25.57649167 8.52549722 1.02 0.4468 Duncan's Multiple Range Test for variable: SỐ NHÁNH CỎ VETIVER GIAI ĐOẠN 150 NGÀY NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 8.346847 Number of Means Critical Range 5.772 5.982 6.086 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 23.110 A A 21.667 3 A A 20.943 A A 19.057 The SAS System Analysis of Variance Procedure Dependent Variable: SỐ NHÁNH CỎ VETIVER GIAI ĐOẠN 210 NGÀY Source DF Sum of Squares Mean Square F Value Pr > F Model 492.05170833 98.41034167 6.07 0.0242 Error 97.27431667 16.21238611 Corrected Total 11 589.32602500 R-Square C.V Root MSE YIELD3 Mean 0.834940 9.452340 4.02646074 42.59750000 Source DF Anova SS Mean Square F Value Pr > F REP 19.99095000 9.99547500 0.62 0.0208 TRT 472.06075833 157.35358611 9.71 0.0102 Duncan's Multiple Range Test for variable: SỐ NHÁNH CỎ VETIVER GIAI ĐOẠN 210 NGÀY NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 16.21239 Number of Means Critical Range 8.044 8.337 8.483 Means with the same letter are not significantly different 168 Duncan Grouping A A B A B B C C C Mean 50.777 N TRT 45.833 3 39.613 34.167 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 90 NGÀY (GAM/KHÓM) Source DF Sum of Squares Mean Square F Value Pr > F Model 450.45818333 90.09163667 2.55 0.0432 Error 212.30508333 35.38418056 Corrected Total 11 662.76326667 R-Square C.V Root MSE YIELD10 Mean 0.679667 5.514301 5.94846035 107.87333333 Source DF Anova SS Mean Square F Value Pr > F REP 303.48371667 151.74185833 4.29 0.0497 TRT 146.97446667 48.99148889 1.38 0.3351 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 90 NGÀY (GAM/KHÓM) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 35.38418 Number of Means Critical Range 11.88 12.32 12.53 Means with the same letter are not significantly different Duncan Grouping A A A A A A A Mean 112.027 N TRT 110.463 105.573 3 103.430 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 150 NGÀY (GAM/KHÓM) Source DF Sum of Squares Mean Square F Value Pr > F Model 458.59575833 91.71915167 3.13 0.0408 Error 176.00513333 29.33418889 Corrected Total 11 634.60089167 R-Square C.V Root MSE YIELD11 Mean 0.722652 4.265965 5.41610459 126.96083333 Source DF Anova SS Mean Square F Value Pr > F REP 162.60086667 81.30043333 2.77 0.0404 TRT 295.99489167 98.66496389 3.36 0.0961 169 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 150 NGÀY (GAM/KHÓM) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 29.33419 Number of Means Critical Range 10.82 11.21 11.41 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 132.447 3 A A 130.990 A A 124.083 A A 120.323 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 210 NGÀY (GAM/KHÓM) Source DF Sum of Squares Mean Square F Value Pr > F Model 1751.26530833 350.25306167 5.73 0.0277 Error 366.81438333 61.13573056 Corrected Total 11 2118.07969167 R-Square C.V Root MSE YIELD12 Mean 0.826817 5.104881 7.81893411 153.16583333 Source DF Anova SS Mean Square F Value Pr > F REP 107.28521667 53.64260833 0.88 0.0432 TRT 1643.98009167 547.99336389 8.96 0.0123 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 210 NGÀY (GAM/KHÓM) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 61.13573 Number of Means Critical Range 15.62 16.19 16.47 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 166.683 A A 157.287 3 A A 154.163 A 134.530 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 90 NGÀY (GAM/KHÓM) Source DF Sum of Squares Mean Square F Value Pr > F Model 382.38931667 76.47786333 2.59 0.0388 Error 176.98435000 29.49739167 170 Corrected Total 11 559.37366667 R-Square C.V Root MSE YIELD13 Mean 0.683603 5.270744 5.43115012 103.04333333 Source DF Anova SS Mean Square F Value Pr > F REP 266.74291667 133.37145833 4.52 0.0335 TRT 115.64640000 38.54880000 1.31 0.3558 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 90 NGÀY (GAM/KHÓM) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 29.49739 Number of Means Critical Range 10.85 11.25 11.44 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 106.863 3 A A 104.583 A A 102.263 A A 98.463 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 150 NGÀY (GAM/KHÓM) Source DF Sum of Squares Mean Square F Value Pr > F Model 653.33454167 130.66690833 2.72 0.0276 Error 288.25388333 48.04231389 Corrected Total 11 941.58842500 R-Square C.V Root MSE YIELD14 Mean 0.693864 5.470281 6.93125630 126.70750000 Source DF Anova SS Mean Square F Value Pr > F REP 378.12065000 189.06032500 3.94 0.0309 TRT 275.21389167 91.73796389 1.91 0.2292 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 150 NGÀY (GAM/KHÓM) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 48.04231 Number of Means Critical Range 13.85 14.35 14.60 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 133.047 A A 129.120 3 A A 124.260 A A 120.403 171 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 210 NGÀY (GAM/KHÓM) Source DF Sum of Squares Mean Square F Value Pr > F Model 441.07004167 88.21400833 0.84 0.0370 Error 630.88245000 105.14707500 Corrected Total 11 1071.95249167 R-Square C.V Root MSE YIELD15 Mean 0.411464 6.933890 10.25412478 147.88416667 Source DF Anova SS Mean Square F Value Pr > F REP 65.15461667 32.57730833 0.31 0.0446 TRT 375.91542500 125.30514167 1.19 0.3894 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 210 NGÀY (GAM/KHÓM) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 105.1471 Number of Means Critical Range 20.49 21.23 21.60 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 154.280 A A 151.090 3 A A 146.827 A A 139.340 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 90 NGÀY (GAM/M2) Source DF Sum of Squares Mean Square F Value Pr > F Model 576796.94811667 115359.38962333 51.03 0.0001 Error 13562.46298333 2260.41049722 Corrected Total11 590359.41110000 R-Square C.V Root MSE YIELD1 Mean 0.977027 6.571858 47.54377454 723.44500000 Source DF Anova SS Mean Square F Value Pr > F REP 13651.50455000 6825.75227500 3.02 0.0238 TRT 563145.44356667 187715.14785556 83.04 0.0001 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 90 NGÀY (GAM/M2) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 2260.41 Number of Means Critical Range 95.0 98.4 100.2 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 1034.30 B 784.19 172 C D 633.44 3 441.85 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 150 NGÀY (GAM/M2) Source DF Sum of Squares Mean Square F Value Pr > F Model 710262.35257500 142052.47051500 125.08 0.0001 Error 813.91911667 1135.65318611 Corrected Total 11 717076.27169167 R-Square C.V Root MSE YIELD2 Mean 0.990498 3.975801 33.69945380 847.61416667 Source DF Anova SS Mean Square F Value Pr > F REP 6886.12281667 3443.06140833 3.03 0.0230 TRT 703376.22975833 234458.74325278 206.45 0.0001 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 150 NGÀY (GAM/M2) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 1135.653 Number of Means Critical Range 67.33 69.78 71.00 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 1203.23 B 868.58 C 794.68 3 D 523.96 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 210 NGÀY (GAM/M2) Source DF Sum of Squares Mean Square F Value Pr > F Model 727030.35177500 145406.070355 44.08 0.0001 Error 19792.65611667 3298.77601944 Corrected Total 11 746823.00789167 R-Square C.V Root MSE YIELD3 Mean 0.973498 5.693823 57.43497209 1008.72416667 Source DF Anova SS Mean Square F Value Pr > F REP 8754.08461667 4377.04230833 1.33 0.0333 TRT 718276.26715833 239425.42238611 72.58 0.0001 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG THÂN LÁ CỎ VETIVER GIAI ĐOẠN 210 NGÀY (GAM/M2) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 3298.776 Number of Means Critical Range 114.7 118.9 121.0 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 1345.30 B 1079.14 C 943.72 3 D 666.73 173 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 90 NGÀY (GAM/M2) Source DF Sum of Squares Mean Square F Value Pr > F Model 571150.30607500 114230.06121500 55.78 0.0001 Error 12287.52135000 2047.92022500 Corrected Total11 583437.82742500 R-Square C.V Root MSE YIELD4 Mean 0.978939 6.531589 45.25395259 692.84750000 Source DF Anova SS Mean Square F Value Pr > F REP 10752.57365000 5376.28682500 2.63 0.0117 TRT 560397.73242500 186799.24414167 91.21 0.0001 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 90 NGÀY (GAM/M2) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 2047.92 Number of Means Critical Range 90.41 93.71 95.34 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 1022.63 B 689.24 B B 641.18 3 C 418.33 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 150 NGÀY (GAM/M2) Source DF Sum of Squares Mean Square F Value Pr > F Model 716535.54546667 143307.10909333 64.08 0.0001 Error 13418.15173334 2236.35862222 Corrected Total 11 729953.69720000 R-Square C.V Root MSE YIELD5 Mean 0.981618 5.595210 47.29015354 845.19000000 Source DF Anova SS Mean Square F Value Pr > F REP 19598.69840000 9799.34920000 4.38 0.0071 TRT 696936.84706667 232312.28235556 103.88 0.0001 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 150 NGÀY (GAM/M2) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 2236.359 Number of Means Critical Range 94.48 97.92 99.63 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 1204.03 174 B C D 869.82 774.72 532.19 3 3 The SAS System Analysis of Variance Procedure Dependent Variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 210 NGÀY (GAM/M2) Source DF Sum of Squares Mean Square F Value Pr > F Model 938551.51864166 187710.30372833 40.36 0.0002 Error 27903.80725000 4650.63454167 Corrected Total 11 966455.32589166 R-Square C.V Root MSE YIELD6 Mean 0.971128 6.914907 68.19556101 986.21083333 Source REP TRT DF Anova SS 8231.19321666 930320.32542500 Mean Square F Value 4115.59660833 0.88 310106.77514167 66.68 Pr > F 0.0005 0.0001 Duncan's Multiple Range Test for variable: KHỐI LƯỢNG RỄ CỎ VETIVER GIAI ĐOẠN 210 NGÀY (GAM/M2) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 4650.635 Number of Means Critical Range 136.2 141.2 143.7 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 1393.40 B 1027.78 B B 906.54 3 C 617.12 The SAS System Analysis of Variance Procedure Dependent Variable: HÀM LƯỢNG As TRONG THÂN LÁ CỎ VETVER (150 NGÀY) Source DF Sum of Squares Mean Square F Value Pr > F Model 19.07543333 3.81508667 1.54 0.0045 Error 14.84073333 2.47345556 Corrected Total 11 33.91616667 R-Square C.V Root MSE YIELD7 Mean 0.562429 16.21085 1.57272234 9.70166667 Source DF Anova SS Mean Square F Value Pr > F REP 6.10406667 3.05203333 1.23 0.0157 TRT 12.97136667 4.32378889 1.75 0.2565 Duncan's Multiple Range Test for variable: HÀM LƯỢNG As TRONG THÂN LÁ CỎ VETVER (150 NGÀY) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 2.473456 Number of Means Critical Range 3.142 3.257 3.313 175 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 11.363 3 A A 9.747 A A 9.077 A A 8.620 The SAS System Analysis of Variance Procedure Dependent Variable: HÀM LƯỢNG As TRONG THÂN LÁ CỎ VETVER (210 NGÀY) Source DF Sum of Squares Mean Square F Value Pr > F Model 16.31761667 3.26352333 1.72 0.0229 Error 11.37168333 1.89528056 Corrected Total 11 27.68930000 R-Square C.V Root MSE YIELD8 Mean 0.589311 9.120185 1.37669189 15.09500000 Source DF Anova SS Mean Square F Value Pr > F REP 7.93625000 3.96812500 2.09 0.0243 TRT 8.38136667 2.79378889 1.47 0.3132 Duncan's Multiple Range Test for variable: HÀM LƯỢNG As TRONG THÂN LÁ CỎ VETVER (210 NGÀY) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 1.895281 Number of Means Critical Range 2.750 2.851 2.900 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 16.393 3 A A 15.250 A A 14.487 A A 14.250 The SAS System Analysis of Variance Procedure Dependent Variable: HÀM LƯỢNG As TRONG RỄ CỎ VETVER (150 NGÀY) Source DF Sum of Squares Mean Square F Value Pr > F Model 408.03111667 81.60622333 2.74 0.0263 Error 178.96885000 29.82814167 Corrected Total 11 586.99996667 R-Square C.V Root MSE YIELD10 Mean 0.695113 8.082950 5.46151459 67.56833333 Source DF Anova SS Mean Square F Value Pr > F REP 354.69501667 177.34750833 5.95 0.0377 TRT 53.33610000 17.77870000 0.60 0.6404 176 Duncan's Multiple Range Test for variable: HÀM LƯỢNG As TRONG RỄ CỎ VETVER (150 NGÀY) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 29.82814 Number of Means Critical Range 10.91 11.31 11.51 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 69.967 A A 68.967 3 A A 66.867 A A 64.473 The SAS System Analysis of Variance Procedure Dependent Variable: HÀM LƯỢNG As TRONG RỄ CỎ VETVER (210 NGÀY) Source DF Sum of Squares Mean Square F Value Pr > F Model 226.50631667 45.30126333 1.21 0.4036 Error 223.77285000 37.29547500 Corrected Total 11 450.27916667 R-Square C.V Root MSE YIELD4 Mean 0.503035 5.643223 6.10700213 108.21833333 Source DF Anova SS Mean Square F Value Pr > F REP 140.77921667 70.38960833 1.89 0.2313 TRT 85.72710000 28.57570000 0.77 0.5532 Duncan's Multiple Range Test for variable: HÀM LƯỢNG As TRONG RỄ CỎ VETVER (210 NGÀY) NOTE: This test controls the type I comparisonwise error rate, not the experimentwise error rate Alpha= 0.05 df= MSE= 37.29548 Number of Means Critical Range 12.20 12.65 12.87 Means with the same letter are not significantly different Duncan Grouping Mean N TRT A 111.677 A A 109.927 A A 106.053 A A 105.217 3