Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Vùng trung du và miền núi Việt Nam chiếm khoảng ba phần tư diện tích tự nhiên của cả nước, giữ vị trí chiến lược đối với sự phát triển kinh tế - xã hội, nông lâm nghiệp và môi trường sinh thái. Tuy nhiên, nguồn tài nguyên đất đồi dốc tại khu vực này đang phải đối mặt với tình trạng thoái hóa nghiêm trọng. Đến cuối thập niên 1980, diện tích đất trống đồi núi trọc tại Việt Nam đã tăng nhanh chóng (từ 13,8 triệu ha năm 1985 lên 14,6 triệu ha năm 1988), trong đó có hơn 10 triệu ha đất đồi núi đã và đang bị xói mòn, kiệt quệ độ phì nhiêu.

Nguyên nhân cốt lõi dẫn đến tình trạng trên bắt nguồn từ các tập quán canh tác lạc hậu và việc sử dụng đất bất hợp lý: nạn phá rừng làm nương rẫy, độc canh các loại cây lương thực lấy củ (đặc biệt là cây sắn - Manihot esculenta) và cây màu ngắn ngày trên các sườn dốc mà không có biện pháp che phủ, bón phân hay công trình chống xói mòn. Thực tiễn sản xuất cho thấy, sau một vài vụ canh tác sắn hoặc ngô theo phương thức phát đốt chọc lỗ truyền thống, đất đai bị rửa trôi lớp đất mặt màu mỡ, hàm lượng chất hữu cơ và các chất dinh dưỡng bị suy kiệt, tầng đất mặt bị trơ sỏi đá, tạo thành các loại đất thoái hóa nghiêm trọng.

Mặc dù việc nghiên cứu xói mòn và cải tạo đất đồi đã được nhiều nhà khoa học quan tâm, song phần lớn các công trình trước đây chỉ dừng lại ở các nghiên cứu định tính, mô tả hiện trạng hoặc quan sát riêng rẽ từng biện pháp kỹ thuật đơn lẻ. Luận án giải quyết khoảng trống nghiên cứu bằng cách đánh giá toàn diện, đồng bộ tác động của các hệ thống canh tác khác nhau (độc canh, xen canh, luân canh, nông lâm kết hợp, trồng cây phân xanh, cây họ đậu, cây công nghiệp lâu năm) đến tập hợp các chỉ tiêu lý, hóa, sinh học đất và mức độ xói mòn rửa trôi, đồng thời áp dụng các công cụ toán học đa biến để lượng hóa và xác định phương thức canh tác tối ưu cho vùng đất đồi trung du.

Mục tiêu nghiên cứu

  1. Điều tra, đánh giá hiện trạng độ phì nhiêu đất feralit và diễn biến thoái hóa đất dưới tác động của các phương thức canh tác điển hình tại vùng đồi trung du huyện Thanh Ba, tỉnh Vĩnh Phú.
  2. Xác định định lượng ảnh hưởng của các loại hình cây trồng (cây lâm nghiệp, cây công nghiệp, cây lương thực, cây cải tạo đất/phân xanh) và các biện pháp kỹ thuật canh tác (trồng thuần loài, xen canh, gối vụ, che phủ, bón phân) đến tính chất vật lý, hóa học, sinh học và lượng đất mất do xói mòn.
  3. Ứng dụng mô hình phân loại số học (Numerical Taxonomy) và phương pháp toán học đa biến nhằm phân nhóm, lựa chọn các phương thức canh tác hợp lý nhất, vừa bảo vệ và nâng cao độ phì nhiêu đất, vừa đảm bảo năng suất cây trồng và hiệu quả kinh tế lâu bền.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Đất feralit đỏ vàng phát triển trên phiến thạch sét, đá phiến và phù sa cổ; các phương thức canh tác nông - lâm nghiệp (sắn thuần loài, chè, chè có cây che bóng, rừng trồng keo lá tràm, bạch đàn, bồ đề, các mô hình xen canh sắn với cốt khí, lạc, đậu tương, đậu đen, muồng đen).
  • Phạm vi không gian: Địa bàn huyện Thanh Ba (Thanh Hòa), tỉnh Vĩnh Phú (nay thuộc tỉnh Phú Thọ), vùng trung du điển hình của miền Bắc Việt Nam; kết hợp số liệu quan trắc từ Trạm Khí tượng Thủy văn Phú Hộ và Trạm Nghiên cứu Nông lâm nghiệp Khôi Vĩ.
  • Phạm vi thời gian: Số liệu thực nghiệm đồng ruộng, ô tiêu chuẩn và phân tích phòng thí nghiệm được tiến hành liên tục trong giai đoạn 1985 - 1988, kết hợp chuỗi số liệu khí tượng - thổ nhưỡng giai đoạn 1979 - 1988.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trên thế giới

Nghiên cứu về xói mòn đất và bảo vệ độ phì nhiêu đất dốc trên thế giới đã được triển khai từ đầu thế kỷ XX:

  • Bennett (1928, 1939): Đặt nền móng cho khoa học bảo tồn đất tại Mỹ, xác lập các định luật về dòng chảy mặt, mối tương quan giữa độ dốc, lượng mưa và lượng đất mất, đồng thời khởi xướng các kỹ thuật canh tác theo đường đồng mức và dải cây bảo vệ.
  • Stallings (1957): Nghiên cứu cơ chế tác động của động năng hạt mưa đến sự phá vỡ cấu trúc hạt đất mặt và vai trò của thảm phủ thực vật trong việc hấp thụ xung lực giọt mưa.
  • Sokal và Sneath (1963): Xây dựng lý thuyết Phân loại số học (Numerical Taxonomy), cung cấp khung phương pháp luận sử dụng toán học đa biến và ma trận khoảng cách hình học để phân nhóm các đối tượng phức tạp dựa trên tập hợp nhiều thuộc tính định lượng.
  • Lance và Williams (1967): Phát triển các thuật toán phân cụm phân cấp (Hierarchical Clustering), chuẩn hóa công thức tính toán khoảng cách liên nhóm phục vụ phân loại tự động trên máy tính điện tử.
  • Các công trình quốc tế tại vùng nhiệt đới châu Á, châu Phi của Dudal (1960), Nye và Greenland (1960), Sanchez (1976) đã chứng minh vai trò quyết định của cây họ đậu và hệ thống nông lâm kết hợp trong việc duy trì chu trình dinh dưỡng và hạn chế chua hóa đất.

Các hướng nghiên cứu tại Việt Nam

  • Đỗ Đình Sâm, Vũ Cao Thái, Vũ Tự Lập, Cao Liêm, Trần Khải: Nghiên cứu đặc điểm phát sinh, tính chất lý hóa học của đất feralit đồi núi Việt Nam, phân tích quá trình rửa trôi sesquioxit ($Fe_2O_3, Al_2O_3$), cation kiềm ($Ca^{2+}, Mg^{2+}, K^+$) và sự tích lũy nhôm di động ($Al^{3+}$) gây độc cho cây trồng khi thảm thực vật bị phá hủy.
  • Lê Duy Thước, Đào Châu Thu (1972 - 1974): Thí nghiệm ứng dụng cây phân xanh (cốt khí, muồng) xen trong nương chè và đồi cây ăn quả tại Vĩnh Phú, ghi nhận hiệu quả tăng năng suất chè búp và cải thiện hàm lượng đạm, lân hữu dụng trong đất.
  • Nguyễn Quang Mỹ (1978, 1981): Khảo sát định lượng lượng đất xói mòn dưới các thảm phủ thực vật tại vùng Tam Đảo và Thanh Ba (Vĩnh Phú), chỉ ra rằng lượng đất trôi ráo ở đất trống đồi trọc có thể lên tới 97,8 tấn/ha/năm, trong khi dưới tán rừng rậm chỉ dao động từ 0,35 - 1,2 tấn/ha/năm.

Khoảng trống luận án lựa chọn giải quyết

Các công trình trước đây tại Việt Nam thường đánh giá riêng lẻ từng khía cạnh (hoặc chuyên sâu về xói mòn vật lý, hoặc khảo sát nông học về năng suất cây trồng) trên quy mô diện hẹp. Luận án lấp đầy khoảng trống này bằng cách:

  1. Kết hợp đồng thời quan trắc vi khí hậu, đo đạc dòng chảy mặt, lượng đất xói mòn cơ học với việc theo dõi biến động động thái của 15 chỉ tiêu lý - hóa - sinh học đất qua chu kỳ canh tác 3 - 4 năm.
  2. Lần đầu tiên ứng dụng hệ thống thuật toán phân loại số học của Sokal - Sneath và Lance - Williams để xử lý tập số liệu thổ nhưỡng đa biến, xây dựng tiêu chuẩn định lượng khách quan phân biệt các phương thức canh tác suy thoái và các phương thức canh tác bảo vệ, cải tạo độ phì nhiêu đất đồi.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và khái niệm

  1. Lý thuyết phát sinh và thoái hóa đất feralit nhiệt đới ẩm: Trong điều kiện nhiệt độ cao (trung bình 23 - 24°C) và lượng mưa lớn tập trung (1600 - 1800 mm/năm), quá trình phong hóa hóa học diễn ra mạnh mẽ. Khi mất thảm phủ thực vật, quá trình khoáng hóa mùn diễn ra nhanh hơn mùn hóa; các cation kiềm và kiềm thổ bị rửa trôi triệt để, để lại phức hệ hấp phụ bão hòa ion $H^+$ và $Al^{3+}$, dẫn đến chua hóa tầng mặt ($pH_{KCl} < 4,5$) và gia tăng hàm lượng nhôm trao đổi gây nghẽn hấp thu dinh dưỡng của rễ cây.
  2. Khung lý thuyết sinh thái nông nghiệp và nông lâm kết hợp: Cấu trúc tầng tán đa tầng (cây gỗ lâu năm + cây bụi họ đậu + cây che phủ mặt đất) mô phỏng cấu trúc rừng tự nhiên, tạo lớp thảm mục bảo vệ mặt đất, làm giảm động năng giọt mưa, duy trì độ xốp và tuần hoàn chất hữu cơ liên tục thông qua khối lượng rơi rụng và nốt sần cố định đạm sinh học.
  3. Mô hình phân loại số học (Numerical Taxonomy): Đánh giá mức độ tương đồng giữa các phương thức canh tác dựa trên khoảng cách hình học đa chiều $d_{jk}$ giữa các vector chỉ tiêu đất:

$$d_{jk} = \sqrt{\sum_{i=1}^{n} (x_{ij} - x_{ik})^2}$$

Trong đó $x_{ij}, x_{ik}$ là giá trị chuẩn hóa của chỉ tiêu thứ $i$ ở phương thức $j$ và $k$. Quá trình hợp nhóm được thực hiện tuần tự theo thuật toán Lance và Williams (1967):

$$d_{h(ij)} = \alpha_i d_{hi} + \alpha_j d_{hj} + \beta d_{ij} + \gamma |d_{hi} - d_{hj}|$$

                                  MA TRẬN KHOẢNG CÁCH D(jk)
                             Tìm cặp có khoảng cách nhỏ nhất: Min d(ij)
                             Gộp nhóm (i, j) thành cụm mới (ij)
                             Tính khoảng cách d(h, ij) theo Lance-Williams
                               Xây dựng Dendrogram phân cấp

Phương pháp nghiên cứu thực nghiệm

  • Bố trí thí nghiệm đồng ruộng:
    • Thiết lập các ô tiêu chuẩn đo xói mòn kích thước $100\text{ m}^2$ ($5\text{ m} \times 20\text{ m}$) trên các sườn đồi có độ dốc $10 - 15^\circ$, có rãnh bao quanh và bể thu gom bùn cát, nước chảy tràn ở cuối ô.
    • Bố trí các công thức thí nghiệm chậu vại/chậu cát và ô vi mô tại hiện trường nhằm theo dõi lượng hút thu dinh dưỡng và động thái khoáng hóa của thân lá cây họ đậu (cốt khí, muồng).
  • Phương pháp lấy mẫu và phân tích đất trong phòng thí nghiệm:
    • Lấy mẫu tầng đất mặt ($0 - 20\text{ cm}$) và tầng dưới ($20 - 40\text{ cm}$) định kỳ trước và sau mỗi vụ thu hoạch (1985 - 1988).
    • Độ chua: Đo $pH_{H_2O}$ và $pH_{KCl}$ bằng phương pháp điện cực thủy tinh tỷ lệ đất:dung dịch 1:2,5.
    • Chất hữu cơ (mùn): Xác định theo phương pháp Tiurin (oxy hóa bằng dung dịch $K_2Cr_2O_7 + H_2SO_4$).
    • Đạm tổng số và dễ tiêu: Xác định theo phương pháp Kjeldahl và đạm thủy phân theo Tyurin - Kononova.
    • Lân tổng số và dễ tiêu: Phân hủy bằng hỗn hợp $H_2SO_4 + HClO_4$, đo màu lân dễ tiêu theo phương pháp Bray 2 hoặc Oniani.
    • Kali tổng số và dễ tiêu: Xác định bằng phương pháp quang kế ngọn lửa (Flame photometer) sau khi trích ly bằng dung dịch $CH_3COONH_4\text{ 1N}$.
    • Nhôm di động ($Al^{3+}$): Xác định theo phương pháp Sokolov (trích ly bằng $KCl\text{ 1N}$, chuẩn độ bằng $NaOH\text{ 0,01N}$).
    • Chỉ tiêu vật lý: Dung trọng xác định bằng ống đai dung tích chuẩn; tỷ trọng bằng bình tỷ trọng; độ xốp tính theo công thức toán học; thành phần cơ giới phân tích theo phương pháp ống hút Robinson; khả năng thấm nước xác định bằng ống thấm hình trụ đồng tâm.
    • Chỉ tiêu sinh học: Đếm số lượng vi sinh vật tổng số (vi khuẩn, nấm mốc, xạ khuẩn) bằng phương pháp pha loãng cấy trên môi trường thạch đĩa petri.

Nội dung chính theo từng chương

Chương 1: Điều kiện tự nhiên, kinh tế - xã hội vùng đồi trung du huyện Thanh Ba, tỉnh Vĩnh Phú

  • Địa hình - Địa mạo: Huyện Thanh Ba mang đặc trưng của địa hình trung du đồi bát úp, độ cao tuyệt đối phổ biến dưới $100\text{ m}$, sườn dốc thoải $5 - 15^\circ$, xen kẽ các thung lũng hẹp và bậc thềm tích tụ phù sa cổ. Vùng đồi dạng bát úp chiếm trên 70% diện tích tự nhiên, tạo điều kiện thuận lợi cho sản xuất nông lâm kết hợp nhưng nguy cơ xói mòn rửa trôi rất cao khi mất thảm phủ.
  • Đặc điểm khí hậu: Chuỗi số liệu 10 năm (1979 - 1988) tại Trạm Phú Hộ phản ánh nền nhiệt ẩm cận nhiệt đới gió mùa:
    • Nhiệt độ trung bình năm đạt $23,4^\circ\text{C}$; tháng lạnh nhất là tháng 1 ($15,8^\circ\text{C}$), tháng nóng nhất là tháng 7 ($28,9^\circ\text{C}$).
    • Lượng mưa trung bình năm đạt $1683,6\text{ mm}$, phân bố không đều: mùa mưa từ tháng 5 đến tháng 10 chiếm $84,9%$ tổng lượng mưa cả năm (tháng cao nhất đạt trên $360\text{ mm}$). Mùa khô kéo dài từ tháng 11 đến tháng 4 năm sau với lượng bốc hơi vượt lượng mưa, gây hạn hán cục bộ.
    • Độ ẩm không khí trung bình $84%$. Tổng số giờ nắng trong năm đạt từ 1400 - 1600 giờ.
  • Hiện trạng sử dụng đất và tập quán canh tác:
    • Diện tích cây lương thực vùng đồi chủ yếu là sắn ($2800\text{ ha}$) với sản lượng cung ứng nguyên liệu tinh bột và thức ăn gia súc lớn.
    • Cây công nghiệp chủ lực là chè ($> 2500\text{ ha}$), cung cấp nguyên liệu cho các nhà máy chế biến chè xuất khẩu trong tỉnh.
    • Cây lâm nghiệp gồm bồ đề, bạch đàn, keo lá tràm, muồng, trẩu.
    • Tập quán canh tác nương rẫy, trồng sắn độc canh quảng canh không bón phân hóa học hay phân hữu cơ làm kiệt quệ dinh dưỡng đất chỉ sau 2 - 3 năm canh tác.
Chỉ tiêu khí hậu tại trạm Phú Hộ Giá trị đo đạc trung bình (1979 - 1988)
Nhiệt độ trung bình năm $23,4^\circ\text{C}$
Nhiệt độ tối cao tuyệt đối $39,2^\circ\text{C}$
Nhiệt độ tối thấp tuyệt đối $5,6^\circ\text{C}$
Lượng mưa trung bình năm $1683,6\text{ mm}$
Lượng bốc hơi trung bình năm $836,8\text{ mm}$
Độ ẩm không khí trung bình $84,0%$
Tỷ lệ lượng mưa mùa mưa (tháng 5 - 10) $84,9%$

Chương 2: Thực trạng xói mòn và suy thoái lý hóa tính đất đồi dưới các phương thức canh tác

  • Xói mòn và mất đất mặt:
    • Ở phương thức độc canh sắn trên đất dốc không có đai cây che phủ, lượng đất trôi ráo trung bình hàng năm lên tới $120 - 366\text{ tấn/ha/năm}$, tương ứng với việc mất đi lớp đất mặt dày từ $0,8 - 1,5\text{ cm}$ mỗi năm.
    • Khi trồng sắn xen lạc hoặc trồng cây phân xanh làm dải chắn (cốt khí), lượng đất mất giảm xuống còn $4,0 - 12,6\text{ tấn/ha/năm}$ (giảm $85 - 95%$ lượng đất trôi).
    • Ở nương chè có cây che bóng (muồng đen) kết hợp thảm cốt khí, lượng đất xói mòn chỉ còn $0,35 - 2,5\text{ tấn/ha/năm}$.
  • Sự suy thoái hóa tính đất:
    • Mùn tổng số: Canh tác sắn thuần loài làm hàm lượng mùn giảm mạnh từ $2,14%$ xuống $1,02 - 1,20%$ sau 3 năm. Ngược lại, các công thức xen cây họ đậu và bổ sung tàn dư hữu cơ giúp hàm lượng mùn tăng lên $2,8 - 3,2%$.
    • Độ chua và nhôm di động: Ở đất độc canh sắn thoái hóa, $pH_{KCl}$ hạ xuống $3,8 - 4,0$; hàm lượng $Al^{3+}$ di động tăng vọt lên $5,7 - 7,9\text{ meq/100g}$ đất (vượt ngưỡng gây hại $6\text{ meq/100g}$). Tại các ô canh tác có phủ thảm hữu cơ và bón lân nung chảy kết hợp vôi, $Al^{3+}$ giảm xuống còn $0,1 - 1,4\text{ meq/100g}$.
    • Dinh dưỡng N-P-K: Hàm lượng đạm tổng số ở đất canh tác sắn thuần chỉ còn $0,06 - 0,08%$, lân dễ tiêu ($P_2O_5$) tụt xuống mức vết ($< 1,5\text{ mg/100g}$ đất), kali trao đổi suy giảm nghiêm trọng do bị rửa trôi qua tầng thẩm lậu sâu.
Phương thức canh tác khảo sát Lượng đất xói mòn (tấn/ha/năm) Mùn tầng mặt (%) $pH_{KCl}$ $Al^{3+}$ di động (meq/100g) $P_2O_5$ dễ tiêu (mg/100g)
1. Đất đồi trọc bỏ hoang $184,6$ $1,02$ $3,85$ $7,86$ $0,85$
2. Sắn độc canh (thuần loài) $268,5$ $1,15$ $3,90$ $6,45$ $1,20$
3. Sắn xen lạc + đậu đỗ $18,4$ $2,45$ $4,45$ $2,15$ $3,80$
4. Sắn xen băng cốt khí $8,5$ $2,78$ $4,60$ $1,35$ $4,50$
5. Chè thuần loài $33,0$ $2,10$ $4,10$ $4,20$ $2,90$
6. Chè + Muồng đen + Cốt khí $2,1$ $3,15$ $4,85$ $0,85$ $6,40$
7. Rừng Keo lá tràm 3 năm tuổi $4,2$ $2,90$ $4,70$ $1,10$ $4,10$

Chương 3: Tác dụng cải tạo, bồi dưỡng độ phì nhiêu đất của các hệ thống nông lâm kết hợp

  • Khả năng cung cấp sinh khối và tuần hoàn dinh dưỡng từ cây phân xanh:
    • Cây cốt khí (Tephrosia candida) khi trồng thành băng bảo vệ có khả năng sinh trưởng mạnh, cung cấp từ $20 - 35\text{ tấn}$ cành lá tươi/ha/năm, tương đương với $4 - 7\text{ tấn}$ chất khô hữu cơ trả lại cho đất.
    • Phân tích sinh hóa thân lá cốt khí cho thấy lượng dinh dưỡng bổ sung hàng năm đạt: $90 - 150\text{ kg N/ha}$, $15 - 25\text{ kg } P_2O_5\text{/ha}$, $40 - 70\text{ kg } K_2O\text{/ha}$.
    • Cây muồng đen (Cassia siamea) trồng làm cây che bóng trong nương chè cung cấp lượng lá rụng định kỳ $7 - 12\text{ tấn/ha/năm}$, tạo lớp phủ hữu cơ dày hạn chế bốc thoát hơi nước trong mùa khô và cải thiện dung trọng đất (giảm từ $1,32\text{ g/cm}^3$ xuống $1,15\text{ g/cm}^3$).
  • Cải thiện tính chất vật lý đất:
    • Các công thức nông lâm kết hợp làm tăng độ xốp tổng số của đất lên $52 - 57%$ (so với $38 - 42%$ ở đất độc canh thoái hóa).
    • Tốc độ thấm nước của đất dưới thảm cây phân xanh và nương chè nhiều tầng tán tăng gấp 2,5 - 4 lần so với đất sắn độc canh, làm giảm lưu lượng dòng chảy mặt và hạn chế hiện tượng đóng váng bề mặt.
  • Tác động đến năng suất cây trồng:
    • Trồng sắn xen cốt khí hoặc xen lạc không làm giảm năng suất củ sắn mà còn duy trì ổn định ở mức $16 - 20\text{ tấn củ tươi/ha}$ qua nhiều năm liên tục, đồng thời thu thêm $0,8 - 1,2\text{ tấn hạt lạc/ha}$.
    • Nương chè có cốt khí và muồng đen cho năng suất búp tươi tăng $25 - 40%$ so với chè độc canh không che bóng, phẩm chất búp chè được nâng cao rõ rệt.
                         HIỆU QUẢ CẢI TẠO ĐẤT CỦA CÂY PHÂN XANH

Chương 4: Ứng dụng toán học tối ưu hóa và phân loại số học xác định phương thức canh tác hợp lý

  • Thiết lập ma trận số học và phân nhóm phương thức canh tác:
    • Luận án chuẩn hóa tập hợp các chỉ tiêu (Mùn, N tổng số, $P_2O_5$ dễ tiêu, $K_2O$ trao đổi, $pH_{KCl}$, $Al^{3+}$ di động, Dung trọng, Lượng đất xói mòn) của 5 nhóm phương thức canh tác chính ($PT_1$ đến $PT_5$).
    • Thiết lập ma trận tương quan và khoảng cách hình học:
Khoảng cách ($d_{jk}$) $PT_1$ (Trọc) $PT_2$ (Sắn thuần) $PT_3$ (Sắn xen lạc) $PT_4$ (Chè + Cốt khí) $PT_5$ (Rừng Keo)
$PT_1$ $0,00$ $0,42$ $0,94$ $1,00$ $0,97$
$PT_2$ $0,42$ $0,00$ $0,61$ $0,97$ $0,87$
$PT_3$ $0,94$ $0,61$ $0,00$ $0,31$ $0,52$
$PT_4$ $1,00$ $0,97$ $0,31$ $0,00$ $0,18$
$PT_5$ $0,97$ $0,87$ $0,52$ $0,18$ $0,00$
  • Tiến trình hợp nhóm theo thuật toán phân cấp:
    1. Bước 1: Cặp $(PT_4, PT_5)$ có khoảng cách nhỏ nhất ($d_{45} = 0,18$), được gộp thành cụm $(PT_{4-5})$.
    2. Bước 2: Cụm $(PT_{4-5})$ kết hợp với $PT_3$ ở mức khoảng cách $d = 0,31$, tạo thành nhóm các phương thức canh tác cải tạo và bảo vệ đất tối ưu (Nông lâm kết hợp, xen canh họ đậu).
    3. Bước 3: Cặp $(PT_1, PT_2)$ có khoảng cách $d_{12} = 0,42$, gộp thành nhóm các phương thức suy thoái độ phì nghiêm trọng (Đất trống trọc và sắn độc canh).
    4. Bước 4: Hai nhánh chính phân tách hoàn toàn ở ngưỡng khoảng cách $d > 0,85$, chứng minh tính đối lập tuyệt đối về mặt nông học và thổ nhưỡng giữa canh tác quảng canh phá hủy thảm phủ và canh tác nông lâm kết hợp bảo vệ đất.
Dendrogram phân loại số học các phương thức canh tác (Sokal & Sneath):

Mức khoảng cách (d)
0.0           0.2           0.4           0.6           0.8           1.0

Kết quả và những đóng góp mới

Đóng góp về mặt lý luận khoa học

  1. Lượng hóa quy luật thoái hóa đất đồi feralit: Xác định định lượng tốc độ suy giảm các hợp phần độ phì (mùn, đạm, lân, độ chua, dung trọng) và chỉ ra mối tương quan tỷ lệ thuận giữa hàm lượng nhôm di động ($Al^{3+}$) với mức độ mất đất do xói mòn khi độc canh sắn trên sườn dốc.
  2. Phương pháp luận phân loại số học trong thổ nhưỡng: Lần đầu tiên tại Việt Nam, lý thuyết Phân loại số học (Numerical Taxonomy) và các thuật toán phân cụm đa biến (Lance & Williams) được áp dụng thành công để phân loại, đánh giá khách quan các hệ thống canh tác nông lâm kết hợp dựa trên tập số liệu thổ nhưỡng thực nghiệm, khắc phục tính chủ quan của các phương pháp đánh giá định tính trước đây.

Đóng góp về mặt thực tiễn và giải pháp đề xuất

  1. Xác lập các mô hình canh tác tối ưu cho vùng trung du Vĩnh Phú:
    • Mô hình Chè + Muồng đen + Cốt khí: Giảm xói mòn $98%$, tăng hàm lượng mùn tầng mặt lên trên $3,0%$, nâng cao năng suất búp tươi $25 - 40%$.
    • Mô hình Sắn xen Lạc / Cốt khí theo đường đồng mức: Giảm lượng đất mất từ $268,5\text{ tấn/ha/năm}$ xuống dưới $18,4\text{ tấn/ha/năm}$, duy trì năng suất sắn ổn định $16 - 20\text{ tấn/ha}$ và tạo thêm nguồn thu từ cây đậu đỗ.
  2. Kiến nghị phân vùng kỹ thuật canh tác đất dốc:
    • Sườn dốc $< 15^\circ$: Cho phép canh tác cây lương thực và cây công nghiệp ngắn ngày nhưng bắt buộc phải áp dụng biện pháp xen canh cây họ đậu, tạo băng cốt khí chắn đất theo đường đồng mức và bón lân nung chảy kết hợp tàn dư hữu cơ.
    • Sườn dốc $> 15^\circ$: Nghiêm cấm độc canh cây ngắn ngày; chuyển đổi toàn bộ sang trồng chè nhiều tầng tán, nông lâm kết hợp hoặc trồng rừng bảo vệ đất (keo lá tràm, bồ đề, muồng đen).

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án

  1. Thời gian quan trắc thực địa: Thời gian theo dõi thực nghiệm đồng ruộng (1985 - 1988) kéo dài 4 năm, tương đối ngắn so với chu kỳ biến đổi lâu dài của các quá trình hình thành và thoái hóa thổ nhưỡng nhiệt đới.
  2. Hạn chế về công cụ phân tích thời kỳ 1990: Các chỉ tiêu vi lượng (Bo, Mo, Zn, Cu) và động thái chi tiết của các phân đoạn mùn (axit humic, axit fulvic) chưa được định lượng đầy đủ do điều kiện trang thiết bị phân tích còn thiếu thốn.
  3. Phạm vi địa bàn: Các ô thí nghiệm tập trung chủ yếu trên nhóm đất feralit phát triển trên đá phiến thạch sét tại huyện Thanh Ba, chưa bao quát hết các nhóm đất feralit trên đá biến chất hay đá vôi của toàn vùng trung du miền núi phía Bắc.

Hướng nghiên cứu tiếp

  • Mở rộng chuỗi quan trắc dài hạn (10 - 20 năm) nhằm đánh giá sự tích lũy hữu cơ bền vững và chu trình carbon trong các hệ thống nông lâm kết hợp.
  • Ứng dụng mô hình toán học đa biến mở rộng trên các vùng sinh thái đất dốc khác như Tây Bắc, Tây Nguyên và Bắc Trung Bộ.

Giá trị tham khảo

  • Đối với nghiên cứu sinh, học viên và giảng viên ngành Khoa học Đất, Nông học, Lâm nghiệp: Luận án cung cấp hệ thống số liệu thực nghiệm gốc chi tiết về xói mòn, biến động lý hóa tính đất feralit và là tài liệu tham khảo mẫu mực về phương pháp ứng dụng toán học đa biến (Taxonomy số học) trong phân tích số liệu nông nghiệp.
  • Đối với các nhà hoạch định chính sách và quản lý nông nghiệp vùng đồi núi: Cung cấp cơ sở khoa học để xây dựng các quy hoạch sử dụng đất dốc bền vững, ban hành các quy chuẩn kỹ thuật bắt buộc về trồng băng cây cải tạo đất và xóa bỏ tập quán độc canh sắn trên đất dốc.

Câu hỏi thường gặp

1. Tác hại nghiêm trọng nhất của việc trồng sắn độc canh trên đất đồi dốc Thanh Ba là gì?

Trồng sắn độc canh không che phủ làm lượng đất xói mòn lên tới $120 - 366\text{ tấn/ha/năm}$ (làm mất $0,8 - 1,5\text{ cm}$ tầng đất mặt mỗi năm), mùn giảm dưới $1,2%$, dung trọng đất tăng cao làm đất bị nén chặt, và hàm lượng $Al^{3+}$ di động tăng vọt lên $5,7 - 7,9\text{ meq/100g}$ đất, gây chua hóa nặng và đầu độc hệ rễ cây trồng.

2. Cây cốt khí (Tephrosia candida) có vai trò cụ thể như thế nào trong việc cải tạo đất đồi?

Cốt khí sinh trưởng nhanh trên đất đồi chua nghèo dinh dưỡng, cung cấp $20 - 35\text{ tấn}$ sinh khối cành lá tươi/ha/năm (tương đương trả lại cho đất $90 - 150\text{ kg N/ha}$, $15 - 25\text{ kg } P_2O_5\text{/ha}$, $40 - 70\text{ kg } K_2O\text{/ha}$), che phủ giảm xói mòn trên $85%$, tăng mùn tầng mặt thêm $0,45 - 0,77%$ và giảm mạnh hàm lượng nhôm di động.

3. Mô hình toán học nào được tác giả sử dụng để tối ưu hóa việc phân loại các phương thức canh tác?

Tác giả sử dụng mô hình Phân loại số học (Numerical Taxonomy) của Sokal và Sneath (1963) kết hợp thuật toán tính ma trận khoảng cách hình học và phân cụm phân cấp của Lance và Williams (1967) được lập trình xử lý trên máy tính điện tử.

4. Lượng mưa tại khu vực nghiên cứu phân bố như thế nào và ảnh hưởng ra sao đến quá trình thoái hóa đất?

Theo số liệu 10 năm tại trạm Phú Hộ, tổng lượng mưa trung bình năm đạt $1683,6\text{ mm}$, nhưng $84,9%$ lượng mưa tập trung dồn dập vào mùa mưa (tháng 5 đến tháng 10). Cường độ mưa lớn kết hợp độ dốc địa hình tạo nên động năng dòng chảy mặt rất mạnh, gây xói mòn cơ học khốc liệt trên các nương đồi thiếu thảm phủ thực vật.

Kết luận

Luận án của nghiên cứu sinh (Hà Nội, 1990) đã giải quyết toàn diện vấn đề đánh giá tác động của các phương thức canh tác đến độ phì nhiêu đất đồi feralit tại huyện Thanh Ba, tỉnh Vĩnh Phú thông qua các dữ liệu thực nghiệm đồng ruộng và phòng thí nghiệm tin cậy. Nghiên cứu đã chứng minh một cách định lượng tác hại phá hủy đất nghiêm trọng của tập quán độc canh sắn và khẳng định vai trò vượt trội của các mô hình nông lâm kết hợp, nương chè có cây che bóng và xen canh cây họ đậu trong việc bảo vệ, phục hồi độ phì nhiêu đất. Việc ứng dụng thành công mô hình toán học phân loại số học Sokal - Sneath không chỉ mang lại giá trị lý luận học thuật sâu sắc mà còn cung cấp cơ sở khoa học vững chắc cho chiến lược phát triển nông lâm nghiệp bền vững trên đất đồi dốc Việt Nam.