Tổng quan nghiên cứu

Khu Bảo tồn Thiên nhiên Bán đảo Sơn Trà thuộc thành phố Đà Nẵng sở hữu vị trí địa lý đắc địa khi cách trung tâm thành phố 10 km về phía Đông Bắc, đóng vai trò như lá phổi xanh điều hòa không khí và bức bình phong sinh thái vững chắc chắn bão gió ven biển. Về mặt địa lý sinh vật học, bán đảo Sơn Trà là điểm giao thoa độc đáo giữa hai luồng thực vật Bắc và Nam, tạo nên một hệ sinh thái rừng nhiệt đới ven biển phong phú với 20 con suối tự nhiên quanh năm cung cấp nước ngọt cho cư dân bán đảo. Tổng diện tích đất lâm nghiệp của khu bảo tồn đạt 4.195,70 ha, chiếm 69,1% diện tích tự nhiên của toàn quận Sơn Trà. Nơi đây từng ghi nhận 985 loài thực vật bậc cao có mạch thuộc 483 chi và 143 họ, trong đó có 22 loài thực vật quý hiếm nằm trong Sách đỏ Việt Nam như Gụ lau, Chò chai, Re hương và Cẩm lai Bà Rịa.

Mặc dù giữ vị trí chiến lược, tính đa dạng thực vật thân gỗ tại Sơn Trà đang đối mặt với nhiều áp lực nghiêm trọng từ sự xâm lấn của các loài dây leo, biến đổi thảm thực vật sau khai thác và tác động mạnh mẽ từ tốc độ đô thị hóa, du lịch sinh thái. Để thiết lập cơ sở dữ liệu định lượng phục vụ công tác quản lý, nghiên cứu này được triển khai nhằm đánh giá toàn diện cấu trúc thành phần loài, xác định các chỉ số sinh thái học cốt lõi như mức độ ưu thế, độ phong phú và mối quan hệ phân cụ giữa các quần thể thực vật thân gỗ qua các dạng sinh cảnh rừng khác nhau. Luận văn cung cấp hệ thống số liệu khoa học tin cậy, đóng góp các chỉ số sinh thái quan trọng giúp cơ quan quản lý bảo tồn bền vững nguồn gen thực vật nguy cấp và tối ưu hóa quy hoạch cảnh quan thiên nhiên tại Sơn Trà.

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 sinh thái học quần xã và Công ước về Đa dạng sinh học (CBD 1992), tiếp cận đa dạng sinh học ở cả ba cấp độ: đa dạng di truyền, đa dạng loài và đa dạng hệ sinh thái. Nghiên cứu tập trung khai thác các bậc cấu trúc đa dạng alpha (tính đa dạng trong cùng một sinh cảnh) và đa dạng beta (mức độ thay đổi thành phần loài giữa các sinh cảnh liền kề).

Để định lượng hóa cấu trúc quần xã thực vật thân gỗ, luận văn vận dụng hệ thống các chỉ số toán sinh thái tiêu chuẩn quốc tế:

  • Chỉ số mức độ phong phú loài Margalef ($d = (S - 1) / \ln N$): Đánh giá độ giàu loài dựa trên tương quan giữa tổng số loài ($S$) và tổng số cá thể ($N$).
  • Chỉ số đa dạng Shannon-Wiener ($H = - \sum p_i \log_2 p_i$): Phản ánh lượng thông tin sinh thái và độ bất định của cấu trúc loài trong quần xã.
  • Chỉ số độ đồng đều Pielou ($J' = H / \log_2 S$): Đo lường mức độ đồng đều trong phân bố cá thể giữa các loài, dao động trong khoảng từ 0 đến 1.
  • Chỉ số mức độ chiếm ưu thế Simpson ($Cd = \sum (n_i / N)^2$): Biểu thị xác suất hai cá thể chọn ngẫu nhiên thuộc cùng một loài, làm nổi bật vai trò của các loài chủ chốt.
  • Chỉ số giá trị quan trọng (Importance Value Index - $IVI = \text{Mật độ tương đối} + \text{Tần suất tương đối} + \text{Độ che phủ tương đối}$): Tổng điểm đạt 300, xác định vị thế sinh thái và độ ưu thế tuyệt đối của từng loài cây gỗ.
  • Tỷ số $A/F$ (Abundance/Frequency): Phân loại quy luật phân bố không gian của các loài thành phân bố ngẫu nhiên ($A/F < 0,025$), phân bố chuẩn ($0,025 \le A/F \le 0,05$) hoặc phân bố cụm ($A/F > 0,05$).

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

Nguồn dữ liệu thực nghiệm được thu thập thông qua phương pháp điều tra lâm học kết hợp hệ thống ô tiêu chuẩn (OTC) định vị bằng thiết bị GPS cầm tay và bản đồ số MapInfo.

Về cỡ mẫu và thiết kế điều tra: Tổng cộng 12 ô tiêu chuẩn được thiết lập trên thực địa với diện tích 500 m² mỗi ô (kích thước chuẩn $10\text{ m} \times 50\text{ m}$ hoặc tương đương), nâng tổng diện tích điều tra chuyên sâu lên 6.000 m². Cấu trúc mỗi ô tiêu chuẩn 500 m² được phân tầng thành 3 cấp tiểu ô:

  • Ô cấp A: Gồm 5 tiểu ô kích thước $10\text{ m} \times 10\text{ m}$ (100 m²) để đo đếm toàn bộ cây gỗ có đường kính thân ở vị trí $D_{1.3\text{m}} > 10\text{ cm}$.
  • Ô cấp B: Gồm 5 tiểu ô kích thước $5\text{ m} \times 5\text{ m}$ (25 m²) đo đếm cây bụi và cây gỗ nhỏ có đường kính $D_{1.3\text{m}} < 10\text{ cm}$ và chiều cao lớn hơn 1,3 m.
  • Ô cấp C: Các tiểu ô kích thước $2\text{ m} \times 2\text{ m}$ (4 m²) đo đếm cây tái sinh có chiều cao từ 0,3 m đến dưới 1,3 m, phân loại rõ tái sinh chồi hay hạt.

Phương pháp chọn mẫu: Nghiên cứu áp dụng phương pháp lấy mẫu phân tầng ngẫu nhiên qua 5 dạng sinh cảnh đại diện tại Sơn Trà: sinh cảnh rừng tự nhiên (OTC 1, 2, 3, 4, 8), sinh cảnh rừng tự nhiên Chò đen (OTC 9, 12), sinh cảnh rừng trồng (OTC 6, 7), sinh cảnh trảng cỏ (OTC 5) và sinh cảnh đất trống phục hồi (OTC 10, 11).

Lý do lựa chọn phương pháp: Thiết kế ô tiêu chuẩn phân tầng cho phép thu thập dữ liệu cấu trúc thảm thực vật từ tầng vượt tán đến tầng tái sinh dưới tán mà không bỏ sót sinh khối. Việc kết hợp phân tích đường cong Đa dạng - Ưu thế (D-D curve) và phân tích cụm (Cluster analysis) theo ma trận tương đồng Bray-Curtis giúp lượng hóa chính xác cấu trúc sinh thái, loại bỏ tính chủ quan trong phân loại học truyền thống.

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

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

Thứ nhất, mức độ đa dạng loài và số lượng cá thể có sự biến động sâu sắc giữa các dạng sinh cảnh. Thành phần loài ($S$) trong 12 ô tiêu chuẩn dao động từ 4 đến 39 loài/OTC (trung bình 20 loài/OTC). Mật độ cá thể cây gỗ ($N$) biến thiên mạnh từ 11 đến 190 cá thể/500 m² (trung bình 88 cá thể/OTC). Trạng thái rừng tự nhiên đạt chỉ số đa dạng cao nhất tại OTC 1 với $S = 39$ loài, $N = 135$ cá thể, chỉ số $d = 7,75$ và $H = 4,76$. Ngược lại, sinh cảnh trảng cỏ tại OTC 5 ghi nhận độ giàu loài thấp nhất với $S = 4$ loài, $N = 11$ cá thể, $d = 1,25$ và $H = 1,62$.

Thứ hai, xác lập trật tự chiếm ưu thế của các loài thực vật thân gỗ qua chỉ số giá trị quan trọng ($IVI$). Xét trên toàn bộ khu hệ khảo sát, loài Chò đen (Parashorea stellata) giữ vị trí ưu thế số một với chỉ số $IVI = 35,38/300$. Đứng ở các vị trí tiếp theo là loài Xoài cuống dài (Mangifera laurina) với $IVI = 14,42$ và Mộc (Planchonella obovata) với $IVI = 10,77$. Ở các trạng thái rừng tự nhiên hỗn giao, mức độ ưu thế phân bố khá đồng đều giữa các loài tầng nhỡ và tầng ưu thế sinh thái.

Thứ ba, sự sai khác cấu trúc ưu thế cục bộ tại sinh cảnh rừng trồng và rừng Chò đen thuần loại. Tại sinh cảnh rừng tự nhiên Chò đen (OTC 9 và OTC 12), loài Chò đen chiếm ưu thế áp đảo với giá trị $IVI$ đạt tới 131,0/300 trong tổng số 300 điểm, dẫn đến chỉ số ưu thế Simpson tăng vọt lên mức $Cd = 0,493$ và độ đồng đều giảm xuống $J' = 0,47$. Tương tự, tại sinh cảnh rừng trồng (OTC 6 và OTC 7), hai loài chiếm vị trí thống trị tuyệt đối là Bạch đàn trắng (Eucalyptus camaldulensis) với $IVI = 166,3/300$ và Hà nu (Ixonanthes reticulata) với $IVI = 104,7/300$.

Thứ tư, xác định quy luật phân bố không gian và mức độ tương đồng giữa các loài. Kết quả phân tích tỷ số $A/F$ chỉ ra rằng có 93/94 loài cây gỗ khảo sát có dạng phân bố cụm ($A/F > 0,05$), điển hình là Chò đen có giá trị $A/F = 2,333$, chứng tỏ điều kiện lập địa và môi trường vi khí hậu tại Sơn Trà có tính ổn định cao. Duy nhất loài Cọ mai nhám (Colona thorelii) thể hiện dạng phân bố ngẫu nhiên với $A/F = 0,03$. Khi tiến hành phân tích cụm ở mức tương đồng 20%, hệ thực vật được chia thành 9 nhóm loài riêng biệt; ở mức tương đồng 50%, các loài được gom thành 18 nhóm cấu trúc sinh thái rõ rệt.

Thảo luận kết quả

Đường cong Đa dạng - Ưu thế (Diversity - Dominance curve) của thảm thực vật cây gỗ tại sinh cảnh rừng tự nhiên Sơn Trà tuân theo mô hình phân bố Log-chuẩn (Log-normal distribution series). Điều này chứng minh rằng tại các sinh cảnh rừng tự nhiên ít bị xáo trộn, thảm thực vật cây gỗ không có hiện tượng một vài loài lấn át hoàn toàn nguồn tài nguyên, tính cạnh tranh sinh thái diễn ra hài hòa và thảm rừng đạt độ bão hòa cao.

Trái lại, đường cong D-D tại sinh cảnh rừng tự nhiên Chò đen, rừng trồng bạch đàn và trảng cỏ lại thể hiện dạng phân bố Hình học (Geometric distribution series) với đường dốc gần như thẳng đứng. Đồ thị này phản ánh sự suy giảm tính đa dạng sinh học cục bộ khi một loài duy nhất chiếm giữ phần lớn sinh khối và không gian dinh dưỡng, làm giảm chỉ số độ đồng đều $J'$ xuống dưới mức 0,50 và cản trở sự sinh trưởng của các loài cây bản địa dưới tán.

Dữ liệu nghiên cứu được biểu diễn trực quan qua hệ thống biểu đồ đường phân tán thể hiện biến động chỉ số $H$ và $Cd$ trên cùng trục hoành 12 ô tiêu chuẩn, kết hợp bảng ma trận phân cụ sinh thái học. Cách trình bày này giúp làm nổi bật mối tương quan nghịch giữa chỉ số đa dạng Shannon-Wiener và chỉ số ưu thế Simpson, tạo công cụ trực quan hỗ trợ các nhà quản lý nhận diện nhanh chóng những khu vực rừng đang bị thoái hóa hoặc mất cân bằng sinh thái.

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

Dựa trên các kết quả phân tích định lượng về cấu trúc thảm thực vật thân gỗ, nghiên cứu đưa ra 4 nhóm giải pháp kỹ thuật và quản lý cụ thể:

Thứ nhất, thực hiện điều tiết cấu trúc tán rừng và tỉa thưa sinh thái. Chi cục Kiểm lâm thành phố Đà Nẵng phối hợp cùng Ban Quản lý Khu bảo tồn tiến hành các biện pháp lâm sinh điều chỉnh mật độ tại 4.195,70 ha đất lâm nghiệp, đặc biệt là các khu vực rừng trồng bạch đàn (nơi Bạch đàn trắng chiếm $IVI = 166,3/300$) và các ô tiêu chuẩn Chò đen có mật độ quá dày. Mục tiêu là mở các khoảng trống tán nhỏ nhằm hạ tỷ lệ che phủ của loài chiếm ưu thế xuống dưới 35%, kích thích quá trình tái sinh tự nhiên của các loài cây gỗ bản địa có giá trị cao. Thời gian triển khai từ quý I năm 2024 đến hết năm 2026.

Thứ hai, thiết lập mạng lưới ô định vị giám sát đa dạng sinh thái dài hạn. Ban Quản lý Khu bảo tồn Thiên nhiên Sơn Trà cần duy trì và chuẩn hóa 12 ô tiêu chuẩn cố định diện tích 500 m² đã thiết lập trong luận văn. Định kỳ 12 tháng một lần, lực lượng kỹ thuật đo đếm lại các chỉ số $S, N, H, Cd$ và tỷ lệ cây tái sinh để theo dõi biến động sinh thái, đảm bảo duy trì chỉ số đa dạng loài $H \ge 3,50$ tại các phân khu phục hồi sinh thái.

Thứ ba, khoanh vùng bảo vệ nguyên vị (In-situ) đối với 22 loài thực vật quý hiếm nằm trong Sách đỏ Việt Nam. Sở Nông nghiệp và Phát triển Nông thôn thành phố Đà Nẵng cần ban hành đề án bảo tồn nguồn gen đối với các loài như Gụ lau, Chò chai, Re hương và Cẩm lai Bà Rịa. Đặt mục tiêu nâng tỷ lệ sống sót của cây con tái sinh tự nhiên thuộc nhóm nguy cấp này lên trên 60% trong giai đoạn 2024–2027 thông qua các biện pháp dọn dây leo xâm lấn Bim Bim và thiết lập đai bảo vệ cục bộ.

Thứ tư, quy hoạch phân khu du lịch sinh thái gắn liền với bảo vệ 20 lưu vực suối tự nhiên. Ủy ban nhân dân quận Sơn Trà và các đơn vị lữ hành phải tuân thủ nghiêm ngặt ranh giới bảo vệ nghiêm ngặt, nghiêm cấm xây dựng hạ tầng bê tông hóa trong bán kính 100 m tính từ các dòng suối chính như Suối Đá, Suối Heo nhằm bảo toàn độ ẩm và điều kiện vi khí hậu của thảm thực vật thân gỗ.

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

Luận văn thạc sĩ này là tài liệu chuyên khảo mang giá trị học thuật và thực tiễn cao cho 4 nhóm đối tượng trọng tâm:

Thứ nhất, Ban Quản lý Khu bảo tồn Thiên nhiên Sơn Trà và Hạt Kiểm lâm liên quận Sơn Trà - Ngũ Hành Sơn: Sử dụng toàn bộ hệ thống số liệu định lượng về 12 ô tiêu chuẩn và danh lục cây gỗ ưu thế để xây dựng phương án quản lý rừng bền vững trên diện tích 4.195,70 ha rừng đặc dụng.

Thứ hai, các nhà khoa học, nghiên cứu sinh và học viên cao học chuyên ngành Sinh thái học, Quản lý Tài nguyên Rừng: Kế thừa khung phương pháp luận về tính toán chỉ số $IVI$, phân tích chuỗi phân bố đường cong D-D và xử lý phân cụ tương đồng loài bằng công cụ toán sinh thái ứng dụng cho các hệ sinh thái rừng nhiệt đới ven biển.

Thứ ba, các cơ quan quản lý nhà nước như Sở Nông nghiệp và Phát triển Nông thôn, Sở Tài nguyên và Môi trường thành phố Đà Nẵng: Ứng dụng các chỉ số đa dạng loài ($H = 3,23; Cd = 0,209$) làm tiêu chí nền tảng để thẩm định báo cáo đánh giá tác động môi trường (ĐTM) cho các dự án phát triển kinh tế, giao thông và du lịch tại bán đảo.

Thứ tư, các tổ chức phi chính phủ và quỹ bảo tồn quốc tế (như WWF, IUCN, GreenViet): Tham khảo dữ liệu phân bố của 22 loài thực vật Sách đỏ để xây dựng các đề xuất tài trợ dự án phục hồi rừng đầu nguồn và giáo dục môi trường cộng đồng tại miền Trung Việt Nam.

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

Chỉ số Shannon-Wiener ($H$) trung bình của thảm thực vật thân gỗ tại Sơn Trà đạt bao nhiêu và có ý nghĩa gì? Giá trị $H$ trung bình trên toàn bộ 12 ô tiêu chuẩn đạt 3,23 (tính theo logarit cơ số 2), biến thiên từ 1,62 ở trảng cỏ đến 4,76 ở rừng tự nhiên. Giá trị này khẳng định hệ thực vật thân gỗ Sơn Trà có cấu trúc thành phần loài phong phú và hệ sinh thái đang trong giai đoạn phát triển ổn định.

Loài thực vật thân gỗ nào chiếm giá trị ưu thế sinh thái cao nhất tại Khu Bảo tồn Thiên nhiên Sơn Trà? Loài Chò đen (Parashorea stellata) chiếm ưu thế tuyệt đối toàn khu hệ với chỉ số $IVI = 35,38/300$. Tại các sinh cảnh rừng Chò đen đặc trưng, loài này có thể đạt $IVI$ lên tới 131,0/300, đóng vai trò cây ưu thế tầng vượt tán định hình sinh cảnh.

Đường cong Đa dạng - Ưu thế (D-D curve) dạng Log-chuẩn tại rừng tự nhiên Sơn Trà phản ánh điều gì? Đường cong phân bố dạng Log-chuẩn tại 6 ô rừng tự nhiên phản ánh một quần xã sinh thái có tính cạnh tranh tài nguyên đồng đều, không có loài nào chiếm ưu thế tuyệt đối làm lấn át các loài khác, biểu thị cấu trúc rừng đã đạt trạng thái cân bằng sinh thái tự nhiên.

Quy mô và phương pháp bố trí ô tiêu chuẩn trong luận văn được thực hiện như thế nào? Nghiên cứu bố trí 12 ô tiêu chuẩn với diện tích 500 m² mỗi ô (tổng diện tích 6.000 m²). Trong mỗi ô lớn lồng ghép 5 tiểu ô $10\text{ m} \times 10\text{ m}$ đo cây gỗ lớn, 5 tiểu ô $5\text{ m} \times 5\text{ m}$ đo cây bụi và các ô $2\text{ m} \times 2\text{ m}$ đo cây tái sinh.

Tại sao cần thực hiện các biện pháp lâm sinh điều chỉnh đối với rừng trồng Bạch đàn trắng tại Sơn Trà? Tại sinh cảnh rừng trồng, loài Bạch đàn trắng chiếm $IVI = 166,3/300$, tạo ra đường phân bố hình học dốc đứng. Cần tỉa thưa để giảm thiểu sự lấn át dinh dưỡng và ánh sáng, tạo không gian cho các loài cây gỗ bản địa tái sinh dưới tán.

Kết luận

  • Luận văn đã định lượng hóa toàn diện cấu trúc đa dạng thực vật thân gỗ tại Bán đảo Sơn Trà thông qua hệ thống 12 ô tiêu chuẩn diện tích 500 m²/ô đại diện cho 5 sinh cảnh sinh thái.
  • Thiết lập bộ chỉ số sinh thái học tiêu chuẩn với số loài trung bình đạt 20 loài/ô, mật độ trung bình 88 cá thể/ô, chỉ số Shannon-Wiener trung bình $H = 3,23$ và chỉ số Simpson $Cd = 0,209$.
  • Nhận diện Chò đen là loài ưu thế sinh thái quan trọng nhất với $IVI = 35,38/300$, đồng thời chỉ ra nguy cơ mất cân bằng đa dạng loài tại các khu vực rừng trồng bạch đàn với $IVI = 166,3/300$.
  • Chứng minh tính ổn định sinh thái cao của khu bảo tồn khi có tới 93/94 loài cây gỗ phân bố không gian theo dạng cụm ($A/F > 0,05$) và xác lập 9 nhóm phân cụ loài ở mức tương đồng 20%.
  • Cung cấp luận cứ khoa học thực chứng để khoanh vùng bảo vệ nghiêm ngặt 22 loài thực vật quý hiếm nằm trong Sách đỏ Việt Nam trên tổng diện tích 4.195,70 ha đất lâm nghiệp.

Nghiên cứu đã đóng góp cơ sở dữ liệu toán sinh thái chuẩn mực, phục vụ đắc lực cho công tác quản lý tài nguyên và định hướng bảo tồn đa dạng sinh thái tại Bán đảo Sơn Trà. Để tiếp tục duy trì tính bền vững, các chương trình quan trắc sinh thái học định kỳ cần được duy trì liên tục trong giai đoạn 2024–2030 trên toàn bộ các phân khu bảo tồn. Hãy tham khảo và trích dẫn chi tiết toàn văn luận văn thạc sĩ khoa học để ứng dụng hiệu quả các phương pháp phân tích đa biến và giải pháp kỹ thuật lâm sinh vào thực tiễn quản lý tài nguyên thiên nhiên.