Tổng quan nghiên cứu

Tình trạng phú dưỡng tại các thủy vực nội địa đô thị đang là thách thức môi trường cấp bách tại Thủ đô Hà Nội. Quận Hoàng Mai nằm ở vùng trũng phía Nam thành phố với tổng diện tích tự nhiên 4.032 ha, sở hữu mạng lưới mặt nước rộng trên 400 ha bao gồm 14 hồ và đầm trũng lớn nhỏ. Trong đó, ba hồ giữ vai trò điều hòa trọng yếu nhất là hồ Yên Sở với diện tích 130 ha, hồ Linh Đàm với 75 ha và hồ Định Công với 25 ha. Do tiếp nhận trực tiếp lượng lớn nước thải sinh hoạt chưa qua xử lý từ các sông thoát nước chính như sông Tô Lịch, sông Sét, sông Lừ và sông Kim Ngưu (với chỉ số nhu cầu oxy sinh học BOD5 vượt quy chuẩn từ 6 đến 12 lần), hệ sinh thái các hồ này đang suy thoái nghiêm trọng, xuất hiện hiện tượng nở hoa tảo và cá chết hàng loạt.

Mục tiêu cốt lõi của nghiên cứu là xác định mối tương quan giữa chỉ số trạng thái phú dưỡng (TSI) và đặc tính phổ phản xạ quang học từ dữ liệu vệ tinh Landsat 8 OLI, từ đó xây dựng mô hình viễn thám bán thực nghiệm để định lượng, lập bản đồ và mô phỏng biến động không gian - thời gian của hiện tượng phú dưỡng tại ba hồ trọng điểm giai đoạn 2013 - 2017. Kết quả phân tích hiện trường ghi nhận hàm lượng Chlorophyll-a trung bình đạt 100,14 µg/L (biến thiên từ 28,5 đến 304,3 µg/L) và chỉ số TSI trung bình ở mức 76 (dao động từ 63 đến 87), đưa các hồ vào ngưỡng từ phú dưỡng đến siêu phú dưỡng. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp giải pháp quan trắc diện rộng với độ tin cậy đạt hệ số xác định r2 = 0,86, giúp các cơ quan quản lý cắt giảm chi phí và thời gian giám sát chất lượng nước mặt so với phương pháp đo đạc truyền thố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 được xây dựng trên nền tảng lý thuyết trạng thái phú dưỡng của Carlson (1977) kết hợp với bảng phân loại mở rộng của Carlson và Simpson (1996). Chỉ số TSI phân chia mức độ dinh dưỡng của vực nước thành bốn cấp độ: nghèo dinh dưỡng (TSI dưới 30, Chl-a dưới 0,95 µg/L), dinh dưỡng trung bình (TSI từ 30 đến 50, Chl-a từ 0,95 đến 7,3 µg/L), phú dưỡng (TSI từ 50 đến 70, Chl-a từ 7,3 đến 56 µg/L) và siêu phú dưỡng (TSI từ 70 đến trên 80, Chl-a từ 56 đến trên 155 µg/L). Trong đó, Chlorophyll-a (Chl-a) được xác định là chỉ thị sơ cấp trực tiếp phản ánh sinh khối thực vật nổi và sức khỏe sinh thái hồ.

Cơ sở quang học viễn thám nước nội địa dựa trên quy luật tương tác bức xạ điện từ của các hợp phần quang học hoạt tính: Chl-a hấp thụ mạnh bức xạ mặt trời ở dải sóng xanh lam (450 - 475 nm) và dải sóng đỏ (670 nm), đồng thời phản xạ mạnh tại dải sóng xanh lục (550 - 560 nm) và cận hồng ngoại (gần 700 nm). Việc sử dụng tỷ số giữa các kênh phổ nhìn thấy giúp triệt tiêu nhiễu xạ khí quyển và phản xạ bề mặt, tạo tiền đề vững chắc cho việc định lượng Chl-a và TSI.

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

Nguồn dữ liệu sử dụng bao gồm 12 cảnh ảnh vệ tinh Landsat 8 OLI (Path/Row: 127/45, múi chiếu WGS84 UTM 48N, độ phân giải không gian 30 m x 30 m, độ phủ mây dưới 20%) thu thập từ Cục Khảo sát Địa chất Hoa Kỳ giai đoạn từ tháng 6/2013 đến tháng 8/2017. Quy trình xử lý ảnh trên phần mềm ENVI 5.3 bao gồm: chuyển đổi giá trị số (DN) sang bức xạ phổ TOA, hiệu chỉnh khí quyển bằng phương pháp loại trừ điểm đen (Dark Subtraction) căn chỉnh theo dữ liệu phổ thực địa, và chiết tách mặt nước tự động thông qua tỷ số kênh SWIR và Green (kênh 6 và kênh 3).

Để xây dựng và kiểm chứng mô hình, tác giả thực hiện 6 đợt khảo sát thực địa đồng bộ với thời gian vệ tinh quét qua (ngày 01/06/2016, 28/09/2016, 01/04/2017, 16/04/2017, 04/06/2017, 07/08/2017). Tổng cỡ mẫu thu thập là 48 mẫu nước mặt tại tầng sâu 0 - 25 cm (phù hợp với độ sâu đĩa Secchi luôn dưới 25 cm), phân bố tại hồ Linh Đàm (18 mẫu), hồ Yên Sở (17 mẫu) và hồ Định Công (13 mẫu). Phương pháp chọn mẫu áp dụng nguyên tắc phân tầng không gian từ khu vực nước trong nhất, khu vực giữa hồ đến các điểm xả thải ven bờ có độ đục cao. Dữ liệu phổ phản xạ mặt nước được đo trực tiếp bằng máy quang phổ hiện trường GER 1500 (dải bước sóng 350 - 1050 nm, độ phân giải phổ 1,5 nm). Hàm lượng Chl-a được phân tích trong phòng thí nghiệm bằng phương pháp chiết xuất acetone 90% đo quang phổ trên máy Hach DR 5000 theo chuẩn APHA. Dữ liệu được phân tích tương quan Pearson và hồi quy đơn biến trên phần mềm IBM SPSS Statistics 20 ở độ tin cậy 95%.

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

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

Thứ nhất, kết quả phân tích 48 mẫu thực địa khẳng định cả ba hồ đều ở trạng thái phú dưỡng nghiêm trọng đến siêu phú dưỡng. Hàm lượng Chl-a trung bình tại các đợt khảo sát đạt 100,14 µg/L, giá trị cực đại lên tới 304,3 µg/L tại hồ Linh Đàm vào ngày 01/06/2016. Chỉ số TSI trung bình đạt mức 76, trong đó hồ Linh Đàm có mức độ ô nhiễm cao nhất với TSI trung bình dao động từ 73 đến 83, hồ Định Công đạt 67 đến 74 và hồ Yên Sở đạt 68 đến 73.

Thứ hai, nghiên cứu xác lập được mối tương quan tuyến tính chặt chẽ giữa chỉ số TSI và tỷ số kênh phổ b3/b2 (kênh xanh lục 561 nm / kênh xanh lam 482 nm) của ảnh Landsat 8 OLI với hệ số tương quan Pearson R = 0,90 và hệ số xác định r2 = 0,76 (sai số chuẩn 2,2). Phương trình toán học chiết tách có dạng hàm tuyến tính: TSI = 50,07 * (b3/b2) - 0,43. Kiểm chứng độc lập giữa TSI tính từ ảnh và TSI thực nghiệm đạt hệ số xác định r2 = 0,86 với sai số trung bình tuyệt đối chỉ 1,1 (tương đương sai lệch 1,5% so với giá trị đo đạc trực tiếp).

Thứ ba, diễn biến theo thời gian từ năm 2013 đến năm 2017 cho thấy chỉ số TSI có xu hướng gia tăng liên tục qua các năm. Điển hình tại hồ Linh Đàm vào tháng 6, giá trị TSI trung bình tăng từ 52,46 năm 2013 lên 56,81 năm 2014, đạt 57,79 năm 2016 và chạm ngưỡng 60,15 năm 2017 (tăng 14,6%). Tại hồ Yên Sở, TSI tăng từ 50,36 (2013) lên 58,95 (2017), tương đương mức tăng 17,1%. Trong cùng một năm, các tháng mùa khô cuối năm (tháng 10 - 12) luôn có TSI cao hơn từ 1,5 đến 3,5 đơn vị so với các tháng đầu hè do lưu lượng nước bổ cập giảm và nồng độ chất ô nhiễm tích tụ.

Thứ tư, về mặt không gian, chỉ số TSI phân bố không đồng đều trên mặt hồ. Các giá trị siêu phú dưỡng (TSI từ 70 đến trên 80) tập trung đậm đặc tại khu vực ven bờ, xung quanh các cửa cống xả thải ngầm từ khu dân cư (như khu đô thị bán đảo Linh Đàm, khu dân cư Định Công) và giảm dần về phía tâm hồ với giá trị dao động từ 50 đến 60.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự gia tăng chỉ số TSI xuất phát từ áp lực đô thị hóa và gia tăng dân số cơ học nhanh chóng tại quận Hoàng Mai. Thống kê giai đoạn 2011 - 2016 cho thấy dân số toàn quận tăng từ 344.785 người lên 379.483 người, với tỷ lệ tăng dân số cơ học năm 2016 đạt đỉnh 10,3%. Diện tích đất xây dựng đô thị mở rộng nhanh làm thu hẹp không gian tự nhiên, trong khi hạ tầng xử lý nước thải sinh hoạt không theo kịp, dẫn đến hàng chục nghìn mét khối nước thải giàu Nitơ và Photpho đổ thẳng vào các hồ.

Khi biểu diễn trực quan qua các bản đồ chuyên đề phân mảnh mật độ (density slicing) trong ENVI và ArcGIS, không gian phú dưỡng được minh họa rõ nét bằng các thang màu chuyển đổi từ vàng nhạt (TSI 50 - 60) ở giữa hồ sang đỏ đậm (TSI 70 - 80) tại các mép bờ tiếp giáp cống xả. Biểu đồ phân tán (scatter plot) giữa tỷ số b3/b2 và TSI thực địa thể hiện sự tập trung cao của các điểm quan sát quanh đường hồi quy tuyến tính, khẳng định độ vững của mô hình. Kết quả này hoàn toàn tương đồng với các nghiên cứu quốc tế của Lim và Choi (2015) khi ứng dụng Landsat 8 đánh giá hồ nội địa Hàn Quốc, cũng như nghiên cứu của nhóm tác giả tại Hồ Tây và hồ Hoàn Kiếm (2016 - 2017), chứng minh tỷ số Green/Blue là chỉ thị quang học tối ưu nhất cho các vực nước đô thị giàu dinh dưỡng.

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

Nhằm kiểm soát hiện tượng phú dưỡng và nâng cao hiệu quả quản lý chất lượng nước các hồ đô thị, nghiên cứu đề xuất 4 nhóm giải pháp trọng tâm:

  1. Nâng cấp và tách biệt hoàn toàn hệ thống thu gom nước thải sinh hoạt: Tiến hành xây dựng hệ thống cống bao quanh hồ Linh Đàm và Định Công để thu gom 100% nước thải sinh hoạt về Nhà máy xử lý nước thải Yên Sở trước khi xả ra môi trường. Mục tiêu giảm 45 - 50% tải lượng Nitơ và Photpho tích tụ trong nước hồ; lộ trình thực hiện từ năm 2024 đến năm 2026; cơ quan chủ trì là UBND quận Hoàng Mai phối hợp với Sở Xây dựng Hà Nội.
  2. Thiết lập quy trình giám sát chất lượng nước định kỳ bằng công nghệ viễn thám: Tích hợp ảnh vệ tinh Landsat 8/9 và Sentinel-2 vào quy trình quan trắc chuẩn với tần suất cập nhật 16 ngày/lần. Mục tiêu cắt giảm 60% chi phí và nhân lực quan trắc thủ công, đồng thời cảnh báo sớm hiện tượng bùng phát tảo nở hoa; triển khai từ quý I/2025; đơn vị thực hiện là Chi cục Bảo vệ Môi trường - Sở Tài nguyên và Môi trường Hà Nội.
  3. Cải tạo thủy vực, nạo vét trầm tích và phát triển vùng đệm sinh thái: Định kỳ nạo vét lớp bùn đáy tích tụ kim loại nặng và chất hữu cơ yếm khí, kết hợp trồng các bè thủy sinh nổi (cây ngổ trâu, thủy trúc) chiếm 10 - 15% diện tích mặt hồ nhằm hấp thụ dinh dưỡng dư thừa. Duy trì độ sâu trung bình đáy hồ trên 2,0 m; kế hoạch thực hiện hàng năm giai đoạn 2024 - 2027; do Ban Quản lý Dự án Đầu tư Xây dựng quận Hoàng Mai đảm trách.
  4. Chuẩn hóa cơ sở dữ liệu GIS và hoàn thiện quy chế pháp lý bảo vệ hành lang hồ: Ban hành quy định xử phạt nghiêm khắc các cơ sở kinh doanh, hộ dân xả rác và nước thải trực tiếp xuống lòng hồ; số hóa 100% hồ điều hòa lên hệ thống cơ sở dữ liệu GIS dùng chung phục vụ điều hành đô thị thông minh; hoàn thành trước tháng 12/2025; cơ quan quản lý là UBND thành phố Hà Nội và Sở Thông tin và Truyền thông.

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

Nội dung và phương pháp tiếp cận của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước về môi trường và đô thị: Chi cục Bảo vệ Môi trường, Sở Xây dựng và UBND các quận huyện có thể ứng dụng trực tiếp quy trình viễn thám để theo dõi diễn biến tài nguyên nước mặt trên phạm vi 4.032 ha, phục vụ công tác quy hoạch thoát nước và phê duyệt dự án phát triển đô thị.
  2. Các nhà khoa học, giảng viên và học viên cao học: Nhóm nghiên cứu thuộc các chuyên ngành Địa chất Môi trường, Quản lý Tài nguyên, Viễn thám và Hệ thông tin địa lý (GIS) có thể khai thác bộ dữ liệu thực nghiệm 48 mẫu nước và phương pháp giải đoán phổ mặt nước làm tài liệu tham khảo cho các công trình nghiên cứu mở rộng.
  3. Các doanh nghiệp tư vấn môi trường và công nghệ xử lý nước: Đơn vị thiết kế công trình thủy lợi, trạm xử lý nước thải đô thị có thể sử dụng dữ liệu biến động TSI giai đoạn 2013 - 2017 để tính toán tải lượng ô nhiễm thiết kế, đánh giá hiệu quả sau xử lý của các giải pháp cải tạo hồ.
  4. Các tổ chức phi chính phủ và cộng đồng dân cư địa phương: Các tổ chức xã hội hoạt động trong lĩnh vực bảo vệ nguồn nước có cơ sở dữ liệu tin cậy để nâng cao nhận thức cộng đồng, thúc đẩy các sáng kiến bảo vệ hành lang xanh và giám sát xã hội đối với các nguồn xả thải.

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

  1. Tại sao tỷ số kênh 3 trên kênh 2 (Green/Blue) của Landsat 8 lại tối ưu nhất để ước tính chỉ số TSI? Tỷ số b3/b2 phản ánh chính xác đặc tính quang học của diệp lục: Chl-a hấp thụ mạnh ánh sáng xanh lam (kênh 2: 482 nm) và phản xạ cực đại ở ánh sáng xanh lục (kênh 3: 561 nm). Nghiên cứu thực nghiệm chứng minh tỷ số này đạt hệ số tương quan cao nhất R = 0,90 và r2 = 0,76, vượt trội hơn hẳn các kênh đơn lẻ.

  2. Mức độ phú dưỡng của các hồ tại quận Hoàng Mai giai đoạn 2013 - 2017 diễn biến thế nào? Cả ba hồ đều ở tình trạng phú dưỡng đến siêu phú dưỡng với TSI trung bình đạt 76 và Chl-a trung bình đạt 100,14 µg/L. Xu hướng ô nhiễm tăng dần theo thời gian, điển hình hồ Linh Đàm tăng từ mức TSI 52,46 (tháng 6/2013) lên 60,15 (tháng 6/2017), đe dọa trực tiếp đến hệ sinh thái thủy vực.

  3. Ưu điểm nổi bật của việc sử dụng ảnh Landsat 8 so với phương pháp đo đạc truyền thống là gì? Ảnh Landsat 8 cung cấp góc nhìn không gian diện rộng bao quát toàn bộ diện tích mặt nước, có tính kế thừa lịch sử đa thời gian với chu kỳ 16 ngày và hoàn toàn miễn phí. Phương pháp này giúp giảm trên 60% kinh phí khảo sát thực địa mà vẫn đảm bảo độ tin cậy r2 = 0,86.

  4. Nguyên nhân chính dẫn đến sự phân bố TSI không đều giữa ven bờ và tâm hồ là gì? Sự chênh lệch bắt nguồn từ các điểm xả thải cục bộ ven bờ. Tại các cửa cống ngầm tiếp nhận nước thải sinh hoạt chưa xử lý, nồng độ Nitơ và Photpho đậm đặc khiến Chl-a bùng phát đạt TSI từ 70 đến 87, sau đó lan tỏa và pha loãng dần ra khu vực giữa hồ (TSI 50 - 60).

  5. Thuật toán bán thực nghiệm xây dựng trong luận văn có thể áp dụng cho các hồ khác tại Việt Nam không? Hoàn toàn có thể áp dụng cho các hồ đô thị nội địa có đặc tính quang học tương đồng, hàm lượng chất rắn lơ lửng và độ sâu đĩa Secchi dưới 25 cm. Tuy nhiên, cần tiến hành đo đạc phổ thực địa để hiệu chỉnh hệ số phương trình hồi quy phù hợp với nền khí quyển địa phương.

Kết luận

  • Luận văn đã làm sáng tỏ hiện trạng ô nhiễm phú dưỡng nghiêm trọng tại ba hồ trọng điểm quận Hoàng Mai (Linh Đàm, Yên Sở, Định Công) với hàm lượng Chl-a trung bình 100,14 µg/L và TSI trung bình 76.
  • Xây dựng thành công mô hình toán học ước tính chỉ số TSI từ tỷ số kênh phổ b3/b2 của ảnh vệ tinh Landsat 8 OLI đạt độ chính xác cao với hệ số xác định r2 = 0,86 và sai số kiểm chứng chỉ 1,1.
  • Phân tích chi tiết quy luật biến động không gian và thời gian của TSI từ 12 cảnh ảnh vệ tinh giai đoạn 2013 - 2017, khẳng định mối liên hệ trực tiếp giữa áp lực gia tăng dân số cơ học (10,3%) với tốc độ suy thoái nguồn nước.
  • Cung cấp cơ sở khoa học vững chắc giúp chuyển đổi phương thức quản lý chất lượng nước từ quan trắc thủ công gián đoạn sang công nghệ viễn thám bán tự động liên tục, tiết kiệm ngân sách nhà nước.
  • Đề xuất hệ thống giải pháp toàn diện bao gồm tách dòng thu gom nước thải, nạo vét đáy hồ định kỳ và số hóa GIS phục vụ phát triển đô thị bền vững.

Kế hoạch hành động tiếp theo yêu cầu các cơ quan chức năng hoàn thiện tích hợp thuật toán viễn thám vào cổng thông tin quản lý môi trường đô thị trong giai đoạn 2024 - 2025. Hãy chủ động ứng dụng công nghệ viễn thám và dữ liệu không gian ngay hôm nay để bảo vệ và phục hồi các hồ đô thị, xây dựng Thủ đô Hà Nội xanh, thông minh và phát triển bền vững.