Tổng quan luận án

Luận án tiến sĩ lâm nghiệp với đề tài "Đặc điểm lâm học và sinh thái quần thể loài Thông 5 lá (Pinus dalatensis Ferré) ở Tây Nguyên" do nghiên cứu sinh Lê Cảnh Nam thực hiện tại Viện Khoa học Lâm nghiệp Việt Nam (thuộc khóa đào tạo 25/2013), dưới sự hướng dẫn khoa học của GS. Bảo Huy (người hướng dẫn 1) và TS. Nguyễn Thành Mến (người hướng dẫn 2). Công trình thuộc chuyên ngành Lâm sinh, mã số 9 62 02 05.

Tính cấp thiết và khoảng trống nghiên cứu được xác định xuất phát từ vị trí phân loại và tình trạng bảo tồn của đối tượng nghiên cứu. Cây Thông 5 lá (Pinus dalatensis Ferré), được nhà thực vật học người Pháp Ferré mô tả và công bố năm 1960, là loài thực vật đặc hữu theo nghĩa rộng của dãy Trường Sơn, đồng thời là loài cây gỗ quý hiếm thuộc nhóm IIA theo Nghị định số 06/2019/NĐ-CP của Chính phủ và xếp vào nhóm Sắp nguy cấp (NT - Near Threatened) trong Danh lục Đỏ của Liên minh Bảo tồn Đa dạng Sinh học Quốc tế (IUCN) năm 2019. Đây là loài cây gỗ lớn thường xanh, có thể đạt chiều cao 30–40 m và đường kính ngang ngực lên đến 2,5 m, phân bố trong kiểu rừng hỗn giao cây lá rộng lá kim ở các cao nguyên Langbiang, Ban Mê Thuột, Pleiku và Kon Tum. Trên phạm vi toàn cầu cũng như tại Việt Nam, Thông 5 lá hiện chỉ còn ghi nhận dưới 10 địa điểm phân bố với số lượng thường dưới 100 cây trưởng thành trong mỗi vùng do sự suy giảm môi trường sống. Các nghiên cứu trước đây chủ yếu tập trung vào phân loại thực vật, mô tả hình thái, nhân giống hoặc ghi nhận vùng phân bố mới một cách riêng lẻ. Nguồn dữ liệu định lượng có hệ thống về cấu trúc lâm phần, đặc tính tái sinh, quy luật sinh trưởng gắn với biến động khí hậu và yêu cầu sinh thái của quần thể Thông 5 lá phục vụ công tác quy hoạch bảo tồn còn thiếu hụt.

Mục tiêu nghiên cứu của luận án được xác lập cụ thể như sau:

  1. Mục tiêu chung: Góp phần cung cấp cơ sở, thông tin dữ liệu khoa học về lĩnh vực lâm học và sinh thái rừng phục vụ bảo tồn và phát triển bền vững loài và quần thể Thông 5 lá quý hiếm ở Tây Nguyên.
  2. Mục tiêu cụ thể:
    • Mục tiêu cụ thể 1: Định lượng được các đặc điểm lâm học chính của cá thể và quần thể Thông 5 lá bao gồm mô phỏng cấu trúc, tái sinh, sinh trưởng và tăng trưởng dưới ảnh hưởng của nhân tố khí hậu và môi trường phục vụ việc áp dụng các biện pháp lâm sinh và bảo tồn.
    • Mục tiêu cụ thể 2: Xác định được các nhân tố sinh thái chủ yếu ảnh hưởng đến mật độ và lập bản đồ mật độ quần thể loài Thông 5 lá phục vụ cho quản lý, xử lý lâm sinh, bảo tồn và phát triển quần thể.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Loài Thông 5 lá (Pinus dalatensis Ferré), thuộc chi Pinus, họ Thông (Pinaceae), bộ Thông (Pinales), lớp Thông (Pinopsida), ngành Hạt trần (Pinophyta).
  • Phạm vi không gian: Các lâm phần có phân bố loài Thông 5 lá thuộc lâm phận quản lý của 3 Vườn quốc gia (VQG) tại khu vực Tây Nguyên, bao gồm VQG Bidoup - Núi Bà (tỉnh Lâm Đồng), VQG Chư Yang Sin (tỉnh Đắk Lắk) và VQG Kon Ka Kinh (tỉnh Gia Lai).
  • Nguồn kế thừa: Luận án có kế thừa một phần số liệu điều tra sinh thái loài tại VQG Chư Yang Sin từ đề tài khoa học công nghệ cấp tỉnh Đắk Lắk "Điều tra phân bố, sinh thái một số loài thực vật thân gỗ quí hiếm phục vụ công tác bảo tồn nguồn gen tại tỉnh Đăk Lăk".

Tổng quan tài liệu và vị trí của luận án

Tổng quan tài liệu trong luận án phân tích các công trình nghiên cứu trên thế giới và tại Việt Nam qua các nhánh chuyên môn:

Về cấu trúc quần thể thực vật rừng và mô hình hóa cấu trúc:

  • Trên thế giới, các nghiên cứu nền tảng về cấu trúc rừng mưa nhiệt đới bắt đầu từ mô tả định tính của Richards (1952), Baur (1976), Odum (1971), sau đó chuyển dần sang lượng hóa bằng mô hình toán thống kê (Zheng và Zhou, 2010). Chỉ số quan trọng loài (IV%) được Curtis và McIntosh (1950) thiết lập dựa trên tổ hợp ba đại lượng tỉ lệ mật độ (N%), tỉ lệ tiết diện ngang (G%) và tần suất xuất hiện (F%). Mô hình hóa phân bố số cây theo cấp đường kính (N/D) được Rollet (1971), Baker và cs, Barbosa và cs áp dụng qua các hàm Weibull, Mayer, Hyperbol, Poisson, Logarit chuẩn và họ Pearson. Đối với cấu trúc tầng rừng (N/H), trắc diện đứng của Davis và Richards (1933–1943) được thay thế dần bằng các phân bố xác suất như Weibull, Poisson. Cấu trúc mặt bằng được lượng hóa thông qua phân bố khoảng cách cá thể và tiêu chuẩn kiểm định $U$ của Clark và Evans (1954).
  • Tại Việt Nam, Trần Ngũ Phương (1970) và Thái Văn Trừng (1978, 1999) đặt nền móng phân loại thảm thực vật và mô tả 5 tầng rừng nhiệt đới. Việc toán học hóa cấu trúc lâm phần được phát triển bởi Đồng Sĩ Hiền (1974), Nguyễn Văn Trương (1973, 1983), Trần Văn Con (1991, 2001, 2011), Bảo Huy (1997, 2017a), Nguyễn Hải Tuất và cs (2006). Nhiều nghiên cứu cụ thể đã kiểm chứng tính phù hợp của hàm phân bố Weibull và phân bố khoảng cách đối với các trạng thái rừng tự nhiên (Nguyễn Thành Mến, 2004; Đỗ Thị Hà và Bùi Thị Hằng, 2010; Lê Cảnh Nam và cs, 2012, 2016; Bùi Mạnh Hưng và Lê Xuân Trường, 2017).

Về sinh thái rừng và quan hệ sinh thái loài:

  • Phương pháp phân tích quan hệ liên loài được Stephen và Garray (1986), Nguyễn Hải Tuất (1990) và Bảo Huy (1997, 2017a) chuẩn hóa để xác định tương tác hỗ trợ, cạnh tranh hay độc lập giữa các loài cây ưu thế. Vấn đề mô hình hóa mối quan hệ giữa biến số thực vật và tổ hợp nhân tố môi trường ($y_i = f(x_j) + \varepsilon$) được Bảo Huy (2015a, 2017a) phát triển qua phương pháp mã hóa biến theo gradient để nhận diện các nhân tố sinh thái chi phối.

Về vòng năm (Dendrochronology) và sinh trưởng cây rừng theo khí hậu:

  • Nghiên cứu niên đại thực vật được khởi xướng bởi Douglass (1901) và Fritts (1987). Jacoby và D’Arrigo (1990), Buckley và cs (1995), Pumijumnong và cs đã chứng minh cây gỗ nhiệt đới Đông Nam Á hình thành vòng năm rõ rệt. Cook (1985) chuẩn hóa bề rộng vòng năm ($Z_t$) nhằm tách biệt biến động sinh học do tuổi cây khỏi tác động của khí hậu. Phương pháp ước lượng hợp lý cực đại phi tuyến (Nonlinear Maximum Likelihood) có trọng số và thẩm định chéo K-Fold được ứng dụng để mô hình hóa sinh trưởng (Swanson và cs, 2011; Pinheiro và cs, 2014; Huy và cs, 2016).
  • Tại Việt Nam, các nghiên cứu vòng năm của Sano và cs, Buckley và cs trên loài Pơmu (Fokienia hodginsii), Hansen và cs trên loài Thiết sam giả lá ngắn (Pseudotsuga sinensis), Phạm Trọng Nhân và cs (2011) trên loài Thông 3 lá (Pinus kesiya), Nguyễn Thị Oanh và cs (2015) và Nguyễn Văn Thiết (2016) trên loài Du sam (Keteleeria evelyniana) đã làm rõ tương quan giữa sinh trưởng vòng năm với nhiệt độ, lượng mưa và dao động khí hậu ENSO.

Về ứng dụng Hệ thống thông tin địa lý (GIS):

  • Aronoff (1989), Franklin (2001), Campbell và cs phân tích vai trò của GIS trong tích hợp dữ liệu không gian. Tại Việt Nam, phương pháp lồng ghép mô hình hồi quy đa biến vào các lớp dữ liệu không gian GIS để phân cấp thích nghi và lập bản đồ mật độ thực vật được thực hiện bởi Bảo Huy (2009, 2013, 2015a, 2017a), Đặng Ngọc Quốc Hưng và Hồ Đắc Thái Hoàng (2009), Phùng Văn Khen và Phạm Trịnh Hùng (2017).

Về các công trình nghiên cứu riêng trên loài Thông 5 lá:

  • Trên thế giới, loài Pinus dalatensis thuộc phân chi Strobus, nhánh Strobi (Richard và cs). Stefanie và Bond (2000) giải phẫu hình thái lá; Thomas và cs (2006) ghi nhận phân bố mới tại Phou Ak (Lào); Averyanov và cs định danh biến chủng; Thomas và Phan (2013) đánh giá hiện trạng bảo tồn.
  • Tại Việt Nam, sau công bố của Ferré (1960), Businsky (1999, 2004) phân chia thành 3 phân loài: Pinus dalatensis var. dalatensis, var. bidoupnensis và subsp. procera. Farjon (2002), Nguyễn Hoàng Nghĩa (2004), Nguyễn Đức Tố Lưu và Thomas (2004), Phan Kế Lộc và cs (2010, 2011, 2013) khảo sát vùng phân bố tại Tây Nguyên, Quảng Bình, Thừa Thiên Huế. Nguyễn Hoàng Nghĩa (2004) và Phí Hồng Hải (2011) cho rằng quần thể Thông 5 lá có nguy cơ thoái hóa do vắng bóng cây cấp kính nhỏ; trong khi Lê Văn Vinh (2013) và Trang (2011) phát hiện Thông 5 lá là loài ưa sáng mạnh, tái sinh tốt ở các lỗ trống tán ($100–200\text{ m}^2$) và bìa rừng. Phong và cs phân tích cấu trúc di truyền quần thể.

Khoảng trống nghiên cứu được luận án giải quyết: Các công trình trước đây chỉ tiếp cận từng khía cạnh hình thái, phân loại hoặc ghi nhận hiện trạng cục bộ. Luận án giải quyết khoảng trống về một nghiên cứu tổng thể, định lượng trên toàn bộ các vùng phân bố chính ở Tây Nguyên bao gồm: cấu trúc không gian 3 chiều và tổ thành loài; tương tác sinh thái liên loài; mô hình hóa toán học tác động của tổ hợp sinh thái (độ cao, tầng dày đất, lượng mưa) đến mật độ cây; phản ứng của bề rộng vòng năm chuẩn hóa đối với nhiệt độ và lượng mưa theo chuỗi thời gian; hệ thống mô hình sinh trưởng đường kính và chiều cao có thẩm định sai số K-Fold; cùng hệ thống cơ sở dữ liệu bản đồ GIS cấp mật độ phục vụ bảo tồn.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án áp dụng kết hợp các khung lý thuyết lâm học định lượng, sinh thái học thực nghiệm, niên đại học thực vật và phân tích không gian GIS:

Khung lý thuyết và mô hình tính toán

  • Cấu trúc tổ thành loài và loài ưu thế: Sử dụng chỉ số giá trị quan trọng (IV%) của Curtis và McIntosh (1950): $$\text{IV%} = \frac{\text{N%} + \text{G%} + \text{F%}}{3}$$ Trong đó: N% là tỉ lệ mật độ loài so với tổng mật độ; G% là tỉ lệ tiết diện ngang ($G = \frac{\pi \cdot D^2}{40000}$) loài so với tổng tiết diện ngang; F% là tỉ lệ tần số xuất hiện của loài trong các ô đo đếm. Đối với tầng cây tái sinh, IV% được tính bằng trung bình cộng của N% và F%.
  • Mô phỏng cấu trúc phân bố $N/D$ và $N/H$: Dữ liệu thực nghiệm được chia cấp kính cự ly 10 cm (cấp nhỏ nhất 6–16 cm, trung bình 11 cm; 16–26 cm, trung bình 21 cm...) và cấp chiều cao cự ly 2 m (cấp nhỏ nhất 5–7 m, trung bình 6 m; 7–9 m, trung bình 8 m...). Kiểm định tính đồng nhất giữa các ô tiêu chuẩn bằng tiêu chuẩn Chi bình phương ($\chi^2$). Thử nghiệm 4 dạng phân bố lý thuyết:
    1. Mô hình Mayer: $\ln(N) = \alpha + \beta \cdot D$
    2. Phân bố khoảng cách (Distance distribution): $P(x) = (1 - \theta) \cdot (1 - \lambda) \cdot \lambda^{x-1}$ (với $x \ge 1$)
    3. Phân bố hình học (Geometric distribution): $P(x) = \theta \cdot (1 - \theta)^x$
    4. Phân bố Weibull: Hàm mật độ $f(x) = \frac{c}{b} \cdot \left(\frac{x - x_{\min}}{b}\right)^{c-1} \cdot \exp\left[-\left(\frac{x - x_{\min}}{b}\right)^c\right]$ Độ phù hợp của mô hình lý thuyết với dữ liệu quan sát được kiểm định bằng tiêu chuẩn $\chi^2$ ở mức ý nghĩa $P < 0,05$.
  • Cấu trúc mặt bằng: Sử dụng tiêu chuẩn $U$ của Clark và Evans (1954): $$U = \frac{\bar{r} - E(r)}{\sigma_r} = \frac{\bar{r} - \frac{1}{2\sqrt{\rho}}}{\frac{0,26136}{\sqrt{n \cdot \rho}}}$$ Trong đó $\bar{r}$ là khoảng cách trung bình từ cây chọn ngẫu nhiên đến cây gần nhất, $\rho$ là mật độ cây trên một đơn vị diện tích ($m^2$), $n$ là dung lượng mẫu ($n > 30$). Quy tắc phân loại: $U < -1,96$ (phân bố cụm); $-1,96 \le U \le 1,96$ (phân bố ngẫu nhiên); $U > 1,96$ (phân bố cách đều).
  • Lý thuyết niên đại thực vật và chuẩn hóa vòng năm: Chuỗi bề rộng vòng năm thô ($Z_r$) thu thập từ lõi khoan được loại trừ xu thế sinh trưởng tự nhiên theo tuổi để chuyển thành chỉ số bề rộng vòng năm chuẩn hóa ($Z_t$) theo thuật toán của Cook (1985). Chuỗi $Z_t$ được đối chiếu với chuỗi số liệu nhiệt độ ($T$) và lượng mưa ($P$) trung bình tháng trong nhiều năm bằng hồi quy phi tuyến có trọng số.
  • Mô hình hóa sinh trưởng cây rừng: Sử dụng phương pháp Maximum Likelihood phi tuyến có trọng số và đưa biến ngẫu nhiên theo vùng phân bố sinh thái vào mô hình (Random effects). Thử nghiệm các dạng phương trình sinh trưởng đường kính theo tuổi ($D/A$) như Mitscherlich ($D = D_m \cdot (1 - e^{-a \cdot A})$ với đường kính giới hạn $D_m = 300\text{ cm}$), mô hình quan hệ chiều cao - đường kính ($H/D$) dạng Power ($H = a \cdot D^b$) và mô hình tỉ lệ tăng trưởng đường kính ($P_d = a \cdot D^b$). Quy trình lựa chọn và kiểm tra mô hình áp dụng kỹ thuật thẩm định chéo K-Fold (chia 10 phần dữ liệu ngẫu nhiên để đánh giá sai số độc lập).
  • Mô hình hóa sinh thái và liên kết GIS: Thiết lập hàm tương quan phi tuyến giữa cấp mật độ cá thể ($N$) với 3 biến môi trường độc lập: Độ cao so với mực nước biển (DC, m), Tầng dày đất (TDD, cm) và Lượng mưa trung bình năm (P, mm/năm). Kết quả mô hình được dùng làm thuật toán chồng xếp các lớp dữ liệu không gian trong phần mềm GIS để phân vùng mật độ.

Phương pháp thu thập dữ liệu thực tế

  • Dữ liệu cấu trúc lâm phần: Lập 17 ô tiêu chuẩn (ÔTC) điển hình tạm thời diện tích $2.500\text{ m}^2$ ($50\text{ m} \times 50\text{ m}$) tại các lâm phần ít bị tác động, gồm 6 ÔTC tại VQG Bidoup - Núi Bà, 6 ÔTC tại VQG Chư Yang Sin và 5 ÔTC tại VQG Kon Ka Kinh. Mỗi ÔTC $2.500\text{ m}^2$ được chia thành 25 ô đơn vị $100\text{ m}^2$ ($10\text{ m} \times 10\text{ m}$) để đo đếm toàn diện đường kính $D \ge 6\text{ cm}$, chiều cao vút ngọn $H$ và khoảng cách đến cây gần nhất. Đo tái sinh trên 5 ô dạng bản $4\text{ m}^2$ ($2\text{ m} \times 2\text{ m}$) ở 4 góc và tâm ô tiêu chuẩn đối với các cá thể có $H \ge 0,5\text{ m}$ và $D < 6\text{ cm}$.
  • Dữ liệu sinh thái và mật độ quần thể: Bố trí 19 điểm nghiên cứu diện tích $1\text{ km}^2$, trong mỗi điểm thiết lập các ô mẫu phân tầng $1.000\text{ m}^2$ để điều tra số lượng cây gỗ, cây tái sinh Thông 5 lá và đo đạc các yếu tố lập địa (độ cao, độ dốc, hướng phơi, tầng dày đất, lượng mưa).
  • Dữ liệu vòng năm: Sử dụng ống khoan tăng trưởng Haglöf khoan lấy mẫu lõi gỗ ngang tâm thân cây tại vị trí $1,3\text{ m}$. Các mẫu lõi khoan được dán nẹp gỗ, đánh bóng bề mặt, định tuổi chéo (cross-dating) và đo bề rộng từng vòng năm bằng kính hiển vi quang học độ phóng đại 40X kết hợp hệ thống bàn trượt đo Velmex kết nối máy tính.
Hạng mục điều tra Phương pháp và dung lượng mẫu Địa điểm thực hiện
Ô tiêu chuẩn cấu trúc lâm phần 17 ÔTC ($2.500\text{ m}^2$/ô, chia 25 ô nhỏ $100\text{ m}^2$) 6 ô Bidoup - Núi Bà, 6 ô Chư Yang Sin, 5 ô Kon Ka Kinh
Ô dạng bản điều tra tái sinh 85 ô dạng bản ($4\text{ m}^2$/ô, 5 ô/ÔTC cấu trúc) Bố trí tại 4 góc và tâm của 17 ÔTC $2.500\text{ m}^2$
Điểm nghiên cứu phân bố và sinh thái 19 điểm nghiên cứu ($1\text{ km}^2$/điểm) và ô mẫu $1.000\text{ m}^2$ 3 vùng phân bố tại Đắk Lắk, Lâm Đồng, Gia Lai
Mẫu lõi khoan vòng năm cây gỗ Khoan tăng trưởng Haglöf, đo kính hiển vi 40X và máy Velmex Các cây mẫu Thông 5 lá tại 3 Vườn quốc gia
Chuỗi số liệu khí hậu đối chứng Nhiệt độ ($T$) và lượng mưa ($P$) trung bình tháng nhiều năm Trạm khí tượng Bidoup (38 năm), Chư Yang Sin (32 năm), Kon Ka Kinh (32 năm)

Nội dung chính theo từng chương

Mở đầu

Phần mở đầu (3 trang) trình bày lý do chọn đề tài, tính cấp thiết của việc bảo tồn loài Thông 5 lá (Pinus dalatensis Ferré), xác định mục tiêu nghiên cứu tổng quát và mục tiêu cụ thể, đối tượng, phạm vi nghiên cứu, ý nghĩa khoa học và thực tiễn của luận án.

Chương 1: Tổng quan vấn đề nghiên cứu

Chương 1 (24 trang) tổng thuật có hệ thống các công trình nghiên cứu trong và ngoài nước liên quan đến cấu trúc quần thể thực vật rừng nhiệt đới, sinh thái rừng và quan hệ liên loài, phương pháp phân tích vòng năm cây rừng gắn với dao động khí hậu, kỹ thuật mô hình hóa sinh trưởng sinh trắc học và ứng dụng công nghệ thông tin địa lý (GIS). Chương này cũng phân tích toàn diện lịch sử nghiên cứu riêng về loài Thông 5 lá và kết thúc bằng phần thảo luận làm rõ các vấn đề học thuật cần được giải quyết trong luận án.

Chương 2: Nội dung, phương pháp và đặc điểm khu vực nghiên cứu

Chương 2 (25 trang) chi tiết hóa 6 nội dung nghiên cứu khoa học, xây dựng khung phương pháp tiếp cận định lượng toàn diện từ rút mẫu ngoại nghiệp, kỹ thuật xử lý số liệu toán thống kê ($\chi^2$, kiểm định $U$, phân bố Weibull, phân bố khoảng cách), phương pháp phân tích niên đại thực vật (Dendrochronology), thuật toán thẩm định chéo sai số K-Fold đối với các mô hình sinh trưởng phi tuyến, đến quy trình xây dựng lớp bản đồ mật độ trong môi trường GIS. Chương này cũng mô tả các điều kiện lập địa, khí hậu, thủy văn và thảm thực vật tự nhiên của 3 khu vực nghiên cứu tại VQG Bidoup - Núi Bà, VQG Chư Yang Sin và VQG Kon Ka Kinh.

Chương 3: Kết quả nghiên cứu và thảo luận

Chương 3 (78 trang) là nội dung trung tâm của luận án, trình bày và thảo luận các phát hiện thực nghiệm thông qua 32 bảng số liệu và 55 hình minh họa:

  1. Cấu trúc lâm phần có Thông 5 lá phân bố:
  • Thành phần loài: Trong kiểu rừng hỗn giao lá rộng lá kim, Thông 5 lá là một trong những loài ưu thế về sinh khối và tầng tán, xuất hiện cùng các loài lá kim và lá rộng đặc trưng như Thông 2 lá dẹt (Pinus krempfii), Pơmu (Fokienia hodginsii), Hồng tùng (Dacrydium elatum). Chỉ số quan trọng IV% của Thông 5 lá dao động tùy vùng sinh thái. Về tầng cây tái sinh, IV% của Thông 5 lá đạt giá trị thấp trong các lâm phần khép tán, phản ánh đặc tính không chịu bóng dưới tán rừng rậm.
  • Cấu trúc số cây theo cấp kính (N/D): Kiểm định $\chi^2$ cho thấy các dãy phân bố N/D của lâm phần có sự đồng nhất trong từng nhóm trạng thái. Dạng phân bố N/D chung của lâm phần tuân theo quy luật giảm liên tục, mô phỏng tốt nhất bằng hàm phân bố khoảng cách và phân bố Weibull. Đối với riêng quần thể Thông 5 lá, phân bố N/D thể hiện dạng có đỉnh lệch hoặc phân bố giảm gián đoạn, trong đó tỷ lệ cá thể tập trung nhiều ở các cấp kính trung bình và lớn ($D > 40\text{ cm}$), số lượng cá thể ở cấp kính nhỏ ($D < 16\text{ cm}$) rất hạn chế trong lâm phần rừng kín.
  • Cấu trúc số cây theo cấp chiều cao (N/H): Phân bố N/H của lâm phần và riêng loài Thông 5 lá có dạng phân bố một đỉnh lệch trái đến gần đối xứng, được mô phỏng phù hợp nhất bằng hàm phân bố Weibull ($P_{\text{value}} < 0,05$ qua kiểm định $\chi^2$). Cây Thông 5 lá chiếm lĩnh tầng vượt tán và tầng ưu thế sinh thái của rừng với chiều cao phổ biến từ 20 m đến trên 30 m.
  • Cấu trúc mặt bằng: Kết quả tính toán chỉ số $U$ theo cự ly khoảng cách giữa các cây cho thấy: đối với toàn bộ lâm phần, cây gỗ phân bố theo kiểu ngẫu nhiên đến cụm nhẹ; riêng đối với quần thể loài Thông 5 lá, giá trị $U < -1,96$ ở hầu hết các ô tiêu chuẩn, khẳng định loài này có cấu trúc mặt bằng dạng cụm trên mặt đất rừng.
  1. Ảnh hưởng của nhân tố sinh thái đến mật độ quần thể Thông 5 lá:
  • Luận án đã xác định được mô hình hồi quy phi tuyến biểu diễn mối quan hệ giữa cấp mật độ Thông 5 lá ($N$, cây/ha) với 3 nhân tố sinh thái chi phối gồm: Độ cao so với mực nước biển (DC, m), Tầng dày đất (TDD, cm) và Lượng mưa trung bình năm (P, mm/năm).
  • Ma trận phân tích 27 tổ hợp lập địa từ 3 biến sinh thái cho thấy: Cấp mật độ cao ($N > 100\text{ cây/ha}$) hình thành tại các khu vực có độ cao trên $1.400\text{ m}$, lượng mưa trung bình năm cao ($P > 2.000\text{ mm/năm}$) và tầng đất mỏng đến trung bình trên các sườn dông núi cao. Cấp mật độ trung bình ($51–100\text{ cây/ha}$) và thấp ($10–50\text{ cây/ha}$) xuất hiện ở các đai cao trung bình, trong khi mật độ hiếm ($< 10\text{ cây/ha}$) xuất hiện ở các vùng giáp ranh giới hạn sinh thái dưới $1.000\text{ m}$.
  1. Mối quan hệ sinh thái loài giữa Thông 5 lá với các loài ưu thế:
  • Phân tích bảng tiếp liên thực nghiệm sinh thái học chỉ ra các mối quan hệ hỗ trợ, cạnh tranh và độc lập giữa Thông 5 lá với các loài cây gỗ đi kèm trong quần xã. Ở tầng cây gỗ, Thông 5 lá có mối quan hệ tương hỗ với một số loài lá kim chịu hạn lạnh trên sườn dông và độc lập hoặc cạnh tranh nhẹ về không gian dinh dưỡng với các loài cây lá rộng bản địa thuộc họ Dẻ (Fagaceae), họ Re (Lauraceae). Ở tầng cây tái sinh, Thông 5 lá chịu sự cạnh tranh gay gắt về ánh sáng từ thảm tươi, cây bụi và cây tái sinh lá rộng tiên phong.
  1. Bề rộng vòng năm và sinh trưởng theo khí hậu:
  • Dữ liệu vòng năm chuẩn hóa ($Z_t$) phản ánh nhịp điệu sinh trưởng đường kính hàng năm và có tương quan rõ nét với các biến động khí hậu vùng:
    • Tại VQG Bidoup - Núi Bà: Chỉ số $Z_t$ tương quan thuận có ý nghĩa với nhiệt độ trung bình tháng 6 ($T_6$) và lượng mưa tháng 11 ($P_{11}$).
    • Tại VQG Chư Yang Sin: $Z_t$ tương quan nghịch với tích nhiệt độ tháng 3 và tháng 4 ($T_3 \times T_4$) – giai đoạn cao điểm mùa khô hạn tại Tây Nguyên.
    • Tại VQG Kon Ka Kinh: $Z_t$ tương quan nghịch với nhiệt độ tháng 4 ($T_4$).
  • Quy luật sinh thái học rút ra: Sự gia tăng nhiệt độ trong mùa mưa thúc đẩy cường độ quang hợp và làm tăng sinh trưởng đường kính của Thông 5 lá; ngược lại, sự gia tăng nhiệt độ trong các tháng mùa khô gây căng thẳng hạn sinh lý, dẫn đến sự suy giảm bề rộng vòng năm.
  1. Hệ thống mô hình sinh trưởng và tăng trưởng:
  • Quan hệ chiều cao - đường kính: Mô hình Power dạng $H = a_i \cdot D^b$ được xác lập với các bộ tham số $a_i$ phân biệt theo từng vùng sinh thái (Bidoup - Núi Bà, Chư Yang Sin, Kon Ka Kinh), giải thích tốt mối quan hệ tăng trưởng $H$ theo $D$.
  • Sinh trưởng đường kính theo tuổi ($D/A$): Mô hình Mitscherlich với giới hạn đường kính $D_m = 300\text{ cm}$ có dạng $D = 300 \cdot (1 - e^{-a_i \cdot A})$ cho kết quả tối ưu nhất qua kiểm định sai số K-Fold. Tham số $a_i$ thay đổi phản ánh tốc độ sinh trưởng đường kính khác nhau giữa 3 vùng phân bố.
  • Tỉ lệ tăng trưởng đường kính ($P_d$): Mô hình Power $P_d = a \cdot D^{b_i}$ mô tả chính xác sự giảm dần của tỉ lệ tăng trưởng đường kính khi kích thước cây tăng lên.
  1. Cơ sở dữ liệu GIS và bản đồ phân bố mật độ:
  • Chuyển hóa mô hình sinh thái đa biến vào hệ thống GIS, luận án đã xây dựng bản đồ số phân vùng mật độ loài Thông 5 lá tại 3 VQG và toàn vùng Tây Nguyên theo 4 cấp: Mật độ cao ($> 100\text{ cây/ha}$), Mật độ trung bình ($51–100\text{ cây/ha}$), Mật độ thấp ($10–50\text{ cây/ha}$) và Hiếm/ít gặp ($< 10\text{ cây/ha}$). Bản đồ tích hợp đầy đủ số hiệu tiểu khu rừng và ranh giới quản lý lâm phần.

Kết luận, tồn tại và kiến nghị

Phần này (3 trang) đúc kết các kết luận khoa học chính, nêu rõ các điểm tồn tại về phạm vi thu thập số liệu và đưa ra các kiến nghị kỹ thuật quản lý rừng, bảo tồn nguồn gen.

Kết quả và những đóng góp mới

Đóng góp mới về mặt khoa học và lý luận

  1. Luận án đã xây dựng thành công mô hình toán học dự đoán mật độ cá thể Thông 5 lá (Pinus dalatensis Ferré) dưới tác động đồng thời của 3 nhân tố sinh thái chính: Độ cao so với mực nước biển (DC), Tầng dày đất (TDD) và Lượng mưa trung bình năm (P), định lượng được điều kiện lập địa tối ưu hình thành nên các quần thể mật độ cao.
  2. Công trình đã làm sáng tỏ cơ chế phản ứng sinh trưởng của loài Thông 5 lá đối với nhân tố khí hậu thông qua phân tích chuỗi bề rộng vòng năm chuẩn hóa ($Z_t$): Gia tăng nhiệt độ trong mùa mưa có tác dụng kích thích sinh trưởng, trong khi gia tăng nhiệt độ trong mùa khô hạn gây ức chế và làm giảm sinh trưởng đường kính của loài.
  3. Luận án đã thiết lập và thẩm định chéo sai số độc lập (K-Fold cross-validation) hệ thống mô hình sinh trưởng đường kính ($D/A$), chiều cao ($H/D$) và tỉ lệ tăng trưởng đường kính ($P_d/D$) có xét đến hiệu ứng ngẫu nhiên của các vùng sinh thái khác nhau ở Tây Nguyên.

Đóng góp mới về mặt thực tiễn và giải pháp đề xuất

  1. Xác lập cơ sở dữ liệu GIS và hệ thống bản đồ phân bố cấp mật độ Thông 5 lá chi tiết đến từng tiểu khu rừng tại VQG Bidoup - Núi Bà, VQG Chư Yang Sin và VQG Kon Ka Kinh, làm công cụ phục vụ việc tuần tra, giám sát và phân vùng bảo tồn nghiêm ngặt loài cây quý hiếm.
  2. Đề xuất các biện pháp kỹ thuật lâm sinh ứng dụng:
    • Xúc tiến tái sinh tự nhiên: Khắc phục tình trạng thiếu hụt cây tái sinh dưới tán rừng rậm bằng giải pháp tạo các lỗ trống tán nhân tạo ($100–200\text{ m}^2$) với độ tàn che thích hợp từ $0,3$ đến $0,5$ ở các vùng gần cây mẹ thành thục, kết hợp luỗng phát thực bì, dây leo, cây bụi cạnh tranh để tạo điều kiện ánh sáng cho cây tái sinh Thông 5 lá sinh trưởng.
    • Gây trồng phục hồi quần thể: Sử dụng bản đồ phân vùng sinh thái thích hợp (kết hợp các điều kiện tối ưu về độ cao, tầng dày đất và lượng mưa) để lựa chọn địa bàn trồng bổ sung, làm giàu rừng bằng Thông 5 lá theo phương thức hỗn giao với các loài cây lá rộng và lá kim bản địa tương hợp sinh thái.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của công trình:

  • Đối tượng nghiên cứu của luận án chỉ tập trung trong phạm vi 3 Vườn quốc gia thuộc khu vực Tây Nguyên (Bidoup - Núi Bà, Chư Yang Sin, Kon Ka Kinh), chưa bao quát được các điểm phân bố tự nhiên khác ngoài vùng Tây Nguyên như khu vực dãy núi Mang (Thừa Thiên Huế), Quảng Bình hoặc các quần thể phân bố tại Lào.
  • Chuỗi dữ liệu khí hậu sử dụng để phân tích tương quan vòng năm phụ thuộc vào chuỗi thời gian quan trắc thực tế của các trạm khí tượng thủy văn lân cận khu vực nghiên cứu (từ 32 đến 38 năm), chưa đủ dài để tái lập toàn diện lịch sử cổ khí hậu hàng thế kỷ.

Hướng nghiên cứu tiếp theo:

  • Mở rộng phạm vi nghiên cứu lâm học và sinh thái quần thể Thông 5 lá sang các vùng sinh cảnh Bắc Trung Bộ và các khu bảo tồn lân cận trên dãy Trường Sơn.
  • Tiếp tục thu thập chuỗi niên đại vòng năm dài hơn trên các cây cổ thụ thành thục nhằm giải mã chi tiết tác động của các hiện tượng khí hậu cực đoan chu kỳ dài như ENSO, hạn hán lịch sử và biến đổi khí hậu toàn cầu đối với động thái quần thể loài Thông 5 lá.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị ứng dụng cho các nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và sinh viên các chuyên ngành Lâm sinh, Quản lý tài nguyên rừng, Sinh thái học thực vật và Thực vật học tham khảo về phương pháp tiếp cận định lượng trong cấu trúc rừng, kỹ thuật xử lý số liệu vòng năm (Dendrochronology), thuật toán thẩm định mô hình K-Fold và tích hợp mô hình sinh thái vào GIS.
  • Cán bộ kỹ thuật, ban quản lý các Vườn quốc gia, Khu bảo tồn thiên nhiên tại Tây Nguyên và vùng duyên hải miền Trung sử dụng hệ thống bản đồ mật độ, các tiêu chuẩn lập địa và biện pháp kỹ thuật mở tán ($100–200\text{ m}^2$, độ tàn che $0,3–0,5$) để tổ chức bảo tồn nguyên vị (in-situ) và gây trồng phục hồi nguồn gen Thông 5 lá.
  • Các nhà hoạch định chính sách lâm nghiệp và bảo tồn đa dạng sinh học trong việc xây dựng danh mục loài nguy cấp, thiết lập hành lang bảo tồn sinh thái và hoạch định chiến lược quản lý rừng thích ứng với biến đổi khí hậu.

Câu hỏi thường gặp

1. Ba nhân tố sinh thái chủ đạo nào quyết định cấp mật độ phân bố của Thông 5 lá trong mô hình dự đoán của luận án?

Luận án xác định 3 nhân tố sinh thái chủ đạo chi phối cấp mật độ ($N$) của Thông 5 lá là: Độ cao so với mực nước biển (DC), Tầng dày đất (TDD) và Lượng mưa trung bình năm (P). Trong đó, mật độ cao ($> 100\text{ cây/ha}$) hình thành tại các khu vực có độ cao trên $1.400\text{ m}$, lượng mưa trên $2.000\text{ mm/năm}$ và tầng đất mỏng đến trung bình trên sườn dông núi cao.

2. Sự biến động của nhiệt độ và lượng mưa ảnh hưởng cụ thể như thế nào đến bề rộng vòng năm chuẩn hóa ($Z_t$) của Thông 5 lá?

Kết quả phân tích hồi quy phi tuyến chuỗi niên đại vòng năm chỉ ra: Tại VQG Bidoup - Núi Bà, $Z_t$ tương quan thuận với nhiệt độ tháng 6 ($T_6$) và lượng mưa tháng 11 ($P_{11}$); tại VQG Chư Yang Sin, $Z_t$ tương quan nghịch với nhiệt độ tháng 3 và tháng 4 ($T_3 \times T_4$); tại VQG Kon Ka Kinh, $Z_t$ tương quan nghịch với nhiệt độ tháng 4 ($T_4$). Khái quát lại, tăng nhiệt độ trong mùa mưa thúc đẩy sinh trưởng, còn tăng nhiệt độ trong mùa khô hạn gây ức chế sinh trưởng đường kính của Thông 5 lá.

3. Mô hình toán học nào được lựa chọn để mô tả quy luật sinh trưởng đường kính theo tuổi ($D/A$) của Thông 5 lá?

Mô hình sinh trưởng Mitscherlich có dạng $D = 300 \cdot (1 - e^{-a_i \cdot A})$ (với đường kính giới hạn tối đa $D_m$ được cố định ở mức $300\text{ cm}$) là mô hình tối ưu nhất sau khi thực hiện thẩm định chéo sai số K-Fold. Tham số $a_i$ của mô hình được ước lượng phân biệt cho từng vùng phân bố (Bidoup - Núi Bà, Chư Yang Sin, Kon Ka Kinh) để phản ánh tốc độ sinh trưởng đặc trưng của từng vùng sinh thái.

4. Tại sao cây tái sinh Thông 5 lá thường vắng bóng dưới tán rừng kín và giải pháp lâm sinh khắc phục là gì?

Thông 5 lá là loài cây có đặc tính ưa sáng mạnh ngay từ giai đoạn cây con, do đó hạt rụng dưới tán rừng rậm có độ tàn che lớn không thể phát triển thành cây non triển vọng. Để xúc tiến tái sinh tự nhiên, luận án đề xuất giải pháp kỹ thuật mở tán tạo các khoảng trống có diện tích $100–200\text{ m}^2$ với độ tàn che duy trì từ $0,3$ đến $0,5$ quanh khu vực có cây mẹ gieo giống, đồng thời dọn dẹp thực bì cạnh tranh.

5. Đặc điểm cấu trúc mặt bằng của quần thể Thông 5 lá tại Tây Nguyên tuân theo quy luật nào?

Thông qua việc tính toán khoảng cách từ cây chọn ngẫu nhiên đến cây gần nhất theo tiêu chuẩn kiểm định $U$ của Clark và Evans ($n > 30$), luận án xác định giá trị $U < -1,96$ ở hầu hết các ô tiêu chuẩn, chứng minh các cá thể Thông 5 lá phân bố theo kiểu cụm trên mặt đất rừng tự nhiên.

Kết luận

Luận án của tác giả Lê Cảnh Nam là công trình nghiên cứu định lượng có hệ thống về đặc điểm lâm học, quy luật cấu trúc, động thái sinh trưởng vòng năm và sinh thái quần thể loài Thông 5 lá (Pinus dalatensis Ferré) tại khu vực Tây Nguyên. Nghiên cứu đã xác lập các mô hình toán học giải thích ảnh hưởng của tổ hợp sinh thái (độ cao, tầng dày đất, lượng mưa) đến mật độ cá thể, làm rõ tác động hai chiều của biến đổi nhiệt độ theo mùa đối với sinh trưởng đường kính và xây dựng hệ thống mô hình sinh trưởng chuẩn hóa qua kiểm định K-Fold. Các kết quả về giải pháp kỹ thuật mở tán xúc tiến tái sinh cùng hệ thống cơ sở dữ liệu bản đồ GIS cung cấp luận cứ khoa học và công cụ thực tiễn phục vụ công tác quản lý, bảo tồn nguyên vị và phát triển bền vững nguồn gen Thông 5 lá quý hiếm.