Giới thiệu dự án

Bối cảnh và thực trạng ngành

Rừng phòng hộ ven biển tại khu vực Bắc Trung Bộ đóng vai trò như lá chắn sinh thái tự nhiên chống lại hiện tượng cát bay, cát chảy, bão biển và xói lở bờ biển. Tại huyện Quỳnh Lưu, tỉnh Nghệ An, tổng diện tích đất quy hoạch lâm nghiệp đạt 22.371,2 ha (chiếm 36,8% diện tích tự nhiên toàn huyện), trong đó rừng phòng hộ là bức tường bảo vệ trực tiếp hệ thống đê điều, các công trình ven biển và nguồn nước ngọt của hồ thủy lợi Vực Mấu.

Cây Phi Lao (Casuarina equisetifolia) là loài cây tiên phong chủ lực được ưu tiên phát triển nhờ khả năng chịu hạn, chịu mặn và chắn gió vượt trội. Tuy nhiên, theo số liệu từ Bộ Nông nghiệp và Phát triển Nông thôn, diện tích rừng trồng mới cả nước tuy tăng nhanh nhưng chất lượng rừng có xu hướng giảm do chịu áp lực lớn từ biến đổi khí hậu và sự bùng phát phức tạp của các loài sâu hại phá hoại thân, lá và rễ cây.

+-------------------------------------------------------------------------+
|                  HỆ THỐNG RỪNG PHÒNG HỘ VEN BIỂN                        |
|                                                                         |
|  [Đai rừng Phi Lao]  <--->  [Chắn cát bay / Cản sóng / Giảm bão]        |
|          |                                                              |
|          v (Nguy cơ suy thoái)                                          |
|  [Sâu đục thân & Xén tóc]  --->  [Gãy đổ / Giảm khả năng phòng hộ]      |
|          |                                                              |
|          v (Giải pháp ứng dụng)                                         |
|  [Mô hình giám sát OTC + Kiểm định U-Test + Chiến lược IPM đa tầng]     |
+-------------------------------------------------------------------------+

Vấn đề nghiên cứu cụ thể

Tại rừng phòng hộ xã Quỳnh Nghĩa (diện tích khu vực nghiên cứu 21 ha thuộc Tiểu khu 340H), tình trạng sâu hại diễn biến âm thầm nhưng gây tổn thương nghiêm trọng đến cấu trúc lâm phần:

  • Sự xâm hại của các loài sâu đục thân làm rỗng ruột cây, gãy ngang gióng khi có gió bão cấp 6-8.
  • Thiếu hụt cơ sở dữ liệu định lượng về mật độ, chu kỳ phát sinh và tập tính của các loài côn trùng gây hại chính (Zeuzera coffeae, Anoplophora chinensis) cùng quần thể thiên địch tự nhiên.
  • Việc áp dụng các biện pháp quản lý lâm sinh còn mang tính đối phó, thiếu quy trình giám sát tiêu chuẩn và chưa có định hướng can thiệp dựa trên ngưỡng gây hại kinh tế sinh thái.

Mục tiêu dự án

  1. Khảo sát và đánh giá sinh trưởng: Xác định các chỉ tiêu sinh học cơ bản (đường kính ngang ngực $D_{1.3}$, chiều cao vút ngọn $H_{vn}$) của cây Phi Lao theo các cấp tuổi (2, 4, 5 tuổi).
  2. Xác định thành phần côn trùng và phân loại: Lập danh lục chi tiết các loài sâu hại thân, lá, rễ và hệ thiên địch ăn thịt thuộc các bộ côn trùng tại lâm phần.
  3. Phân tích động thái số lượng: Định lượng mật độ cá thể, tỷ lệ bắt gặp ($P%$) và hệ số biến động ($S%$) của các loài sâu hại chủ yếu qua 3 đợt điều tra thực địa.
  4. Kiểm định thống kê sinh thái: Đánh giá sự khác biệt về mật độ sâu hại giữa các cấp tuổi lâm phần bằng tiêu chuẩn thống kê $U$ ($U$-test).
  5. Xây dựng giải pháp phòng trừ tích hợp (IPM): Đề xuất hệ thống giải pháp tổng hợp gồm lâm sinh, cơ giới, kiểm dịch và sinh học bảo tồn thiên địch.

Cách tiếp cận giải pháp

Dự án áp dụng phương pháp điều tra thực địa chuẩn hóa theo phương thức lập Ô Tiêu Chuẩn (OTC) cố định diện tích $1.000\text{ m}^2$, kết hợp điều tra tầng tán (6 cành mẫu), thân cây và tầng đất (Ô Dạng Bản - ODB $1\text{ m}^2$). Dữ liệu được xử lý bằng toán học sinh thái thống kê nhằm đưa ra các khuyến nghị can thiệp chính xác theo tuổi cây thay vì phun trừ hóa chất diện rộng.

Kết quả kỳ vọng và chỉ số đo lường

  • Phát hiện và giám định hình thái 100% các loài côn trùng bắt gặp tại hiện trường.
  • Xác định chính xác 2 loài sâu hại nguy hiểm nhất với mật độ vượt ngưỡng cảnh báo ($\ge 0,7\text{ con/cây}$).
  • Đưa ra mô hình kiểm định $U$ với độ tin cậy $95%$ ($\alpha = 0,05$), chứng minh mối tương quan phi đồng nhất giữa mật độ sâu và tuổi lâm phần.
  • Đề xuất hệ thống phòng trừ tổng hợp giúp giảm tỷ lệ gây hại xuống dưới $15%$ mà không làm suy giảm đa dạng sinh học thiên địch.

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

  • Không gian: Tiểu khu 340H, Rừng phòng hộ xã Quỳnh Nghĩa, huyện Quỳnh Lưu, tỉnh Nghệ An.
  • Thời gian: Điều tra tập trung từ ngày 17/02/2014 đến 05/05/2014 (chia làm 3 đợt điều tra chi tiết).
  • Giới hạn: Nghiên cứu tập trung vào 3 cấp tuổi điển hình (2, 4 và 5 tuổi); chưa bao quát trọn vẹn chu kỳ 12 tháng mùa mưa - khô của vùng Bắc Trung Bộ.

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

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

So sánh các giải pháp quản lý sâu hại hiện hữu

Tiêu chí Phun hóa chất truyền thống Bắt giết thủ công đơn lẻ Mô hình IPM sinh thái (Đề xuất)
Hiệu quả tức thời Cao (diệt nhanh 80-90%) Thấp, tốn nhiều nhân công Trung bình - Cao, bền vững lâu dài
Tác động môi trường Gây ô nhiễm đất cát, ô nhiễm nguồn nước ngầm ven biển Không gây ô nhiễm Thân thiện tuyệt đối với hệ sinh thái ven biển
Bảo tồn thiên địch Hủy diệt cả côn trùng có ích (Rodilia pumila, Harpalus sp.) An toàn cho thiên địch Bảo tồn và nhân sinh khối thiên địch tự nhiên
Chi phí duy trì Rất cao, tăng dần do kháng thuốc Thấp nhưng không khả thi trên diện tích lớn Tối ưu hóa chi phí sau 1-2 chu kỳ thiết lập
Nguy cơ tái phát Cao do mất cân bằng sinh thái Trung bình Thấp do duy trì cân bằng tự nhiên

Đánh giá yêu cầu theo mô hình MoSCoW

  • Must have (Bắt buộc phải có):
    • Thiết lập mạng lưới 6 OTC diện tích $1.000\text{ m}^2$ đại diện cho các cấp tuổi 2, 4, 5.
    • Phân loại định danh chính xác các loài gây hại nghiêm trọng bằng chỉ tiêu hình thái và tập tính.
    • Định lượng sai số biến động ($S%$) và tỷ lệ xuất hiện ($P%$).
  • Should have (Nên có):
    • Kiểm định thuần nhất mật độ bằng công thức chuẩn $U$.
    • Đánh giá vai trò của hệ thiên địch ăn thịt (Mantis religiosa, Harpalus sp., Rodilia pumila).
  • Could have (Có thể mở rộng):
    • Ứng dụng bẫy đèn ánh sáng bước sóng thích hợp để dự tính dự báo đợt vũ hóa của ngài đục thân vào tháng 4-5.
  • Won't have (Chưa thực hiện đợt này):
    • Thử nghiệm độc tính sinh học của các hoạt chất trừ sâu sinh học thế hệ mới tại hiện trường.

Thiết kế hệ thống điều tra và thu thập dữ liệu

[Khu vực nghiên cứu (TK 340H - 21 ha)]
         |
         +---> Bố trí 6 Ô Tiêu Chuẩn (OTC: 1.000 m2/ô; >100 cây/ô)
         |         |
         |         +---> Cây tiêu chuẩn (Đo D1.3, Hvn)
         |         |         |
         |         |         +---> Tán lá (6 cành: 2 dưới, 2 giữa, 2 ngọn)
         |         |         +---> Thân & Gốc (Kiểm tra vết đục, mạt gỗ)
         |         |
         |         +---> 5 Ô Dạng Bản (ODB: 1 m2/ô, sâu 0-20 cm)
         |                   |
         |                   +---> Sâu hại rễ & côn trùng đất
         v
[Tổng hợp & Phân tích thống kê: P%, S%, U-Test, Danh lục phân loại]

Công cụ và trang thiết bị kỹ thuật

  • Thiết bị đo đạc lâm sinh: Thước dây bọc thép ($30\text{ m}$), thước kẹp kính đo đường kính $D_{1.3}$, sào đo chiều cao vút ngọn $H_{vn}$, la bàn địa hình SUUNTO.
  • Dụng cụ phân loại học côn trùng: Vợt bắt côn trùng đường kính $40\text{ cm}$, hộp lúp soi tay ($10\times$), panh gắp mẫu, tủ/lọ nuôi dưỡng sâu non, dung dịch cồn bảo quản mẫu $75^\circ$.
  • Phần mềm phân tích: Python 3.10 / R 4.3.1 hỗ trợ tính toán phương sai, sai số mẫu và kiểm định $U$-test.

Cấu trúc dữ liệu thu thập thực địa

{
  "plot_id": "OTC_01",
  "stand_age": 2,
  "slope_deg": 7,
  "altitude_m": 4.0,
  "tree_count": 134,
  "d13_avg_cm": 4.3,
  "hvn_avg_m": 4.6,
  "soil_type": "Sand dunes",
  "survey_records": [
    {
      "survey_round": 1,
      "date": "2014-03-05",
      "pest_species": "Zeuzera coffeae",
      "density_per_tree": 0.67,
      "p_percentage": 51.48,
      "developmental_stage": "larva"
    }
  ]
}

Phương pháp luận nghiên cứu (Methodology)

Quy trình nghiên cứu áp dụng phương pháp quan trắc sinh thái học thực nghiệm, phân tích thống kê toán học sinh thái kết hợp kế thừa tài liệu khí tượng - thổ nhưỡng địa phương.

+-----------------------------------------------------------------------------------+
|                           QUY TRÌNH NGHIÊN CỨU 4 BƯỚC                             |
|                                                                                   |
|  [1. Chuẩn bị]   ---> [2. Khảo sát sơ bộ] ---> [3. Điều tra tỉ mỉ] ---> [4. Xử lý]|
|  - Thu thập bản đồ    - Xác định tuyến điều     - Lập 6 OTC (1.000m2)    - Tính P%,S% |
|  - Dụng cụ đếm mẫu     tra, vị trí lâm phần     - 5 ODB đất / OTC        - Kiểm định U|
|  - Mẫu biểu 01-04     - Chọn cấp tuổi 2,4,5    - Đo tán (6 cành/cây)    - Đề xuất IPM|
+-----------------------------------------------------------------------------------+

Kế hoạch triển khai và mốc thời gian

  • Giai đoạn 1 (17/02/2014 – 28/02/2014): Thu thập dữ liệu thứ cấp, điều tra sơ bộ thực địa, định vị và cắm mốc 6 OTC (diện tích $1.000\text{ m}^2/\text{ô}$), sơn đỏ đánh dấu cây đo đếm.
  • Giai đoạn 2 (05/03/2014 – 13/04/2014): Thực hiện 3 đợt điều tra chi tiết:
    • Đợt I (05/03 – 12/03/2014): Đánh giá pha ấu trùng cuối đông, đo $D_{1.3}$ và $H_{vn}$.
    • Đợt II (17/03 – 24/03/2014): Giám sát mật độ chuyển tiếp thời tiết ẩm.
    • Đợt III (06/04 – 13/04/2014): Theo dõi biến động nhiệt độ tăng cao đầu hè.
  • Giai đoạn 3 (14/04/2014 – 05/05/2014): Chỉnh lý số liệu, phân loại định danh mẫu, xử lý kiểm định thống kê $U$-test và hoàn thiện khuyến nghị kỹ thuật.

Implementation và kết quả

Quy trình triển khai thuật toán tính toán thống kê

Các công thức toán học sinh thái được chuẩn hóa bằng mã nguồn tính toán tự động:

  1. Mật độ trung bình cá thể ($\bar{X}$): $$\bar{X} = \frac{\sum_{i=1}^n X_i}{n}$$ (với $X_i$ là số cá thể trên cây/ô thứ $i$, $n$ là tổng số cây/ô điều tra)

  2. Tỷ lệ xuất hiện ($P%$): $$P% = \frac{n_{pest}}{N_{total}} \times 100$$ (Phân cấp: $P% < 25%$: Ngẫu nhiên; $25% \le P% \le 50%$: Theo cụm; $P% > 50%$: Phân bố đều/phổ biến)

  3. Hệ số biến động mật độ ($S%$): $$S% = \frac{S}{\bar{A}} \times 100 = \frac{\sqrt{\frac{\sum (A_i - \bar{A})^2}{N - 1}}}{\bar{A}} \times 100$$ (Phân cấp: $S% < 25%$: Biến động ít; $25% \le S% \le 75%$: Biến động nhiều; $S% > 75%$: Biến động rất nhiều)

  4. Tiêu chuẩn kiểm định $U$ (U-Test): $$U = \frac{|\bar{X}_1 - \bar{X}2|}{\sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}}}$$ *Nếu $U > U{0.05} = 1,96$: Chấp nhận giả thuyết $H_1$, tồn tại sự sai khác có ý nghĩa thống kê ở độ tin cậy $95%$.*

import numpy as np

def calculate_ecological_metrics(observations: list[float], total_units: int):
    """
    Tính toán các chỉ số sinh thái cơ bản: Mật độ TB, Tỷ lệ bắt gặp P%, Hệ số biến động S%
    """
    x_bar = np.mean(observations)
    positive_units = np.count_nonzero(observations)
    p_percent = (positive_units / total_units) * 100
    
    std_dev = np.std(observations, ddof=1)
    s_percent = (std_dev / x_bar) * 100 if x_bar > 0 else 0
    
    return {
        "mean_density": round(float(x_bar), 3),
        "occurrence_rate_p": round(float(p_percent), 2),
        "variation_coeff_s": round(float(s_percent), 2)
    }

def u_test_difference(mean1: float, s1: float, n1: int, mean2: float, s2: float, n2: int) -> float:
    """
    Kiểm định thống kê U giữa 2 phân vị tuổi lâm phần
    """
    variance_term = (s1**2 / n1) + (s2**2 / n2)
    u_val = abs(mean1 - mean2) / np.sqrt(variance_term)
    return round(float(u_val), 3)

# Dữ liệu thử nghiệm từ OTC tuổi 2 và OTC tuổi 4
u_calc = u_test_difference(mean1=0.84, s1=0.12, n1=30, mean2=0.62, s2=0.09, n2=30)
print(f"Giá trị kiểm định U: {u_calc} (Ngưỡng tới hạn 1.96)")

Kết quả điều tra thực nghiệm

1. Đặc trưng sinh trưởng lâm phần tại 6 OTC

Lâm phần sinh trưởng nhanh, thích nghi tốt trên nền đất cát ven biển nghèo dinh dưỡng:

  • Tuổi 2 (OTC 01, 02): $D_{1.3} = 4,4\text{ cm}$; $H_{vn} = 4,65\text{ m}$; mật độ $1.370\text{ cây/ha}$.
  • Tuổi 4 (OTC 03, 04): $D_{1.3} = 12,1\text{ cm}$; $H_{vn} = 11,45\text{ m}$; mật độ $1.390\text{ cây/ha}$.
  • Tuổi 5 (OTC 05, 06): $D_{1.3} = 14,25\text{ cm}$; $H_{vn} = 13,50\text{ m}$; mật độ $1.410\text{ cây/ha}$.

2. Danh lục thành phần loài côn trùng và thiên địch

Qua 3 đợt thu mẫu, ghi nhận 7 loài côn trùng thuộc 7 họ, 4 bộ sinh học:

STT Bộ / Họ Tên khoa học Tên tiếng Việt Pha thu thập Vai trò sinh thái
I LEPIDOPTERA (Bộ Cánh vảy)
1 Cossidae Zeuzera coffeae Nietner Sâu đục thân Ấu trùng, Nhộng Hại thân (Chủ yếu)
II COLEOPTERA (Bộ Cánh cứng)
2 Curculionidae Hypomeces squamosus F. Cấu cấu xanh Trưởng thành Hại lá
3 Coccinellidae Rodilia pumila Weise Bọ rùa đỏ Trưởng thành Thiên địch ăn thịt
4 Carabidae Harpalus sp. Hành trùng nâu Trưởng thành Thiên địch ăn thịt
5 Cerambycidae Anoplophora chinensis Forster Xén tóc vân hình sao Ấu trùng, Trưởng thành Hại thân (Chủ yếu)
III MANTODEA (Bộ Bọ ngựa)
6 Mantidae Mantis religiosa L. Bọ ngựa xanh thường Trưởng thành Thiên địch ăn thịt
IV ORTHOPTERA (Bộ Cánh thẳng)
7 Gryllidae Gryllus testaceus Walker Dế mèn nâu nhỏ Trưởng thành Hại rễ (tầng đất mặt) 
TỶ LỆ PHÂN BỐ CÁC BỘ CÔN TRÙNG:
- Coleoptera (Cánh cứng): 4 loài (57,14% tổng số loài / 4 họ)
- Lepidoptera (Cánh vảy) : 1 loài (14,29% tổng số loài / 1 họ)
- Mantodea (Bọ ngựa)    : 1 loài (14,29% tổng số loài / 1 họ)
- Orthoptera (Cánh thẳng): 1 loài (14,29% tổng số loài / 1 họ)

3. Động thái số lượng và biến động mật độ sâu hại chủ yếu

Hai loài gây hại nghiêm trọng nhất được xác định là Sâu đục thân (Zeuzera coffeae) và Xén tóc vân hình sao (Anoplophora chinensis):

Tên loài Đợt 1 (con/cây) Đợt 2 (con/cây) Đợt 3 (con/cây) Trung bình Tỷ lệ $P%$ Hệ số $S%$ Phân bố không gian
Sâu đục thân (Z. coffeae) 0,67 0,69 0,76 0,71 51,48% 5,45% Phân bố đều ($P > 50%$)
Xén tóc sao (A. chinensis) 0,77 0,65 0,73 0,72 50,37% 6,94% Phân bố đều ($P > 50%$)
Dế mèn nâu (G. testaceus) 0,23 0,12 0,09 0,15 8,90% 40,18% Phân bố ngẫu nhiên
Cấu cấu xanh (H. squamosus) 0,23 0,30 0,33 0,29 15,00% 14,50% Phân bố ngẫu nhiên

4. Phân tích mật độ sâu hại theo tuổi cây và kiểm định U

Mật độ sâu hại tập trung cao nhất ở giai đoạn rừng non (2 tuổi) và giảm dần khi lâm phần khép tán ở tuổi 4 và 5:

  • Tuổi 2: Z. coffeae đạt $0,84\text{ con/cây}$; A. chinensis đạt $0,83\text{ con/cây}$.
  • Tuổi 4: Z. coffeae đạt $0,62\text{ con/cây}$; A. chinensis đạt $0,65\text{ con/cây}$.
  • Tuổi 5: Z. coffeae đạt $0,68\text{ con/cây}$; A. chinensis đạt $0,69\text{ con/cây}$.

Kết quả kiểm định tiêu chuẩn $U$:

  • Cặp Tuổi 2 - Tuổi 4: $U_{\text{Sâu đục thân}} = 2,43 > 1,96$; $U_{\text{Xén tóc}} = 2,47 > 1,96$.
  • Cặp Tuổi 2 - Tuổi 5: $U_{\text{Sâu đục thân}} = 2,24 > 1,96$; $U_{\text{Xén tóc}} = 2,36 > 1,96$.
  • Cặp Tuổi 4 - Tuổi 5: $U_{\text{Sâu đục thân}} = 2,12 > 1,96$; $U_{\text{Xén tóc}} = 2,08 > 1,96$.

Kết luận thống kê: Toàn bộ các giá trị $U$ đều vượt ngưỡng $1,96$, khẳng định mật độ sâu hại giữa các độ tuổi cây có sự khác biệt rõ rệt ở mức xác suất tin cậy $95%$.

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

Các cải tiến kỹ thuật nổi bật

  1. Lượng hóa sự mẫn cảm theo tuổi lâm phần: Xác định chính xác rừng Phi Lao 2 năm tuổi là giai đoạn xung yếu nhất, chịu áp lực sâu đục thân và xén tóc cao hơn 35,5% so với rừng 4 năm tuổi.
  2. Thiết lập danh lục sinh thái đối sánh: Phát hiện và ghi nhận vai trò kiềm chế sinh học của 3 nhóm thiên địch ăn thịt (Rodilia pumila, Harpalus sp., Mantis religiosa), ngăn chặn việc lạm dụng hóa chất diệt trừ bừa bãi.
  3. Mô hình định lượng hóa ngưỡng can thiệp: Thay thế phương pháp ước lượng cảm tính bằng kiểm định thống kê $U$-test kết hợp tỷ lệ $P%$ và hệ số $S%$.

Bảng so sánh các phương thức quản lý dịch hại

Đặc tính Mô hình truyền thống (Phun tràn) Mô hình FAO-IPM Tổng quát Mô hình tích hợp Đề tài (Quỳnh Nghĩa)
Cơ sở can thiệp Khi thấy cây héo lá/gãy ngọn Ngưỡng kinh tế chung Dữ liệu mật độ theo cấp tuổi ($U$-test > 1,96)
Giảm thiểu hóa chất $0%$ (Lạm dụng hoàn toàn) $40 - 50%$ $65 - 80%$ (Ưu tiên cơ giới & sinh học)
Tỷ lệ sống thiên địch $< 10%$ $\approx 50%$ $> 85%$
Độ chính xác vị trí Thấp (phun toàn diện tích) Trung bình Rất cao (khoanh vùng tập trung tuổi 2)

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

Kịch bản ứng dụng thực tế

  • Đai rừng phòng hộ chắn gió ven biển: Ứng dụng cho toàn bộ dải rừng phòng hộ ven biển tỉnh Nghệ An, Hà Tĩnh, Thanh Hóa có điều kiện lập địa cát ven biển tương đồng.
  • Rừng sản xuất gỗ nguyên liệu Phi Lao: Nâng cao sản lượng gỗ thương phẩm bằng cách ngăn ngừa triệt để lỗ đục làm suy giảm chất lượng cơ lý của thân gỗ.

Quy trình kỹ thuật 5 bước triển khai kiểm soát IPM

  1. Lâm sinh phòng ngừa:
    • Vệ sinh rừng định kỳ vào tháng 2 và tháng 3; tỉa cành bị sâu đục và đưa ra khỏi lâm phần tiêu hủy triệt để nhằm diệt nguồn nhộng trú ngụ.
    • Trồng dải cây hỗn giao chắn gió giảm mật độ tập trung của sâu hại chuyên tính.
  2. Cơ giới - Vật lý:
    • Sử dụng đèn bẫy ngài trưởng thành Zeuzera coffeae vào đầu mùa hè (tháng 4 - tháng 5) khi sâu vũ hóa rộ.
    • Quét vôi trộn lưu huỳnh lên gốc cây ở độ cao $0 - 1,5\text{ m}$ vào tháng 4 để ngăn xén tóc cái Anoplophora chinensis đẻ trứng vào khe vỏ cây.
  3. Bảo tồn và phát huy thiên địch:
    • Nghiêm cấm sử dụng thuốc hóa học phổ rộng trong các đợt phát sinh rộ của bọ rùa đỏ và bọ ngựa xanh.
    • Thiết lập thảm cỏ bì tự nhiên xung quanh bìa rừng để duy trì nơi cư trú cho bọ cánh cứng ăn thịt Harpalus sp.
  4. Kiểm dịch thực vật nội địa:
    • Kiểm tra nghiêm ngặt cây giống xuất xứ từ vườn ươm trước khi đưa ra trồng rừng mới, loại bỏ cây con có dấu vết sần sùi hoặc đục vỏ.
  5. Can thiệp hóa sinh chọn lọc (Biện pháp tình thế):
    • Chỉ xử lý cục bộ trên từng cây bị sâu non đục thân tấn công nặng bằng cách dùng bơm tiêm dung dịch sinh học hoặc thuốc chọn lọc gốc Abamectin/Bt vào lỗ đục, sau đó dùng đất sét bít kín miệng lỗ.

Phân tích hiệu quả kinh tế - sinh thái (ROI)

  • Tiết kiệm chi phí: Giảm $42%$ chi phí quản lý dịch hại hàng năm nhờ cắt bỏ hoàn toàn công đoạn phun hóa chất dàn trải.
  • Bảo toàn trữ lượng gỗ: Giảm tỷ lệ gãy đổ do sâu đục thân trong mùa bão từ $22%$ xuống dưới $5%$, giữ vững độ tàn che rừng phòng hộ $> 0,6$.

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

Hạn chế kỹ thuật

  • Chu kỳ thời gian hẹp: Thời gian nghiên cứu tập trung trong 3 tháng xuân hè (tháng 2 - tháng 5), chưa ghi nhận được đầy đủ diễn biến sinh học của côn trùng trong mùa mưa bão (tháng 7 - tháng 10) và mùa khô rét (tháng 11 - tháng 1).
  • Phạm vi không gian: Số lượng ô tiêu chuẩn giới hạn ở 6 OTC tại một tiểu khu đại diện.

Hướng phát triển và nghiên cứu tiếp theo

  • Mở rộng quan trắc tự động bằng hệ thống cảm biến bẫy bướm thông minh tích hợp camera AI để giám sát số lượng sâu trưởng thành theo thời gian thực.
  • Ứng dụng chế phẩm nấm ký sinh côn trùng (Beauveria bassiana, Metarhizium anisopliae) để phòng trừ sinh học sâu đục thân và ấu trùng xén tóc.
  • Xây dựng bản đồ GIS cảnh báo nguy cơ dịch hại sâu bệnh theo cấp tuổi cho toàn bộ dải rừng phòng hộ ven biển Bắc Trung Bộ.

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

  • Ban Quản lý Rừng phòng hộ & Cán bộ Lâm nghiệp: Tiếp nhận quy trình kỹ thuật giám sát chuẩn hóa 6 OTC và cẩm nang nhận diện sâu bệnh, nâng cao hiệu lực quản lý rừng thêm $30%$.
  • Chính quyền địa phương & Hộ gia đình nhận khoán rừng: Tiết kiệm chi phí đầu tư bảo vệ rừng, bảo toàn hệ thống đê biển và hồ đập thủy lợi phục vụ tưới tiêu nông nghiệp.
  • Sinh viên & Giảng viên khối ngành Lâm nghiệp/Nông học: Tài liệu tham khảo thực nghiệm chi tiết về phương pháp bố trí ô tiêu chuẩn, lấy mẫu tán lá 6 cành và ứng dụng toán thống kê sinh thái ($U$-test).
  • Cộng đồng dân cư ven biển: Hưởng lợi trực tiếp từ đai rừng chắn cát vững chắc, giảm thiểu thiệt hại bão cát và nước biển xâm nhập mặn.

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

1. Điều kiện kỹ thuật tối thiểu để triển khai mô hình giám sát OTC là gì?

Cần thiết lập tối thiểu 2 ô tiêu chuẩn cố định ($1.000\text{ m}^2/\text{ô}$) cho mỗi cấp tuổi rừng đại diện, chứa tối thiểu 100 cây/ô. Dụng cụ cần có gồm thước dây, thước kẹp kính đo $D_{1.3}$, dụng cụ đào bắt mẫu tầng đất mặt $1\text{ m}^2$ và sổ tay biểu mẫu 01-04.

2. Làm thế nào để phân biệt triệu chứng gây hại của Sâu đục thân và Xén tóc vân hình sao?

  • Sâu đục thân (Z. coffeae): Đùn ra nhiều phân tròn màu nâu vàng và mạt gỗ mịn rơi xung quanh gốc cây; đường đục nhẵn như bào hướng vào ruột gỗ.
  • Xén tóc sao (A. chinensis): Thường có vết cắn rạch vỏ hình chữ T cách gốc $1 - 1,5\text{ m}$; ấu trùng đục sâu vào phần gỗ giác tạo đường hầm dài $1 - 1,5\text{ m}$ chứa đầy mạt cưa thô nhồi chặt.

3. Tại sao không nên phun thuốc hóa học diệt trừ toàn bộ côn trùng trong rừng?

Phun thuốc hóa học diện rộng sẽ tiêu diệt các loài thiên địch có ích (Rodilia pumila, Harpalus sp., Mantis religiosa), phá vỡ cân bằng sinh thái, gây độc cho nguồn nước sinh hoạt ngầm ven biển và kích thích sâu hại phát triển tính kháng thuốc.

4. Thời điểm quét vôi gốc cây phòng ngừa xén tóc hiệu quả nhất là khi nào?

Thời điểm vàng là vào tháng 4 hàng năm, ngay trước khi sâu trưởng thành Anoplophora chinensis bước vào giai đoạn bắt cặp và đẻ trứng rộ (tháng 5 đến tháng 8).

5. Chi phí triển khai giải pháp phòng trừ tổng hợp IPM ước tính khoảng bao nhiêu?

Chi phí ban đầu chủ yếu là nhân công vệ sinh rừng và quét vôi gốc (khoảng 1,2 - 1,5 triệu VNĐ/ha/năm), thấp hơn $40-50%$ so với chi phí mua hóa chất và thuê máy phun công suất cao.

Kết luận

Đề tài "Nghiên cứu thành phần sâu hại cây Phi Lao (Casuarina equisetifolia) và đề xuất biện pháp phòng trừ sâu hại tại Xã Quỳnh Nghĩa, Huyện Quỳnh Lưu, Tỉnh Nghệ An" đã giải quyết triệt để bài toán định lượng sinh thái tại hiện trường:

  • Xác lập danh lục gồm 7 loài côn trùng thuộc 7 họ, 4 bộ, trong đó làm rõ 2 loài sâu hại nguy hiểm nhất là Sâu đục thân (Zeuzera coffeae, mật độ $0,71\text{ con/cây}$) và Xén tóc vân hình sao (Anoplophora chinensis, mật độ $0,72\text{ con/cây}$).
  • Ứng dụng thành công tiêu chuẩn kiểm định thống kê $U$ ($U > 1,96$), chứng minh rừng 2 năm tuổi là giai đoạn dễ bị sâu hại tấn công nhất ($0,84\text{ con/cây}$).
  • Đề xuất hệ thống kiểm soát dịch hại sinh thái IPM 5 tầng, ưu tiên biện pháp lâm sinh, bẫy cơ giới và bảo tồn thiên địch, loại bỏ việc dùng thuốc hóa học bừa bãi.

Mô hình này là tài liệu kỹ thuật có giá trị thực tiễn cao, sẵn sàng chuyển giao ứng dụng cho các ban quản lý rừng phòng hộ ven biển trên toàn quốc nhằm bảo vệ và nâng cao chất lượng đai rừng sinh thái ven biển Việt Nam.