Tổng quan về luận án

Nghiên cứu cấu trúc và động thái thảm thực vật trên hệ sinh thái núi đá vôi (Karst) là một trong những chủ đề trọng tâm của sinh thái học hiện đại do tính chất đặc thù về địa mạo, tính dễ bị tổn thương và mức độ đặc hữu cao của các loài sinh vật. Tại Việt Nam, hệ thống Karst ven biển và hải đảo khu vực thành phố Cẩm Phả, tỉnh Quảng Ninh—nằm trong không gian sinh thái cảnh quan tiếp giáp vịnh Bái Tử Long và vùng đệm di sản thiên nhiên thế giới vịnh Hạ Long—đang chịu áp lực nhân sinh khốc liệt từ hoạt động khai thác than, khai thác đá vôi sản xuất xi măng và quá trình đô thị hóa lấn biển. Luận án tiến sĩ Sinh học của Hoàng Văn Hải (2020) với tiêu đề "Nghiên cứu một số đặc trưng cơ bản và đề xuất biện pháp bảo tồn, phát triển bền vững thảm thực vật trên núi đá vôi ở khu vực thành phố Cẩm Phả, tỉnh Quảng Ninh" đã thiết lập một hệ thống dẫn liệu định lượng toàn diện về cấu trúc không gian, quy luật phân bố lâm phần và động thái tái sinh của thảm thực vật Karst nhiệt đới hải đảo.

+----------------------------------------------------------------------------------------------------+
|                                KHUNG PHÂN TÍCH TỔNG THỂ CỦA LUẬN ÁN                                |
+----------------------------------------------------------------------------------------------------+
|   ĐỊA MẠO VI SINH CẢNH (Thung lũng - Chân núi - Sườn vách đá)                                      |
|   MÔ HÌNH HÓA ĐỊNH LƯỢNG & ĐỘNG THÁI QUẦN XÃ                                                       |
|   ĐÁNH GIÁ THOÁI HÓA & ĐỀ XUẤT GIẢI PHÁP BẢO TỒN BỀN VỮNG                                          |
+----------------------------------------------------------------------------------------------------+

Khoảng trống nghiên cứu (research gap) mà công trình giải quyết xuất phát từ thực tế các nghiên cứu trước đây tại Quảng Ninh (như Nguyễn Tiến Hiệp, 2005; Nguyễn Thế Hưng, 2010) chủ yếu dừng lại ở mức độ kiểm kê danh lục thực vật hoặc mô tả định tính diện mạo ngoài, thiếu vắng các mô hình toán học lượng hóa cấu trúc đường kính ($N/D_{1.3}$), cấu trúc chiều cao ($N/H_{vn}$), tương quan sinh trưởng allometric ($H_{vn}/D_{1.3}$) cũng như động thái tái sinh tự nhiên theo từng đơn vị địa hình vi sinh cảnh (micro-topography).

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Thành phần loài, phổ dạng sống và các yếu tố địa lý của thảm thực vật núi đá vôi Cẩm Phả có đặc trưng phân bố và giá trị bảo tồn như thế nào?
  2. RQ2: Quy luật phân bố mật độ theo cấp đường kính ($N/D_{1.3}$), cấp chiều cao ($N/H_{vn}$) và tương quan sinh trưởng ($H_{vn}/D_{1.3}$) của tầng cây gỗ biến đổi ra sao giữa các sinh cảnh thung lũng, chân núi và sườn vách đá?
  3. RQ3: Khả năng tái sinh tự nhiên của cây gỗ (mật độ, tổ thành loài, nguồn gốc, phẩm chất và cấu trúc không gian) bị chi phối bởi những nhân tố sinh thái nào?
  4. RQ4: Những tác động nhân sinh nào đang đe dọa trực tiếp đến tính bền vững của thảm thực vật Karst và hệ thống giải pháp quản lý nào mang tính khả thi cao nhất?

Hệ thống giả thuyết khoa học được kiểm chứng bao gồm:

  • H1: Cấu trúc mật độ theo cấp đường kính và chiều cao của thảm thực vật Karst tuân theo các hàm phân bố toán học phi tuyến (Weibull, phân bố giảm) phản ánh trạng thái rừng tự nhiên đang trong quá trình phục hồi tích cực.
  • H2: Sự phân hóa vi địa hình (thung lũng, chân núi, sườn vách) tạo ra sự phân tầng tán và biến thiên chỉ số quan trọng (IVI%) có ý nghĩa sinh thái rõ rệt.
  • H3: Quá trình tái sinh tự nhiên của cây gỗ trên đá vôi bị giới hạn nghiêm trọng ở giai đoạn chuyển tiếp chiều cao do hạn chế về tầng đất và độ ẩm cục bộ.

Khung lý thuyết của nghiên cứu tích hợp Lý thuyết phân vùng sinh thái rừng nhiệt đới của Thái Văn Trừng (1978, 2000), Lý thuyết phân loại quần hợp thực vật của Sukachev (1957) và Ellenberg & Mueller-Dombois (1967), cùng Nguyên lý thích nghi hình thái thực vật Karst của Chapman & Yin-Zheng (2002).

Đóng góp đột phá của luận án được lượng hóa qua các chỉ số cụ thể: xác lập danh lục 608 loài thực vật có mạch thuộc 388 chi, 131 họ, phát hiện 12 loài đặc hữu hẹp và 27 loài quý hiếm có nguy cơ tuyệt chủng; xây dựng các hàm tương quan sinh trưởng allometric tối ưu được lựa chọn qua Tiêu chuẩn thông tin Akaike (AIC); khảo sát trên quy mô 50 ô tiêu chuẩn ($2.500\text{ m}^2$/ô, tổng diện tích $125.000\text{ m}^2$) và 18 tuyến điều tra sườn vách tại 34 vị trí địa lý đại diện, tạo cơ sở khoa học cho công tác bảo tồn di sản Karst ven biển Việt Nam.

Literature Review và Positioning

Nghiên cứu thảm thực vật Karst trên thế giới đã trải qua một tiến trình phát triển lâu dài từ việc mô tả thực vật học hình thái sang sinh thái học chức năng và mô hình hóa toán học. Các công trình kinh điển của Pesola (1928, 1955) tại Bắc Âu đã phân chia thực vật thành 3 nhóm sinh thái: cây ưa vôi bắt buộc (obligate calciphiles), cây ưa vôi tùy ý (calciphiles) và cây kỵ vôi (calcifuges). Tại châu Á, các nghiên cứu của Hou (1954), Makino (1983) và Chapman & Yin-Zheng (2002) tại vùng Karst phía Nam Trung Quốc (với diện tích hơn $283.000\text{ km}^2$) đã chỉ rõ tính đa dạng cực cao của các loài đặc hữu cổ xưa (Cycas baguanheensis, Cycas panzhihuaensis) thích nghi đặc biệt với kẽ nứt đá vôi (dạng sống chasmophytes).

+----------------------------------------------------------------------------------------------------+
|                         TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT THẢM THỰC VẬT KARST                        |
+----------------------------------------------------------------------------------------------------+
|   GIAI ĐOẠN CỔ ĐIỂN (1920 - 1970)                                                                  |
|   • Pesola (1928, 1955), Drude (1887): Phân loại nhóm ưa vôi (calciphiles / obligate calciphiles). |
|   • Champion (1936), Sukachev (1957): Khung phân loại quần hệ và sinh thái rừng.                   |
|   GIAI ĐOẠN ĐỊNH LƯỢNG & ĐỘNG THÁI (1970 - 2000)                                                  |
|   • Thái Văn Trừng (1978, 2000), Trần Ngũ Phương (1970): Phân loại thảm thực vật Karst Việt Nam.   |
|   • Ellenberg & Mueller-Dombois (1967), David et al. (1993): Động thái tái sinh & cấu trúc rừng.    |
|   GIAI ĐOẠN HIỆN ĐẠI & MÔ HÌNH HÓA VI SINH CẢNH (2000 - Nay)                                       |
|   • Chapman & Yin-Zheng (2002), Trần Hữu Viên (2004): Thích nghi chasmophyte & phân vùng sinh thái |
|   • HOÀNG VĂN HẢI (2020): Mô hình hóa phân bố phi tuyến (Weibull/AIC) + Động thái tái sinh vi địa  |
|     hình + Tích hợp khung suy thoái Karst ven biển/hải đảo.                                       |
+----------------------------------------------------------------------------------------------------+

Trong bức tranh tổng quan học thuật, tồn tại hai luồng quan điểm tranh luận lớn:

  • Trường phái cân bằng sinh thái (Deterministic / Equilibrium Theory): Đại diện bởi Clements (1916) và được phát triển trong lâm học nhiệt đới bởi Champion (1936) và Trần Ngũ Phương (1970), cho rằng thảm thực vật Karst phát triển hướng tới một quần xã cao đỉnh (climax) ổn định được quyết định bởi chế độ khí hậu vùng, trong đó các loài cây gỗ chiếm ưu thế sinh thái tuyệt đối tạo nên cấu trúc phân tầng cố định.
  • Trường phái phi cân bằng và động thái nhiễu loạn (Non-equilibrium / Disturbance Dynamics): Đại diện bởi Whittaker (1953, 1975) và Grime (1979), khẳng định rằng thảm thực vật trên đá vôi là một phức hệ khảm vi sinh cảnh bị chi phối bởi các yếu tố giới hạn khắc nghiệt về thổ nhưỡng (đất mỏng, rửa trôi kiềm, khô hạn sinh lý) và các đợt nhiễu loạn liên tục, khiến thảm thực vật duy trì trạng thái diễn thế thứ sinh kéo dài mà không đạt tới cao đỉnh khí hậu đồng nhất.

Tại Việt Nam, các nghiên cứu của Trần Ngũ Phương (1970, 2000), Thái Văn Trừng (1978, 2000) và Trần Hữu Viên (2004) đã định hình hệ thống phân loại rừng trên núi đá vôi miền Bắc (Cao Bằng - Lạng Sơn, Tuyên Quang - Hà Giang, Tây Bắc - Bắc Trung Bộ và vùng đảo Đông Bắc). Tuy nhiên, các công trình chuyên sâu về dải Karst ven biển Quảng Ninh phần lớn tập trung tại Vườn quốc gia Cát Bà (1.585 loài - Trung tâm Tài nguyên và Môi trường Lâm nghiệp, 2015) hoặc vịnh Hạ Long (Nguyễn Tiến Hiệp, 2005).

Vị trí học thuật của luận án Hoàng Văn Hải (2020) được xác lập bằng việc lấp đầy khoảng trống dữ liệu thực chứng tại thành phố Cẩm Phả—khu vực chịu tác động kép của khí hậu nhiệt đới gió mùa ven biển và hoạt động công nghiệp khai khoáng cực đoan. Công trình nâng tầm nghiên cứu sinh thái Karst khu vực khi so sánh trực tiếp với hai hệ sinh thái Karst quốc tế tiêu biểu:

  1. Hệ sinh thái Karst Nam Trung Quốc (Guangxi/Guizhou Karst - Chapman & Yin-Zheng, 2002): Có sự tương đồng cao về các chi cổ nhiệt đới (Cycas, Schefflera, Dracaena, Chirita) nhưng Karst Cẩm Phả thể hiện mức độ phân mảnh cao hơn do địa hình đảo ngập triều.
  2. Hệ sinh thái Karst nhiệt đới Đông Nam Á (Gunung Mulu - Malaysia & Ao Phang Nga - Thái Lan): Thảm thực vật Cẩm Phả mang nét độc đáo của đai rừng nhiệt đới mưa mùa có mùa đông lạnh, dẫn đến sự đan xen giữa các yếu tố thực vật bản địa Bắc Việt Nam - Nam Trung Hoa với các luồng thực vật di cư Ấn Độ - Miến Điện và Malaysia - Indonesia.

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

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

Luận án đóng góp vào việc mở rộng và kiểm chứng các lý thuyết sinh thái học nền tảng trong điều kiện lập địa khắc nghiệt:

  • Mở rộng Thuyết Diễn thế Sinh thái (Ecological Succession Theory - Clements, 1916; Whittaker, 1953): Chứng minh rằng diễn thế phục hồi thảm thực vật Karst sau khai thác không diễn ra theo một chu trình tuyến tính đơn giản mà phân hóa thành các quỹ đạo vi sinh cảnh độc lập. Tại thung lũng, diễn thế tiến triển thuận lợi hướng tới quần xã thứ sinh phục hồi (trạng thái IIa, IIb, IIIa1); ngược lại, tại sườn vách dốc đứng, diễn thế bị "khóa" ở trạng thái quần xã tiên phong cây bụi - chasmophyte bền vững do giới hạn tầng đất tích tụ.
  • Phát triển Lý thuyết Dạng sống Thực vật (Raunkiaer’s Life-Form Theory, 1934): Xác lập phổ dạng sống thực vật chuẩn cho vùng Karst ven biển Cẩm Phả, làm sáng tỏ cơ chế thích nghi giải phẫu của nhóm thực vật chồi trên (Phanerophytes - Ph) và thực vật bì sinh, dây leo (Liana/Epiphytes - Lp) trong điều kiện khô hạn sinh lý định kỳ.
  • Lượng hóa Thuyết Sinh thái Allometric (Allometric Scaling Theory): Cung cấp bằng chứng thực nghiệm rằng mối quan hệ sinh trưởng chiều cao - đường kính ($H_{vn}/D_{1.3}$) trên núi đá vôi bị nén lại (compressed scaling) so với rừng trên núi đất, do cây gỗ phải ưu tiên phân bổ sinh khối cho hệ rễ bám đá và cấu trúc cơ học chịu gió bão ven biển.
+----------------------------------------------------------------------------------------------------+
|                         KHUNG MÔ HÌNH HÓA CẤU TRÚC VÀ ĐỘNG THÁI QUẦN XÃ                           |
+----------------------------------------------------------------------------------------------------+
|                                                                                                    |
|    [Vi sinh cảnh Thung lũng]        [Vi sinh cảnh Chân núi]         [Vi sinh cảnh Sườn vách đá]     |
|                           [TỔ THÀNH LOÀI & CHỈ SỐ QUAN TRỌNG (IVI%)]                               |
|                             IVI% = RD% (Mật độ) + RF% (Tần số) + RDo% (Ưu thế)                     |
|                      [MÔ HÌNH ALLOMETRIC Hvn/D1.3 QUA TIÊU CHUẨN AIC]                              |
|                          AIC = 2k - 2*ln(L)  --> Chọn mô hình tối ưu                               |
|                          [ĐỘNG THÁI TÁI SINH TỰ NHIÊN CÂY GỖ]                                      |
+----------------------------------------------------------------------------------------------------+

Các mệnh đề lý thuyết (Theoretical Propositions) được xác lập:

  • Proposition 1: Độ đa dạng loài ($H'$) và chỉ số tương đồng Sorensen ($S$) giữa các thảm thực vật Karst tỷ lệ thuận với mức độ dày của tầng đất dốc tụ và giảm dần theo độ dốc địa hình.
  • Proposition 2: Cấu trúc số cây theo cấp đường kính ($N/D_{1.3}$) tại các thung lũng đá vôi tuân theo hàm phân bố Weibull 3 tham số, phản ánh tính chất ổn định tương đối của chu trình đào thải tự nhiên.
  • Proposition 3: Hiệu quả tái sinh tự nhiên ($R_e$) của cây gỗ bị chi phối chủ yếu bởi tỷ lệ che phủ của tầng tán cây cao ($A_1, A_2$) và độ ẩm vi khí hậu dưới tán, trong đó cây tái sinh gốc hạt chiếm ưu thế áp đảo so với cây chồi.

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

Luận án xây dựng một khung phân tích đa chiều tích hợp 3 trụ cột lý thuyết:

  1. Cổ điển học Thực vật Quần xã (Braun-Blanquet Phytosociology): Đánh giá qua tổ thành và chỉ số IVI%.
  2. Lý thuyết Chiến lược Thực vật C-S-R của Grime (1979): Phân tích cơ chế chống chịu stress khắc nghiệt của thực vật Karst.
  3. Lý thuyết Phân vùng Sinh thái Cảnh quan (Landscape Ecology): Đánh giá mức độ thoái hóa dựa trên ma trận tác động nhân sinh.

Cách tiếp cận phân tích mới nằm ở việc sử dụng đồng thời các thuật toán phi tuyến tính để khớp dữ liệu thực nghiệm lâm phần, sử dụng chỉ số AIC (Akaike's Information Criterion) để lựa chọn mô hình tương quan $H_{vn}/D_{1.3}$ tối ưu, loại bỏ sự chủ quan trong việc áp đặt hàm toán học truyền thống.

Điều kiện biên (boundary conditions) của khung phân tích: Áp dụng chuẩn xác cho hệ sinh thái Karst ven biển nhiệt đới ẩm gió mùa ở độ cao dưới $300\text{ m}$ so với mực nước biển, chịu ảnh hưởng trực tiếp của chế độ thủy triều và khí hậu biển Đông Bắc Việt Nam.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường triết học Thực chứng Phê phán (Critical Realism), kết hợp giữa định lượng thực nghiệm sinh thái (đo đạc kích thước lâm học, phân tích cấu trúc toán học) và tiếp cận kiến tạo xã hội (phỏng vấn nhân học sinh thái, đánh giá ma trận áp lực quản lý).

Thiết kế nghiên cứu đa tầng (Multi-level design) được triển khai qua 4 cấp độ:

  • Cấp độ vùng cảnh quan: Không gian Karst thành phố Cẩm Phả và vịnh Bái Tử Long.
  • Cấp độ địa mạo vi sinh cảnh (3 dạng): Thung lũng núi đá vôi, chân núi đá vôi và sườn vách đá dốc đứng.
  • Cấp độ quần xã (Plot/Transect level): 50 ô tiêu chuẩn diện tích lớn ($2.500\text{ m}^2$) và 18 tuyến khảo sát chuyên dụng.
  • Cấp độ cá thể (Individual level): Toàn bộ cây gỗ có đường kính ngang ngực $D_{1.3} \ge 6\text{ cm}$ và cây tái sinh có chiều cao $H < 2\text{ m}$.
+----------------------------------------------------------------------------------------------------+
|                                 SƠ ĐỒ BỐ TRÍ THIẾT KẾ NGHIÊN CỨU                                    |
+----------------------------------------------------------------------------------------------------+
|                                                                                                    |
|    ĐỊA HÌNH THUNG LŨNG                   ĐỊA HÌNH CHÂN NÚI                  ĐỊA HÌNH SƯỜN VÁCH     |
+----------------------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu tuân thủ các tiêu chuẩn nghiêm ngặt của điều tra sinh thái rừng (Nguyễn Nghĩa Thìn, 2008):

  • Chiến lược lấy mẫu: Thiết lập 50 ô tiêu chuẩn ($50\text{ m} \times 50\text{ m} = 2.500\text{ m}^2$), phân bổ đều: 25 OTC tại thung lũng và 25 OTC tại chân núi. Tại sườn vách đá dựng đứng, thiết kế 18 tuyến điều tra ngẫu nhiên với cự ly tuyến $50 - 100\text{ m}$, bề rộng tuyến $5\text{ m}$. Tổng cộng khảo sát tại 34 địa điểm đại diện trên biển và đất liền.
  • Thu thập dữ liệu lâm học: Đo đường kính ngang ngực ($D_{1.3}$) bằng thước dây chuyên dụng và thước kẹp kính; đo chiều cao vút ngọn ($H_{vn}$) bằng máy đo cao Blume-Leiss; thu thập mẫu tiêu bản thực vật bậc cao theo quy chuẩn quốc tế, định danh qua bộ sách Cây cỏ Việt Nam (Phạm Hoàng Hộ, 1991, 1999) và Thực vật chí Việt Nam.
  • Điều tra cây tái sinh: Trong mỗi OTC $2.500\text{ m}^2$, bố trí 5 ô dạng bản (ODB) diện tích $4\text{ m}^2$ ($2\text{ m} \times 2\text{ m}$) tại 4 góc và tâm ô để đo đếm số lượng, đánh giá nguồn gốc (hạt/chồi), phẩm chất (A: tốt, B: trung bình, C: xấu) và phân cấp chiều cao ($H < 0,5\text{ m}$; $0,5\text{ m} \le H \le 1,0\text{ m}$; $1,0\text{ m} < H \le 1,5\text{ m}$; $H > 1,5\text{ m}$).
  • Tam giác hóa dữ liệu (Triangulation): Kết hợp kiểm chứng chéo giữa số liệu đo đếm thực địa, ảnh viễn thám Landsat/vệ tinh Google Earth, hồ sơ kiểm lâm địa phương và phỏng vấn cấu trúc xã hội học (KIP - Key Informant Panel).

Data và phân tích

Các công cụ toán học và phần mềm thống kê chuyên dụng (SPSS, Statistica, Excel Data Analysis Toolpak) được vận dụng để giải mã cấu trúc dữ liệu:

  • Chỉ số tương đồng Sorensen: $$S = \frac{2C}{A + B}$$ (Trong đó $A, B$ là số loài của hai quần xã so sánh; $C$ là số loài xuất hiện đồng thời ở cả hai khu hệ).
  • Chỉ số Giá trị Quan trọng (Importance Value Index - IVI%): $$IVI% = \frac{RD% + RF% + RDo%}{3}$$ (Với $RD%$ là Mật độ tương đối, $RF%$ là Tần số xuất hiện tương đối, $RDo%$ là Độ tàn che/Ưu thế tương đối).
  • Mô phỏng quy luật phân bố cấu trúc lâm phần ($N/D_{1.3}$ và $N/H_{vn}$): Kiểm định độ phù hợp của 3 hàm phân bố:
    • Hàm khoảng cách (Distance distribution)
    • Hàm phân bố giảm (Exponential decay)
    • Hàm phân bố Weibull 3 tham số: $$f(x) = \frac{c}{b} \left(\frac{x - a}{b}\right)^{c-1} e^{-\left(\frac{x - a}{b}\right)^c}$$
  • Mô hình hóa tương quan Allometric ($H_{vn}/D_{1.3}$): Kiểm tra 4 dạng hàm toán học:
    • Hàm Logarit: $H = a + b \ln(D)$
    • Hàm Power: $H = a D^b$
    • Hàm Compound: $H = a b^D$
    • Hàm S-curve: $H = e^{a + b/D}$
  • Lựa chọn mô hình tối ưu qua Tiêu chuẩn thông tin Akaike (AIC): $$AIC = 2k - 2\ln(L)$$ (Mô hình có giá trị AIC nhỏ nhất và hệ số xác định $R^2$ cao nhất được lựa chọn làm phương trình tương quan chuẩn).

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

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

Kết quả phân tích định lượng từ luận án đem lại 5 phát hiện mang tính đột phá cho sinh thái học Karst Việt Nam:

+----------------------------------------------------------------------------------------------------+
|                         TỔNG HỢP CÁC KẾT QUẢ ĐỘT PHÁ CỦA LUẬN ÁN                                   |
+----------------------------------------------------------------------------------------------------+
|   1. ĐA DẠNG KHU HỆ THỰC VẬT                                                                       |
|      • 608 loài thực vật có mạch, thuộc 388 chi, 131 họ (Ngành Magnoliophyta chiếm 95,72%).        |
|      • 12 loài đặc hữu hẹp (Cycas tropophylla, Chirita drakei, Schefflera alongensis...).          |
|      • 27 loài quý hiếm (Sách Đỏ Việt Nam / IUCN: Paphiopedilum concolor, Garcinia fagraeoides).   |
|                                                                                                    |
|   2. PHỔ DẠNG SỐNG RAUNKIAER                                                                       |
|      • SB = 78,12% Ph + 1,97% Ch + 2,14% H + 1,48% Cr + 1,81% Th + 14,48% Lp                       |
|      • Phanerophytes (Cây chồi trên) chiếm ưu thế áp đảo; Liana (Dây leo) đạt tỷ lệ cao đặc thù.  |
|                                                                                                    |
|   3. QUY LUẬT CẤU TRÚC ĐƯỜNG KÍNH & CHIỀU CAO (N/D1.3, N/Hvn)                                      |
|      • Thung lũng: Hàm Weibull mô phỏng chuẩn xác nhất cấu trúc N/D1.3 (dạng giảm lệch trái).       |
|      • Tương quan Allometric Hvn/D1.3: Hàm Power và Logarit đạt AIC tối ưu (R2 = 0,68 - 0,84).    |
|                                                                                                    |
|   4. ĐỘNG THÁI TÁI SINH TỰ NHIÊN CÂY GỖ                                                            |
|      • Nguồn gốc: Cây tái sinh từ hạt chiếm 76,4% - 82,1%; Tái sinh từ chồi chiếm 17,9% - 23,6%.  |
|      • Phẩm chất: Loại A (Tốt) chiếm 58,2% - 64,5%; Điểm nghẽn ở cấp chiều cao H > 1,5m.          |
|                                                                                                    |
|   5. LƯỢNG HÓA MỨC ĐỘ THOÁI HÓA                                                                    |
|      • 47,2% diện tích khảo sát ở mức thoái hóa trung bình - nặng do áp lực nổ mìn, than & đô thị.|
+----------------------------------------------------------------------------------------------------+
  1. Đặc trưng khu hệ thực vật phong phú và mức độ đặc hữu cao: Đã ghi nhận 608 loài thực vật có mạch thuộc 388 chi, 131 họ của 5 ngành thực vật bậc cao (Magnoliophyta chiếm ưu thế tuyệt đối với 582 loài, tương ứng 95,72%; Polypodiophyta có 18 loài; Lycopodiophyta có 4 loài; Psilotophyta có 2 loài; Gymnospermae có 2 loài). Đã phát hiện 12 loài đặc hữu hẹp của vùng Karst Hạ Long - Bái Tử Long như Chirita drakei, Chirita halongensis, Paraboea halongensis, Livistona halongensis, Cycas tropophylla, Schefflera alongensis, cùng 27 loài quý hiếm có tên trong Sách Đỏ Việt Nam (2007) và Danh lục Đỏ IUCN như Lan hài đốm (Paphiopedilum concolor - EN), Huyết giác (Dracaena cambodiana - VU), Trai lý (Garcinia fagraeoides - EN).
  2. Cấu trúc dạng sống phản ánh tính thích nghi khô hạn sinh lý: Phổ dạng sống thực vật chuẩn xác định theo công thức Raunkiaer là: $$SB = 78,12% Ph + 1,97% Ch + 2,14% H + 1,48% Cr + 1,81% Th + 14,48% Lp$$ Sự thống trị của nhóm chồi trên ($Ph = 78,12%$) kết hợp với tỷ lệ cao của dây leo thân gỗ ($Lp = 14,48%$) chứng minh thảm thực vật Karst Cẩm Phả là dạng rừng kín thường xanh mưa ẩm nhiệt đới chịu hạn trên đá xốp.
  3. Sự phân hóa không gian vi sinh cảnh sâu sắc:
    • Thung lũng: Cấu trúc 3-4 tầng rõ rệt, tầng vượt tán ($A_1$) đạt chiều cao $> 20\text{ m}$ gồm Sấu (Dracontomelon duperreanum), Phay (Duabanga grandiflora), Sâng (Pometia pinnata); tầng ưu thế sinh thái ($A_2$) gồm các loài họ Dẻ (Fagaceae), Long não (Lauraceae).
    • Chân núi: Cấu trúc 2-3 tầng, chiếm ưu thế bởi Ruối ô rô (Streblus ilicifolius), Mạy tèo (Streblus macrophyllus), Lát hoa (Chukrasia tabularis).
    • Sườn vách đá: Tán cây bị phá vỡ hoàn toàn, chỉ còn 1 tầng cây bụi và thảm bám vách đá gồm các chasmophyte điển hình (Euphorbia antiquorum, Dracaena cambodiana, Ficus alongensis).
  4. Quy luật cấu trúc toán học lâm phần: Mô hình hóa toán học chứng minh phân bố số cây theo cấp đường kính ($N/D_{1.3}$) ở thung lũng tuân theo hàm phân bố Weibull có dạng đường cong giảm (inverted-J shape), chỉ ra rằng lâm phần có mật độ cây gỗ nhỏ chiếm ưu thế và số cây giảm dần đều đặn khi đường kính tăng—dấu hiệu điển hình của rừng tự nhiên có khả năng tự điều chỉnh tốt. Phương trình tương quan $H_{vn}/D_{1.3}$ theo hàm Power ($H = a D^b$) và hàm Logarit đạt độ tin cậy cao nhất ($p < 0,001$, $R^2 > 0,72$, chỉ số AIC thấp nhất so với các dạng hàm khác).
  5. Động thái tái sinh và "nút thắt cổ chai" sinh thái: Mật độ cây tái sinh đạt mức trung bình khá ($2.800 - 4.500\text{ cây/ha}$). Cây tái sinh có nguồn gốc từ hạt chiếm áp đảo ($76,4% - 82,1%$), tỷ lệ cây có phẩm chất tốt (loại A) đạt trên $58%$. Tuy nhiên, có sự suy giảm nghiêm trọng về số lượng cá thể khi chuyển từ cấp chiều cao I ($H < 0,5\text{ m}$) lên cấp chiều cao IV ($H > 1,5\text{ m}$) (giảm từ $45,2%$ xuống chỉ còn $11,3%$), chứng minh sự thiếu hụt độ ẩm và dinh dưỡng khoáng trong kẽ đá là rào cản chọn lọc sinh thái lớn nhất đối với sự thành thục của cây con.

Implications đa chiều

  • Tiến bộ lý thuyết: Cung cấp khung lý thuyết hoàn chỉnh về cơ chế thích nghi và phân hóa cấu trúc lâm phần Karst ven biển, làm cơ sở bổ sung vào giáo trình Sinh thái học thực vật và Lâm học nhiệt đới.
  • Đổi mới phương pháp luận: Chuẩn hóa quy trình ứng dụng kết hợp giữa điều tra OTC diện tích lớn, tuyến vách đá hiểm trở và phương pháp chọn mô hình toán học qua tiêu chuẩn AIC, có thể nhân rộng cho các khu bảo tồn Karst dọc bờ biển miền Trung và vịnh Thái Lan.
  • Ứng dụng thực tiễn & Quản lý: Xác định danh mục 40 loài cây bản địa có khả năng tái sinh và bám đá tốt (như Streblus macrophyllus, Garcinia fagraeoides, Chukrasia tabularis, Ficus microcarpa) để phục vụ trực tiếp cho các dự án trồng rừng phục hồi cảnh quan mỏ than và bãi thải đá vôi sau khai thác tại Quảng Ninh.
  • Khuyến nghị chính sách: Đề xuất quy hoạch phân vùng bảo vệ nghiêm ngặt đối với 12 loài đặc hữu hẹp, thiết lập hành lang sinh thái kết nối các đảo đá vôi vịnh Bái Tử Long với di sản vịnh Hạ Long, chấm dứt việc cấp phép các mỏ đá vôi xi măng xâm lấn vùng đệm cảnh quan.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu khách quan:

  1. Rào cản tiếp cận không gian (Logistical Access Boundary): Do địa hình núi đá vôi trên biển Cẩm Phả có độ dốc thẳng đứng hiểm trở, phụ thuộc chặt chẽ vào con nước thủy triều và các rạn đá ngầm nguy hiểm, số lượng điểm khảo sát sườn vách đá bị giới hạn ở 34 vị trí và 18 tuyến điều tra—chưa bao phủ được toàn bộ các đảo biệt lập ngoài khơi xa.
  2. Hạn chế về dữ liệu chuỗi thời gian (Temporal Cross-sectional Constraint): Các số liệu cấu trúc và tái sinh được thu thập tại các thời điểm cắt ngang (cross-sectional snapshot), thiếu vắng các ô định vị sinh thái dài hạn (LTER - Long Term Ecological Research plots) để theo dõi động thái sinh trưởng từng cá thể qua nhiều thập kỷ.
  3. Thiếu vắng dữ liệu sinh thái phân tử và vi sinh vật: Nghiên cứu tập trung vào sinh thái quần xã hình thái, chưa phân tích cấu trúc di truyền quần thể của các loài đặc hữu hẹp (Cycas tropophylla, Chirita drakei) cũng như vai trò của hệ vi sinh vật nội cộng sinh (endophytic microbiome) trong rễ cây bám đá vôi.

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

  • Hướng 1 (LTER Karst Network): Thiết lập mạng lưới các ô tiêu chuẩn vĩnh cửu $1\text{ ha}$ theo tiêu chuẩn CTFS-ForestGEO tại vùng Karst vịnh Bái Tử Long để giám sát dòng carbon và động thái suy thoái lâm phần dưới tác động của biến đổi khí hậu.
  • Hướng 2 (Trait-based Ecology): Phân tích các đặc tính chức năng lá (Specific Leaf Area - SLA, thế nước mô lá, hàm lượng canxi mô bào) nhằm giải mã cơ chế sinh lý thích nghi của thực vật ưa vôi bắt buộc.
  • Hướng 3 (Conservation Genomics): Ứng dụng giải trình tự gen thế hệ mới (NGS) để đánh giá mức độ đa dạng di truyền và nguy cơ suy thoái cận huyết của các loài nguy cấp có kích thước quần thể nhỏ trên các đảo đá vôi.
  • Hướng 4 (Advanced Remote Sensing): Sử dụng thiết bị bay không người lái (UAV-LiDAR) kết hợp ảnh siêu phổ (Hyperspectral) để quét cấu trúc 3D tán rừng trên các vách đá vôi dựng đứng không thể tiếp cận bằng phương pháp thủ công.

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

Công trình của Hoàng Văn Hải (2020) tạo ra những ảnh hưởng sâu rộng trên nhiều bình diện:

  • Tác động học thuật (Academic Impact): Tạo lập bộ cơ sở dữ liệu thực vật học Karst ven biển chuẩn xác bậc nhất cho vùng Đông Bắc Bộ, được dự báo là nguồn trích dẫn nền tảng cho các nghiên cứu về bảo tồn đa dạng sinh học vịnh Bắc Bộ và các luận án chuyên ngành Thực vật học, Sinh thái học.
  • Chuyển đổi thực tiễn ngành Lâm nghiệp và Môi trường: Cung cấp cơ sở khoa học trực tiếp cho Tập đoàn Công nghiệp Than - Khoáng sản Việt Nam (TKV) và Tổng công ty Đông Bắc trong việc lựa chọn tập đoàn cây bản địa để hoàn nguyên môi trường các bãi thải mỏ than quanh khu vực Cẩm Phả.
  • Tác động chính sách công: Là tài liệu khoa học cốt lõi phục vụ Sở Nông nghiệp và Phát triển Nông thôn tỉnh Quảng Ninh, Ban Quản lý Vườn quốc gia Bái Tử Long và UBND thành phố Cẩm Phả trong việc điều chỉnh Quy hoạch bảo tồn đa dạng sinh học tỉnh Quảng Ninh tầm nhìn 2030, phục vụ việc duy trì tính toàn vẹn của kỳ quan thiên nhiên thế giới.
  • Giá trị kinh tế - xã hội: Bảo vệ thảm thực vật Karst chính là bảo vệ nguồn sinh thủy nước ngầm ngầm Karst cung cấp cho sinh hoạt dân sinh thành phố Cẩm Phả, đồng thời gìn giữ cảnh quan du lịch sinh thái biển đảo bền vững, giảm thiểu nguy cơ sạt lở đá và xói mòn ven biển.

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

Hệ thống kết quả nghiên cứu mang lại giá trị ứng dụng cụ thể cho từng nhóm đối tượng:

  • Nghiên cứu sinh và Giảng viên chuyên ngành Sinh thái / Lâm học: Tiếp cận bộ dữ liệu định lượng chuẩn xác về cấu trúc lâm phần, phương pháp mô hình hóa toán học ($N/D, N/H, H/D$) qua chỉ số AIC và kỹ thuật phân tích dạng sống Raunkiaer.
  • Các nhà thực vật học và phân loại học: Khai thác danh lục 608 loài với tọa độ phân bố, dữ liệu hình thái của 12 loài đặc hữu và 27 loài quý hiếm để phục vụ công tác định danh và lưu trữ nguồn gen.
  • Ban quản lý rừng đặc dụng và Vườn quốc gia: Ứng dụng 5 nhóm giải pháp quản trị, quy trình khoanh nuôi xúc tiến tái sinh tự nhiên và danh mục loài cây tái sinh mục tiêu cho hệ sinh thái núi đá vôi.
  • Các nhà hoạch định chính sách và Quản lý đô thị tỉnh Quảng Ninh: Sử dụng bản đồ phân vùng thoái hóa và ma trận phân tích SWOT để ban hành các quy chế quản lý xây dựng công trình, khai thác khoáng sản bền vững tại vùng đệm vịnh Bái Tử Long.

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

+----------------------------------------------------------------------------------------------------+
|                               GIẢI MÃ CÁC VẤN ĐỀ HỌC THUẬT CHUYÊN SÂU                               |
+----------------------------------------------------------------------------------------------------+
|                                                                                                    |
|  1. ĐÓNG GÓP LÝ THUYẾT ĐỘC ĐÁO NHẤT                                                                |
|     --> Mở rộng Lý thuyết Allometric & Diễn thế sinh thái vi sinh cảnh Karst.                      |
|                                                                                                    |
|  2. ĐỔI MỚI PHƯƠNG PHÁP NGHIÊN CỨU                                                                 |
|     --> Kết hợp OTC diện tích lớn (2500m2) + Tuyến vách đá + Lựa chọn mô hình phi tuyến qua AIC.   |
|                                                                                                    |
|  3. PHÁT HIỆN GÂY BẤT NGỜ NHẤT TỪ DỮ LIỆU                                                          |
|     --> Mật độ cây tái sinh gốc hạt cao (76-82%) nhưng nghẽn chuyển tiếp sinh trưởng H > 1,5m.     |
|                                                                                                    |
|  4. QUY TRÌNH TÁI LẬP THỰC NGHIỆM (REPLICATION PROTOCOL)                                           |
|     --> Hệ thống tọa độ chuẩn, thiết kế ODB chuẩn mực, công thức toán học tường minh.              |
|                                                                                                    |
|  5. CHƯƠNG TRÌNH NGHIÊN CỨU 10 NĂM TIẾP THEO                                                       |
|     --> Tích hợp UAV-LiDAR 3D + Giải trình tự gen NGS + Ô vĩnh cửu LTER 1ha.                       |
|                                                                                                    |
+----------------------------------------------------------------------------------------------------+

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng Thuyết Phân vùng Sinh thái Thảm thực vật Karst thông qua lượng hóa toán học mối quan hệ giữa vi địa hình (thung lũng - chân núi - sườn vách) với cấu trúc allometric ($H_{vn}/D_{1.3}$) và động thái tái sinh. Luận án đã bác bỏ giả định cho rằng rừng trên đá vôi có cấu trúc ngẫu nhiên vô trật tự, chứng minh rằng cấu trúc lâm phần Karst vẫn tuân theo các quy luật phân bố xác suất nghiêm ngặt (hàm Weibull và hàm lũy thừa) nhưng bị nén biên độ sinh trưởng bởi giới hạn đất và nước.

2. Phương pháp nghiên cứu có gì đổi mới so với các công trình trước đây?

So với nghiên cứu của Nguyễn Tiến Hiệp (2005) (chỉ điều tra mô tả tuyến định tính) hay Trần Hữu Viên (2004) (chỉ khảo sát tổng quát trên diện rộng), luận án đã:

  • Chuẩn hóa kích thước ô tiêu chuẩn lên $2.500\text{ m}^2$ ($50\text{ m} \times 50\text{ m}$) phù hợp với độ phân tán cao của cây gỗ Karst.
  • Thiết kế hệ thống 18 tuyến điều tra chuyên biệt trên sườn vách đá dốc đứng.
  • Ứng dụng Tiêu chuẩn thông tin Akaike (AIC) để thẩm định tính ưu việt của các mô hình toán học sinh thái thay vì chỉ dựa vào hệ số $R^2$ đơn thuần.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất từ dữ liệu?

Phát hiện bất ngờ nhất là tỷ lệ tái sinh tự nhiên từ hạt đạt rất cao ($76,4% - 82,1%$) ngay cả ở những khu vực rừng bị tác động mạnh (trạng thái IIa, IIb). Điều này đi ngược lại định kiến sinh thái học trước đây cho rằng trên núi đá vôi khô hạn, cây rừng chủ yếu tái sinh chồi vô tính từ gốc cây bị chặt. Tuy nhiên, dữ liệu cũng chỉ ra nghịch lý sinh thái: tỷ lệ sống sót để vượt qua cấp chiều cao $1,5\text{ m}$ lại cực kỳ thấp do hiện tượng "chết khát sinh lý" vào mùa khô khi rễ cây chưa kịp đâm sâu vào các kẽ nứt có tích tụ mùn.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Toàn bộ quy trình được mô tả chi tiết tại Chương 2 và Phụ lục 2:

  • Tọa độ địa lý và sơ đồ bố trí 50 OTC và 18 tuyến điều tra tại 34 điểm.
  • Quy cách kích thước OTC ($50\text{ m} \times 50\text{ m}$) và bố trí 5 ô dạng bản ($2\text{ m} \times 2\text{ m}$).
  • Bảng tính mẫu $f_i$, $X_i$ và các tham số khởi tạo cho các phương trình hồi quy phi tuyến trong môi trường SPSS/Excel.

5. Chương trình nghiên cứu 10 năm tiếp theo (10-Year Research Agenda) được phác thảo như thế nào?

Chương trình 10 năm tập trung vào 3 trọng tâm:

  1. Năm 1-3: Số hóa toàn bộ dữ liệu lâm phần bằng công nghệ viễn thám LiDAR và thiết lập 5 ô định vị LTER vĩnh cửu.
  2. Năm 4-7: Triển khai nghiên cứu Functional Genomics trên 12 loài đặc hữu hẹp để xác định các gen thích nghi chịu canxi cao và chịu hạn.
  3. Năm 8-10: Hoàn thiện mô hình kinh tế sinh thái chi trả dịch vụ môi trường rừng Karst (PES) và quy chuẩn kỹ thuật phục hồi rừng Karst quốc gia.

Kết luận

Luận án tiến sĩ của Hoàng Văn Hải (2020) là một công trình khoa học công phu, chuẩn mực và có giá trị học thuật cao trong lĩnh vực Sinh thái học thực vật Karst nhiệt đới. Những kết luận cốt lõi bao gồm:

  1. Xác lập hệ thống dẫn liệu khu hệ thực vật toàn diện: Ghi nhận 608 loài thực vật có mạch, 12 loài đặc hữu hẹp và 27 loài quý hiếm tại vùng Karst thành phố Cẩm Phả.
  2. Định lượng phổ dạng sống thích nghi: Chứng minh tính ưu thế tuyệt đối của nhóm chồi trên ($Ph = 78,12%$) và tỷ lệ cao của dây leo thân gỗ ($Lp = 14,48%$), phản ánh cấu trúc rừng kín thường xanh mưa ẩm nhiệt đới chịu hạn.
  3. Mô hình hóa thành công cấu trúc lâm phần: Xác định hàm Weibull và các phương trình Allometric Power/Logarit (thẩm định qua chỉ số AIC) là các công cụ toán học tối ưu mô tả quy luật phân bố $N/D_{1.3}$ và tương quan $H_{vn}/D_{1.3}$.
  4. Giải mã động thái và điểm nghẽn tái sinh: Phát hiện ưu thế của cây tái sinh gốc hạt ($>76%$) cùng sự suy giảm mật độ nghiêm trọng ở cấp chiều cao $> 1,5\text{ m}$ do hạn chế lập địa.
  5. Hoạch định 5 nhóm giải pháp bảo tồn đồng bộ: Kết hợp giữa nâng cao nhận thức, tăng cường năng lực tuần tra kiểm lâm, hoàn thiện cơ chế chính sách vùng đệm di sản, phát triển kinh tế sinh thái và áp dụng các giải pháp lâm sinh phục hồi rừng bản địa.

Công trình không chỉ góp phần thúc đẩy sự phát triển của sinh thái học quần xã thực vật tại Việt Nam mà còn là tài liệu tham khảo thiết yếu cho các chiến lược bảo tồn đa dạng sinh học và phát triển bền vững không gian kinh tế - sinh thái vịnh Bái Tử Long và vịnh Hạ Long.