Giới thiệu dự án

Nghiên cứu đặc tính vật lý đất đóng vai trò cốt lõi trong công tác quản lý tài nguyên rừng bền vững và quy hoạch lâm sinh. Các nghiên cứu quốc tế từ Tổ chức Nông - Lương Liên Hợp Quốc (FAO) và Trung tâm Lâm nghiệp Quốc tế (CIFOR) chỉ ra rằng việc chuyển đổi rừng tự nhiên sang rừng trồng thuần loài cây mọc nhanh (như Thông, Bạch đàn, Keo) với chu kỳ ngắn có thể làm suy giảm 15–30% độ phì nhiêu và làm suy thoái lý tính đất tầng mặt nếu không có biện pháp giám sát lập địa định kỳ.

+-----------------------------------------------------------------------------+
|                      BỐI CẢNH & VẤN ĐỀ NGHIÊN CỨU                           |
+-----------------------------------------------------------------------------+
| Rừng trồng Thông mã vĩ (Pinus massoniana Lamb) tại núi Nhật Lệ, Thanh Hóa   |
|                                                                             |
| [Phương pháp truyền thống]                     [Phương pháp đề xuất]        |
| - Đào phẫu diện, lấy mẫu ống dung trọng 100cm3 - Đo đâm xuyên nhanh         |
| - Sấy khô kiệt 105-110°C (6-8 giờ)             - Dụng cụ Daiki Cone 5553    |
| - Đun sôi bình tỷ trọng picnometer             - Áp dụng mô hình hồi quy    |
| -> Tốn 24-48h/mẫu, chi phí cao, phá hủy cấu trúc -> Đo đạc tức thì (3 phút) |
+-----------------------------------------------------------------------------+

Vấn đề kỹ thuật lớn nhất đặt ra là việc xác định độ xốp của đất ($P$) theo phương pháp cổ điển đòi hỏi quy trình phức tạp trong phòng thí nghiệm: sấy mẫu xác định dung trọng ($D$) và đun sôi trong bình tỷ trọng ($d$) để loại bỏ bọt khí mao quản. Phương pháp này tiêu tốn từ 24–48 giờ cho mỗi đợt phân tích, chi phí hóa chất và thiết bị cao, không đáp ứng được yêu cầu đánh giá nhanh trên quy mô diện rộng.

Dự án tập trung giải quyết bài toán thiết lập mô hình toán học thực nghiệm biểu diễn mối quan hệ giữa độ chặt đâm xuyên ($H_{cone}$), độ ẩm tự nhiên ($W$) và độ xốp ($P$) của đất dưới tán rừng Thông mã vĩ (Pinus massoniana Lamb) tại núi Nhật Lệ, xã Quảng Vinh, huyện Quảng Xương, tỉnh Thanh Hóa.

Mục tiêu dự án

  1. Xác định chính xác các chỉ số lý tính cơ bản của đất tầng mặt (0–20 cm): độ chặt đâm xuyên ($H_{cone}$), độ ẩm tự nhiên ($W$), dung trọng ($D$), tỷ trọng ($d$) và độ xốp ($P$).
  2. Thiết lập phương trình hồi quy tương quan đơn biến giữa độ chặt và độ ẩm; giữa độ chặt và độ xốp; giữa độ ẩm và độ xốp.
  3. Xây dựng mô hình hồi quy đa biến biểu diễn mối quan hệ đồng thời giữa độ ẩm, độ chặt với độ xốp đất.
  4. Kiểm định mức độ phù hợp và độ tin cậy thống kê của mô hình thông qua hệ số xác định ($R^2$) và phân tích phương sai ANOVA ($p$-value $< 0.05$).
  5. Xây dựng bảng tra cứu nhanh độ xốp đất dựa trên giá trị độ chặt đo trực tiếp tại hiện trường.

Phạm vi và giới hạn

  • Phạm vi không gian: Rừng trồng Thông mã vĩ trên đất feralit cát pha/thịt nhẹ tại khu vực núi Nhật Lệ, xã Quảng Vinh, huyện Quảng Xương, tỉnh Thanh Hóa.
  • Phạm vi đối tượng: Tầng đất mặt (độ sâu 0–20 cm) chịu tác động trực tiếp của lớp thảm mục và rễ cây hấp thu.
  • Giới hạn kỹ thuật: Mô hình hồi quy phụ thuộc vào điều kiện thổ nhưỡng đồng nhất của khu vực nghiên cứu; cần hiệu chỉnh hệ số khi áp dụng trên các nhóm đất có thành phần cơ giới nặng (đất sét) hoặc hàm lượng hữu cơ biến động lớn.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí Phương pháp cổ điển (Picnometer + Tủ sấy) Phương pháp thước rơi tự do (GS. Vương Văn Quỳnh) Phương pháp đo đâm xuyên hình nón Daiki Cone
Nguyên lý Cân đo khối lượng sấy khô kiệt và thể tích bị chiếm chỗ trong nước cất Thả thanh sắt 10mm, dài 1m rơi tự do từ độ cao 50cm ($X = 0.52 \cdot ST$) Đo lực kháng xuyên cơ học bằng mũi nón nén lò xo (Daiki Push-Cone 5553)
Thời gian thực hiện 24–48 giờ/mẫu 5–10 phút/ô tiêu chuẩn 1–3 phút/vị trí đo
Độ chính xác ($R^2$) Chuẩn tuyệt đối (Lab Benchmark) $R^2 \approx 0.35 - 0.42$ $R^2 = 0.5176$ (Mô hình bậc 2)
Độ tiện dụng thực địa Kém (phải lấy mẫu đất nguyên thổ đóng hộp mang về) Trung bình (thiết bị cồng kềnh, dễ lệch phương thẳng đứng) Cao (thiết bị cầm tay nhỏ gọn, hiển thị số đo trực tiếp)
Chi phí thiết bị Cao (Tủ sấy, cân 4 số, bình tỷ trọng Glassco) Rất thấp (Gia công cơ khí thủ công) Trung bình (Thiết bị tiêu chuẩn công nghiệp)

Phân loại yêu cầu kỹ thuật (MoSCoW)

  • Must-have (Bắt buộc):
    • Đo đạc chính xác độ đâm xuyên đất tại hiện trường bằng dụng cụ Daiki Push-Cone 5553 với đơn vị mm.
    • Phân tích độ ẩm ($W$), dung trọng ($D$), tỷ trọng ($d$) trong phòng thí nghiệm theo tiêu chuẩn TCVN/FAO.
    • Kiểm định ANOVA với mức ý nghĩa $\alpha = 0.05$ (Sig $< 0.05$).
  • Should-have (Nên có):
    • Mô hình hóa phi tuyến bậc 2 cho cặp biến độ chặt – độ xốp để tối ưu hệ số $R^2$.
    • Lọc bỏ sai số bất thường (outliers) trước khi phân tích hồi quy.
  • Could-have (Có thể có):
    • Xây dựng bảng tra cứu hai chiều (tra cứu theo cả độ ẩm và độ chặt).
  • Won't-have (Chưa thực hiện đợt này):
    • Tự động hóa truyền dữ liệu không dây từ đầu dò đất vào hệ thống cơ sở dữ liệu GIS thời gian thực.
+-----------------------------------------------------------------------------+
|                      SƠ ĐỒ QUY TRÌNH THỰC NGHIỆM                            |
+-----------------------------------------------------------------------------+
| [105 Điểm Lưới 5x5m]                                                       |
|        |                                                                    |
|        +---> Đo độ chặt tại hiện trường (Daiki Push-Cone 5553) -> [H_cone]  |
|        |                                                                    |
|        +---> Lấy mẫu nguyên thổ (Ống dung trọng 100cm3)                    |
|                    |                                                        |
|                    v                                                        |
|        [Phân tích Phòng Lab]                                                |
|        +--- Tủ sấy Memmert UNB-200 (105-110°C) -> Độ ẩm W (%), Dung trọng D |
|        +--- Bình tỷ trọng Glassco + Bếp sỏi    -> Tỷ trọng d (g/cm3)        |
|                    |                                                        |
|                    v                                                        |
|        [Xử lý Thống kê & Mô hình hóa]                                       |
|        +--- Tính độ xốp: P = (1 - D/d) * 100                                |
|        +--- Lọc Outliers (N=105 -> N=100)                                   |
|        +--- SPSS 20.0 / Python Scipy: Hồi quy Tuyến tính & Phi tuyến        |
|                    |                                                        |
|                    v                                                        |
|        [Bảng tra thực nghiệm P = f(H_cone)]                                 |
+-----------------------------------------------------------------------------+

Thông số kỹ thuật thiết bị và công cụ

  • Thiết bị hiện trường: Dụng cụ đo độ chặt đâm xuyên hình nón Daiki Push-Cone 5553 (Daiki Rika Kogyo Co., Nhật Bản), dải đo 0–40 mm. Ống dung trọng thép không gỉ thể tích chuẩn $V = 100 \text{ cm}^3$.
  • Thiết bị phân tích Lab:
    • Cân phân tích điện tử Ohaus Pioneer PA213 (độ chính xác $0.0001\text{ g}$).
    • Tủ sấy đối lưu nhiệt Memmert UNB-200 (nhiệt độ sấy ổn định $105^\circ\text{C} - 110^\circ\text{C}$).
    • Bình hút ẩm chân không Model A036-01.
    • Bình đo tỷ trọng định mức Glassco tiêu chuẩn ISO 3507.
  • Phần mềm tính toán: IBM SPSS Statistics phiên bản 20.0, Microsoft Excel 2016, Python 3.10 (thư viện Pandas, SciPy, NumPy, Matplotlib).

Implementation và kết quả

Quy trình phân tích thực nghiệm

Quy trình thu thập và phân tích mẫu đất được thực hiện nghiêm ngặt qua 3 giai đoạn:

  1. Giai đoạn thu thập mẫu hiện trường:

    • Thiết lập mạng lưới ô mẫu kích thước $5\text{ m} \times 5\text{ m}$ dưới tán rừng Thông mã vĩ núi Nhật Lệ.
    • Thu thập 105 mẫu đất tại các giao điểm lưới vào cùng khung giờ buổi sáng để đồng nhất điều kiện vi khí hậu.
    • Tại mỗi điểm: Cắm thiết bị Daiki Push-Cone 5553 theo phương thẳng đứng, đọc độ nén lò xo ($H_{cone}$ tính bằng mm). Cách vị trí đo 10 cm, đóng ống dung trọng $100\text{ cm}^3$ ngập hoàn toàn, đào cắt phẳng 2 đầu, niêm phong trong túi zip chuyên dụng.

  2. Giai đoạn phân tích phòng thí nghiệm:

    • Độ ẩm ($W$) và Dung trọng ($D$): Cân $20\text{ g}$ đất tươi ($M_2$), sấy trong tủ Memmert UNB-200 ở $105^\circ\text{C} - 110^\circ\text{C}$ liên tục 6–8 giờ. Chuyển vào bình hút ẩm A036-01 trong 20 phút, cân lại khối lượng ($M_3$). Lặp lại chu kỳ sấy 2 giờ đến khi khối lượng không đổi ($\Delta m \le 0.0001\text{ g}$). $$\text{Dung trọng } D = \frac{M_{kk}}{V} = \frac{M_3 - M_1}{100} \quad (\text{g/cm}^3)$$ $$\text{Độ ẩm } W = \frac{M_2 - M_{kk}}{M_{kk}} \times 100 \quad (%)$$
    • Tỷ trọng ($d$): Đất sau xử lý rây qua mắt sàng $1\text{ mm}$, cân $10\text{ g}$ đất khô kiệt cho vào bình tỷ trọng Glassco chứa $1/3$ thể tích nước cất, đun sôi trên bếp sỏi 5 phút để đẩy toàn bộ bọt khí mao quản. Làm nguội, châm nước cất đến vạch mức, cân khối lượng tổng hợp. $$\text{Tỷ trọng } d = \frac{m_0}{(m_1 - m_2) + m_0} \quad (\text{g/cm}^3)$$
    • Độ xốp tổng số ($P$): $$P = \left(1 - \frac{D}{d}\right) \times 100 \quad (%)$$

  3. Thuật toán xử lý và phân tích hồi quy (Python/SPSS):

import numpy as np
import pandas as pd
from scipy import stats

def analyze_soil_metrics(data_path: str):
    # Đọc dữ liệu thực nghiệm 105 mẫu
    df = pd.read_csv(data_path)
    
    # 1. Tính toán độ xốp P (%)
    df['P_porosity'] = (1.0 - (df['D_bulk_density'] / df['d_specific_gravity'])) * 100.0
    
    # 2. Lọc bỏ Outliers bằng Z-score (|z| < 3.0) -> N = 100
    z_scores = np.abs(stats.zscore(df[['H_cone', 'W_moisture', 'P_porosity']]))
    clean_df = df[(z_scores < 3.0).all(axis=1)].copy()
    
    # 3. Khớp mô hình hồi quy bậc 2: P = a*H^2 + b*H + c
    x = clean_df['H_cone'].values
    y = clean_df['P_porosity'].values
    poly_coeffs = np.polyfit(x, y, deg=2)
    p_model = np.poly1d(poly_coeffs)
    
    # 4. Tính toán R-squared
    y_pred = p_model(x)
    ss_res = np.sum((y - y_pred) ** 2)
    ss_tot = np.sum((y - np.mean(y)) ** 2)
    r2_score = 1.0 - (ss_res / ss_tot)
    
    # 5. Phân tích ANOVA
    f_stat, p_val = stats.f_oneway(y, y_pred)
    
    return {
        "Coefficients": poly_coeffs,
        "R2": r2_score,
        "ANOVA_F": f_stat,
        "ANOVA_p_val": p_val,
        "Clean_Sample_Size": len(clean_df)
    }

Kết quả phân tích thống kê và kiểm định

Sau khi loại bỏ 5 giá trị ngoại biên (outliers), bộ dữ liệu $N = 100$ mẫu hợp lệ cho kết quả kiểm định thống kê như sau:

Cặp biến tương quan Dạng mô hình toán học Phương trình thực nghiệm Hệ số $R^2$ Hệ số $R$ Mức ý nghĩa ($p$-value) Đánh giá thống kê
Độ chặt ($x$) - Độ ẩm ($y$) Tuyến tính $y = -0.208x + 38.45$ $0.1056$ $0.325$ $0.001 < 0.05$ Tương quan rất yếu, không phù hợp lập bảng
Độ ẩm ($x$) - Độ xốp ($y$) Tuyến tính $y = 0.048x + 46.12$ $0.0616$ $0.248$ $0.013 < 0.05$ Không có ý nghĩa thực tiễn dự báo
Độ chặt ($x$) - Độ xốp ($y$) Tuyến tính $y = -0.842x + 61.20$ $0.4680$ $0.684$ $0.000 < 0.001$ Tương quan nghịch chặt chẽ
Độ chặt ($x$) - Độ xốp ($y$) Bậc hai (Tối ưu) $y = 0.0185x^2 - 1.342x + 65.51$ $0.5176$ $0.719$ $0.000 < 0.001$ Mô hình tối ưu, đủ tin cậy lập bảng tra
Đa biến ($x_1$: Chặt, $x_2$: Ẩm - $Y$: Xốp) Đa biến tuyến tính $Y = -1.465x_1 + 0.021x_2 + 63.82$ $0.4830$ $0.695$ $0.000 < 0.001$ $x_2$ không cải thiện $R^2$ đáng kể

Theo thang đánh giá đất của Katrinski:

  • Dung trọng đất rừng Thông mã vĩ dao động từ $1.15 - 1.38\text{ g/cm}^3$ (thuộc nhóm đất tơi xốp trung bình, điển hình cho đất trồng trọt và lâm nghiệp).
  • Tỷ trọng đất đạt $2.55 - 2.68\text{ g/cm}^3$ (phản ánh hàm lượng mùn trung bình, đất không bị glêi hóa mạnh).
  • Độ xốp tổng số đạt $42.5% - 56.8%$ (phân cấp đất xốp trung bình đến khá xốp).

Bảng tra thực nghiệm độ xốp đất rừng Thông mã vĩ theo độ chặt

Dựa trên phương trình bậc hai $y = 0.0185x^2 - 1.342x + 65.51$, bảng tra thực nghiệm ngoài hiện trường được xây dựng cho dải độ chặt đâm xuyên $H_{cone}$ từ $9.0\text{ mm}$ đến $22.0\text{ mm}$:

Độ chặt $H_{cone}$ (mm) Độ xốp đất $P$ (%) Phân cấp độ xốp Katrinski Độ chặt $H_{cone}$ (mm) Độ xốp đất $P$ (%) Phân cấp độ xốp Katrinski
9.0 54.93 Khá xốp 16.0 48.78 Xốp trung bình
10.0 53.94 Khá xốp 17.0 48.04 Xốp trung bình
11.0 53.00 Khá xốp 18.0 47.35 Xốp trung bình
12.0 52.09 Khá xốp 19.0 46.70 Xốp trung bình
13.0 51.22 Khá xốp 20.0 46.10 Xốp trung bình
14.0 50.38 Khá xốp 21.0 45.54 Xốp trung bình
15.0 49.56 Xốp trung bình 22.0 45.03 Xốp trung bình

Đổi mới và đóng góp

  1. Đổi mới phương pháp luận đo đạc lập địa:

    • Thay thế hoàn toàn phương pháp rơi tự do thủ công không kiểm soát được góc rơi bằng thiết bị chuẩn hóa cơ học Daiki Push-Cone 5553.
    • Giảm thiểu $85%$ thời gian thu thập và phân tích dữ liệu lý tính đất: từ $36\text{ giờ/lô}$ xuống còn $15\text{ phút/lô}$ điều tra tại hiện trường.

  2. Đóng góp về mặt mô hình toán học:

    • Khám phá tính chất phi tuyến bậc hai giữa lực kháng đâm xuyên và độ xốp tổng số ($R^2 = 0.5176$), chứng minh rằng trong điều kiện đất đồi núi nhiệt đới có thảm mục Thông mã vĩ, độ ẩm tự nhiên tức thời không chi phối đáng kể đến độ xốp cấu trúc đất ($R^2 = 0.0616$).
    • Bác bỏ giả thiết cần bảng tra hai chiều phức tạp, tinh gọn quy trình đánh giá chỉ qua một biến số đầu vào duy nhất là độ đâm xuyên $H_{cone}$.

  3. Đóng góp cơ sở dữ liệu lâm học:

    • Bổ sung chuỗi số liệu $100$ mẫu chuẩn về dung trọng ($D$), tỷ trọng ($d$), độ chặt và độ xốp ($P$) của đất rừng Thông mã vĩ vùng ven biển Thanh Hóa vào cơ sở dữ liệu thổ nhưỡng lâm nghiệp Việt Nam.

Ứng dụng thực tế và triển khai

Tình huống ứng dụng thực tế (Use Cases)

  • Giám sát suy thoái đất rừng sau khai thác: Các công ty lâm nghiệp có thể sử dụng Daiki Push-Cone để kiểm tra mức độ nén dẽ đất do xe cơ giới khai thác gỗ gây ra; nếu độ chặt $> 20\text{ mm}$ (độ xốp $< 46%$), cần tiến hành xới xáo cải tạo trước khi trồng luân kỳ mới.
  • Quy hoạch loài cây trồng lâm nghiệp: Xác định nhanh tầng đất mặt có đủ độ xốp tối thiểu ($> 50%$) cho các loài cây rễ cọc nhạy cảm với độ nén dẽ như Thông, Lát hoa, Lim xanh.
+-----------------------------------------------------------------------------+
|                      QUY TRÌNH TRIỂN KHAI THỰC ĐỊA                          |
+-----------------------------------------------------------------------------+
| BƯỚC 1: Dọn sạch lớp thảm mục mặt đất (bán kính 20cm)                      |
|         Vị trí: Đất rừng Thông mã vĩ hoặc đất lâm nghiệp quy hoạch          |
|                               |                                             |
|                               v                                             |
| BƯỚC 2: Cắm đầu nón Daiki Cone vuông góc 90° xuống mặt đất                  |
|         Ấn đều tay đến khi nón ngập hoàn toàn                               |
|                               |                                             |
|                               v                                             |
| BƯỚC 3: Đọc trị số nén lò xo (mm) trên thân thước đo                        |
|         Ghi nhận giá trị H_cone                                             |
|                               |                                             |
|                               v                                             |
| BƯỚC 4: Tra cứu trực tiếp trên Bảng thực nghiệm                             |
|         -> Xác định ngay Độ xốp P (%) và Phân cấp Katrinski                 |
+-----------------------------------------------------------------------------+

Phân tích chi phí - hiệu quả (Cost-Benefit & ROI)

  • Chi phí phương pháp Lab truyền thống: 100 mẫu $\times$ $150.000\text{ VNĐ/mẫu}$ (tiền điện tủ sấy, hóa chất, khấu hao cân, công phân tích) = $15.000.000\text{ VNĐ}$; thời gian xử lý: 15–20 ngày làm việc.
  • Chi phí phương pháp đo nhanh Daiki + Bảng tra: Chi phí đầu tư thiết bị ban đầu: $8.000.000\text{ VNĐ}$ (sử dụng $> 5\text{ năm}$); chi phí khảo sát: 1 ngày công thực địa ($500.000\text{ VNĐ}$).
  • Điểm hòa vốn (ROI): Hoàn vốn đầu tư thiết bị ngay sau 2 đợt điều tra lập địa quy mô 50 ha; tiết kiệm $90%$ kinh phí phân tích cho các ban quản lý rừng phòng hộ.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Độ nhạy lập địa: Phương trình hồi quy $y = 0.0185x^2 - 1.342x + 65.51$ được xây dựng trên nền đất feralit cát pha dưới tán rừng Thông mã vĩ núi Nhật Lệ; khi áp dụng trên đất sét nặng hoặc đất có lẫn nhiều đá mẹ phong hóa cứng ($> 2\text{ mm}$), đầu đo nón dễ gặp vật cản cục bộ tạo sai số cơ học.
  • Ảnh hưởng của mạng lưới rễ cây: Đất tầng mặt rừng thông có mật độ rễ cám cao có thể làm tăng chỉ số kháng đâm xuyên mà không phản ánh đúng dung trọng thực của khoáng vật đất.

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

  • Mở rộng thu thập mẫu trên các trạng thái rừng khác: Rừng Keo tai tượng (Acacia mangium), Keo lá tràm (Acacia auriculiformis), và rừng tự nhiên phục hồi.
  • Tích hợp cảm biến đo ẩm điện dung (FDR/TDR) trực tiếp vào đầu đo Daiki Cone, kết nối Bluetooth với điện thoại thông minh để tự động tính toán và vẽ bản đồ phân bố độ xốp đất trên nền tảng QGIS.

Đối tượng hưởng lợi

  • Sinh viên & Giảng viên ngành Lâm học, Quản lý tài nguyên rừng: Tiếp cận phương pháp nghiên cứu thực nghiệm chuẩn tắc, bộ dữ liệu thực tế minh họa cho giáo trình Khoa học đất và Đánh giá lập địa.
  • Kỹ sư lâm nghiệp & Cán bộ điều tra quy hoạch rừng: Sở hữu công cụ bảng tra cứu hiện trường tiện lợi, cắt giảm $85%$ thời gian lập hồ sơ đất lập địa phục vụ các dự án trồng rừng thay thế.
  • Ban quản lý rừng phòng hộ & Doanh nghiệp lâm nghiệp: Tối ưu hóa chi phí giám sát chất lượng đất định kỳ, phát hiện sớm hiện tượng nén dẽ đất để có biện pháp can thiệp cơ giới hóa kịp thời.
  • Nhà nghiên cứu thổ nhưỡng & Sinh thái môi trường: Nắm bắt mô hình toán học thực nghiệm về mối quan hệ giữa các tính chất lý tính đất nhiệt đới, làm tiền đề mở rộng cho các hệ sinh thái rừng khác tại Việt Nam.

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

1. Yêu cầu kỹ thuật tối thiểu để triển khai bảng tra độ xốp ngoài thực địa là gì?

Cần trang bị thiết bị đo độ chặt hình nón nén lò xo chuẩn Daiki Push-Cone 5553 (hoặc các dòng penetrometer có thông số kỹ thuật tương đương: góc mũi nón $30^\circ$, tiết diện đáy nón chuẩn), thước chia vạch mm và bảng tra thực nghiệm đi kèm. Vị trí đo phải được gạt sạch lớp thảm mục và lá thông khô bề mặt.

2. Tại sao nghiên cứu loại bỏ yếu tố độ ẩm khi xây dựng bảng tra độ xốp?

Kết quả hồi quy đa biến và kiểm định ANOVA cho thấy hệ số tương quan giữa độ ẩm và độ xốp rất thấp ($R^2 = 0.0616$). Khi đưa độ ẩm vào mô hình đa biến ($Y = -1.465x_1 + 0.021x_2 + 63.82$), hệ số xác định tổng thể $R^2$ chỉ đạt $0.483$, không cao hơn mô hình đơn biến bậc hai của độ chặt ($R^2 = 0.5176$). Do đó, việc loại bỏ biến độ ẩm giúp tinh giản quy trình đo mà vẫn đảm bảo độ chính xác thống kê.

3. Phương pháp này có thể áp dụng cho đất trồng lúa nước hoặc đất nông nghiệp không?

Không nên áp dụng trực tiếp bảng tra này cho đất nông nghiệp hoặc lúa nước. Đất lúa nước có tầng đế cày nén chặt nhân tạo và trạng thái bão hòa nước thường xuyên, làm thay đổi hoàn toàn cơ chế kháng xuyên cơ học. Muốn áp dụng cần lấy mẫu thực nghiệm riêng để xác lập phương trình hồi quy đặc thù cho từng loại đất.

4. Xử lý thế nào khi mũi đo đâm xuyên gặp sỏi đá hoặc rễ cây lớn?

Nếu trong quá trình ấn đầu đo, kim chỉ thị tăng đột ngột hoặc cảm nhận gặp vật cản cứng (sỏi, đá dăm, rễ cây lớn $> 5\text{ mm}$), người đo phải hủy kết quả lần đo đó, di chuyển điểm đo sang vị trí bên cạnh cách $15–20\text{ cm}$ và tiến hành đo lại để tránh tạo ra giá trị ngoại biên (outlier).

5. Chi phí trang bị dụng cụ Daiki Push-Cone 5553 và thời gian hoàn vốn?

Một bộ dụng cụ Daiki Push-Cone 5553 nhập khẩu có giá dao động từ 7.500.000 – 9.000.000 VNĐ. Với đơn vị điều tra đất lâm nghiệp trung bình 100 mẫu/năm, chi phí tiết kiệm được từ việc không phải sấy mẫu và phân tích tỷ trọng trong phòng lab giúp hoàn vốn đầu tư chỉ sau 1 đến 2 dự án khảo sát.

Kết luận

Đề tài khóa luận đã giải quyết thành công bài toán đánh giá nhanh tính chất vật lý đất tầng mặt dưới tán rừng Thông mã vĩ (Pinus massoniana Lamb) tại núi Nhật Lệ, huyện Quảng Xương, tỉnh Thanh Hóa. Bằng việc kết hợp thực nghiệm hiện trường với phân tích chuẩn tắc trong phòng thí nghiệm trên hệ thống thiết bị tiêu chuẩn (Daiki Cone 5553, Memmert UNB-200, Ohaus PA213, Glassco), nghiên cứu đã xác lập mô hình hồi quy phi tuyến bậc hai tối ưu:

$$y = 0.0185x^2 - 1.342x + 65.51 \quad (R^2 = 0.5176; \ p < 0.001)$$

Mô hình chứng minh lực kháng đâm xuyên là chỉ số tin cậy để dự báo độ xốp tổng số của đất, cho phép thay thế quy trình phân tích dung trọng - tỷ trọng kéo dài 24–48 giờ bằng thao tác đo đạc 3 phút tại hiện trường. Bảng tra cứu thực nghiệm thu được là công cụ hữu ích, đóng góp trực tiếp vào công tác quản lý lập địa, đánh giá thoái hóa đất và quy hoạch phát triển rừng bền vững tại khu vực Bắc Trung Bộ.