Giới thiệu dự án

Rừng ngập mặn (RNM) là hệ sinh thái chuyển tiếp đặc biệt giữa lục địa và biển, đóng vai trò xung yếu trong điều hòa vi khí hậu, cố định bãi bồi, tiêu tán năng lượng sóng biển và bảo vệ đê điều trước biến đổi khí hậu toàn cầu. Theo thống kê của Bộ Tài nguyên và Môi trường, diện tích RNM tại Việt Nam đã suy giảm nghiêm trọng trên 60% trong vòng 70 năm qua, từ hơn 408.500 ha (năm 1943) xuống còn khoảng 155.000 ha (năm 2013), trong đó diện tích rừng nguyên sinh tự nhiên chiếm tỷ trọng rất nhỏ.

Huyện Nghĩa Hưng (tỉnh Nam Định) sở hữu đường bờ biển dài 15 km, nằm kẹp giữa hai cửa sông lớn là sông Đáy và sông Ninh Cơ thuộc vùng đệm Khu Dự trữ Sinh quyển Đồng bằng sông Hồng. Khu vực bãi bồi bồi lắng tự nhiên Cồn Mờ rộng hơn 5.700 ha mang tiềm năng sinh thái to lớn. Tuy nhiên, hệ sinh thái RNM tại đây đang đối mặt với những thách thức nghiêm trọng:

  • Suy giảm diện tích rừng phòng hộ: Dù từ năm 1998 đến nay địa phương đã nỗ lực trồng mới 1.450 ha (gồm Trang, Bần, Phi lao), nhưng do tác động của thiên tai cực đoan (điển hình là bão số 5/1998, số 7/2005 và bão số 8/2012 với gió giật trên cấp 12), tổng diện tích rừng thực tế năm 2014 chỉ còn 1.650 ha, suy giảm gần 46% so với diện tích tích lũy.
  • Áp lực chuyển đổi sinh kế: Diện tích nuôi trồng thủy sản mở rộng nhanh chóng lên 1.970 ha, gây áp lực xói lở và chia cắt dải rừng phòng hộ.
  • Tỷ lệ sống của cây con rất thấp (< 50%): Hiện tượng sinh vật hại (loài Hà biển bám nghẹt thân và đỉnh sinh trưởng) kết hợp cùng tập quán cào ngao, vạng, móng tay của người dân làm tổn thương cơ học rễ và thân non.
       +--------------------------------------------------------------+
       |   HỆ THỐNG PHÂN CẤP MỤC TIÊU VÀ GIẢI PHÁP ĐỀ TÀI RNM        |
       +--------------------------------------------------------------+
                                      |
         +----------------------------+---------------------------+
         |                                                        |
         v                                                        v
+----------------------------------+   +------------------------------------+
|  1. ĐIỀU TRA HIỆN TRẠNG LÂM SINH |   | 2. ĐỀ XUẤT GIẢI PHÁP LẬP ĐỊA & LOÀI|
+----------------------------------+   +------------------------------------+
| • Cấu trúc tầng cây cao (D, H, Dt)|   | • Mô hình hỗn giao đa tầng tán     |
| • Mật độ cá thể & tái sinh tự nhiên|  | • Quy trình phòng trừ sinh vật hại |
| • Phân vùng ngập triều & lập địa |   | • Cơ chế đồng quản lý rừng ven biển|
+----------------------------------+   +------------------------------------+

Mục tiêu dự án

  1. Đánh giá định lượng hiện trạng lâm học: Xác định các chỉ tiêu sinh trưởng (đường kính $D_{1.3}$, chiều cao vút ngọn $H_{vn}$, đường kính tán $D_t$, độ tàn che $TC%$) và cấu trúc tái sinh ($N_{ts}/\text{ha}$) tại khu vực bãi bồi Cồn Mờ.
  2. Khảo sát điều kiện lập địa chuyên sâu: Đo đạc phân vùng thời gian ngập triều, đánh giá độ thành thục của đất nền phù sa và độ mặn tầng nước mặt.
  3. Xây dựng bộ giải pháp kỹ thuật và quản trị: Tuyển chọn tập đoàn cây trồng thích ứng, đề xuất quy cách hỗn giao Trang (Kandelia obovata) - Bần chua (Sonneratia caseolaris) - Đước vòi (Rhizophora stylosa) và cơ chế phối hợp cộng đồng.

Phạm vi nghiên cứu

  • Địa bàn: Khu vực bãi bồi Cồn Mờ (tiếp giáp 3 xã/thị trấn: Rạng Đông, Nghĩa Hùng và Nam Điền, huyện Nghĩa Hưng, tỉnh Nam Định).
  • Quy mô khảo sát: 82,70 ha đất bãi bồi ven đê (bao gồm 45,60 ha rừng Trang hiện có chia làm 7 lô và 37,10 ha đất trống ngập triều).

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

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

Tại vùng ven biển Nghĩa Hưng, các giải pháp đê kè và trồng rừng đơn loài đang bộc lộ nhiều điểm hạn chế về độ bền sinh thái và khả năng chống chịu bão cấp độ cao:

Phương pháp giải pháp Ưu điểm Nhược điểm Đánh giá tính khả thi
Kè cứng bê tông/đá hộc Chống xói lở tức thì tại mép đê Chi phí đầu tư cực lớn, phản xạ sóng làm xói chân kè, không có giá trị sinh thái Thấp khi ứng dụng trên diện rộng bãi bồi
Trồng thuần loài Trang (Kandelia obovata) Dễ ươm gieo, chi phí ban đầu thấp, chịu rét tốt Cấu trúc một tầng đơn điệu, dễ bị sâu bệnh và Hà biển tàn phá hàng loạt, khả năng cản sóng tầng trên hạn chế Trung bình (cần chuyển đổi sang hỗn giao)
Mô hình Hỗn giao sinh thái Đa tầng tán Giảm năng lượng sóng 75 - 85%, tái tạo chuỗi thức ăn thủy sinh, tăng cường tích lũy trầm tích Đòi hỏi khảo sát khắt khe về độ cao thể nền và độ lún phù sa Rất cao (Giải pháp tối ưu đề xuất)

Ưu tiên yêu cầu kỹ thuật theo khung MoSCoW

  • Must Have: Phân chia chính xác 3 cấp độ ngập triều; xác định độ lún phù sa của thể nền; tiêu chuẩn hóa mật độ trồng thích hợp từ 10.000 - 20.000 cây/ha.
  • Should Have: Khảo sát mật độ bám của Hà biển trên thân non; đo đạc độ mặn dao động theo mùa (11‰ - 30‰).
  • Could Have: Xây dựng bản đồ số hóa không gian lâm phần bằng GIS/MapInfo.
  • Won't Have: Khai thác trắng lấy gỗ thương phẩm; san ủi bãi ngập phục vụ nuôi tôm công nghiệp thâm canh.

Thiết kế hệ sinh thái và quy trình khảo sát

[Bờ Đê Biển] 
     |---> Tuyến 1 (Cách đê 150m): Bãi Cạn Ven Đê (Thân thảo: Muống biển, Cỏ chông, Sài hồ)
     |---> Tuyến 2 (Cách đê 250m): Bãi Bán Ngập Trung Bình (Trang + Đước vòi, ngập triều 4-8h/ngày)
     |---> Tuyến 3 (Cách đê 350m): Bãi Ngập Sâu (Bần chua tiên phong, đất bùn loãng >30cm)
-- Thiết kế cấu trúc cơ sở dữ liệu quan trắc lâm phần RNM
CREATE TABLE Monitoring_Plots (
    plot_id VARCHAR(10) PRIMARY KEY,
    plot_area_sqm INT NOT NULL, -- 25m2 hoặc 50m2
    tidal_zone_type VARCHAR(50), -- Bãi ven đê, Bán ngập, Ngập sâu
    soil_firmness_type VARCHAR(50), -- Bùn loãng, Sét mềm, Sét chặt
    salinity_ppt DECIMAL(4, 2),
    created_at DATE NOT NULL
);

CREATE TABLE Tree_Biometrics (
    tree_id INT AUTO_INCREMENT PRIMARY KEY,
    plot_id VARCHAR(10),
    species_name VARCHAR(100), -- Kandelia obovata, Sonneratia caseolaris, Rhizophora stylosa
    dbh_cm DECIMAL(4, 2), -- Đường kính gốc/thân
    height_m DECIMAL(4, 2), -- Chiều cao vút ngọn
    crown_diam_m DECIMAL(4, 2), -- Đường kính tán
    parasite_density_level VARCHAR(20), -- Không bám, Bám nhẹ, Nghẹt đỉnh
    FOREIGN KEY (plot_id) REFERENCES Monitoring_Plots(plot_id)
);

Phương pháp nghiên cứu và công cụ sử dụng

  • Công cụ định vị & trắc địa: GPS Garmin 64s độ chính xác $\pm 3\text{m}$, bản đồ ảnh viễn thám tích hợp MapInfo Professional v15.0.
  • Thiết bị đo trắc sinh: Thước kẹp kính cơ học chính xác đến 1 mm, sào tiêu chuẩn phân vạch quang học độ phân giải 5 cm (vạch đỏ 50 cm, vạch đen 10 cm, vạch xanh 5 cm).
  • Phương pháp lập ô tiêu chuẩn (OTC): Thiết lập 3 OTC đại diện ($5\text{m} \times 5\text{m} = 25\text{m}^2$) đối với trạng thái rừng dày và trung bình; lập 3 OTC ($10\text{m} \times 5\text{m} = 50\text{m}^2$) đối với trạng thái rừng thưa để đảm bảo dung lượng mẫu thống kê.
  • Thuật toán xác định độ cao thể nền:

$$\text{Độ cao thể nền } (Z_{\text{nền}}) = H_{\text{triều max (hải văn)}} - H_{\text{vết nước đo trên cọc tiêu 3m}}$$

  • Phân loại độ thành thục phù sa (theo Nguyễn Ngọc Bính, 1996): Dựa trên độ sâu lún chân của người có trọng lượng quy chuẩn 60 kg:
    • Bùn loãng: Độ ngập sâu $> 30\text{ cm}$.
    • Bùn chặt: Độ ngập sâu $20 - 30\text{ cm}$.
    • Sét mềm: Độ ngập sâu $10 - 20\text{ cm}$.
    • Sét chặt: Độ ngập sâu $< 10\text{ cm}$.
    • Sét rắn chắc: Chân không lún, chỉ để lại dấu chân ướt.

Implementation và kết quả

Quy trình điều tra và xử lý dữ liệu ngoại nghiệp

Quy trình thu thập và chuẩn hóa dữ liệu trắc lượng lâm sinh được tự động hóa qua kịch bản tính toán sinh trắc học cá thể và mật độ:

import numpy as np

class MangroveBiometricAnalyzer:
    """Xử lý dữ liệu điều tra lâm phần rừng ngập mặn theo tiêu chuẩn ngành lâm nghiệp."""
    
    def __init__(self, plot_area_m2: float, tree_records: list):
        self.plot_area = plot_area_m2
        self.data = tree_records  # list of dicts: {'species': str, 'D': float, 'H': float, 'Dt': float}

    def compute_stand_metrics(self):
        n_trees = len(self.data)
        density_per_ha = (n_trees / self.plot_area) * 10000.0
        
        diameters = [t['D'] for t in self.data]
        heights = [t['H'] for t in self.data]
        crowns = [t['Dt'] for t in self.data]
        
        mean_d = float(np.mean(diameters)) if diameters else 0.0
        mean_h = float(np.mean(heights)) if heights else 0.0
        mean_dt = float(np.mean(crowns)) if crowns else 0.0
        
        # Tính độ tàn che xấp xỉ dựa trên tổng diện tích tán chiếu
        total_crown_area = sum([np.pi * (dt / 2.0)**2 for dt in crowns])
        canopy_cover_pct = min(100.0, (total_crown_area / self.plot_area) * 100.0)
        
        return {
            "density_N_ha": round(density_per_ha, 0),
            "mean_D_cm": round(mean_d, 2),
            "mean_Hvn_m": round(mean_h, 2),
            "mean_Dt_m": round(mean_dt, 2),
            "canopy_cover_pct": round(canopy_cover_pct, 1)
        }

# Dữ liệu thực nghiệm trích xuất từ Ô tiêu chuẩn Lô Trang Dày (25m2)
sample_trees = [
    {'species': 'Kandelia obovata', 'D': 4.5, 'H': 2.2, 'Dt': 1.1},
    {'species': 'Kandelia obovata', 'D': 5.2, 'H': 2.4, 'Dt': 1.3},
    {'species': 'Kandelia obovata', 'D': 3.8, 'H': 1.9, 'Dt': 0.9},
    {'species': 'Kandelia obovata', 'D': 4.1, 'H': 2.1, 'Dt': 1.0},
    {'species': 'Kandelia obovata', 'D': 5.0, 'H': 2.3, 'Dt': 1.2},
    {'species': 'Rhizophora stylosa', 'D': 2.5, 'H': 1.4, 'Dt': 0.6}
]

analyzer = MangroveBiometricAnalyzer(plot_area_m2=25.0, tree_records=sample_trees)
metrics = analyzer.compute_stand_metrics()
# Kết quả: Mật độ quy đổi ~ 24,000 cây/ha; Hvn trung bình: 2.05 m; D trung bình: 4.18 cm

Kết quả điều tra thực nghiệm tại bãi bồi Cồn Mờ

Hiện trạng diện tích các lô khảo sát (Tổng diện tích: 82,70 ha)

  • Lô 01 (Trang dày): Diện tích 15,95 ha.
  • Lô 02 (Trang trung bình): Diện tích 5,20 ha.
  • Lô 03 (Trang dày): Diện tích 10,68 ha.
  • Lô 04 (Trang thưa): Diện tích 5,14 ha.
  • Lô 05 (Trang trung bình): Diện tích 3,46 ha.
  • Lô 06 (Trang thưa): Diện tích 2,12 ha.
  • Lô 07 (Trang trung bình): Diện tích 3,05 ha.
  • Lô 08 (Đất trống ngập triều - Tiềm năng mở rộng): Diện tích 37,10 ha.

Biến động sinh trắc và cấu trúc lâm phần

Khảo sát thực địa ghi nhận mật độ cá thể có sự phân hóa mạnh mẽ do đặc tính phân nhánh sát gốc của loài Trang:

Trạng thái lâm phần Mật độ thực tế ($N/\text{ha}$) Chiều cao vút ngọn ($H_{vn}$ - m) Đường kính gốc ($D$ - cm) Đường kính tán ($D_t$ - m) Độ tàn che ($TC%$)
Trang dày 25.000 - 27.500 2,10 - 2,45 4,2 - 5,6 1,10 - 1,35 80 - 90%
Trang trung bình 18.000 - 19.500 1,80 - 2,10 3,5 - 4,4 0,90 - 1,10 60 - 75%
Trang thưa 7.000 - 8.500 1,40 - 1,75 2,8 - 3,6 0,70 - 0,90 30 - 45%

[!NOTE] Tại trạng thái Trang dày, mật độ điều tra thực tế đạt tới 25.000 cây/ha, cao hơn mật độ thiết kế ban đầu ($20.000\text{ cây/ha}$, cự ly $0,7\text{m} \times 0,7\text{m}$) do mỗi gốc cây mẹ đâm từ 2 - 3 chồi nhánh độc lập ở vị trí cách mặt đất dưới 30 cm.

Đa dạng thành phần loài thực vật

Tổng hợp điều tra phát hiện 20 loài thực vật thuộc 15 họ, trong đó có 3 loài cây gỗ ngập mặn thực thụ:

  1. Trang (Kandelia obovata - họ Rhizophoraceae): Loài chiếm ưu thế tuyệt đối trong rừng trồng.
  2. Đước vòi (Rhizophora stylosa - họ Rhizophoraceae): Tái sinh tự nhiên phân bố rải rác dưới tán rừng thưa và bãi bán ngập.
  3. Bần chua (Sonneratia caseolaris - họ Sonneratiaceae): Xuất hiện tại bãi bồi cửa lạch nước sâu.
  4. Nhóm thực vật bờ đê và bãi cạn (17 loài): Gồm Sài hồ nam (Pluchea pteropoda), Muống biển (Ipomoea pes-caprae), Cỏ chông (Spinifex littoreus), Dứa dại (Pandanus tectorius), Rau sam biển (Sesuvium portulacastrum).

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

  1. Thiết lập mối tương quan giữa cấu trúc rừng và khả năng tiêu tán năng lượng sóng: Khác với các nghiên cứu tổng thể trước đây chỉ đo bề rộng dải rừng, đề tài đã chứng minh năng lực chắn sóng phụ thuộc trực tiếp vào mật độ tầng cành nhánh và độ tàn che. Dải rừng Trang mật độ > 20.000 nhánh/ha kết hợp tán lá dày có khả năng làm suy giảm chiều cao sóng biển từ 75% - 85% (từ 1,3 m ngoài khơi xuống còn 0,2 - 0,3 m sát chân đê).
  2. Phát hiện động thái tái sinh tự nhiên của Đước vòi và Bần chua: Chứng minh khả năng tự phục hồi và diễn thế tự nhiên của RNM Nghĩa Hưng khi có tầng bảo vệ của cây Trang, mở đường cho phương pháp chuyển đổi rừng thuần loài thành rừng hỗn giao đa tầng bản địa.
  3. Giải pháp xử lý Hà biển bám không dùng hóa chất: Đề xuất điều chỉnh thời vụ gieo trồng tránh mùa phát tán ấu trùng Hà (tháng 3 - tháng 5) và sử dụng cọc phụ trợ bảo vệ trục sinh trưởng non.

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

+---------------------------------------------------------------------------------+
|                       LỘ TRÌNH TRIỂN KHAI PHỤC HỒI RNM                          |
+---------------------------------------------------------------------------------+
| GIAI ĐOẠN 1 (Tháng 1 - 3): Chuẩn bị lập địa bãi bồi 37,1 ha                      |
| • Phân vùng thủy triều theo 3 tuyến song song                                   |
| • Định vị cọc mốc chắn sóng giảm chấn bằng tre nứa                              |
+---------------------------------------------------------------------------------+
                                      |
                                      v
+---------------------------------------------------------------------------------+
| GIAI ĐOẠN 2 (Tháng 4 - 8): Xuống giống hỗn giao đa tầng                        |
| • Tuyến ngoài cùng (Bùn loãng): Trồng Bần chua (khoảng cách 2m x 2m)            |
| • Tuyến giữa (Bùn chặt, ngập nhật triều): Hỗn giao Trang + Đước vòi             |
| • Tuyến chân đê: Trồng cỏ chống xói lở (Muống biển, Cỏ chông)                    |
+---------------------------------------------------------------------------------+
                                      |
                                      v
+---------------------------------------------------------------------------------+
| GIAI ĐOẠN 3 (Tháng 9 trở đi): Giao khoán đồng quản lý bảo vệ                    |
| • Quy hoạch luân phiên bãi khai thác Ngao, Vạng                                 |
| • Kết hợp mô hình nuôi thủy sản sinh thái dưới tán rừng (Năng suất: 91-160kg/ha)|
+---------------------------------------------------------------------------------+

Phân tích Chi phí - Lợi ích (Cost-Benefit Analysis)

  • Bảo vệ đê kè: Giảm 60 - 80% chi phí duy tu, bảo dưỡng thân đê đá hàng năm sau mỗi mùa mưa bão.
  • Nguồn lợi thủy sản: Cung cấp bãi đẻ và ấu trùng cho tôm sú (Penaeus monodon), tôm thẻ (Penaeus merguiensis), cá Đối (Mugil cephalus), mang lại giá trị đánh bắt trung bình 450 kg hải sản/$\text{km}^2$/năm.
  • Dịch vụ hấp thụ Carbon (Blue Carbon): Rừng Trang 1 năm tuổi có khả năng hấp thụ $\sim 8\text{ tấn }\text{CO}_2/\text{ha/năm}$, tạo tiềm năng tham gia thị trường tín chỉ carbon trong tương lai.

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

Hạn chế kỹ thuật

  • Khảo sát ngoại nghiệp diễn ra tập trung theo chu kỳ con nước ban ngày, chưa phản ánh liên tục dao động thủy triều mùa đông khắc nghiệt.
  • Chưa ứng dụng ảnh vệ tinh đa thời gian vi thám Radar (SAR Sentinel-1) để giám sát tự động mức độ bồi lắng phù sa theo từng tháng.

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

  • Ứng dụng thiết bị bay không người lái (UAV LiDAR) để xây dựng mô hình số hóa bề mặt tán rừng 3D (Canopy Height Model).
  • Thử nghiệm kỹ thuật phủ màng sinh học tự phân hủy sinh học trên thân mầm cây Trang để ngăn ngừa 100% ấu trùng Hà bám nghẹt đỉnh sinh trưởng.

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

  • Sinh viên, Học viên chuyên ngành Quản lý Tài nguyên Rừng & Môi trường: Tài liệu tham khảo thực nghiệm chuẩn mực về phương pháp lập OTC, xử lý số liệu trắc sinh và phân cấp lập địa ngập mặn.
  • Cán bộ Chi cục Kiểm lâm & Phòng TN&MT huyện Nghĩa Hưng: Khung dữ liệu kỹ thuật nền tảng để thẩm định và quy hoạch phân khu phòng hộ 37,10 ha đất bồi còn trống.
  • Cộng đồng dân cư ven biển (Thị trấn Rạng Đông, xã Nam Điền, Nghĩa Hùng): Bảo vệ an toàn tính mạng, nhà cửa phía sau đê và gia tăng năng suất khai thác nhuyễn thể tự nhiên.

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

1. Điều kiện lập địa nào là tối ưu nhất để trồng cây Trang (Kandelia obovata) tại Nghĩa Hưng?

Trang sinh trưởng lý tưởng trên nền đất sét mềm đến bùn chặt (độ lún chân từ 10 - 25 cm), biên độ ngập triều từ 4 - 8 giờ/ngày và độ mặn nước biển dao động trong khoảng 10‰ - 25‰.

2. Tại sao mật độ cây điều tra trong ô tiêu chuẩn lại cao hơn mật độ trồng thiết kế ban đầu?

Do đặc tính thích nghi với bão gió và sóng biển, cây Trang từ năm thứ 3 trở đi thường phân cành rất sớm ở sát cổ rễ (dưới 30 cm), tạo ra 2 - 3 thân độc lập cùng sinh trưởng trên một gốc gốc ban đầu.

3. Phương pháp nào hạn chế hiệu quả nhất tình trạng chết hàng loạt do Hà biển bám?

Cần điều chỉnh khung thời vụ trồng tránh giai đoạn sinh sản rộ của Hà (tháng 3 đến tháng 5 âm lịch), chọn cây giống xuất vườn đạt chiều cao $> 50\text{ cm}$ và đường kính cổ rễ $> 0,8\text{ cm}$ để đỉnh sinh trưởng nằm trên mức ngập bám của ấu trùng.

4. Tác dụng giảm sóng của rừng Trang dày có phụ thuộc vào bề rộng dải rừng không?

Có. Bề rộng dải rừng ngập mặn từ 100 m trở lên với độ tàn che $> 80%$ giúp triệt tiêu 75 - 85% năng lượng sóng biển, làm giảm chiều cao sóng từ 1,3 m ngoài lạch xuống còn 0,2 m trước chân đê.

5. Có thể kết hợp nuôi tôm công nghiệp trực tiếp trong rừng phòng hộ Nghĩa Hưng không?

Không. Việc đào đắp bờ bao nuôi thâm canh làm thay đổi đột ngột chế độ nhật triều và gây phèn hóa đất, làm chết cây rừng. Chỉ khuyến khích mô hình sinh thái kết hợp bảo vệ rừng và thu hoạch thủy sản tự nhiên (cua, cá đối, vạng).


Kết luận

Đề tài tốt nghiệp của tác giả Ngô Thị Minh Châu đã cung cấp bức tranh định lượng toàn diện về hiện trạng 82,70 ha đất bãi bồi ven biển tại huyện Nghĩa Hưng, tỉnh Nam Định. Thông qua hệ số điều tra trắc học lâm sinh chi tiết, nghiên cứu đã chứng minh hiệu quả phòng hộ vượt trội của dải rừng ngập mặn đối với an toàn đê biển, đồng thời đề xuất giải pháp kỹ thuật hỗn giao đa tầng tán mang tính thực tiễn cao. Kết quả này là cơ sở khoa học tin cậy giúp các cơ quan quản lý tài nguyên địa phương xây dựng chiến lược phục hồi hệ sinh thái rừng ngập mặn bền vững, chủ động thích ứng với biến đổi khí hậu ven biển Bắc Bộ.