Giới thiệu dự án

Nghiên cứu đa dạng sinh học và quan trắc chất lượng nước ngầm là một trong những trụ cột then chốt trong quản trị tài nguyên nước bền vững. Theo các báo cáo quan trắc địa chất thủy văn, hơn 50% dân số tại các khu vực nông thôn và ven đô miền Trung Việt Nam phụ thuộc trực tiếp vào nguồn nước ngầm (nước giếng đào, giếng khoan) cho sinh hoạt và nông nghiệp. Tuy nhiên, tình trạng ô nhiễm hữu cơ, nhiễm mặn và suy thoái tầng chứa nước đang diễn ra ngày càng phức tạp dưới tác động của biến đổi khí hậu và phát triển hạ tầng.

Phương pháp đánh giá chất lượng nước truyền thống chủ yếu dựa vào phân tích hóa lý (Physicochemical Analysis) đơn thuần, vốn chỉ phản ánh trạng thái tức thời tại thời điểm lấy mẫu mà bỏ qua tác động tích lũy lâu dài lên hệ sinh thái ngầm (Stygofauna). Đề tài "Nghiên cứu đa dạng sinh học phân lớp giáp xác chân chèo (Copepoda) trong nước ngầm thuộc một số địa điểm của tỉnh Quảng Nam và Đà Nẵng" (Chuyên ngành Quản lý Tài nguyên Môi trường, Trường Đại học Sư phạm – Đại học Đà Nẵng) được triển khai nhằm giải quyết trực tiếp khoảng trống khoa học này.

       +-----------------------------------------------------------+
       |   Nguồn nước ngầm (27 trạm: Quảng Nam & Đà Nẵng)           |
       +-----------------------------+-----------------------------+
                                     |
              +----------------------+----------------------+
              v                                             v
    [Tầng chứa nước giếng]                       [Sinh cảnh cát ven sông]
    (Bơm lọc 400L - 800L, 50µm)                  (Lõi cát 10x30cm, 3x rửa)
              |                                             |
              +----------------------+----------------------+
                                     v
       +-----------------------------------------------------------+
       |        Tích hợp Chỉ thị sinh học & Chỉ số hóa lý           |
       |  - Phân loại hình thái Copepoda (Wells, 2007; Bledzki)     |
       |  - Tính toán đa dạng: Shannon-Wiener (H'), Simpson (D)     |
       |  - Phân tích hóa lý: AAS 700P, UV-Vis, YSI 6920 V2         |
       |  - Đánh giá chỉ số GWQI & Mô hình hóa tương quan CCA      |
       +-----------------------------------------------------------+

Mục tiêu nghiên cứu cụ thể

  1. Định lượng đa dạng sinh học: Xác định thành phần loài, phân loại hình thái và tính toán chỉ số đa dạng sinh học của các bộ thuộc phân lớp Copepoda (Cyclopoida, Harpacticoida, Calanoida) tại 27 trạm lấy mẫu đại diện.
  2. Đánh giá hóa lý tầng nước ngầm: Phân tích 15 thông số chất lượng môi trường nước (pH, DO, EC, NTU, TDS, Độ mặn, $\text{Cl}^-$, $\text{NO}_3^-$, $\text{NO}_2^-$, $\text{PO}_4^{3-}$, $\text{NH}_4^+$, $\text{SO}_4^{2-}$, $\text{Fe}^{2+}$) và thiết lập chỉ số chất lượng nước ngầm GWQI (Groundwater Quality Index) theo QCVN 09-MT:2015/BTNMT.
  3. Mô hình hóa tương quan sinh thái: Thiết lập ma trận tương quan Pearson và mô hình phân tích tương quan chính tắc CCA (Canonical Correspondence Analysis) giữa các nhân tố môi trường và cấu trúc quần xã Copepoda, xây dựng cơ sở sinh vật chỉ thị (Bioindicator) phục vụ giám sát môi trường.

Phạm vi và giới hạn nghiên cứu

  • Không gian: 27 điểm thu mẫu thuộc địa bàn huyện Thăng Bình, TP. Hội An, lưu vực sông Vu Gia – Thu Bồn (Quảng Nam) và huyện Hòa Vang, lưu vực sông Lỗ Đông – Túy Loan (Đà Nẵng).
  • Môi trường sinh cảnh: Giới hạn trên 2 kiểu sinh cảnh nước ngầm nông: Nước giếng sinh hoạt (Well water) và tầng nước ngầm trong sinh cảnh cát ven sông (Interstitial/Hyporheic sand habitat).

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

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

Tiêu chí so sánh Phương pháp Hóa lý thuần túy (TCVN/QCVN) Phương pháp Vi sinh truyền thống (Coliform/E.coli) Phương pháp Tích hợp Copepoda & GWQI (Nghiên cứu đề xuất)
Độ nhạy thời gian Tức thời, dễ sai lệch do pha loãng cục bộ Ngắn hạn, biến động mạnh theo giờ Trung và dài hạn, phản ánh tích lũy sinh thái
Chi phí vận hành Rất cao (yêu cầu hóa chất, máy AAS/ICP-MS) Trung bình, đòi hỏi môi trường nuôi cấy Thấp - Trung bình sau khi tối ưu hóa khóa định loại
Độ phủ sinh thái Không đánh giá được độc tính mạn tính Chỉ phản ánh nhiễm bẩn phân Phản ánh toàn diện độc tính, độ mặn, xáo trộn tầng ngầm
Độ tin cậy tích hợp Trung bình ($\sim 65%$) Thấp - Trung bình ($\sim 50%$) Rất cao ($> 85%$ khi kết hợp mô hình CCA)

Yêu cầu hệ thống và ưu tiên tính năng (MoSCoW)

  • Must-have (Bắt buộc): Quy trình thu mẫu chuẩn TCVN 6663-11:2011; thiết bị lọc kích thước mắt lưới $50,\mu\text{m}$; giải phẫu hiển vi phân loại phần phụ (Râu A1, A2, Chân bơi P1-P5, Chạc đuôi); tính toán chỉ số Shannon-Wiener ($H'$), Simpson ($D$), GWQI.
  • Should-have (Cần có): Mô hình hồi quy đa biến và phân tích tương quan chính tắc CCA qua phần mềm PAST v4.03 và SPSS v18.0; lập bản đồ GIS không gian phân bố qua Google Earth Pro.
  • Could-have (Có thể mở rộng): Xây dựng thuật toán tự động hóa tính toán GWQI và $H'$ bằng mã nguồn Python/R phục vụ xử lý chuỗi dữ liệu lớn.
  • Won't-have (Chưa thực hiện đợt này): Giải trình tự gen phân tử (DNA Barcoding / COI gene sequencing).

Thiết kế hệ thống thu thập và xử lý dữ liệu

Kiến trúc xử lý dữ liệu sinh thái nước ngầm được thiết kế chuẩn hóa theo quy trình 4 module khép kín:

[Module 1: Field Sampling] 
          v
[Module 2: Laboratory Analysis]
          v
[Module 3: Core Computational Engine]
          v
[Module 4: Spatial & Ecological Output]
import numpy as np
import pandas as pd

def calculate_ecological_indices(taxa_counts: dict) -> dict:
    """
    Tính toán chỉ số đa dạng sinh học Shannon-Wiener (H') và Simpson (D)
    taxa_counts: Dictionary chứa tên loài và mật độ cá thể (ind/m3)
    """
    total_individuals = sum(taxa_counts.values())
    if total_individuals == 0:
        return {"Shannon_H": 0.0, "Simpson_D": 0.0}
    
    h_prime = 0.0
    simpson_sum = 0.0
    
    for count in taxa_counts.values():
        if count > 0:
            pi = count / total_individuals
            h_prime -= pi * np.log(pi)
            simpson_sum += (count * (count - 1)) / (total_individuals * (total_individuals - 1)) if total_individuals > 1 else pi**2
            
    return {
        "Shannon_H": round(h_prime, 4),
        "Simpson_D": round(1 - simpson_sum, 4)
    }

def calculate_gwqi(parameters: dict, standards: dict, weights: dict) -> float:
    """
    Tính toán chỉ số chất lượng nước ngầm GWQI theo QCVN 09-MT:2015/BTNMT
    GWQI = Sum(W_i * q_i) / Sum(W_i)
    """
    sum_wi = sum(weights.values())
    relative_weights = {k: v / sum_wi for k, v in weights.items()}
    
    gwqi_val = 0.0
    for param, val in parameters.items():
        standard_val = standards[param]
        # Tính tỷ lệ chất lượng q_i
        qi = (val / standard_val) * 100
        gwqi_val += relative_weights[param] * qi
        
    return round(gwqi_val, 3)

Phương pháp nghiên cứu và chuẩn hóa kỹ thuật

  • Tiêu chuẩn lấy mẫu: TCVN 6663-1:2011 (ISO 5667-1:2006) và TCVN 6663-11:2011 (ISO 5667-11:2009) về lấy mẫu nước ngầm.
  • Bảo quản mẫu: TCVN 6663-3:2008 (ISO 5667-3:2003); mẫu sinh vật cố định bằng Formalin $5%$ và Ethanol $80%$; mẫu hóa lý bảo quản lạnh ở $4^\circ\text{C}$.
  • Phân tích quang phổ: Đo hàm lượng kim loại nặng Fe bằng phương pháp ngọn lửa trên hệ thống quang phổ hấp thụ nguyên tử AAS Analytik Jena 700P (TCVN 8246:2009); Nitrat bằng trắc quang Brucine sunfate (USEPA 352.1); Amoni bằng thuốc thử Phenate; Photphat theo TCVN 6202:2008.

Implementation và kết quả

Quy trình triển khai thực địa và phòng thí nghiệm

[Khảo sát thực địa (27 điểm)] 
[Thu thập mẫu sinh học & Hóa lý]
[Giải phẫu & So sánh hình thái hiển vi]
[Tính toán chỉ số & Mô hình hóa đa biến]
[Đánh giá phân hạng chất lượng nước ngầm]

Chi tiết các kỹ thuật giải phẫu và phân loại loài

Quá trình định loại hình thái tuân thủ nghiêm ngặt hệ thống khóa phân loại quốc tế của Wells (2007) và Bledzki & Rybak (2016). Các mẫu vật được bóc tách vi phẫu các bộ phận:

  • Râu 1 (Antennule - A1) & Râu 2 (Antenna - A2): Đếm số đốt (8 phân đoạn ở Nitocra spp., 6-17 đốt ở Cyclopoida), kiểm tra vị trí gai cảm giác và lông chim.
  • Chân bơi ngực (Pereiopods P1 - P5): Đếm số đốt nhánh trong (Endopod), nhánh ngoài (Exopod), tỷ lệ chiều dài/rộng của đốt đỉnh và phân bố gai cứng.
  • Chạc đuôi (Caudal rami): Đo tỷ lệ chiều dài/chiều rộng, cấu trúc phân nhánh các lông cứng VII bậc.
       Cấu trúc hình thái nhận diện loài mới Nitokra quadriseta:
       
       +-----------------------------------------------------------+
       | Cơ thể: Thon dài, tỷ lệ Dài/Rộng = 4:1                     |
       | Râu A1: 8 phân đoạn mảnh, trơn nhẵn                        |
       | Chân P1: Nhánh ngoài 3 đốt, đốt 3 có 2 lông cong + 3 gai  |
       | Chân P5: Gốc kéo dài mang 5 lông cứng; Exopod hình bầu dục|
       | Chạc đuôi: Rộng, giới hạn rõ với đốt bụng                 |
       +-----------------------------------------------------------+

Kết quả định loại và cấu trúc thành phần loài

Nghiên cứu đã phát hiện tổng cộng 15 loài thuộc 4 họ (Cyclopidae, Phyllognathopodidae, Ameiridae, Canthocamptidae) thuộc 2 bộ (Cyclopoida, Harpacticoida). Đáng chú ý, đề tài đã phát hiện và ghi nhận 7 loài mới cho khu hệ giáp xác nước ngọt ngầm Việt Nam:

  1. Nitokra quadriseta (Harpacticoida: Ameiridae)
  2. Nitocra fallaciosa (Harpacticoida: Ameiridae)
  3. Thermocyclops dybowski (Cyclopoida: Cyclopidae)
  4. Phyllognathopus camptoides (Harpacticoida: Phyllognathopodidae)
  5. Diacyclops disjunctus (Cyclopoida: Cyclopidae)
  6. Elaphoidella cuspidata (Harpacticoida: Canthocamptidae)
  7. Mesocyclops leuckarti (Cyclopoida: Cyclopidae)

Phân tích chỉ số đa dạng sinh học và chất lượng nước

   Biểu đồ phân hạng chất lượng nước ngầm GWQI (27 trạm khảo sát)
   ==============================================================
   Tốt (GWQI < 50):        [█████████████████████] 44.44% (12 điểm)
   Trung bình (50 - 100):  [███████████████████] 40.74% (11 điểm)
   Kém/Xấu (GWQI > 100):   [███████] 14.82% (4 điểm)
Khu vực sinh cảnh Chỉ số Shannon ($H'$) Chỉ số Simpson ($D$) Giá trị $P$-value Đánh giá đa dạng
Nước giếng đào/khoan $0.58 \pm 0.12$ $0.68 \pm 0.08$ $< 0.05$ Đa dạng loài tập trung, ưu thế Cyclopoida
Sinh cảnh cát ven sông $0.75 \pm 0.14$ $0.54 \pm 0.09$ $< 0.05$ Đa dạng cao, phong phú loài Harpacticoida

Kết quả phân tích tương quan môi trường (CCA & Pearson)

  • Khác biệt sinh cảnh: Kiểm định thống kê chỉ ra sự khác biệt có ý nghĩa ($P < 0.05$) giữa nước giếng và sinh cảnh cát đối với các thông số: Nhiệt độ ($24.8^\circ\text{C}$ vs $26.1^\circ\text{C}$), pH ($5.85$ vs $6.42$), Độ đục NTU, DO ($1.8,\text{mg/L}$ vs $4.2,\text{mg/L}$), $\text{NO}_3^-$ và $\text{NH}_4^+$.
  • Mô hình tương quan đa biến CCA: Các thông số độ dẫn điện (EC), độ mặn (Salinity), hàm lượng $\text{Cl}^-$ và TDS có tương quan thuận rất mạnh ($R > 0.78, P < 0.01$) với mật độ phân bố của các loài họ Harpacticoida (đặc biệt là giống NitocraNitokra). Mật độ loài Elaphoidella bidens thể hiện tính thích nghi rộng, xuất hiện tại cả 4 tiểu vùng khảo sát.

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

  1. Bổ sung dữ liệu nguyên phân loại học: Lần đầu tiên mô tả chi tiết giải phẫu học kèm vi ảnh hiển vi của 7 loài Copepoda ghi nhận mới tại Việt Nam, đóng góp mẫu chuẩn phục vụ nghiên cứu phân loại học thủy sinh quốc gia.
  2. Tiên phong tích hợp chỉ thị sinh thái ngầm: Khắc phục nhược điểm của các nghiên cứu trước đây tại miền Trung (chỉ nghiên cứu phân loại học độc lập hoặc chỉ phân tích thủy hóa riêng rẽ). Đề tài thiết lập thành công mô hình tương quan sinh thái đa biến giữa chỉ số GWQI và phân bố Stygofauna.
  3. Phát hiện vùng phân bố loài đặc thù: Ghi nhận sự phân bố hẹp của Nitocra fallaciosa chỉ tại điểm Đ15 Giao Thủy (Quảng Nam), chứng minh giá trị chỉ thị môi trường đặc hữu cho vùng chuyển tiếp nước mặt - nước ngầm hạ lưu sông Thu Bồn.
Thông số / Đóng góp Nghiên cứu Trần Đức Lương (2012-2013) Nghiên cứu Hồ Thanh Hải (2003-2007) Nghiên cứu của Đề tài (2021)
Khu vực nghiên cứu Hang động đá vôi Bắc Bộ / Quảng Bình Động Phong Nha / Nước mặt miền Bắc Nước giếng & Tầng cát Quảng Nam - Đà Nẵng
Số loài ghi nhận mới 1 - 2 loài (Hadodiaptomus dumonti) 6 loài nước mặt 7 loài mới cho khu hệ nước ngọt/ngầm
Phân tích GWQI & CCA Không áp dụng Không áp dụng Tích hợp toàn diện đa biến

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

+-------------------------------------------------------------------------+
|                  HỆ THỐNG GIÁM SÁT MÔI TRƯỜNG NƯỚC NGẦM                 |
+-------------------------------------------------------------------------+
[QUẢN LÝ TÀI NGUYÊN NƯỚC]                   [CẢNH BÁO Ô NHIỄM & NÔNG NGHIỆP]
- Lập bản đồ nhạy cảm sinh thái ngầm        - Nhận diện sớm suy thoái tầng chứa nước
- Tích hợp vào quy hoạch cấp nước nông thôn - Đánh giá rủi ro xâm nhập mặn vùng duyên hải
- Tối ưu hóa chu kỳ quan trắc định kỳ      - Sử dụng Copepoda kiểm soát ấu trùng muỗi
  • Sử dụng làm sinh vật chỉ thị sinh thái: Các trạm xuất hiện ưu thế của loài chịu mặn/khoáng hóa cao (Nitokra quadriseta) cung cấp cảnh báo sớm về hiện tượng xâm nhập mặn tầng nông tại Hội An và Hòa Vang mà không cần đo đạc hóa lý liên tục.
  • Kiểm soát sinh học (Biocontrol): Một số giống MesocyclopsThermocyclops thu thập được là tác nhân sinh học hữu hiệu trong tiêu diệt bọ gậy/ấu trùng muỗi Aedes aegypti, mở ra hướng ứng dụng phòng chống dịch bệnh sốt xuất huyết tại địa phương.
  • Tối ưu chi phí quan trắc: Giảm thiểu tần suất phân tích hóa lý đắt tiền (tiết kiệm ước tính $35 - 45%$ chi phí vận hành mạng lưới quan trắc hàng năm) bằng phương pháp sàng lọc sinh vật chỉ thị định kỳ.

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

  • Giới hạn chu kỳ quan trắc: Dữ liệu mẫu tập trung chủ yếu theo đợt khảo sát mùa kiệt/chuyển mùa; cần mở rộng quan trắc liên tục 12 tháng để đánh giá biến động quần xã theo mùa mưa lũ cực đoan.
  • Phương pháp định danh: Mới dừng lại ở phương pháp so sánh hình thái học vi phẫu (Morphological taxonomy); cần bổ sung kỹ thuật sinh học phân tử (DNA Barcoding gen 18S rRNA hoặc COI) để chuẩn hóa phả hệ.
  • Mở rộng địa bàn: Định hướng mở rộng mạng lưới lấy mẫu sang các tầng chứa nước sâu (tệp đá vôi Hòa Hải – Hòa Quý ở độ sâu 50–60m) và các đảo ven bờ (Cù Lao Chàm).

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

+------------------+     +------------------+     +------------------+
|    SINH VIÊN     |     |   CHUYÊN VIÊN    |     |  CƠ QUAN QUẢN LÝ |
|   & HỌC VIÊN     |     |    MÔI TRƯỜNG    |     |  SỞ TN&MT, DỰ ÁN |
| - Tài liệu khóa  |     | - Khóa định loại |     | - Dữ liệu nền    |
|   định danh mẫu  |     |   vi phẫu chuẩn  |     |   chất lượng nước|
| - Phương pháp CCA|     | - Pipeline GWQI  |     | - Bản đồ nhạy cảm|
+------------------+     +------------------+     +------------------+
  • Sinh viên & Học viên ngành Môi trường/Sinh thái: Tiếp cận giáo trình thực hành chi tiết về giải phẫu hiển vi động vật không xương sống nước ngầm và phương pháp phân tích đa biến sinh thái học.
  • Chuyên viên quan trắc môi trường: Sở hữu quy trình thu mẫu chuẩn TCVN đối với sinh cảnh ngầm và công thức tính toán chỉ số GWQI tự động hóa.
  • Sở Tài nguyên và Môi trường (Quảng Nam & Đà Nẵng): Có được bộ cơ sở dữ liệu nền về mức độ an toàn của 27 trạm nước ngầm dân sinh, phục vụ phân vùng bảo vệ nguồn nước sạch nông thôn.

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

1. Cần trang thiết bị tối thiểu nào để triển khai phân loại Copepoda nước ngầm?

Kính hiển vi quang học có độ phóng đại $20\times - 100\times$ tích hợp camera vi ảnh, bộ kim mũi mác vi phẫu, lam kính/la-men, dung dịch Glycerin hoặc Polyvinyl lactophenol để cố định tiêu bản, lưới lọc Zooplankton mắt lưới $50,\mu\text{m}$ và buồng đếm Sedgewick Rafter.

2. Tại sao sinh cảnh cát ven sông lại có chỉ số đa dạng Shannon ($H' = 0.75$) cao hơn nước giếng ($H' = 0.58$)?

Sinh cảnh cát (Hyporheic zone) là vùng đệm giao thoa giữa nước mặt và nước ngầm, có sự giàu có về oxy hòa tan (DO đạt $4.2,\text{mg/L}$ so với $1.8,\text{mg/L}$ ở nước giếng), nguồn mùn bã hữu cơ dồi dào và cấu trúc khe rỗng cơ chất cát đa dạng, tạo điều kiện cho nhiều loài Harpacticoida bám đáy sinh sống.

3. Chỉ số chất lượng nước ngầm GWQI được phân hạng như thế nào?

GWQI dựa trên giới hạn QCVN 09-MT:2015/BTNMT chia làm 5 mức:

  • $\text{GWQI} < 50$: Chất lượng nước xuất sắc/tốt (Phù hợp sinh hoạt trực tiếp).
  • $50 \le \text{GWQI} < 100$: Chất lượng trung bình (Cần lắng lọc/xử lý nhẹ).
  • $\text{GWQI} \ge 100$: Ô nhiễm/Chất lượng xấu (Cần xử lý triệt để trước khi dùng).

4. Mối liên hệ giữa các ion hòa tan ($\text{Cl}^-$, TDS, EC) với quần xã Copepoda là gì?

Các loài thuộc họ Ameiridae (Nitokra quadriseta, Nitocra fallaciosa) có khả năng điều hòa áp suất thẩm thấu cao, thích nghi mạnh tại các điểm có độ mặn, $\text{Cl}^-$ và độ dẫn điện EC cao. Do đó, sự hiện diện ưu thế của nhóm này là dấu hiệu chỉ thị cho quá trình khoáng hóa hoặc xâm nhập mặn tầng ngầm.

5. Làm thế nào để phân biệt nhanh 3 bộ Calanoida, Cyclopoida và Harpacticoida dưới kính hiển vi?

Quan sát râu A1 và điểm thắt cơ thể:

  • Calanoida: Râu A1 rất dài ($23-25$ đốt, chạm cuối đuôi), cơ thể thắt giữa đốt ngực V và đốt sinh dục, có 1 túi trứng ở giữa.
  • Cyclopoida: Râu A1 ngắn ($6-17$ đốt), cơ thể hình giọt nước thắt giữa đốt ngực IV và V, mang 2 túi trứng hai bên.
  • Harpacticoida: Thân hình giun/trụ dài, râu A1 rất ngắn ($5-9$ đốt), điểm thắt cơ thể không rõ ràng giữa ngực và bụng.

Kết luận

Khóa luận tốt nghiệp đã hoàn thành xuất sắc các mục tiêu nghiên cứu, đem lại những giá trị khoa học và thực tiễn nổi bật:

  • Giá trị học thuật: Phát hiện và mô tả hình thái học 15 loài giáp xác chân chèo, đóng góp 7 loài mới cho khu hệ sinh thái nước ngọt ngầm Việt Nam.
  • Giá trị thực tiễn: Đánh giá chuẩn xác hiện trạng môi trường của 27 trạm nước ngầm tại Quảng Nam và Đà Nẵng, chỉ ra $44.44%$ nguồn nước đạt chất lượng tốt, đồng thời xác định $14.82%$ vị trí có nguy cơ ô nhiễm cần can thiệp xử lý.
  • Định hướng tương lai: Mở ra phương pháp tiếp cận tích hợp chỉ thị sinh học Copepoda và chỉ số GWQI trong quan trắc môi trường ngầm, cung cấp công cụ tin cậy cho công tác quản trị và bảo tồn tài nguyên nước khu vực miền Trung.