Nghiên Cứu, Phân Tích Chỉ Số Đa Dạng Sinh Học Thực Vật Thân Gỗ Tại Khu Bảo Tồn Thiên Nhiên Sơn Trà, TP. Đà Nẵng
Tóm tắt nghiên cứu
Nghiên cứu tập trung giải quyết câu hỏi then chốt: Cấu trúc quần xã và mức độ đa dạng sinh học của hệ thực vật thân gỗ tại Khu bảo tồn thiên nhiên (BTTN) Sơn Trà hiện trạng ra sao dưới tác động cộng hưởng của các yếu tố sinh thái và áp lực nhân sinh? Để trả lời vấn đề này, đề tài đã ứng dụng phương pháp tiếp cận định lượng sinh thái học quần xã thông qua điều tra hệ thống 12 ô tiêu chuẩn ($500\text{ m}^2/\text{ô}$) đại diện cho 5 dạng sinh cảnh: rừng tự nhiên, rừng tự nhiên Chò đen, rừng trồng, trảng cỏ và đất trống.
Kết quả đã ghi nhận 96 loài thực vật thân gỗ thuộc 46 họ. Phân tích định lượng chỉ ra trật tự ưu thế rõ nét với loài Chò đen (Parashorea stellata Kurz) giữ vị trí thống lĩnh ($\text{IVI} = 35{,}38$). Chỉ số Shannon-Wiener ($H'$) dao động rộng từ $1{,}62$ đến $4{,}76$, phản ánh sự phân hóa mạnh mẽ giữa các sinh cảnh nguyên sinh và thoái hóa. Đồng thời, chỉ số chiếm ưu thế Simpson ($C_d$) biến thiên từ $0{,}051$ đến $0{,}499$, cảnh báo nguy cơ suy giảm tính bền vững của các quần xã chịu xáo trộn. Hệ số tương đồng Sorensen ($SI = 0{,}04 - 0{,}50$) chứng minh sự cách biệt sâu sắc về thành phần loài qua các trạng thái sinh cảnh. Nghiên cứu cung cấp cơ sở dữ liệu định lượng chuẩn xác, làm tiền đề khoa học cho công tác quy hoạch phân khu chức năng, kiểm soát sinh vật xâm hại và thực thi chính sách bảo tồn tài nguyên gen thực vật đặc hữu tại khu vực Duyên hải Nam Trung Bộ.
Bối cảnh và tầm quan trọng
Việt Nam được cộng đồng quốc tế công nhận là một trong những trung tâm đa dạng sinh học (ĐDSH) hàng đầu thế giới với hệ sinh thái nhiệt đới phong phú, là sinh cảnh của khoảng 10% tổng số loài chim và thú trên hành tinh. Tuy nhiên, trước áp lực từ tăng trưởng dân số, quá trình đô thị hóa nhanh chóng và tình trạng khai thác quá mức, ĐDSH nước ta đang đối mặt với nguy cơ suy thoái nghiêm trọng với hàng trăm loài thực vật và động vật nằm trong Sách Đỏ và Danh lục đỏ IUCN.
Bán đảo Sơn Trà (TP. Đà Nẵng) giữ vai trò như một "lá phổi xanh" và bức tường thành sinh thái che chắn cho đô thị loại I Đà Nẵng, sở hữu hệ thực vật nhiệt đới giao thoa giữa các luồng di lưu thực vật bản địa Bắc Việt Nam - Trung Hoa, luồng Malaysia - Indonesia và Ấn Độ - Myanmar. Mặc dù đã có các khảo sát định tính về thành phần loài, song một khoảng trống nghiên cứu (research gap) lớn tồn tại suốt nhiều thập kỷ: sự thiếu hụt các đánh giá định lượng chuyên sâu về cấu trúc quần xã thực vật thân gỗ thông qua các chỉ số sinh thái học toán học. Hầu hết dữ liệu trước đây dừng lại ở danh lục kiểm kê sơ bộ, chưa phân tích mối tương quan cạnh tranh ổ sinh thái (ecological niche), kiểu phân bố không gian và trật tự ưu thế của các loài chủ chốt.
Trong bối cảnh Luật Đa dạng sinh học ngày càng đòi hỏi các chương trình bảo tồn loài mục tiêu phải dựa trên bằng chứng khoa học cụ thể, cùng với sự bùng nổ của hạ tầng du lịch nghỉ dưỡng và các nguy cơ sinh học (như sự phát triển bùng phát của dây leo Bìm bìm), việc định lượng hóa ĐDSH là nhiệm vụ cấp thiết. Kết quả nghiên cứu không chỉ thiết lập hệ quy chiếu sinh thái chuẩn xác cho Bán đảo Sơn Trà mà còn mang lại mô hình tham chiếu thực tiễn cho hệ thống rừng đặc dụng duyên hải miền Trung - Tây Nguyên.
Methodology và approach
Nghiên cứu ứng dụng phương pháp tiếp cận tổ hợp giữa sinh thái học thực địa (Phyto-sociological survey), phân tích thống kê toán học và ứng dụng công nghệ viễn thám không gian:
1. Thiết kế lấy mẫu thực địa (Sampling Design)
Hệ thống 12 ô tiêu chuẩn (OTC) với diện tích $500\text{ m}^2/\text{ô}$ được thiết kế ngẫu nhiên tầng hóa trên 5 trạng thái sinh cảnh tiêu biểu: rừng tự nhiên (OTC 1, 2, 3, 4, 8), rừng tự nhiên Chò đen (OTC 9, 12), rừng trồng (OTC 6, 7), trảng cỏ (OTC 5) và đất trống thoái hóa (OTC 10, 11).
Trong mỗi OTC, cấu trúc phân tầng được chia thành các ô lồng nhau (nested sub-plots):
- Ô A ($10\text{m} \times 10\text{m} = 100\text{ m}^2$): Điều tra toàn bộ tầng cây gỗ có đường kính ngang ngực ($D_{1.3}$) $\ge 10\text{ cm}$.
- Ô B ($5\text{m} \times 5\text{m} = 25\text{ m}^2$): Đo đếm tầng cây gỗ nhỏ có $D_{1.3} < 10\text{ cm}$.
- Ô C ($2\text{m} \times 2\text{m} = 4\text{ m}^2$): Đánh giá lớp cây tái sinh có chiều cao từ $0{,}3\text{ m}$ đến $1{,}3\text{ m}$, phân loại tái sinh chồi hay hạt.
2. Bộ chỉ số định lượng sinh thái (Quantitative Indices)
- Chỉ số giá trị quan trọng ($\text{IVI}$ - Importance Value Index): Tính toán theo công thức của Mishra (1968):
$$\text{IVI} = \text{RD} + \text{RF} + \text{RBA}$$
Trong đó: $\text{RD}$ là Mật độ tương đối (%), $\text{RF}$ là Tần suất xuất hiện tương đối (%), và $\text{RBA}$ là Tổng tiết diện thân tương đối (%).
- Dạng phân bố không gian ($A/F$): Tỷ số giữa Độ phong phú ($A$) và Tần suất ($F$). Nhận diện kiểu phân bố: Liên tục ($A/F < 0{,}025$), Ngẫu nhiên ($0{,}025 \le A/F \le 0{,}05$), hoặc Cụm/Lan truyền ($A/F > 0{,}05$).
- Chỉ số đa dạng loài Shannon-Wiener ($H'$):
$$H' = -\sum_{i=1}^{n} \left(\frac{N_i}{N}\right) \log_2 \left(\frac{N_i}{N}\right)$$
- Chỉ số mức độ chiếm ưu thế Simpson ($C_d$):
$$C_d = \sum_{i=1}^{n} \left(\frac{N_i}{N}\right)^2$$
- Chỉ số tương đồng Sorensen ($SI$):
$$SI = \frac{2C}{A + B}$$
- Phân tích đường cong Đa dạng - Ưu thế ($D\text{-}D\text{ curve}$): Đánh giá sự chia sẻ ổ sinh thái thông qua các mô hình phân phối chuẩn logarit (Log-normal series) và chuỗi hình học (Geometric series).
Phát hiện chính
Quá trình điều tra và tính toán toán học đã xác lập 4 phát hiện cốt lõi về hệ sinh thái thực vật thân gỗ tại BTTN Sơn Trà:
1. Tính phong phú taxon và sự vượt trội của các họ ưu thế
Trong tổng số $6.000\text{ m}^2$ diện tích khảo sát thuộc 12 OTC, đã ghi nhận 96 loài thực vật thân gỗ thuộc 46 họ. Trong đó, các họ chiếm ưu thế tuyệt đối về mật độ cá thể là Họ Thầu dầu (Euphorbiaceae) và Họ Long não (Lauraceae) cùng chiếm $9{,}4%$; Họ Sim (Myrtaceae) chiếm $7{,}3%$; các họ Dẻ (Fagaceae), Dầu (Dipterocarpaceae), Dâu tằm (Moraceae) và Na (Annonaceae) cùng đóng góp $5{,}2%$.
2. Trật tự ưu thế IVI và sự phân bố ổ sinh thái
Loài Chò đen (Parashorea stellata Kurz) khẳng định vai trò là loài ưu thế sinh thái quan trọng nhất với $\text{IVI} = 35{,}38$, mật độ tương đối $\text{RD} = 14{,}6%$, tần suất tương đối $\text{RF} = 25%$, và mật độ bình quân đạt $58{,}33\text{ cá thể/OTC}$. Xếp sau là Xoài cuống dài ($\text{IVI} = 14{,}42$) và Mộc ($\text{IVI} = 10{,}77$). Cấu trúc IVI cho thấy quần xã chưa bị độc tôn bởi một loài duy nhất (tổng giá trị chia sẻ giữa các tầng tán), tạo nên tính ổn định tương đối của rừng nhiệt đới thường xanh.
3. Sự biến thiên mạnh mẽ của chỉ số đa dạng loài ($H'$) và độ tập trung ($C_d$)
- Chỉ số Shannon ($H'$): Đạt giá trị rất cao tại sinh cảnh Rừng tự nhiên và Rừng tự nhiên Chò đen ($H' = 3{,}61 - 4{,}76$), tương đương với các khu rừng mưa nhiệt đới giàu tính đa dạng tại Đông Nam Á. Tuy nhiên, chỉ số này sụt giảm nghiêm trọng tại các sinh cảnh chịu tác động: Đất trống ($2{,}62 - 2{,}82$), Rừng trồng ($1{,}86 - 2{,}60$), và thấp nhất ở Trảng cỏ ($1{,}62$).
- Chỉ số Simpson ($C_d$): Biến thiên từ $0{,}051$ đến $0{,}499$. Giá trị $C_d$ tăng cao ở các trạng thái thoái hóa phản ánh sự nghèo nàn về thành phần loài và sự thống trị cục bộ của một vài loài cây tiên phong chịu hạn, báo hiệu suy giảm đa dạng sinh thái.
Chỉ số H' (Shannon-Wiener) qua các Sinh cảnh
4. Dạng phân bố không gian ($A/F$) và phân tích đường cong $D\text{-}D$
Phần lớn các loài cây gỗ tại Sơn Trà có tỷ số $A/F > 0{,}05$, xác nhận dạng phân bố lan truyền/theo cụm (contagious pattern) – hình thái thích nghi điển hình của thực vật rừng tự nhiên giúp tăng khả năng thụ phấn và chống chịu gió bão ven biển.
Đường cong $D\text{-}D$ tại sinh cảnh rừng tự nhiên có dạng chuỗi logarit - bình thường (Log-normal distribution series), chứng minh các loài cùng chia sẻ nguồn sống hài hòa trong ổ sinh thái. Ngược lại, tại sinh cảnh rừng trồng và trảng cỏ, đường cong biến dạng thành dạng chuỗi hình học (Geometric distribution series) dốc đứng, phản ánh mức độ cạnh tranh thấp, tính đa dạng suy kiệt và ổ sinh thái còn nhiều khoảng trống dễ bị xâm thực bởi thực vật ngoại lai.
Đóng góp khoa học
Nghiên cứu mang lại những bước tiến quan trọng cả về mặt lý luận học thuật lẫn giá trị thực tiễn:
- Về mặt lý luận và học thuật: Đây là công trình đầu tiên ứng dụng tổ hợp đầy đủ các chỉ số định lượng sinh thái ($\text{IVI}, H', C_d, SI, A/F, D\text{-}D\text{ curve}$) để giải phẫu cấu trúc thảm thực vật thân gỗ tại BTTN Sơn Trà. Đề tài mở rộng phương pháp luận nghiên cứu cấu trúc rừng bán đảo nhiệt đới ven biển, bổ sung tư liệu giá trị cho sinh thái học rừng Việt Nam.
- Công bố khoa học quốc tế và trong nước: Các phát hiện của đề tài đã được công bố qua 02 bài báo khoa học chuẩn mực: 01 bài trên Tạp chí Khoa học Lâm nghiệp (Viện Khoa học Lâm nghiệp Việt Nam) và 01 bài trên Tạp chí Khoa học Lâm nghiệp của Trường Đại học Tổng hợp Lâm nghiệp Quốc gia Saint-Petersburg (Liên bang Nga).
- Ứng dụng thực tiễn trong công tác quản lý: Dữ liệu nghiên cứu chuyển giao trực tiếp cho Ban Quản lý Khu BTTN Sơn Trà và Chi cục Kiểm lâm Đà Nẵng, hỗ trợ:
- Xác định ranh giới và bảo vệ nghiêm ngặt 2.593,1 ha rừng tự nhiên - nơi tập trung các loài có giá trị bảo tồn cao như Chò đen, Dầu rái, Chò nâu.
- Định hướng chiến lược xử lý sinh học ngăn chặn loài dây leo Bìm bìm xâm lấn tại các sinh cảnh rừng phục hồi ($IIA$) và rừng nghèo kiệt.
- Lựa chọn các loài bản địa ưu thế ($\text{IVI}$ cao) như Chò đen, Sao đen, Dầu con quay để xúc tiến tái sinh tự nhiên và làm giàu rừng thay thế rừng trồng thuần loài keo, bạch đàn.
Đối tượng quan tâm
Báo cáo nghiên cứu là tài liệu tham khảo khoa học hữu ích cho nhiều nhóm chuyên môn:
- Các nhà nghiên cứu và viện chuyên ngành: Tiếp cận bộ dữ liệu thực địa tin cậy để so sánh, đối chiếu diễn thế sinh thái rừng ven biển qua các thời kỳ và các vùng sinh thái tương đồng.
- Nhà hoạch định chính sách và Quản lý bảo tồn: Sở Tài nguyên & Môi trường, Ban Quản lý rừng đặc dụng áp dụng kết quả để thiết kế các trạm kiểm lâm cơ động, phân vùng bảo vệ nghiêm ngặt và lập phương án chia sẻ lợi ích sinh thái với cộng đồng dân cư Thọ Quang.
- Đơn vị phát triển du lịch bền vững: Cung cấp thông tin sinh thái phục vụ công tác thiết kế các tuyến du lịch diễn giải môi trường, nâng cao ý thức bảo vệ tài nguyên rừng của du khách.
- Cộng đồng học thuật, sinh viên: Nguồn học liệu tham khảo chuẩn mực trong các học phần Sinh thái học thực vật, Quản lý môi trường và Đánh giá tác động môi trường (ĐTM).
FAQ (Câu hỏi thường gặp)
1. Phát hiện quan trọng nhất của công trình nghiên cứu này là gì?
Trả lời: Phát hiện then chốt là việc xác lập định lượng loài Chò đen (Parashorea stellata) là loài ưu thế tuyệt đối ($\text{IVI} = 35{,}38$) trong cấu trúc rừng tự nhiên Sơn Trà, cùng sự phân hóa sâu sắc của chỉ số đa dạng Shannon ($H'$) từ mức rất cao ở rừng tự nhiên ($4{,}76$) xuống mức suy kiệt ở trảng cỏ ($1{,}62$). Điều này minh chứng rừng tự nhiên Sơn Trà có giá trị bảo tồn sinh thái ngoại hạng nhưng cực kỳ nhạy cảm trước các tác động nhân sinh.
2. Phương pháp nghiên cứu trong đề tài có điểm gì đột phá so với các khảo sát trước đây?
Trả lời: Đề tài đã chuyển hóa phương thức khảo sát truyền thống (thuần túy lập danh lục kiểm kê) sang đánh giá định lượng toán sinh thái toàn diện. Việc kết hợp hệ thống ô tiêu chuẩn tầng hóa (Quadrat lồng ghép $100\text{ m}^2 - 25\text{ m}^2 - 4\text{ m}^2$), định vị vệ tinh VN-2000 và tính toán đồng thời 5 nhóm chỉ số ($\text{IVI}, H', C_d, SI, A/F$) cùng đường cong $D\text{-}D$ mang lại độ chuẩn xác và tin cậy khoa học cao.
3. Kết quả nghiên cứu có thể tổng quát hóa (generalize) cho các khu vực khác không?
Trả lời: Phương pháp luận và quy trình tính toán chỉ số hoàn toàn có thể nhân rộng chuẩn hóa cho toàn bộ hệ sinh thái rừng nhiệt đới gió mùa tại Việt Nam. Về mặt dữ liệu thực chứng, kết quả có giá trị tham chiếu cao đối với hệ thống rừng đặc dụng ven biển và hải đảo thuộc khu vực Duyên hải Miền Trung và Tây Nguyên có điều kiện lập địa trên nền đá Granit tương đồng.
4. Những định hướng nghiên cứu tiếp theo (next steps) cần triển khai là gì?
Trả lời: Cần thiết lập các ô định vị sinh thái cố định dài hạn (Permanent Sampling Plots) để quan trắc diễn thế rừng nhiều năm; mở rộng điều tra đa dạng di truyền các loài nguy cấp có tên trong Sách Đỏ như Chò đen, Gụ lau, Re hương; và nghiên cứu cơ chế tương tác sinh thái giữa thảm thực vật thân gỗ với quần thể linh trưởng đặc hữu Voọc chà vá chân nâu (Pygathrix nemaeus).
5. Kết quả này ứng dụng cụ thể như thế nào trong thực tiễn lâm sinh?
Trả lời: Dữ liệu $\text{IVI}$ và tầng tái sinh là cơ sở để chọn lọc các loài cây bản địa có sức sống và độ che phủ cao (Chò đen, Sao đen, Dầu rái) phục vụ trồng phục hồi rừng sau khai thác và thay thế dần diện tích rừng trồng keo lai ngoại lai, đồng thời khoanh vùng ưu tiên xử lý các điểm nóng xâm lấn của dây leo Bìm bìm.
Kết luận
Công trình nghiên cứu của nhóm tác giả do TS. Phạm Thị Kim Thoa chủ nhiệm đã hoàn thành xuất sắc mục tiêu định lượng hóa bức tranh đa dạng sinh học thực vật thân gỗ tại Khu Bảo tồn Thiên nhiên Sơn Trà. Bằng hệ phương pháp khoa học chuẩn mực, đề tài đã chứng minh Sơn Trà là một kho tàng ĐDSH quý báu với 96 loài cây gỗ thuộc 46 họ, sở hữu cấu trúc rừng tự nhiên giàu có nhưng đang đứng trước nguy cơ phân mảnh sinh cảnh do áp lực phát triển kinh tế - du lịch.
Để giữ vững tính toàn vẹn sinh thái của "viên ngọc xanh" Đà Nẵng, việc cấp bách hiện nay là thực thi nghiêm ngặt quy chế bảo tồn phân khu 2.593,1 ha, ứng dụng các giải pháp lâm sinh phục hồi rừng bản địa dựa trên dữ liệu $\text{IVI}$, và tăng cường giám sát khoa học định kỳ. Các cơ quan quản lý, viện nghiên cứu và cộng đồng cần chung tay hành động vì sự phát triển bền vững của tài nguyên thiên nhiên Sơn Trà.