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" được thực hiện bởi nghiên cứu sinh Lê Cảnh Nam, 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), thuộc chuyên ngành Lâm sinh (mã số: 9 62 02 05) tại Viện Khoa học Lâm nghiệp Việt Nam, bảo vệ năm 2020.

Tính cấp thiết và khoảng trống nghiên cứu

Loài Thông 5 lá hay Thông 5 lá đà lạt (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 cây lá kim đặc hữu theo nghĩa rộng của dãy Trường Sơn và là loài cây gỗ lớn quý hiếm có giá trị khoa học và bảo tồn cao. Loài cây này được xếp vào Nhóm IIA theo Nghị định số 06/2019/NĐ-CP của Chính phủ Việt Nam và thuộc 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, 2019). Trên quy mô toàn cầu và toàn quốc, Thông 5 lá chỉ còn phân bố ở dưới 10 địa điểm chia cắt, với số lượng cá thể trưởng thành tại mỗi quần thể thường dưới 100 cây.

Khoảng trống nghiên cứu được tác giả xác định là: các công trình trước đây chủ yếu dừng lại ở mức độ phân loại thực vật, mô tả hình thái bên ngoài, ghi nhận điểm phân bố mới và khảo sát nhân giống đơn lẻ. Các dữ liệu định lượng có hệ thống về lâm học (như quy luật phân bố số cây theo cấp kính, cấp chiều cao, cấu trúc không gian mặt bằng, đặc điểm tái sinh tự nhiên), sinh thái quần thể (mối quan hệ tương tác sinh thái giữa Thông 5 lá với các loài cây ưu thế cùng xuất hiện, tổ hợp các nhân tố sinh thái quy định mật độ), và phản ứng sinh trưởng bề rộng vòng năm đối với biến động khí hậu theo chuỗi thời gian chưa được nghiên cứu toàn diện để làm cơ sở cho công tác quy hoạch bảo tồn và xử lý lâm sinh.

Mục tiêu nghiên cứu

  1. Mục tiêu chung: Cung cấp cơ sở lý luận và hệ thống thông tin, dữ liệu khoa học định lượng về lĩnh vực lâm học và sinh thái rừng phục vụ công tác bảo tồn và phát triển bền vững loài cùng các quần thể Thông 5 lá quý hiếm ở khu vực Tây Nguyên.
  2. Mục tiêu cụ thể:
    • Mục tiêu 1: Định lượng 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 2: Xác định 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). Cây gỗ lớn thường xanh có chiều cao đạt 30–40 m, đường kính ngang ngực lên đến 2,5 m.
  • Phạm vi không gian: Các lâm phần rừng tự nhiên hỗn giao lá rộng - lá kim có phân bố tự nhiên tập trung Thông 5 lá thuộc lâm phận quản lý của 3 Vườn quốc gia (VQG) tại Tây Nguyên: 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).
  • Phạm vi thời gian: Số liệu ngoại nghiệp và mẫu khoan vòng năm được thu thập trong giai đoạn từ năm 2017 đến năm 2018, đồng thời kế thừa một phần dữ liệu sinh thái tại VQG Chư Yang Sin từ đề tài khoa học công nghệ cấp tỉnh Đắk Lắk.

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

Các hướng nghiên cứu trong và ngoài nước được điểm luận

  1. Cấu trúc quần thể thực vật rừng: Richards (1952), Baur (1976), Odum (1971) đã đặt nền tảng mô tả tổ thành loài, dạng sống và phân tầng rừng mưa nhiệt đới. Curtis và McIntosh (1950) đề xuất chỉ số quan trọng IV% (Importance Value Index) để xác định loài ưu thế. Mô phỏng toán học cấu trúc không gian ba chiều được phát triển qua các hàm phân bố số cây theo cấp kính $N/D$ và cấp chiều cao $N/H$ như hàm Meyer (Meyer, 1943), hàm Weibull (Bailey, 1973; Rollet, 1971; Felfili, 1997; Baker et al., 2005; Barbosa et al., 2018), phân bố Khoảng cách, phân bố Poisson và Lognormal. Cấu trúc mặt bằng được Jayaraman (1999) phân loại thành ba dạng (ngẫu nhiên, cụm, đều) và được lượng hóa thông qua tiêu chuẩn $U$ của Clark và Evans (1954). Tại Việt Nam, Trần Ngũ Phương (1970), Thái Văn Trừng (1978, 1999), Nguyễn Văn Trương (1973, 1983), Đồng Sĩ Hiền (1974), Phùng Ngọc Lan (1986), Trần Văn Con (1991, 2001, 2011), Bảo Huy (1997, 2014, 2017a, 2017c), Võ Đại Hải (2014) và Bùi Mạnh Hưng, Lê Xuân Trường (2017) đã ứng dụng các mô hình phân bố xác suất trong phân tích cấu trúc rừng tự nhiên.
  2. Tái sinh rừng: Phùng Ngọc Lan (1986) và Thái Văn Trừng (1978) nhấn mạnh vai trò quyết định của ánh sáng đối với tái sinh rừng nhiệt đới. Các giải pháp phục hồi rừng thông qua tái sinh tự nhiên và nhân tạo được nghiên cứu bởi Lieth và Mooney (1991), Kamo et al. (2002), Lamb et al. (2005), Phạm Ngọc Thường (2003), Bùi Chính Nghĩa (2012), và Bảo Huy (2015a, 2015b).
  3. Mối quan hệ sinh thái loài và mô hình hóa sinh thái: Stephen và Garray (1986), Bùi Đoàn (2001), và Bảo Huy (2015a, 2017a) sử dụng phương pháp thực nghiệm sinh thái học và bảng phân tích 2x2 để xác định mức độ liên kết, tương tác sinh thái giữa các loài cây ưu thế. Mô hình hóa mối quan hệ phi tuyến giữa các nhân tố môi trường (độ cao, đất, lượng mưa) với mật độ và phân bố loài được chuẩn hóa bởi Bảo Huy (2009, 2014a, 2017a).
  4. Khoa học niên đại cây rừng (Dendrochronology) và mô hình sinh trưởng: Khởi xướng từ Andrew Ellicott Douglass (Fritts, 1987), nghiên cứu vòng năm được mở rộng sang vùng nhiệt đới qua các công trình của Jacoby và D'Arrigo (1990), Buckley et al. (1995, 2010, 2017), Biondi (1999), Cook (1985 - chuẩn hóa chỉ số $Z_t$), Pastur et al. (2008), và Barbosa et al. (2018). Tại Việt Nam, Sano et al. (2009) và Buckley et al. (2010, 2017) nghiên cứu vòng năm cây Pơ mu; Hansen et al. (2017) nghiên cứu Thiết sam giả lá ngắn; Phạm Trọng Nhân và cs (2011) nghiên cứu Thông 3 lá; Nguyễn Thị Oanh và cs (2015) nghiên cứu Pơ mu ở Kon Tum; Nguyễn Văn Thiết (2016) nghiên cứu Du sam. Phương pháp mô hình hóa phi tuyến Maximum Likelihood có trọng số và tác động ngẫu nhiên (Random effects) kết hợp kiểm định chéo (Cross-validation) được phát triển bởi Swanson et al. (2011), Pinheiro et al. (2014), Archontoulis và Miguez (2015), và Huy et al. (2016a, 2016b, 2019).
  5. Ứng dụng GIS trong quản lý lâm nghiệp: Aronoff (1989), Franklin (2001), Bảo Huy (2009, 2013, 2014), Đặng Ngọc Quốc Hưng và Hồ Đắc Thái Hoàng (2009), và Trịnh Thị Ngoan và cs (2018) ứng dụng GIS để tích hợp các biến môi trường và xây dựng bản đồ phân bố thích nghi của loài.
  6. Nghiên cứu về Thông 5 lá (Pinus dalatensis Ferré): Trên thế giới, các tác giả Richard et al. (2004), Stefanie và Bond (2000), Thomas et al. (2007), Businsky (1999, 2004, 2010), và Averyanov et al. (2014) chủ yếu tập trung mô tả phân loại học, giải phẫu hình thái lá, hạt và công bố vùng phân bố mới tại Lào (Khammoune, Xe Kong). Tại Việt Nam, các nghiên cứu của Farjon (2002), Phan Kế Lộc và cs (2010, 2011, 2013), Nguyễn Hoàng Nghĩa (2004), Phí Hồng Hải (2011), Lê Văn Vinh (2013), Trang (2011), Racz và Huyen (2007), Lương Viết Hùng và cs (2010), Lê Cảnh Nam và cs (2016), cùng Phong et al. (2014) khảo sát vùng phân bố, cấu trúc cục bộ và đa dạng di truyền.

Khoảng trống luận án lựa chọn giải quyết

Luận án giải quyết khoảng trống về sự thiếu hụt các mô hình định lượng cấu trúc lâm phần có Thông 5 lá phân bố, đặc điểm cấu trúc riêng của loài, bản chất mối quan hệ tương tác sinh thái giữa Thông 5 lá với các loài ưu thế, xác định tổ hợp nhân tố sinh thái quy định mật độ, quy luật sinh trưởng đường kính và phản ứng của bề rộng vòng năm đối với khí hậu theo chuỗi thời gian nhiều năm, cùng việc xây dựng cơ sở dữ liệu không gian GIS phục vụ quản lý và bảo tồn tại 3 Vườn quốc gia trọng điểm ở Tây Nguyên.

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

Khung lý thuyết và khái niệm sử dụng

  • Lý thuyết hệ sinh thái và phân tầng rừng nhiệt đới: Tiếp cận theo quan điểm sinh thái học của Odum (1971) và Thái Văn Trừng (1978, 1999), coi lâm phần hỗn giao lá rộng - lá kim là một hệ thống phân tầng chức năng với các nhóm loài ưu thế.
  • Chỉ số quan trọng loài (IV%): Xác định vị trí sinh thái của các loài cây gỗ và cây tái sinh theo Curtis và McIntosh (1950) cùng Narayan và Anshumali (2015): $$\text{IV% (cây gỗ)} = \frac{N% + G% + F%}{3}$$ $$\text{IV% (cây tái sinh)} = \frac{N% + F%}{2}$$ Trong đó: $N%$ là tỷ lệ mật độ loài, $G%$ là tỷ lệ tiết diện ngang tầng cây gỗ, $F%$ là tần suất xuất hiện loài qua hệ thống ô tiêu chuẩn.
  • Lý thuyết cấu trúc phân bố kích thước và không gian mặt bằng: Kiểm định tính đồng nhất của các dãy phân bố thực nghiệm bằng tiêu chuẩn $\chi^2$. Thử nghiệm 4 hàm phân bố xác suất lý thuyết: hàm Meyer, hàm Khoảng cách (Distance), hàm Hình học (Geometric) và hàm Weibull (đánh giá mức độ phù hợp ở ngưỡng xác suất $P < 0,05$). Kiểm tra dạng phân bố mặt bằng ngẫu nhiên, cụm hay cách đều bằng tiêu chuẩn $U$ của Clark và Evans (1954) dựa trên khoảng cách tới cây gần nhất: $$U = \frac{\bar{r} - E(r)}{\sigma_r}$$ (Nếu $-1,96 \le U \le 1,96$: phân bố ngẫu nhiên; $U > 1,96$: phân bố đều; $U < -1,96$: phân bố cụm).
  • Lý thuyết niên đại học thực vật (Dendrochronology): Tách ảnh hưởng của tuổi cây khỏi bề rộng vòng năm ($Z_r$) để thu được chuỗi chỉ số bề rộng vòng năm chuẩn hóa ($Z_t$) không thứ nguyên theo phương pháp của Cook (1985), làm đại lượng phản ánh biến động khí hậu.
  • Mô hình hóa sinh trắc rừng hiện đại: Sử dụng phương pháp ước lượng cực đại hợp lý phi tuyến có trọng số (Weighted Non-linear Maximum Likelihood), tích hợp biến ngẫu nhiên theo vùng phân bố sinh thái (Random effects), và thẩm định chéo sai số khách quan bằng phương pháp $K$-Fold cross-validation (chia $1/10$ dữ liệu kiểm định độc lập).

Phương pháp nghiên cứu thực tế và thu thập dữ liệu

  • Thiết kế mẫu điều tra cấu trúc lâm phần: Thiết lập 17 ô tiêu chuẩn (ÔTC) điển hình 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 Bidoup - Núi Bà, 6 ÔTC tại Chư Yang Sin và 5 ÔTC tại Kon Ka Kinh.
    • Tầng cây gỗ: Mỗi ÔTC $2.500\text{ m}^2$ chia thành 25 ô phụ $100\text{ m}^2$ ($10\text{ m} \times 10\text{ m}$) để đo toàn diện các cây có đường kính ngang ngực $D \ge 6\text{ cm}$, đo chiều cao $H\text{ (m)}$, xác định tên loài và cự ly đến cây gần nhất (đối với Thông 5 lá đo 2 cự ly: đến cây Thông 5 lá gần nhất và cây gỗ gần nhất).
    • Tầng cây tái sinh: Trên mỗi ÔTC $2.500\text{ m}^2$ bố trí 5 ô dạng bản $4\text{ m}^2$ ($2\text{ m} \times 2\text{ m}$) tại 4 góc và tâm ô để điều tra cây tái sinh triển vọng có $H \ge 0,5\text{ m}$ và $D < 6\text{ cm}$.
  • Thiết kế mẫu điều tra sinh thái và mật độ: Thiết lập 19 điểm nghiên cứu diện tích $1\text{ km}^2$ đại diện tại 3 vùng (Bidoup - Núi Bà: 5 điểm/50 ô; Chư Yang Sin: 4 điểm/40 ô; Kon Ka Kinh: 10 điểm/100 ô). Tại mỗi điểm $1\text{ km}^2$, bố trí 2 tuyến song song cách nhau $500\text{ m}$, mỗi tuyến đặt 5 ô mẫu tròn phân tầng diện tích $1.000\text{ m}^2$ ($R = 17,84\text{ m}$) và ô phụ đo tái sinh $100\text{ m}^2$ ($R = 5,64\text{ m}$). Tổng số ô mẫu điều tra thực tế là 173 ô tròn $1.000\text{ m}^2$, trong đó có 38 ô điều tra toàn bộ các loài cây gỗ đi kèm.
  • Dữ liệu vòng năm và thiết bị phân tích: Sử dụng khoan tăng trưởng Thụy Điển Haglöf để lấy mẫu lõi gỗ cây Thông 5 lá còn sống mà không chặt hạ. Lõi gỗ được dán lên máng gỗ cố định, đánh bóng bề mặt bằng giấy nhám mịn, soi định tuổi chéo và đo độ rộng từng vòng năm bằng bàn đo vi sai Velmex kết hợp kính hiển vi phóng đại 40X kết nối phần mềm máy tính chuyên dụng.
  • Xử lý dữ liệu không gian: Sử dụng hệ thống thông tin địa lý (GIS) để số hóa, nội suy và chồng xếp các lớp dữ liệu sinh thái nhằm chiết xuất ranh giới tiểu khu rừng và phân cấp mật độ quần thể.
Nội dung điều tra Số lượng mẫu thực tế Vùng Bidoup - Núi Bà Vùng Chư Yang Sin Vùng Kon Ka Kinh
Ô tiêu chuẩn cấu trúc ($2.500\text{ m}^2$) 17 ÔTC 6 ÔTC 6 ÔTC 5 ÔTC
Ô mẫu sinh thái & mật độ ($1.000\text{ m}^2$) 173 ô mẫu 48 ô mẫu 37 ô mẫu 88 ô mẫu
Điểm nghiên cứu quy mô không gian 19 điểm ($1\text{ km}^2$) 5 điểm 4 điểm 10 điểm
Phương pháp lấy mẫu vòng năm Khoan tăng trưởng Haglöf Thu thập mẫu lõi gỗ Thu thập mẫu lõi gỗ Thu thập mẫu lõi gỗ

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

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

Chương này tổng thuật toàn diện cơ sở lý luận lâm học nhiệt đới, cấu trúc lâm phần, sinh thái cá thể và quần thể rừng, khoa học niên đại vòng năm và ứng dụng GIS trong lâm nghiệp từ các công trình quốc tế và trong nước. Qua phân tích tổng quan, tác giả khẳng định các nghiên cứu trước về Thông 5 lá tại Việt Nam và Đông Dương chỉ tập trung vào phân loại học thực vật (Ferré, Businsky, Phan Kế Lộc) hoặc khảo sát hiện trạng định tính. Các luận điểm về diễn thế cho rằng quần thể Thông 5 lá bị "thoái hóa" do thiếu vắng cây con dưới tán rừng khép tán được phân tích lại: thực chất Thông 5 lá là loài ưa sáng mạnh ở giai đoạn non, đòi hỏi không gian mở tán hoặc lỗ trống rừng ($100–200\text{ m}^2$, độ tàn che $0,3–0,5$) để tái sinh tự nhiên thành công, chứ không phải loài suy thoái sinh học tự nhiê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 trình bày chi tiết 6 nội dung nghiên cứu, khung phương pháp thu thập mẫu ngoại nghiệp trên 17 ÔTC ($2.500\text{ m}^2$) và 173 ô mẫu tròn ($1.000\text{ m}^2$) phân bố tại 3 Vườn quốc gia đại diện cho 3 tiểu vùng sinh thái Tây Nguyên:

  • VQG Bidoup - Núi Bà (Lâm Đồng): Đại diện cho khối núi cao Langbiang, khí hậu á nhiệt đới núi cao, độ ẩm lớn.
  • VQG Chư Yang Sin (Đắk Lắk): Đại diện cho vùng núi cao phía Nam cao nguyên Đắk Lắk với địa hình chia cắt mạnh mẽ.
  • VQG Kon Ka Kinh (Gia Lai): Đại diện cho khu vực chuyển tiếp Bắc Tây Nguyên - Kon Tum, chịu ảnh hưởng của các chế độ khí hậu đặc thù.

Chương này cũng thiết lập toàn bộ quy trình toán học: kiểm định $\chi^2$, thuật toán ước lượng tham số hàm Weibull, tiêu chuẩn $U$ phân bố mặt bằng, quy trình chuẩn hóa chuỗi vòng năm $Z_t$, các phương trình phi tuyến ước lượng chiều cao $H/D$, sinh trưởng đường kính $D/A$, tỷ lệ tăng trưởng $P_d/D$ và quy trình xây dựng lớp bản đồ chuyên đề trong GIS.

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

1. Đặc điểm cấu trúc lâm phần có phân bố Thông 5 lá

  • Cấu trúc tổ thành loài: Tầng cây gỗ rừng tự nhiên hỗn giao lá rộng - lá kim có sự tham gia của Thông 5 lá cùng các loài lá kim quý hiếm khác như Thông 2 lá dẹt (Pinus krempfii), Pơ mu (Fokienia hodginsii), Thông 3 lá (Pinus kesiya), Hồng tùng (Dacrydium elatum) và các loài cây lá rộng bản địa. Giá trị $IV%$ khẳng định Thông 5 lá là một trong những loài ưu thế tầng tán chính trong các lâm phần nghiên cứu.
  • Cấu trúc cây tái sinh: Tái sinh tự nhiên của Thông 5 lá dưới tán rừng rậm rạp có độ tàn che cao là rất thấp; cây con tái sinh chỉ ghi nhận tập trung ở bìa rừng, đường dông, vách đá hở hoặc các khoảng trống mở tán.
  • Phân bố số cây theo cấp kính ($N/D$) và cấp chiều cao ($N/H$): Phân bố $N/D$ của toàn lâm phần phản ánh cấu trúc rừng tự nhiên nhiều tuổi. Đối với riêng loài Thông 5 lá, phân bố $N/D$ và $N/H$ biểu thị dạng phân bố có đỉnh (lệch trái hoặc tiệm cận chuẩn), mô tả chính xác nhất qua hàm phân bố xác suất Weibull.
  • Cấu trúc mặt bằng: Kiểm định tiêu chuẩn $U$ cho thấy phần lớn cây gỗ toàn lâm phần và riêng loài Thông 5 lá có kiểu phân bố cụm ($U < -1,96$) hoặc ngẫu nhiên ($-1,96 \le U \le 1,96$), phản ánh đặc tính sử dụng không gian dinh dưỡng tập trung theo vi địa hình thích hợp.

2. Ảnh hưởng của nhân tố sinh thái đến phân bố mật độ Thông 5 lá

Nghiên cứu xác định được 3 nhân tố sinh thái chủ đạo chi phối đến cấp mật độ phân bố ($N$, cây/ha) của Thông 5 lá:

  1. Độ cao so với mặt nước biển ($DC$, m).
  2. Tầng dày đất ($TDD$, cm).
  3. Lượng mưa trung bình năm ($P$, mm/năm).

Mô hình phi tuyến đa biến tổ hợp từ 27 tổ hợp biến sinh thái đã giải thích và dự báo chính xác các cấp mật độ của Thông 5 lá trên thực địa.

Nhân tố sinh thái Ký hiệu Vai trò ảnh hưởng đến phân bố mật độ
Độ cao địa hình $DC$ Ảnh hưởng đến nhiệt độ không khí, độ ẩm và sự phân hóa đai thực vật
Tầng dày đất $TDD$ Ảnh hưởng đến khả năng neo bám cơ giới và dinh dưỡng rễ cây gỗ lớn
Lượng mưa bình quân $P$ Chi phối cán cân ẩm độ và nhịp sinh học mùa của quần thể

3. Mối quan hệ sinh thái giữa Thông 5 lá với các loài ưu thế trong quần xã

Phân tích bảng kết hợp thực nghiệm sinh thái học chỉ ra các mối quan hệ hỗ trợ hoặc cạnh tranh giữa Thông 5 lá (cả tầng cây gỗ và lớp cây tái sinh) với các loài cây lá rộng và lá kim ưu thế cùng xuất hiện, làm cơ sở lựa chọn danh mục loài đi kèm khi phục hồi sinh thái.

4. Bề rộng vòng năm và sinh trưởng cây cá thể dưới ảnh hưởng của khí hậu

  • Tác động của khí hậu đến chỉ số vòng năm chuẩn hóa ($Z_t$):
    • Tại VQG Bidoup - Núi Bà: $Z_t$ tương quan thuận chặt chẽ với nhiệt độ tháng 6 ($T_6$) và lượng mưa tháng 11 ($P_{11}$) theo các phương trình hồi quy có trọng số $Z_t = f(T_6)$ và $Z_t = f(P_{11})$.
    • Tại VQG Chư Yang Sin: $Z_t$ tương quan nghịch với nhiệt độ mùa khô gồm tháng 3 và tháng 4 qua mô hình $Z_t = f(T_3, T_4)$ trong chuỗi số liệu 32 năm (1980–2011).
    • Tại VQG Kon Ka Kinh: $Z_t$ tương quan nghịch với nhiệt độ tháng 4 ($T_4$) qua mô hình $Z_t = f(T_4)$ giai đoạn 1980–2011.
    • Quy luật chung: Sự gia tăng nhiệt độ trong mùa mưa thúc đẩy sinh trưởng vòng năm của Thông 5 lá; ngược lại, nhiệt độ tăng cao trong mùa khô hạn gây stress nước và làm suy giảm sinh trưởng đường kính.
  • Hệ thống mô hình sinh trưởng và tăng trưởng:
    • Mô hình quan hệ chiều cao - đường kính ($H/D$): Dạng hàm Power $H = a_i \times D^b$, tham số $a_i$ phân hóa theo 3 vùng sinh thái.
    • Mô hình sinh trưởng đường kính theo tuổi ($D/A$): Hàm Mitscherlich có dạng: $$D = 300 \times \left(1 - e^{-a_i \times A}\right)$$ (với đường kính giới hạn tối đa $D_m = 300\text{ cm}$, tham số $a_i$ thay đổi theo vùng phân bố).
    • Mô hình tỷ lệ tăng trưởng đường kính ($P_d/D$): Hàm Power $P_d = a \times D^{b_i}$.
    • Toàn bộ các mô hình đều được thẩm định sai số khách quan qua phương pháp $K$-Fold cross-validation.
Mô hình sinh trắc học Dạng phương trình Ý nghĩa và tham số
Tương quan $H/D$ $H = a_i \times D^b$ Ước tính chiều cao $H\text{ (m)}$ theo đường kính $D\text{ (cm)}$
Sinh trưởng $D/A$ $D = 300 \times \left(1 - e^{-a_i \times A}\right)$ Đường cong Mitscherlich ước tính đường kính theo tuổi $A\text{ (năm)}$
Tăng trưởng $P_d/D$ $P_d = a \times D^{b_i}$ Ước tính tỷ lệ tăng trưởng đường kính $P_d\text{ (%/năm)}$ theo $D$

5. Bản đồ và cơ sở dữ liệu GIS về phân bố mật độ và sinh thái Thông 5 lá

Luận án tích hợp các phương trình toán học và điều kiện địa hình để xuất bản hệ thống bản đồ phân cấp mật độ Thông 5 lá chi tiết đến từng tiểu khu rừng tại 3 Vườn quốc gia:

  • VQG Bidoup - Núi Bà: Mật độ cấp Thấp ($10–50\text{ cây/ha}$) và cấp Hiếm ($<10\text{ cây/ha}$).
  • VQG Chư Yang Sin: Mật độ cấp Trung bình ($51–100\text{ cây/ha}$), cấp Thấp ($10–50\text{ cây/ha}$) và cấp Hiếm ($<10\text{ cây/ha}$).
  • VQG Kon Ka Kinh: Mật độ cấp Cao ($>100\text{ cây/ha}$), cấp Trung bình ($51–100\text{ cây/ha}$), cấp Thấp ($10–50\text{ cây/ha}$) và cấp Hiếm ($<10\text{ cây/ha}$).

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

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

  1. Xây dựng mô hình dự đoán mật độ quần thể sinh thái: Thiết lập thành công mô hình toán học lượng hóa mối quan hệ phi tuyến giữa cấp mật độ phân bố của Thông 5 lá với tổ hợp 3 nhân tố sinh thái chủ đạo: độ cao so với mặt nước biển, độ dày tầng đất và lượng mưa trung bình năm.
  2. Làm rõ quy luật phản ứng sinh trưởng đối với khí hậu: Phát hiện và lượng hóa tác động trái ngược của nhiệt độ theo mùa đối với bề rộng vòng năm chuẩn hóa ($Z_t$): gia tăng nhiệt độ trong mùa mưa thúc đẩy sinh trưởng tế bào tầng sinh sinh dưỡng gỗ, trong khi gia tăng nhiệt độ trong mùa khô hạn gây hạn sinh lý và kìm hãm sinh trưởng đường kính cây Thông 5 lá.
  3. Thiết lập hệ thống mô hình sinh trưởng đường kính chuẩn hóa: Áp dụng phương pháp ước lượng Maximum Likelihood phi tuyến có xét ảnh hưởng ngẫu nhiên của vùng sinh thái (Random effects) và thẩm định chéo sai số độc lập $K$-Fold để xây dựng chuỗi mô hình quan hệ $H/D$, sinh trưởng tích lũy $D/A$ (theo hàm Mitscherlich) và tỷ lệ tăng trưởng $P_d/D$ cho 3 vùng sinh thái Tây Nguyên.

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

  1. Cung cấp bản đồ và cơ sở dữ liệu không gian GIS: Cung cấp hệ thống bản đồ số và dữ liệu phân bố mật độ Thông 5 lá chi tiết theo tiểu khu rừng tại VQG Bidoup - Núi Bà, VQG Chư Yang Sin và VQG Kon Ka Kinh, phục vụ trực tiếp cho công tác tuần tra, cắm mốc ranh giới bảo tồn nghiêm ngặt nguồn gen quý hiếm.
  2. Đề xuất các biện pháp kỹ thuật lâm sinh:
    • Xúc tiến tái sinh tự nhiên: Mở các khoảng trống tán rừng nhân tạo với diện tích thích hợp ($100–200\text{ m}^2$, độ tàn che giảm xuống $0,3–0,5$) xung quanh cây mẹ thành thục để tạo điều kiện ánh sáng cho hạt nảy mầm và cây con sinh trưởng.
    • Trồng bổ sung và làm giàu rừng: Chọn lựa lập địa theo đúng 3 tiêu chí sinh thái (độ cao, độ dày tầng đất, lượng mưa) đã được mô hình hóa; kết hợp trồng xen với các loài cây lá rộng ưu thế có quan hệ sinh thái tương hỗ.

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

  • Hạn chế:
    • Phạm vi nghiên cứu thực địa chỉ tập trung tại 3 Vườn quốc gia lớn ở Tây Nguyên, chưa bao phủ toàn bộ các dải phân bố tự nhiên nhỏ lẻ khác của Thông 5 lá tại Việt Nam (như khu vực Núi Mang - VQG Bạch Mã, Thừa Lưu - Huế, Ngọc Linh - Kon Tum) và khu vực Nam Lào.
    • Số lượng ô mẫu điều tra sinh thái thực tế đạt 173 ô mẫu tròn (so với thiết kế ban đầu 190 ô) do địa hình núi cao hiểm trở và chia cắt phức tạp cản trở việc tiếp cận một số tuyến mẫu.
    • Số liệu khí hậu sử dụng trong phân tích chuỗi thời gian dựa trên các trạm khí tượng thủy văn lân cận khu vực nghiên cứu, chưa có trạm vi khí hậu đo đạc tự động liên tục ngay dưới tán lâm phần Thông 5 lá.
  • Hướng nghiên cứu tiếp:
    • Mở rộng nghiên cứu đối sánh sinh thái học quần thể và cấu trúc di truyền học phân tử của Thông 5 lá trên toàn bộ phạm vi phân bố địa lý dọc dãy Trường Sơn.
    • Thiết lập các ô định vị theo dõi diễn thế sinh thái lâu dài, nghiên cứu sâu về sinh lý học hạt giống, khả năng cộng sinh nấm rễ (Mycorrhiza) và thử nghiệm các mô hình trồng rừng thực nghiệm làm giàu rừng bằng Thông 5 lá.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị tham khảo chuyên môn cao đối với:

  • Nghiên cứu sinh, học viên cao học, sinh viên ngành Lâm sinh, Quản lý tài nguyên rừng và Sinh thái học: Tham khảo phương pháp luận kết hợp giữa điều tra sinh trắc học rừng, khoa học niên đại vòng năm (Dendrochronology) không chặt hạ cây, xử lý mô hình phi tuyến có trọng số và ứng dụng GIS trong sinh thái quần thể thực vật.
  • Ban quản lý các Vườn quốc gia và Khu bảo tồn thiên nhiên (đặc biệt là VQG Bidoup - Núi Bà, Chư Yang Sin, Kon Ka Kinh): Ứng dụng trực tiếp hệ thống cơ sở dữ liệu GIS và bản đồ phân cấp mật độ theo tiểu khu để quy hoạch các phân khu bảo vệ nghiêm ngặt nguồn gen Thông 5 lá.
  • Các cơ quan quản lý lâm nghiệp và nhà hoạch định chính sách: Cung cấp cơ sở khoa học định lượng để xây dựng quy trình kỹ thuật lâm sinh, hướng dẫn bảo tồn và nhân rộng các loài cây lá kim bản địa quý hiếm nằm trong Danh lục Đỏ tại Tây Nguyên.

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

1. Ba nhân tố sinh thái chủ đạo nào quy định mật độ phân bố của quần thể Thông 5 lá ở Tây Nguyên?

Qua phân tích đa biến tổ hợp từ 27 tổ hợp biến sinh thái, luận án chỉ ra 3 nhân tố sinh thái chính ảnh hưởng quyết định đến mật độ phân bố của Thông 5 lá là: độ cao so với mặt nước biển ($DC$), độ dày tầng đất ($TDD$), và lượng mưa trung bình năm ($P$).

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

Chỉ số $Z_t$ phản ánh rõ nét tác động khí hậu theo mùa: sự gia tăng nhiệt độ trong mùa mưa (ví dụ nhiệt độ tháng 6 tại Bidoup - Núi Bà) có tác dụng thúc đẩy sinh trưởng đường kính thân cây; ngược lại, sự gia tăng nhiệt độ trong mùa khô hạn (nhiệt độ tháng 3, tháng 4 tại Chư Yang Sin và Kon Ka Kinh) làm gia tăng mức độ thiếu hụt ẩm, gây kìm hãm sinh trưởng vòng năm của cây.

3. Phương trình toán học nào được lựa chọn để mô phỏng sinh trưởng đường kính ($D$) theo tuổi ($A$) của Thông 5 lá?

Luận án sử dụng hàm sinh trưởng Mitscherlich có dạng: $$D = 300 \times \left(1 - e^{-a_i \times A}\right)$$ với đường kính giới hạn tối đa được cố định $D_m = 300\text{ cm}$ và tham số tốc độ sinh trưởng $a_i$ phân hóa riêng biệt cho từng vùng phân bố sinh thái (Bidoup - Núi Bà, Chư Yang Sin, Kon Ka Kinh), được kiểm định sai số bằng phương pháp $K$-Fold cross-validation.

4. Tại sao cây tái sinh Thông 5 lá hầu như không xuất hiện dưới tán rừng rậm rạp?

Thông 5 lá là loài cây có đặc tính sinh thái ưa sáng mạnh ngay từ giai đoạn non. Dưới tán rừng rậm rạp có độ tàn che lớn, thiếu ánh sáng tán xạ và trực xạ, hạt không thể phát triển thành cây con triển vọng. Cây tái sinh chỉ xuất hiện dày đặc tại các khu vực mở tán, lỗ trống rừng có diện tích từ $100–200\text{ m}^2$ (độ tàn che $0,3–0,5$), ven đường dông hoặc sườn núi đá hở.

5. Cấp mật độ phân bố của Thông 5 lá khác nhau như thế nào giữa ba Vườn quốc gia nghiên cứu?

Theo bản đồ số hóa GIS:

  • VQG Bidoup - Núi Bà chỉ xuất hiện cấp mật độ Thấp ($10–50\text{ cây/ha}$) và cấp Hiếm ($<10\text{ cây/ha}$).
  • VQG Chư Yang Sin có cấp Trung bình ($51–100\text{ cây/ha}$), cấp Thấp và cấp Hiếm.
  • VQG Kon Ka Kinh là khu vực có mật độ phong phú nhất, xuất hiện cả 4 cấp gồm cấp Cao ($>100\text{ cây/ha}$), Trung bình, Thấp và Hiếm.

Kết luận

Luận án tiến sĩ của tác giả Lê Cảnh Nam là một công trình nghiên cứu khoa học định lượng toàn diện về đặc điểm lâm sinh học, cấu trúc không gian và sinh thái quần thể của loài Thông 5 lá (Pinus dalatensis Ferré) tại khu vực Tây Nguyên. Tác giả đã làm sáng tỏ tổ hợp 3 nhân tố sinh thái chủ đạo quy định mật độ phân bố, chứng minh quy luật tác động trái ngược của nhiệt độ mùa mưa và mùa khô tới bề rộng vòng năm, đồng thời thiết lập hệ thống phương trình sinh trưởng đường kính chuẩn hóa theo hàm Mitscherlich cho 3 vùng sinh thái. Các bản đồ chuyên đề GIS và đề xuất mở tán xúc tiến tái sinh tự nhiên từ luận án cung cấp luận cứ kỹ thuật trực tiếp cho công tác bảo tồn bền vững loài cây hạt trần quý hiếm này tại Việt Nam.