Tổng quan luận án

Luận án tiến sĩ "Nghiên cứu mối liên hệ giữa thảm thực vật rừng với dòng chảy và chất lượng nước trong lưu vực sông Đăkbla, tỉnh Kon Tum" được tác giả Nguyễn Tấn Liêm thực hiện tại Trường Đại học Lâm nghiệp vào năm 2014, thuộc chuyên ngành Lâm sinh (mã số 62.05), dưới sự hướng dẫn khoa học của GS. Nguyễn Ngọc Lung và PGS. Phùng Văn Khoa.

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

Rừng và nước là hai tài nguyên giữ vai trò cấu thành cơ bản của cảnh quan địa lý, có mối quan hệ tương tác phức tạp với khí hậu, địa hình, thổ nhưỡng và hoạt động nhân sinh. Sông Đăkbla là một trong hai nhánh chính của hệ thống sông Sê San thuộc địa bàn hai tỉnh Kon Tum và Gia Lai, sở hữu diện tích lưu vực 3.060 km². Lưu vực này đóng vai trò cung cấp nguồn nước sinh hoạt, sản xuất cho cư dân địa phương và là nguồn nước chủ yếu vận hành các công trình thủy điện quốc gia trên hệ thống sông Sê San gồm thủy điện Yaly, Sê San 3, Sê San 3A và Sê San 4. Địa hình lưu vực sông Đăkbla có đặc trưng cao, dốc, chia cắt mạnh, chế độ mưa lớn nhưng phân bố không đều theo mùa, dẫn đến nguy cơ hạn hán trong mùa khô, lũ lụt trong mùa mưa, cùng hiểm họa xói mòn, sạt lở đất gây bồi lắng các lòng hồ thủy điện.

Văn bản luận án khẳng định khoảng trống nghiên cứu thực tế: trước thời điểm đề tài triển khai, tại khu vực lưu vực sông Đăkbla vẫn chưa có công trình nghiên cứu khoa học nào được tiến hành một cách sâu rộng, có hệ thống nhằm làm sáng tỏ quy luật tương tác giữa thảm thực vật rừng với dòng chảy và chất lượng nước. Do đó, việc xây dựng cơ sở lý luận và thực tiễn để phục vụ công tác quản lý, sử dụng bền vững hai nguồn tài nguyên này tại địa phương là yêu cầu cấp thiết.

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

Mục tiêu nghiên cứu của luận án được văn bản phân định thành hai cấp độ cụ thể:

  1. Mục tiêu nghiên cứu chung: Phân tích và đánh giá khả năng điều tiết nước trong dòng chảy sông suối và hạn chế xói mòn đất của thảm thực vật rừng, bước đầu đề xuất quy mô diện tích và chất lượng rừng cần thiết nhằm nâng cao sản lượng và chất lượng nước trong dòng chảy sông Đăkbla.
  2. Mục tiêu nghiên cứu cụ thể: Xây dựng các mô hình toán mô phỏng mối liên hệ giữa thảm thực vật rừng và các yếu tố lập địa với dòng chảy và chất lượng nước trong lưu vực sông Đăkbla làm cơ sở tính toán diện tích rừng cần thiết trong lưu vực.

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

  • Đối tượng nghiên cứu: Các mối liên hệ giữa đặc trưng của thảm thực vật rừng (bao gồm độ che phủ rừng - $CP$ và tỷ lệ diện tích rừng giàu và trung bình - $P_{gtb}$) cùng các nhân tố lập địa (lượng mưa, độ dốc, chỉ số hình dạng lưu vực) với các đặc trưng dòng chảy (lưu lượng dòng chảy - $Q$, mô đuyn dòng chảy - $M$, hệ số dòng chảy - $\eta$) và chất lượng nước (được đại diện bởi mô đuyn bùn cát lơ lửng mùa lũ - $M_{bclu}$).
  • Phạm vi không gian: Toàn bộ phạm vi lưu vực sông Đăkbla thuộc địa bàn tỉnh Kon Tum với quy mô diện tích tự nhiên là 3.060 km² cùng các lưu vực/tiểu lưu vực hợp thành.
  • Phạm vi thời gian và tư liệu: Dữ liệu quan trắc khí tượng, thủy văn và hiện trạng rừng tập trung phân tích thực nghiệm trong giai đoạn từ năm 2011 đến năm 2013, kết hợp chuỗi dữ liệu bản đồ và không gian viễn thám tại thời điểm nghiên cứu.

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

Các hướng nghiên cứu trên thế giới

Tác giả hệ thống hóa các hướng nghiên cứu thủy văn rừng trên thế giới qua ba quy mô tiếp cận chính:

  1. Quy mô ô thí nghiệm: Tập trung trên diện tích nhỏ nhằm nhận diện cơ chế của từng quá trình vật lý thủy văn cơ bản, làm tiền đề ngoại suy lên quy mô lớn hơn.
  2. Quy mô khu rừng: Vận dụng nguyên lý và phương trình cân bằng nước để theo dõi lượng giáng thủy qua các thành phần: lượng nước giữ lại trên tán và bốc hơi, lượng nước chảy men thân cây, lượng nước lọt qua tán xuống đất; trên mặt đất phân chia thành lượng hấp thụ qua thảm mục, bốc hơi mặt đất, dòng chảy tràn bề mặt, dòng chảy ngấm vào đất hình thành nước ngầm nông hoặc ngầm sâu, và lượng thoát hơi qua rễ cây.
  3. Quy mô lưu vực: Tiếp cận theo hai phương pháp gồm nghiên cứu lưu vực cặp đôi (paired watershed, áp dụng diện tích $< 10\text{ km}^2$ theo Zhang và cộng sự, 2004) và nghiên cứu từng lưu vực riêng rẽ (single watershed, áp dụng cho lưu vực quy mô lớn).

Về mô hình hóa dòng chảy lưu vực, luận án phân tích các nhóm mô hình phân bố không gian tự nhiên và mô hình tổng cộng (như BASINS, SWAT, PERFECT, SoilWat, WMS). Tác giả dẫn chứng tài liệu của FAO và CIFOR (2005) về lịch sử "lý thuyết miếng bọt biển" hình thành từ thế kỷ XIX (được Myers nhắc lại năm 1986) cho rằng đất rừng tích nước mùa mưa để nhả ra mùa khô, nhưng lý thuyết này bị tranh luận từ năm 1920.

Các nghiên cứu thực nghiệm được tác giả điểm lại cho thấy nhiều nhận định trái chiều:

  • Bosch (1982) tổng kết 94 lưu vực chỉ ra: biến đổi 10% độ che phủ rừng thông và bạch đàn làm thay đổi trung bình 40 mm sản lượng nước năm; rừng gỗ cứng lá rộng rụng lá biến đổi 10% độ che phủ làm thay đổi 25 mm; rừng cây bụi và thảm cỏ thay đổi 10 mm.
  • Bruijnzeel (2004) chỉ ra thảm thực vật rừng giảm dòng chảy mặt và lọc sạch nước, nhưng hoạt động trồng rừng lại có thể mang tác động ngược lại.
  • Vidan Sahin và Michael Hall (1996) phân tích 145 thí nghiệm cho thấy giảm 10% che phủ rừng lá kim làm tăng 20–25 mm sản lượng nước, rừng rụng lá tăng 17–19 mm, rừng bạch đàn tăng 6 mm; ngược lại phục hồi 10% che phủ cây bụi làm giảm 5 mm sản lượng nước.
  • Kathleen A. Farley và cộng sự (2005) chứng minh chuyển đổi trảng cỏ và cây bụi sang rừng trồng làm giảm sản lượng dòng chảy năm 44% và 31% (bạch đàn giảm 75%, thông giảm 40%). Ge Sun và cộng sự (2005) ghi nhận trồng rừng ở vùng đất bán khô hạn và nhiệt đới thoái hóa tại Trung Quốc làm giảm dòng chảy từ 50 đến 300 mm/năm (giảm tương đối 20–40%). P. Versfeld (2007) ghi nhận trồng rừng tại Nam Phi giảm 3,2% dòng chảy năm và 7,8% dòng chảy kiệt. Egginton tổng kết hơn 500 lưu vực toàn cầu cho thấy trên 70% trường hợp trồng rừng trên đồng cỏ, cây bụi làm giảm 38% dòng chảy năm.
  • Ngược lại, Bishop và cộng sự nghiên cứu tại lưu vực Koga (266 km², Ethiopia) ghi nhận khi che phủ rừng giảm từ 16% (năm 1957) xuống 1% (năm 1986), chuỗi thủy văn từ năm 1960 đến 2002 không ghi nhận biến đổi rõ rệt trong chế độ dòng chảy.
  • Về lũ lụt, báo cáo của FAO và CIFOR (2005) dẫn các nghiên cứu của Gilmour (1987), Hamilton (1987), Bruijnzeel (1990, 2004), Calder (2000), Hofer và Messerli (2006) kết luận: rừng có ảnh hưởng rất hạn chế đối với các trận lũ lụt lớn trên diện rộng ở hạ lưu, do khi mưa lớn kéo dài, đất rừng đã bão hòa nước thì dòng chảy tràn vẫn phát sinh.
  • Về chất lượng nước, thí nghiệm tại lưu vực Houay Pano (Lào) ghi nhận chất bồi lắng từ nương rẫy trồng lúa cao gấp 3 lần so với rừng tre nứa và gấp 9 lần so với bãi cỏ; Van Dijk và cộng sự (2007), Zhang và cộng sự (2007) chỉ ra mất rừng làm gia tăng nồng độ muối trong dòng chảy; Phùng Văn Khoa (2006) chỉ ra lượng mưa, đá mẹ và thảm thực vật rừng chi phối rõ rệt sản lượng ion trong nước tại các lưu vực ở Hoa Kỳ.

Các hướng nghiên cứu tại Việt Nam

Tại Việt Nam, nghiên cứu thủy văn rừng phát triển từ thập niên 1970 và mở rộng từ thập niên 1990:

  • Ở quy mô khu rừng: Nguyễn Tử Siêm và Thái Phiên (1995) nghiên cứu xói mòn đất dưới canh tác nông nghiệp và rừng trồng; Võ Đại Hải (1996) nghiên cứu cấu trúc rừng phòng hộ (giảm độ tàn che từ 0,7–0,8 xuống 0,3–0,4 làm dòng chảy mặt tăng 30,4% ở rừng tự nhiên và 33,8% ở rừng le); Nguyễn Ngọc Lung và Võ Đại Hải (1997) xác định tán rừng ngăn cản 5,7–11,6% lượng mưa, đất rừng thấm 88,2–92,5% tổng lượng mưa, xây dựng bảng tra hệ số thảm thực vật $C$; Phạm Văn Điển (1998, 2006) tại hồ Hòa Bình xác định tỷ lệ giữ nước trên tán từ 2,91–18,55%, dòng chảy mặt từ 104,7–574,7 mm/ha/năm (hệ số 5,2–28,7%), thiết lập chỉ số nuôi dưỡng nguồn nước $U = (GT + CP + TM) / (K \cdot S)$; Phùng Văn Khoa (1997) nghiên cứu rừng thông đuôi ngựa (dòng chảy mặt chiếm 3–5%, thoát hơi nước chiếm 30–40%, bốc hơi đất chiếm 30–35% lượng mưa).
  • Ở quy mô lưu vực: Nguyễn Viết Phổ (1992), Vũ Văn Tuấn (1977, 1981, 1982); Ngô Đình Quế và cộng sự (2003–2004) tại Miền Trung và Tây Nguyên; Võ Đại Hải và Ngô Đình Quế (2006) tại sông Rào Nậy (giảm 11% rừng làm tăng đỉnh lũ 7,6%, giảm lưu lượng kiệt 16%) và sông Pô Kô (tăng 29% che phủ rừng làm giảm đỉnh lũ 13%); Vũ Tấn Phương, Phạm Thị Hương Lan và Nguyễn Thị Hải (2007) tại sông Chảy, sông Bồ, sông Ba (tăng 20–25% che phủ rừng làm giảm dòng chảy mùa lũ 7,5–9,3%, tăng dòng chảy mùa kiệt 8,5–11,2%, giảm xói mòn 7,2–38%); Lương Thị Vân (2001) tại sông Lại và sông Côn; Vương Văn Quỳnh và cộng sự (2007) đưa ra phương pháp rừng quy đổi ($RQĐ$), chỉ số giữ nước $K$ và công thức tính tỷ lệ $RQĐ$ cần thiết:

$$\text{Tỷ lệ } RQĐ \text{ cần thiết (%)} = \frac{4500 \cdot \Delta H_{bq}}{F_{lv}}$$

(trong đó $\Delta H_{bq}$ là độ chênh cao bình quân lưu vực tính bằng mét, $F_{lv}$ là diện tích lưu vực tính bằng hecta).

Ứng dụng mô hình: Vũ Văn Tuấn và Phạm Thị Hương Lan (2004) dùng SWAT tại Mù Cang Chải; Nguyễn Thị Bích và Nguyễn Kiên Dũng (2005) dùng SWAT tính bùn cát tại sông Sê San; Dương Văn Khánh (2001) dùng mô hình WMS tại lưu vực sông Zagodzonka (Ba Lan).

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

Luận án chỉ ra các hạn chế trong nghiên cứu trước đây: số lượng công trình ở quy mô lưu vực tại Việt Nam còn ít; thời gian nghiên cứu ngắn; các mô hình vật lý phân bố như SWAT đòi hỏi hệ thống dữ liệu chuỗi thời gian khí tượng, thủy văn và thổ nhưỡng liên tục, hoàn chỉnh — điều kiện mà trạm quan trắc tại lưu vực sông Đăkbla chưa thể đáp ứng. Hơn nữa, các nghiên cứu chưa kết hợp đồng thời ảnh hưởng của độ che phủ ($CP$) lẫn chất lượng rừng (tỷ lệ rừng giàu và trung bình - $P_{gtb}$) trong mối tương tác với các yếu tố hình thái và lập địa đến dòng chảy và bùn cát lơ lửng. Luận án lựa chọn giải quyết khoảng trống này bằng phương pháp nghiên cứu lưu vực riêng rẽ kết hợp giải đoán viễn thám, GIS và mô hình hóa thống kê đa biế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

Nghiên cứu dựa trên nguyên lý cân bằng nước lưu vực và lý thuyết thủy văn rừng. Lưu vực được tiếp cận như một hệ thống tự nhiên đón nhận giáng thủy và tập trung dòng chảy về điểm thoát chung.

  • Dòng chảy lưu vực gồm dòng chảy mặt (surface flow), dòng chảy sát mặt đất (subsurface flow) và dòng chảy ngầm (groundwater flow); phân theo mùa gồm dòng chảy năm, dòng chảy lũ và dòng chảy kiệt.
  • Chất lượng nước: Luận án tập trung đánh giá thông qua lượng bùn cát lơ lửng mùa lũ. Căn cứ theo Nguyễn Thanh Sơn (2003), dòng chảy bùn cát gồm phù sa lơ lửng, phù sa đáy và chất hòa tan, trong đó phù sa lơ lửng chiếm tỷ lệ trên 80%.
Ký hiệu Đại lượng và đơn vị đo Ý nghĩa trong luận án
$F_{lv}$ Diện tích lưu vực ($\text{km}^2$ hoặc $\text{ha}$) Quy mô không gian đón nhận giáng thủy
$K_d$ Chỉ số hình dạng lưu vực Tỷ lệ hình học ảnh hưởng tập trung dòng chảy
$S_{bq}$ Độ dốc bình quân lưu vực (độ hoặc %) Độ dốc địa hình bề mặt lưu vực
$X_{bq}, X_{năm}, X_{lũ}$ Lượng mưa bình quân, năm, mùa lũ ($\text{mm}$) Lượng nước rơi bề mặt theo chu kỳ
$T_{năm}$ Tổng lượng nước mưa năm ($\text{m}^3$ hoặc triệu $\text{m}^3$) Nguồn cấp nước đầu vào của lưu vực
$CP$ Độ che phủ rừng (%) Tỷ lệ phần trăm diện tích có rừng che phủ
$P_{gtb}$ Tỷ lệ diện tích rừng giàu và trung bình (%) Phản ánh cấu trúc và chất lượng thảm rừng
$Q, Q_{năm}$ Lưu lượng dòng chảy bình quân, bình quân năm ($\text{m}^3/\text{s}$) Thể tích nước chảy qua mặt cắt trong 1 giây
$M, M_{năm}, M_{lũ}$ Mô đuyn dòng chảy bình quân năm, mùa lũ ($\text{l}/\text{s}\cdot\text{km}^2$) Suất sinh dòng chảy trên một đơn vị diện tích
$\eta, \eta_{năm}$ Hệ số dòng chảy năm Tỷ số giữa lớp dòng chảy và lượng mưa
$M_{bclu}$ Mô đuyn bùn cát lơ lửng mùa lũ ($\text{tấn}/\text{km}^2$ hoặc $\text{g}/\text{s}\cdot\text{km}^2$) Mức độ xói mòn và chất lượng nước mùa lũ

Phương pháp nghiên cứu thực tế

  • Phương pháp thu thập thông tin và dữ liệu không gian:
    • Ứng dụng công nghệ viễn thám: Sử dụng ảnh vệ tinh Landsat đa thời kỳ, xử lý trên phần mềm ENVI để giải đoán hiện trạng rừng và phân loại các trạng thái rừng.
    • Ứng dụng hệ thống thông tin địa lý (GIS): Sử dụng phần mềm ArcGIS, số hóa bản đồ địa hình, xây dựng mô hình số hóa độ cao DEM (Digital Elevation Model) nhằm trích xuất ranh giới lưu vực, diện tích ($F_{lv}$), chiều dài ($L$), chiều rộng bình quân ($B$), độ dốc bình quân ($S_{bq}$), mật độ lưới sông suối và chỉ số hình dạng lưu vực ($K_d$).
    • Thu thập dữ liệu thứ cấp: Tiếp nhận số liệu quan trắc khí tượng (lượng mưa) và thủy văn (lưu lượng, mô đuyn dòng chảy, bùn cát lơ lửng) giai đoạn 2011–2013 từ Trung tâm Khí tượng Thủy văn tỉnh Kon Tum; số liệu hiện trạng và phân loại rừng từ Chi cục Lâm nghiệp, Chi cục Kiểm lâm tỉnh Kon Tum và Công ty tư vấn lâm nghiệp Chinh Nguyên.
  • Phương pháp phân tích và xử lý thống kê:
    • Sử dụng phần mềm SPSS (Statistical Package for the Social Sciences) thực hiện phân tích tương quan và hồi quy đa biến.
    • Khảo sát các dạng hàm quan hệ toán học: Tuyến tính (Linear), Parabol bậc hai (Quadratic), Parabol bậc ba (Cubic), Logarit (Logarithmic), Nghịch đảo (Inverse), Hàm lũy thừa (Power), Hàm phức hợp (Compound), Hàm S dạng:

$$\ln Y = B_0 + \frac{B_1}{X} \quad \text{hoặc} \quad Y = \exp\left(B_0 + \frac{B_1}{X}\right)$$

  • Đánh giá mức độ phù hợp và ý nghĩa thống kê của các mô hình bằng phân tích phương sai ANOVA thông qua các đại lượng: Tổng bình phương hồi quy ($SS_{Regression}$), Tổng bình phương sai số ($SS_{Residual}$), Tổng bình phương chung ($SS_{Total}$), hệ số tương quan ($R$), hệ số xác định ($R^2$) và mức ý nghĩa ($p\text{-value} / \text{Sig.}$).

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

Mở đầu

Trình bày bối cảnh chung về vị trí chiến lược của sông Đăkbla đối với phát triển kinh tế, cấp nước sinh hoạt và an ninh năng lượng cho chuỗi thủy điện sông Sê San. Xác định mục tiêu tổng quát, mục tiêu cụ thể, đối tượng nghiên cứu, ý nghĩa khoa học trong việc bổ sung luận cứ cho thủy văn rừng Tây Nguyên và ý nghĩa thực tiễn đối với công tác quy hoạch, quản lý rừng phòng hộ đầu nguồn tại tỉnh Kon Tum.

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

Hệ thống hóa toàn diện các thành tựu và tranh luận học thuật về thủy văn rừng trên phạm vi quốc tế và Việt Nam. Điểm lại các nghiên cứu ở ba cấp độ (ô thí nghiệm, khu rừng, lưu vực), phân tích các quan điểm trái chiều về ảnh hưởng của độ che phủ rừng đến sản lượng nước, dòng chảy kiệt, dòng chảy lũ và xói mòn đất. Luận giải chi tiết lý do không áp dụng mô hình thủy văn SWAT cho lưu vực sông Đăkbla do hạn chế về mật độ trạm quan trắc và chuỗi số liệu đất đai, khí văn đầu vào, từ đó xác lập tính tất yếu của việc lựa chọn phương pháp nghiên cứu thống kê kết hợp viễn thám và GIS.

Chương 2: Nội dung và phương pháp nghiên cứu

Trình bày cấu trúc ba nội dung nghiên cứu lớn:

  1. Phân tích các đặc trưng cơ bản về địa hình, hình thái, thổ nhưỡng, khí hậu và thảm thực vật rừng tại các lưu vực nghiên cứu.
  2. Thiết lập và phân tích mô hình liên hệ giữa thảm thực vật rừng (độ che phủ $CP$, tỷ lệ rừng giàu và trung bình $P_{gtb}$) cùng các yếu tố lập địa với các đại lượng dòng chảy ($Q_{năm}$, $M_{năm}$, $M_{lũ}$, $\eta_{năm}$) và chất lượng nước ($M_{bclu}$).
  3. Đề xuất quy mô diện tích và chất lượng rừng cần thiết nhằm nâng cao hiệu quả điều tiết nước và giảm thiểu xói mòn.

Quy định chi tiết các bước giải đoán ảnh Landsat, quy trình xử lý không gian bằng ArcGIS và ENVI, công thức xác định các thông số hình thái lưu vực, phương pháp lấy mẫu đo đạc thủy văn và quy trình phân tích hồi quy, kiểm định phương sai bằng phần mềm SPSS.

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

Chương này tập trung toàn bộ các kết quả phân tích số liệu thực tế giai đoạn 2011–2013 được trình bày qua hệ thống bảng biểu và đồ thị:

  • Đặc trưng địa mạo, thổ nhưỡng và thảm thực vật: Xác định vị trí, ranh giới 3.060 km² của sông Đăkbla trong hệ thống sông Kon Tum; phân tích diện tích, độ dốc, chỉ số hình dạng các lưu vực; diện tích các loại đất; thống kê lượng mưa 2011–2013; phân loại diện tích các trạng thái rừng tự nhiên (giàu, trung bình, nghèo), rừng trồng, đất trống; xác định độ che phủ rừng ($CP$) và tỷ lệ diện tích rừng giàu, trung bình ($P_{gtb}$).
  • Đặc trưng chế độ dòng chảy và bùn cát: Tổng hợp lưu lượng bình quân năm ($Q_{năm}$), mùa lũ ($Q_{lũ}$), mùa cạn ($Q_{cạn}$); lượng bùn cát lơ lửng mùa lũ; tổng lượng dòng chảy, mô đuyn dòng chảy năm ($M_{năm}$), mô đuyn mùa lũ ($M_{lũ}$), hệ số dòng chảy năm ($\eta_{năm}$) và mô đuyn bùn cát lơ lửng ($M_{bclu}$).
  • Mối liên hệ giữa độ che phủ rừng ($CP$) với dòng chảy và bùn cát:
    • Thiết lập biểu đồ thăm dò và bảng tương quan giữa $CP$ với $Q_{năm}$, $M_{năm}$, $M_{lũ}$, $\eta_{năm}$, $M_{bclu}$.
    • Phân tích phương sai ANOVA và xác định hệ số phương trình hồi quy đơn biến giữa $CP$ với từng biến phụ thuộc.
    • Thiết lập mô hình hồi quy đa biến: phương trình kết hợp giữa tổng lượng mưa năm và độ che phủ rừng $(T_{năm}, CP)$ với $Q_{năm}$; mô hình mở rộng $(T_{năm}, CP, S_{bq}, K_d)$ với $Q_{năm}$; mô hình giữa lượng mưa năm và độ che phủ $(X_{năm}, CP)$ với $M_{năm}$; mô hình giữa lượng mưa lũ và độ che phủ $(X_{lũ}, CP)$ với $M_{lũ}$; mô hình giữa $(X_{lũ}, CP)$ với $M_{bclu}$.
  • Mối liên hệ giữa chất lượng rừng ($P_{gtb}$) với dòng chảy:
    • Thăm dò và thiết lập phương trình hồi quy đơn biến giữa $P_{gtb}$ với $Q_{năm}$, $M_{năm}$, $M_{lũ}$.
    • Phân tích phương sai và lập mô hình hồi quy đa biến: giữa $(T_{năm}, P_{gtb})$ với $Q_{năm}$; giữa $(X_{năm}, P_{gtb})$ với $M_{năm}$; giữa $(X_{lũ}, P_{gtb})$ với $M_{lũ}$.
  • Mô phỏng biến động và đề xuất quy mô, chất lượng rừng:
    • Lập bảng tính toán kết quả biến động mô đuyn dòng chảy năm ($M_{năm}$) theo sự biến đổi kết hợp của lượng mưa bình quân năm ($X_{năm}$) và độ che phủ rừng ($CP$).
    • Lập bảng biến động mô đuyn bùn cát lơ lửng mùa lũ ($M_{bclu}$) theo sự thay đổi của lượng mưa lũ ($X_{lũ}$) và độ che phủ rừng ($CP$).
    • Lập bảng biến động mô đuyn dòng chảy năm ($M_{năm}$) theo sự biến đổi của lượng mưa và tỷ lệ diện tích rừng giàu, trung bình ($P_{gtb}$).
    • Đánh giá năng lực điều tiết nước và chống xói mòn của hiện trạng rừng, từ đó lượng hóa mức độ che phủ tối thiểu và tỷ lệ chất lượng rừng cần duy trì cho lưu vực.

Kết luận và kiến nghị

Tổng kết các phát hiện định lượng rút ra từ các mô hình toán học; khẳng định vai trò điều tiết của thảm thực vật rừng đối với chế độ dòng chảy và bùn cát lơ lửng trên lưu vực sông Đăkbla; đưa ra các kiến nghị cụ thể cho cơ quan quản lý lâm nghiệp tỉnh Kon Tum trong công tác khoanh nuôi, bảo vệ và phục hồi chất lượng rừng đầu nguồn.

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. Lượng hóa mối quan hệ tương tác phức hợp: Luận án đã xác lập một cách định lượng mối liên hệ giữa các chỉ tiêu đặc trưng của thảm thực vật rừng (cả về số lượng là độ che phủ $CP$, lẫn chất lượng là tỷ lệ rừng giàu và trung bình $P_{gtb}$) với các yếu tố hình thái lưu vực ($S_{bq}, K_d$) và điều kiện giáng thủy ($T_{năm}, X_{năm}, X_{lũ}$) tác động đến các đặc trưng dòng chảy ($Q_{năm}, M_{năm}, M_{lũ}, \eta_{năm}$) và chất lượng nước ($M_{bclu}$) tại vùng Tây Nguyên.
  2. Xây dựng hệ thống mô hình toán học thực nghiệm: Thông qua phân tích phương sai ANOVA và kiểm định thống kê trên SPSS, luận án đã xác lập hệ thống các phương trình hồi quy đơn biến và đa biến (tuyến tính và phi tuyến). Các mô hình này cung cấp công cụ toán học mô phỏng sự biến động của mô đuyn dòng chảy năm, mô đuyn dòng chảy mùa lũ và mô đuyn bùn cát lơ lửng khi có sự thay đổi của độ che phủ rừng và chất lượng rừng trong các kịch bản lượng mưa khác nhau.

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

  1. Đánh giá thực trạng rừng đầu nguồn: Cung cấp cơ sở dữ liệu số hóa hoàn chỉnh về hiện trạng phân bố các trạng thái rừng, độ dốc, đặc trưng hình thái lưới sông suối lưu vực sông Đăkbla dựa trên công nghệ GIS và viễn thám.
  2. Cơ sở xác định diện tích và chất lượng rừng phòng hộ: Các bảng tính toán biến động mô đuyn dòng chảy và mô đuyn bùn cát lơ lửng là căn cứ khoa học trực tiếp giúp các nhà quản lý lâm nghiệp và chính quyền tỉnh Kon Tum xác định quy mô diện tích rừng cần bảo vệ, định hướng tỷ lệ diện tích rừng giàu và trung bình ($P_{gtb}$) tối ưu.
  3. Bảo đảm an toàn cho các công trình thủy điện: Kết quả tính toán về sự suy giảm mô đuyn bùn cát lơ lửng khi tăng độ che phủ rừng là luận cứ thực tiễn để xây dựng chính sách chi trả dịch vụ môi trường rừng, giảm bồi lắng lòng hồ và duy trì nguồn nước ổn định cho các nhà máy thủy điện Yaly, Sê San 3, Sê San 3A và Sê San 4.

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

Hạn chế được ghi nhận từ văn bản

  • Độ dài chuỗi số liệu quan trắc: Chuỗi số liệu khí tượng và thủy văn phục vụ tính toán thực nghiệm trong luận án chỉ giới hạn trong giai đoạn 3 năm (từ năm 2011 đến năm 2013).
  • Hệ thống quan trắc và nguồn dữ liệu đầu vào: Mạng lưới trạm khí tượng thủy văn trên địa bàn lưu vực sông Đăkbla còn thưa thớt, lạc hậu; các nghiên cứu cơ bản về đất rừng còn nhiều hạn chế, dẫn đến việc thiếu hụt chuỗi dữ liệu liên tục kéo dài. Do đó, luận án không thể áp dụng các mô hình vật lý phân bố chi tiết cao như mô hình SWAT.
  • Phạm vi chỉ tiêu chất lượng nước: Luận án mới chỉ tập trung phân tích một chỉ tiêu chất lượng nước duy nhất là mô đuyn bùn cát lơ lửng trong mùa lũ ($M_{bclu}$), chưa có điều kiện khảo sát toàn diện các chỉ tiêu hóa lý, dinh dưỡng hoặc nồng độ ion hòa tan trong nước sông suối.

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

  • Mở rộng quan trắc và thu thập chuỗi dữ liệu thủy văn, khí tượng và bùn cát liên tục qua nhiều năm để kiểm định và nâng cao độ chính xác của các tham số mô hình hồi quy.
  • Tiếp tục nghiên cứu bổ sung các chỉ tiêu chất lượng nước khác (như thành phần hóa học, độ dẫn điện, độ mặn, các hợp chất dinh dưỡng $N, P$) trong mối liên hệ với các trạng thái thảm thực vật rừng.
  • Khi điều kiện mạng lưới quan trắc và cơ sở dữ liệu đất đai được hoàn thiện, tiến hành thử nghiệm ứng dụng các mô hình thủy văn phân bố không gian (như SWAT) để so sánh và đối chứng kết quả với mô hình hồi quy thực nghiệm.

Giá trị tham khảo

Luận án mang lại giá trị tham khảo học thuật và ứng dụng thực tiễn cho nhiều nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và giảng viên: Cung cấp tài liệu tham khảo chuyên sâu trong các chuyên ngành Lâm sinh, Thủy văn rừng, Quản lý tài nguyên rừng và Môi trường. Đáng tham khảo nhất là quy trình kết hợp công nghệ viễn thám Landsat (xử lý trên ENVI), hệ thống thông tin địa lý ArcGIS và phân tích thống kê đa biến trên SPSS để giải quyết bài toán thủy văn rừng ở quy mô lưu vực trong điều kiện thiếu chuỗi số liệu quan trắc dài hạn.
  • Nhà quy hoạch và cơ quan quản lý nhà nước: Sở Nông nghiệp và Phát triển nông thôn, Chi cục Kiểm lâm, Chi cục Lâm nghiệp tỉnh Kon Tum và Gia Lai có thể khai thác các bảng mô phỏng biến động $M_{năm}$ và $M_{bclu}$ để xác định phân vùng phòng hộ đầu nguồn xung yếu, lập quy hoạch diện tích rừng phòng hộ và đề xuất giải pháp khoanh nuôi nâng cao chất lượng rừng tự nhiên.
  • Đơn vị quản lý khai thác thủy điện: Các ban quản lý công trình thủy điện trên dòng sông Sê San (Yaly, Sê San 3, 3A, 4) có cơ sở khoa học để đánh giá nguy cơ bồi lắng lòng hồ theo mức độ che phủ rừng thượng nguồn, làm căn cứ thực hiện chính sách chi trả dịch vụ môi trường rừng đối với chức năng điều tiết nguồn nước và hạn chế xói mòn.

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

1. Vì sao luận án không sử dụng mô hình SWAT để mô phỏng thủy văn lưu vực sông Đăkbla?

Mô hình SWAT (Soil and Water Assessment Tools) là mô hình vật lý phân bố đòi hỏi các dữ liệu đầu vào không gian và thuộc tính rất chi tiết, bao gồm nhiều loại bản đồ số hóa cùng số liệu khí tượng, dòng chảy được quan trắc liên tục trong thời gian dài. Tại lưu vực sông Đăkbla, hệ thống quan trắc khí tượng thủy văn còn thiếu thốn và các nghiên cứu về đất rừng còn hạn chế, không đáp ứng đủ chuỗi dữ liệu đầu vào để vận hành mô hình. Do đó, tác giả lựa chọn phương pháp phân tích thống kê tương quan - hồi quy kết hợp công nghệ viễn thám và GIS thay cho SWAT.

2. Các đặc trưng của thảm thực vật rừng được lượng hóa bằng những chỉ tiêu nào trong mô hình?

Thảm thực vật rừng trong luận án được lượng hóa qua hai chỉ tiêu cơ bản:

  1. Độ che phủ rừng ($CP$): Thể hiện tỷ lệ phần trăm diện tích có rừng che phủ trên tổng diện tích lưu vực nghiên cứu.
  2. Tỷ lệ diện tích rừng giàu và trung bình ($P_{gtb}$): Phản ánh chất lượng và trạng thái thảm thực vật rừng tự nhiên có trữ lượng cao, quyết định trực tiếp đến dung tích chứa nước của đất và khả năng cản dòng chảy mặt.

3. Vì sao luận án chọn mô đuyn bùn cát lơ lửng mùa lũ để đánh giá chất lượng nước?

Dòng chảy bùn cát trong sông suối bao gồm ba thành phần: phù sa lơ lửng, phù sa đáy và vật chất hòa tan. Theo cơ sở lý thuyết được tác giả viện dẫn từ Nguyễn Thanh Sơn (2003), phù sa lơ lửng là thành phần chiếm tỷ lệ áp đảo (trên 80% tổng lượng dòng chảy bùn cát). Trong điều kiện giới hạn về nguồn lực nghiên cứu, chỉ tiêu mô đuyn bùn cát lơ lửng mùa lũ ($M_{bclu}$) phản ánh trực tiếp mức độ xói mòn đất trên bề mặt lưu vực và mức độ bồi lắng hồ chứa thủy điện, do đó được chọn làm đại lượng đặc trưng cho chất lượng nước.

4. Báo cáo của FAO và CIFOR (2005) được luận án trích dẫn nêu nhận định gì về vai trò của rừng đối với lũ lụt lớn?

Theo tổng kết của FAO và CIFOR (2005) trong luận án, rừng có ảnh hưởng rất hạn chế đối với các trận lũ lụt lớn ở vùng hạ lưu, đặc biệt là các trận lũ xảy ra trên diện rộng. Rừng và đất rừng chỉ có khả năng làm giảm dòng chảy tràn trong các trận mưa nhỏ ở quy mô cục bộ. Khi lượng mưa lớn và kéo dài, tầng đất rừng đạt trạng thái bão hòa hoàn toàn, không thể thấm hút thêm nước, dẫn đến hiện tượng nước chảy tràn trên bề mặt đất tương tự như các vùng đất không có rừng.

5. Những yếu tố lập địa và hình thái lưu vực nào được đưa vào các phương trình hồi quy đa biến?

Bên cạnh các chỉ tiêu thảm thực vật rừng ($CP, P_{gtb}$), các phương trình hồi quy đa biến trong luận án tích hợp các nhân tố: tổng lượng nước mưa năm trên lưu vực ($T_{năm}$), lượng mưa bình quân năm ($X_{năm}$), lượng mưa mùa lũ ($X_{lũ}$), độ dốc bình quân lưu vực ($S_{bq}$) và chỉ số hình dạng lưu vực ($K_d$).

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Tấn Liêm đã hoàn thành việc đánh giá thực nghiệm mối quan hệ giữa thảm thực vật rừng và các yếu tố lập địa với dòng chảy và bùn cát lơ lửng tại lưu vực sông Đăkbla giai đoạn 2011–2013. Bằng việc kết hợp công nghệ viễn thám, GIS và công cụ phân tích hồi quy đa biến, công trình đã xây dựng hệ thống các hàm toán học định lượng vai trò của độ che phủ và chất lượng rừng đối với khả năng điều tiết nguồn nước mùa lũ, dòng chảy năm và giảm thiểu xói mòn đất. Kết quả nghiên cứu là tài liệu khoa học hữu ích, đóng góp luận cứ thực tiễn cho công tác quy hoạch lâm nghiệp, bảo vệ rừng đầu nguồn và duy trì an toàn vận hành các hồ chứa thủy điện trên lưu vực sông Sê San.