PHẦN MỞ ĐẦU Sự cần thiết luận án Hiện nước ta, keo (Acacia) loài chủ lực trồng rừng công nghiệp cung cấp nguyên liệu cho ngành chế biến gỗ nước xuất Sự phát triển rừng trồng công nghiệp có đóng góp đáng kể phát triển kinh tế xã hội, góp phần giảm nhập gỗ đáp ứng nhu cầu chế biến hàng năm tăng cao Tuy nhiên, khuynh hướng suy giảm suất rừng qua chu kỳ (CK) kinh doanh mối lo ngại doanh nghiệp người trồng rừng không nước mà nhiều quốc gia giới, mà nguyên nhân quan trọng quản lý lập địa thiếu bền vững trồng rừng Kết nghiên cứu mạng lưới dự án Trung tâm Nghiên cứu Lâm nghiệp Quốc tế (CIFOR) thực 16 nước vùng Nhiệt đới Á nhiệt đới cho thấy: việc quản lý hợp lý vật liệu hữu sau khai thác, kiểm soát thảm thực bì sử dụng phân bón phù hợp có tác dụng tích cực đến độ phì đất suất rừng trồng qua chu kỳ kinh doanh (Nambiar, 1996) Một nghiên cứu bản, định lượng, có sở khoa học để xác định nguyên nhân suy giảm suất rừng trồng từ khía cạnh lập địa Việt Nam tản mạn, nghiên cứu sử dụng vật liệu hữu sau khai thác rừng để trả lại chất hữu cho đất chu kỳ sau mẻ chưa nghiên cứu sâu rộng, biện pháp canh tác truyền thống phát, đốt, dọn thực bì cày xới để trồng rừng Do vậy, tác giả thực luận án “Nghiên cứu ảnh hưởng vật liệu hữu để lại sau khai thác đến độ phì đất suất rừng trồng Keo tràm chu kỳ sau Phú Bình, tỉnh Bình Dương làm sở đánh giá toàn diện diễn biến độ phì đất suất rừng qua chu kỳ kinh doanh Mục tiêu nghiên cứu luận án Đánh giá mức độ ảnh hưởng vật liệu hữu để lại sau khai thác (VLHCSKT) đến độ phì đất lượng hóa mối quan hệ giữ vật liệu hữu để lại sau khai thác đến sinh trưởng, suất, sinh khối dinh dưỡng rừng trồng Keo tràm qua chu kỳ kinh doanh Xác định nguồn cung cấp, khả tích lũy cân dinh dưỡng thông qua để lại VLHCSKT làm sở cho đề xuất kỹ thuật quản lý lập địa trồng rừng trồng Keo tràm Ý nghĩa khoa học thực tiễn luận án * Ý nghĩa khoa học: góp phần làm sáng tỏ sở khoa học kinh doanh rừng trồng bền vững Việt Nam Đó tìm mối quan hệ biện chứng dinh dưỡng đất trồng với suất sinh học thực vật loài Keo tràm Từ tác động kỹ thuật đề xuất nghiên cứu trì vật liệu hữu sau khai thác rừng, bổ sung dinh dưỡng phân Lân phát huy có hiệu việc cải thiện độ phì đất suất rừng trồng Keo tràm * Ý nghĩa thực tiễn: Góp phần làm thay đổi nhận thức người trồng rừng, dần xóa bỏ tập quán canh tác truyền thống thiếu bền vững, là: lạm dụng giới làm đất, phát dọn, đốt, lấy VLHCSKT trước trồng rừng, cày xới hàng năm chăm sóc phòng chống cháy rừng Những đóng góp luận án - Đã xác định vai trò quản lý VLHCSKT đến cải thiện dinh dưỡng đất, nâng cao suất rừng trồng khả tự cân dinh dưỡng rừng trồng Keo tràm qua chu kỳ kinh doanh - Bước đầu nghiên cứu số mắt xích quan trọng chu trình dinh dưỡng rừng trồng Keo tràm, làm sở cho việc quản lý kinh doanh rừng trồng bền vững Việt Nam Đối tượng địa điểm nghiên cứu 5.1 Đối tượng nghiên cứu Đối tượng nghiên cứu đất tán rừng rừng trồng Keo tràm loài đánh giá qua chu kỳ kinh doanh, cụ thể sau: +) Chu kỳ 1: Rừng trồng Keo tràm trồng từ hạt, mật độ trồng 833 cây/ha, rừng trồng năm 1995 khai thác năm 2002 (7 năm) +) Chu kỳ 2: Rừng trồng thí nghiệm loài Keo tràm trồng từ hạt với nguồn giống tuyển chọn, mật độ trồng 1.667 cây/ha, rừng trồng năm 2002 khai thác năm 2008 (6 năm) +) Chu kỳ 3: Rừng trồng thí nghiệm loài Keo tràm dòng AA1 AA9 giống tiến kỹ thuật công nhận, mật độ trồng 1.667 cây/ha, rừng trồng năm 2008 theo dõi đến năm 2013 (5 năm) Vật liệu hữu sau khai thác: Trong nghiên cứu hiểu toàn cành, nhánh rừng trồng có đường kính < cm, tất bụi, thảm tươi vật rơi rụng tán rừng sau chặt hạ cắt ngắn có chiều dài từ 0,5 - 1m, rải toàn diện tích ô thí nghiệm, không đốt, không cày xới rừng chăm sóc năm đầu thông qua kiểm soát cỏ dại thuốc diệt cỏ Round-up với liều lượng lít/ha 5.2 Địa điểm nghiên cứu Nghiên cứu thực trạm Thực nghiệm lâm nghiệp Phú Bình, thuộc Viện Khoa học Lâm nghiệp Nam Bộ, nằm địa bàn xã Tam Lập, huyện Phú Giáo, tỉnh Bình Dương có vị trí tọa độ địa lý: 10 52’ 12” đến 110 30’ vĩ độ Bắc 1060 20’ đến 1070 06’ kinh độ Giới hạn phạm vi nghiên cứu + Chỉ đánh giá tính chất đất tầng từ - 10cm 10 - 20cm tầng đất dễ bị thay đổi tác động môi trường biện pháp kỹ thuật lâm sinh tác động + Chỉ nghiên cứu loại đất xám phù sa cổ trạm Thực nghiệm Lâm nghiệp Phú Bình, huyện Phú Giáo, tỉnh Bình Dương mà chưa có điều kiện nghiên cứu nhóm đất khác nhiều vùng sinh thái khác + Chỉ nghiên cứu số mắt xích quan trọng chu trình dinh dưỡng rừng trồng Keo tràm mà không tham vọng thiết lập chu trình dinh dưỡng Bố cục luận án Luận án gồm 129 trang, có 35 bảng 46 hình Kết cấu luận án gồm phần sau đây: Phần mở đầu Chương I: Tổng quan vấn đề nghiên cứu Chương II: Nội dung phương pháp nghiên cứu Chương III: Điều kiện tự nhiên khu vực nghiên cứu Chương IV: Kết nghiên cứu thảo luận Kết luận, tồn kiến nghị Tài liệu tham khảo: 89 tài liệu (40 tài liệu tiếng việt 49 tài liệu tiếng anh) Phần phụ lục gồm 22 trang với 35 bảng số liệu loại Chương 1: TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU Hiện diện tích rừng trồng Keo nước khoảng 1,1 triệu ha, rừng trồng Keo tràm khoảng gần 400.000 có xu hướng ngày tăng Bên cạnh lợi ích kinh tế, kinh doanh rừng trồng thiếu bền vững suất rừng giảm dần qua chu kỳ kinh doanh Phương thức canh tác truyền thống cày, đốt hay bỏ vật liệu sau khai thác cho nguyên nhân gây ảnh hưởng xấu đến độ phì làm xói mòn đất Nghiên cứu biện pháp quản lý lập địa thực nước có Lâm nghiệp tiên tiến từ 30 năm trước tiếp tục nghiên cứu ngày Các kết nghiên cứu cho thấy biện pháp quản lý lập địa việc để lại vật liệu hữu sau khai thác làm tăng độ phì đất suất rừng trồng số loài mọc nhanh như: Thông, Bạch đàn Keo tai tượng, … Tuy nhiên, Việt Nam nghiên cứu mẻ kết đạt khiêm tốn Tổng quan nghiên cứu nước giới biện pháp kỹ thuật bảo vệ nâng cao độ phì đất tập trung vào nội dung như: (i) Nghiên cứu quản lý vật liệu hữu sau khai thác rừng; (ii) Nghiên cứu quản lý dinh dưỡng đất qua bón phân cho rừng trồng (iii) Nghiên cứu quản lý thực vật qua kiểm soát thực vật cạnh tranh tán rừng Nghiên cứu quản lý VLHCSKT phần nghiên cứu quản lý lập địa Đối với loài Keo tràm có nghiên cứu nước đề cập đến, loại trừ dự án CIFOR thực từ năm 2002 – 2007 Bình Dương Do vậy, việc tiếp tục phát triển nghiên cứu chu kỳ cần thiết khẳng định vai trò quản lý VLHCSKT đánh giá, so sánh với chu kỳ kinh doanh rừng trước Hơn nữa, nghiên cứu ảnh hưởng để lại VLHCSKT rừng trồng Việt Nam nói chung rừng trồng Keo tràm nói riêng chưa nghiên cứu toàn diện độ phì đất, sinh trưởng, suất sinh khối rừng liên tục qua nhiều chu kỳ kinh doanh Chương 2: NỘI DUNG VÀ PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Nội dung nghiên cứu - Nghiên cứu ảnh hưởng vật liệu hữu để lại đến độ phì đất qua chu kỳ kinh doanh - Nghiên cứu ảnh hưởng vật liệu hữu để lại sau khai thác đến sinh trưởng, trữ lượng, suất sinh khối rừng qua chu kỳ kinh doanh - Nghiên cứu tích lũy chất dinh dưỡng khả cân dinh dưỡng rừng trồng Keo tràm chu kỳ 2.2 Phương pháp nghiên cứu 2.2.1 Phương pháp bố trí thí nghiệm • Kế thừa số liệu đất, sinh trưởng, suất sinh khối rừng trồng thí nghiệm chu kỳ & từ kết nghiên cứu dự án “Quản lý lập địa suất rừng trồng Keo Việt Nam” thực từ 2002 – 2007 tổ chức CIFOR tài trợ • Phương pháp bố trí thí nghiệm chu kỳ Thí nghiệm bố trí rừng trồng chu kỳ theo khối ngẫu nhiên đầy đủ gồm công thức lần lặp lại Phương pháp bố trí thí nghiệm dựa theo phương pháp nghiên cứu quản lý lập địa rừng trồng Nhiệt đới Trung tâm Nghiên cứu Lâm nghiệp Quốc tế (CIFOR), (Tiarks cs., 1998) sau: Fl: Lấy toàn vật liệu hữu sau khai thác chu kỳ Fm: Để lại vật liệu hữu sau khai thác chu kỳ chu kỳ Fh: Để lại gấp đôi vật liệu hữu sau khai thác chu kỳ 1; để lại VLHCSKT rừng chu kỳ đồng thời bón lót 300g super phosphate 16,5% P205/cây tương đương 36 kgP/ha Tổng diện tích thí nghiệm ha, gồm phần: Ô thí nghiệm gồm công thức thí nghiệm bố trí lần lặp lại Tổng diện tích 15 ô thí nghiệm 17.280 m2, diện tích ô 1.152 m (12 hàng x 16 cây/hàng), diện tích ô đo đếm 576 m (8 hàng x 12 cây/hàng = 96 cây), diện tích vùng đệm 576 m2 gồm hàng bao xung quanh để hạn chế tác động yếu tố bên làm ảnh hưởng tới kết nghiên cứu luận án Ngoài ra, bố trí ô phụ thiết kế để lấy mẫu sinh khối hàng năm với diện tích 720 m2 • Phương pháp phân tích xử lý số liệu - Phân tích tính chất vật lý, hóa học đất, mẫu VLHCSKT, vật rụng, phân hủy thực vật, sinh khối thành phần chất dinh dưỡng thực vật thực phòng phân tích môn sinh thái môi trường rừng thuộc Viện Khoa học Lâm nghiệp Nam Bộ - Phân tích vi sinh vật (VSV) tổng số VSV phân giải lâm thực phòng thí nghiệm vi sinh Viện Khoa học Nông nghiệp Miền Nam - Số liệu thu thập xử lý theo phương pháp thống kê sinh học với phần mềm Genstat version 12.0 Excel 7.0 để tính toán Sử dụng trắc nghiệm Duncan để xếp hạng tiêu đánh giá theo mức độ khác biệt có ý nghĩa mặt thống kê 2.2.2 Phương pháp thu thập số liệu phân tích theo nội dung nghiên cứu • Phương pháp nghiên cứu độ phì đất qua chu kỳ kinh doanh - Lấy mẫu đất định kỳ vào tháng hàng năm để phân tích tính chất vật lý hóa học: Sử dụng khoan tay, ô thí nghiệm lấy từ 10 điểm dải ô với độ sâu tầng đất từ - 10cm 10 - 20cm Những mẫu độ sâu tầng đất trộn lại thành mẫu hỗn hợp Từ mẫu hỗn hợp, lấy hai mẫu phụ ngẫu nhiên phơi khô không khí Một nửa mẫu phụ dùng cho phân tích nửa lại lưu trữ Xác định dung trọng đất ống kim loại hình trụ đường kính 10cm, thể tích 1dm lấy mẫu phần đất không bị tác động từ điểm ô thí nghiệm tầng phát sinh – 10cm từ 10 - 20cm - Mẫu đánh giá động vật đất lấy từ công thức x lần lặp x lần/năm = 18 mẫu/năm, diện tích ô mẫu x 1m sâu 0,1m Đất lấy toàn làm nhỏ để đếm toàn động vật đất có kích thước > 0,2cm Mẫu phân tích VSV tổng số VSV phân giải lân lấy từ công thức x lần lặp x vị trí lấy khác x lần/năm = 54 mẫu/năm, diện tích ô mẫu 10 x 10 x 10cm Vị trí lấy mẫu cách gốc - 30cm, 30 - 80cm 80 150cm Mẫu sau lấy bảo tươi thùng đá đem phân tích phòng thí nghiệm Thời gian lấy mẫu định kỳ lần/năm vào mùa khô (tháng 3) mùa mưa (tháng hàng năm) - Phương pháp phân tích mẫu đất: thực theo qui trình phòng thí nghiệm (sổ tay phân tích đất, nước, phân bón trồng Viện Nông hóa thổ nhưỡng, 1998) có tham khảo tài liệu từ trung tâm thông tin đất Quốc tế ISRIS/FAO” tác giả van Reeuwijk,1995 - Phương pháp phân tích vi sinh vật tổng số (CFU/g) áp dụng tiêu chuẩn TCVN 4881-1989 TCVN 4884:2001 “Vi sinh vật học – hướng dẫn chung định lượng vi sinh vật – kỹ thuật đếm khuẩn lạc 300C” - Phương pháp phân tích vi sinh vật phân giải lân (CFU/g) áp dụng TCVN 6167:1996 “Phân bón vi sinh vật phân giải hợp chất phốt khó tan Phosphat-solubling microbial fertilizer” - Cách tính thành phần dinh dưỡng đất: áp dụng công thức X (kg) = BD (kg/m3) x 10.000 (m2) x SD (m) x TPDD (mg/kg) Trong đó: X thành phần dinh dưỡng đất cần tìm; BD dung trọng theo tầng đất; 10.000m2 diện tích cho 1ha; SD: độ sâu tầng đất TPDD thành phần dinh dưỡng chất X đất • Phương pháp nghiên cứu sinh trưởng, trữ lượng, suất sinh khối mặt đất rừng qua chu kỳ kinh doanh - Đánh giá sinh trưởng rừng: + Số liệu sinh trưởng, trữ lượng, suất, sinh khối rừng chu kỳ 1, đề tài luận án kế thừa số liệu dự án CIFOR (2002 - 2007) làm sở để so sánh đánh giá qua chu kỳ kinh doanh + Đo rừng trồng thí nghiệm cuối chu kỳ trước khai thác, đo sinh trưởng rừng chu kỳ định kỳ 1lần/năm, đo định vị toàn số công thức thí nghiệm theo lần lặp, tiêu gồm: Tỷ lệ sống (%), Hvn (m), D1.3 (cm) - Xác định suất, trữ lượng rừng sinh khối rừng: Ở CK2 trước khai thác rừng trồng CK3: Đo tiêu D 1.3, Hvn công thức thí nghiệm giải tích 15 tiêu chuẩn đại diện cho cấp kính khác để xác định trữ lượng lâm phần Trong 15 mẫu chọn để lấy mẫu theo phận như: thân, vỏ, cành, hoa phân tích thành phần dinh dưỡng - Tính toán sinh khối: Từ 15 tiêu chuẩn giải tích hàng năm tiến hành lập phương trình tương quan đường kính trọng lượng khô cá thể có dạng hàm Mũ: Y= a Xb • Phương pháp nghiên cứu tích lũy chất dinh dưỡng khả cân dinh dưỡng rừng trồng chu kỳ Xác định tốc độ phân huỷ vật liệu hữu sau khai thác Vật liệu lấy từ rừng khai thác chu kỳ gồm: Lá; Cành < cm; Cành từ - 5cm cắt thành đoạn dài 1cm Khối lượng mẫu khô 25g cho vào túi lưới nylon kích thước 20 x 20cm có mắt lưới rộng 2mm Mẫu đặt rải ô thí nghiệm rừng trồng chu kỳ 3, với tổng số 210 mẫu (5 lặp x 14 lần x loại vật liệu) Thời gian khảo sát 14 lần vào thời điểm tháng thứ: 0, ½, 1, 2, 3, 4, 5, 7, 9, 12, 15, 18, 21, 24 Mỗi lần thu mẫu x loại vật liệu = 15 mẫu, rửa để lấy phần vật liệu chưa phân hủy mang sấy khô nhiệt độ 760C trọng lượng không đổi cân để xác định khối lượng khô đánh giá lượng trình phân hủy - Xác định nguồn dinh dưỡng có khả bổ sung cho đất sau: + Xác định lượng VLHCSKT việc tính trọng lượng sinh khối toàn cành, có đường kính < 5cm cây, hoa, quả,… CK2 lấy mẫu phân tích thành phần dinh dưỡng loại VLHC để lại sau khai thác + Xác định sinh khối bụi, thảm tươi thảm mục tán rừng việc bố trí đặt 25 lưới ô vuông kích thước 1m (1 x 1m) rải ô thí nghiệm để phân loại thu gom tất thành phần : bụi, cỏ, cành, lá, (tươi khô) Mẫu sau thu sấy khô nhiệt độ 76 0C 48 đến trọng lượng không đổi để tính sinh khối khô phân tích hàm lượng chất dinh dưỡng + Xác định lượng vật rụng chu kỳ thu mẫu vật rụng hàng năm bẫy lượng rơi có kích thước 1m2 (1x1m) công thức Fm, gồm bẫy x lặp = 25 Thời gian thu mẫu từ bắt đầu thu tuổi đến tuổi (2 lần/tháng vào ngày 15 30) Xác định trọng lượng khô sấy mẫu 76 0C 48 trọng lượng không đổi Hàng năm, mẫu trộn lấy để phân tích thành chất dinh dưỡng gồm: N, P, K, Ca Mg - Xác định nguồn dinh dưỡng tích lũy + Lấy mẫu phân tích: Để xác định nguồn dinh dưỡng tích lũy rừng thí nghiệm chu kỳ thông qua giải tích 15 tiêu chuẩn/năm rải theo cấp kính để xác định sinh khối tươi, khô phân tích hàm lượng chất dinh dưỡng tích lũy + Phương pháp phân tích mẫu thực vật: áp dụng theo phân tích đất Quốc tế tác giả L.P van Reeuwijk, 1995 (ISRIC), cụ thể sau: N phương pháp Kjeldahl; P phương pháp so màu; K phương pháp quang kế lửa; Ca Mg phương pháp hấp thụ nguyên tử + Tính toán số liệu: Từ tổng sinh khối theo phận thân công thức thí nghiệm tính toán cho thành phần chất dinh dưỡng qui đổi 01 theo công thức: X (kg/ha) = ∑[SKBPi (tấn/ha) x NC (kg/ha] Trong đó: X thành phần dinh dưỡng cần tìm; SKBPi tổng sinh khối phận cây; Bpi phận thân, cành lớn, cành nhỏ, vỏ cây; NC hàm lượng dinh dưỡng X theo phận Chương 3: ĐIỀU KIỆN TỰ NHIÊN – KINH TẾ XÃ HỘI KHU VỰC NGHIÊN CỨU 3.1 Điều kiện tự nhiên Khu vực nghiên cứu trạm thực nghiệm lâm nghiệp Phú Bình, thuộc Viện Khoa học Lâm nghiệp Nam Bộ, phần đất công ty TNHH MTV Lâm nghiệp Bình Dương quản lý thuộc huyện Phú Giáo, tỉnh Bình Dương Địa hình dạng đồi núi thấp, độ dốc trung bình từ - 0, độ cao trung bình 20 - 25m so với mực nước biển Khí hậu đặc trưng vùng miền Đông Nam Bộ nhiệt đới gió mùa, có mùa mưa mùa khô rõ rệt, lượng mưa trung bình năm biến động từ biến động từ 1.671,3 - 2.319,7 mm Đất thuộc loại đất xám vàng phù sa cổ (Acrisol chromic), mẫu chất đá phiến Đất có cấu trúc thịt sét pha cát tầng mặt khoảng - 20cm, tầng đất dày, địa hình tương đối bằng, bề mặt đất bị xói mòn Phản ứng môi trường đất trung bình giao động pH = 4,0 không thay đổi nhiều tầng Độ chua trao đồi pHkcl biến động 3,98 - 4,0 cho thấy đất chua Hàm lượng cation trao đổi Ca+2 mức trung bình (0,38 - 0,40 cmol/kg), K+ Mg+2 Hàm lượng dinh dưỡng nghèo, chất hữu Cacbon biến động khoảng 2%, N < 0,2%, lân dễ tiêu từ - mg/kg 3.2 Điều kiện kinh tế - xã hội Khu vực nghiên cứu luận án có lợi vị trí địa lý có giao thông lại thuận lợi, nằm khu vực tiếp giáp khu đô thị khu công nghiệp phát triển với khu vực nông nghiệp nông thôn Mật độ dân cư trung bình 151 người/km2, nguồn lao động dồi dào, có kinh nghiệm tổ chức sản xuất kinh doanh nông nghiệp Đời sống hộ gia đình chủ yếu trồng công nghiệp (cao su, điều) phát triển chăn nuôi gia súc, gia cầm Tuy nhiên, thu nhập bình quân đầu người mức trung bình đến so với toàn tỉnh nên áp lực vào tài nguyên rừng không lớn, công tác quản lý bảo vệ rừng thuận lợi để lại vật liệu hữu sau khai thác rừng Chương 4: KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 4.1.1 Ảnh hưởng để lại VLHCSKT đến độ phì đất 4.1.2 Tính chất vật lý đất a) Chỉ tiêu dung trọng đất Dung trọng tiêu quan trọng để đánh giá độ tơi xốp đất Kết đánh giá diễn biến dung trọng đất tầng đất từ – 10cm 10 – 20cm qua chu kỳ kinh doanh thể hình 4.1 4.2 đây: Hình 4.1; 4.2 Biến động dung trọng đất tầng đất 0-10cm 10-20cm Như vậy, dung trọng đất cải thiện qua chu kỳ kinh doanh rừng trồng Keo tràm tầng đất Giữa công thức có khác biệt rõ rệt, cụ thể đầu chu kỳ dung trọng đất 1,41 g/cm đến giai đoạn tuổi rừng trồng chu kỳ (sau 11 năm) công thức Fh có dung trọng nhỏ (1,23 g/cm3), công thức Fm (1,28 g/cm 3) cao công thức Fl (1,32 g/cm3) Điều chứng tỏ độ xốp đất cải thiện để lại vật liệu hữu sau khai thác qua chu kỳ kinh doanh b) Chỉ tiêu thành phần giới đất Thành phần giới đất tỷ lệ phần trăm cấp hạt có kích thước khác đoàn lạp đất trạng thái bị phá hủy Thành phần giới đặc trưng cho kết cấu đất biến đổi thời gian ngắn Kết nghiên cứu luận án đánh giá công thức Fh có cường độ tác động mạnh thông qua việc để lại VLHCSKT chu kỳ sau, kết tổng hợp bảng 4.1 đây: Bảng 4.1 Chỉ tiêu thành phần giới đất qua chu kỳ kinh doanh TPCG CK1 (2002) CK2 (2008) CK3 (2013) Tầng đất từ - 10cm (tỉ lệ % cấp hạt) Tầng đất từ 10 - 20cm (tỉ lệ % cấp hạt) Sét (