Tổng quan về luận án

Hệ sinh thái rừng (HSTR) trên núi đá vôi (karst) tại Vườn Quốc gia (VQG) Cát Bà là đại diện đặc thù và có giá trị bảo tồn bậc nhất cho thảm thực vật nhiệt đới biển đảo miền Bắc Việt Nam. Với tổng diện tích tự nhiên 17.362,96 ha, trong đó đất quy hoạch lâm nghiệp chiếm 10.443,6 ha (60,2%), diện tích rừng tự nhiên đạt 5.891,69 ha (chiếm 56,6% đất lâm nghiệp). Tuy nhiên, hiện trạng bảo tồn đang đối mặt với thách thức to lớn khi rừng nguyên sinh chỉ còn duy trì 1.056,27 ha (17,3% diện tích rừng), trong khi rừng thứ sinh nghèo và nghèo kiệt chiếm tới 4.856,09 ha, tương đương gần 70% tổng diện tích rừng tự nhiên (VQG Cát Bà, 2020). Thảm thực vật rừng tự nhiên thuộc kiểu rừng kín lá rộng thường xanh mưa ẩm nhiệt đới chịu áp lực phân hóa mạnh mẽ do đặc thù lập địa xương xẩu, tầng đất mỏng, thẩm thấu cao và tính dễ tổn thương sinh thái sâu sắc. Khi bị suy thoái bởi khai thác kiệt hoặc nương rẫy, khả năng tự phục hồi của rừng karst diễn ra vô cùng chậm chạp và khó khăn so với rừng trên núi đất.

Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại ở việc thiếu vắng các cơ sở khoa học định lượng toàn diện về cấu trúc lâm phần, quy luật tái sinh tự nhiên và động thái phục hồi lý - hóa tính đất trong mối tương tác đa chiều với các nhân tố sinh thái - xã hội (SES - Social-Ecological Systems). Các can thiệp kỹ thuật lâm sinh trước đây phần lớn mang tính khái quát, thiếu cơ sở phân cấp lập địa và vi khí hậu vi mô trên nền địa hình chia cắt phức tạp. Luận án tiến sĩ lâm nghiệp của nghiên cứu sinh Lê Hồng Liên dưới sự hướng dẫn khoa học của PGS. Bùi Thế Đồi (Trường Đại học Lâm nghiệp, 2024) đã giải quyết triệt để khoảng trống này thông qua việc thiết lập hệ thống câu hỏi và giả thuyết nghiên cứu tường minh:

  1. Câu hỏi nghiên cứu 1 (RQ1): Đặc trưng cấu trúc tầng cây cao, lớp cây tái sinh tự nhiên và tính chất thổ nhưỡng biến đổi như thế nào giữa các trạng thái diễn thế từ rừng nguyên sinh giàu đến các kiểu phụ thứ sinh thoái hóa?
    Giả thuyết 1 (H1): Mức độ suy thoái do con người tỷ lệ thuận với sự suy giảm chỉ số đa dạng loài, phá vỡ cấu trúc không gian nhiều tầng và làm suy thoái nghiêm trọng hàm lượng mùn hữu cơ cũng như khả năng giữ nước của đất karst.
  2. Câu hỏi nghiên cứu 2 (RQ2): Mối quan hệ tương tác giữa các nhân tố sinh thái (địa hình, đất, cấu trúc lâm phần) và nhân tố kinh tế - xã hội (cự ly tiếp cận, áp lực sinh kế của cộng đồng bản địa) chi phối động thái phục hồi rừng ra sao?
    Giả thuyết 2 (H2): Các yếu tố xã hội và cự ly nhân tác đóng vai trò biến điều chỉnh quyết định quỹ đạo phục hồi hơn là các biến số lập địa tự nhiên đơn lẻ.
  3. Câu hỏi nghiên cứu 3 (RQ3): Căn cứ phân chia mức độ phục hồi sinh thái nào cho phép xây dựng gói giải pháp lâm sinh tích hợp chuẩn xác cho từng quần xã thực vật rừng (QXTVR)?
    Giả thuyết 3 (H3): Mô hình đánh giá phục hồi đa chỉ số chuẩn quốc tế kết hợp phân tích đa biến cho phép tối ưu hóa việc phân vùng khoanh nuôi bảo vệ, xúc tiến tái sinh và trồng bổ sung loài bản địa mục đích.

Dựa trên khung lý thuyết diễn thế sinh thái nhiệt đới (Ecological Succession Theory), nguyên lý lâm sinh gần với tự nhiên (Close-to-nature Silviculture) và khung đánh giá phục hồi rừng của Holl & Aide (2011), nghiên cứu được triển khai trên phạm vi toàn diện bao gồm cả vùng lõi (phân khu bảo vệ nghiêm ngặt, phục hồi sinh thái) và vùng đệm VQG Cát Bà. Công trình thu thập dữ liệu qua 2 đợt khảo sát thực địa nghiêm ngặt (2017 - 2018 và 2021 - 2022), phân tích 13 tuyến khảo sát, 63 ô tiêu chuẩn điển hình (1.000 m²/ô), 16 ô tiêu chuẩn đánh giá mô hình lâm sinh, 21 mẫu đất tầng mặt và phỏng vấn 102 chủ thể quản lý cùng hộ gia đình địa phương. Ý nghĩa đột phá của luận án là xác lập thành công cơ sở khoa học định lượng cho 7 kiểu phụ rừng và 21 QXTVR, định hình bản đồ phân bố và cung cấp hệ thống giải pháp kỹ thuật phục hồi thích ứng sinh thái - xã hội bền vững cho quần đảo Cát Bà.

Literature Review và Positioning

Tổng quan y văn quốc tế ghi nhận hệ sinh thái karst bao phủ khoảng 12% diện tích bề mặt đất liền toàn cầu và 20% diện tích đất không đóng băng của Trái đất (Ford và Williams, 2007). Khu vực Đông Nam Á và Nam Trung Quốc tập trung hơn 800.000 km² núi đá vôi với giá trị bảo tồn và tỷ lệ đặc hữu sinh thái rất cao (Williams, 2008). Các nghiên cứu nền tảng của Wang và cộng sự (2019) chỉ rõ rằng thảm thực vật karst phân bố nghiêm ngặt theo vi địa hình từ chân núi (piedmont) lên vách dốc và đỉnh đá, hình thành các "ốc đảo sinh thái" với tầng đất mỏng, nghèo dinh dưỡng và chịu khô hạn sinh lý theo mùa. Trong các nghiên cứu so sánh quốc tế điển hình, Chen và cộng sự (2022) khi phân tích rừng nhiệt đới trên núi đá vôi tại Xishuangbanna (Trung Quốc) đã phát hiện ra hiện tượng quần hợp sinh thái với tính đa dạng loài tổng thể thấp hơn rừng núi đất liền kề, nhưng bị chi phối bởi các loài ưu thế thích nghi cao như Cleistanthus sumatranus và Lasiococca comberi, song song với tỷ lệ rất lớn của các loài cây quý hiếm phân bố cục bộ. Hiện tượng quần xã đơn ưu (monodominant forest) cũng được ghi nhận tương tự tại rừng nhiệt đới núi đá Zaire với loài Gilbertiodendron dewevrei (van der Velden và cộng sự, 2014). Tại miền Nam Trung Quốc, các công trình phục hồi sa mạc hóa đá vôi (rocky desertification) của Tong và cộng sự (2018) cùng Yirdaw và cộng sự (2019) nhấn mạnh rằng việc suy thoái thảm thực vật karst bắt nguồn từ sự kết hợp giữa biến đổi khí hậu khắc nghiệt và áp lực dân số, đòi hỏi sự can thiệp tích hợp giữa sinh thái phục hồi và tái thiết sinh kế cộng đồng.

Trong bức tranh học thuật đó, tồn tại hai cuộc tranh luận lý thuyết lớn chưa có hồi kết:

  1. Tranh luận về chiến lược phục hồi thụ động (passive restoration) đối lập phục hồi chủ động (active restoration): Một trường phái sinh thái học bảo tồn cho rằng các lập địa karst nhạy cảm cần được bảo vệ tuyệt đối để diễn thế tự nhiên diễn ra tự điều chỉnh nhằm bảo tồn toàn vẹn tính đặc hữu và cấu trúc vi mô bản địa (Lamb và Gilmour, 2003; Holl và Aide, 2011). Ngược lại, các nhà lâm sinh học ứng dụng lập luận rằng do tầng đất mặt bị xói mòn mạnh và nguồn hạt giống của các loài cây gỗ lớn mục đích bị suy kiệt, việc chỉ dựa vào tái sinh tự nhiên sẽ dẫn đến bẫy suy thoái (degradation trap) kéo dài hàng thập kỷ, do đó bắt buộc phải có các biện pháp làm giàu rừng và trồng bổ sung loài bản địa tiên phong (Tong và cộng sự, 2018).
  2. Tranh luận về bảo tồn đa dạng sinh học đối lập cải thiện trữ lượng kinh tế: Long và cộng sự (2021) cùng Triệu Thái Hưng (2021) khi nghiên cứu tại Cát Bà đã chỉ ra một mâu thuẫn thực tiễn: mô hình khoanh nuôi xúc tiến tái sinh có trồng bổ sung giúp gia tăng nhanh trữ lượng và mật độ loài mục đích, nhưng lại có xu hướng làm sụt giảm độ phong phú thành phần loài tổng thể (chỉ đạt 13 - 15 loài/quần xã) so với mô hình khoanh nuôi bảo vệ tự nhiên (16 - 20 loài/quần xã).

Tại Việt Nam, các công trình kinh điển của Thái Văn Trừng (1978) đã phân loại rừng núi đá vôi vào "kiểu phụ thổ nhưỡng kiệt nước trên đất đá vôi xương xẩu (Đk)", trong khi Trần Ngũ Phương (1970) xếp vào đai rừng nhiệt đới mưa mùa lá rộng thường xanh ưu thế Nghiến (Excentrodendron tonkinense). Tiếp nối các nghiên cứu của Trần Hữu Viên (2004, 2005), Trần Văn Con (2006, 2015) và Phạm Xuân Hoàn (2002) về lâm sinh "gần với tự nhiên", luận án của Lê Hồng Liên đã định vị chính xác vị thế học thuật của mình. Công trình đã vượt qua cách tiếp cận mô tả tĩnh hay điều tra kiểm kê đơn thuần của các nghiên cứu trước đó tại Cát Bà (như Hoàng Văn Thập, 2008) để tiến tới một khung phân tích định lượng động thái lâm phần, liên kết ma trận cấu trúc - tái sinh - lập địa với hệ thống sinh thái - xã hội dưới chuẩn mực đánh giá phục hồi quốc tế.

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

                        HỆ THỐNG SINH THÁI - XÃ HỘI (SES)
[ÁP LỰC KINH TẾ - XÃ HỘI]                                     [ĐIỀU KIỆN LẬP ĐỊA VÀ KHÍ HẬU]
- Lịch sử canh tác nương rẫy                                 - Tỷ lệ đá lộ đầu, độ dày tầng đất
- Cự ly tiếp cận, áp lực sinh kế                             - Độ dốc, vi địa hình núi đá vôi
                  [TIẾN TRÌNH DIỄN THẾ VÀ TƯƠNG TÁC ĐA BIẾN]
             (Phân tích PCA ma trận sinh thái - trắc lượng đa cấp)
[CẤU TRÚC LÂM PHẦN TẦNG CÂY CAO]                             [ĐỘNG THÁI LỚP CÂY TÁI SINH]
- Mật độ (N), Tiết diện ngang (G)                            - Mật độ, phân cấp phẩm chất (TSTV)
- Trữ lượng (M), Chỉ số quan trọng (IVI)                     - Tỷ lệ chồi/hạt, chỉ số đa dạng Rényi
                 [ĐÁNH GIÁ MỨC ĐỘ PHỤC HỒI & PHÂN CẤP CAN THIỆP]
                 (Chuẩn hóa 5 sao McDonald et al. 2016 trên Radar)
[PHÂN KHU BẢO VỆ NGHIÊM NGẶT]   [PHÂN KHU PHỤC HỒI SINH THÁI]      [VÙNG ĐỆM TRONG & NGOÀI]
  Phục hồi thụ động (I.Đk1,        Khoanh nuôi xúc tiến tái           Khoanh nuôi xúc tiến tái
     I.Np1-1, I.Np1-2)            sinh có trồng bổ sung (I.Np1-3)    sinh + Trồng làm giàu (I.Np1-4,
                                                                          I.Np2-1, I.Np2-2)

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

Luận án mở rộng và làm sâu sắc lý thuyết diễn thế sinh thái (Ecological Succession Theory) trong sinh thái học nhiệt đới bằng việc chứng minh rằng: "các quy luật thể hiện trong cấu trúc, tổ thành rừng hiện tại là sự phản chiếu của tổng hợp các tác động trong quá khứ". Công trình cụ thể hóa hệ thống phân loại của Thái Văn Trừng (1978) thành mô hình phân loại chi tiết bao gồm 7 kiểu phụ rừng và 21 QXTVR điển hình tại VQG Cát Bà, chia thành 2 nhóm lớn:

  • Kiểu phụ thổ nhưỡng kiệt nước nguyên sinh ít bị tác động (ký hiệu I.Đk1, làm đối chứng chuẩn sinh thái).
  • 6 kiểu phụ thứ sinh nhân tác phản ánh các nấc thang suy thoái diễn thế khác nhau:
    • I.Np1-1: Thứ sinh nhân tác trên đất đá vôi xương xẩu (trạng thái trung bình);
    • I.Np1-2: Thứ sinh nhân tác phục hồi sau khai thác chọn (rừng nghèo);
    • I.Np1-3: Thứ sinh nhân tác phục hồi sau khai thác mạnh (rừng nghèo);
    • I.Np1-4: Thứ sinh nhân tác phục hồi sau khai thác kiệt (vùng đệm trong);
    • I.Np2-1: Thứ sinh nhân tác thoái hóa chân núi sau khai thác (rừng nghèo kiệt);
    • I.Np2-2: Thứ sinh nhân tác thoái hóa sau nương rẫy (rừng nghèo kiệt).

Về mặt lý thuyết lâm sinh học, nghiên cứu đã củng cố lý luận "lâm học gần với tự nhiên" của Phạm Xuân Hoàn (2002) bằng cách xác lập các thông số định lượng về cấu trúc tổ thành, đường kính ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), và chỉ số quan trọng (IVI) của các loài cây ưu thế bản địa (như Chẹo tía, Côm tầng, Lòng mang, Trâm núi, Ô rô, Sấu, Hồng tùng, Kháo vàng). Luận án khẳng định sự chuyển đổi mô hình (paradigm shift) từ quản lý lâm phần đơn tầng, thuần loại sang điều tiết lâm phần đa tầng, hỗn loài, khác tuổi dựa trên cơ chế tự tái sinh của các loài cây mục đích trong các hốc đất karst.

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

Khung phân tích của luận án tích hợp liên ngành giữa 3 trụ cột lý thuyết:

  1. Lý thuyết phục hồi thảm thực vật (Lamb & Gilmour, 2003; Holl & Aide, 2011) với sự phân định ranh giới giữa 3 hành động: cải tạo (reclamation), khôi phục (restoration) và phục hồi (rehabilitation).
  2. Khung phân tích Hệ thống Sinh thái - Xã hội (Social-Ecological Systems - SES): Đặt cấu trúc thảm thực vật trong mối tương tác động với các biến số kinh tế - xã hội (lịch sử canh tác, mật độ dân cư, mức độ phụ thuộc lâm sản ngoài gỗ, chính sách giao khoán rừng).
  3. Trắc lượng sinh thái định lượng đa biến (Multivariate Quantitative Ecology): Sử dụng Phân tích thành phần chính (PCA) để xử lý ma trận đa chiều giữa 16 biến sinh thái - lập địa - xã hội với các kiểu phụ rừng, kết hợp Phân cấp cụm kết tụ (AHC) theo phương pháp liên kết Ward và khoảng cách Euclidean nhằm phân nhóm khách quan mức độ suy thoái và tương đồng.

Ranh giới áp dụng (Boundary conditions): Khung phân tích được xác lập tối ưu cho kiểu phụ thổ nhưỡng kiệt nước núi đá vôi nhiệt đới ven biển và hải đảo thuộc vùng khí hậu nhiệt đới gió mùa Đông Bắc Bộ, nơi có lượng mưa thay đổi theo mùa rõ rệt và chịu tác động kép của địa mạo karst và chế độ hải đảo.

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng lập trường nhận thức luận Thực chứng định lượng (Positivism) kết hợp Hiện thực phê phán (Critical Realism) trong khuôn khổ tiếp cận hệ thống sinh thái xã hội. Thiết kế nghiên cứu đa cấp độ (multi-level design) được triển khai qua các cấp:

  • Cấp độ vĩ mô: Toàn bộ diện tích 17.362,96 ha VQG Cát Bà (phân khu bảo vệ nghiêm ngặt, phục hồi sinh thái, vùng đệm trong và ngoài).
  • Cấp độ trung mô: 7 kiểu phụ rừng và 21 QXTVR điển hình.
  • Cấp độ vi mô: Hệ thống ô tiêu chuẩn (OTC) cố định và tạm thời, ô thứ cấp và các điểm lấy mẫu đất đại diện.

Dung lượng mẫu được xác định chính xác: 13 tuyến khảo sát thực địa chính (với các tuyến phụ 2 - 12 km), 63 ô tiêu chuẩn điển hình diện tích 1.000 m² (50 x 20 m), 16 ô tiêu chuẩn đánh giá 4 mô hình kỹ thuật lâm sinh (4 mô hình x 4 ô lặp lại), 21 mẫu đất tầng mặt (0 - 30 cm) phân tích 10 chỉ tiêu lý hóa, và 102 đối tượng tham gia phỏng vấn bán cấu trúc (gồm cán bộ quản lý VQG, kiểm lâm địa bàn 5 trạm, đại diện 6 cộng đồng và 60 hộ gia đình nhận khoán bảo vệ rừng).

Quy trình nghiên cứu rigorous

Thiết kế cấu trúc ô tiêu chuẩn kế thừa và cải tiến từ phương pháp ô mẫu biến đổi của Whittaker (Stohlgren và cộng sự, 1995) nhằm thích ứng hoàn hảo với địa hình núi đá dốc:

  • Ô lớn A (1.000 m² - 50 x 20 m): Điều tra toàn bộ cây gỗ có $D_{1.3} \ge 6\text{ cm}$. Đo đếm $D_{1.3}$ bằng thước đo dây chuyên dụng (chính xác đến mm), $H_{vn}$ bằng thước đo cao laser Vertex (chính xác đến 0,5 m), giám định tên loài và phân cấp chất lượng cây thành 3 cấp (tốt/triển vọng, xấu/không triển vọng, trung bình). Độ tàn che được đo đạc chính xác bằng phương pháp phân tích ảnh chụp tán vuông góc tại 5 điểm (tâm và 4 góc cách tâm 5 m).
  • Ô thứ cấp B (100 m² tại tâm ô A): Điều tra các cá thể cây gỗ nhỏ có $D_{1.3} < 6\text{ cm}$ và $H_{vn} > 2\text{ m}$.
  • 20 ô dạng bản C (bố trí phân tán trong ô A): Đo đếm cây tái sinh có $D_{1.3} < 6\text{ cm}$ và/hoặc $H_{vn} \le 2\text{ m}$, phân loại theo nguồn gốc phát sinh (tái sinh hạt hoặc chồi), đánh giá phẩm chất và mô tả tỷ lệ che phủ của tầng cây bụi, thảm tươi.

Tam giác đạc dữ liệu (Data triangulation) được thực hiện chặt chẽ qua việc đối chiếu giữa số liệu kiểm đếm thực địa, dữ liệu bản đồ tài nguyên rừng kế thừa từ đề tài cấp Nhà nước ĐTĐL.CN-25/17 của Viện Nghiên cứu Lâm sinh, và dữ liệu điều tra xã hội học từ 102 phiếu phỏng vấn. Hệ thực vật được định danh theo hệ thống phân loại của APG và Danh lục các loài thực vật Việt Nam (Nguyễn Tiến Bân, 2003 - 2005), phân loại dạng sống theo Raunkiaer (1934) với 14 nhóm, và phân loại công dụng theo 11 nhóm chuyên biệt. Độ tin cậy và tính hợp thức của mẫu thổ nhưỡng được bảo đảm thông qua quy chuẩn phân tích mẫu đất hỗn hợp tại phòng thí nghiệm chuyên ngành Đại học Lâm nghiệp theo hệ thống TCVN hiện hành (TCVN 6860:2001, TCVN 5979:2007, TCVN 4050-85, TCVN 6498:1999, TCVN 8567:2010).

Data và phân tích

Số liệu lâm học và đa dạng sinh thái được xử lý trên nền tảng ngôn ngữ lập trình thống kê R (R version 4.0.0) với các gói phần mềm chuyên dụng vegan và factoextra:

  • Tính toán các chỉ số đa dạng loài cấu trúc: Chỉ số Margalef ($d_1$), Menhinick ($d_2$), Simpson ($D$), Shannon ($H'$), chỉ số tương đồng thành phần loài Sorensen/Jaccard ($SI$).
  • Đặc biệt, luận án áp dụng chuỗi chỉ số đa dạng Rényi ($H_\alpha$) với tham số bậc $\alpha \in [0, \infty)$ để vẽ phổ đa dạng sinh học liên tục cho cả tầng cây cao và lớp cây tái sinh, tránh các sai lệch cố hữu của từng chỉ số đơn lẻ.
  • Phân tích tương quan không định hướng PCA trên ma trận đa biến sinh thái - xã hội nhằm trích xuất các giá trị Eigenvalues, tỷ lệ phương sai giải thích ($r^2$) và hệ số tương quan ($r$) trên không gian 2 chiều ($Dim_1, Dim_2$).
  • Áp dụng phương pháp phân cấp cụm kết tụ (AHC) theo thuật toán Ward để phân định ranh giới các cụm lâm phần tương đồng.
  • Đánh giá mức độ phục hồi rừng dựa trên khung tiêu chuẩn 5 sao quốc tế của McDonald và cộng sự (2016). Luận án chia các chỉ số lâm học và xã hội theo hàm tứ phân vị (Quartiles: $<25%, 25-50%, 50-75%, >75%$) tương ứng với thang điểm từ 1 đến 5 sao, biểu diễn trực quan trên đồ thị mạng nhện (Radar chart). Ngưỡng phục hồi sinh thái bền vững đạt chuẩn được xác lập khi toàn bộ các chỉ tiêu thành phần đạt tối thiểu 4 điểm (4 sao).

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

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

  1. Xác lập quy luật phân hóa cấu trúc lâm phần giữa các kiểu phụ rừng:
    Rừng nguyên sinh đối chứng ($I.Đk1$) duy trì mật độ tầng cây cao trung bình 820 - 950 cây/ha, tổng tiết diện ngang $G = 32,4 - 38,6\text{ m}^2/\text{ha}$, trữ lượng $M = 240 - 285\text{ m}^3/\text{ha}$, với các loài ưu thế đạt IVI cao như Lòng mang đa hình (Pterospermum diversifolium), Chẹo tía (Engelhardia chrysolepis), Nhội (Bischofia javanica), Sấu (Dracontomelon duperreanum). Ngược lại, các kiểu phụ thoái hóa sau nương rẫy ($I.Np2-2$) và thoái hóa chân núi ($I.Np2-1$) bị suy giảm nghiêm trọng: mật độ chỉ đạt 280 - 410 cây/ha, $G = 8,2 - 14,5\text{ m}^2/\text{ha}$, $M = 28 - 55\text{ m}^3/\text{ha}$, lâm phần bị đơn giản hóa cấu trúc tầng tán, ưu thế thuộc về các loài cây tiên phong rụng lá hoặc cây bụi gỗ nhỏ như Vả (Ficus auriculata), Bời lời lá tròn (Litsea glutinosa), Sảng nhung (Sterculia velutina).
  2. Quy luật động thái tái sinh tự nhiên và bẫy nghèo kiệt:
    Mật độ cây tái sinh dưới tán dao động từ 4.200 đến 12.500 cây/ha. Tuy nhiên, tỷ lệ cây tái sinh có triển vọng (TSTV - thân thẳng, khỏe, không sâu bệnh) ở các kiểu phụ thoái hóa ($I.Np1-4, I.Np2-1, I.Np2-2$) chỉ chiếm dưới 22 - 30%, phần lớn là cây tái sinh chồi kém bền vững sau các chu kỳ chặt lặp lại. Ở các vách đá dốc và đỉnh núi, tỷ lệ cây mạ chết trong mùa khô đạt trên 65% do kiệt nước cục bộ, chứng minh rằng sự có mặt của lớp phủ thảm tươi và độ tàn che tầng cây cao ($>0,6$) là điều kiện tiên quyết duy trì độ ẩm hạt nảy mầm.
  3. Phát hiện phản trực giác về tác động của mô hình làm giàu rừng đến đa dạng sinh thái:
    So sánh 16 ô tiêu chuẩn mô hình lâm sinh ghi nhận phát hiện quan trọng: mô hình khoanh nuôi xúc tiến tái sinh có trồng bổ sung (KNXTTS-TBS) nâng cao nhanh trữ lượng và mật độ loài mục đích (Lát hoa, Re hương, Nghiến, Kim giao), nhưng lại làm suy giảm đáng kể chỉ số đa dạng Shannon ($H' = 1,82 - 2,15$) so với mô hình khoanh nuôi bảo vệ thụ động ($H' = 2,85 - 3,42$). Sự can thiệp luỗng phát dây leo bụi rậm không đúng kỹ thuật đã vô tình loại bỏ ngân hàng cây tái sinh bản địa của các loài cây tầng dưới.
  4. Động thái suy thoái lý - hóa tính đất karst:
    Phân tích 21 mẫu đất chỉ ra rằng đất dưới rừng nguyên sinh ($I.Đk1$) có hàm lượng mùn hữu cơ (OM) rất cao ($8,4 - 12,2%$), đạm tổng số $N = 0,38 - 0,52%$, CEC đạt $32 - 45\text{ meq}/100\text{g}$ đất, độ xốp $>58%$. Tuy nhiên, ở kiểu phụ thoái hóa sau nương rẫy ($I.Np2-2$), tầng đất mặt bị xói mòn rửa trôi nghiêm trọng: OM sụt giảm còn $2,1 - 3,4%$, CEC giảm còn $12 - 18\text{ meq}/100\text{g}$, dung trọng tăng cao ($1,28 - 1,42\text{ g/cm}^3$) làm đất bị nén chặt, cản trở rễ cây bám vào các khe nứt đá vôi.
  5. Vai trò chi phối của nhân tố kinh tế - xã hội (Kết quả PCA và AHC):
    Trục $Dim_1$ trong phân tích PCA (giải thích $>48,5%$ tổng phương sai) phản ánh tương quan nghịch rất chặt giữa khoảng cách tới khu dân cư/đường giao thông với mức độ suy thoái cấu trúc lâm phần ($r = -0,78, p < 0,001$). Các QXTVR phân bố gần các thung lũng có dân cư sinh sống chịu tần suất tác động chặt trộm gỗ và khai thác lâm sản ngoài gỗ $>5\text{ điểm/km}$, khiến quỹ đạo phục hồi bị chệch hướng hoàn toàn so với diễn thế tự nhiên.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm vững chắc xác lập mô hình diễn thế phục hồi rừng karst nhiệt đới ven biển, chứng minh sự tương tác phi tuyến tính giữa vi địa hình, tính chất thổ nhưỡng và áp lực nhân sinh trong việc định hình cấu trúc quần thể thực vật.
  • Về phương pháp luận: Chuẩn hóa quy trình điều tra OTC Whittaker biến đổi trên địa hình karst dốc đứng và hoàn thiện khung đánh giá phục hồi 5 sao McDonald áp dụng cho hệ thống rừng đặc dụng Việt Nam.
  • Về kỹ thuật lâm sinh ứng dụng: Đề xuất gói giải pháp kỹ thuật theo từng nhóm đối tượng:
    • Đối với nhóm kiểu phụ $I.Đk1, I.Np1-1, I.Np1-2$ (vùng lõi, bảo vệ nghiêm ngặt): Áp dụng giải pháp bảo tồn nguyên vị, bảo vệ nghiêm ngặt, nghiêm cấm mọi tác động cơ học, phát huy tối đa năng lực tự phục hồi.
    • Đối với kiểu phụ $I.Np1-3, I.Np1-4$ (rừng nghèo, vùng đệm trong): Áp dụng khoanh nuôi xúc tiến tái sinh kết hợp trồng bổ sung (mật độ 300 - 500 cây/ha) bằng các loài bản địa ưu thế chịu hạn và ưa canxi như Lát hoa (Chukrasia tabularis), Trám đen (Canarium pimela), Re hương (Cinnamomum parthenoxylon), Nghiến (Excentrodendron tonkinense), Kim giao (Nageia fleuryi).
    • Đối với kiểu phụ thoái hóa $I.Np2-1, I.Np2-2$ (chân núi, nương rẫy cũ): Cần thực hiện các biện pháp cải tạo lý hóa tính đất ban đầu bằng cây tiên phong cải tạo đất (cây họ Đậu), kết hợp gieo hạt thẳng vào các hốc đá có đất trước khi đưa cây gỗ mục đích vào trồng.
  • Về chính sách và quản lý: Thiết lập cơ chế đồng quản lý rừng (co-management) giữa Ban quản lý VQG Cát Bà và cộng đồng cư dân 6 xã vùng đệm; lồng ghép sinh kế bảo tồn với du lịch sinh thái và phát triển chuỗi giá trị lâm sản ngoài gỗ (cây dược liệu, mật ong hoa rừng tự nhiên) nhằm giảm thiểu áp lực xâm lấn tài nguyên.

Limitations và Future Research

  1. Hạn chế về dữ liệu chuỗi thời gian định vị (Longitudinal time-series limitation): Do điều kiện nghiên cứu luận án, dữ liệu được thu thập qua 2 đợt khảo sát tĩnh (2017 - 2018 và 2021 - 2022). Mặc dù phương pháp thay thế thời gian bằng không gian (space-for-time substitution) đã được áp dụng để phân tích các pha diễn thế, việc thiếu hệ thống ô tiêu chuẩn định vị giám sát liên tục 10 - 20 năm là một hạn chế trong việc nắm bắt chi tiết tốc độ tăng trưởng đường kính và sinh khối thực tế.
  2. Hạn chế về quan trắc vi khí hậu liên tục: Nghiên cứu chưa thiết lập được hệ thống trạm cảm biến đo tự động liên tục các chỉ số bốc thoát hơi nước, độ ẩm không khí và nhiệt độ bề mặt tán tại các vi địa hình khác nhau (thung lũng so với đỉnh núi dốc).
  3. Phạm vi phân tích sinh học phân tử và khu hệ vi sinh vật đất: Luận án tập trung chuyên sâu vào thực vật bậc cao có mạch và các chỉ tiêu lý hóa tính đất cơ bản, chưa đi sâu phân tích hệ vi sinh vật rễ nấm (mycorrhizae) cộng sinh trên nền đá vôi – một nhân tố tiềm năng chi phối khả năng hấp thu dinh dưỡng của cây tái sinh.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Thiết lập mạng lưới 20 ô tiêu chuẩn định vị vĩnh cửu (Permanent Sampling Plots - PSPs) kích thước 1 ha theo tiêu chuẩn CTFS-ForestGEO tại VQG Cát Bà để theo dõi biến động động thái rừng karst ven biển dưới tác động của biến đổi khí hậu toàn cầu và nước biển dâng.
  • Nghiên cứu cơ chế thích ứng sinh lý học (ecophysiological traits) và động thái trao đổi nước - carbon của các loài cây đặc hữu quý hiếm (Excentrodendron tonkinense, Nageia fleuryi, Livistona halongensis) trong điều kiện khô hạn cực đoan mùa đông.
  • Ứng dụng công nghệ viễn thám LiDAR đa thời gian kết hợp máy bay không người lái (UAV) và mô hình học máy (Machine Learning) để lượng hóa trữ lượng carbon và cấu trúc tán rừng 3D trên toàn bộ quần thể đảo Cát Bà.
  • Nghiên cứu cơ chế cộng sinh nấm rễ và hệ gen chức năng của thảm thực vật karst nhằm sản xuất các chế phẩm vi sinh phục vụ ươm tạo cây con bản địa phục hồi sinh thái.

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

  • Tác động học thuật: Luận án đóng góp hệ thống dữ liệu định lượng gốc phong phú cho ngành Lâm sinh và Sinh thái học nhiệt đới Việt Nam, làm tài liệu tham khảo chuẩn mực cho các chương trình đào tạo sau đại học tại Trường Đại học Lâm nghiệp và các viện nghiên cứu sinh thái. Dự kiến các công bố quốc tế trích xuất từ luận án sẽ đóng góp vào cơ sở dữ liệu Scopus/ISI về bảo tồn hệ sinh thái karst Đông Nam Á.
  • Tác động quản lý và chuyển đổi ngành lâm nghiệp: Cung cấp cơ sở khoa học trực tiếp giúp Ban Quản lý VQG Cát Bà điều chỉnh "Phương án quản lý rừng bền vững giai đoạn 2021 - 2030" và triển khai chuẩn xác các dự án phục hồi rừng ven biển tại thành phố Hải Phòng.
  • Tác động chính sách quốc gia và quốc tế: Cung cấp luận cứ khoa học thực tiễn phục vụ công tác quản lý và bảo tồn Di sản Thiên nhiên Thế giới Vịnh Hạ Long - Quần đảo Cát Bà (được UNESCO công nhận năm 2023), đồng thời hỗ trợ triển khai hiệu quả Thông tư 29/2018/TT-BNNPTNT và Thông tư 17/2022/TT-BNNPTNT của Bộ Nông nghiệp và PTNT về quy định các biện pháp kỹ thuật lâm sinh.
  • Lợi ích xã hội: Tạo cơ sở triển khai các mô hình sinh kế gắn với chi trả dịch vụ môi trường rừng (PES) và du lịch sinh thái dựa vào cộng đồng cho hơn 60 hộ dân nhận khoán bảo vệ rừng tại 6 xã thuộc huyện Cát Hải, nâng cao thu nhập và giảm nghèo bền vững.

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

  • Nghiên cứu sinh và Giảng viên Lâm nghiệp: Tiếp cận hệ thống phương pháp luận sinh thái định lượng tiên tiến (PCA, AHC, chỉ số Rényi, tiêu chuẩn 5 sao McDonald) và các khoảng trống nghiên cứu chuyên sâu về sinh thái học thực vật karst.
  • Nhà khoa học Lâm sinh và Bảo tồn: Kế thừa tập dữ liệu chuẩn về 7 kiểu phụ rừng và danh lục thực vật chuyên khảo để phát triển các nghiên cứu chuyên sâu về động thái rừng và thích ứng biến đổi khí hậu.
  • Ban Quản lý các Vườn Quốc gia và Khu bảo tồn thiên nhiên có núi đá vôi: Áp dụng trực tiếp gói hướng dẫn kỹ thuật lâm sinh phân tầng để quy hoạch phân vùng bảo tồn, thiết kế mô hình khoanh nuôi xúc tiến tái sinh và lựa chọn đúng danh mục loài cây bản địa trồng bổ sung.
  • Cơ quan Hoạch định Chính sách (Cục Lâm nghiệp, Sở NN&PTNT các tỉnh có rừng Karst): Cơ sở dữ liệu khoa học thực chứng để phân bổ nguồn vốn đầu tư công cho các chương trình trồng rừng phòng hộ đặc dụng và bảo tồn đa dạng sinh học.
  • Cộng đồng địa phương vùng đệm Cát Bà: Hưởng lợi từ sự cải thiện dịch vụ hệ sinh thái nguồn nước ngầm, cơ hội việc làm thông qua các dự án trồng rừng bổ sung và hưởng lợi từ du lịch sinh thái 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 luận án là gì và đã mở rộng lý thuyết nào?
Đóng góp lý thuyết độc đáo nhất là việc tích hợp thành công Lý thuyết Diễn thế Sinh thái học (Ecological Succession Theory) với Khung Hệ sinh thái - Xã hội (SES) để giải mã cơ chế phục hồi rừng trên nền địa chất karst ven biển. Luận án đã mở rộng hệ thống phân loại thảm thực vật rừng Việt Nam của Thái Văn Trừng (1978) bằng việc thiết lập chi tiết 7 kiểu phụ rừng và 21 QXTVR định lượng, cung cấp bằng chứng thực nghiệm khẳng định rằng mức độ phục hồi cấu trúc và đa dạng loài của rừng karst bị điều tiết đồng thời bởi vi địa hình thổ nhưỡng và lịch sử tác động nhân sinh trong quá khứ.

2. Điểm đổi mới phương pháp luận nổi bật khi so sánh với ít nhất 2 nghiên cứu tiền nhiệm?
So với các nghiên cứu trước đây tại Cát Bà của Hoàng Văn Thập (2008) chỉ dừng lại ở điều tra mô tả cấu trúc tĩnh đơn giản theo phương pháp lâm học truyền thống, hay nghiên cứu của Long và cộng sự (2021) chỉ đánh giá các chỉ số đa dạng đơn lẻ ($H', D$), luận án đã có 3 đột phá phương pháp luận:

  • Cải tiến phương pháp ô tiêu chuẩn biến đổi của Whittaker (Stohlgren et al., 1995) thiết lập cấu trúc ô phân cấp 3 tầng (A: 1.000 m², B: 100 m², C: 20 ô dạng bản) đo đếm toàn diện từ cây gỗ lớn đến ngân hàng cây tái sinh.
  • Sử dụng chuỗi đa dạng bậc cao Rényi ($H_\alpha$) kết hợp trắc lượng đa biến PCA và AHC trên môi trường R Studio để phân cụm khách quan các trạng thái rừng.
  • Ứng dụng khung đánh giá phục hồi 5 sao quốc tế của McDonald và cộng sự (2016) chuẩn hóa bằng hàm tứ phân vị trên đồ thị Radar.

3. Phát hiện nào gây bất ngờ nhất và được chứng minh bằng dữ liệu nào?
Phát hiện bất ngờ và có ý nghĩa thực tiễn sâu sắc nhất là: mô hình khoanh nuôi xúc tiến tái sinh có trồng bổ sung (KNXTTS-TBS) – vốn luôn được kỳ vọng là tối ưu hóa phục hồi – lại làm suy giảm độ phong phú thành phần loài tổng thể của lâm phần ($H' = 1,82 - 2,15$, chỉ có 13 - 15 loài tham gia tổ thành) so với mô hình khoanh nuôi bảo vệ tự nhiên ($H' = 2,85 - 3,42$, có 16 - 20 loài). Dữ liệu điều tra 16 ô tiêu chuẩn mô hình chứng minh rằng các tác động kỹ thuật luỗng phát dây leo bụi rậm thiếu kiểm soát đã triệt hạ vô tình các loài cây bản địa tầng dưới và làm đảo lộn cấu trúc tái sinh tự nhiên.

4. Luận án có cung cấp quy trình nhân rộng (Replication Protocol) rõ ràng không?
Luận án cung cấp quy trình nghiên cứu và ứng dụng thực tiễn hoàn chỉnh gồm 5 bước chuẩn hóa:

  1. Thu thập dữ liệu thứ cấp và giải đoán ảnh viễn thám phân vùng kiểu phụ rừng;
  2. Thiết lập tuyến khảo sát và hệ thống OTC Whittaker cải tiến (1.000 m²);
  3. Điều tra thu thập đồng thời 4 nhóm biến số: lâm học tầng cây cao, tái sinh lớp dưới, lý hóa tính đất và phỏng vấn cấu trúc kinh tế - xã hội;
  4. Xử lý đa biến PCA, AHC và tính toán chuỗi đa dạng Rényi trên phần mềm R;
  5. Đánh giá mức độ phục hồi theo thang 5 sao McDonald để đề xuất giải pháp kỹ thuật lâm sinh phù hợp.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao?
Chương trình nghiên cứu 10 năm tập trung vào 4 hướng chiến lược:

  • Xây dựng mạng lưới ô định vị vĩnh cửu 1 ha giám sát sinh khối và carbon dài hạn;
  • Nghiên cứu sinh thái sinh lý học thực vật chịu hạn và thích ứng biến đổi khí hậu của các loài đặc hữu;
  • Ứng dụng công nghệ viễn thám LiDAR/UAV và trí tuệ nhân tạo giám sát cấu trúc tán 3D;
  • Phân tích hệ vi sinh vật đất và nấm rễ cộng sinh trên núi đá vôi phục vụ công nghệ sinh học lâm nghiệp.

Kết luận

  1. Xác lập cơ sở khoa học định lượng toàn diện: Luận án đã làm sáng tỏ cấu trúc sinh thái, động thái tái sinh tự nhiên, đa dạng loài và đặc tính lý - hóa tính đất cho 7 kiểu phụ rừng và 21 QXTVR tại VQG Cát Bà.
  2. Chứng minh quy luật tương tác sinh thái - xã hội: Bằng phân tích PCA đa biến ($Dim_1 > 48,5%$ phương sai), công trình chứng minh cự ly nhân tác và áp lực sinh kế của cộng đồng đóng vai trò quyết định quỹ đạo phục hồi của thảm thực vật karst.
  3. Đột phá trong phương pháp đánh giá phục hồi: Ứng dụng thành công chuẩn mực phục hồi quốc tế 5 sao của McDonald và cộng sự (2016) trên đồ thị mạng nhện đa chiều, lượng hóa chính xác khoảng cách sinh thái giữa các trạng thái thứ sinh so với rừng nguyên sinh chuẩn đối chứng ($I.Đk1$).
  4. Hệ thống giải pháp lâm sinh thích ứng phân tầng: Đề xuất chính xác danh mục loài cây bản địa mục đích (Lát hoa, Re hương, Nghiến, Kim giao, Trám đen) và quy trình can thiệp lâm sinh riêng biệt cho từng phân khu chức năng (bảo tồn thụ động tại vùng lõi; khoanh nuôi xúc tiến tái sinh có trồng bổ sung tại vùng đệm trong; cải tạo đất và gieo hạt thẳng tại vùng nương rẫy thoái hóa).
  5. Mở ra các hướng nghiên cứu học thuật mới: Khởi xướng các luồng nghiên cứu chuyên sâu về giám sát sinh khối vĩnh cửu dài hạn, sinh thái sinh lý thực vật karst thích ứng khô hạn cực đoan và ứng dụng viễn thám LiDAR cấu trúc tán tại các khu bảo tồn thiên nhiên biển đảo Việt Nam.