Tổng quan về luận án

Hệ sinh thái rừng trên núi đá vôi (karst) tại Vườn Quốc gia (VQG) Cát Bà là một trong những cảnh quan sinh thái biển đảo độc đáo bậc nhất Việt Nam và khu vực Đông Nam Á, mang giá trị bảo tồn đa dạng sinh học toàn cầu. Tuy nhiên, các cảnh quan rừng karst có tính nhạy cảm sinh thái cực kỳ cao do lớp phủ thổ nhưỡng mỏng, tốc độ xói mòn vật lý vượt xa tốc độ phong hóa hình thành đất và khả năng giữ nước kém. Bối cảnh thực tế cho thấy: "Tổng diện tích của VQG là 17.362,96 ha, trong đó diện tích quy hoạch cho đất lâm nghiệp là 10.443,6 ha, chiếm 60,2% tổng diện tích tự nhiên. Rừng tự nhiên (chủ yếu là rừng trên núi đá vôi) có 5.891,69 ha chiếm 56,6% diện tích đất lâm nghiệp; trong đó: 1.056,27 ha là rừng nguyên sinh chiếm 17,3% tổng diện tích rừng; rừng thứ sinh 4.856,09 ha, chủ yếu là rừng nghèo và nghèo kiệt, chiếm gần 70% tổng diện tích rừng tự nhiên (VQG Cát Bà, 2020)". Sự suy thoái nghiêm trọng do khai thác kiệt quệ và nương rẫy trong quá khứ đã làm đứt gãy cấu trúc quần xã, khiến nhiều loài thực vật quý hiếm như Nghiến (Excentrodendron tonkinense), Đinh (Markhamia stipulata), Kim giao (Nageia fleuryi), Cọ hạ long (Livistona halongensis) rơi vào tình trạng bị đe dọa tuyệt chủng cục bộ.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tiễn: Các công trình lâm học trước đây tại Việt Nam (Thái Văn Trừng, 1978; Trần Hữu Viên, 2005; Phan Thanh Lâm, 2016) chủ yếu dừng lại ở mức độ mô tả định tính thảm thực vật hoặc thống kê tài nguyên đơn thuần, chưa giải mã được động thái phục hồi rừng theo các đơn vị vi lập địa karst kết hợp sinh thái định lượng. Đồng thời, việc ứng dụng các biện pháp lâm sinh truyền thống thường tách rời động thái thảm thực vật khỏi bối cảnh tương tác sinh thái - xã hội (Social-Ecological Systems - SES), dẫn đến các can thiệp phục hồi thiếu tính thích ứng.

Để giải quyết khoảng trống này, luận án tập trung vào ba câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H):

  • RQ1: Cấu trúc tầng cây cao, động thái tái sinh tự nhiên và đặc tính lý - hóa học đất biến đổi như thế nào giữa các kiểu phụ rừng thứ sinh nhân tác so với trạng thái rừng nguyên sinh đối chứng trên núi đá vôi?
    • H1: Quá trình suy thoái nhân tác làm suy giảm phi tuyến tính độ phì đất và mật độ cây tái sinh có triển vọng (TSTV), làm thay đổi căn bản tổ thành loài tầng cây cao từ nhóm loài ưu thế đỉnh cực sang nhóm tiên phong ưa sáng.
  • RQ2: Các nhân tố sinh thái - xã hội (khoảng cách dân sinh, áp lực khai thác lâm sản ngoài gỗ, địa hình, thổ nhưỡng) tương tác và chi phối quá trình phục hồi rừng tự nhiên như thế nào?
    • H2: Cường độ tác động nhân sinh kết hợp với tỷ lệ đá lộ đầu và độ dày tầng đất là những biến số chính kiểm soát quỹ đạo diễn thế và mức độ tương đồng quần xã.
  • RQ3: Giải pháp kỹ thuật lâm sinh và khung quản trị nào tối ưu hóa khả năng tự phục hồi của từng kiểu phụ rừng karst?
    • H3: Phân vùng phục hồi dựa trên hệ thống tiêu chí định lượng 5 sao quốc tế (SER) cho phép thiết lập các biện pháp can thiệp có mục tiêu (khoanh nuôi xúc tiến tái sinh tự nhiên so với khoanh nuôi có trồng bổ sung loài bản địa), vượt trội hơn các can thiệp đồng nhất.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa Lý thuyết diễn thế sinh thái (Hobbs & Walker, 2007), Thuyết lập địa sinh thái rừng nhiệt đới (Thái Văn Trừng, 1978; Holl & Aide, 2011) và Khung phục hồi sinh thái đa chuẩn mực (McDonald et al., 2016). Luận án tạo bước đột phá khi phân loại chi tiết 7 kiểu phụ rừng và 21 quần xã thực vật rừng (QXTVR) điển hình, lượng hóa toàn diện các chỉ số đa dạng (Margalef, Menhinick, Simpson, Shannon-Wiener, Rényi) và cơ cấu lại hệ thống giải pháp kỹ thuật lâm sinh trên quy mô 63 ô tiêu chuẩn (OTC) định lượng, 16 OTC mô hình thử nghiệm và 102 đối tượng phỏng vấn xã hội học trong hai giai đoạn khảo sát (2017 - 2018 và 2021 - 2022).

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về hệ sinh thái karst phản ánh một trường nghiên cứu phức tạp với nhiều dòng lý thuyết đan xen:

[Tiếp cận Thổ nhưỡng - Địa mạo Karst]
(Ford & Williams, 2007; Gutiérrez et al., 2014)
[Thích nghi Sinh thái & Cấu trúc Quần xã]
(Denslow, 1987; Wang et al., 2019; Liu et al., 2021)
[Tranh luận Phục hồi: Thụ động vs Chủ động]    [Động thái Lâm sinh Karst Việt Nam]
(Lamb & Gilmour, 2003; Holl & Aide, 2011)      (Thái Văn Trừng, 1978; Trần Hữu Viên, 2005)
        [Định vị Luận án: Sinh thái Định lượng & SES]
        (Lê Hồng Liên, 2024 - VQG Cát Bà)
  1. Dòng nghiên cứu địa mạo - thủy văn và độ nhạy cảm thổ nhưỡng karst: Ford & Williams (2007) cùng Gutiérrez et al. (2014) đã chứng minh rằng các cảnh quan carbonate toàn cầu (chiếm 12% diện tích lục địa) có đặc tính thủy văn karst phân mảnh, lớp đất mặt mỏng, tốc độ xói mòn vật lý cao hơn nhiều so với tốc độ hình thành đất. Khi thảm phủ rừng bị phá hủy, sự xói mòn rửa trôi chất hữu cơ diễn ra triệt để, tạo ra hiện tượng "sa mạc hóa đá" (rocky desertification).
  2. Dòng nghiên cứu cấu trúc quần xã và thích nghi sinh thái học: Liu et al. (2021) và Wang et al. (2019) chỉ ra rằng thực vật núi đá vôi phát triển các cơ chế chịu hạn sinh lý, rễ bám sâu vào khe nứt đá và thích nghi với môi trường giàu ion Ca2+, Mg2+. Chen et al. (2022) khi khảo cứu rừng mưa nhiệt đới Tây Song Bản Nạp (Xishuangbanna, Trung Quốc) nhấn mạnh sự tồn tại của các quần hợp sinh thái đơn ưu xen lẫn tỷ lệ lớn các loài hiếm (rare species), nơi 1 - 2 loài ưu thế chiếm tới 90% số cá thể tầng cây cao (van der Velden et al., 2014).
  3. Dòng nghiên cứu kỹ thuật lâm sinh và phục hồi rừng: Tiếp cận của Lamb & Gilmour (2003) phân định rõ ba cấp độ: Cải tạo (Reclamation), Khôi phục (Restoration) và Phục hồi (Rehabilitation). Holl & Aide (2011) thiết lập khung căn cứ lựa chọn giữa phục hồi thụ động (dựa vào tái sinh tự nhiên) và phục hồi chủ động (can thiệp lâm sinh).

Trong y văn tồn tại hai cuộc tranh luận lý thuyết lớn:

  • Tranh luận 1: Phục hồi thụ động (Passive) đối đầu Phục hồi chủ động (Active). Các nhà bảo tồn truyền thống ưu tiên khoanh nuôi thuần túy để duy trì tính tự nhiên nguyên vẹn. Ngược lại, các nhà lâm sinh ứng dụng cho rằng đối với lập địa karst thoái hóa nghiêm trọng bị kiệt nguồn hạt giống, phục hồi tự nhiên có thể mất hàng thế kỷ hoặc bị bế tắc ở giai đoạn trảng cây bụi thoái hóa, đòi hỏi phải trồng làm giàu hoặc gieo hạt thẳng (Hoàng Kim Ngũ & Bùi Thế Đồi, 2002; Trần Hữu Viên, 2004).
  • Tranh luận 2: Tính cân bằng sinh thái của quần hợp đơn ưu so với đa dạng loài cao. Một số quan điểm cho rằng rừng karst tự nhiên phải đạt độ đa dạng loài cực đại (Denslow, 1987), trong khi các bằng chứng thực nghiệm tại Xishuangbanna (Chen et al., 2022) và Cúc Phương (Nguyễn Bá Thụ, 1995) chứng minh rừng karst trưởng thành thường hình thành các ưu hợp cục bộ với chỉ số đa dạng loài tổng thể thấp hơn rừng núi đất nhưng tính đặc hữu và độ thích nghi lại cao hơn rõ rệt.

Luận án định vị chính xác tại điểm giao thoa giữa sinh thái học định lượng và lâm sinh học cảnh quan đảo. Bằng cách so sánh đối chứng với các nghiên cứu quốc tế tại Tây Song Bản Nạp (Chen et al., 2022) và chương trình phục hồi karst quy mô lớn tại Tây Nam Trung Quốc (Tong et al., 2018; Ouyang et al., 2019), công trình tại VQG Cát Bà đóng góp một mô hình phân tích toàn diện: không chỉ xem xét lâm phần như một thực thể sinh học độc lập mà tích hợp chặt chẽ các biến số nhân sinh, thể chế quản lý và lịch sử xáo trộn đất đai vào cấu trúc mô hình phân cấp cụm và phân tích thành phần chính.

Đó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 Lý thuyết diễn thế sinh thái (Ecological Succession Theory - Clements, 1916; Hobbs & Walker, 2007) và Thuyết phân loại thảm thực vật phát sinh sinh thái của Thái Văn Trừng (1978) trong điều kiện lập địa karst biển đảo đặc thù thông qua các luận điểm khoa học có tính hệ thống:

                     [QUỸ ĐẠO DIỄN THẾ & PHỤC HỒI RỪNG KARST]
                     
    - Độ phì đất hồi phục            - Điều chỉnh vi khí hậu             - OM% < 2.5%
    - IVI loài mục đích tăng         - Trồng bổ sung loài bản địa        - TSTV thiếu hụt
    - Tầng tán phân hóa              - Bảo vệ nghiêm ngặt                - Tỷ lệ đá lộ đầu cao
  • Mở rộng lý thuyết diễn thế thoái hóa và phục hồi: Xác lập mô hình động thái biến đổi cấu trúc giữa 7 kiểu phụ rừng (từ kiểu phụ thổ nhưỡng kiệt nước ít bị tác động I.Đk1 đến kiểu phụ thứ sinh sau nương rẫy I.Np2-2). Luận án chứng minh rằng trên nền đá vôi, diễn thế phục hồi không diễn ra theo đường tuyến tính đơn giản mà phụ thuộc vào các ngưỡng suy thoái sinh thái (ecological thresholds) của tầng đất mặt. Khi hàm lượng mùn (OM%) và dung tích hấp thu cation (CEC) suy giảm dưới ngưỡng tới hạn, quá trình diễn thế tự nhiên bị phong tỏa (arrested succession), đòi hỏi các tác động hỗ trợ lâm sinh có chủ đích.
  • Mô hình hóa lý thuyết qua các mệnh đề khoa học (Propositions):
    • Mệnh đề P1: Mức độ tương đồng thành phần loài (SI) giữa tầng cây cao và lớp cây tái sinh tỷ lệ nghịch với cường độ xáo trộn nhân tác trong quá khứ.
    • Mệnh đề P2: Trong các kiểu phụ rừng thoái hóa nặng, sự hiện diện của tầng cây bụi thảm tươi dày đặc đóng vai trò như một "bộ lọc sinh thái" cản trở sự nảy mầm của các loài cây gỗ bản địa tái sinh từ hạt.
    • Mệnh đề P3: Cấu trúc không gian và tổ thành loài tầng cây cao tại các vi lập địa karst chịu sự kiểm soát đồng thời của độ dốc/tỷ lệ đá lộ đầu (nhân tố vô sinh) và khoảng cách tới ranh giới khu dân cư/đường giao thông (nhân tố hữu sinh nhân tác).

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

Luận án thiết lập một khung phân tích liên ngành tích hợp 3 trụ cột lý thuyết: (1) Sinh thái học thảm thực vật định lượng, (2) Thổ nhưỡng học rừng karst, và (3) Hệ thống Sinh thái - Xã hội (SES).

Khung phân tích này kết hợp chặt chẽ việc phân loại thảm thực vật theo tiếp cận phát sinh của Thái Văn Trừng (1978), phân cấp trạng thái rừng theo Thông tư 33/2018/TT-BNNPTNT và Loeschau (1963) với các phương pháp toán sinh thống kê đa biến tiên tiến.

Điều kiện biên (Boundary conditions): Khung phân tích được xác lập cho hệ sinh thái rừng kín lá rộng thường xanh mưa ẩm nhiệt đới phát triển trên núi đá vôi thuộc đai độ cao dưới 500 m so với mực nước biển, trong vùng sinh thái biển đảo Đông Bắc Việt Nam, chịu tác động của khí hậu nhiệt đới gió mùa có mùa đông lạnh khô và lượng mưa trung bình năm khoảng 1.700 - 1.800 mm.

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

Thiết kế nghiên cứu

Luận án vận hành trên nền tảng triết học Hiện thực Phê phán (Critical Realism) kết hợp Thực chứng luận Định lượng (Empirical Positivism). Nghiên cứu thừa nhận các quy luật cấu trúc sinh thái tồn tại khách quan nhưng biểu hiện hiện tại của chúng là kết quả tích hợp lịch sử của các can thiệp xã hội. Do đó, thiết kế nghiên cứu đa tầng bậc (Multi-level mixed design) được áp dụng:

[Khảo sát Tuyến: 13 Tuyến (30m x Chiều dài tuyến, đệm 10m)]
[Hệ thống Ô tiêu chuẩn Whittaker Cải tiến: 63 OTC (1.000 m²)]
             [Phân tích Lý - Hóa Đất]              [Khảo sát Xã hội học Lâm nghiệp]
             - 21 mẫu hỗn hợp (0-30 cm)            - 102 phiếu phỏng vấn bán cấu trúc
             - 10 chỉ tiêu lý hóa tiêu chuẩn       - 4 nhóm đối tượng liên quan

Cụ thể, quy mô mẫu nghiên cứu bao gồm:

  • Tầng vĩ mô (Tuyến điều tra): 13 tuyến khảo sát thực địa chính (bề rộng 30 m, mở rộng quan sát mỗi bên 10 m, với các tuyến phụ 2 - 12 km phân nhánh mỗi 1,5 km), cắt qua mọi dạng địa hình từ chân núi, sườn dốc đến đỉnh núi cao nhất.
  • Tầng trung mô (Ô tiêu chuẩn thảm thực vật): 63 OTC điển hình tạm thời (kích thước 1.000 m², 50 m x 20 m) phân bổ đều cho 7 kiểu phụ rừng (mỗi kiểu phụ gồm 3 QXTVR x 3 OTC lặp lại). Ngoài ra, thiết lập 16 OTC (1.000 m²) để đánh giá 4 mô hình kỹ thuật lâm sinh thử nghiệm (4 biện pháp x 4 OTC/biện pháp).
  • Tầng vi mô (Ô thứ cấp lồng ghép theo Whittaker): Trong mỗi OTC 1.000 m² (Ô A), lập 01 ô thứ cấp B (100 m² tại tâm) để đo đếm cây gỗ nhỏ ($D_{1.3} < 6\text{ cm}$, $H_{vn} > 2\text{ m}$) và 20 ô dạng bản C (4 m², tổng $80\text{ m}^2$) phân bố đều để đo đếm cây tái sinh ($D_{1.3} < 6\text{ cm}$ và/hoặc $H_{vn} \le 2\text{ m}$).
  • Tầng xã hội học: 102 đối tượng phỏng vấn bán cấu trúc (2 cán bộ quản lý VQG, 10 cán bộ kiểm lâm từ 5 trạm bảo vệ rừng, 30 đại diện cộng đồng từ 6 xã, và 60 hộ gia đình nhận khoán bảo vệ rừng tại 6 xã vùng lõi và vùng đệm).
  • Tầng thổ nhưỡng: 21 mẫu đất hỗn hợp tầng mặt (0 - 30 cm) thu thập tại 5 điểm (4 góc và tâm ô) của 21 QXTVR điển hình.

Quy trình nghiên cứu rigorous

Quy trình thu thập và chuẩn hóa dữ liệu tuân thủ các tiêu chuẩn kỹ thuật nghiêm ngặt:

Hạng mục nghiên cứu Quy chuẩn / Dụng cụ / Phương pháp áp dụng Tiêu chuẩn kiểm soát chất lượng
Đo trắc lượng cây gỗ Thước dây chuyên dụng (đo D1.3, mm), Máy đo cao laser Vertex (đo Hvn, độ chính xác 0,5 m) Kiểm định sai số thiết bị trước mỗi đợt khảo sát
Xác định độ tàn che Ảnh chụp tán thẳng đứng tại 5 điểm/OTC (tâm và 4 hướng Đông, Tây, Nam, Bắc cách 5 m) Chuẩn hóa ảnh quang học số hóa
Phân loại dạng sống & Công dụng Phổ dạng sống Raunkiaer (1934) cải tiến cho Việt Nam (Nguyễn Nghĩa Thìn, 2007: 14 dạng sống); Phân loại giá trị sử dụng 11 nhóm Đối chiếu Sách Đỏ Việt Nam (2007), Nghị định 06/2019/NĐ-CP, Nghị định 84/2021/NĐ-CP, CITES, IUCN (2020)
Phân tích lý tính đất Tỷ trọng ($d$) bình Picnome; Dung trọng ($D$) ống dung tích tâm OTC (TCVN 6860:2001); Cấp hạt 3 cấp FAO (TCVN 8567:2010) Thực hiện tại Phòng Thí nghiệm Thổ nhưỡng Đại học Lâm nghiệp
Phân tích hóa tính đất pHH2O máy pH-meter (TCVN 5979:2007); Mùn OM% Chiurin (TCVN 4050-85); Đạm tổng số $N%$ Kjeldahl (TCVN 6498:1999); Lân tổng số $P_2O_5%$ so màu (TCN 453:2001); Kali tổng số $K_2O%$ quang kế ngọn lửa (TCN 454:2001); $Ca^{2+}$, $Mg^{2+}$ Complexon; CEC Amoni-Axetat (TCVN 6646:2000) Chuẩn mẫu đối chứng nội bộ và lặp lại 3 lần cho mỗi chỉ tiêu

Tam giác đạc (Triangulation) được thiết lập giữa ba nguồn dữ liệu: số liệu sinh thái cấu trúc thực địa, số liệu phân tích lý - hóa thổ nhưỡng trong phòng thí nghiệm, và dữ liệu điều tra xã hội học lâm nghiệp, đảm bảo tính hợp lệ kiến tạo (construct validity) và độ tin cậy ngoại tại (external validity).

Data và phân tích

Toàn bộ dữ liệu định lượng được số hóa và xử lý trên môi trường R Studio (R phiên bản 4.0):

  • Phân tích Thành phần Chính không định hướng (PCA): Sử dụng R package vegan để phân tích ma trận tương quan đa biến tuyến tính giữa Ma trận 1 (7 kiểu phụ rừng và 21 QXTVR cho cả tầng cây cao và lớp cây tái sinh) và Ma trận 2 gồm 16 biến môi trường - lâm học: Đường kính ($D_{1.3}$), Chiều cao ($H_{vn}$), Tiết diện ngang ($G$), Trữ lượng ($M$), Độ cao tuyệt đối, Độ dốc, Tỷ lệ đá lẫn, Tỷ lệ đá lộ đầu, Mùn (OM%), Độ xốp, Độ tàn che, Độ che phủ thảm tươi, Tỷ lệ cây mục đích, Mật độ cây tái sinh, Số lượng cây TSTV, và Tần suất tác động cộng đồng. Mức độ giải thích phương sai được đánh giá qua các giá trị Eigenvalues và hệ số tương quan $r, r^2$ trên không gian hai chiều của Trục 1 và Trục 2.
  • Phân cấp cụm kết tụ (AHC - Agglomerative Hierarchical Clustering): Ứng dụng R package Factoextra dựa trên ma trận khoảng cách Euclid (Euclidean distance) và thuật toán liên kết phương sai tối thiểu Ward (Ward’s minimum variance method) để phân nhóm khách quan các kiểu phụ rừng và quần xã thực vật theo mức độ tương đồng sinh thái.
  • Hệ thống đánh giá phục hồi 5 sao SER: Vận dụng khung tiêu chuẩn phục hồi quốc tế (McDonald et al., 2016) trên biểu đồ Radar. Các chỉ số lâm học (mật độ, $D_{1.3}, G, M$, số loài, Margalef, Menhinick, Simpson, Shannon-Wiener của tầng cây cao và tầng tái sinh) cùng chỉ số tác động xã hội được chuẩn hóa thành thang điểm từ 1 đến 5 bằng hàm phân vị Quartile ($<25% = 1-2$ điểm; $25-50% = 3$ điểm; $50-75% = 4$ điểm; $>75% = 5$ điểm). Mức phục hồi đạt chuẩn sinh thái khi tất cả các chỉ số đạt tối thiểu 4 điểm (ngưỡng hài hòa tự nhiên).

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

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

Nghiên cứu mang lại 5 phát hiện khoa học mang tính đột phá với các cứ liệu thực chứng định lượng:

  1. Phân hóa sâu sắc giữa 7 kiểu phụ rừng và 21 QXTVR trên nền xáo trộn nhân sinh: Luận án đã làm rõ sự suy thoái nghiêm trọng của rừng tự nhiên karst tại Cát Bà khi trạng thái rừng nguyên sinh giàu (I.Đk1) chỉ còn chiếm 17,3% (1.056,27 ha), trong khi rừng thứ sinh nghèo và nghèo kiệt (I.Np1-2, I.Np1-3, I.Np1-4, I.Np2-1, I.Np2-2) chiếm tới gần 70% tổng diện tích rừng tự nhiên. Tầng cây cao ở các trạng thái suy thoái bị mất hoàn toàn các loài cây gỗ lớn ưu thế thuộc họ Fabaceae, Dipterocarpaceae, Meliaceae, thay thế bằng các loài tiên phong chịu hạn như Ô rô (Acanthus ebracteatus), Ba bét (Mallotus floribundus), Côm tầng, Chẹo tía.
  2. Sự đứt gãy động thái giữa tầng cây cao và lớp cây tái sinh: Tồn tại sự phân kỳ rõ rệt giữa mật độ hạt nảy mầm và mật độ cây tái sinh có triển vọng (TSTV). Mặc dù mật độ tái sinh tự nhiên ban đầu ở các kiểu phụ thoái hóa (I.Np1-3, I.Np1-4) có thể đạt từ 3.500 đến trên 5.000 cây/ha, tỷ lệ cây TSTV ($H_{vn} > 1\text{ m}$, thân thẳng, không sâu bệnh) chỉ đạt dưới 25 - 30%. Hiện tượng thiếu hụt tái sinh nghiêm trọng đối với các loài cây bản địa quý hiếm (Nghiến, Kim giao, Lát hoa) đã được ghi nhận trên khắp các lâm phần vùng đệm do thiếu hụt nguồn giống mẹ và điều kiện vi khí hậu khắc nghiệt.
  3. Quy luật suy thoái thổ nhưỡng karst theo cường độ xáo trộn: Phân tích 21 mẫu đất khẳng định tính chất đất biến động tiêu cực rõ rệt dọc theo gradient suy thoái:
    • Kiểu phụ nguyên sinh đối chứng I.Đk1 duy trì hàm lượng mùn tầng mặt cao (OM > 8,5%), đạm tổng số $N > 0,35%$, dung tích hấp thu $CEC > 35\text{ cmol}^+/\text{kg}$.
    • Kiểu phụ thoái hóa sau nương rẫy I.Np2-2 và chân núi I.Np2-1 bị suy giảm mùn nghiêm trọng (OM < 2,5%), đất bị nén chặt, dung trọng tăng cao ($D > 1,35\text{ g/cm}^3$), độ xốp giảm xuống dưới 42%, làm mất khả năng giữ ẩm mao dẫn trong mùa khô.
  4. Phát hiện nghịch lý trong mô hình phục hồi nhân tạo: Kết quả điều tra thực chứng trên các mô hình lâm sinh khẳng định một kết quả phản trực giác quan trọng: Can thiệp trồng bổ sung mật độ cao hoặc làm giàu rừng bằng cây sinh trưởng nhanh có thể đẩy nhanh độ tàn che tán rừng tạm thời nhưng lại làm suy giảm tính đa dạng loài thực vật bản địa tái sinh dưới tán. Cụ thể: "Số loài cây trong mỗi QXTV của các mô hình phục hồi biến động từ 13 - 20 loài, trong đó số loài ở mô hình khoanh nuôi bảo vệ dao động từ 16 - 20 loài, cao hơn hẳn ở mô hình khoanh nuôi xúc tiến tái sinh có trồng bổ sung, chỉ từ 13 - 15 loài (Triệu Thái Hưng, 2021)".
  5. Cơ chế tác động phi đối xứng của hệ sinh thái - xã hội: Kết quả phân tích PCA chứng minh các áp lực nhân sinh (thu hái dược liệu, bẫy bắt động vật, lấy củi, du lịch sinh thái tự phát) có tần suất xuất hiện cao nhất ($>5\text{ điểm/km}$) tại các kiểu phụ rừng vùng đệm ngoài (I.Np2-1, I.Np2-2) và vùng đệm trong (I.Np1-4). Tuy nhiên, địa hình karst hiểm trở (vách đá dựng đứng, độ dốc $>35^\circ$) đóng vai trò như một rào cản vật lý tự nhiên, bảo toàn được cấu trúc của kiểu phụ nguyên sinh I.Đk1 tại vùng lõi bảo vệ nghiêm ngặt.
       [KẾT QUẢ ĐÁNH GIÁ PHỤC HỒI RỪNG KARST CÁT BÀ - THANG ĐIỂM 5 SAO SER]
                                Mật độ cây gỗ
                                     5 ★
                                  Mật độ TSTV
                                      1 ★
  [Chú giải: I.Đk1 đạt chuẩn 5★ toàn diện; I.Np1 đạt 3-4★; I.Np2 suy thoái nghiêm trọng 1-2★]

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung luận cứ thực chứng khẳng định mô hình diễn thế phục hồi trên karst chịu sự kiểm soát kép (Dual-control mechanism) của vi lập địa thổ nhưỡng và áp lực kinh tế - xã hội vùng đệm, hoàn thiện khung lý thuyết lâm học cận tự nhiên cho vùng karst nhiệt đới.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình ứng dụng kết hợp giữa phân tích đa biến PCA, phân cấp cụm AHC Ward và thang đánh giá phục hồi 5 sao SER, cung cấp bộ công cụ định lượng có thể chuyển giao cho hệ thống 34 vườn quốc gia tại Việt Nam.
  • Về ứng dụng thực tiễn: Định hình chính xác danh mục loài cây mục đích bản địa phù hợp cho phục hồi rừng karst Cát Bà: Lát hoa (Chukrasia tabularis), Re hương (Cinnamomum parthenoxylon), Nghiến (Excentrodendron tonkinense), Kim giao (Nageia fleuryi), Trám đen (Canarium pimela), kết hợp các loài lâm sản ngoài gỗ có giá trị cao như Mây nếp, Song mật, Cọc giậu.
  • Về chính sách và quản trị: Cung cấp cơ sở khoa học trực tiếp để thể chế hóa các quy định tại Thông tư 29/2018/TT-BNNPTNT và Thông tư 17/2022/TT-BNNPTNT của Bộ Nông nghiệp và Phát triển nông thôn. Luận án đề xuất lộ trình tích hợp: chuyển đổi từ giao khoán bảo vệ rừng hình thức sang cơ chế đồng quản lý rừng (Co-management), chia sẻ lợi ích từ dịch vụ môi trường rừng và du lịch sinh thái carbon thấp.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 hạn chế khoa học xuất phát từ điều kiện khách quan:

  1. Thiết kế chuỗi thời gian thay thế (Space-for-time substitution): Do hạn chế về thời gian đào tạo tiến sĩ, nghiên cứu sử dụng phương pháp điều tra đối chứng không gian giữa các trạng thái thoái hóa thay vì theo dõi liên tục trên hệ thống ô định vị dài hạn (Permanent Sample Plots - PSP) qua nhiều thập kỷ.
  2. Khoảng trống về vi sinh vật và nấm rễ nội cộng sinh (Mycorrhizae): Chưa phân tích sâu vai trò của hệ vi sinh vật đất và nấm rễ trong việc hỗ trợ cây gỗ bản địa hấp thụ dinh dưỡng trên các kẽ đá vôi nghèo kiệt.
  3. Biến số khí hậu cực đoan: Chưa lượng hóa được tác động tức thời của các hiện tượng thời tiết cực đoan (bão biển nhiệt đới mạnh, hạn hán kéo dài do El Niño) đến tỷ lệ chết của cây con mới tái sinh.
  4. Phạm vi không gian: Dữ liệu tập trung chủ yếu tại đảo Cát Bà, cần thận trọng khi ngoại suy tuyệt đối cho các hệ sinh thái karst nội địa vùng núi cao Tây Bắc (Hà Giang, Sơn La) nơi có biên độ nhiệt độ mùa đông khắc nghiệt hơn.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Thiết lập mạng lưới ô định vị PSP $1\text{ ha}$ theo tiêu chuẩn quốc tế CTFS-ForestGEO trên núi đá vôi Cát Bà để giám sát động thái diễn thế 10 - 20 năm.
  • Nghiên cứu cơ chế sinh lý thích ứng chịu hạn và mạng lưới nấm rễ mycorhizae của các loài cây gỗ cực kỳ nguy cấp (Nghiến, Kim giao).
  • Đánh giá khả năng hấp thụ và lưu giữ carbon (Carbon sequestration) của thảm thực vật và đất karst trong thị trường tín chỉ carbon rừng.
  • Ứng dụng công nghệ viễn thám LiDAR đa thời gian kết hợp trí tuệ nhân tạo (AI) để giải đoán sinh khối và theo dõi biến động cấu trúc rừng karst tự động.

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

  • Tác động học thuật: Cung cấp cơ sở dữ liệu định lượng gốc về 7 kiểu phụ rừng và 21 QXTVR karst, đóng góp nguồn trích dẫn quan trọng cho các nghiên cứu lâm sinh nhiệt đới, sinh thái học bảo tồn và phục hồi cảnh quan rừng (FLR) tại khu vực Châu Á - Thái Bình Dương.
  • Tác động ngành lâm nghiệp: Cung cấp sổ tay hướng dẫn kỹ thuật phục hồi rừng trên núi đá vôi cho Chi cục Kiểm lâm Hải Phòng và Ban quản lý VQG Cát Bà, trực tiếp nâng cao tỷ lệ sống của cây trồng phục hồi từ $<50%$ lên $>75%$ nhờ chọn đúng vi lập địa và loài cây bản địa.
  • Tác động chính sách: Cung cấp luận cứ vững chắc cho việc thực hiện Quy hoạch tổng thể bảo tồn đa dạng sinh học VQG Cát Bà đến năm 2030, hỗ trợ thực thi cam kết quốc gia về phục hồi các hệ sinh thái bị suy thoái theo Khung Đa dạng Sinh học Toàn cầu Kunming-Montreal (GBF).
  • Tác động xã hội và sinh kế: Định hướng phát triển các mô hình trồng cây dược liệu bản địa (Hà thủ ô đỏ, Bình vôi, Cốt toái bổ) dưới tán rừng khoanh nuôi, tạo nguồn sinh kế bền vững cho hơn 1.500 hộ dân vùng đệm, giảm thiểu áp lực khai thác rừng trái phép.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận sinh thái định lượng kết hợp R-scripts (PCA, AHC Ward) và quy trình điều tra ô lồng Whittaker cải tiến trên địa hình hiểm trở.
  • Nhà quản lý Khu bảo tồn & Vườn Quốc gia: Nhận chuyển giao bản đồ phân vùng mức độ phục hồi và các gói giải pháp kỹ thuật lâm sinh được may đo cho từng kiểu phụ rừng.
  • Các nhà hoạch định chính sách Lâm nghiệp: Sở hữu bộ cứ liệu thực chứng để điều chỉnh chính sách chi trả dịch vụ môi trường rừng (DVMTR) và định mức kinh tế kỹ thuật khoanh nuôi xúc tiến tái sinh trên đất núi đá.
  • Cộng đồng dân cư địa phương: Hưởng lợi từ sự gia tăng dịch vụ hệ sinh thái (nguồn nước ngầm karst trong sạch, giảm thiểu xói mòn lũ quét) và cơ hội phát triển du lịch sinh thái dựa vào cộng đồng.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào? Trả lời: Đóng góp độc đáo nhất là việc tích hợp thành công Thuyết diễn thế sinh thái (Hobbs & Walker, 2007) với Thuyết thảm thực vật sinh thái phát sinh (Thái Văn Trừng, 1978) trên nền tảng lập địa karst biển đảo. Luận án đã mở rộng lý thuyết khi chứng minh rằng: trên lập địa karst cực kỳ khắc nghiệt, sự phục hồi của thảm thực vật không đơn thuần diễn ra theo trật tự thời gian mà chịu sự kiểm soát mang tính quyết định của "ngưỡng thổ nhưỡng mao dẫn" và "rào cản nhân sinh". Việc lượng hóa thành công 7 kiểu phụ và 21 QXTVR đã chuyển đổi cách phân loại lâm học từ mô tả hình thái tĩnh sang mô hình động lực học sinh thái thích ứng.

2. Đột phá về mặt phương pháp luận so với các nghiên cứu trước đây là gì? Trả lời: So với các nghiên cứu lâm học truyền thống tại Việt Nam (như Hoàng Văn Thập, 2008 chỉ thống kê mô tả mật độ và công thức tổ thành cơ bản; Phan Thanh Lâm, 2016 chỉ phân loại theo quy phạm QPN 684), luận án tạo đột phá qua ba điểm:

  • Ứng dụng hệ thống ô lồng đa quy mô Whittaker (1.000 m² - 100 m² - 4 m²) thu thập đồng thời 3 tầng cấu trúc thực vật.
  • Sử dụng toán tin sinh học đa biến (PCA trên R vegan và AHC Ward trên R Factoextra) xử lý đồng thời 16 biến sinh thái - xã hội.
  • Ứng dụng khung chuẩn mực quốc tế 5 sao SER (McDonald et al., 2016) trên biểu đồ Radar để chuẩn hóa mức độ phục hồi rừng.

3. Phát hiện nào trong luận án gây bất ngờ nhất và được giải thích bằng dữ liệu thực nghiệm ra sao? Trả lời: Phát hiện bất ngờ nhất là: Mô hình khoanh nuôi xúc tiến tái sinh có trồng bổ sung lại có số lượng loài tái sinh trong quần xã thấp hơn mô hình khoanh nuôi bảo vệ thuần túy (13 - 15 loài so với 16 - 20 loài, trích dẫn số liệu kiểm chứng từ Triệu Thái Hưng, 2021). Dữ liệu thực nghiệm giải thích rằng: hoạt động can thiệp lâm sinh (dọn thảm tươi, xới xáo đất để trồng dặm) trên nền đất karst mỏng đã vô tình phá vỡ cấu trúc vi sinh cảnh tự nhiên và tầng rêu giữ ẩm, đồng thời các loài cây trồng bổ sung cạnh tranh ánh sáng và dinh dưỡng, ức chế sự nảy mầm của các loài cây bản địa thứ sinh khác.

4. Quy trình nghiên cứu có cung cấp đầy đủ giao thức để lặp lại (Replication protocol) không? Trả lời: Hoàn toàn đầy đủ và minh bạch. Luận án mô tả chi tiết:

  • Tọa độ và sơ đồ 13 tuyến khảo sát, 63 OTC (50 x 20 m) và 16 OTC mô hình.
  • Tiêu chuẩn kích thước đo đếm cây ($D_{1.3} \ge 6\text{ cm}$, cây nhỏ $D_{1.3} < 6\text{ cm}$ và $H_{vn} > 2\text{ m}$, cây tái sinh $H_{vn} \le 2\text{ m}$).
  • Quy trình phân tích đất tuân thủ tuyệt đối các Tiêu chuẩn Quốc gia (TCVN 6860:2001, TCVN 5979:2007, TCVN 4050-85, TCVN 6498:1999).
  • Mã lệnh phân tích PCA, AHC và Quartile tính điểm SER trên phần mềm R 4.0 được chuẩn hóa, cho phép các nhà khoa học khác tái lập hoàn toàn quy trình phân tích.

5. Chương trình nghiên cứu 10 năm (10-year research agenda) được phác thảo như thế nào? Trả lời: Chương trình 10 năm tập trung vào 3 trụ cột:

  • Giai đoạn 2024 - 2027: Xây dựng 03 ô định vị dài hạn PSP diện tích $1\text{ ha}$ tại 3 phân khu chức năng của VQG Cát Bà; giải mã hệ microbiome và nấm rễ rhexiphytic trên nền đá vôi.
  • Giai đoạn 2027 - 2030: Đánh giá sinh khối ngầm và động thái tích lũy carbon của rừng phục hồi karst dưới tác động của biến đổi khí hậu toàn cầu.
  • Giai đoạn 2030 - 2034: Ứng dụng công nghệ AI và Digital Twin (Bản sao số sinh thái) để mô phỏng và dự báo rủi ro suy thoái rừng karst biển đảo trên toàn tuyến duyên hải Bắc Bộ.

Kết luận

  1. Xác lập cơ sở khoa học toàn diện: Luận án đã xây dựng thành công luận cứ khoa học có hệ thống và thực chứng vững chắc cho công tác phục hồi rừng tự nhiên trên núi đá vôi tại VQG Cát Bà thông qua sự kết hợp giữa sinh thái học định lượng và hệ thống kinh tế - xã hội nhân văn.
  2. Phân loại chính xác 7 kiểu phụ rừng và 21 QXTVR: Định danh và lượng hóa chi tiết các đặc trưng lâm học, trắc lượng lâm phần, đặc điểm lý hóa tính đất và mức độ đa dạng thực vật cho 7 kiểu phụ rừng karst (từ kiểu phụ nguyên sinh I.Đk1 đến kiểu phụ thoái hóa sau nương rẫy I.Np2-2).
  3. Đổi mới phương pháp luận đánh giá phục hồi: Ứng dụng thành công ma trận tương quan đa biến PCA, mô hình phân cấp cụm kết tụ AHC Ward và thang đo 5 sao quốc tế SER để lượng hóa chính xác mức độ phục hồi sinh thái của từng trạng thái rừng.
  4. Giải pháp lâm sinh "may đo" theo vi lập địa: Đề xuất hệ thống giải pháp kỹ thuật lâm sinh phân hóa: ưu tiên bảo vệ nghiêm ngặt đối với kiểu phụ giàu I.Đk1; khoanh nuôi xúc tiến tái sinh tự nhiên thuần túy đối với kiểu phụ trung bình và nghèo (I.Np1-1, I.Np1-2); áp dụng khoanh nuôi có trồng bổ sung loài cây bản địa mục đích (Lát hoa, Re hương, Kim giao) kết hợp cải tạo đất đối với các kiểu phụ thoái hóa nặng (I.Np1-3, I.Np1-4, I.Np2-1, I.Np2-2).
  5. Mô hình đồng quản lý và phát triển sinh kế bền vững: Hoàn thiện cơ chế phối hợp đa bên giữa Ban quản lý VQG, Kiểm lâm và cộng đồng 6 xã vùng đệm, gắn kết hoạt động bảo tồn rừng karst với chi trả dịch vụ môi trường rừng và khai thác bền vững lâm sản ngoài gỗ.
  6. Mở ra các hướng nghiên cứu sinh thái karst tiên phong: Đặt nền móng vững chắc cho các nghiên cứu tiếp theo về sinh thái học chức năng (Functional trait ecology), động thái carbon rừng karst biển đảo và quản trị thích ứng hệ sinh thái trước bối cảnh biến đổi khí hậu toàn cầu.