Bối cảnh và vấn đề nghiên cứu

Loài Thông đỏ bắc (Taxus chinensis (Pilg.) Rehder), thuộc họ Thông đỏ (Taxaceae), là loài thực vật hạt trần bản địa quý hiếm, phân bố rải rác tại một số tỉnh miền núi phía Bắc Việt Nam như Lào Cai, Sơn La, Hà Giang, Tuyên Quang, Cao Bằng, Hòa Bình và Thanh Hóa. Về mặt giá trị sinh học và kinh tế, loài cây này không chỉ cung cấp nguồn gỗ có giác lõi phân biệt, màu hồng thẫm, thớ mịn, chịu nước dùng trong mỹ nghệ và xây dựng, mà còn chứa các hoạt chất thứ cấp đặc biệt có giá trị dược liệu cao. Cụ thể, các hợp chất nhóm taxoid gồm 7-xylosyl-10-deacetylbaccatine III và 10-deacetylbaccatine III chiết xuất từ vỏ và lá của loài đã được chứng minh là tiền chất quan trọng để bán tổng hợp taxol - hoạt chất sử dụng trong phác đồ điều trị các bệnh ung thư.

Mặc dù có giá trị lớn, các quần thể Thông đỏ bắc trong tự nhiên đang đứng trước nguy cơ suy thoái và tuyệt chủng nghiêm trọng. Các quần thể tự nhiên thường có quy mô nhỏ, bị chia cắt không gian mạnh, khả năng tái sinh tự nhiên từ hạt rất hạn chế do nón cái chỉ phát triển đơn độc một hạt với áo hạt bao quanh, đồng thời chịu áp lực lớn từ nạn khai thác gỗ trái phép, bóc đẽo vỏ làm dược liệu, cháy rừng và sự thu hẹp sinh cảnh do canh tác nương rẫy. Hiện nay, loài này đã được xếp vào Phụ lục II của Công ước CITES, hạng Nguy cấp (EN) theo Danh lục đỏ của IUCN (2017), hạng Sẽ nguy cấp (VU) trong Sách Đỏ Việt Nam và thuộc Nhóm IIA theo Nghị định số 32/2006/NĐ-CP của Chính phủ.

Khu bảo tồn thiên nhiên Hang Kia – Pà Cò (huyện Mai Châu, tỉnh Hòa Bình) được xác định là khu vực có tính đa dạng sinh học cao với 750 loài thực vật bậc cao có mạch, trong đó có 41 loài nằm trong danh mục bị đe dọa tuyệt chủng. Sinh cảnh đặc trưng tại đây là kiểu phụ rừng kín thường xanh trên núi đá vôi á nhiệt đới, nơi lưu trữ các loài lá kim cổ như Thông Pà Cò (Pinus kwangtungensis), Pơ mu (Fokienia hodginsii), Dẻ tùng (Amentotaxus argotaenia) và Thông đỏ bắc (Taxus chinensis). Tuy nhiên, số lượng cá thể Thông đỏ bắc ghi nhận tại đây còn rất ít, phân bố rải rác trên các đỉnh núi cao hiểm trở, khả năng tái sinh tự nhiên thấp và chưa có quy hoạch phân vùng bảo tồn chuyên biệt. Xuất phát từ yêu cầu thực tiễn đó, đề tài nghiên cứu đã đặt ra các mục tiêu và nhiệm vụ cụ thể:

  1. Mục tiêu chung: Xây dựng cơ sở khoa học phục vụ công tác quy hoạch, định hướng không gian và đề xuất giải pháp trồng bảo tồn loài Thông đỏ bắc tại huyện Mai Châu, tỉnh Hòa Bình cũng như toàn bộ khu vực miền Bắc Việt Nam.
  2. Mục tiêu cụ thể:
    • Đánh giá đặc tính sinh học, sinh thái học và hiện trạng phân bố của quần thể Thông đỏ bắc tại Khu bảo tồn thiên nhiên Hang Kia – Pà Cò.
    • Ứng dụng mô hình sinh cảnh trên nền tảng hệ thống thông tin địa lý (GIS) để xác định các khu vực ưu tiên trồng bảo tồn loài.
    • Đề xuất hệ thống các giải pháp quản lý, kỹ thuật lâm sinh và kinh tế - xã hội nhằm bảo tồn và phát triển loài bền vững.
  3. Đối tượng và phạm vi nghiên cứu:
    • Đối tượng nghiên cứu: Loài Thông đỏ bắc (Taxus chinensis (Pilg.) Rehder) và sinh cảnh sống tự nhiên của loài.
    • Phạm vi không gian: Nghiên cứu thực địa chuyên sâu tại Khu bảo tồn thiên nhiên Hang Kia – Pà Cò (bao gồm địa giới hành chính của các xã Hang Kia, Pà Cò, Tân Sơn, Bao La, Piềng Vế, Cun Pheo, Nà Mèo thuộc huyện Mai Châu, tỉnh Hòa Bình) và mở rộng phạm vi dự báo sinh cảnh trên toàn địa bàn các tỉnh miền Bắc Việt Nam.
    • Phạm vi thời gian và dữ liệu: Dữ liệu ngoại nghiệp thu thập từ giai đoạn 2020–2021; dữ liệu viễn thám sử dụng ảnh vệ tinh Landsat 8 chụp ngày 28/06/2020 cùng các chuỗi số liệu khí hậu và hiện trạng rừng giai đoạn 2015–2020.

Cơ sở lý thuyết và phương pháp

Khung lý thuyết của nghiên cứu dựa trên quan điểm sinh thái học cảnh quan và mô hình hóa sinh cảnh (Habitat Modeling). Sinh cảnh của một loài thực vật được nhìn nhận là một thực thể đa chiều, cấu thành từ sự tương tác đồng thời giữa các yếu tố môi trường phi sinh học (địa hình, khí hậu, ánh sáng) và sinh học (lớp phủ thực vật). Việc phân tích thích nghi sinh cảnh trong đề tài sử dụng phương pháp thống kê đa biến Khoảng cách Mahalanobis ($D^2$). Đây là mô hình dự báo sinh thái dựa trên dữ liệu chỉ có điểm xuất hiện loài (presence-only model), cho phép định lượng mức độ sai khác giữa tổ hợp điều kiện môi trường tại một vị trí không gian bất kỳ với giá trị môi trường trung bình tối ưu được tính toán từ các tọa độ có sự hiện diện thực tế của loài.

Khoảng cách Mahalanobis ($D^2$) cho mỗi pixel trên bản đồ raster được xác định theo công thức ma trận:

$$D^2 = (x - m)^T \cdot C^{-1} \cdot (x - m)$$

Trong đó:

  • $x$ là véc-tơ các giá trị môi trường tại ô lưới (pixel) cần đánh giá.
  • $m$ là véc-tơ giá trị trung bình của các biến môi trường tại các điểm ghi nhận sự xuất hiện của Thông đỏ bắc trong thực địa.
  • $C^{-1}$ là nghịch đảo của ma trận hiệp phương sai giữa các biến môi trường độc lập tại các điểm xuất hiện loài.
  • $T$ là ký hiệu chuyển vị ma trận.

Giá trị $D^2$ tỷ lệ nghịch với mức độ thích hợp sinh cảnh: giá trị $D^2$ càng nhỏ thì điều kiện môi trường tại ô lưới càng gần với điều kiện môi trường chuẩn của loài, đồng nghĩa với tiềm năng phân bố và khả năng sống sót của loài càng cao. Nhằm trực quan hóa và lượng hóa xác suất thích nghi sinh thái về khoảng giá trị từ 0 đến 1, các giá trị $D^2$ được chuyển đổi sang phân phối Khi bình phương (Chi-square distribution, $p$-value). Sau đó, công cụ phân loại tự nhiên (Natural Breaks - Jenks) trong GIS được áp dụng để phân chia không gian thành 3 cấp độ ưu tiên trồng bảo tồn: Ưu tiên cao, Ưu tiên trung bình và Ưu tiên thấp.

Phương pháp nghiên cứu kết hợp giữa điều tra ngoại nghiệp và phân tích không gian nội nghiệp:

  • Kế thừa dữ liệu: Tiếp nhận bản đồ số độ cao (DEM), số liệu khí hậu trung bình năm (nhiệt độ, lượng mưa) từ nguồn WorldClim, cơ sở dữ liệu phân bố loài toàn cầu (GBIF) và các báo cáo hiện trạng tài nguyên rừng, kinh tế - xã hội của Ban quản lý Khu bảo tồn và Ủy ban nhân dân huyện Mai Châu.
  • Điều tra ngoại nghiệp: Thiết lập 4 tuyến khảo sát thực địa bằng thiết bị định vị GPS cầm tay, bao gồm: Tuyến đồi Hang Dơi; Tuyến Đỉnh Hang Kia; Tuyến đồi sau bể nước xã Hang Kia; và Tuyến khu đồi thuộc xã Pà Cò. Nội dung điều tra tập trung định vị tọa độ cá thể, đo đếm đường kính thân, chiều cao vút ngọn, đánh giá tầng cây tái sinh, tổ thành loài đi kèm và phỏng vấn cán bộ kiểm lâm địa bàn cùng người dân bản địa.
  • Xử lý nội nghiệp: Toàn bộ dữ liệu môi trường được chuẩn hóa về định dạng raster, độ phân giải không gian $30 \times 30\text{ m}$, quy chiếu về hệ tọa độ WGS 1984 UTM Zone 48N.
Tên biến môi trường Ký hiệu Mô tả đặc tính Nguồn cung cấp dữ liệu
Độ cao Elevation Độ cao tuyệt đối so với mực nước biển (m) Mô hình số độ cao DEM ($30 \times 30\text{ m}$) từ Bandolamnghiep.com
Độ dốc Slope Độ dốc sườn núi (độ) Tính toán phái sinh từ DEM
Hướng dốc Aspect Hướng phơi của sườn núi ($0^\circ - 360^\circ$) Tính toán phái sinh từ DEM
Bức xạ mặt trời Solar Insolation Chỉ số bức xạ/phơi sáng địa hình ($WH/m^2$) Tính toán phái sinh từ DEM qua phân tích bức xạ
Lượng mưa Rainfall Lượng mưa trung bình năm (mm) Dữ liệu khí hậu WorldClim.org
Nhiệt độ Temperature Nhiệt độ trung bình năm ($^\circ\text{C}$) Dữ liệu khí hậu WorldClim.org
Chỉ số thực vật NDVI Mức độ tươi tốt và mật độ tán che thảm thực vật Tính toán từ ảnh vệ tinh Landsat 8 (ngày chụp 28/06/2020)

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

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

Chương này tổng kết có hệ thống cơ sở phân loại học thực vật hạt trần và vị trí của chi Taxus trong ngành Thông (Pinophyta). Tác giả dẫn chứng tài liệu phân loại học quốc tế xác định chi Taxus gồm 24 loài trên thế giới; trong đó loài Taxus chinensis được Pilger và Rehder mô tả chính thức năm 1919 dựa trên mẫu chuẩn thu thập tại Trung Quốc. Về nghiên cứu ứng dụng dược liệu, chương tổng quan phân tích sâu các công bố của Trinh Thi Thuy và cộng sự (2005) chứng minh sự hiện diện của hoạt chất nhóm taxane tại Hòa Bình.

Tại Việt Nam, các dẫn liệu của Phan Kế Lộc, Nguyễn Tiến Hiệp chỉ rõ ngành Thông bản địa có 33 loài thuộc 19 chi, 5 họ, với trên 90% số loài đang bị đe dọa tuyệt chủng cấp quốc gia. Về mặt công nghệ, chương phân tích tiến trình phát triển của GIS từ các hệ thống máy tính lớn thập niên 1960 đến xu thế phân tích không gian nguồn mở và dữ liệu viễn thám độ phân giải cao hiện nay. Tác giả chỉ ra khoảng trống nghiên cứu rõ nét: mô hình khoảng cách Mahalanobis trên thế giới và tại Việt Nam (như nghiên cứu cây Mai dương của Lê Thái Sơn, 2019) phần lớn chỉ được ứng dụng để cảnh báo thực vật ngoại lai xâm hại; việc áp dụng mô hình này vào công tác xác định không gian ưu tiên bảo tồn loài thực vật quý hiếm bản địa là hướng tiếp cận hoàn toàn mới.

Chương 2: Đối tượng, mục tiêu, nội dung và phương pháp nghiên cứu

Nội dung chương trình bày chi tiết về đối tượng nghiên cứu là loài Thông đỏ bắc tại Khu bảo tồn thiên nhiên Hang Kia – Pà Cò. Tác giả cụ thể hóa hệ thống mục tiêu, thiết kế cấu trúc 5 nội dung nghiên cứu gắn liền với khung phương pháp luận GIS. Chương này giải thích cặn kẽ thuật toán chuyển đổi không gian đa chiều của Mahalanobis, quy trình tiền xử lý 7 lớp bản đồ chuyên đề raster có kích thước ô lưới $30 \times 30\text{ m}$, và việc ứng dụng phần mềm GIS kết hợp gói công cụ Land Facet Corridor Tools của Jenness nhằm tự động hóa quá trình tính toán nghịch đảo ma trận hiệp phương sai từ các điểm phân bố thực tế. Phương pháp điều tra thực địa theo 4 tuyến chính và phương pháp điều tra xã hội học thông qua phỏng vấn bán cấu trúc cũng được thiết lập cụ thể tại chương này.

Chương 3: Đặc điểm tự nhiên, kinh tế - xã hội khu vực nghiên cứu

Chương 3 lượng hóa toàn bộ điều kiện tự nhiên và kinh tế - xã hội của Khu bảo tồn thiên nhiên Hang Kia – Pà Cò, đơn vị quản lý có diện tích quy hoạch 5.252,98 ha nằm trên địa bàn 7 xã vùng cao huyện Mai Châu (ranh giới tọa độ từ $20^\circ 41'$ đến $20^\circ 46'$ vĩ độ Bắc; $104^\circ 05'$ đến $105^\circ 01'$ kinh độ Đông). Địa hình chia cắt hiểm trở với độ cao trung bình từ 800 đến 900 m, đỉnh cao nhất đạt 1.500 m, độ dốc sườn núi phổ biến từ $28^\circ$ đến $32^\circ$. Về địa chất, khu bảo tồn có hai dạng đá chính là trầm tích hóa học (đá vôi tạo đất feralit màu nâu dầy trên 120 cm) và trầm tích cơ học (đá phiến sét phong hóa thành đất feralit trên sa thạch). Khí hậu mang tính chất á nhiệt đới núi cao, nhiệt độ trung bình $21^\circ\text{C}$, lượng mưa bình quân 1.929 mm/năm (85% lượng mưa tập trung vào mùa mưa từ tháng 5 đến tháng 10).

Loại đất / Trạng thái rừng Diện tích (ha) Tỷ lệ (%)
I. Đất lâm nghiệp 5.244,07 99,83
1. Đất có rừng tự nhiên 5.221,38 99,40
- Rừng gỗ giàu 584,50 11,13
- Rừng gỗ trung bình 860,71 16,39
- Rừng gỗ nghèo 2.890,22 55,02
- Rừng gỗ phục hồi 885,95 16,86
2. Đất trống, đồi núi trọc quy hoạch lâm nghiệp 22,69 0,43
II. Đất khác 8,91 0,17
Tổng diện tích tự nhiên 5.252,98 100,00

Về tài nguyên sinh vật, khu hệ thực vật ghi nhận 750 loài thuộc 140 họ (trong đó có 35 loài trong Sách Đỏ Việt Nam 2007; 16 loài theo Nghị định 32/2006/NĐ-CP); hệ động vật có 451 loài với 11 loài quý hiếm đặc hữu. Về nhân sinh, khu vực có 12.207 nhân khẩu, dân tộc Mường chiếm 45,9%, còn lại là người H'Mông, Thái, Kinh và Tày. Tỷ lệ lao động nông nghiệp chiếm tới 86,9% (3.420 người trên tổng số 3.935 lao động); hệ thống thủy lợi gần như không có, canh tác phụ thuộc nước trời, dẫn tới đời sống nhân dân còn nhiều khó khăn và tạo sức ép trực tiếp lên tài nguyên rừng.

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

Chương này công bố toàn bộ dữ liệu định lượng thu được từ nghiên cứu. Về hình thái, loài cây tại khu vực nghiên cứu là cây thân gỗ nhỡ thường xanh, chiều cao thực tế đo được từ 8 đến 15 m, đường kính ngang ngực ($D_{1.3}$) dao động từ 15 đến 65 cm; vỏ ngoài màu nâu đỏ bong mảng nhỏ; lá hình dải mác cong lưỡi liềm dài 1,5 - 2,2 cm, rộng 3 mm; cây đơn tính khác gốc, nón cái đơn độc mang hạt hình trứng $8 \times 5\text{ mm}$ bao quanh bởi áo hạt đỏ. Khảo sát thực địa dọc 4 tuyến chỉ ghi nhận duy nhất 14 cá thể trưởng thành tại 7 bản: Pà Cò Lớn, Pà Háng Lớn (xã Pà Cò) và Hang Kia 2, Thung Ẳng, Thung Mặn, Thung Mài, Pà Khôm (xã Hang Kia), đều nằm ở đai cao từ 1.200 đến 1.500 m trên các sườn và đỉnh núi đá vôi. Tác giả chỉ ra rằng không tìm thấy cây cổ thụ lớn, tầng cây tái sinh cực kỳ thưa thớt, sinh trưởng còi cọc và bị cạnh tranh ánh sáng mạnh mẽ.

Phân tích mô hình sinh cảnh Mahalanobis cho thấy sự thích hợp của Thông đỏ bắc liên quan chặt chẽ với điều kiện vi khí hậu mát mẻ, lượng mưa cao và địa hình sườn dốc hiểm trở trên núi đá vôi. Bản đồ phân phối giá trị Khi bình phương ($p$-value) phản ánh không gian ưu tiên trồng bảo tồn loài được phân thành 3 cấp:

  • Khu vực ưu tiên cao: Tập trung tại các đỉnh dông cao thuộc xã Hang Kia và xã Pà Cò, nơi có điều kiện địa hình, tầng tán rừng thường xanh á nhiệt đới tương đồng nhất với 14 điểm phát hiện cá thể mẹ.
  • Khu vực mở rộng toàn miền Bắc: Mô hình định vị các dải sinh cảnh thích hợp cao chạy dọc theo các khối núi đá vôi hiểm trở thuộc vùng Tây Bắc và Đông Bắc (như Lào Cai, Hà Giang, Sơn La, Cao Bằng, Thanh Hóa).

Các mối đe dọa được tác giả tổng hợp từ hồ sơ kiểm lâm gồm tình trạng lấn chiếm rừng làm nương rẫy, khai thác trộm gỗ quý và hoạt động chăn thả gia súc tự do trong vùng lõi khu bảo tồn. Từ đó, chương đề xuất hệ thống 3 nhóm giải pháp bảo tồn toàn diện: tổ chức quản lý bảo vệ nghiêm ngặt quần thể hiện có, nâng cao sinh kế vùng đệm và ứng dụng kỹ thuật nhân giống sinh dưỡng bằng giâm hom cành (sử dụng chất điều hòa sinh trưởng IBA/IAA nồng độ 1.000 ppm) để tạo nguồn cây giống phục vụ trồng phục hồi sinh thái.

Kết quả và đóng góp

Kết quả nghiên cứu chính

  1. Dữ liệu phân bố và hiện trạng quần thể: Khóa luận đã xác định chính xác vị trí tọa độ của 14 cá thể Thông đỏ bắc trưởng thành còn sót lại trong tự nhiên tại 7 bản thuộc 2 xã Hang Kia và Pà Cò. Toàn bộ các cá thể này đều nằm ở đai cao từ 1.200 m đến 1.500 m trên sườn núi đá vôi dốc đứng. Quần thể phản ánh tình trạng cấu trúc suy giảm: không có cá thể đường kính lớn trên 65 cm, đường kính thân phổ biến từ 15 đến 65 cm, chiều cao 8 đến 15 m. Khả năng tái sinh tự nhiên từ hạt rất thấp, cây mầm dễ bị thối hỏng hoặc bị dẫm đạp, chết do thiếu hụt ánh sáng dưới tán rừng kín.
  2. Xây dựng thành công mô hình sinh cảnh số: Đề tài đã tích hợp đồng bộ 7 biến môi trường tự nhiên (Elevation, Slope, Aspect, Solar Insolation, Temperature, Rainfall, NDVI) vào công cụ tính toán Khoảng cách Mahalanobis ($D^2$) và chuyển đổi thành công sang bản đồ phân phối Khi bình phương ($p$-value). Kết quả phân loại tự nhiên Jenks đã xác lập bản đồ phân vùng ưu tiên trồng bảo tồn Thông đỏ bắc cho 2 phạm vi: nội tại Khu bảo tồn thiên nhiên Hang Kia – Pà Cò và không gian địa lý mở rộng trên toàn miền Bắc Việt Nam.
  3. Lượng hóa các mối đe dọa sinh cảnh: Nghiên cứu đã làm rõ các nguyên nhân trực tiếp đe dọa sự tồn vong của loài tại Hang Kia – Pà Cò, bao gồm sự cạnh tranh sinh thái tự nhiên của các loài đi kèm, nguy cơ suy giảm sinh cảnh do xâm lấn đất rừng làm nương rẫy, chăn thả gia súc tự do và tình trạng khai thác lâm sản trái phép.

Hệ thống giải pháp đề xuất

  • Nhóm giải pháp tổ chức - quản lý: Tăng cường lực lượng kiểm lâm và tổ bảo vệ rừng cộng đồng tuần tra nghiêm ngặt tại 7 bản có phân bố các cá thể mẹ; cắm mốc ranh giới bảo vệ tuyệt đối các cá thể trưởng thành; xử lý triệt để các hành vi phá rừng làm nương rẫy và chăn thả gia súc trong phân khu bảo vệ nghiêm ngặt.
  • Nhóm giải pháp kinh tế - xã hội: Giao khoán bảo vệ rừng gắn với chi trả dịch vụ môi trường rừng cho các hộ đồng bào dân tộc H'Mông, Mường; hỗ trợ phát triển sinh kế phi nông nghiệp, thâm canh lúa nước và cây ngô nhằm ổn định lương thực, giảm thiểu tối đa sự phụ thuộc sinh kế vào việc khai thác tài nguyên rừng.
  • Nhóm giải pháp kỹ thuật lâm sinh: Thu thập vật liệu giống (cành bánh tẻ) từ 14 cây mẹ hiện có để tiến hành giâm hom nhân giống sinh dưỡng bằng chất kích thích ra rễ (IBA hoặc IAA nồng độ 1.000 ppm) theo quy trình kỹ thuật đã được kiểm chứng; tiến hành trồng phục hồi, bổ sung nguồn giống vào đúng các vùng không gian được mô hình GIS chỉ định là "vùng ưu tiên cao" nhằm bảo đảm tỷ lệ sống và khả năng thích nghi cao nhất.

Đóng góp khoa học của đề tài

Khóa luận đã tạo ra một đóng góp mới về mặt phương pháp luận tại Việt Nam: chuyển hướng ứng dụng của mô hình sinh cảnh dựa trên thuật toán Khoảng cách Mahalanobis từ lĩnh vực quản lý sinh vật ngoại lai xâm hại sang lĩnh vực bảo tồn loài thực vật rừng nguy cấp, quý hiếm. Nghiên cứu cung cấp bản đồ số phân vùng ưu tiên có tọa độ chuẩn xác, giúp tối ưu hóa nguồn lực kinh phí và nhân lực cho công tác quy hoạch bảo tồn chuyển vị (ex-situ) và bảo tồn tại chỗ (in-situ) trên diện rộng mà phương pháp điều tra thực địa truyền thống không thể bao quát hết.

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

  • Hạn chế của nghiên cứu: Do điều kiện địa hình vùng núi đá vôi Hang Kia – Pà Cò bị chia cắt hiểm trở, nhiều vách đá dựng đứng khó tiếp cận và thời gian thực hiện đề tài tốt nghiệp có hạn, tác giả chỉ mới khảo sát được 4 tuyến điều tra ngoại nghiệp chính, chưa thể rà soát toàn bộ diện tích các đỉnh núi đá cao của khu bảo tồn. Số lượng mẫu điểm xuất hiện loài đưa vào mô hình còn khiêm tốn (14 điểm cá thể), hoàn toàn phụ thuộc vào việc tìm kiếm thực địa. Mặt khác, các lớp dữ liệu môi trường raster (như dữ liệu khí hậu WorldClim và mô hình DEM) có độ phân giải $30 \times 30\text{ m}$, do đó chưa phản ánh được hết các biến thiên cực nhỏ của điều kiện vi khí hậu cục bộ dưới các khe nứt, vách đá vôi ngách.
  • Hướng nghiên cứu tiếp theo:
    • Mở rộng phạm vi điều tra ngoại nghiệp tại các khu vực được mô hình dự báo thuộc cấp "Ưu tiên cao" ở các tỉnh lân cận như Sơn La, Lào Cai, Hà Giang và Thanh Hóa để kiểm chứng và hiệu chỉnh độ chính xác của mô hình sinh cảnh miền Bắc.
    • Tích hợp thêm các biến môi trường có độ chi tiết cao hơn, đặc biệt là các biến về tính chất lý hóa thổ nhưỡng núi đá vôi và dữ liệu vi khí hậu thu nhận trực tiếp từ các trạm đo chuyên dụng.
    • Thiết lập các ô tiêu chuẩn thực nghiệm trồng bảo tồn loài Thông đỏ bắc bằng cây con nhân giống từ hom tại các khu vực có giá trị $p$-value cao nhất để theo dõi, đánh giá tỷ lệ sống, khả năng sinh trưởng và hoàn thiện quy trình kỹ thuật lâm sinh phục hồi quần thể.

Giá trị tham khảo

Bản khóa luận tốt nghiệp là tài liệu tham khảo chuyên khảo có giá trị thiết thực đối với:

  • Cơ quan quản lý và lực lượng kiểm lâm: Ban quản lý Khu bảo tồn thiên nhiên Hang Kia – Pà Cò cùng các Chi cục Kiểm lâm vùng Tây Bắc có thể sử dụng trực tiếp bản đồ phân cấp ưu tiên để xây dựng phương án tuần tra chuyên đề, khoanh vùng bảo vệ nghiêm ngặt các dông núi có Thông đỏ bắc và lập dự án trồng phục hồi loài.
  • Cán bộ nghiên cứu và sinh viên: Sinh viên, học viên cao học thuộc các ngành Quản lý tài nguyên rừng, Quản lý tài nguyên thiên nhiên, Lâm học, Khoa học môi trường và Công nghệ địa không gian (GIS) có thể tham khảo phương pháp luận kết hợp thuật toán Mahalanobis với phần mềm GIS (Land Facet Corridor Tools) để xử lý bài toán mô hình hóa sinh cảnh loài khi chỉ có dữ liệu điểm xuất hiện (presence-only data).
  • Các dự án bảo tồn nguồn gen: Cung cấp thông tin thực chứng về phân bố, đặc tính hình thái học, cấu trúc tầng cây đi kèm và các dẫn liệu thử nghiệm giâm hom lâm sinh phục vụ các chương trình bảo tồn nguồn gen cây thuốc quý hiếm tại Việt Nam.

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

1. Tại sao loài Thông đỏ bắc (Taxus chinensis) lại có giá trị dược liệu đặc biệt quan trọng?

Vỏ và lá của cây Thông đỏ bắc chứa các hợp chất taxoid quan trọng, cụ thể là 7-xylosyl-10-deacetylbaccatine III và 10-deacetylbaccatine III (đã được Trinh Thi Thuy và cộng sự phân lập từ mẫu cây tại Hòa Bình năm 2005). Đây là các tiền chất trực tiếp phục vụ quy trình bán tổng hợp hoạt chất taxol – một trong những loại dược chất chính được sử dụng trong các phác đồ hóa trị liệu điều trị các bệnh ung thư trên thế giới.

2. Có bao nhiêu cá thể Thông đỏ bắc trưởng thành được phát hiện qua điều tra thực địa tại Khu bảo tồn thiên nhiên Hang Kia – Pà Cò?

Khảo sát thực địa dọc theo 4 tuyến điều tra tại khu bảo tồn đã ghi nhận được chính xác 14 cá thể Thông đỏ bắc trưởng thành. Các cây này phân bố rải rác trên sườn dông núi đá vôi ở đai cao từ 1.200 đến 1.500 m, thuộc địa bàn 7 bản: bản Pà Cò Lớn, Pà Háng Lớn (xã Pà Cò) và bản Hang Kia 2, Thung Ẳng, Thung Mặn, Thung Mài, Pà Khôm (xã Hang Kia). Toàn bộ cá thể có chiều cao 8 - 15 m và đường kính thân $15 - 65\text{ cm}$.

3. Mô hình sinh cảnh trong đề tài sử dụng bao nhiêu biến môi trường và gồm những biến cụ thể nào?

Mô hình sử dụng 07 biến môi trường dạng raster với độ phân giải không gian $30 \times 30\text{ m}$ quy chiếu UTM Zone 48N, bao gồm:

  1. Độ cao tuyệt đối (Elevation - trích xuất từ DEM).
  2. Độ dốc (Slope - tính từ DEM).
  3. Hướng dốc (Aspect - tính từ DEM).
  4. Chỉ số phơi sáng địa hình (Solar Insolation - tính từ DEM).
  5. Lượng mưa trung bình năm (Rainfall - từ WorldClim).
  6. Nhiệt độ trung bình năm (Temperature - từ WorldClim).
  7. Chỉ số thực vật (NDVI - tính toán từ ảnh viễn thám Landsat 8 chụp tháng 6/2020).

4. Vì sao tác giả lựa chọn phương pháp Khoảng cách Mahalanobis ($D^2$) thay vì các phương pháp thống kê đơn biến?

Các kỹ thuật thống kê đơn biến không phản ánh được bản chất đa chiều và mối tương quan phức tạp giữa các yếu tố môi trường cấu thành nên sinh cảnh sống của loài. Phép đo Khoảng cách Mahalanobis ($D^2$) cho phép tính toán khoảng cách đa chiều dựa trên ma trận hiệp phương sai của các biến môi trường tại các điểm loài xuất hiện thực tế (presence-only data). Nhờ đó, mô hình lượng hóa chính xác mức độ sai khác điều kiện sống ở mọi vị trí không gian so với điều kiện tối ưu của loài mà không đòi hỏi số liệu về các điểm vắng mặt loài (absence data) - vốn rất khó xác định chính xác đối với các loài quý hiếm, phân bố hẹp.

5. Những yếu tố nào đang trực tiếp đe dọa quần thể Thông đỏ bắc tại khu vực nghiên cứu?

Sự tồn vong của loài bị đe dọa bởi cả nhân tố nội tại lẫn tác động nhân sinh:

  • Nhân tố tự nhiên: Quần thể rất nhỏ (chỉ 14 cây), khả năng tái sinh tự nhiên từ hạt rất kém, cây con sinh trưởng chậm và bị chèn ép ánh sáng bởi các loài cây khác.
  • Nhân tố con người: Tình trạng chặt trộm lấy gỗ do gỗ có thớ mịn và chịu nước; bóc vỏ làm thuốc; lấn chiếm đất rừng núi đá vôi để canh tác nương rẫy; và hoạt động chăn thả gia súc tự do gây dẫm đạp làm chết cây tái sinh.

Kết luận

Khóa luận tốt nghiệp đã làm rõ thực trạng báo động của quần thể Thông đỏ bắc (Taxus chinensis) tại Khu bảo tồn thiên nhiên Hang Kia – Pà Cò với chỉ 14 cá thể trưởng thành còn sót lại và khả năng tái sinh tự nhiên suy giảm nghiêm trọng. Bằng việc ứng dụng thành công kỹ thuật mô hình hóa sinh cảnh đa chiều dựa trên hàm Khoảng cách Mahalanobis và phân phối Khi bình phương trong môi trường GIS, nghiên cứu đã phân lập được các không gian thích vi khí hậu tối ưu để phân cấp ưu tiên trồng bảo tồn loài tại khu bảo tồn cũng như mở rộng dự báo cho toàn miền Bắc Việt Nam. Kết quả đạt được cung cấp luận cứ khoa học tin cậy và công cụ bản đồ số hóa thiết thực, phục vụ trực tiếp cho công tác quy hoạch lâm sinh, bảo tồn chuyển vị và phục hồi nguồn gen thực vật nguy cấp có giá trị dược liệu cao này.