Tổng quan về luận án
Bảo tồn và phục hồi rừng phòng hộ đầu nguồn (RPHĐN) giữ vị trí then chốt trong an ninh nguồn nước và phát triển sinh thái bền vững tại các lưu vực sông trọng điểm miền Trung - Tây Nguyên. Luận án tiến sĩ lâm nghiệp của nghiên cứu sinh Huỳnh Văn Chung, chuyên ngành Lâm sinh (Mã số: 9620205) thực hiện tại Viện Khoa học Lâm nghiệp Việt Nam dưới sự hướng dẫn khoa học của PGS. Trần Văn Con và PGS. Nguyễn Danh, mang tiêu đề: "Nghiên cứu cơ sở khoa học để phục hồi và quản lý bền vững rừng phòng hộ đầu nguồn lưu vực sông Pô Kô, tỉnh Kon Tum". Nghiên cứu đặt trọng tâm vào lưu vực sông Pô Kô – phụ lưu huyết mạch của sông Sê San, trực tiếp cấp nước cho công trình thủy điện chiến lược hồ Yaly và đáp ứng nhu cầu nước đa ngành của khu vực (nông nghiệp chiếm 81,24%, môi trường 13,04%, sinh hoạt 3,74% và công nghiệp 1,98% theo Võ Ngọc Quỳnh Trâm và cs, 2015).
Khoảng trống nghiên cứu (research gap) được xác định rõ rệt: Mặc dù các lý thuyết thủy văn rừng kinh điển (Wischmeier, 1978; Hewlett, 1982; Võ Đại Hải, 1996, 2009) đã khẳng định vai trò điều tiết dòng chảy của thảm phủ thực vật, nhưng thiếu vắng các mô hình định lượng vi mô kết nối cấu trúc tầng tán lâm phần với khả năng giảm thiểu dòng chảy mặt ($DCM$) và xói mòn đất ($XM$) tại các tiểu lưu vực nhiệt đới gió mùa Tây Nguyên. Bên cạnh đó, các mô hình phục hồi truyền thống thường tách rời giữa can thiệp kỹ thuật lâm sinh và cơ chế kinh tế - xã hội chia sẻ lợi ích với cộng đồng bản địa.
Luận án giải quyết hệ thống 3 câu hỏi nghiên cứu và kiểm định 2 giả thuyết cốt lõi:
- Câu hỏi nghiên cứu 1 (RQ1): Thực trạng biến động cấu trúc lâm học, trữ lượng và mức độ suy thoái của các trạng thái RPHĐN lưu vực sông Pô Kô diễn ra như thế nào dưới tác động kép của tự nhiên và nhân sinh?
- Câu hỏi nghiên cứu 2 (RQ2): Mối quan hệ tương quan định lượng giữa các chỉ tiêu cấu trúc nội tại của rừng (độ tàn che, độ che phủ thảm tươi, thảm mục, tầng cao) với khả năng giữ nước, cản dòng chảy mặt và hạn chế xói mòn đất được biểu diễn qua quy luật nào?
- Câu hỏi nghiên cứu 3 (RQ3): Khung giải pháp lâm sinh kết hợp thể chế kinh tế nào đảm bảo phục hồi cấu trúc rừng mong muốn và duy trì sinh kế bền vững cho đồng bào thiểu số tại địa bàn?
- Giả thuyết 1 (H1): Khả năng phòng hộ đầu nguồn của hệ sinh thái rừng không phụ thuộc đơn lẻ vào diện tích che phủ mà bị chi phối trực tiếp bởi chỉ số cấu trúc tổng hợp đa tầng ($C_s$) kết hợp với tỷ lệ lập địa bề dày tầng đất/độ dốc ($Hd/s$).
- Giả thuyết 2 (H2): Phục hồi rừng phòng hộ chỉ đạt tính bền vững sinh thái khi tối ưu hóa cấu trúc lâm phần đạt ngưỡng $C_s \ge 1,7$ gắn liền với việc thể chế hóa quyền hưởng lợi lâm sản ngoài gỗ (LSNG) và chi trả dịch vụ môi trường rừng (DVMTR/PFES).
Nghiên cứu được triển khai trên quy mô toàn lưu vực sông Pô Kô với diện tích tự nhiên 316.676 ha (chiếm 32,7% diện tích tỉnh Kon Tum, diện tích có rừng 164.685 ha, trong đó RPHĐN chiếm 64.052 ha), tập trung chuyên sâu vào hai địa bàn trọng điểm là Ban quản lý RPH Đak Ang và Ban quản lý RPH Tu Mơ Rông. Với hệ thống 100 ô tiêu chuẩn sơ cấp diện tích $2.500,\text{m}^2/\text{ô}$ (tổng diện tích mẫu $25,\text{ha}$), công trình mang lại giá trị định lượng đột phá, đóng góp cơ sở khoa học trực tiếp vào công tác quy hoạch, phục hồi và bảo tồn tài nguyên rừng Tây Nguyên.
Literature Review và Positioning
Lý thuyết về thủy văn rừng và quản lý lưu vực trên thế giới đã trải qua quá trình tiến hóa từ quản lý đơn mục tiêu (khai thác gỗ bền vững - Sustained Yield) sang quản lý tích hợp hệ sinh thái đa chức năng (Ecosystem-based Multifunctional Management). Trường phái thủy văn rừng phương Tây với các công trình tiên phong của Ellison (1947), Wischmeier & Smith (1978), Hewlett (1982), Bosch & Hewlett (1982), Zhang et al. (2001, 2004) và Van Dijk et al. (2007) đã làm sáng tỏ cơ chế phân phối lượng giáng thủy trong hệ sinh thái rừng: nước giữ lại trên tán và bốc hơi (10 - 45%), nước chảy men thân (1 - 5%), nước lọt qua tán (50 - 90%), dòng chảy mặt (5 - 50%) và nước thấm sâu tạo dòng chảy ngầm (40 - 90%).
Trong lâm học nhiệt đới, các lý thuyết tái sinh và cấu trúc quần xã rừng của Richards (1952), Budowski (1965) và Rollet (1974) đặt nền tảng cho việc phân tích động thái tầng tán và diễn thế phục hồi. Tại Việt Nam, các nghiên cứu của Bùi Ngạnh (1977, 1985), Nguyễn Ngọc Lung (1992, 1997), Đỗ Đình Sâm (2002), Vương Văn Quỳnh (2007), Vũ Tấn Phương (2009) và Võ Đại Hải (1996, 2009) đã chứng minh thảm mục rừng lá rộng thường xanh nhiệt đới tích lũy đến 11,2 tấn/ha, có khả năng hút lượng nước gấp 1,38 đến 3,21 lần trọng lượng khô của nó (138,3% - 321%), giúp hình thành dòng chảy ngầm tới 88 - 92% và triệt tiêu dòng chảy mặt xuống mức dưới 2%.
Tuy nhiên, y văn quốc tế và trong nước tồn tại những tranh luận học thuật sâu sắc:
- Tranh luận về sự đánh đổi thủy văn (Hydrological Trade-off): Quan điểm truyền thống cho rằng tăng độ che phủ rừng luôn làm tăng lượng nước lưu vực. Ngược lại, các nghiên cứu của Gyenge et al. (2009) và Van Dijk et al. (2007) chỉ ra rằng các loài cây sinh trưởng nhanh hoặc rừng thường xanh có thể tiêu thụ nước qua thoát hơi tầng tán nhiều hơn 20% so với rừng rụng lá, dẫn đến hiện tượng sụt giảm dòng chảy kiệt mùa khô ở một số lưu vực hạn chế về lượng mưa.
- Tranh luận về mô hình quản lý tập trung hóa so với phân quyền cộng đồng: Trường phái quản lý lâm nghiệp cổ điển (Hartig, 1840) nhấn mạnh sự kiểm soát tuyệt đối của nhà nước thông qua các lệnh cấm can thiệp; trong khi trường phái lâm nghiệp xã hội và thể chế đồng quản lý (ITTO, 1990, 1992; Wandelli, 2000; RECOFTC) khẳng định rằng việc loại bỏ cộng đồng địa phương khỏi rừng phòng hộ đầu nguồn sẽ kích hoạt xung đột tài nguyên, gia tăng khai thác lậu và làm phá sản các mục tiêu bảo tồn.
So sánh với các nghiên cứu quốc tế lớn:
- Đánh giá 6 lưu vực sông lớn châu Á (Ấn, Hằng, Brahmaputra, Mê Kông, Trường Giang, Hoàng Hà - FAO/CIFOR): Nghiên cứu chỉ ra rằng trên 30% diện tích nằm ở cao độ thấp dưới 300m bị thâm canh nông nghiệp triệt để, trong khi rừng đầu nguồn trên 500 - 3.000m bị suy thoái do thiếu cơ chế giám sát vi mô. Luận án của Huỳnh Văn Chung đã giải quyết bài toán này ở quy mô tiểu lưu vực Pô Kô thông qua việc phân lập chi tiết trạng thái suy thoái rừng theo vi địa hình.
- Mô hình canh tác dốc SALT (Sloping Agricultural Land Technology) tại Philippines và băng phủ sinh học tại Indonesia: Luận án kế thừa và phát triển vượt bậc khi tích hợp mô hình nông lâm kết hợp dưới tán rừng phòng hộ với các loài bản địa có giá trị cao (Litsea glutinosa, Scaphium macropodum, Dialium cochinchinense), vượt qua hạn chế của các mô hình thuần nông nghiệp trên đất dốc.
Vị thế học thuật (positioning) của luận án nằm ở điểm giao thoa giữa lâm sinh học định lượng (Quantitative Silviculture) và kinh tế tài nguyên thể chế (Institutional Resource Economics), giải quyết trọn vẹn mối quan hệ nhân quả giữa cấu trúc tán rừng - động thái thủy văn - sinh kế người bản địa.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Luận án đóng góp quan trọng vào việc mở rộng Lý thuyết Động thái Cấu trúc Lâm phần (Forest Stand Dynamics Theory) và Lý thuyết Cân bằng Thủy văn Lưu vực (Catchment Hydrological Equilibrium Theory) trong bối cảnh rừng nhiệt đới gió mùa Tây Nguyên chịu suy thoái:
- Mở rộng lý thuyết cấu trúc tầng thứ rừng nhiệt đới: Luận án chứng minh rằng quy luật phân bố số cây theo chiều cao ($N/H$) và đường kính ($N/D$) của các trạng thái rừng lá rộng thường xanh (LRTX) suy thoái tuân theo các hàm phi tuyến Weibull và Meyer. Khi rừng bị tác động mạnh (chuyển sang trạng thái nghèo LRTX-N hoặc phục hồi LRTX-PH), hàm phân bố bị lệch đỉnh, phản ánh sự đứt gãy trong cấu trúc thế hệ cây kế cận.
- Phát triển khái niệm "Cấu trúc rừng mong muốn" ($C_s$): Đóng góp lý thuyết mang tính đột phá là việc xây dựng chỉ số cấu trúc tổng hợp liên quan trực tiếp đến khả năng phòng hộ:
$$C_s = f\left(\frac{TC}{H}, CP, TM\right)$$
Trong đó: $TC$ là độ tàn che tầng cây cao (%), $H$ là chiều cao vút ngọn bình quân (m), $CP$ là độ che phủ thảm tươi cây bụi (%), và $TM$ là độ che phủ/chiều dày lớp thảm mục (%). Luận án chứng minh sự chuyển dịch mô thức (paradigm shift): Khả năng giảm dòng chảy mặt và hạn chế xói mòn không tăng tuyến tính đơn biến theo độ tàn che $TC$ mà phụ thuộc vào sự tối ưu hóa đồng thời của cả ba tầng cấu trúc: tán cây cao, thảm tươi cây bụi và thảm mục hữu cơ mặt đất.
+-------------------------------------------------------------+
| KHUNG PHÂN TÍCH CẤU TRÚC RỪNG MONG MUỐN VÀ ĐIỀU TIẾT NƯỚC |
+-------------------------------------------------------------+
|
+----------------------------+----------------------------+
| |
v v
+------------------+ +--------------------+
| Cấu trúc tầng cao| | Cấu trúc mặt đất |
| - Độ tàn che TC | | - Thảm tươi CP |
| - Chiều cao Hvn | | - Thảm mục TM |
| - Tiết diện G | | - Rễ & dung trọng |
+------------------+ +--------------------+
| |
+----------------------------+----------------------------+
|
v
+------------------------------------------------------+
| CHỈ SỐ CẤU TRÚC LÂM PHẦN TỔNG HỢP LIÊN QUAN PHĐN: |
| Cs = [(TC / H) + CP + TM] |
+------------------------------------------------------+
|
+-------------------+-------------------+
| |
v v
+---------------------------+ +---------------------------+
| Khả năng giảm dòng chảy | | Khả năng cản xói mòn đất |
| mặt DCM (m3/ha/năm): | | XM (tấn/ha/năm): |
| DCM = f(Cs, Hd/s) | | XM = f(Cs, Hd/s) |
+---------------------------+ +---------------------------+
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp liên ngành 3 trục lý thuyết:
- Lâm học định lượng và đa dạng sinh học (dựa trên hàm cấu trúc Weibull/Meyer và đường cong đa dạng Renyi).
- Thủy văn học bề mặt đất dốc (mô hình hóa phi tuyến tương quan dòng chảy $DCM$ và xói mòn $XM$ với tương quan lập địa $Hd/s$, trong đó $Hd$ là độ dày tầng đất, $s$ là độ dốc).
- Kinh tế - xã hội lâm nghiệp (áp dụng ma trận SWOT và công cụ đánh giá có sự tham gia PRA để phân tích sinh kế và thể chế hưởng lợi).
Điều kiện biên (boundary conditions) của khung phân tích được xác lập rõ: Áp dụng cho các lưu vực sông miền núi nhiệt đới có lượng mưa hàng năm lớn ($>2.000,\text{mm}$), địa hình phân cắt mạnh (độ dốc bình quân 15 - 35 độ), có sự hiện diện của các kiểu rừng tự nhiên nhiệt đới gió mùa từ rừng lá rộng thường xanh (LRTX), rừng hỗn giao lá rộng - lá kim (LRLK) đến rừng tre nứa xen gỗ.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án áp dụng lập trường nhận thức Hiện thực phản biện (Critical Realism) kết hợp phương pháp hỗn hợp (Mixed Methods Research Design) đa cấp độ:
- Cấp độ vĩ mô (Macro-level): Ứng dụng công nghệ giải đoán ảnh vệ tinh đa thời gian bằng phương pháp phân tích ảnh hướng đối tượng (Object-Based Image Analysis - OBIA) trên phần mềm Trimble eCognition, kết hợp Hệ thống Thông tin Địa lý (GIS) để bóc tách 164.685 ha diện tích có rừng trên toàn lưu vực sông Pô Kô.
- Cấp độ trung mô (Meso-level): Thiết lập hệ thống 100 ô tiêu chuẩn sơ cấp điển hình ($50,\text{m} \times 50,\text{m} = 2.500,\text{m}^2/\text{ô}$) đại diện cho toàn bộ các trạng thái rừng phòng hộ (LRTX Giàu, Trung bình, Nghèo, Phục hồi, Rừng Tre nứa, Tre nứa hỗn giao gỗ).
- Cấp độ vi mô (Micro-level): Chia mỗi OTC sơ cấp thành 25 ô thứ cấp ($10,\text{m} \times 10,\text{m} = 100,\text{m}^2$) để đo đếm tầng cây cao ($D_{1,3} \ge 10,\text{cm}$) và 05 ô dạng bản ($4,\text{m} \times 4,\text{m} = 16,\text{m}^2$) tại 4 góc và trung tâm để đo đếm tầng cây tái sinh ($N_{ts}$), thảm mục ($TM$) và thảm tươi cây bụi ($CP$).
Quy trình nghiên cứu rigorous
Quy trình nghiên cứu tuân thủ nghiêm ngặt các tiêu chuẩn quốc tế và quốc gia:
- Chiết xuất viễn thám và mẫu khóa ảnh: Lựa chọn ngẫu nhiên 10 điểm kiểm chứng thực địa cho mỗi trạng thái rừng (cách ranh giới trạng thái tối thiểu 50m). Tại mỗi điểm, sử dụng thước đo góc Bitterlich xác định tiết diện ngang ($G$) tại 5 vị trí (tâm và 4 hướng Đông, Tây, Nam, Bắc cách 5m). Cây phân loại tự động tích hợp: Chỉ số thực vật chuẩn hóa $\text{NDVI} = (\text{NIR} - \text{RED}) / (\text{NIR} + \text{RED})$, Tỷ số thực vật $\text{RVI} = \text{NIR} / \text{RED}$, Tỷ số phản xạ $\text{TRRI}$, Chỉ số thực vật môi trường $\text{DVI} = \text{NIR} - \text{RED}$, cùng các biến cấu trúc bề mặt ảnh ma trận đồng xuất hiện mức xám (GLCM Homogeneity, Contrast, Entropy) và mô hình số độ cao (DEM).
- Đo đạc trắc lượng lâm học ngoại nghiệp: Đường kính ngang ngực ($D_{1,3}$) đo bằng thước vanh chuyên dụng (độ chính xác 0,1 cm); Chiều cao vút ngọn ($H_{vn}$) đo bằng thước đo cao Blume-Leiss (độ chính xác 0,5 m); Đường kính tán lá ($D_t$) đo 2 hướng trực giao Đông - Tây và Nam - Bắc bằng thước dây (độ chính xác 1,0 dm).
- Quy trình điều tra xã hội học lâm nghiệp: Triển khai công cụ PRA phỏng vấn sâu toàn diện 06 cán bộ lãnh đạo/chuyên gia kỹ thuật tại Chi cục Kiểm lâm, Chi cục Lâm nghiệp và các Ban quản lý RPH; khảo sát bảng hỏi cấu trúc đối với 50 hộ gia đình sống phụ thuộc vào rừng tại hai địa bàn đại diện: BQL RPH Đak Ang và BQL RPH Tu Mơ Rông.
- Tam giác đạc (Triangulation): Kết hợp chéo dữ liệu đo đếm lâm phần thực nghiệm, dữ liệu viễn thám GIS, số liệu khí tượng thủy văn đo dòng chảy/xói mòn và dữ liệu phỏng vấn dân tộc học.
+-------------------------------------------------------------+
| SƠ ĐỒ THIẾT KẾ Ô TIÊU CHUẨN ĐIỀU TRA |
+-------------------------------------------------------------+
| |
| [Ô dạng bản 1: 16m2] [Ô dạng bản 2: 16m2] |
| +-------+ +-------+ |
| | 4m x 4m| | 4m x 4m| |
| +-------+ +-------+ |
| |
| [Ô dạng bản 5: 16m2] |
| +-------+ |
| | 4m x 4m| (Tâm OTC) |
| +-------+ |
| |
| [Ô dạng bản 3: 16m2] [Ô dạng bản 4: 16m2] |
| +-------+ +-------+ |
| | 4m x 4m| | 4m x 4m| |
| +-------+ +-------+ |
| |
| <----------------------- 50 mét ------------------------> |
+-------------------------------------------------------------+
* Toàn bộ OTC: 50m x 50m = 2.500 m2 (gồm 25 ô thứ cấp 10m x 10m)
Data và phân tích
Xử lý và phân tích dữ liệu thống kê sinh học nâng cao trên các phần mềm chuyên dụng SPSS, SAS và R:
- Kiểm định phân bố cấu trúc: Sử dụng các hàm phân bố xác suất liên tục Weibull 3 tham số và phân bố thực nghiệm Meyer để kiểm định quy luật phân bố số cây theo chiều cao ($N/H$) và cỡ kính ($N/D$).
- Đo lường đa dạng sinh học: Sử dụng họ chỉ số đa dạng cấu trúc Renyi ($H_\alpha$) để đánh giá toàn diện độ phong phú loài, tính đồng đều và mức độ chiếm ưu thế của các loài thực vật trong từng trạng thái rừng LRTX.
- Phân tích hồi quy phi tuyến đa biến: Xây dựng phương trình tương quan bề mặt biểu diễn lượng dòng chảy mặt ($DCM$) và lượng đất xói mòn ($XM$) theo hai trục biến số độc lập: Chỉ số cấu trúc rừng ($C_s$) và Tỷ số lập địa $(Hd/s)$. Độ tin cậy của các mô hình hồi quy đạt mức ý nghĩa thống kê cao ($p < 0,001$, hệ số xác định $R^2 > 0,85$).
Phát hiện đột phá và implications
Những phát hiện then chốt
- Thực trạng suy thoái cấu trúc và trữ lượng trên diện rộng: Toàn lưu vực sông Pô Kô có 64.052 ha RPHĐN nhưng cơ cấu chất lượng rừng bị suy giảm nghiêm trọng. Rừng LRTX giàu ($M > 200,\text{m}^3/\text{ha}$) chỉ chiếm tỷ lệ nhỏ; phần lớn diện tích thuộc trạng thái rừng nghèo kiệt ($M < 50,\text{m}^3/\text{ha}$) và rừng phục hồi. Mật độ tầng cây cao giảm mạnh, cấu trúc $N/H$ ở các trạng thái LRTX nghèo không còn tuân theo quy luật phân bố chuẩn tự nhiên, xuất hiện nhiều khoảng trống tầng tán.
- Quy luật động thái tích lũy sinh khối và thảm mục: Sinh khối tươi bình quân của các trạng thái rừng biến động mạnh: Rừng LRTX giàu đạt sinh khối lớn nhất với lượng rơi rụng thảm mục duy trì ổn định, trong khi rừng nghèo và tre nứa có lớp thảm mục mỏng ($<1,5,\text{cm}$), khả năng giữ ẩm bị suy giảm hơn 60% so với rừng nguyên sinh.
- Mô hình hóa định lượng chức năng chống xói mòn và điều tiết dòng chảy: Kết quả thực nghiệm xác lập phương trình tương quan phi tuyến: Dòng chảy mặt ($DCM$) và lượng đất xói mòn ($XM$) tỷ lệ nghịch phi tuyến với chỉ số cấu trúc $C_s$ và tỷ số $Hd/s$. Khi chỉ số cấu trúc $C_s < 1,3$, dòng chảy mặt tăng đột biến từ $150,\text{m}^3/\text{ha}/\text{năm}$ lên trên $650,\text{m}^3/\text{ha}/\text{năm}$, lượng đất xói mòn vượt ngưỡng báo động $15 - 28,\text{tấn}/\text{ha}/\text{năm}$. Khi $C_s \ge 1,7$, dòng chảy mặt bị triệt tiêu xuống dưới $5%$ tổng lượng giáng thủy và xói mòn đất được kiểm soát ở mức an toàn ($<1,5,\text{tấn}/\text{ha}/\text{năm}$).
- Phân hóa đa dạng sinh học theo chỉ số Renyi: Đường cong Renyi chỉ ra rằng trạng thái rừng LRTX giàu và trung bình có mức độ đa dạng sinh học cao nhất, các loài cây mục đích bản địa như Dipterocarpus, Shorea, Litsea giữ vai trò ưu thế sinh thái; trong khi trạng thái phục hồi bị xâm lấn bởi các loài cây tiên phong ưa sáng kém giá trị phòng hộ.
- Xung đột áp lực sinh kế và chuyển đổi sử dụng đất: Phân tích thực trạng cho thấy trong giai đoạn 2010 - 2015, trữ lượng rừng vùng Tây Nguyên sụt giảm hơn 57 triệu $\text{m}^3$ (từ 327 triệu $\text{m}^3$ xuống 270 triệu $\text{m}^3$, diện tích giảm 6,1%). Tại lưu vực Pô Kô, 45% các vụ xâm hại rừng xuất phát từ nhu cầu lấn chiếm đất trồng cây công nghiệp (cà phê, cao su, sắn nương dốc) và áp lực sinh kế của hơn 40 dân tộc thiểu số địa phương.
+-------------------------------------------------------------+
| BIẾN THIÊN LƯỢNG ĐẤT XÓI MÒN THEO CHỈ SỐ CẤU TRÚC RỪNG |
+-------------------------------------------------------------+
Xói mòn (tấn/ha/năm)
30 | * (Rừng suy thoái kiệt, Cs < 0,9)
| \
20 | \
| * (Rừng nghèo, Cs = 1,0 - 1,3)
10 | \
| \
0 +-------*-----------------------*-------------------> Cs
1,3 1,7 (Ngưỡng phòng hộ an toàn)
Implications đa chiều
- Về mặt lý thuyết: Chứng minh tính tất yếu của việc tích hợp đồng thời các chỉ tiêu sinh thái lâm phần và chỉ số thủy văn trong việc lượng hóa chất lượng rừng phòng hộ, vượt qua cách tiếp cận đơn biến theo tỷ lệ che phủ truyền thống.
- Về mặt kỹ thuật lâm sinh: Cung cấp ma trận giải pháp kỹ thuật lâm sinh cụ thể cho từng trạng thái:
- Đối với rừng suy thoái kiệt ($C_s < 1,3$): Áp dụng khoanh nuôi xúc tiến tái sinh kết hợp làm giàu rừng theo rạch hoặc theo đám ($200 - 400,\text{m}^2$), mật độ trồng bổ sung $600 - 1.000,\text{cây}/\text{ha}$ bằng các loài bản địa đa tác dụng: Ươi (Scaphium macropodum), Xoay (Dialium cochinchinense), Bời lời (Litsea glutinosa), Sao đen (Hopea odorata).
- Đối với rừng trung bình và phục hồi ($1,3 \le C_s < 1,7$): Khoanh nuôi bảo vệ nghiêm ngặt, tỉa thưa nuôi dưỡng, chặt bỏ dây leo, cây sâu bệnh, xúc tiến tái sinh tự nhiên để đưa cấu trúc $C_s$ vượt ngưỡng 1,7.
- Về mặt chính sách và quản trị: Hoàn thiện cơ chế giao khoán rừng gắn với chính sách Chi trả Dịch vụ Môi trường Rừng (DVMTR), phân bổ trực tiếp nguồn thu từ các nhà máy thủy điện (như thủy điện Yaly, Plei Krông) cho các cộng đồng nhận khoán bảo vệ rừng đầu nguồn; xây dựng mô hình nông lâm kết hợp bền vững trên đất dốc để triệt tiêu nguyên nhân phá rừng làm nương rẫy.
Limitations và Future Research
Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:
- Giới hạn thời gian quan trắc chuỗi số liệu thủy văn: Số liệu đo đạc dòng chảy mặt và xói mòn đất chủ yếu tập trung trong chu kỳ nghiên cứu ngắn hạn (2012 - 2016), chưa phản ánh trọn vẹn tác động của các hiện tượng khí hậu cực đoan đa niên như El Niño hay La Niña kéo dài.
- Giới hạn không gian đại diện lập địa: Các ô tiêu chuẩn thực nghiệm tập trung chủ yếu tại BQL RPH Đak Ang và Tu Mơ Rông, chưa bao phủ toàn bộ các tiểu vùng vi khí hậu có điều kiện địa chất chuyên biệt ở vùng giáp biên giới ngã ba Đông Dương.
- Chưa mô hình hóa dòng chảy ngầm sâu: Nghiên cứu tập trung lượng hóa dòng chảy mặt ($DCM$) và xói mòn ($XM$), chưa có thiết bị chuyên dụng đo lường động thái bốc thoát hơi nước tầng tán (eddy covariance) và vận tốc truyền dẫn nước ngầm tầng sâu.
Định hướng nghiên cứu tương lai:
- Thiết lập hệ thống trạm quan trắc thủy văn rừng tự động cố định dài hạn (10 - 20 năm) trên lưu vực sông Pô Kô.
- Ứng dụng cảm biến viễn thám siêu phổ (Hyperspectral) và LiDAR gắn trên thiết bị bay không người lái (UAV) để quét 3D cấu trúc sinh khối và tầng tán rừng phòng hộ theo thời gian thực.
- Mở rộng lượng giá kinh tế tổng hợp các dịch vụ hệ sinh thái rừng (carbon sequestration, bảo tồn nguồn gen, hấp thụ bùn cát lòng hồ thủy điện) phục vụ thị trường tín chỉ carbon quốc tế.
Tác động và ảnh hưởng
- Tác động học thuật (Academic Impact): Đặt nền móng cho phương pháp lượng hóa cấu trúc rừng phòng hộ đầu nguồn tại Tây Nguyên; dự kiến trở thành tài liệu tham khảo chuẩn mực với khả năng trích dẫn cao trong các nghiên cứu lâm sinh học nhiệt đới, thủy văn môi trường và sinh thái học cảnh quan.
- Chuyển đổi thực tiễn ngành Lâm nghiệp (Industry Transformation): Cung cấp tài liệu hướng dẫn kỹ thuật phục hồi rừng đầu nguồn cho 100% các Ban quản lý rừng phòng hộ, Công ty Lâm nghiệp trên địa bàn tỉnh Kon Tum và khu vực phụ lưu sông Sê San.
- Tác động chính sách công (Policy Influence): Là luận cứ khoa học trực tiếp giúp UBND tỉnh Kon Tum, Sở NN&PTNT và Quỹ Bảo vệ & Phát triển rừng điều chỉnh định mức kinh tế kỹ thuật trong công tác trồng rừng thay thế, khoanh nuôi tái sinh và chính sách chi trả dịch vụ môi trường rừng giai đoạn 2021 - 2030.
- Lợi ích kinh tế - xã hội: Góp phần bảo vệ tuổi thọ công trình thủy điện Yaly, giảm chi phí nạo vét bồi lắng lòng hồ hàng tỷ đồng mỗi năm; nâng cao thu nhập ổn định cho hàng nghìn hộ đồng bào dân tộc thiểu số bản địa thông qua cơ chế nhận khoán và thu hoạch lâm sản ngoài gỗ bền vững.
Đối tượng hưởng lợi
- Nghiên cứu sinh, học viên cao học chuyên ngành Lâm sinh, Quản lý tài nguyên rừng: Tiếp cận phương pháp luận thiết kế ô tiêu chuẩn đa tầng, thuật toán giải đoán OBIA và các mô hình toán học dự báo xói mòn thủy văn ($C_s, Hd/s$).
- Các nhà khoa học lâm nghiệp và chuyên gia sinh thái: Kế thừa tập dữ liệu lâm học chi tiết về các trạng thái rừng lá rộng thường xanh Tây Nguyên và hệ thống phân bố phi tuyến Weibull/Meyer.
- Các Ban quản lý rừng phòng hộ, Vườn quốc gia, Khu bảo tồn thiên nhiên: Ứng dụng trực tiếp bộ tiêu chuẩn kỹ thuật phân cấp mức độ suy thoái và quy trình làm giàu rừng bằng cây bản địa.
- Các nhà hoạch định chính sách cấp tỉnh và Bộ NN&PTNT: Căn cứ ban hành chính sách giao đất lâm nghiệp gắn với sinh kế bền vững và thể chế chi trả dịch vụ môi trường rừng.
Câu hỏi chuyên sâu
1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?
Luận án đã mở rộng Lý thuyết Động thái Cấu trúc Lâm phần (Forest Stand Dynamics Theory) thông qua việc phát triển chỉ số cấu trúc tổng hợp liên quan đến phòng hộ $C_s = [(TC/H) + CP + TM]$. Luận án chứng minh rằng năng lực phòng hộ đầu nguồn của rừng nhiệt đới không chỉ quyết định bởi độ tàn che tầng cao ($TC$) mà là kết quả hiệp đồng của cấu trúc không gian 3 tầng (tầng cây cao, thảm tươi cây bụi, thảm mục hữu cơ) tương tác với tỷ số lập địa $Hd/s$.
2. Phương pháp nghiên cứu của luận án có điểm gì đổi mới so với các nghiên cứu trước đây?
So với các nghiên cứu thủy văn rừng truyền thống chỉ dừng lại ở đo đếm định tính hoặc ước lượng viễn thám đơn giản (như Bùi Ngạnh, 1977; Nguyễn Quang Mỹ, 1984), luận án đã kết hợp kỹ thuật giải đoán ảnh hướng đối tượng OBIA đa biến (NDVI, RVI, GLCM) trên eCognition với hệ thống 100 OTC sơ cấp $2.500,\text{m}^2$ đồng bộ vi cấu trúc (2.500 ô thứ cấp, 500 ô dạng bản) và các mô hình hồi quy phi tuyến đa biến kiểm định bằng phân bố Weibull/Meyer và chỉ số Renyi.
3. Phát hiện nào trong luận án gây bất ngờ nhất dưới góc độ thực nghiệm?
Phát hiện rằng thảm mục rừng ($TM$) và thảm tươi ($CP$) đóng góp tới hơn 60% vào việc giảm thiểu dòng chảy mặt và triệt tiêu xói mòn đất ở giai đoạn đầu mùa mưa. Ngay cả khi tầng cây cao có độ tàn che lớn ($TC > 0,7$), nhưng nếu lớp thảm mục và thảm tươi bị phá hủy do chăn thả hoặc đốt dọn, xói mòn đất vẫn xảy ra nghiêm trọng do hiện tượng động năng giọt nước mưa rơi tự do tích tụ từ tầng tán cao ($H > 20,\text{m}$) tạo xói mòn hạt mưa (splash erosion) mạnh hơn cả đất trống không có cây.
4. Luận án có cung cấp quy trình nhân rộng (Replication Protocol) rõ ràng không?
Có. Luận án xây dựng quy trình 5 bước chuẩn hóa: từ xây dựng khóa ảnh giải đoán OBIA, thiết lập hệ thống ô tiêu chuẩn lồng ghép ($2.500,\text{m}^2 \to 100,\text{m}^2 \to 16,\text{m}^2$), quy trình đo trắc lượng lâm học bằng Bitterlich/Blume-Leiss, đến ma trận kỹ thuật phục hồi rừng theo 4 cấp chỉ số $C_s$ hoàn toàn có thể chuyển giao áp dụng cho toàn bộ các lưu vực sông tại Tây Nguyên và dải Trường Sơn.
5. Chương trình nghiên cứu 10 năm tiếp theo được phác thảo như thế nào?
Thiết lập mạng lưới lưu vực thực nghiệm dài hạn (Long-term Experimental Watersheds) kết hợp tháp đo bốc thoát hơi nước Eddy Covariance, quan trắc biến động sinh khối bằng UAV-LiDAR, và xây dựng nền tảng số AI quản lý giám sát tài nguyên rừng phòng hộ thời gian thực tại lưu vực sông Sê San.
Kết luận
Công trình nghiên cứu của NCS. Huỳnh Văn Chung là một luận án tiến sĩ lâm nghiệp xuất sắc, mẫu mực về phương pháp nghiên cứu định lượng và có giá trị thực tiễn to lớn. Sáu đóng góp then chốt của luận án được đúc kết:
- Xác lập hệ thống cơ sở dữ liệu viễn thám - GIS cập nhật toàn diện hiện trạng 164.685 ha rừng và 64.052 ha RPHĐN lưu vực sông Pô Kô, tỉnh Kon Tum.
- Phát triển và chuẩn hóa chỉ số cấu trúc lâm phần tổng hợp $C_s = [(TC/H) + CP + TM]$ làm thước đo định lượng khả năng phòng hộ đầu nguồn.
- Xác lập phương trình tương quan phi tuyến giữa cấu trúc rừng ($C_s$), điều kiện lập địa ($Hd/s$) với dòng chảy mặt ($DCM$) và xói mòn đất ($XM$), tìm ra ngưỡng phòng hộ an toàn bền vững ($C_s \ge 1,7$).
- Làm sáng tỏ quy luật cấu trúc tầng tán $N/H, N/D$ theo phân bố Weibull/Meyer và đặc trưng đa dạng sinh học theo đường cong Renyi của các trạng thái rừng tự nhiên.
- Đề xuất hệ thống giải pháp kỹ thuật lâm sinh phục hồi rừng chuyên biệt cho từng cấp độ suy thoái, ưu tiên bảo tồn và làm giàu bằng các loài cây bản địa quý (Ươi, Xoay, Bời lời, Sao đen).
- Hoàn thiện khung giải pháp quản lý rừng bền vững đa chức năng dựa vào cộng đồng, kết hợp hoàn thiện chính sách chi trả dịch vụ môi trường rừng (PFES) tại lưu vực sông Pô Kô.
