Tổng quan nghiên cứu

Sóng nhiệt biển đang trở thành một trong những hiểm họa sinh thái nghiêm trọng nhất đối với đa dạng sinh học đại dương trong bối cảnh biến đổi khí hậu toàn cầu. Tại khu vực biển Đông Nam Á, các dữ liệu vệ sinh học ghi nhận tần suất và thời gian kéo dài của các đợt sóng nhiệt biển đã tăng từ 5 đến 9 ngày mỗi thập kỷ trong giai đoạn từ năm 1982 đến 2018, với số ngày có nhiệt độ bề mặt nước biển vượt ngưỡng 34°C dao động từ 40 đến 60 ngày mỗi năm. Hiện tượng gia tăng nhiệt độ cực đoan này thường vượt quá giới hạn chống chịu sinh lý tối ưu của các loài sinh vật phù du ven bờ, vốn là mắt xích nền tảng duy trì năng suất sơ cấp và thứ cấp của toàn bộ hệ sinh thái ven biển nhiệt đới.

Nghiên cứu của luận văn tập trung giải quyết bài toán cấp thiết: Đánh giá tác động độc lập và tương tác giữa sóng nhiệt biển cực đoan cùng áp lực săn mồi phi tiêu thụ từ ấu trùng cá lên sự sinh trưởng, phát triển và tính mềm dẻo xuyên thế hệ của loài giáp xác chân mái chèo Pseudodiaptomus incisus. Đây là loài sinh vật phù du chiếm ưu thế tuyệt đối tại các đầm phá ven biển miền Trung Việt Nam, đóng góp từ 97% đến 98% khẩu phần thức ăn tự nhiên của ấu trùng nhiều loài cá kinh tế giá trị cao như cá chẽm và cá bớp. Mục tiêu cụ thể của đề tài là phân tích tỷ lệ sống, kích thước bọc trứng, tỷ lệ nở, ấu trùng nauplii tích lũy và lượng viên phân thải ra qua hai thế hệ liên tiếp F1 và F2.

Toàn bộ quá trình thực nghiệm được triển khai nghiêm ngặt tại Phòng thí nghiệm Nghiên cứu Nuôi trồng Thủy sản, Trung tâm Nghiên cứu Thủy sản Nhiệt đới Cam Ranh thuộc Đại học Nha Trang từ ngày 27/07/2020 đến ngày 08/09/2020, kết hợp xử lý và phân tích thống kê chuyên sâu tại Khoa Sinh học, Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội từ tháng 09/2021 đến tháng 11/2023. Kết quả nghiên cứu cung cấp hệ thống dữ liệu định lượng quan trọng, giúp các nhà khoa học và nhà quản lý môi trường biển dự báo chính xác mức độ suy giảm năng suất sinh học của mạng lưới thức ăn thủy sinh dưới các kịch bản biến đổi khí hậu khắc nghiệt.

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

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng của ba khung lý thuyết sinh thái học tiến hóa và sinh lý học hiện đại:

Thứ nhất là Lý thuyết Tính mềm dẻo xuyên thế hệ (Transgenerational Plasticity - TGP). Lý thuyết này giải thích cơ chế mà điều kiện môi trường bất lợi do thế hệ bố mẹ trải qua có thể điều chỉnh và nâng cao khả năng thích ứng của thế hệ con cái trong cùng điều kiện đó thông qua các biến đổi biểu sinh phi di truyền (epigenetics), chẳng hạn như quá trình methyl hóa gen hoặc tối ưu hóa nguồn dưỡng chất noãn hoàng truyền lại từ mẹ.

Thứ hai là Khung lý thuyết Tương tác đa áp lực sinh thái (Multiple Stressors Framework). Khung lý thuyết này phân tích tác động đồng thời của áp lực phi sinh học (nhiệt độ cực đoan) và áp lực sinh học (mùi hóa học từ động vật ăn thịt), phân định rõ các dạng phản ứng theo cơ chế cộng gộp (additive), đối kháng (antagonistic) hoặc hiệp đồng (synergistic).

Thứ ba là Lý thuyết Ngân sách năng lượng sinh học động (Dynamic Energy Budget), chỉ ra sự đánh đổi năng lượng nghiêm trọng giữa việc duy trì sự sống cơ bản và các chức năng sinh sản khi sinh vật đối mặt với tình trạng sốc nhiệt kéo dài.

Nghiên cứu vận dụng 4 khái niệm cốt lõi:

  • Sóng nhiệt biển (MHW): Thời kỳ nhiệt độ nước biển tăng cao bất thường vượt quá bách phân vị thứ 90 của chuỗi số liệu 30 năm và kéo dài trên 5 ngày liên tục.
  • Tín hiệu săn mồi từ cá (FPC / Kairomones): Các hợp chất hóa học do động vật săn mồi tiết ra môi trường nước, tạo ra áp lực đe dọa phi tiêu thụ đối với con mồi.
  • Hiệu ứng bố mẹ (Parental Effects): Sự truyền tải trực tiếp tài nguyên, chất lượng noãn hoàng từ thế hệ F1 sang thế hệ F2.
  • Chỉ số viên phân (Faecal Pellets): Chỉ số phản ánh tốc độ bắt mồi, tiêu hóa và chuyển hóa vật chất của giáp xác nổi.

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

Nguồn mẫu sinh vật khởi đầu (thế hệ F0) gồm các cá thể Pseudodiaptomus incisus trưởng thành mang bọc trứng được thu thập tự nhiên tại các đầm ven biển Vịnh Cam Ranh bằng lưới vớt động vật nổi chuyên dụng có kích thước mắt lưới 200 µm, ở điều kiện nhiệt độ môi trường 27°C đến 28°C và độ mặn 35 PSU. Cỡ mẫu nghiên cứu được thiết kế với độ lặp lại cao: thế hệ F1 bố trí 4 nghiệm thức trực giao (2 mức nhiệt độ 30°C và 34°C × 2 mức tín hiệu cá săn mồi gồm có hoặc không có FPC) với 10 bình lặp lại cho mỗi nghiệm thức, mật độ 180 đến 240 ấu trùng nauplii trong mỗi bình 1,2 lít. Thế hệ F2 mở rộng thành 16 nghiệm thức trực giao kết hợp toàn diện giữa lịch sử phơi nhiễm của F1 và môi trường sống của F2, mỗi nghiệm thức gồm 10 lần lặp lại độc lập. Thức ăn tiêu chuẩn là vi tảo Isochrysis galbana được duy trì ổn định ở mật độ từ 30.000 đến 33.000 tế bào/mL (tương đương 800 đến 850 µg Carbon/lít).

Phương pháp phân tích dữ liệu được thực hiện trên môi trường phần mềm R phiên bản 4.x với các mô hình tuyến tính tổng quát (GLM) và mô hình hỗn hợp (GLMER). Tỷ lệ sống của cá thể đực và cái được phân tích qua mô hình hồi quy Quasibinomial với kiểm định F nhằm kiểm soát triệt để hiện tượng phân tán quá mức (overdispersion). Tỷ lệ nở của trứng được mô hình hóa theo phân phối Binomial với kiểm định Chi-square. Kích thước bọc trứng và số lượng nauplii nở trên mỗi bọc được phân tích bằng mô hình hồi quy Poisson và mô hình nhị thức âm (Negative Binomial GLM). Các chỉ số tích lũy nauplii và tổng lượng viên phân sau 5 ngày được đánh giá qua mô hình Quasipoisson. Việc lựa chọn các mô hình phi tuyến tính này là bắt buộc do dữ liệu đếm và tỷ lệ sinh học không tuân theo phân phối chuẩn, giúp đảm bảo giá trị AIC nhỏ nhất và hạn chế tối đa sai số dự báo.

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

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

Phân tích số liệu thực nghiệm qua hai thế hệ giáp xác đã chứng minh 4 phát hiện then chốt:

Thứ nhất, ở thế hệ F1, sóng nhiệt biển 34°C trực tiếp gây tổn thương nghiêm trọng đến khả năng duy trì quần thể của P. incisus. Tỷ lệ sống của cả cá thể đực và cái đều sụt giảm từ 17% đến 18% so với điều kiện đối chứng 30°C (p < 0,001). Số lượng ấu trùng nauplii nở thành công từ mỗi bọc trứng sụt giảm kỷ lục tới 62%, trong khi tổng lượng nauplii tích lũy sau 5 ngày giảm 27% và lượng viên phân thải ra giảm 28%.

Thứ hai, tác động của tín hiệu cá săn mồi (FPC) phụ thuộc chặt chẽ vào nền nhiệt độ. Tại nhiệt độ tối ưu 30°C, FPC kích thích phản ứng phòng vệ thích nghi giúp giáp xác tăng khoảng 15% lượng nauplii tích lũy và tăng số lượng viên phân thải ra. Tuy nhiên, khi nhiệt độ tăng lên 34°C, tác động kích thích sinh sản của FPC hoàn toàn biến mất do sinh vật bị giới hạn nghiêm trọng về mặt năng lượng duy trì.

Thứ ba, ở thế hệ F2, nghiên cứu phát hiện bằng chứng rõ nét về tính mềm dẻo xuyên thế hệ (TGP) đối với tỷ lệ sống. Các cá thể F2 có bố mẹ F1 từng trải qua sóng nhiệt biển (nghiệm thức F1-MHW × F2-MHW) đạt tỷ lệ sống tương đương với nhóm đối chứng sống ở 30°C, cao hơn khoảng 10% so với nhóm F2 lần đầu tiên tiếp xúc với sóng nhiệt biển.

Thứ tư, sóng nhiệt biển từ thế hệ bố mẹ tạo ra hiệu ứng tiêu cực kéo dài (Parental effect). Nhóm F2 sống ở nhiệt độ đối chứng 30°C nhưng có bố mẹ từng chịu sốc nhiệt ở F1 bị suy giảm tới 22% kích thước bọc trứng và giảm 33% số nauplii nở ra từ mỗi bọc.

Thảo luận kết quả

Cơ chế sinh học giải thích cho sự suy giảm nghiêm trọng ở thế hệ F1 bắt nguồn từ tình trạng mất cân bằng năng lượng tế bào. Khi nhiệt độ vượt ngưỡng tối ưu lên 34°C, nhu cầu oxy tế bào tăng vọt vượt quá khả năng cung cấp của hệ tuần hoàn mở, đồng thời ức chế quá trình tổng hợp phân tử ATP và gây biến tính các cấu trúc đại phân tử. Lượng viên phân thải ra giảm 28% đến 32% phản ánh rõ rệt sự suy thoái trong hoạt động lọc mồi và tiêu hóa vi tảo. Mối tương quan thuận chặt chẽ giữa lượng viên phân và số lượng nauplii tích lũy khẳng định rằng sự thiếu hụt năng lượng nạp vào đã trực tiếp tước đoạt nguồn năng lượng phân bổ cho việc tạo trứng và nuôi dưỡng phôi.

Sự phục hồi hoàn toàn tỷ lệ sống ở thế hệ F2 thông qua tính mềm dẻo xuyên thế hệ (TGP) là một cơ chế thích ứng tiến hóa quan trọng. Sự thích nghi nhanh chóng này đạt được thông qua việc tái lập trình biểu sinh và chọn lọc tự nhiên các kiểu gen chịu nhiệt từ thế hệ F1. Tuy nhiên, một sự đánh đổi sinh học (trade-off) sâu sắc đã diễn ra: trong khi tỷ lệ sống được bảo tồn hoàn hảo, khả năng phục hồi của các chỉ tiêu sinh sản chỉ đạt từ 5% đến 12%, và lượng nauplii tích lũy ở thế hệ F2 vẫn thấp hơn khoảng 50% so với quần thể đối chứng. Tín hiệu cá săn mồi chỉ tạo ra biến thiên dưới 5% trong tổng thể các phản ứng, chứng minh sóng nhiệt biển là yếu tố chi phối hoàn toàn sự tồn vong của loài.

Các dữ liệu thực nghiệm đa nhân tố này khi được biểu diễn trực quan trên hệ tọa độ trực giao của phần mềm R thông qua biểu đồ hộp (box-plots) và biểu đồ tương tác đa chiều (interaction plots) đã phản ánh nhất quán xu hướng phân hóa sinh lý rõ rệt giữa các nghiệm thức.

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

Dựa trên các bằng chứng thực nghiệm thu được, luận văn đưa ra 4 nhóm giải pháp mang tính hành động cao:

Thứ nhất, xây dựng hệ thống giám sát và cảnh báo sớm sóng nhiệt biển tại các vùng sinh cảnh đầm phá ven biển miền Trung. Các trạm quan trắc hải văn cần thiết lập ngưỡng kích hoạt cảnh báo khi nhiệt độ nước biển duy trì trên 33°C liên tục quá 3 ngày, giúp các cơ sở bảo tồn và ngư dân có phương án ứng phó trước 7 đến 10 ngày. Trách nhiệm thực hiện thuộc về Viện Hải dương học phối hợp cùng Trung tâm Dự báo Khí tượng Thủy văn Quốc gia, lộ trình triển khai từ năm 2024 đến 2026.

Thứ hai, cải tiến quy trình kỹ thuật sản xuất giống thức ăn tươi sống trong nuôi trồng thủy sản thông qua giải pháp thuần dưỡng nhiệt xuyên thế hệ. Các trại sản xuất giống cần áp dụng kỹ thuật tôi luyện nhiệt độ tăng dần 1°C mỗi 12 giờ cho quần thể bố mẹ Pseudodiaptomus incisus để khai thác tính mềm dẻo TGP, hướng tới mục tiêu duy trì tỷ lệ sống của quần thể giống trên 85% trong các đợt nắng nóng mùa hè. Trách nhiệm chuyển giao công nghệ do các Viện Nghiên cứu Nuôi trồng Thủy sản đảm trách trong giai đoạn 2024-2025.

Thứ ba, nâng cao hàm lượng dinh dưỡng và chất lượng thức ăn vi tảo trong các giai đoạn sốc nhiệt. Các đơn vị nuôi trồng cần duy trì mật độ vi tảo Isochrysis galbana ổn định ở mức tối thiểu 30.000 đến 35.000 tế bào/mL, bổ sung các axit béo không no cao phân tử (PUFA) để bù đắp mức thâm hụt từ 28% đến 32% hiệu suất tiêu hóa của giáp xác khi nhiệt độ nước tăng cao. Chủ thể thực hiện là các kỹ sư thủy sản tại các trang trại nuôi trồng từ quý II năm 2024.

Thứ tư, tích hợp chỉ số nhạy cảm của các mắt xích sinh vật phù du vào quy hoạch không gian biển và đánh giá rủi ro sinh thái tổng thể. Cục Thủy sản và Bộ Tài nguyên và Môi trường cần đưa thông số biến động sinh khối giáp xác nổi vào bộ chỉ thị quan trắc sức khỏe hệ sinh thái biển ven bờ, hoàn thiện khung tiêu chuẩn kỹ thuật trong giai đoạn 2025-2027.

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

Nội dung luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cho 4 nhóm đối tượng:

Nhóm 1: Các nhà nghiên cứu và giảng viên chuyên ngành Sinh thái học biển và Thủy sinh học. Tài liệu cung cấp cơ sở dữ liệu thực nghiệm chuẩn mực về tính mềm dẻo xuyên thế hệ (TGP) dưới tác động của đa áp lực sinh thái, làm tài liệu tham khảo chất lượng cao cho các bài giảng chuyên đề sau đại học và các công trình công bố quốc tế.

Nhóm 2: Kỹ sư công nghệ nuôi trồng thủy sản và quản lý trại sản xuất giống hải sản. Luận văn cung cấp các thông số kỹ thuật chi tiết về chu kỳ phát triển 12 giai đoạn biến thái, kích thước cơ thể 680 đến 750 µm của con cái trưởng thành, cùng giải pháp tối ưu hóa sinh khối Pseudodiaptomus incisus làm thức ăn sống cao cấp cho ấu trùng cá chẽm và cá biển.

Nhóm 3: Các cơ quan quản lý nhà nước về môi trường, biển và hải đảo. Báo cáo cung cấp luận cứ khoa học định lượng giúp các nhà hoạch định chính sách xây dựng các chiến lược bảo tồn đa dạng sinh học biển, đánh giá mức độ tổn thương của nguồn lợi thủy sản ven bờ dưới tác động của biến đổi khí hậu toàn cầu.

Nhóm 4: Học viên cao học và nghiên cứu sinh khối ngành Sinh học, Môi trường và Thủy sản. Công trình là hình mẫu phương pháp luận chuẩn mực về thiết kế thí nghiệm trực giao 16 nghiệm thức phức tạp, quy trình xử lý thống kê nâng cao với các mô hình GLM, GLMER và kỹ thuật kiểm soát phần dư trong phần mềm R.

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

Sóng nhiệt biển 34°C gây tổn hại cụ thể như thế nào đối với giáp xác Pseudodiaptomus incisus?

Sóng nhiệt biển 34°C làm suy giảm 17% đến 18% tỷ lệ sống của cả con đực và con cái ở thế hệ F1. Đặc biệt, năng suất sinh sản bị tê liệt nghiêm trọng khi số lượng ấu trùng nauplii nở từ mỗi bọc trứng giảm tới 62%, đồng thời lượng thức ăn tiêu thụ đo qua viên phân giảm 28% do mất cân bằng năng lượng trao đổi chất.

Tính mềm dẻo xuyên thế hệ (TGP) có giúp loài phục hồi hoàn toàn trước sóng nhiệt biển không?

TGP giúp phục hồi hoàn toàn tỷ lệ sống của thế hệ F2 về mức tương đương với nhiệt độ tối ưu 30°C (tăng 10% so với nhóm chịu sốc nhiệt đơn lẻ). Tuy nhiên, năng suất sinh sản không thể phục hồi hoàn toàn, số lượng ấu trùng tích lũy ở F2 vẫn giảm khoảng 50% do sinh vật phải ưu tiên năng lượng cho sự sống còn.

Áp lực từ tín hiệu cá săn mồi (FPC) tương tác như thế nào với sóng nhiệt biển?

Ở nhiệt độ đối chứng 30°C, FPC kích thích phản ứng thích nghi giúp tăng khoảng 15% lượng nauplii tích lũy. Nhưng khi sóng nhiệt biển 34°C xuất hiện, tác động tích cực này hoàn toàn bị triệt tiêu vì sinh vật cạn kiệt năng lượng dự trữ, cho thấy sóng nhiệt biển là yếu tố chi phối tuyệt đối.

Hiệu ứng bố mẹ (Parental effects) kéo dài sang thế hệ sau ra sao nếu môi trường sống trở lại bình thường?

Khi thế hệ bố mẹ F1 chịu sốc nhiệt 34°C nhưng thế hệ con F2 được sống ở nhiệt độ bình thường 30°C, hiệu ứng tiêu cực từ bố mẹ vẫn làm giảm 22% kích thước bọc trứng và giảm 33% số nauplii nở ra ở F2 do nguồn dinh dưỡng noãn hoàng từ mẹ bị suy giảm nghiêm trọng.

Nghiên cứu sử dụng những công cụ phân tích thống kê nào để đảm bảo tính chuẩn xác?

Nghiên cứu ứng dụng toàn diện ngôn ngữ thống kê R phiên bản 4.x, sử dụng các mô hình phi tuyến tính bao gồm Quasibinomial GLM cho tỷ lệ sống, Binomial GLM cho tỷ lệ nở, Poisson và Negative Binomial GLM cho kích thước bọc trứng, cùng Quasipoisson GLM cho sản lượng nauplii và viên phân, tối ưu hóa qua tiêu chí thông tin AIC.

Kết luận

  • Luận văn xác định định lượng sóng nhiệt biển 34°C làm giảm 17% đến 18% tỷ lệ sống và làm suy giảm 62% số lượng ấu trùng nauplii nở ra ở thế hệ F1 của loài Pseudodiaptomus incisus.
  • Nghiên cứu chứng minh tính mềm dẻo xuyên thế hệ (TGP) giúp phục hồi hoàn toàn 100% tỷ lệ sống của thế hệ F2 nhưng năng suất sinh sản vẫn bị suy giảm khoảng 50% do sự đánh đổi năng lượng sinh học.
  • Đề tài chỉ ra hiệu ứng bố mẹ tiêu cực kéo dài làm giảm 22% kích thước bọc trứng của thế hệ con cái ngay cả khi thế hệ con được chuyển về môi trường nhiệt độ thuận lợi 30°C.
  • Tác phẩm làm sáng tỏ sự tương tác sinh thái đa yếu tố, khẳng định sóng nhiệt biển có vai trò áp đảo hoàn toàn so với áp lực hóa học từ động vật ăn thịt (chiếm hơn 95% mức độ ảnh hưởng).
  • Kế hoạch nghiên cứu tiếp theo giai đoạn 2024-2027 cần tập trung giải mã cơ chế methyl hóa ADN tế bào và kiểm chứng mô hình sinh thái thực địa tại các đầm phá Vịnh Cam Ranh.

Hãy liên hệ trực tiếp với tác giả hoặc Khoa Sinh học, Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội để tiếp cận toàn văn công trình nghiên cứu và khai thác bộ dữ liệu thực nghiệm có giá trị này phục vụ công tác nghiên cứu và sản xuất thủy sản bền vững.