Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu đang gia tăng áp lực nặng nề lên các hệ sinh thái đất ngập nước ven biển, trong đó mực nước biển toàn cầu đã tăng với tốc độ khoảng 3,1 mm/năm trong giai đoạn 1993 - 2003 và có xu hướng tiếp tục gia tăng. Tại Việt Nam, nhiệt độ trung bình năm trong 50 năm qua đã tăng 0,7°C và mực nước biển dâng khoảng 20 cm. Theo các kịch bản biến đổi khí hậu của Bộ Tài nguyên và Môi trường, nếu mực nước biển dâng 1 m vào cuối thế kỷ 21, hơn 10% diện tích đồng bằng sông Hồng và 39% diện tích đồng bằng sông Cửu Long có nguy cơ bị ngập chìm, tác động trực tiếp đến đời sống của trên 9% dân số miền Bắc.

Trong bối cảnh đó, Vườn Quốc gia Xuân Thủy (tỉnh Nam Định) – khu Ramsar đầu tiên của Việt Nam và Đông Nam Á – đóng vai trò là vùng đệm sinh thái trọng yếu bảo vệ bờ biển trước thiên tai, sóng bão. Tuy nhiên, hệ sinh thái rừng ngập mặn tại đây đang đứng trước nguy cơ suy thoái nghiêm trọng do sự thay đổi bất thường của độ ngập nước triều, dòng chảy phù sa và sự thu hẹp không gian thích ứng sinh học. Vấn đề nghiên cứu trọng tâm được đặt ra là làm thế nào để định lượng và dự báo chính xác động thái phát triển, diễn thế của các loài cây ngập mặn trước tác động của chế độ ngập biến đổi.

Mục tiêu cụ thể của luận văn tập trung vào việc xác định mối quan hệ tương tác giữa đặc điểm môi trường vật lý, địa hình với thảm thực vật; phân tích các yếu tố chi phối chế độ ngập; và ứng dụng công cụ mô hình hóa toán học hướng cá thể nhằm mô phỏng diễn biến không gian, cấu trúc rừng ngập mặn theo các kịch bản nước biển dâng. Phạm vi nghiên cứu được triển khai tại vùng lõi Vườn Quốc gia Xuân Thủy thuộc tọa độ từ 20°10’ đến 20°15’ vĩ độ Bắc và 106°20’ đến 106°32’ kinh độ Đông, trải dài trên địa bàn 5 xã vùng đệm huyện Giao Thủy với tổng diện tích tự nhiên 40,18 km². Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng cho việc hoạch định 2.821 ha đất có rừng và đất bãi bồi, hỗ trợ cơ quan quản lý đưa ra giải pháp bảo tồn thích ứng kịp thời.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng lý thuyết diễn thế sinh thái rừng ngập mặn vùng cửa sông ven biển nhiệt đới và lý thuyết mô hình hóa hướng cá thể (Individual-Based Model - IBM). Diễn thế sinh thái tại khu vực được phân định thành hai tuyến nguyên sinh đặc trưng: diễn thế bãi bồi ven biển với vai trò tiên phong của quần xã Mắm biển (Avicennia marina) và diễn thế vùng nước lợ cửa sông với vai trò tiên phong của Bần chua (Sonneratia caseolaris). Các lý thuyết này giải thích quy luật bồi tụ nền đáy, sự tích lũy phù sa của hệ rễ hô hấp và quá trình thay thế loài theo độ cao địa hình cũng như thời gian ngập triều.

Mô hình nghiên cứu kế thừa cấu trúc từ các mô hình sinh thái rừng kinh điển trên thế giới như JABOWA, FORMAN, KiWi và phát triển từ mô hình CGMM (Cần Giờ Mangrove Model). Các khái niệm then chốt trong nghiên cứu bao gồm:

  • Độ ngập triều (Inundation regime): Biến môi trường quyết định phản ứng sinh lý, hô hấp của bộ rễ và khả năng nảy mầm của trụ mầm.
  • Vùng ảnh hưởng cạnh tranh (Field of neighborhood): Không gian dinh dưỡng và ánh sáng mà mỗi cá thể chiếm giữ, quyết định tỷ lệ sống sót và tốc độ tăng trưởng đường kính, chiều cao cây.
  • Biên độ sinh thái độ mặn: Ngưỡng chịu đựng nồng độ muối của thực vật ngập mặn, tối ưu trong khoảng 10‰ đến 25‰.
  • Cấu trúc loài thực thụ và loài gia nhập: Phân loại 17 loài thực vật thuộc 11 họ xuất hiện tại vùng lõi vườn quốc gia.

Phương pháp nghiên cứu

Nghiên cứu sử dụng hệ thống dữ liệu đa nguồn kết hợp giữa khảo sát thực địa, số liệu khí tượng thủy văn trạm Hòn Dáu và kịch bản phát thải trung bình B2 của Bộ Tài nguyên và Môi trường.

Quy trình chọn mẫu và cỡ mẫu điều tra hiện trường được thiết lập chặt chẽ:

  • Cỡ mẫu và bố trí: Thiết lập 3 tuyến băng điều tra đại diện cắt ngang từ bờ ra mép ngoài bãi triều, mỗi băng bố trí 3 ô tiêu chuẩn (OTC), tổng cộng có 9 OTC khảo sát. OTC số 1 của mỗi băng có diện tích 100 m² (10m x 10m) trên các dải rừng trồng giai đoạn 1997 - 1998 với tầng mùn 5 - 15 cm; các OTC số 2 và số 3 có diện tích 25 m² (5m x 5m) được bố trí tại các khu vực rừng trồng 1999 - 2000 (tầng mùn dày 20 - 40 cm) và khu vực tái sinh tự nhiên (tầng mùn 10 - 25 cm).
  • Phương pháp thu thập: Tiến hành đo đếm định kỳ 2 đợt thực địa vào tháng 4 và tháng 10 năm 2012 để đánh giá các chỉ tiêu lâm học: đường kính thân (D1.3), chiều cao vút ngọn (H), đường kính tán lá và mật độ cây tái sinh.
  • Lý do lựa chọn phương pháp phân tích: Phương pháp mô hình hóa hướng cá thể (IBM) lập trình trên môi trường phần mềm MATLAB được lựa chọn thay vì các mô hình biến tổng hợp truyền thống. Lý do là IBM cho phép theo dõi tương tác vi mô phi tuyến tính giữa từng cây riêng lẻ trong không gian 2 chiều, phản ánh chính xác sự cạnh tranh tài nguyên và phản ứng sinh thái của từng loài khi độ ngập thay đổi từng centimét theo chu kỳ nhật triều.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình điều tra thực địa và mô phỏng toán học đã mang lại các phát hiện cốt lõi về hiện trạng và động thái sinh thái rừng ngập mặn Vườn Quốc gia Xuân Thủy:

  • Sự phân hóa cấu trúc nền đất và thảm thực vật: Toàn bộ diện tích đất vùng lõi gồm 1.593 ha đất mặn sú vẹt đước gley sâu (chiếm 56,4% diện tích đất khảo sát), 659 ha đất mặn nhiều gley sâu (chiếm 23,4%) và 569 ha đất cát biển điển hình (chiếm 20,2%). Thảm thực vật ghi nhận 17 loài thuộc 11 họ, trong đó các loài chiếm ưu thế tuyệt đối là Trang (Kandelia candel), Sú (Aegiceras corniculatum) và Bần chua (Sonneratia caseolaris).
  • Tương quan giữa độ dày tầng mùn và tái sinh tự nhiên: Tại các OTC số 1 ở địa hình cao (0,9 - 1,2 m), tầng mùn mỏng (5 - 15 cm), rừng đã khép tán trên 90% nhưng mật độ cây tái sinh rất thấp, dưới 2 cây/m². Ngược lại, ở OTC số 3 với tầng mùn 10 - 25 cm chịu ngập triều định kỳ, số lượng cây con tái sinh tăng gấp 3,5 lần, tạo thành cấu trúc rừng nhiều tầng tán phong phú.
  • Dự báo biến động diện tích dưới tác động của nước biển dâng: Theo kịch bản mực nước biển dâng thêm 30 cm vào năm 2050 và 65 cm vào năm 2100, chế độ ngập triều tại các bãi bồi sẽ tăng từ 2 đến 5 giờ/ngày. Mô phỏng chỉ ra rằng dải rừng tiên phong Mắm biển và Bần chua sẽ buộc phải dịch chuyển lùi sâu vào nội địa từ 150 m đến 300 m. Nếu gặp vật cản đê điều kiên cố ở phía trong, diện tích rừng ngập mặn ngoài đê có nguy cơ suy thoái cục bộ từ 25% đến 35% do hiện tượng ngập úng kéo dài vượt ngưỡng chịu đựng sinh lý.
  • Áp lực từ hoạt động kinh tế xã hội vùng đệm: Khảo sát tại 5 xã vùng đệm cho thấy 69,45% trong tổng số 23.350 lao động hoạt động trong lĩnh vực nông - ngư nghiệp. Gần 500 ha bãi cát pha đã bị người dân chuyển đổi tự phát sang nuôi thả nhuyễn thể (vạng, ngao), gây suy giảm thảm thực vật tự nhiên và cản trở dòng bồi lắng phù sa.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự thay đổi cấu trúc rừng là do chế độ nhật triều thuần nhất với biên độ trung bình 1,7 - 1,9 m (lớn nhất đạt 4,0 m) kết hợp với lưu lượng 114 tỷ m³ nước và 115 triệu tấn phù sa hàng năm từ sông Hồng. Khi mực nước biển dâng cao, thời gian ngập úng kéo dài làm giảm nồng độ oxy trong đất, khiến bộ rễ hô hấp của cây Trang và Vẹt dù bị nghẹt khí, làm suy giảm tốc độ quang hợp ước tính khoảng 30% đến 40%.

So sánh với các nghiên cứu tại rừng ngập mặn Cần Giờ (Nam Bộ), tốc độ tăng trưởng đường kính cây ngập mặn tại Xuân Thủy đạt trung bình 0,4 - 0,6 cm/năm, thấp hơn khoảng 30% so với khu vực phía Nam (0,6 - 0,8 cm/năm). Điều này phản ánh rõ nét tác động bất lợi của mùa đông lạnh miền Bắc (nhiệt độ có thời điểm xuống dưới 10°C) kết hợp với hiện tượng ngọt hóa theo mùa vào thời kỳ lũ sông Hồng.

Dữ liệu mô phỏng từ phần mềm MATLAB được trực quan hóa qua biểu đồ phân bố đường kính thân theo nhóm tuổi (1 đến 15 năm) và các bản đồ không gian thể hiện rõ 3 đai ngập triều: đai dương không ngập (cao trình 1,2 - 1,5 m), đai ngập theo chu kỳ (0,9 - 1,2 m) và đai ngập thường xuyên (0,5 - 0,9 m). Bảng quy hoạch không gian đến các năm 2030, 2050 và 2100 cung cấp bức tranh chi tiết về sự chuyển dịch sinh khối, làm rõ hiện tượng "kẹt sinh thái" khi thảm rừng bị ép chặt giữa mực nước biển dâng và hệ thống đê biển kiên cố.

Đề xuất và khuyến nghị

Dựa trên kết quả mô hình hóa động thái rừng, luận văn đề xuất 4 nhóm giải pháp chiến lược nhằm bảo tồn và phát triển bền vững Vườn Quốc gia Xuân Thủy:

  1. Quy hoạch hành lang thích ứng sinh thái ven biển: Ban Quản lý Vườn Quốc gia Xuân Thủy chủ trì, phối hợp với Sở Nông nghiệp và Phát triển nông thôn tỉnh Nam Định rà soát và mở rộng tối thiểu 200 ha vùng đệm bãi bồi ngoài đê trước năm 2030. Giải pháp này nhằm tạo không gian mở cho rừng ngập mặn tự di cư lấn biển và tiến dần vào đất liền khi mực nước biển dâng thêm 17 - 24 cm.
  2. Trồng phục hồi và làm giàu rừng bằng tập đoàn cây bản địa chịu ngập: Chi cục Kiểm lâm tỉnh Nam Định triển khai kế hoạch trồng bổ sung 150 ha rừng ngập mặn hỗn loài (ưu tiên Bần chua và Mắm biển tại các bãi bồi thấp ven sông Trà, sông Vọp) trong giai đoạn 2026 - 2035, đặt mục tiêu nâng tỷ lệ sống của cây con tái sinh lên trên 85%.
  3. Nạo vét, khơi thông dòng chảy thủy văn liên thông: Công ty Khai thác Công trình Thủy lợi huyện Giao Thủy phối hợp chính quyền địa phương hoàn thành việc nạo vét, thông dòng 12 km luồng lạch trên sông Vọp và sông Trà trước năm 2028. Hoạt động này nhằm điều hòa độ mặn ổn định trong ngưỡng 10‰ đến 22‰, xóa bỏ tình trạng ngọt hóa hoặc mặn hóa cực đoan từng vùng.
  4. Kiểm soát hoạt động nuôi trồng thủy sản và chuyển đổi sinh kế: Ủy ban nhân dân huyện Giao Thủy cùng Ủy ban nhân dân 5 xã vùng đệm thực hiện chuyển đổi 500 ha vùng nuôi ngao, vạng tự phát sang mô hình nuôi sinh thái kết hợp bảo vệ rừng trước năm 2030; đồng thời đào tạo nghề phi nông nghiệp nhằm giảm tỷ lệ lao động khai thác tự nhiên từ 69,45% xuống dưới 45%.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả mô hình hóa của luận văn là tài liệu tham chiếu giá trị cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước và Ban Quản lý các Khu bảo tồn thiên nhiên: Luận văn cung cấp phương pháp luận và dữ liệu dự báo không gian giúp các nhà hoạch định chính sách xây dựng phương án quy hoạch bờ biển, thiết kế đê điều mềm thích ứng biến đổi khí hậu trên quy mô hơn 7.000 ha đất ngập nước ven biển miền Bắc.
  2. Giảng viên, nghiên cứu sinh và học viên chuyên ngành Khoa học Môi trường, Lâm nghiệp: Tài liệu cung cấp quy trình chi tiết về xây dựng mô hình toán học hướng cá thể (IBM) trên nền tảng MATLAB, phục vụ giảng dạy và nghiên cứu mở rộng mô hình sinh thái cho 28 tỉnh thành ven biển Việt Nam.
  3. Các tổ chức phi chính phủ và đơn vị tài trợ quốc tế (IUCN, WWF, Quỹ Môi trường Toàn cầu): Các số liệu định lượng về sinh khối, cấu trúc loài và ngưỡng ngập triều là cơ sở vững chắc để xây dựng các dự án trồng rừng ngập mặn dựa vào cộng đồng (CbCR) với ngân sách từ hàng trăm nghìn đến hàng triệu USD.
  4. Chính quyền địa phương và cộng đồng doanh nghiệp vùng đệm: Phân tích thực trạng kinh tế xã hội giúp chính quyền địa phương xây dựng cơ chế đồng quản lý tài nguyên, định hướng các hộ gia đình khai thác thủy sản bền vững, bảo vệ nguồn thu hàng trăm tỷ đồng mỗi năm từ kinh tế biển.

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

Mô hình hướng cá thể (IBM) mang lại ưu thế gì vượt trội so với các mô hình sinh thái truyền thống?
Mô hình IBM theo dõi trực tiếp từng cá thể cây với các thuộc tính riêng biệt như đường kính thân, chiều cao vút ngọn và bán kính tán thay vì lấy giá trị trung bình toàn diện tích. Điều này giúp phản ánh chính xác sự cạnh tranh không gian giữa các cá thể trong 9 ô tiêu chuẩn khảo sát khi độ ngập thay đổi từng centimét.

Nước biển dâng ảnh hưởng trực tiếp nhất đến rừng ngập mặn Xuân Thủy qua cơ chế nào?
Nước biển dâng làm gia tăng thời gian ngập triều từ 2 đến 5 giờ mỗi ngày, khiến đất bãi bồi rơi vào trạng thái yếm khí kéo dài. Hiện tượng này làm yếu cấu trúc bám của bộ rễ hô hấp ở cây Trang và Sú, làm giảm 30% đến 40% khả năng trao đổi khí và dẫn đến nguy cơ bật gốc khi có sóng lớn.

Tại sao khu vực cửa sông Trà lại hình thành quần xã Bần chua ưu thế thay vì cây Trang hay Vẹt dù?
Cửa sông Trà nhận nguồn nước ngọt và phù sa trực tiếp từ sông Hồng với lưu lượng bình quân 114 tỷ m³/năm, tạo ra môi trường nước lợ có độ mặn thấp (9‰ - 17‰ vào mùa lũ). Bần chua là loài cây tiên phong có biên độ chịu lợ tốt và hệ rễ cọc phát triển mạnh mẽ trên nền bùn lỏng sâu.

Vì sao các loài cây ngập mặn phong phú ở miền Nam như Đước đôi, Dừa nước không phát triển tự nhiên tại Xuân Thủy?
Do đặc thù khí hậu vịnh Bắc Bộ có mùa đông lạnh với nhiệt độ thấp nhất xuống tới 6,8°C và gió mùa Đông Bắc mạnh. Các loài như Đước đôi (Rhizophora apiculata) hay Dừa nước đòi hỏi nhiệt độ trung bình năm trên 24°C không ổn định và không chịu được các đợt rét đậm kéo dài trên 30 ngày.

Hoạt động nuôi vạng và ngao ở vùng đệm tác động tiêu cực thế nào đến hệ sinh thái rừng?
Việc bao bãi nuôi thả tự phát trên diện tích gần 500 ha đã làm thay đổi dòng chảy tự nhiên của lạch triều và cản trở quá trình phát tán trụ mầm của cây ngập mặn. Đồng thời, chất thải hữu cơ và khí sinh học tích tụ ở đáy đầm nuôi làm suy thoái chất lượng đất bãi triều xung quanh.

Kết luận

  • Luận văn đã thiết lập thành công mô hình toán học hướng cá thể (IBM) kết hợp phần mềm MATLAB để mô phỏng tương tác động thái của thảm thực vật rừng ngập mặn Vườn Quốc gia Xuân Thủy dưới tác động của độ ngập nước triều.
  • Xác định rõ mối tương quan giữa 3 đai cao trình địa hình (0,5 m đến 1,5 m), 3 nhóm đất chính và biên độ mặn dao động từ 9‰ đến 32‰ với sự phân bố của 17 loài thực vật ngập mặn.
  • Định lượng hóa nguy cơ suy giảm 25% đến 35% diện tích rừng ngoài đê vào năm 2100 dưới kịch bản nước biển dâng trung bình B2 (dâng 65 cm) nếu không có hành lang di cư sinh thái nội địa.
  • Đóng góp bộ công cụ dự báo định lượng phục vụ công tác quy hoạch phân bố không gian và quản lý bền vững cho 2.821 ha diện tích vùng lõi khu Ramsar Xuân Thủy.
  • Đề xuất lộ trình hành động liên ngành giai đoạn 2026 - 2035 tích hợp giữa bảo tồn đa dạng sinh học, điều tiết thủy văn và chuyển đổi sinh kế bền vững cho 49.224 người dân tại 5 xã vùng đệm.

Nghiên cứu mở ra hướng đi tiên phong trong việc ứng dụng công nghệ mô hình hóa sinh thái vào quản lý tài nguyên thiên nhiên tại Việt Nam. Các cơ quan quản lý và viện nghiên cứu cần tiếp tục ứng dụng, chuẩn hóa dữ liệu viễn thám độ phân giải cao để hoàn thiện các kịch bản thích ứng biến đổi khí hậu trong tương lai.