STT Thông số Đơn vị Phương pháp xác định Số hiệu tiêu chuẩn, phương pháp
1 Tỉ trọng g/cm3
Sấy ở 105oC cho đến khi mẫu đất
khơng cịn thay đổi trọng lượng, cân khối lượng khô
2 Ẩm độ %
Sấy ở 105oC cho đến khi mẫu đất
khơng cịn thay đổi trọng lượng, cân khối lượng khô
3 pH
Trích bằng nước cất theo tỷ lệ 1:2.5 (đất: nước). sử dụng máy đo pH WalkLAB TI9000
TCVN 5979:2007 (ISO 10390:2005); TCVN
4402:1987; TCVN 4401:1987
4 CHC % Phương pháp Walkley Black TCVN 6642:2000;
TCVN 6644:2000
5 TN mg/L Phương pháp Kjeldahl TCVN 6645:2000 (ISO
13878:1998)
6 TP mg/L Phương pháp so màu TCVN 8940:2011
3.3.2.4 Phương pháp thu mẫu và phân tích các chỉ tiêu về cá a. Phương pháp thu mẫu cá
Việc lựa chọn phương pháp khảo sát, thu mẫu cũng như xử lý mẫu được thực hiện theo “Hướng dẫn điều tra đa dạng sinh học cá” được ban hành kèm theo Công văn 2149/TCMTBTĐDSH ngày 14 tháng 9 năm 2016 của Tổng cục Môi trường.
Tuyến/điểm thu mẫu: Đặt ngư cụ, tổ chức đánh bắt cá theo cấp kênh đối với vùng
trồng Tràm, Keo lai và thu cá trên các kênh bao, gần vị trí các cống tiêu thốt nước ở vùng đệm cũng như vùng lõi.
- Dùng ngư cụ đánh bắt chủ động: Dùng lưới cước chặn hay khoanh vùng đoạn kênh 30 – 50 m sau đó dùng lưới kéo (lưới cước và lưới chày) để kéo cá trên tồn bộ diện tích kênh được khoanh vùng với sự hỗ trợ của ngư dân, ngư cụ này bắt được cá mồi, cá thương phẩm ở cả phân bố ở 3 tầng nước.
- Dùng ngư cụ đánh bắt thụ động: Dùng lưới quế (cở mắc 3 cm), lưới 3 màng và lưới 1 màng (4,5 cm) bắt được cá thương phẩm với nhiều loại kích cở khác nhau ở cả 3 tầng nước; ma trận (lú) được đặt cắt ngang toàn bộ bề mặt đáy kênh và chạy dọc đáy kênh, chủ yếu bắt cả tầng đáy với độ cao lưới khoảng 28 – 30 cm tính từ đáy kênh.
b. Phương pháp xử lý mẫu cá
- Đối với mẫu cá có kích thước lớn có thể định danh và xác định khối lượng được tại hiện trường sau đó thả cá về tự nhiên.
- Đối với mẫu cần phân tích lại: mẫu sau khi thu được rửa sạch bằng nước ngọt để mẫu được sạch, đồng thời loại bỏ các vi sinh vật bám vào mẫu. Sau đó mẫu được phân biệt theo từng loài khác nhau và mẫu thu được ghi nhãn để đánh dấu.
- Mẫu thu được giữ trong nước sạch có sục khí để giữ sống và giữ lạnh đối với cá thể chết để đưa về phân tích.
- Mẫu sau khi thu được phân loại sơ bộ và chụp ảnh tại thủy vực thu mẫu. Sau đó, cố định mẫu trong dung dịch formol 5 - 10% và đưa về phân tích tại phịng thí nghiệm đa dạng sinh học, Khoa Mơi trường & TNTN, Trường Đại học Cần Thơ.
Sau khi phân tích, mẫu được bảo quản trong cồn 70% để lưu trữ tại phịng thí nghiệm (Đa dạng sinh học) của Bộ môn Quản lý môi trường & TNTN, khoa Môi Trường & TNTN, trường Đại học Cần Thơ như là mẫu để đối chiếu.
c. Phương pháp định danh
Tất cả những mẫu vật thu thực tế được xử lý, định danh tên và sắp xếp hệ thống phân loại theo Fishbase (Froese & Pauly, 2019). Bên cạnh đó, việc định danh lồi cịn được dựa trên các tài liệu ở các nghiên cứu trước về thành phần lồi cá ở lưu vực sơng Mekong hay ở Đồng bằng sông Cửu Long như Khoa & Hương. (1993), Định và ctv. (2013) và Yên và ctv. (1992).
3.3.2.5 Phương pháp phỏng vấn nông hộ
Nghiên cứu tiến hành phỏng vấn ngẫu nhiên 90 hộ dân đại diện cho các mơ hình canh tác (Keo lai, Tràm và Lúa 2 vụ) để thu thập các thơng tin cần thiết liên quan. Cụ thể, có 30 hộ dân trên một mơ hình canh tác và được lặp lại trên 03 mơ hình. Biểu mẫu phiếu phỏng vấn với các nội dung chính như thơng tin chung nông hộ, hiện trạng canh tác, kĩ thuật canh tác, đa dạng cá và nguyên nhân biến động nguồn lợi cá.
3.3.3 Phương pháp xử lý số liệu 3.3.3.1 Phương pháp so sánh 3.3.3.1 Phương pháp so sánh
Các kết quả sau khi phân tích được lưu trữ trong phần mềm Microsoft Excel (Microsoft Crop., Washington, USA), điều này có thể tạo điều kiện thuận lợi cho việc phân tích thống kê. Dữ liệu chất lượng đất và nước được so sánh và phân tích theo cấp tuổi, tầng phèn và mùa do lượng mưa có thể dẫn đến sự thay đổi lớn trong quá trình nghiên cứu.
Đối với mẫu đất: Các chỉ tiêu sau khi phân tích được so sánh đối chiếu với các thang đánh giá chất lượng đất. Thang đánh giá này được trình bày chi tiết trong Chương 2. Tổng quan tài liệu.
Đối với mẫu nước: Các chỉ tiêu chất lượng môi trường nước sau khi phân tích được so sánh với quy chuẩn chất lượng mơi trường nước mặt QCVN 08-MT:2015/BTNMT cột A1 và A2 (Bộ Tài nguyên và Môi trường, 2015) đối với chất lượng nước sử dụng cho mục đích cấp nước sinh hoạt (sau khi áp dụng xử lý thông thường), bảo tồn động thực vật thủy sinh và các mục đích khác. Giá trị giới hạn của các chỉ tiêu phân tích được trình bày trong Bảng 3.3. Các chỉ tiêu chất lượng môi trường nước vượt giới hạn cho
phép ở cột A1, A2 được xem là khơng cịn phù hợp cho đời sống thủy sinh vật và có thể ảnh hưởng bất lợi đến cá tại khu vực nghiên cứu.
Bảng 3.3 Giá trị giới hạn của một số thông số đánh giá chất lượng môi trường nước STT Thông số Đơn vị Giá trị giới hạn A1 Giá trị giới hạn A2
1 pH - 6 - 8,5 6 - 8,5 2 EC mS/cm - - 3 N-NH4+ mg/L 0,3 0,3 4 Al3+ mg/L - - 5 Fe3+ mg/L 0,5 1 6 COD mg/L 10 15 7 BOD5 mg/L 4 6 8 DO mg/L ≥ 6 ≥ 5 9 N-NO3- mg/L 2 5
3.3.3.2 Phương pháp phân tích thống kê đa biến a. Phân tích sự khác biệt
Sự khác biệt trung bình về đặc tính mơi trường đất, nước và cá ở các mơ hình theo mùa (mùa mưa và mùa khơ), theo tầng phèn (phèn nông và phèn sâu), và theo cấp độ tuổi được đánh giá sử dụng phương pháp thống kê Independent Sample T-Test và phân tích phương sai một nhân tố (One - way ANOVA). Trong đó, thử nghiệm của Duncan được sử dụng để xác định sự khác biệt đáng kể về chất lượng đất và nước giữa các mơ hình trên cùng một tầng phèn hoặc trong cùng một mùa với độ tin cậy 95% (mức ý nghĩa p < 0,05). Bên cạnh đó, phân tích tương quan (Pearson) được sử dụng để mô tả mối quan hệ và sự phụ thuộc lẫn nhau giữa các chỉ tiêu chất lượng môi trường đất, nước và đa đạng cá tại khu vực nghiên cứu. Kết quả tương quan có thể cho biết mối liên hệ giữa các thơng số chất lượng đất, nước tại các mơ hình, tương quan của các thơng số môi trường đất với các thông số chất lượng môi trường nước. Các mối tương quan được coi là có ý nghĩa khi p < 0,05 và hệ số tương quan (r) chỉ ra mức độ tương quan của các thông số. Hệ số tương quan Pearson (r) thể hiện mối liên hệ giữa các biến, có giá trị trong khoảng từ -1 đến +1. Tương quan thuận là hai biến tăng hoặc giảm cùng nhau (r > 0); và tương quan nghịch là một biến giảm, một biến tăng (r < 0) (Gazzaz và ctv., 2012). Những mối tương quan được quan tâm khi giá trị tuyệt đối của hệ số tương quan r > 0,5 và mức ý nghĩa thống kê p < 0,05 (Feher và ctv., 2016). Tất cả phân tích trên được thực hiện bằng phần mềm thống kê IBM SPSS 20.0 Windows và Statgraphics Centurion XVI (Statgraphics Technologies,Virginia, USA). Cơng thức sau được sử dụng để tính tốn mối tương quan Pearson.
r = ∑n (Xi− X̅) × (Yi− Y̅) i=1 √∑n (Xi− X̅)2 i=1 × √∑n (Yi − Y̅)2 i=1 Trong đó:
r = hệ số tương quan Pearson giữa tham số X và tham số Y. n = số lần quan sát.
Xi = giá trị của X (cho lần quan sát thứ i). Yi = giá trị của Y (cho lần quan sát thứ i).
b. Phân tích cụm
Phân tích cụm (Cluster analysis, CA) được ứng dụng để phân nhóm nguồn nước theo không gian và thời gian (Feher et al., 2016; Chounlamany et al., 2017). Có tổng cộng 9 thơng số chất lượng nước và 6 thông số chất lượng đất được sử dụng trong phân tích cụm; các thơng số này được sử dụng tích hợp để xem xét tính tương đồng chất lượng đất và nước tại một vị trí nhất định so với một vị trí khác. Những vị trí thu mẫu có đặc tính ơ nhiễm tương đồng sẽ được nhóm vào cùng một nhóm, các đặc tính ơ nhiễm khác nhau sẽ được nhóm vào một nhóm khác. Việc phân tích cụm được tiến hành theo phương pháp của Ward (Salah et al., 2006) và được trình bày dưới dạng cây cấu trúc hay dendogram (Feher et al., 2016; Chounlamany et al., 2017). Các sơ đồ có thể giúp xác định được số lượng nhóm vị trí có đặc điểm tương đồng. Trong nghiên cứu này, phân tích cụm được sử dụng để gom nhóm những vị trí có chất lượng nước, đất ở các nghiệm thức tại khu vực nghiên cứu.
c. Phân tích thành phần chính
Tương tự phân tích cụm, các biến đầu vào của phân tích cũng được sử dụng tương tự, với 9 thông số chất lượng nước và 6 thơng số chất lượng đất. Phân tích nhân tố chính (Principal Component Analysis, PCA) dùng để rút trích thơng tin quan trọng từ bộ số liệu ban đầu (Feher et al., 2016; Chounlamany et al., 2017). Phương pháp xoay trục được sử dụng trong PCA là Varimax, mỗi biến số liệu ban đầu sẽ được xếp vào một thành phần chính (PC) và mỗi PC sẽ đại diện cho một nhóm nhỏ các biến ban đầu (Feher et al., 2016). Kết quả của PCA là giảm bớt những biến số liệu ban đầu khơng có đóng góp quan trọng vào sự biến động của số liệu. Tương quan giữa PCs và các biến số liệu ban đầu được thể hiện thông qua hệ số tương quan gia trọng (loading) (Feher et al., 2016). Trị tuyệt đối của hệ số tương quan gia trọng lớn hơn 0.75 có nghĩa là mối tương quan chặt giữa thành phần chính và thơng số chất lượng nước. từ 0.75-0.50 là tương quan trung bình và 0.5-0.3 là tương quan yếu. Kết quả phân tích thành phần chính giúp xác định yếu tố chính gây ảnh hưởng đến chất lượng đất, nước và dự đoán các ảnh hưởng đến cá tại khu vực nghiên cứu. Phân tích CA và PCA được tiến hành bằng cách sử dụng
phần mềm Statgraphics Centurion XVI (Statgraphics Technologies,Virginia, USA). Việc tính tốn PC có thể được dựa theo cơng thức sau:
Zij = ai1x1j + ai2x2j + ai3x3j + … … … . . + aimxmj
Trong đó z là của điểm thành phần, a là tải thành phần, x biểu thị giá trị đo được của biến, i là số thành phần, j là số mẫu và m là tổng số biến.
d. Phân tích biệt số
Phương pháp phân tích sự khác biệt (Discriminant Analysis - DA) là kỹ thuật phân tích dữ liệu khi biến phụ thuộc là biến phân loại và biến độc lập là biến định lượng. Giúp xác định biến độc lập nào là nguyên nhân chính gây ra sự khác biệt giữa các nhóm. Chính vì vậy, nghiên cứu đã sử dụng phân tích biệt số đã được áp dụng nhằm xác định các thơng số có ý nghĩa và quan trọng nhất góp phần vào sự biến đổi chất lượng nước và đất giữa hai mùa và hai tầng phèn trong nghiên cứu (Samsudin et al., 2019). Trong nghiên cứu này, biến phụ thuộc là mùa và tầng phèn; trong khi, biến độc lập là các chỉ tiêu chất lượng nước (9 thông số) và đất (6 thông số). Từ đó, tìm ra thơng số chính và có ý nghĩa dẫn đế sự sự khác biệt nhất giữa hai mùa mưa và mùa khô, tầng phèn nông và phèn sâu (Samsudin et al., 2019). Phân tích DA được thực hiện bằng phần mềm thống kê mềm Statgraphics Centurion XVI (Statgraphics Technologies,Virginia, USA).
e. Phân tích mối liên hệ mơi trường và đa dạng cá
Phân tích BIO-ENV: Cho phép so sánh các ma trận khoảng cách hoặc đồng dạng giữa hai bộ dữ liệu có chung mẫu hoặc biến. Thường trong việc khám phá mối tương quan giữa các biến môi trường và các chỉ tiêu sinh học. Trong trường hợp này, ma trận đồng dạng của số liệu sinh học là cố định, trong khi các tập hợp con của các biến môi trường được sử dụng để tính tốn ma trận đồng dạng mơi trường. Sau đó, một hệ số tương quan (hệ số tương quan Spearman) được tính giữa hai ma trận và tập hợp con tốt nhất (BEST) của các biến mơi trường sau đó có thể được xác định và tiếp tục kiểm tra hoán vị để xác định tầm quan trọng của các biến môi trường (nước) đối với các chỉ tiêu sinh học. Cụ thể, 9 thông số chất lượng nước và cấu trúc thành phần loài cá đã được sử dụng để thực hiện phân tích trong nghiên cứu hiện tại. Phân tích BIO-ENV được tiến hành bằng cách sử dụng phần mềm Primer 5.2 for Windows (PRIMER-E Ltd, Plymouth, UK).
CHƯƠNG 4
KẾT QUẢ VÀ THẢO LUẬN
4.1 Khảo sát cấu trúc, kỹ thuật trồng rừng Keo lai, Tràm và lúa 2 vụ tác động của chúng đến môi trường đất, nước và cá tự nhiên ở vùng nghiên cứu chúng đến môi trường đất, nước và cá tự nhiên ở vùng nghiên cứu
4.1.1 Khảo sát cấu trúc, kỹ thuật lên liếp trồng Keo lai và Tràm
Các thông số liên quan đến việc thiết kế và cải tạo khu vực trồng Keo lai và Tràm được trình bày trong Bảng 4.1 và Bảng 4.2. Kết quả cho thấy khu vực trồng Tràm tại VQG U Minh Hạ có chiều dài liếp nhỏ nhất 300 m, lớn nhất 1000 m; trung bình 827,3 m. Chiều rộng liếp nhỏ nhất được ghi nhận khoảng 4 m, lớn nhất 18 m; trung bình 11,3 m; chiều cao liếp thấp nhất 0,3 m, lớn nhất 1,5 m; trung bình 0,6 m. Chiều dài mương có kích thước trung bình bằng với chiều dài liếp, trung bình khoảng 834 m. Trong khi đó, chiều rộng mương nhỏ nhất 0,6 m, lớn nhất 15 m; trung bình 3,6 m. Theo báo cáo của Xuân (1993) chiều rộng mặt mương giữa các liếp canh tác khoảng 4 – 5 m đối với mơ hình trồng Tràm tại vùng Đồng bằng sơng Cửu Long; điều này có thể thấy kỹ thuật trồng có sự chênh lệch đáng kể.
Bảng 4.1 Kích thước mặt liếp và kênh mương ở mơ hình Tràm
Giá trị
Kích thước mặt liếp Kích thước kênh/mương Chiều dài (m) Chiều rộng (m) Chiều cao (m) Chiều dài (m) Chiều rộng (m) Nhỏ nhất 300 4 0,3 300 0,6 Lớn nhất 1000 18 1,5 1000 15 Trung bình 827,3 11,3 0,6 834 3,6 Trung vị 955 12 0,6 955 3
Đối với khu vực Keo lai khu vực trồng Keo lai có kích thước chiều dài liếp nhỏ nhất 100 m, lớn nhất 1000 m; trung bình 717 m; chiều rộng liếp nhỏ nhất 3 m, lớn nhất 17m; trung bình 10,1m; chiều cao liếp thấp nhất 0,4 m, lớn nhất 1,5 m; trung bình 0,9 m. Tương tự với kỹ thuật lên liếp trồng Tràm thì kích thước chiều dài liếp và chiều dài kênh/mương tương tự nhau, đạt giá trị trung bình 717 m. Ở mơ hình trồng Keo lai, chiều rộng mương nhỏ nhất 2,5 m, lớn nhất 15 m; trung bình 4,6 m (Bảng 4.2). Tuy nhiên, kết quả khảo sát đã ghi nhận có sự chênh lệch chiều cao liếp so với đề xuất của Trung tâm Giống nông nghiệp Cà Mau lớn hơn 0,3 m đối với Tràm và 1 m đối với trồng Keo lai. Yếu tố độ cao địa hình tại các khu vực trồng được xem là nguyên nhân dẫn đến sự chênh lệch này.