Tổng quan nghiên cứu

Hệ sinh thái rừng ngập mặn Cần Giờ có tổng diện tích tự nhiên đạt 704,22 km2, chiếm gần một phần ba diện tích Thành phố Hồ Chí Minh, với tỷ lệ che phủ của rừng và đất rừng đạt trên 54%. Được UNESCO công nhận là Khu Dự trữ Sinh quyển Thế giới vào năm 2000, khu vực này giữ vai trò là lá phổi xanh và quả thận sinh thái lọc sạch chất thải từ lưu vực sông Đồng Nai và sông Sài Gòn. Tuy nhiên, áp lực nhân sinh từ các vùng nuôi trồng thủy sản tập trung và hoạt động xả thải công nghiệp dọc lưu vực sông Thị Vải đang đe dọa trực tiếp đến chất lượng môi trường nền đáy. Trong bối cảnh đó, nhóm động vật đáy không xương sống cỡ trung bình, đặc biệt là tuyến trùng sống tự do (Nematoda) chiếm 90% đến 95% tổng số lượng cá thể sinh vật đáy, trở thành công cụ quan trắc sinh học lý tưởng.

Mục tiêu cốt lõi của nghiên cứu là xác định đặc điểm cấu trúc thành phần loài, phân bố mật độ và các chỉ số đa dạng sinh học của quần xã tuyến trùng sống tự do, đồng thời áp dụng phương pháp đường cong ABC và chỉ số môi trường W để đánh giá hiện trạng xáo trộn môi trường trầm tích. Phạm vi nghiên cứu được thực hiện tại ba phân khu đặc trưng gồm vùng lõi bảo tồn, vùng nuôi trồng thủy sản và vùng công nghiệp ven sông Thị Vải qua hai đợt khảo sát thực địa đại diện cho mùa khô (tháng 4 năm 2016) và mùa mưa (tháng 9 năm 2016). Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp bộ cơ sở dữ liệu định lượng gồm 80 loài tuyến trùng và các ngưỡng biến động mật độ từ 166 đến 1395 cá thể trên 10 cm2, đóng góp công cụ định lượng hỗ trợ các cơ quan quản lý thiết lập chiến lược bảo tồn và giám sát sức khỏe sinh thái bền vững.

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 vận dụng hệ thống phân loại tuyến trùng hiện đại của De Ley (2004) kết hợp cấu trúc phân loại dưới họ của Lorenzen (1994) và các khóa định loại hình thái học của Warwick (1998) nhằm định danh chính xác các đơn vị phân loại đến mức giống và loài. Bên cạnh đó, mô hình đánh giá xáo trộn sinh thái thông qua đường cong tích lũy số lượng - sinh khối ABC (Abundance Biomass Comparison) và chỉ số thống kê W của Warwick (1986) cùng Clarke & Warwick (2001) được sử dụng để phát hiện trạng thái ô nhiễm môi trường nền đáy. Khi môi trường ổn định, đường sinh khối nằm phía trên đường số lượng cá thể (chỉ số W dương từ 0 đến +1); ngược lại, khi chịu tác động ô nhiễm nặng, đường số lượng sẽ vượt lên trên đường sinh khối (chỉ số W âm từ -1 đến 0).

Nền tảng lý thuyết còn tích hợp các khái niệm then chốt gồm: Tuyến trùng sống tự do (Free-living marine nematodes), Sinh vật chỉ thị môi trường (Bioindicators), Động vật đáy cỡ trung bình (Meiobenthos có kích thước rây lọc từ 38 đến 1000 µm), Độ phong phú loài Margalef, Độ đồng đều Pielou và Đa dạng sinh học Shannon-Wiener.

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

Nguồn dữ liệu sơ cấp được thu thập từ 18 mẫu trầm tích tầng mặt tại 3 khu vực nghiên cứu đại diện cho các mức độ tác động nhân sinh khác nhau: Vùng lõi tự nhiên (tọa độ 10°30'17''B - 106°55'44''Đ), Vùng nuôi trồng thủy sản (tọa độ 10°29'39''B - 106°51'41''Đ) và Vùng công nghiệp ven sông Thị Vải (tọa độ 10°31'53''B - 106°59'36''Đ). Cỡ mẫu tại mỗi khu vực gồm 3 lần lặp lại độc lập cho mỗi đợt khảo sát vào tháng 4 và tháng 9 năm 2016. Phương pháp chọn mẫu tầng ngẫu nhiên được áp dụng khi triều rút bằng ống trụ trong suốt có đường kính 3,5 cm, dài 40 cm cắm sâu 10 cm vào lớp bùn đáy.

Lý do lựa chọn phương pháp tách chiết tuyển nổi bằng dung dịch keo silica Ludox TM50 (tỷ trọng 1,18 g/l) kết hợp cố định Formalin-Axit Axetic (FAA) và làm trong theo De Grisse cải tiến là giúp thu hồi nguyên vẹn các cá thể giun tròn siêu nhỏ mà không làm biến dạng cấu trúc vi phẫu. Dữ liệu sau đó được phân tích thống kê đa biến trên phần mềm chuyên dụng PRIMER phiên bản 6, bao gồm phân tích thành phần chính (PCA) cho các biến thủy lý hóa, phân tích độ tương đồng ANOSIM với ma trận khoảng cách Bray-Curtis, kỹ thuật sắp xếp đa chiều không gian 2D-MDS và phân tích cụm thứ bậc Cluster.

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

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

Nghiên cứu đã ghi nhận tổng cộng 80 loài và dạng loài tuyến trùng sống tự do thuộc 20 họ và 6 bộ tại rừng ngập mặn Cần Giờ. Trong đó, bộ Monhysterida chiếm ưu thế tuyệt đối về thành phần loài với 25 loài (chiếm 31,25%), tiếp theo là bộ Araeolaimida với 15 loài (18,75%), bộ Enoplida với 12 loài (15,00%), bộ Desmoscolecida với 10 loài (12,50%), bộ Plectida với 9 loài (11,25%) và bộ Chromadorida với 9 loài (11,25%).

Mật độ cá thể trung bình ghi nhận sự biến động rõ rệt giữa hai mùa và giữa các phân khu khảo sát. Tại vùng nuôi trồng thủy sản, mật độ tuyến trùng vào mùa mưa tăng vọt lên mức 1395 ± 34 cá thể trên 10 cm2, cao gấp gần 4 lần (tăng trưởng 296,3%) so với mật độ mùa khô là 352 ± 68 cá thể trên 10 cm2. Tại vùng lõi bảo tồn, mật độ tăng từ 172 ± 83 cá thể trên 10 cm2 trong mùa khô lên 398 ± 20 cá thể trên 10 cm2 trong mùa mưa. Trái lại, tại khu vực chịu ảnh hưởng công nghiệp, mật độ giảm nhẹ từ 200 ± 59 cá thể trên 10 cm2 vào mùa khô xuống 166 ± 24 cá thể trên 10 cm2 vào mùa mưa.

Cấu trúc ưu thế của quần xã có sự phân hóa theo mùa: Vào mùa khô, 7 họ tuyến trùng chiếm tới 92% tổng số cá thể, dẫn đầu là Desmodoridae (23%) và Chromadoridae (18%), với loài ưu thế nổi bật là Metachromadora minor (15%). Sang mùa mưa, cấu trúc dịch chuyển sang họ Comesomatidae (24%) và Linhomoeidae (23%), trong đó các loài thuộc giống Paracomesoma chiếm ưu thế áp đảo. Kiểm định ANOSIM khẳng định sự khác biệt về cấu trúc quần xã giữa hai mùa có ý nghĩa thống kê rất cao với hệ số R toàn cầu bằng 1 (mức ý nghĩa p = 0,1%), trong khi sự sai khác giữa ba vùng thu mẫu là không đáng kể (hệ số R đều nhỏ hơn 0,1; p > 15%).

Thảo luận kết quả

Sự biến động mạnh mẽ của mật độ và cấu trúc loài bắt nguồn từ sự thay đổi của các yếu tố thủy lý hóa và cơ giới trầm tích. Vào mùa mưa, lượng nước ngọt khổng lồ từ thượng lưu đổ về làm độ mặn suy giảm từ ngưỡng 23,5‰ - 25,0‰ ở mùa khô xuống còn 15,0‰ - 21,0‰, đồng thời độ pH tăng từ khoảng 6,94 - 6,98 lên mức kiềm nhẹ 7,11 - 8,02. Cùng với đó, tỷ lệ hạt bùn mịn kích thước dưới 63 µm luôn chiếm trên 40% ở tất cả các điểm khảo sát và đạt cực đại 72,2% tại vùng lõi vào mùa mưa. Lượng mùn bã hữu cơ dồi dào từ thức ăn dư thừa và chất thải nuôi tôm tại vùng nuôi trồng đã kích thích sự bùng nổ mật độ của các nhóm tuyến trùng ăn vi khuẩn và mùn bã.

Khi biểu diễn trên đồ thị đa chiều 2D-MDS và biểu đồ phân cụm Cluster, dữ liệu sinh học phân tách thành hai nhóm mùa hoàn toàn riêng biệt tại ngưỡng tương đồng 40%, thể hiện mức độ dị biệt lên tới 60%. Mức độ co cụm điểm dày đặc trong mùa khô phản ánh tính ổn định môi trường cao hơn sự biến động phân tán của mùa mưa. Kết quả này tương đồng với các nghiên cứu sinh thái tại rừng ngập mặn Vịnh Gazi ở Kenya ghi nhận mật độ từ 223 đến 1241 cá thể trên 10 cm2 và vịnh Nha Trang. Đồ thị đường cong k-dominance và đường cong ABC cho thấy sự hiện diện ổn định của các giống chỉ thị bùn mịn chịu mặn như Terschellingia, Paracomesoma và Ptycholaimellus, phản ánh chất lượng nền đáy tại vùng lõi vẫn duy trì trạng thái cân bằng tự nhiên trong khi vùng nuôi trồng và công nghiệp bộc lộ các dấu hiệu xáo trộn sinh thái cục bộ.

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

  1. Thiết lập mạng lưới quan trắc sinh học định kỳ: Sở Tài nguyên và Môi trường Thành phố Hồ Chí Minh phối hợp cùng các viện nghiên cứu chuyên ngành triển khai chương trình quan trắc sinh học nền đáy định kỳ 2 lần mỗi năm (vào tháng 4 và tháng 9), sử dụng bộ chỉ số tuyến trùng (H', Margalef d, đường cong ABC) tại 100% các điểm xả thải công nghiệp và vùng nuôi tôm trọng điểm ven sông Thị Vải, hoàn thành trong giai đoạn 2026-2027.
  2. Kiểm soát và nâng cấp hệ thống xử lý nước thải thủy sản: Ủy ban nhân dân huyện Cần Giờ cùng Phòng Kinh tế địa phương chỉ đạo bắt buộc 100% cơ sở nuôi tôm thương phẩm xây dựng hệ thống ao lắng sinh học đạt chuẩn, giảm tối thiểu 35% hàm lượng chất hữu cơ xả thải trực tiếp ra hệ thống kênh rạch tự nhiên trước quý IV năm 2027.
  3. Tăng cường giám sát liên tỉnh xả thải công nghiệp lưu vực sông: Lực lượng Cảnh sát Môi trường và Thanh tra Môi trường Thành phố Hồ Chí Minh phối hợp tỉnh Đồng Nai lắp đặt hệ thống quan trắc tự động nước thải tại các nhà máy ven sông Thị Vải, đặt mục tiêu duy trì chỉ số sinh thái W của quần xã động vật đáy luôn đạt giá trị dương trên 0,15 từ năm 2026 đến 2030.
  4. Xây dựng ngân hàng dữ liệu số về sinh vật chỉ thị môi trường: Viện Sinh thái và Tài nguyên sinh vật chủ trì đề án số hóa cơ sở dữ liệu mã vạch DNA và hình thái học cho hơn 80 loài tuyến trùng tại Cần Giờ, tích hợp phần mềm cảnh báo sớm suy thoái sinh thái nền đáy rừng ngập mặn trước năm 2028.

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

  1. Các nhà quản lý tài nguyên và cơ quan bảo tồn thiên nhiên: Sử dụng quy trình ứng dụng chỉ số môi trường W và đường cong ABC từ luận văn để xây dựng khung đánh giá sức khỏe sinh thái cho Khu Dự trữ Sinh quyển Cần Giờ và hệ thống rừng phòng hộ ven biển.
  2. Giảng viên, nghiên cứu sinh và sinh viên chuyên ngành Sinh học, Sinh thái học, Thủy sinh vật học: Khai thác chi tiết phương pháp tách mẫu bằng Ludox TM50, kỹ thuật lên tiêu bản Seinhorst và tài liệu mô tả phân loại của 80 loài giun tròn biển phục vụ công tác giảng dạy và nghiên cứu mở rộng.
  3. Các đơn vị tư vấn lập Báo cáo Đánh giá tác động môi trường (ĐTM): Kế thừa bộ số liệu nền về mật độ (166 đến 1395 cá thể trên 10 cm2) và chỉ số đa dạng loài để định lượng tác động sinh thái nền đáy cho các dự án nạo vét luồng hàng hải, phát triển cảng biển và khu công nghiệp ven duyên hải.
  4. Doanh nghiệp và hợp tác xã nuôi trồng thủy sản ven biển: Ứng dụng các chỉ báo về sự tích tụ mùn bã hữu cơ đáy ao để tối ưu hóa quy trình xiphong, kiểm soát bùn đáy và duy trì cân bằng hệ vi sinh vật có lợi cho tôm cá.

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

  1. Tuyến trùng sống tự do đóng vai trò gì trong việc chỉ thị chất lượng môi trường nền đáy? Tuyến trùng biển có chu kỳ sống ngắn, sinh khối lớn chiếm hơn 90% động vật đáy và phản ứng rất nhạy với các biến đổi lý hóa. Tỷ lệ giữa các giống nhạy cảm và giống chống chịu giúp phát hiện nhanh tình trạng ô nhiễm hữu cơ và kim loại nặng mà các chỉ tiêu hóa lý đơn thuần dễ bỏ sót.

  2. Nguyên nhân nào khiến mật độ tuyến trùng tại vùng nuôi trồng thủy sản tăng vọt vào mùa mưa? Vào mùa mưa, lưu lượng nước ngọt hòa loãng độ mặn kết hợp với nguồn thức ăn dư thừa dồi dào từ bùn thải ao tôm tạo môi trường dinh dưỡng lý tưởng. Điều này thúc đẩy các loài ăn mùn bã và vi khuẩn thuộc họ Comesomatidae bùng nổ số lượng cá thể lên tới 1395 cá thể trên 10 cm2.

  3. Đâu là sự khác biệt cấu trúc cốt lõi của quần xã tuyến trùng giữa hai mùa nghiên cứu? Sự khác biệt mang ý nghĩa thống kê tuyệt đối với hệ số R toàn cầu đạt 1 theo phân tích ANOSIM. Mùa khô ưu thế bởi các họ Desmodoridae (23%) và Chromadoridae (18%), trong khi mùa mưa dịch chuyển hoàn toàn sang ưu thế của họ Comesomatidae (24%) và Linhomoeidae (23%).

  4. Phương pháp tách lọc tuyến trùng bằng dung dịch Ludox TM50 có ưu điểm gì vượt trội? Dung dịch keo Ludox TM50 điều chỉnh tỷ trọng chuẩn 1,18 g/l cho phép tách nổi chọn lọc toàn bộ tuyến trùng ra khỏi trầm tích cát bùn thô mà không làm hư hại cấu trúc giải phẫu của giun tròn, đảm bảo độ chính xác tuyệt đối khi định danh hình thái dưới kính hiển vi quang học.

  5. Đường cong ABC và chỉ số W phản ánh trạng thái xáo trộn sinh thái như thế nào? Đường cong ABC so sánh đường tích lũy sinh khối và mật độ cá thể. Khi môi trường không bị xáo trộn, đường sinh khối nằm trên đường mật độ (chỉ số W dương). Khi ô nhiễm gia tăng, các loài nhỏ có khả năng chịu đựng chiếm ưu thế, đẩy đường mật độ vượt lên trên và chỉ số W chuyển sang âm.

Kết luận

  • Luận văn đã xác lập bộ cơ sở dữ liệu phân loại học toàn diện với 80 loài tuyến trùng biển thuộc 20 họ và 6 bộ tại ba phân khu đặc trưng của rừng ngập mặn Cần Giờ.
  • Chứng minh tính phân hóa mùa sâu sắc trong cấu trúc quần xã thông qua kiểm định thống kê ANOSIM với giá trị R toàn cầu đạt 1 (p = 0,1%).
  • Ghi nhận sự biến động mật độ cá thể mạnh mẽ tại vùng nuôi trồng thủy sản, tăng gần gấp 4 lần từ 352 cá thể lên 1395 cá thể trên 10 cm2 khi chuyển mùa.
  • Khẳng định giá trị ứng dụng thực tiễn của chỉ số W và đường cong ABC trong việc đánh giá nhanh và chính xác mức độ xáo trộn sinh thái nền đáy.
  • Hoàn thiện lộ trình mở rộng quan trắc sinh học định kỳ và số hóa dữ liệu nhận dạng tuyến trùng chỉ thị phục vụ công tác quản lý tài nguyên bền vững đến năm 2030.