Tổng quan luận án

Luận án tiến sĩ với đề tài "Diễn thế nguyên sinh rừng ngập mặn ở Cồn Ông Trang thuộc Khu Dự trữ sinh quyển Mũi Cà Mau" được thực hiện bởi nghiên cứu sinh Lư Ngọc Trâm Anh, dưới sự hướng dẫn khoa học của PGS. Viên Ngọc Nam và PGS. Nguyễn Thị Phương Thảo, chuyên ngành Sinh thái học (mã số 9.20), bảo vệ tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam vào năm 2020.

Tính cấp thiết và khoảng trống nghiên cứu

Hệ sinh thái rừng ngập mặn giữ vai trò sinh thái, môi trường và kinh tế quan trọng ở vùng cửa sông ven biển nhiệt đới, đặc biệt trong việc bồi tụ phù sa, chống xói lở bờ biển, hạn chế xâm nhập mặn, bảo tồn đa dạng sinh học và hấp thụ carbon nhằm giảm thiểu tác động của biến đổi khí hậu. Khu vực Mũi Cà Mau là trung tâm phân bố rừng ngập mặn ngập triều tự nhiên lớn nhất Việt Nam. Trong đó, hệ thống Cồn Ông Trang (gồm Cồn Trong, Cồn Ngoài và Cồn Mới) nằm ở cửa sông Cửa Lớn, là bãi bồi cửa sông hình thành hoàn toàn tự nhiên và chịu sự bảo vệ nghiêm ngặt của Vườn Quốc gia Mũi Cà Mau.

Tuy nhiên, phần lớn các công trình nghiên cứu trước đây tại Việt Nam về thảm thực vật ngập mặn chủ yếu tập trung mô tả định tính thành phần loài, phân loại hình thái hoặc khảo sát lát cắt không gian đơn lẻ tại một thời điểm. Vẫn còn thiếu các nghiên cứu định lượng toàn diện kết hợp giải đoán chuỗi thời gian hình thành đất bồi với khảo sát đa dạng loài, cấu trúc lâm phần, đặc tính thổ nhưỡng, chế độ thủy triều và trữ lượng carbon tích lũy theo từng giai đoạn diễn thế nguyên sinh. Việc nghiên cứu diễn thế nguyên sinh tại một địa bàn không chịu tác động trực tiếp của con người như Cồn Ông Trang là cơ sở khoa học cần thiết để làm sáng tỏ quy luật diễn thế tự nhiên và hỗ trợ công tác quy hoạch, phục hồi rừng.

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

Luận án xác định 3 mục tiêu nghiên cứu cụ thể được đánh số như sau:

  1. Xác định được ranh giới và diện tích các khu vực của Cồn Ông Trang có thời gian hình thành khác nhau.
  2. Phân tích được quá trình diễn thế nguyên sinh của rừng ngập mặn ở Cồn Ông Trang cũng như xây dựng được sơ đồ diễn thế nguyên sinh của khu vực.
  3. Xác định được sự phân bố của các loài thực vật ngập mặn, cấu trúc, đa dạng sinh học và trữ lượng carbon của rừng ngập mặn ở Cồn Ông Trang, đồng thời phân tích được ảnh hưởng của các nhân tố sinh thái cơ bản đến cấu trúc, đa dạng và phân bố thực vật ngập mặn theo không gian và thời gian.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Hệ sinh thái rừng ngập mặn tại 3 cồn ở sông Cửa Lớn gồm Cồn Trong, Cồn Ngoài và Cồn Mới (gọi chung là Cồn Ông Trang) thuộc Vườn Quốc gia Mũi Cà Mau / Khu Dự trữ sinh quyển Mũi Cà Mau.
  • Phạm vi nghiên cứu:
    • Không gian: Khu vực 3 cồn bồi tụ tự nhiên tại Cửa Ông Trang.
    • Thời gian hình thành (các giai đoạn diễn thế): Được phân chia thành 5 giai đoạn theo chuỗi thời gian số hóa bản đồ và ảnh viễn thám, gồm Giai đoạn 1 (2004–2016), Giai đoạn 2 (1992–2004), Giai đoạn 3 (1979–1992), Giai đoạn 4 (1962–1979), và Giai đoạn 5 (trước năm 1962).
    • Các chỉ tiêu môi trường và sinh thái: Thành phần loài thực vật ngập mặn; đường kính ngang ngực ($Dbh$), chiều cao vút ngọn ($H$), mật độ ($N$), tiết diện ngang ($BA$), thể tích thân cây ($V$), chỉ số giá trị quan trọng ($IVI$); độ cao địa hình, chế độ và độ ngập triều; đặc tính thổ nhưỡng gồm thể nền, thành phần cơ giới đất, độ mặn ($Sal$), $pH$, hàm lượng $NPK$ tổng số ở hai tầng đất ($0-20\text{ cm}$ và $20-60\text{ cm}$); sinh khối tươi, sinh khối khô và lượng carbon tích tụ trong sinh khối trên mặt đất ($CAGB$), sinh khối dưới mặt đất ($CBGB$) và trong đất.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước đây được điểm lại trong văn bản

Tác giả đã tổng thuật các hướng nghiên cứu trong nước và quốc tế liên quan đến bốn nhóm vấn đề chính:

Nhóm vấn đề Các công trình và tác giả quốc tế Các công trình và tác giả tại Việt Nam
Thành phần loài và đa dạng sinh học rừng ngập mặn Tomlinson (1986); Spalding và cộng sự (1997); Mandal và Naskar (2008); Hossain và cộng sự (2015); Alappatt (2018); hệ thống phân loại phát sinh loài thực vật có hoa APG IV (2016). Phan Nguyên Hồng và Hoàng Thị Sản (1993); Phan Nguyên Hồng và cộng sự (1999); Huỳnh Quốc Tính và cộng sự (2009); Đặng Văn Sơn (2014); Phạm Văn Ngọt và cộng sự (2013); Vien Ngoc Nam và cộng sự (2014); Phạm Hồng Tính và Mai Sỹ Tuấn (2016).
Cấu trúc rừng và ảnh hưởng của nhân tố sinh thái Chen và Twilley (1998); Calegario và cộng sự (2015); Seedo và cộng sự (2017); Datta và Deb (2017); Duke (2012); Clough (2014); Van Loon và cộng sự (2016). Phùng Trung Ngân và Châu Quang Hiển (1987); Thái Văn Trừng (1998); Nguyễn Hoàng Trí (1999); Đỗ Đình Sâm và cộng sự (2005); Đặng Trung Tấn (2007); Hoàng Văn Thơi và cộng sự (2010, 2012); Lê Tấn Lợi (2011); Nguyễn Văn Tú và Bùi Lai (2012); Võ Ngươn Thảo và cộng sự (2013).
Sinh khối và tích tụ carbon Murdiyarso và cộng sự (2009); Donato và cộng sự (2011); Kauffman và Donato (2012). Viên Ngọc Nam và cộng sự (2009, 2010, 2011, 2016); Cao Huy Bình và Viên Ngọc Nam (2014); N. Tue và cộng sự (2014); Nguyễn Thị Hồng Hạnh và cộng sự (2015, 2017); Lê Tấn Lợi và Lý Hằng Ni (2015); Nguyễn Hà Quốc Tín và Lê Tấn Lợi (2015).
Lý thuyết và mô hình diễn thế sinh thái Clements (1916); Walker và Moral (2003); Ball (1980); Cunha-Lignon và cộng sự (2009). Phan Nguyên Hồng và Hoàng Thị Sản (1993); Đỗ Đình Sâm và cộng sự (2005); Brake và Huijgevoort (2008); Nguyễn Hà Quốc Tín và cộng sự (2014).

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

Mặc dù Cồn Ngoài hoặc Cồn Trong đã từng được khảo sát sơ bộ trong một số báo cáo ngắn hạn (như Brake và Huijgevoort, 2008; Nguyễn Hà Quốc Tín và cộng sự, 2014, 2015), nhưng các nghiên cứu này chỉ khảo sát đơn lẻ, thiếu tính hệ thống trên toàn bộ 3 cồn và chưa làm rõ quá trình diễn thế nguyên sinh gắn với trục thời gian hình thành địa mạo bãi bồi. Luận án giải quyết khoảng trống này bằng cách kết hợp phương pháp không gian thay thế thời gian (space-for-time substitution) dựa trên ảnh viễn thám và bản đồ lịch sử, phân tích định lượng đồng thời cả cấu trúc thực vật, quan hệ đa biến của các yếu tố sinh thái thổ nhưỡng - thủy triều và trữ lượng carbon trong các bể chứa theo 5 giai đoạn hình thành.

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

Cơ sở lý thuyết và khung phân loại

  • Lý thuyết diễn thế nguyên sinh: Dựa trên quan điểm diễn thế sinh thái cổ điển của Clements (1916) và Walker & Moral (2003), xác định diễn thế nguyên sinh là quá trình phát triển của hệ sinh thái trên nền bãi bồi mới tạo lập chưa từng có thực vật sinh sống, chuyển dịch từ quần xã tiên phong đến quần xã ổn định.
  • Khung phân loại thực vật: Luận án áp dụng hệ thống phân loại APG IV (2016), theo đó chi Mấm (Avicennia) được chuyển từ họ Mấm (Avicenniaceae) / họ Cỏ roi ngựa (Verbenaceae) sang họ Ô rô (Acanthaceae); chi Bần (Sonneratia) được xếp thành một phân họ trong họ Bằng lăng (Lythraceae).
  • Phân vùng chế độ ngập triều: Áp dụng phân loại 4 vùng ngập triều theo Đỗ Đình Sâm và cộng sự (2005): vùng ngập triều thường xuyên, vùng ngập triều thấp, vùng ngập triều trung bình, và vùng ngập bởi triều cao.

Phương pháp nghiên cứu thực tế

Luận án phối hợp phương pháp tiếp cận không gian thay thế thời gian, điều tra thực địa theo ô tiêu chuẩn và phân tích thống kê định lượng:

  1. Xác định các giai đoạn hình thành cồn: Sử dụng các tư liệu bản đồ và ảnh vệ tinh lịch sử qua các thời kỳ (trước 1962, 1962, 1979, 1992, 2004, 2016) để số hóa, xác định ranh giới và tính toán diện tích tăng thêm của Cồn Trong, Cồn Ngoài và Cồn Mới qua 5 giai đoạn diễn thế ($gd1$ đến $gd5$).
  2. Điều tra lâm học bằng ô tiêu chuẩn (OTC): Thiết lập các tuyến điều tra cắt ngang các cồn và bố trí hệ thống ô tiêu chuẩn đại diện cho các giai đoạn hình thành đã số hóa. Đo đếm tất cả các cây gỗ có $Dbh \ge 5\text{ cm}$ về đường kính ngang ngực ($Dbh$), chiều cao vút ngọn ($H$), xác định số lượng cá thể để tính mật độ ($N$), tiết diện ngang ($BA$) và thể tích thân cây ($V$).
  3. Tính toán các chỉ số đa dạng và cấu trúc:
    • Chỉ số giá trị quan trọng: $IVI (%) = \frac{RN + RF + RBA}{3}$ (trong đó $RN$ là mật độ tương đối, $RF$ là tần suất xuất hiện tương đối, $RBA$ là tiết diện ngang tương đối).
    • Các chỉ số đa dạng sinh thái: Chỉ số đa dạng Shannon-Wiener ($H'$), chỉ số đồng đều Pielou ($J'$), chỉ số phong phú Margalef ($d'$), chỉ số ưu thế Simpson.
    • Xây dựng phương trình tương quan $H - Dbh$ và phân bố số cây theo cấp đường kính ($N - Dbh$), cấp chiều cao ($N - H$).
  4. Khảo sát môi trường và phân tích thổ nhưỡng:
    • Đo đạc cao trình mặt đất và ghi nhận chế độ ngập triều thực tế đối chiếu số liệu trạm hải văn Năm Căn.
    • Thu thập mẫu đất tại 2 tầng độ sâu: tầng mặt $0 - 20\text{ cm}$ ($T_{0-20}$) và tầng sâu $20 - 60\text{ cm}$ ($T_{20-60}$).
    • Phân tích thành phần cơ giới đất (cát, limon, sét theo tam giác phân loại đất của Bộ Nông nghiệp Hoa Kỳ - USDA), độ $pH$ ($pH_{20}, pH_{60}$), độ mặn đất ($Sal_{20}, Sal_{60}$), hàm lượng dinh dưỡng tổng số ($N_{20}, N_{60}, P_{20}, P_{60}, K_{20}, K_{60}$) và hàm lượng carbon hữu cơ trong đất.
  5. Ước tính sinh khối và trữ lượng carbon:
    • Sử dụng các phương trình tương quan sinh khối (allometric equations) chuyên biệt cho từng loài cây ngập mặn kết hợp tỉ trọng gỗ để tính sinh khối tươi, sinh khối khô trên mặt đất ($AGB$) và sinh khối rễ dưới mặt đất ($BGB$).
    • Tính toán trữ lượng carbon tích tụ trong sinh khối ($CAGB, CBGB$) bằng cách nhân sinh khối khô với hệ số chuyển đổi carbon tương ứng của từng loài.
    • Tính lượng carbon tích tụ trong đất cho 2 tầng đất $0-20\text{ cm}$ và $20-60\text{ cm}$ từ hàm lượng carbon phân tích và dung trọng đất.
  6. Xử lý số liệu và phân tích đa biến:
    • Sử dụng thống kê mô tả (trung bình $\bar{X}$, sai số chuẩn $SE$, độ lệch chuẩn $SD$) và kiểm định ANOVA / so sánh trung bình.
    • Phân tích cụm nhóm loài (Cluster analysis) dựa trên khoảng cách tương đồng Bray-Curtis.
    • Phân tích sắp xếp đa chiều không gian (MDS) để đánh giá sự phân bố loài theo chế độ ngập triều và thể nền.
    • Phân tích thành phần chính (PCA) để xác định mối quan hệ giữa phân bố quần xã thực vật với các biến môi trường thổ nhưỡng.

Nội dung chính theo từng chương

Chương 1. Tổng quan vấn đề nghiên cứu

Chương 1 hệ thống hóa các cơ sở lý luận về hệ sinh thái rừng ngập mặn trên thế giới và tại Việt Nam. Nội dung tập trung phân tích:

  • Đặc điểm thảm thực vật rừng ngập mặn, các nhóm cây ngập mặn chính thức và tham gia; cập nhật vị trí phân loại học của chi Avicennia và Sonneratia theo APG IV (2016).
  • Quy luật phân bố địa lý của rừng ngập mặn toàn cầu (khoảng $30^\circ\text{B} - 30^\circ\text{N}$, tập trung lớn nhất tại vùng Ấn Độ - Tây Thái Bình Dương) và phân vùng 4 khu vực rừng ngập mặn tại Việt Nam (Đông Bắc, Đồng bằng Bắc Bộ, Trung Bộ và Nam Bộ).
  • Các nhân tố sinh thái kiểm soát cấu trúc rừng: chế độ ngập triều, biên độ triều, cao trình địa hình, tính chất cơ lý và hóa học của đất ($pH$, độ mặn, $NPK$).
  • Vai trò tích tụ carbon trong 4 bể chứa (sinh khối trên mặt đất, sinh khối dưới mặt đất, đất và thảm mục/vật rụng), trong đó bể carbon dưới mặt đất và đất ngập mặn chiếm tỉ trọng vượt trội.
  • Lý thuyết về diễn thế sinh thái, các pha diễn thế nguyên sinh từ bãi bồi cửa sông ra biển, và tổng kết các công trình nghiên cứu diễn thế, sinh khối đã thực hiện tại vùng rừng ngập mặn Mũi Cà Mau và Cồn Ông Trang.

Chương 2. Địa điểm và phương pháp nghiên cứu

Chương 2 mô tả chi tiết địa bàn thực nghiệm và quy trình kỹ thuật:

  • Địa điểm nghiên cứu: Cồn Trong, Cồn Ngoài và Cồn Mới thuộc cửa sông Cửa Lớn, xã Viên An, huyện Ngọc Hiển, tỉnh Cà Mau. Địa bàn mang đặc trưng khí hậu nhiệt đới gió mùa cận xích đạo với lượng mưa lớn và chế độ nhật triều không đều của biển Tây kết hợp bán nhật triều biển Đông (theo số liệu trạm Năm Căn).
  • Phương pháp luận phân chia 5 giai đoạn hình thành cồn ($gd1: 2004-2016$, $gd2: 1992-2004$, $gd3: 1979-1992$, $gd4: 1962-1979$, $gd5: \text{trước } 1962$) dựa trên giải đoán chuỗi ảnh viễn thám và bản đồ địa hình.
  • Thiết kế mạng lưới tuyến điều tra và ô tiêu chuẩn ($OTC$) trên cả 3 cồn; quy trình đo đếm lâm học đường kính $Dbh$, chiều cao $H$, lập phương trình tương quan.
  • Quy trình lấy mẫu đất tầng $0-20\text{ cm}$ và $20-60\text{ cm}$, các phương pháp phân tích phòng thí nghiệm xác định thành phần cơ giới, $pH$, độ mặn, $N$, $P$, $K$, carbon hữu cơ.
  • Các công thức toán - sinh thái tính chỉ số $IVI$, độ đa dạng ($H', J', d'$), các phương trình sinh khối $AGB, BGB$, carbon sinh khối, carbon đất và phương pháp phân tích thống kê đa biến (Cluster, MDS, PCA).

Chương 3. Kết quả và thảo luận

Chương 3 trình bày các kết quả thực nghiệm và thảo luận chuyên sâu qua 5 nội dung chính:

3.1. Sự hình thành các cồn ở Cửa Ông Trang theo thời gian

  • Phân tích chuỗi ảnh vệ tinh và bản đồ lịch sử ghi nhận quá trình bồi tụ mở rộng liên tục diện tích đất bồi từ sông Cửa Lớn đổ ra biển.
  • Cồn Trong là cồn hình thành sớm nhất, có đầy đủ 5 giai đoạn diễn thế ($gd1$ đến $gd5$). Tiếp theo là Cồn Ngoài hình thành ở phía ngoài Cồn Trong, và Cồn Mới là bãi bồi mới nhất đang trong quá trình tích tụ phù sa ban đầu.
  • Xác định diện tích tăng thêm và ranh giới cụ thể của từng cồn tương ứng với các mốc thời gian trước năm 1962, 1962–1979, 1979–1992, 1992–2004 và 2004–2016.

3.2. Diễn thế nguyên sinh ở Cồn Trong

  • Môi trường: Đất có xu hướng chuyển dần từ bùn lỏng sét pha limon ở khu vực mới bồi sang đất sét chặt ở các giai đoạn hình thành lâu đời; độ cao địa hình tăng dần từ mép nước bồi tụ vào trung tâm cồn kéo theo số ngày và thời gian ngập triều giảm dần. Độ mặn và hàm lượng dinh dưỡng có sự biến đổi theo độ sâu tầng đất ($0-20\text{ cm}$ so với $20-60\text{ cm}$) và giai đoạn hình thành.
  • Thành phần loài và cấu trúc: Ghi nhận sự thay đổi rõ rệt của chỉ số giá trị quan trọng ($IVI$). Ở giai đoạn mới hình thành ($gd1$), loài Mấm trắng (Avicennia alba) đóng vai trò tiên phong chiếm ưu thế tuyệt đối. Sang các giai đoạn già hơn ($gd2 \to gd5$), Đước đôi (Rhizophora apiculata) tăng nhanh về mật độ và tiết diện ngang, dần thay thế Mấm trắng để trở thành loài ưu thế sinh thái đỉnh cực.
  • Phân tích đa biến: Sơ đồ Cluster và MDS chỉ rõ sự phân nhóm thực vật theo chế độ ngập triều và thể nền đất. Phân tích PCA xác định các biến thổ nhưỡng ($pH$, $Sal$, hàm lượng sét, $N, P, K$) có tương quan chặt chẽ với sự sắp xếp không gian của quần xã thực vật.
  • Tích tụ carbon: Trữ lượng carbon sinh khối ($CAGB, CBGB$) và carbon tích lũy trong đất tầng $0-60\text{ cm}$ tăng dần theo cấp tuổi của giai đoạn diễn thế. Xây dựng sơ đồ mô phỏng diễn thế nguyên sinh hoàn chỉnh tại Cồn Trong.

3.3. Diễn thế nguyên sinh ở Cồn Ngoài

  • Môi trường: Chịu tác động trực tiếp hơn của sóng và thủy triều so với Cồn Trong. Thể nền bãi bồi ven rìa ngoài có độ ngập triều cao, hàm lượng bùn sét biến động.
  • Thành phần loài và cấu trúc: Quy luật thay thế loài diễn ra tương tự Cồn Trong nhưng có sự tham gia đáng kể của các loài đi kèm như Vẹt tách (Bruguiera parviflora), Mấm đen (Avicennia officinalis), Bần trắng (Sonneratia alba). Mấm trắng chiếm ưu thế ở vành đai ngoài (bãi bồi thấp), sau đó chuyển sang dải rừng hỗn giao Mấm trắng - Đước đôi, và đi sâu vào vùng đất cao bên trong là quần xã Đước đôi thuần loại hoặc hỗn giao với Vẹt tách.
  • Cấu trúc và carbon: Các chỉ tiêu $Dbh$, $H$, tiết diện ngang $BA$ và trữ lượng rừng $M$ ($m^3/\text{ha}$) tăng dần từ $gd1$ đến $gd5$. Trữ lượng carbon trong sinh khối cây đứng và trong đất tăng tương ứng theo thời gian hình thành cồn. Xây dựng sơ đồ diễn thế nguyên sinh tại Cồn Ngoài.

3.4. Diễn thế nguyên sinh ở Cồn Mới

  • Đặc điểm sinh cảnh: Cồn Mới là cồn mới xuất hiện gần đây ($gd1$), đại diện cho pha khởi đầu của quá trình diễn thế nguyên sinh trên bãi bồi ngập triều thường xuyên với thể nền bùn lỏng, địa hình thấp.
  • Thực vật và carbon: Đa dạng loài ở mức thấp, thảm thực vật chủ yếu gồm các cây con và lâm phần non của loài tiên phong Mấm trắng (Avicennia alba). Chỉ số $IVI$ của Mấm trắng chiếm vị trí chủ đạo; kích thước $Dbh$ và $H$ nhỏ, mật độ cao nhưng trữ lượng và sinh khối tích tụ còn thấp. Bể chứa carbon đất chiếm tỷ trọng áp đảo so với bể carbon sinh khối trên mặt đất.

3.5. Diễn thế nguyên sinh tổng thể ở Cồn Ông Trang

  • Tổng hợp quy luật phân bố và đa dạng: Tổng kết toàn diện sự biến động cấu trúc thực vật của 3 cồn theo không gian và thời gian. Phân bố thực vật ngập mặn Cồn Ông Trang tuân theo quy luật thay thế liên tục: từ quần thể tiên phong đơn loài Mấm trắng trên đất bùn ngập triều sâu $\to$ quần xã hỗn giao Mấm trắng + Đước đôi $\to$ quần xã Đước đôi chiếm ưu thế cùng Vẹt tách trên nền đất cao, chặt, ít ngập triều.
  • Quy luật cấu trúc: Phân bố số cây theo cấp đường kính ($N - Dbh$) và cấp chiều cao ($N - H$) của lâm phần ở các giai đoạn có dạng phân bố giảm (lệch trái), phản ánh đặc trưng của quần xã rừng tự nhiên đang trong quá trình tái sinh và phát triển liên tục.
  • Tổng hợp trữ lượng carbon: Định lượng toàn diện các bể carbon của Cồn Ông Trang:
    • Sinh khối trên mặt đất ($AGB$) và dưới mặt đất ($BGB$) tăng đều qua các giai đoạn $gd1 \to gd5$.
    • Lượng carbon trong đất tầng $0-60\text{ cm}$ chiếm tỷ trọng lớn trong tổng tích lũy carbon của toàn hệ sinh thái, trong đó tầng $0-20\text{ cm}$ có hàm lượng carbon biến động phụ thuộc vào lượng vật rơi rụng và thảm mục hữu cơ bề mặt.
    • Tổng trữ lượng carbon trung bình của hệ sinh thái rừng ngập mặn Cồn Ông Trang tăng dần theo thời gian hình thành cồn.
  • Sơ đồ diễn thế tổng quát: Luận án hoàn thiện sơ đồ mô hình hóa diễn thế nguyên sinh rừng ngập mặn Cồn Ông Trang phản ánh sự tương tác đồng thời giữa động thái bồi tụ địa mạo, biến đổi lý hóa tính thổ nhưỡng, chế độ ngập triều và sự biến đổi tổ thành loài thực vật.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận

  1. Luận án đã xác định và lượng hóa được mối liên hệ chặt chẽ giữa quá trình diễn thế nguyên sinh của rừng ngập mặn với sự biến đổi các nhân tố sinh thái (chế độ ngập triều, cao trình địa hình, thành phần cơ giới và lý hóa tính của đất) theo chuỗi thời gian hình thành địa mạo bãi bồi tự nhiên.
  2. Xây dựng thành công sơ đồ mô phỏng quá trình diễn thế nguyên sinh tổng thể và chi tiết cho từng cồn tại Cồn Ông Trang, chứng minh quy luật thay thế từ quần xã tiên phong thích nghi bùn lỏng ngập triều (Avicennia alba) tiến đến quần xã đỉnh cực thích nghi nền đất bồi cao, ổn định (Rhizophora apiculata).
  3. Đóng góp dữ liệu định lượng về động thái biến đổi của các chỉ số cấu trúc lâm phần ($N, Dbh, H, BA, V, IVI$), đa dạng sinh học ($H', J', d'$) và các bể chứa carbon ($AGB, BGB, CAGB, CBGB$, carbon đất $0-60\text{ cm}$) theo các giai đoạn diễn thế nguyên sinh tại vùng cửa sông ven biển phía Nam Việt Nam.
  4. Cập nhật và áp dụng hệ thống phân loại hiện đại APG IV (2016) vào việc định danh và đánh giá phát sinh loài thực vật ngập mặn tại hệ sinh thái Mũi Cà Mau.

Đóng góp mới về mặt thực tiễn

  1. Cung cấp cơ sở khoa học định lượng phục vụ công tác quy hoạch, bảo tồn và phát triển rừng ngập mặn tại Vườn Quốc gia Mũi Cà Mau; xác định đúng các loài ưu thế sinh thái tương ứng với từng điều kiện lập địa và giai đoạn bồi tụ.
  2. Cung cấp luận cứ khoa học cho việc lựa chọn cơ cấu loài cây trồng rừng ngập mặn phù hợp: đề xuất các phương thức trồng rừng thuần loài hoặc hỗn giao dựa trên kết quả phân tích mức độ tương đồng và liên kết sinh thái giữa các loài (như các tổ hợp Đước đôi, Mấm trắng, Vẹt tách).
  3. Cung cấp bộ số liệu thực tế về sinh khối và khả năng hấp thụ, tích lũy carbon của rừng ngập mặn Cồn Ông Trang, tạo tiền đề khoa học cho việc xây dựng chính sách chi trả dịch vụ môi trường rừng (PES) và tham gia các cơ chế giảm phát thải khí nhà kính (REDD, REDD+).

Kiến nghị và giải pháp đề xuất

  • Duy trì chế độ bảo vệ nghiêm ngặt các vùng bãi bồi tự nhiên tại Cồn Ông Trang để quá trình diễn thế nguyên sinh diễn ra tự nhiên không bị xáo trộn.
  • Áp dụng các sơ đồ diễn thế và kết quả phân nhóm loài thích nghi theo thể nền/chế độ triều vào các dự án trồng rừng và phục hồi bãi bồi ven biển tại đồng bằng sông Cửu Long.
  • Tiếp tục duy trì các ô tiêu chuẩn định vị để theo dõi diễn biến lâu dài của các bể carbon và cấu trúc rừng dưới tác động của biến đổi khí hậu.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế

  • Luận án sử dụng phương pháp tiếp cận không gian thay thế thời gian (space-for-time substitution) dựa trên các giai đoạn số hóa ảnh viễn thám; phương pháp này có thể chịu ảnh hưởng cục bộ bởi các yếu tố biến động ngẫu nhiên về khí hậu hoặc dòng chảy thủy văn qua các năm.
  • Khảo sát carbon hữu cơ trong đất chỉ thực hiện đến độ sâu $60\text{ cm}$ (chia 2 tầng $0-20\text{ cm}$ và $20-60\text{ cm}$), chưa đánh giá được toàn bộ trữ lượng carbon ở các tầng trầm tích bùn cát sâu hơn ($>100\text{ cm}$).

Hướng nghiên cứu tiếp

  • Tiếp tục theo dõi đo đếm định kỳ trên mạng lưới ô tiêu chuẩn cố định tại 3 cồn nhằm kiểm chứng tốc độ diễn thế thực tế theo thời gian thực.
  • Mở rộng nghiên cứu khoan lấy mẫu trầm tích sâu ($1-3\text{ m}$) để đánh giá toàn diện bể chứa carbon đất và lịch sử bồi lắng trầm tích cửa sông.
  • Nghiên cứu sâu hơn về động thái dòng chảy, sự tương tác giữa sóng biển và rễ cây ngập mặn trong việc giữ phù sa bồi tụ bãi bồi cửa sông Cửa Lớn dưới bối cảnh nước biển dâng.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo có giá trị học thuật và thực tiễn cao cho nhiều nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và sinh viên các ngành Sinh thái học, Lâm nghiệp, Quản lý tài nguyên rừng, Quản lý môi trường và Khoa học biển: Tham khảo phương pháp luận tích hợp viễn thám GIS với điều tra ô tiêu chuẩn lâm học, phương pháp phân tích thống kê sinh thái đa biến (Cluster, MDS, PCA) và phương pháp tính toán sinh khối, carbon rừng.
  • Giảng viên và nhà nghiên cứu sinh học/lâm nghiệp: Sử dụng dữ liệu thực nghiệm về cấu trúc, tổ thành loài, hệ số tương quan $H-Dbh$ và sơ đồ diễn thế nguyên sinh phục vụ giảng dạy các chuyên đề sinh thái học thực vật, sinh thái học rừng ngập mặn.
  • Nhà quản lý, cán bộ kiểm lâm và cơ quan hoạch định chính sách: Ban Quản lý Vườn Quốc gia Mũi Cà Mau, Hạt Kiểm lâm, Sở Nông nghiệp và Phát triển nông thôn các tỉnh ven biển đồng bằng sông Cửu Long sử dụng các dẫn liệu về loài ưu thế, điều kiện lập địa và trữ lượng carbon để xây dựng phương án trồng rừng ven biển, khoanh nuôi xúc tiến tái sinh tự nhiên và triển khai chi trả dịch vụ môi trường rừng.

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

1. Quá trình diễn thế nguyên sinh tại Cồn Ông Trang được chia thành mấy giai đoạn và căn cứ trên cơ sở nào?
Quá trình diễn thế nguyên sinh tại Cồn Ông Trang được phân chia thành 5 giai đoạn theo chuỗi thời gian hình thành diện tích cồn bồi: Giai đoạn 1 ($gd1$: 2004–2016), Giai đoạn 2 ($gd2$: 1992–2004), Giai đoạn 3 ($gd3$: 1979–1992), Giai đoạn 4 ($gd4$: 1962–1979) và Giai đoạn 5 ($gd5$: trước năm 1962). Sự phân chia này căn cứ vào việc giải đoán và số hóa ranh giới bồi tụ cồn từ chuỗi tư liệu bản đồ và ảnh vệ tinh lịch sử.

2. Luận án áp dụng hệ thống phân loại nào cho các loài thực vật ngập mặn chủ yếu?
Luận án áp dụng hệ thống phân loại thực vật APG IV (2016) dựa trên phát sinh loài ở mức phân tử. Theo hệ thống này, chi Mấm (Avicennia) được xếp vào họ Ô rô (Acanthaceae) thay vì họ Mấm (Avicenniaceae) hay họ Cỏ roi ngựa (Verbenaceae); chi Bần (Sonneratia) được xếp vào họ Bằng lăng (Lythraceae) thay vì họ Bần (Sonneratiaceae).

3. Quy luật thay thế loài ưu thế diễn ra như thế nào trong quá trình diễn thế nguyên sinh tại Cồn Ông Trang?
Tại các bãi bồi mới hình thành ($gd1$, Cồn Mới, mép ngoài Cồn Trong và Cồn Ngoài) có địa hình thấp, bùn lỏng và ngập triều thường xuyên, loài Mấm trắng (Avicennia alba) đóng vai trò tiên phong với chỉ số giá trị quan trọng ($IVI$) cao nhất. Khi đất được bồi tụ cao hơn, nén chặt hơn và thời gian ngập triều giảm ($gd2 \to gd5$), Đước đôi (Rhizophora apiculata) xâm lấn, cạnh tranh và thay thế Mấm trắng để trở thành loài chiếm ưu thế tuyệt đối trong quần xã rừng đỉnh cực ổn định, đi kèm với một số loài như Vẹt tách (Bruguiera parviflora).

4. Trữ lượng carbon của rừng ngập mặn Cồn Ông Trang được đánh giá thông qua các bể chứa nào?
Trữ lượng carbon được đánh giá định lượng qua 3 bể chứa chính:

  • Bể carbon trong sinh khối gỗ trên mặt đất ($CAGB$);
  • Bể carbon trong sinh khối gỗ rễ dưới mặt đất ($CBGB$);
  • Bể carbon tích lũy trong đất ở hai tầng độ sâu $0-20\text{ cm}$ và $20-60\text{ cm}$.

5. Các phương pháp thống kê đa biến nào được sử dụng để phân tích quan hệ giữa thảm thực vật và môi trường?
Luận án sử dụng ba phương pháp phân tích đa biến chính: Phân tích cụm (Cluster analysis) dựa trên hệ số tương đồng Bray-Curtis để phân nhóm các loài thực vật; Phân tích sắp xếp đa chiều không gian (MDS) để đánh giá sự phân bố của loài theo chế độ ngập triều và thể nền đất; và Phân tích thành phần chính (PCA) để xác định mức độ tác động của các yếu tố thổ nhưỡng ($pH$, độ mặn, $NPK$, thành phần cơ giới) đến cấu trúc quần xã thực vật.

Kết luận

Luận án của nghiên cứu sinh Lư Ngọc Trâm Anh là công trình nghiên cứu sinh thái học thực nghiệm toàn diện về diễn thế nguyên sinh rừng ngập mặn tại hệ thống Cồn Ông Trang thuộc Khu Dự trữ sinh quyển Mũi Cà Mau. Bằng việc kết hợp phương pháp chuỗi thời gian số hóa viễn thám với điều tra ô tiêu chuẩn lâm học và phân tích thổ nhưỡng đa biến, tác giả đã làm sáng tỏ quy luật thay thế thành phần loài, biến động cấu trúc lâm phần và tích lũy carbon trong các bể chứa qua 5 giai đoạn hình thành cồn bồi. Các sơ đồ diễn thế nguyên sinh và bộ dẫn liệu định lượng được xây dựng trong luận án cung cấp cơ sở khoa học trực tiếp cho công tác bảo tồn, phục hồi rừng ngập mặn và hoạch định chính sách chi trả dịch vụ môi trường rừng tại vùng bờ biển Nam Bộ.