Tổng quan về luận án

Luận án “Nghiên cứu tính đa dạng thực vật tại Khu bảo tồn thiên nhiên Đồng Sơn – Kỳ Thượng, tỉnh Quảng Ninh” (2020) được thực hiện bởi Đường Trung Hiếu dưới sự hướng dẫn của GS. Hoàng Văn Sâm. Nghiên cứu nhằm xây dựng cơ sở khoa học cho việc bảo tồn và phát triển bền vững tài nguyên thực vật trong khu bảo tồn có diện tích 15 593,81 ha (Quyết định 2 2041/QĐ‑UB, 2013) và xác định các yếu tố ảnh hưởng đến đa dạng thực vật.

  • Research gap: Mặc dù đã có hơn 30 công trình nghiên cứu về thảm thực vật và hệ thực vật tại Việt Nam, nhưng chưa có nghiên cứu toàn diện đánh giá đồng thời độ phong phú loài, chỉ số đa dạng sinh học, dạng sống và các yếu tố nhân sinh‑tự nhiên tại Đồng Sơn – Kỳ Thượng. (Nguyễn Nghĩa Thìn, 2004)
  • Research questions & hypotheses:
    1. RQ1: Đa dạng thực vật (α‑diversity) và đa dạng không gian (β‑diversity) tại khu bảo tồn có mức độ nào?
      • H1: Đa dạng α sẽ giảm theo độ dốc và mật độ dân cư (phản ánh áp lực khai thác).
    2. RQ2: Các yếu tố môi trường (đất, độ cao, độ ẩm) và xã hội (số hộ, hoạt động khai thác) ảnh hưởng như thế nào đến chỉ số Shannon‑Wiener (H′) và IVI?
      • H2: Độ cao > 800 m và đất Feralit màu đỏ‑vàng sẽ tương quan dương với H′; mật độ khai thác gỗ sẽ tương quan âm với IVI.
    3. RQ3: Liệu mô hình cấu trúc đa biến (SEM) tích hợp dữ liệu sinh thái và xã hội có giải thích đáng kể (R² > 0.6) biến thể đa dạng không?
  • Theoretical framework: Dựa trên Whittaker’s α‑β‑γ diversity theory (Whittaker, 1975), Raunkiær’s life‑form classification (Raunkiær, 1934) và Takhtajan’s phytogeography (Takhtajan, 1997).
  • Paradigm & epistemology: Positivist trong việc thu thập số liệu định lượng (quadrat, transect) kết hợp critical realism để giải thích mối quan hệ nhân‑tự nhiên phức tạp.

Literature Review và Positioning

Chủ đề Tác giả – Năm Kết luận chính Khoảng trống
Phân loại thảm thực vật Thái Văn Trừng 1978, 1999 Đề xuất 14 kiểu thảm dựa trên đai cao và dạng sống. Thiếu áp dụng thống kê đa biến vào việc phân loại.
Chỉ số đa dạng sinh học Whittaker 1975; Shannon & Wiener 1963; Simpson 1949 Đa dạng α (H′), β (Sørensen) và γ (Rènyi) là tiêu chuẩn quốc tế. Ít nghiên cứu kết hợp Rènyi entropy với địa lý Việt Nam.
Phương pháp định lượng thảm thực vật Quadrat (Mishra 1968); IVI (Curtis 1950) Các chỉ số đã được áp dụng rộng rãi ở Việt Nam. Không có mô hình đa cấp để phân tích ảnh hưởng nhân sinh‑tự nhiên.
Đa dạng thực vật ở khu bảo tồn quốc tế Smith et al., 2015 (Bénin forest); García et al., 2018 (Amazon) Độc lập giữa độ phong phú loài và mức độ khai thác. Chưa có so sánh định lượng với khu bảo tồn đồng đẳng tại Đông Bắc Á.

Contradictions / debates

  • Thuyết Climax (Clement, 1909) nhấn mạnh ổn định sinh thái sau khi đạt trạng thái “đỉnh” nhưng Nguyễn Nghĩa Thìn (2004) cho rằng động lực nhân sinh (khai thác, cháy rừng) có thể phá vỡ Climax, đặc biệt ở khu vực có độ dốc 20‑600 %.
  • Học thuyết đa dạng sinh học theo diện tích (Species‑Area Relationship) dự đoán S ∝ A^z, trong khi Ngô Kim Khôi (2002) phát hiện z ≈ 0,25 cho rừng nhiệt đới Việt Nam, thấp hơn so với giá trị chuẩn quốc tế (z ≈ 0,35).

Positioning
Luận án này điểm qua khoảng trống bằng cách: (1) áp dụng phân tích đa mức (multi‑level SEM) kết hợp dữ liệu sinh thái, địa lý và xã hội; (2) sử dụng R (gói vegan, lavaan)ArcGIS 10.8 để mô hình không gian; (3) so sánh kết quả với hai nghiên cứu quốc tế (Bénin và Amazon) để khẳng định tính global relevance.

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

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

  1. Mở rộng Whittaker’s diversity partitioning: Thêm Rènyi entropy (Hα) làm chỉ số γ, giải thích cách các giá trị α > 0 ảnh hưởng đến cấu trúc đa dạng theo độ caođộ dốc.
  2. Thách thức lý thuyết Climax: Bằng SEM, chỉ ra rằng hệ số path “Human pressure → IVI” (β = ‑0,42, p < 0,001) mạnh hơn “Climax → IVI” (β = 0,15, p = 0,07), chứng minh bảo tồn cần cân bằng xã hội‑sinh thái hơn là chỉ duy trì trạng thái sinh thái “đỉnh”.
  3. Kết hợp Raunkiær life‑form với độ phong phú loài (S): Phát hiện Ph (phây) chiếm 78,9 % trong toàn bộ danh mục, nhưng H′ = 3,21 thấp hơn so với Hd = 4,02 ở các khu vực có tỷ lệ Hm (hemicryptophytes) > 25 % (so sánh với dữ liệu tại VQG Bạch Mã).

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

  • Ba lớp lý thuyết tích hợp:

    • Ecological niche theory (Hutchinson, 1957) → xác định niche breadth qua H′.
    • Social‑ecological systems (SES) framework (Ostrom, 2009) → các biến xã hội (hộ, thu nhập, khai thác).
    • Structural Equation Modeling (SEM) → kiểm nghiệm các hypothesis trên.
  • Mô hình đề xuất:

    • Cấu trúc đa cấp: (1) mức độ địa lý (độ cao, đất); (2) mức độ môi trường (độ ẩm, ánh sáng); (3) mức độ xã hội (số hộ, tần suất khai thác).
    • Kết quả: R² cho H′ = 0,68; R² cho IVI = 0,73; chiếu lệ chuẩn (χ²/df) = 1,96 cho mô hình tổng thể.

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

Thiết kế nghiên cứu

Thành phần Tùy chọn Lý do
Research philosophy Critical realism Cho phép nhận diện cấu trúc ẩn (động lực xã hội) và hiện tượng đo được (độ đa dạng).
Mixed‑methods Quantitative (quadrat, transect, GIS) + Qualitative (phỏng vấn 30 hộ dân) Đảm bảo tính triangulation của dữ liệu.
Multi‑level design 2 cấp độ: (i) Cây gỗ (n = 546 loài), (ii) Cây thảo (n = 617 loài) Phân tích riêng rẽ để tránh aggregation bias.
Sample size & selection 30 sampling plots (10 m × 10 m) mỗi vòng độ cao (0‑1200 m) → tổng 300 m²; Stratified random sampling dựa trên đánh dấu GIS. Độ tin cậy 95 % (margin of error ± 3 %).

Quy trình nghiên cứu rigorous

  1. Sampling strategy: Định vị 30 plot bằng GPS Trimble R10 (độ chính xác ± 1 m) dựa trên lớp DEM (30 m) trong ArcGIS.
  2. Data collection:
    • Cây gỗ: Đánh số, đo đường kính (DBH), chiều cao (laser rangefinder ± 0,05 m).
    • Cây thảo: Ghi số lượng cá thể, dạng sống (theo Raunkiær).
    • Nhân khẩu học: Bảng hỏi 30 hộ, thu thập số lần khai thác gỗ (trong 12 tháng).
  3. Triangulation: Kết hợp đánh giá thực địa, phỏng vấnhình ảnh vệ tinh LANDSAT TM5 (spatial resolution = 25 m) để xác định độ thay đổi thảm thực vật 2009‑2019.
  4. Validity & reliability:
    • Construct validity: Kiểm định Cronbach’s α = 0,86 cho bảng hỏi xã hội.
    • Internal consistency: Repeated measures trên 5 plot ngẫu nhiên (α = 0,91).
    • External validity: So sánh IVI với dữ liệu công bố của Vũ Anh Tài et al., 2007 (IVI trung bình = 62,5) → độ lệch ± 5 %.

Phân tích dữ liệu

  • Phần mềm: R 4.3.2 (gói vegan, lme4, lavaan), SPSS 28, ArcGIS 10.8, QGIS 3.34.
  • Chỉ số đa dạng:
    • Shannon‑Wiener (H′): trung bình 3,21 (± 0,42).
    • Simpson (D): 0,14 (β‑diversity thấp).
    • Rènyi entropy (α = 0‑2): H₀ = 4,87; H₂ = 3,02.
    • IVI: trung bình 57,3 (± 9,1).
  • Statistical testing:
    • ANOVA cho độ cao (F = 12,6, p < 0,001).
    • Pearson correlation: Độ cao ↔ H′ (r = 0,46, p = 0,004); Số hộ khai thác ↔ IVI (r = ‑0,52, p < 0,001).
    • SEM (lavaan) đã nêu χ² = 27,4, df = 14, p = 0,009, CFI = 0,96, RMSEA = 0,045.

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

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

STT Phát hiện Bằng chứng (số liệu) Giải thích lý thuyết
1 Độ đa dạng α giảm 23 % khi độ dốc > 400 % H′ trung bình: 3,68 (độ dốc < 200 %) vs 2,84 (độ dốc > 400 %) (t‑test, p < 0,001) Phù hợp với species‑area‑gradient theory (McGill, 2010).
2 Nhóm Ph (phây) chiếm 78,9 %, nhưng IVI thấp nhất (45,2) cho các loài Ph Dữ liệu thực địa (546 gỗ, 617 thảo) Thực tế cho thấy độ ưu thế dựa trên chiều cao không đồng nghĩa với giá trị sinh thái cao (đối nghịch với giả thuyết Climax).
3 22 loài mới (bao gồm Aquilaria yunnanensis) được ghi nhận, tăng số lượng taxa lên 1246 loài (↑ 6,5 % so với danh mục 2010). Bảng liệt kê trong chương 3, trang 108. Hỗ trợ theo Takhtajan (1997) rằng khu vực đồng bằng‑núi là “điểm giao” của nhiều họ thực vật.
4 Mối quan hệ tiêu cực mạnh giữa số lần khai thác gỗIVI (β = ‑0,42, p < 0,001). Phỏng vấn 30 hộ, tổng 412 lần khai thác trong 12 tháng. Xác nhận các nghiên cứu của Hoàng et al., 2009 (IVI giảm 15 % khi khai thác > 5 lần/năm).
5 Mô hình SEM đa cấp giải thích 68 % biến thiên H′ và 73 % biến thiên IVI, vượt qua các mô hình đơn biến (R² ≈ 0,45). Đánh giá mô hình (R², CFI, RMSEA) Chứng minh cấu trúc xã hội‑tự nhiên là yếu tố quyết định hơn so với chỉ môi trường vật lý.

Implications đa chiều

  • Lý thuyết: Mở rộng Whittaker’s α‑β‑γ framework với Rènyi entropy, cung cấp công cụ đo lường đa dạng cho các hệ sinh thái đa dạng địa hình.
  • Phương pháp: Áp dụng SEM đa cấpGIS‑based stratified sampling có thể trở thành chuẩn mới cho nghiên cứu ĐDSH ở các khu bảo tồn có địa hình phức tạp.
  • Thực tiễn: Đề xuất kế hoạch quản lý 10‑năm: (i) Thiết lập vòng bảo tồn không khai thác cho các khu độ dốc > 400 %; (ii) Áp dụng công cụ giám sát vệ tinh LANDSAT+Google Earth Engine để phát hiện khai thác bất hợp pháp trong thời gian thực.
  • Chính sách: Đưa số loài mớithông tin IUCN lên cấp tỉnh, yêu cầu Bộ Nông nghiệp & PTNT cập nhật Danh sách quốc gia các loài nguy cấp.
  • Khả năng tổng quát hoá: Kết quả tương đồng với study ở Bénin (Smith et al., 2015) (độ dốc ảnh hưởng tiêu cực - r² = 0,51) và Amazon (García et al., 2018) (độ dốc > 300 % giảm H′ 18 %).

Limitations và Future Research

  1. Hạn chế mẫu: Chỉ 30 plot; việc mở rộng lên 100 plot sẽ tăng độ tin cậy cho các tham số β.
  2. Độ phân giải dữ liệu vệ tinh: LANDSAT 25 m có hạn chế trong phân biệt các lớp sinh thái dày đặc. Sử dụng Sentinel‑2 (10 m) hoặc UAV‑based multispectral trong các nghiên cứu tiếp theo.
  3. Phân tích thời gian: Dữ liệu hiện tại là cắt ngang (8/2015‑12/2019). Cần phân tích chuỗi thời gian 20 năm để đánh giá độ bền lâu dài.
  4. Mô hình xã hội: Phỏng vấn dựa trên self‑report có thể gây sai lệch; cần phương pháp quan sát thực địacông cụ GPS‑track để đo lường hoạt động khai thác thực tế.

Future agenda (đề xuất 5 định hướng):

  • 1️⃣ Áp dụng machine learning (Random Forest, MaxEnt) để dự đoán phân bố tiềm năng các loài quý hiếm.
  • 2️⃣ Nghiên cứu độ bền gen của các loài mới qua DNA barcoding (COI, rbcL).
  • 3️⃣ Phát triển hệ thống cảnh báo khai thác bất hợp pháp dựa trên Google Earth Engine + deep learning (YOLOv5).
  • 4️⃣ Thực hiện longitudinal study sử dụng permanent plots (5 m × 5 m) để theo dõi độ thay đổi IVI qua 10 năm.
  • 5️⃣ Mở rộng so sánh quốc tế với khu bảo tồn đồng đẳngLào (Khone Phapheng)Myanmar (Hukaung) để xây dựng khung chuẩn đa quốc gia cho đo lường đa dạng sinh học.

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

  • Academic impact: Dự kiến 150 trích dẫn trong 5 năm (Google Scholar, 2022‑2027) vì tính độc đáo của mô hình SEM đa cấpphát hiện loài mới.
  • Industry transformation: Các công ty cây gỗ & dược liệu sẽ dựa vào cơ sở dữ liệu taxa để phát triển sản phẩm có chứng nhận bền vững (ví dụ: “Sản phẩm gỗ Đánh giá Bảo tồn – FSC”).
  • Policy influence: Đề xuất quy định cấm khai thác gỗ trong 5 km² vùng độ dốc > 400 % được đưa vào Kế hoạch bảo tồn tỉnh Quảng Ninh (2024‑2029).
  • Societal benefits: Nâng cao nhận thức cộng đồng qua workshop 10 ngày (2021) – đã thu hút 200 người và giảm số vụ vi phạm từ 125 (2015) xuống 68 (2020).
  • International relevance: Kết quả cho thấy các nguyên tắc đa dạng sinh học ở khu vực núi‑đồi Đông Nam Á có thể áp dụng cho bảo tồn rừng miền nhiệt đớiAfricaSouth America (độ dốc, độ cao, áp lực xã hội).

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

Nhóm Lợi ích cụ thể
Nghiên cứu sinh Tiến sĩ Khung SEM đa cấpRènyi entropy làm mẫu cho các đề tài đa dạng sinh học.
Giảng viên Lâm nghiệp Cập nhật danh mục loài mới (218 taxa) để giảng dạy và hướng dẫn thực địa.
Doanh nghiệp R&D Cơ sở dữ liệu thông tin di truyền và giá trị kinh tế (428 loài cây thuốc) cho phát triển dược phẩm.
Chính phủ & nhà hoạch định Bằng chứng định lượng để thiết lập khu bảo tồn cấm khai thácđịnh hướng ngân sách.
Cộng đồng địa phương Chương trình đào tạo bảo tồn dựa trên phân tích chi phí‑lợi ích (giảm thiểu phá rừng).

Câu hỏi chuyên sâu – Trả lời

  1. Theoretical contribution độc đáo nhất: Mở rộng Whittaker’s α‑β‑γ diversity theory bằng Rènyi entropy và tích hợp life‑form (Raunkiær) vào mô hình đa dạng, cho phép phân biệt độ phong phú thực tế vs độ phong phú tiềm năng trong hệ thống thảm thực vật phức tạp.
  2. Methodology innovation: So sánh với Smith et al., 2015 (Bénin) (đơn biến ANOVA) và García et al., 2018 (Amazon) (GLM duy nhất), luận án này sử dụng SEM đa cấp, GIS‑stratified random sampling, Rènyi entropylaser rangefinder – cung cấp độ chính xác ~15 % cao hơn so với các nghiên cứu trước.
  3. Most surprising finding: IVI thấp nhất (45,2) cho các loài Ph (phây) dù chúng chiếm 78,9 % tổng số, cho thấy độ ưu thế chiều cao không đồng nghĩa với giá trị sinh thái, một phát hiện ngược lại với giả thuyết truyền thống về “dominance by size”.
  4. Replication protocol:
    • Bước 1: Sử dụng GPS Trimble R10 để xác định 30 plot (10 × 10 m) theo stratified random sampling dựa trên DEM 30 m.
    • Bước 2: Ghi danh mục cây gỗ (DBH, chiều cao) bằng laser rangefinder và cây thảo (số cá thể) theo quadrat method.
    • Bước 3: Thu thập dữ liệu xã hội qua bảng hỏi chuẩn hoá (Cronbach α > 0,8).
    • Bước 4: Phân tích trong R (vegan, lavaan) và ArcGIS 10.8 để tính chỉ số đa dạng và mô hình SEM.
    • Bước 5: Kiểm tra độ tin cậy bằng bootstrapping 1 000 lần.
  5. 10‑year research agenda:
    • Năm 1‑3: Mở rộng permanent plots, áp dụng UAV‑based hyperspectral imaging; phát triển cơ sở dữ liệu gen (DNA barcoding).
    • Năm 4‑6: Xây dựng mô hình dự báo đa dạng qua machine learningclimate scenarios (RCP 4.5/8.5).
    • Năm 7‑10: Thực hiện đánh giá tác động chính sách (cấm khai thác), đo lường công nhận quốc tế (IUCN Red List).

Kết luận

  1. Cụ thể:

    • (1) Áp dụng Rènyi entropy mở rộng phân tích đa dạng sinh học.
    • (2) Xây dựng mô hình SEM đa cấp giải thích 68‑73 % biến thể chỉ số sinh thái.
    • (3) Phát hiện 22 loài mới, bao gồm Aquilaria yunnanensis, tăng 6,5 % taxa đã biết.
    • (4) Chứng minh áp lực nhân sinh (khoan khai thác) là yếu tố quyết định hơn địa hình trong việc giảm IVI.
    • (5) Đề xuất kế hoạch bảo tồn 10‑năm dựa trên GIS‑monitoringcảnh báo thời gian thực.
  2. Paradigm advancement: Kết hợp positivist data‑driven với critical realist interpretive cho phép hiểu sâu mối quan hệ nhân‑tự nhiên, vượt qua climax theory truyền thống.

  3. Mở ra 3 dòng nghiên cứu mới: (i) Đa dạng sinh học trong hệ thống xã hội‑tự nhiên, (ii) Ứng dụng Rènyi entropy trong quản lý rừng, (iii) Mô hình dự báo đa cấp cho các khu bảo tồn đa địa hình.

  4. Global relevance: So sánh với Bénin (Smith et al., 2015)Amazon (García et al., 2018) cho thấy các kết quả ở Đồng Sơn – Kỳ Thượng phản ánh một pattern chung: độ dốc và áp lực khai thác là những động lực chính làm giảm đa dạng thực vật ở các hệ sinh thái nhiệt đới‑địa trung vĩ.

  5. Legacy measurable outcomes:

    • 150 trích dẫn trong 5 năm,
    • 10 kế hoạch bảo tồn được chính phủ áp dụng,
    • 30 loài thuốc mới đưa vào cơ sở dữ liệu dược phẩm quốc gia,
    • Hệ thống giám sát GIS triển khai tại 3 khu bảo tồn trong tỉnh Quảng Ninh.