Tổng quan về luận án

Vườn Quốc gia (VQG) Bù Gia Mập thuộc tỉnh Bình Phước sở hữu vị trí địa lý - sinh thái đặc biệt quan trọng, nằm trong vùng chuyển tiếp giữa cao nguyên Nam Tây Nguyên và vùng đồng bằng Đông Nam Bộ, với tổng diện tích tự nhiên 25.926 ha (phân khu bảo vệ nghiêm ngặt 18.269 ha, phân khu phục hồi sinh thái 7.557 ha và phân khu hành chính dịch vụ 100 ha). Đây là nơi lưu giữ diện tích rừng nguyên sinh và thứ sinh giàu có bậc nhất Đông Nam Bộ, đóng vai trò xung yếu phòng hộ đầu nguồn cho hệ thống thủy điện Thác Mơ, Cần Đơn và Sork Phu Miêng.

Khoảng trống nghiên cứu (research gap) lớn nhất tại VQG Bù Gia Mập trước năm 2019 là hầu hết các công trình nghiên cứu hệ sinh thái thực vật chỉ dừng lại ở phương pháp điều tra mô tả định tính, khảo sát thành phần loài theo tuyến (transect) rời rạc (Viện Sinh học Nhiệt đới, 1997; Phân viện Điều tra Quy hoạch Rừng Nam Bộ, 2004; VQG Bù Gia Mập, 2012). Chưa có bất kỳ công trình nào thiết lập hệ thống ô tiêu chuẩn định vị dài hạn quy chuẩn quốc tế (1 ha) để định lượng hóa cấu trúc lâm phần, đánh giá mối quan hệ tương hỗ giữa tính đa dạng sinh thái của thực vật thân gỗ với quy luật phân bố không gian và động thái tầng tán.

Luận án tiến sĩ lâm sinh của nghiên cứu sinh Vương Đức Hòa với đề tài "Nghiên cứu tính đa dạng thực vật thân gỗ và đặc điểm cấu trúc của một số kiểu rừng tại Vườn Quốc gia Bù Gia Mập, tỉnh Bình Phước" (chuyên ngành Lâm sinh, mã số 9620205, Viện Khoa học Lâm nghiệp Việt Nam) đã giải quyết triệt để khoảng trống khoa học này. Luận án đặt ra hai câu hỏi nghiên cứu cốt lõi:

  1. Mức độ đa dạng alpha ($\alpha$), beta ($\beta$) và vai trò sinh thái của các loài thực vật thân gỗ ($D_{1,3} \ge 10\text{ cm}$) biến đổi như thế nào giữa hai kiểu rừng chính là kiểu rừng kín thường xanh mưa nhiệt đới (Rkx) và kiểu rừng kín nửa thường xanh ẩm nhiệt đới (Rkn) qua các trạng thái phục hồi $IIIA_1, IIIA_2, IIIA_3$?
  2. Quy luật cấu trúc hình thái không gian (phân bố số cây theo đường kính $N/D$, phân bố theo chiều cao $N/H_{vn}$, tương quan $H/D$) và cấu trúc sinh thái (tổ thành loài $IV%$, độ ưu thế họ $FIV$) tuân theo các mô hình hàm toán học nào, và mối quan hệ định lượng giữa đa dạng loài với kích thước cây rừng diễn biến ra sao?

Hai giả thuyết nghiên cứu được kiểm định:

  • Giả thuyết 1 ($H_1$): Các loài thực vật thân gỗ ở các trạng thái rừng $IIIA_1, IIIA_2, IIIA_3$ có vai trò sinh thái và mức độ đóng góp khác nhau rõ rệt vào việc hình thành cấu trúc quần xã thực vật (QXTV) giữa hai kiểu rừng Rkx và Rkn.
  • Giả thuyết 2 ($H_2$): Cấu trúc tổ thành và quy luật phân bố không gian của thực vật thân gỗ bị chi phối trực tiếp bởi sự biến động về mật độ ($N$), tiết diện ngang ($G$) và trữ lượng ($V$) của các nhóm loài ưu thế sinh thái.

Nghiên cứu dựa trên khung lý thuyết sinh thái học phân tầng của Richards (1952), học thuyết phân loại thảm thực vật phát sinh quần thể của Thái Văn Trừng (1978, 1999), kết hợp lý thuyết ngẫu lựu và cạnh tranh ổ sinh thái (Niche Partitioning Theory). Khảo sát thực địa thu thập số liệu trên 6 ô tiêu chuẩn định vị (OTCDV) quy mô 1 ha ($100\text{ m} \times 100\text{ m}$) trong khung thời gian 2015–2016 tại vùng lõi VQG Bù Gia Mập, mang lại giá trị khoa học đột phá trong việc xây dựng cơ sở dữ liệu định lượng phục vụ giám sát đa dạng sinh học lâu dài.

Literature Review và Positioning

Lịch sử nghiên cứu thảm thực vật nhiệt đới chứng kiến sự đối thoại và tranh luận học thuật sâu sắc giữa các trường phái phân loại sinh thái. Trường phái sinh thái học tổ chức quần xã (Organismic / Climax Theory) khởi xướng bởi Clements (1916) và tiếp tục bởi Schimper (1918) xem quần xã thực vật là một "siêu sinh vật" vận động tiến tới trạng thái đỉnh cực ổn định dưới sự kiểm soát của khí hậu. Ngược lại, trường phái cá thể luận và liên tục sinh thái (Individualistic / Continuum Concept) của Ramensky (1938), Gleason (1926) và Whittaker (1953, 1975) lập luận rằng lớp phủ thực vật là sự tập hợp liên tục của các cá thể loài phản ứng độc lập với gradient môi trường sống.

Tại châu Âu, hai hệ thống phân loại song song định hình ngành địa thực vật: phân loại dựa trên loài chỉ thị của Braun-Blanquet (1928) thuộc trường phái Pháp – Thụy Sĩ và phân loại cấu trúc ngoại mạo (Physiognomic) của UNESCO (1973). Tại Việt Nam, hệ thống trạng thái rừng kinh doanh dựa trên trữ lượng của Loeschau (1960) được quy chuẩn hóa trong Quy phạm QPN 6-84, tuy nhiên bộc lộ nhược điểm là thiếu nền tảng sinh thái học phát sinh. Để khắc phục, Thái Văn Trừng (1978, 1999) đã xây dựng hệ sinh thái thực vật phát sinh quần thể phân chia thảm thực vật Việt Nam thành 14 kiểu rừng, khẳng định tính đa dạng cao và tính đặc trưng của các hệ sinh thái rừng nhiệt đới.

Trên bình diện phương pháp luận định lượng, việc đo lường đa dạng và cấu trúc rừng đã có những bước tiến vượt bậc từ các nghiên cứu quốc tế. Curtis & McIntosh (1951) xây dựng chỉ số giá trị quan trọng ($IV%$), Shannon & Wiener (1963) và Simpson (1949) chuẩn hóa thước đo đa dạng alpha, trong khi Condit (1992) và Dallmeier (1992) thiết lập tiêu chuẩn điều tra ô định vị 1 ha chia lưới thứ cấp để giám sát rừng nhiệt đới toàn cầu. Các nghiên cứu quốc tế tiêu biểu như Wending Huang và cộng sự (2003) tại dãy núi East Usambara (Tanzania) đã chứng minh cấu trúc lâm phần và mức độ tác động của con người chi phối trực tiếp đến độ phong phú loài; nghiên cứu của Slik và cộng sự (2003) tại vùng đất thấp Borneo chỉ ra sự tương tác giữa địa hình và lượng mưa tới cấu trúc họ Dầu (Dipterocarpaceae); hay nghiên cứu của Hossain và cộng sự (2007) tại VQG Dulhazara Safari (Bangladesh) định lượng tính đa dạng thực vật thông qua chỉ số Margalef ($d$) và Shannon-Wiener ($H'$).

Công trình của nghiên cứu sinh Vương Đức Hòa định vị chính xác tại giao điểm giữa lý thuyết thảm thực vật phát sinh quần thể Thái Văn Trừng và phương pháp luận định lượng hiện đại bằng ô tiêu chuẩn diện tích lớn (1 ha). Luận án kế thừa các nghiên cứu cấu trúc lâm phần của Vũ Tiến Hinh (2012), Đồng Sỹ Hiền (1974), Trần Văn Con (1991), Bảo Huy (1993), đồng thời vượt qua hạn chế của các khảo sát thực vật định tính trước đó tại VQG Bù Gia Mập, đưa cơ sở dữ liệu lâm sinh vùng Đông Nam Bộ tiệm cận với chuẩn mực sinh thái học thực địa quốc tế.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và kiểm chứng các lý thuyết nền tảng trong lâm học và sinh thái học quần xã:

  • Kiểm chứng học thuyết cấu trúc và phân tầng rừng nhiệt đới của Richards (1952): Luận án chứng minh rằng trong rừng tự nhiên nhiệt đới hỗn loài không đều tuổi, cấu trúc tầng tán không tồn tại ranh giới vật lý cố định mà là sự phân tầng liên tục, tích tụ tán lá phức hợp. Cấu trúc này được lượng hóa chính xác qua mối tương quan phi tuyến giữa đường kính ngang ngực ($D_{1,3}$) và chiều cao vút ngọn ($H_{vn}$).
  • Bổ sung cơ sở lý thuyết tổ thành lâm phần của Curtis & McIntosh (1951) và Mori et al. (1983): Thiết lập mối quan hệ giữa chỉ số giá trị quan trọng loài ($IV%$) và chỉ số giá trị gia đình ($FIV$). Luận án chứng minh sự chiếm lĩnh ưu thế tuyệt đối của các họ thực vật đặc trưng vùng nhiệt đới gió mùa (như Dipterocarpaceae, Fabaceae, Euphorbiaceae, Moraceae) quyết định tính ổn định của toàn bộ cấu trúc lâm phần.
  • Lượng hóa sự thay đổi paradigm từ mô tả định tính sang phân tích hàm phân bố xác suất sinh thái: Luận án mở rộng việc ứng dụng hàm mật độ xác suất Weibull 3 tham số, hàm khoảng cách và hàm giảm của Meyer để mô tả cấu trúc phân bố số cây theo cỡ kính ($N/D$) và chiều cao ($N/H$). Minh chứng dữ liệu thực nghiệm cho thấy phân bố số cây theo cấp đường kính của cả hai kiểu rừng Rkx và Rkn đều có dạng đường cong giảm (dạng chữ J ngược), khẳng định tính chất tái sinh tự nhiên liên tục và khả năng tự phục hồi của lâm phần rừng kín nhiệt đới.

Khung phân tích độc đáo

Khung phân tích của nghiên cứu là sự tích hợp đa tầng giữa 3 trụ cột lý thuyết:

  1. Hệ thống phân loại thảm thực vật phát sinh sinh thái (Thái Văn Trừng, 1999);
  2. Lý thuyết đa dạng sinh học quần xã và chỉ số sinh thái học (Shannon & Wiener, 1963; Simpson, 1949; Pielou, 1966; Whittaker, 1972);
  3. Toán học thống kê sinh thái và phân tích đa biến (Primer 6, phân tích thành phần chính PCA, phân tích đa chiều không gian MDS, đường cong tích lũy loài $K$-Dominance).

Khung phân tích xác định các điều kiện biên (boundary conditions) chặt chẽ: đối tượng khảo sát được giới hạn ở tầng cây gỗ có $D_{1,3} \ge 10\text{ cm}$ và chiều cao từ $3\text{ m}$ trở lên theo phân loại dạng sống cây có chồi trên mặt đất (Raunkiær, 1934; Huxley, 1992); không gian nghiên cứu định vị tại các trạng thái rừng $IIIA_1, IIIA_2, IIIA_3$ ở phân khu bảo vệ nghiêm ngặt và phục hồi sinh thái của VQG Bù Gia Mập; loại trừ các trạng thái rừng chuyển đổi bất định như rừng tre nứa thuần loại hoặc rừng hỗn giao tre nứa ($IIB, IIIA_1+L, IIIA_2+L$).

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo triết lý thực chứng khoa học (positivism), áp dụng phương pháp định lượng sinh thái học thực địa kết hợp điều tra xã hội học lâm nghiệp. Luận án xây dựng thiết kế nghiên cứu đa cấp (multi-level sampling design) trên quy mô diện tích lớn:

  • Địa bàn nghiên cứu: Phân khu bảo vệ nghiêm ngặt và phân khu phục hồi sinh thái VQG Bù Gia Mập, với tọa độ địa lý $12^\circ 08'30''$ đến $12^\circ 17'03''$ vĩ độ Bắc và $107^\circ 03'30''$ đến $107^\circ 14'30''$ kinh độ Đông; độ cao dưới $1.000\text{ m}$ so với mực nước biển.
  • Hệ thống ô mẫu: Thiết lập 6 ô tiêu chuẩn định vị (OTCDV) điển hình, mỗi ô có diện tích chuẩn $1\text{ ha}$ ($100\text{ m} \times 100\text{ m} = 10.000\text{ m}^2$). Trong đó:
    • 3 ô đại diện cho kiểu rừng kín thường xanh mưa nhiệt đới (Rkx) tại các trạng thái: $QXTV\text{ }Rkx-IIIA_3-OII$, $QXTV\text{ }Rkx-IIIA_2-OIII$, $QXTV\text{ }Rkx-IIIA_2-OIV$.
    • 3 ô đại diện cho kiểu rừng kín nửa thường xanh ẩm nhiệt đới (Rkn) tại các trạng thái: $QXTV\text{ }Rkn-IIIA_3-OI$, $QXTV\text{ }Rkn-IIIA_1-OV$, $QXTV\text{ }Rkn-IIIA_1-OVI$.
  • Cấu trúc ô thứ cấp: Mỗi OTCDV 1 ha được chia thành 25 ô thứ cấp kích thước $20\text{ m} \times 20\text{ m}$ ($400\text{ m}^2$) để đo đếm toàn diện cây gỗ trưởng thành ($D_{1,3} \ge 10\text{ cm}$), và bên trong mỗi ô thứ cấp thiết lập các ô tiểu tiêu chuẩn $5\text{ m} \times 5\text{ m}$ ($25\text{ m}^2$) để theo dõi cây tái sinh và thảm tươi.

Quy trình nghiên cứu rigorous

Quy trình thu thập dữ liệu tuân thủ các giao thức chuẩn mực sinh thái rừng nhiệt đới:

  • Giao thức đo đếm ngoại nghiệp: Tiến hành gắn số thẻ nhôm định vị cố định cho từng cá thể cây gỗ $D_{1,3} \ge 10\text{ cm}$, xác định tọa độ không gian $(x, y)$, đo đường kính ngang ngực ($D_{1,3}$) bằng thước kẹp kính/thước dây chuyên dụng với độ chính xác $0,1\text{ cm}$, đo chiều cao vút ngọn ($H_{vn}$) và chiều cao dưới cành ($H_{dc}$) bằng thước đo cao Blume-Leiss.
  • Định danh thực vật: Mẫu thực vật lá, hoa, quả được thu thập nghiêm ngặt theo quy chuẩn bảo tàng thực vật, phối hợp định danh hình thái học với các chuyên gia Viện Sinh thái học Miền Nam và Trường Đại học Nông Lâm TP. Hồ Chí Minh, đối soát danh lục với bộ sách Cây cỏ Việt Nam (Phạm Hoàng Hộ, 1999–2000), Danh lục các loài thực vật Việt Nam (Nguyễn Tiến Bân, 2001–2005) và hệ thống phân loại APG.
  • Đánh giá tình trạng bảo tồn: Tra cứu và phân hạng mức độ nguy cấp dựa trên Sách Đỏ Việt Nam (2007), Danh lục Đỏ IUCN (2017), Nghị định 32/2006/NĐ-CP của Chính phủ và Công ước CITES.
  • Điều tra kinh tế - xã hội vùng đệm: Phỏng vấn bán cấu trúc bảng hỏi tại 3 xã vùng đệm (xã Đắk Ơ, Bù Gia Mập thuộc huyện Bù Gia Mập, tỉnh Bình Phước và xã Quảng Trực thuộc huyện Tuy Đức, tỉnh Đắk Nông) với tổng dân số 5.845 hộ (26.296 nhân khẩu, tỷ lệ dân tộc thiểu số chiếm 56,21%).

Data và phân tích

Số liệu lâm phần được xử lý và phân tích bằng các phần mềm chuyên ngành hiện đại:

  • Chỉ số sinh thái học: Tính toán chỉ số quan trọng loài $IV% = (F% + N% + G%)/3$; chỉ số giá trị gia đình $FIV = RD_{fam}% + RDo_{fam}% + RDi_{fam}%$; chỉ số đa dạng loài Shannon-Wiener ($H' = -\sum p_i \ln p_i$), chỉ số phong phú Margalef ($d = (S-1)/\ln N$), chỉ số đồng đều Pielou ($J' = H'/\ln S$), chỉ số ưu thế Simpson ($Cd = \sum p_i^2$), chỉ số tương đồng sinh thái Sorensen ($SI = 2c/(a+b)$), và chỉ số hiếm $IR = (\sum C_i / S) \times 100$.
  • Mô hình hóa cấu trúc hình thái: Sử dụng các hàm phân bố xác suất phi tuyến trên phần mềm Statistica và phần mềm lâm nghiệp chuyên dụng để nắn các đường cong phân bố thực nghiệm:
    • Hàm phân bố Weibull 3 tham số: $f(x) = \frac{c}{b}\left(\frac{x-a}{b}\right)^{c-1}\exp\left[-\left(\frac{x-a}{b}\right)^c\right]$
    • Hàm Meyer: $N = N_0 \cdot \exp(-b \cdot D)$
    • Phương trình hồi quy phi tuyến tương quan $H/D$: dạng hàm parabol $H = a \cdot D^2 + b \cdot D + c$, hàm logarit $H = a + b \cdot \ln(D)$, hàm lũy thừa $H = a \cdot D^b$.
    • Kiểm định độ phù hợp giữa phân bố thực nghiệm và phân bố lý thuyết bằng tiêu chuẩn $\chi^2$ (Chi-square test) hoặc Kolmogorov-Smirnov ở mức ý nghĩa $\alpha = 0,05$.
  • Phân tích cấu trúc quần xã đa biến: Sử dụng gói phần mềm sinh thái PRIMER 6 (Plymouth Routines In Multivariate Ecological Research) để thiết lập ma trận tương đồng Bray-Curtis, vẽ đồ thị phân nhóm Cluster Analysis, đồ thị sắp xếp đa chiều phi chỉ số (Non-metric Multidimensional Scaling - MDS), phân tích thành phần chính (PCA) và đường cong tích lũy ưu thế $K$-Dominance. Quản lý cơ sở dữ liệu không gian rừng bằng MapInfo và Google Earth Pro.

Phát hiện đột phá và implications

Những phát hiện then chốt

Kết quả nghiên cứu đã làm sáng tỏ cấu trúc toàn diện của hệ sinh thái thực vật thân gỗ VQG Bù Gia Mập với những phát hiện mang tính quy luật sinh thái sâu sắc:

+-----------------------------------------------------------------------------------------+
|                        HỆ THỰC VẬT THÂN GỖ VQG BÙ GIA MẬP                               |
|       (Diện tích tự nhiên: 25.926 ha | Vùng đệm: 18.038 ha, 14.783 người DTTS)          |
+--------------------------------------------+--------------------------------------------+
|  KIỂU RỪNG KÍN THƯỜNG XANH MƯA NHIỆT ĐỚI   | KIỂU RỪNG KÍN NỬA THƯỜNG XANH ẨM NHIỆT ĐỚI |
|                   (Rkx)                    |                   (Rkn)                    |
+--------------------------------------------+--------------------------------------------+
| • Mật độ: 387 - 527 cây/ha                 | • Mật độ: 310 - 451 cây/ha                 |
| • Tiết diện ngang: 27,81 - 39,45 m2/ha     | • Tiết diện ngang: 21,34 - 34,21 m2/ha     |
| • Trữ lượng gỗ: 245,6 - 412,3 m3/ha        | • Trữ lượng gỗ: 189,2 - 356,8 m3/ha        |
| • Chỉ số Shannon H': 3,82 - 4,45           | • Chỉ số Shannon H': 3,41 - 3,98           |
| • Họ ưu thế FIV: Fabaceae, Dipterocarpaceae| • Họ ưu thế FIV: Lythraceae, Datiscaceae   |
| • Quy luật N/D: Tuân theo hàm Weibull (J)  | • Quy luật N/D: Tuân theo hàm Weibull (J)  |
| • Tương quan H/D: Hàm Logarit/Parabol      | • Tương quan H/D: Hàm Logarit/Parabol      |
+--------------------------------------------+--------------------------------------------+

  1. Sự phân hóa đa dạng loài và ưu thế họ sinh thái giữa hai kiểu rừng: Trong kiểu rừng kín thường xanh mưa nhiệt đới (Rkx), tính đa dạng loài ở mức rất cao với chỉ số Shannon-Wiener $H'$ dao động từ $3,82$ đến $4,45$; chỉ số phong phú Margalef $d$ đạt từ $10,21$ đến $14,65$; chỉ số đồng đều $J'$ từ $0,82$ đến $0,91$. Các họ có chỉ số giá trị gia đình ($FIV$) cao nhất thuộc về Fabaceae, Dipterocarpaceae, Euphorbiaceae, Burseraceae, Myrtaceae. Ngược lại, kiểu rừng kín nửa thường xanh ẩm nhiệt đới (Rkn) có chỉ số $H'$ dao động từ $3,41$ đến $3,98$; ưu thế thuộc về các họ rụng lá mùa khô đặc trưng như Lythraceae (với loài Bằng lăng ổi Lagerstroemia calyculata), Datiscaceae (với loài Đăng Tetrameles nudiflora), Combretaceae, Sterculiaceae.

  2. Quy luật cấu trúc mật độ, tiết diện ngang và trữ lượng lâm phần:

    • Tại kiểu rừng Rkx, mật độ cây gỗ dao động từ $387$ đến $527\text{ cây/ha}$, tổng tiết diện ngang ($G$) từ $27,81$ đến $39,45\text{ m}^2/\text{ha}$, và trữ lượng gỗ ($V$) đạt từ $245,6$ đến $412,3\text{ m}^3/\text{ha}$. Cấu trúc tổ thành loài cho thấy sự chi phối của các loài cây gỗ lớn như Vên vên (Anisoptera costata), Dầu rái (Dipterocarpus alatus), Sao đen (Hopea odorata), Gõ đỏ (Afzelia xylocarpa), Cẩm lai (Dalbergia oliveri).
    • Tại kiểu rừng Rkn, mật độ cây gỗ biến động từ $310$ đến $451\text{ cây/ha}$, tổng tiết diện ngang ($G$) từ $21,34$ đến $34,21\text{ m}^2/\text{ha}$, trữ lượng gỗ từ $189,2$ đến $356,8\text{ m}^3/\text{ha}$, với công thức tổ thành ghi nhận Bằng lăng (Lagerstroemia calyculata), Đăng (Tetrameles nudiflora), Râm (Canarium subulatum) chiếm tỷ lệ tổ thành sinh thái từ $35%$ đến $55%$.

  3. Mô hình hóa cấu trúc không gian đường kính ($N/D$) và chiều cao ($N/H_{vn}$): Kiểm định thống kê $\chi^2$ khẳng định phân bố số cây theo cỡ kính ($N/D_{1,3}$) của toàn bộ 6 OTCDV đều thích ứng cao nhất với hàm phân bố Weibull 3 tham số với $R^2 > 0,92$, sai số phương trình $SEE < 5%$, thể hiện cấu trúc giảm liên tục dạng chữ J ngược. Cây rừng tập trung mật độ cao nhất ở cấp kính nhỏ ($D_{1,3} = 10 - 25\text{ cm}$ chiếm trên $65 - 75%$ tổng số cá thể), chứng minh quá trình tái sinh tự nhiên dưới tán rừng diễn ra liên tục, tạo nguồn cây dự trữ bổ cập dồi dào cho tầng cây vượt tán. Phân bố số cây theo chiều cao ($N/H_{vn}$) thích ứng tối ưu với hàm phân bố khoảng cách và hàm Weibull với dạng một đỉnh lệch trái, tập trung chủ yếu ở tầng cao $12 - 20\text{ m}$.

  4. Tương quan sinh trưởng $H_{vn}/D_{1,3}$ và quan hệ đa dạng loài theo cấp kính: Mối quan hệ tương quan sinh trưởng giữa chiều cao vút ngọn và đường kính ngang ngực biểu diễn chính xác nhất qua phương trình hồi quy phi tuyến dạng hàm logarit ($H_{vn} = a + b \cdot \ln(D_{1,3})$) và hàm parabol với hệ số tương quan $R$ đạt từ $0,78$ đến $0,89$ ($p < 0,001$). Đồng thời, phân tích định lượng chỉ ra quy luật: tính đa dạng loài đạt giá trị cực đại ở các cấp đường kính nhỏ và trung bình ($10 - 30\text{ cm}$), sau đó giảm dần theo quy luật hàm mũ ở các cấp đường kính cổ thụ ($D_{1,3} > 60\text{ cm}$) do sự cạnh tranh đào thải tự nhiên về không gian dinh dưỡng và ánh sáng.

  5. Đánh giá nguồn tài nguyên thực vật quý hiếm và mức độ đe dọa: Nghiên cứu ghi nhận hệ thống thực vật nguy cấp, quý hiếm có giá trị bảo tồn toàn cầu cao tại VQG Bù Gia Mập với hàng chục loài nằm trong Sách Đỏ Việt Nam và Danh lục Đỏ IUCN như: Gõ đỏ (Afzelia xylocarpa - EN), Gõ mật (Sindora siamensis - EN), Cẩm lai (Dalbergia oliveri - EN), Giáng hương quả to (Pterocarpus macrocarpus - EN), Vên vên (Anisoptera costata - EN), Dầu cát (Dipterocarpus costatus - CR), Trắc bách diệp (Platycladus orientalis), Thủy tùng (Glyptostrobus pensilis). Chỉ số hiếm $IR$ tại các QXTV đạt từ $81,2%$ đến $94,6%$, phản ánh mức độ đe dọa sinh thái nghiêm trọng nếu sinh cảnh bị phân mảnh.

Implications đa chiều

  • Về mặt lý thuyết lâm học: Cung cấp mô hình toán học định lượng mô tả mối quan hệ giữa cấu trúc lâm phần và đa dạng sinh thái thực vật nhiệt đới vùng Đông Nam Bộ, làm sâu sắc thêm học thuyết diễn thế rừng thứ sinh và cân bằng động của hệ sinh thái rừng rụng lá/thường xanh.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình thiết lập ô tiêu chuẩn định vị 1 ha chia lưới thứ cấp kết hợp phân tích đa biến sinh thái (PRIMER 6, PCA, MDS) và hệ thống thông tin địa lý (GIS/Google Earth), tạo quy trình chuẩn có thể chuyển giao áp dụng cho các khu bảo tồn thiên nhiên và VQG trên toàn quốc.
  • Về mặt quản lý thực tiễn và kỹ thuật lâm sinh:
    • Cung cấp cơ sở khoa học để xác định các loài cây bản địa ưu thế, cây gỗ lớn có giá trị kinh tế - sinh thái cao (như các loài thuộc họ Dầu Dipterocarpaceae, họ Đậu Fabaceae) phục vụ công tác làm giàu rừng, phục hồi sinh thái các diện tích rừng nghèo kiệt ($IIIA_1$) và đất sau nương rẫy.
    • Hướng dẫn các giải pháp khoanh nuôi xúc tiến tái sinh tự nhiên có tác động lâm sinh (tỉa thưa điều chỉnh mật độ cây phi mục đích, mở tán cho cây tái sinh mục đích).
  • Về mặt chính sách bảo tồn: Đề xuất xây dựng quy chế phối hợp liên tỉnh giữa Ban Quản lý VQG Bù Gia Mập (Bình Phước) với Công ty Lâm nghiệp Nam Tây Nguyên, BQL Rừng phòng hộ Thác Mơ (Đắk Nông); triển khai chính sách chi trả dịch vụ môi trường rừng (DVMTR) gắn với giao khoán bảo vệ rừng cho $5.845$ hộ dân vùng đệm (đặc biệt là đồng bào dân tộc S'tiêng, M'nông) nhằm giảm áp lực khai thác lâm sản trái phép.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu xuất phát từ điều kiện thực địa và nguồn lực:

  1. Giới hạn về phạm vi không gian và số lượng ô mẫu: Mặc dù 6 OTCDV quy mô 1 ha cung cấp khối lượng dữ liệu mẫu lớn ($60.000\text{ m}^2$), song nghiên cứu mới chỉ tập trung tại vùng lõi (phân khu bảo vệ nghiêm ngặt và phục hồi sinh thái), chưa bao quát hết các sinh cảnh đất ngập nước ven sông Đắk Huýt hay các đai cao sát ranh giới tỉnh Đắk Nông.
  2. Giới hạn về đối tượng tầng cây: Luận án chỉ giới hạn phân tích định lượng thực vật thân gỗ có đường kính $D_{1,3} \ge 10\text{ cm}$, chưa lượng hóa chi tiết động thái biến động của tầng thảm tươi, dây leo gỗ và các loài thực vật phụ sinh (epiphytes) vốn rất phong phú trong rừng mưa nhiệt đới.
  3. Giới hạn về chuỗi thời gian quan trắc: Dữ liệu điều tra thu thập trong giai đoạn 2015–2016 mang tính chất mặt cắt tĩnh (cross-sectional snapshot), chưa đủ chuỗi dữ liệu dài hạn (5–10 năm) để đánh giá tốc độ tăng trưởng đường kính hàng năm ($MAI_D$), tỷ lệ chết tự nhiên ($mortality\text{ }rate$) và động thái tích lũy sinh khối/carbon sinh thái của từng loài cây theo thời gian.
  4. Giới hạn về tương tác vi khí hậu và thổ nhưỡng: Nghiên cứu chưa đo đạc đồng bộ các thông số lý - hóa tính chuyên sâu của đất (độ dày tầng đất, hàm lượng mùn, đạm, lân, kali) và tiểu khí hậu dưới tán rừng tương ứng với từng ô tiêu chuẩn.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) cần tập trung vào 4 định hướng:

  • Định hướng 1: Tiếp tục duy trì và tái đo đếm định kỳ hệ thống 6 OTCDV 1 ha theo chu kỳ 5 năm một lần để lượng hóa động thái tăng trưởng sinh khối ($AGB$) và khả năng cô lập carbon ($carbon\text{ }sequestration$) của hai kiểu rừng Rkx và Rkn trong bối cảnh biến đổi khí hậu.
  • Định hướng 2: Mở rộng nghiên cứu đa dạng di truyền (genetic diversity) và cấu trúc quần thể bằng chỉ thị phân tử DNA (microsatellite/ISSR) đối với các loài cây cực kỳ nguy cấp như Anisoptera costata, Dalbergia oliveri, Dipterocarpus costatus.
  • Định hướng 3: Ứng dụng công nghệ viễn thám độ phân giải cao LiDAR kết hợp ảnh máy bay không người lái (UAV) để mô hình hóa cấu trúc 3D tán rừng và phân loại thảm thực vật tự động trên toàn bộ diện tích 25.926 ha của Vườn Quốc gia.
  • Định hướng 4: Nghiên cứu thực nghiệm các mô hình làm giàu rừng và phục hồi sinh thái bằng các tổ hợp loài cây bản địa hỗn giao dựa trên công thức tổ thành sinh thái tối ưu đã được luận án xác lập.

Tác động và ảnh hưởng

  • Tác động học thuật (Academic Impact): Công trình là nguồn tài liệu tham khảo nền tảng cho các nhà nghiên cứu sinh thái, lâm học nhiệt đới; cung cấp bộ dữ liệu thực địa tin cậy cho các phân tích tổng hợp (meta-analysis) về cấu trúc rừng Đông Nam Á; ước tính tạo ra tiềm năng trích dẫn khoa học cao trong hệ thống tạp chí lâm nghiệp quốc gia và quốc tế.
  • Đổi mới quản lý ngành lâm nghiệp: Chuyển đổi căn bản phương pháp quản lý bảo tồn tại VQG Bù Gia Mập từ phương thức tuần tra bảo vệ thụ động sang giám sát đa dạng sinh học chủ động dựa trên hệ thống cơ sở dữ liệu số hóa địa không gian GIS.
  • Tác động chính sách cấp vùng và quốc gia: Cung cấp luận cứ khoa học thực chứng phục vụ việc điều chỉnh Quy hoạch bảo tồn đa dạng sinh học tỉnh Bình Phước đến năm 2030, hỗ trợ định giá dịch vụ môi trường rừng đầu nguồn cho lưu vực sông Bé và hệ thống thủy điện Thác Mơ.
  • Lợi ích kinh tế - xã hội: Góp phần định hình các mô hình phát triển sinh kế bền vững dưới tán rừng, giảm thiểu xung đột tài nguyên giữa cộng đồng dân tộc bản địa vùng đệm và Vườn Quốc gia thông qua chương trình đồng quản lý rừng.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học chuyên ngành Lâm sinh, Sinh thái học, Quản lý tài nguyên rừng: Tiếp cận phương pháp luận mẫu mực về thiết kế ô tiêu chuẩn định vị diện tích lớn (1 ha) và quy trình mô hình hóa toán học cấu trúc lâm phần.
  • Các nhà khoa học lâm nghiệp và chuyên gia sinh thái: Kế thừa hệ thống số liệu định lượng chi tiết về thành phần loài, công thức tổ thành, chỉ số đa dạng $H', d, J', FIV$ và quy luật phân bố Weibull của rừng nhiệt đới Đông Nam Bộ.
  • Ban Giám đốc và cán bộ kỹ thuật các Vườn Quốc gia, Khu Bảo tồn Thiên nhiên: Sở hữu bộ công cụ chẩn đoán sức khỏe rừng, quy trình giám sát động thái sinh thái dài hạn và cơ sở dữ liệu định vị để bảo vệ các loài cây nguy cấp, quý hiếm.
  • Các nhà hoạch định chính sách lâm nghiệp (Cục Lâm nghiệp, Sở Nông nghiệp và PTNT): Sử dụng các bằng chứng khoa học thực địa để xây dựng chính sách giao khoán quản lý rừng, chi trả dịch vụ môi trường và quy hoạch lâm nghiệp bền vững.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của nghiên cứu là gì và đã mở rộng lý thuyết nào?

Nghiên cứu đã mở rộng và lượng hóa thành công lý thuyết thảm thực vật sinh thái phát sinh quần thể của Thái Văn Trừng (1978, 1999) bằng cách tích hợp vào mô hình phân tích cấu trúc lâm phần định lượng hiện đại. Đóng góp độc đáo nhất là xác lập được mối tương quan toán học chặt chẽ giữa tính đa dạng loài theo các cấp đường kính ($D_{1,3}$) với hàm mật độ phân bố Weibull 3 tham số. Luận án chỉ ra rằng sự duy trì tính đa dạng loài cực đại ở các cấp kính nhỏ ($10 - 25\text{ cm}$) chính là cơ chế bảo toàn ổ sinh thái (niche partitioning) giúp các hệ sinh thái rừng kín nhiệt đới duy trì trạng thái ổn định động bền vững trước các biến động môi trường.

2. Phương pháp luận của luận án có điểm đổi mới đột phá nào so với các nghiên cứu trước đây tại Việt Nam?

So sánh với các nghiên cứu lâm sinh truyền thống trước đây (thường chỉ sử dụng các ô tiêu chuẩn tạm thời diện tích nhỏ $500 - 2.500\text{ m}^2$ điều tra đơn lẻ):

  • Luận án đã chuẩn hóa phương pháp thiết lập 6 ô tiêu chuẩn định vị dài hạn quy mô $1\text{ ha}$ ($100\text{ m} \times 100\text{ m}$) chia lưới 25 ô thứ cấp ($400\text{ m}^2$) theo tiêu chuẩn giám sát rừng nhiệt đới quốc tế (Condit, 1992; Dallmeier, 1992).
  • Ứng dụng đồng bộ phần mềm thống kê sinh thái đa biến PRIMER 6 (phân tích MDS, PCA, đường cong $K$-Dominance) kết hợp các phần mềm MapInfo và Google Earth để số hóa toàn diện tọa độ và đặc trưng sinh trưởng của từng cá thể cây gỗ, tạo lập nền tảng cơ sở dữ liệu không gian mở cho việc theo dõi diễn thế rừng lâu dài.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực địa là gì?

Phát hiện bất ngờ nhất nằm ở sự biến thiên đối nghịch về mật độ và chỉ số đa dạng loài theo cấp đường kính: Mặc dù các trạng thái rừng $IIIA_3$ có trữ lượng gỗ và tiết diện ngang tầng cây ưu thế áp đảo hoàn toàn ($V$ đạt tới $412,3\text{ m}^3/\text{ha}$), nhưng chỉ số đa dạng loài $H'$ ở tầng cây vượt tán ($D_{1,3} > 60\text{ cm}$) lại suy giảm mạnh và chỉ tập trung vào một số rất ít loài thuộc họ Dầu (Dipterocarpaceae) hoặc họ Bằng lăng (Lythraceae). Điều này phản ánh quy luật cạnh tranh loại trừ khốc liệt ở tầng tán trên cao của rừng nhiệt đới, bác bỏ quan niệm định tính trước đây cho rằng rừng càng già thì độ đa dạng loài ở mọi tầng kích thước đều tăng đều như nhau.

4. Luận án có cung cấp quy trình tái lập (replication protocol) đầy đủ không?

Có. Luận án cung cấp chi tiết toàn bộ quy trình tái lập bao gồm: tọa độ địa lý chính xác của 6 OTCDV; sơ đồ phân chia ô thứ cấp ($20\text{ m} \times 20\text{ m}$) và ô tiểu tiêu chuẩn ($5\text{ m} \times 5\text{ m}$); giao thức gắn biển nhôm định vị cá thể cây gỗ; bảng mã định danh loài; công thức toán học chi tiết tính toán các chỉ số $IV%, FIV, H', d, J', Cd, IR$; và các tham số khởi tạo của hàm phân bố Weibull, Meyer trên phần mềm Statistica và PRIMER 6.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Chương trình nghiên cứu 10 năm tiếp theo tập trung vào:

  • Tái đo đếm hệ thống 6 OTCDV 1 ha vào các năm 2021 và 2026 để xác lập mô hình sinh trưởng thực nghiệm, tỷ lệ tích lũy sinh khối và trữ lượng carbon rừng.
  • Mở rộng mạng lưới ô mẫu 1 ha lên 15 ô bao phủ toàn bộ các đai độ cao và kiểu phụ thổ nhưỡng của VQG Bù Gia Mập.
  • Nghiên cứu giải mã gen quần thể đối với các loài nguy cấp (Gõ đỏ, Cẩm lai, Vên vên) và tích hợp hệ thống cảm biến vi khí hậu tự động dưới tán rừng để đánh giá phản ứng của hệ sinh thái trước hiện tượng biến đổi khí hậu toàn cầu.

Kết luận

Công trình nghiên cứu của tiến sĩ Vương Đức Hòa đại diện cho bước tiến học thuật xuất sắc của ngành Lâm sinh Việt Nam với 5 đóng góp cốt lõi:

  1. Lượng hóa toàn diện tính đa dạng thực vật thân gỗ: Xác lập hệ thống chỉ số đa dạng sinh thái ($H', d, J', Cd, IR$) cho 2 kiểu rừng tiêu biểu (Rkx và Rkn) qua các trạng thái phục hồi $IIIA_1, IIIA_2, IIIA_3$, ghi nhận mức độ đa dạng cao của hệ thực vật vùng chuyển tiếp Nam Tây Nguyên – Đông Nam Bộ.
  2. Xác lập mô hình toán học cấu trúc lâm phần chuẩn xác: Chứng minh sự thích ứng ưu việt của hàm phân bố xác suất Weibull 3 tham số trong việc mô phỏng quy luật phân bố số cây theo đường kính ($N/D_{1,3}$) và hàm logarit cho quan hệ tương quan chiều cao – đường kính ($H_{vn}/D_{1,3}$).
  3. Định vị vai trò sinh thái của các nhóm loài ưu thế: Định lượng chính xác chỉ số quan trọng loài ($IV%$) và chỉ số giá trị gia đình ($FIV$), làm rõ cơ chế chiếm lĩnh ưu thế của các họ Dipterocarpaceae, Fabaceae, Lythraceae trong cấu trúc rừng nhiệt đới.
  4. Chuẩn hóa phương pháp luận ô định vị quy mô lớn (1 ha): Đưa quy trình điều tra định lượng dài hạn kết hợp phân tích đa biến sinh thái (PRIMER 6, PCA, MDS) và công nghệ GIS vào thực tiễn quản lý rừng đặc dụng tại Việt Nam.
  5. Đề xuất hệ thống giải pháp kỹ thuật lâm sinh và chính sách đồng bộ: Cung cấp cơ sở khoa học tin cậy cho công tác làm giàu rừng bản địa, bảo tồn nguồn gen các loài thực vật nguy cấp quý hiếm (Sách Đỏ Việt Nam, IUCN) và phát triển sinh kế bền vững cho cộng đồng vùng đệm VQG Bù Gia Mập.

Công trình tạo nền tảng lý luận và thực tiễn vững chắc, mở ra các hướng nghiên cứu dài hạn về động thái carbon, sinh thái học phục hồi và thích ứng biến đổi khí hậu của rừng nhiệt đới trong khu vực và quốc tế.