Luận án tiến sĩ: Ứng dụng viễn thám và GIS trong quản lý cây công nghiệp huyện Bảo Lâm, Lâm Đồng - Lê Quang Toan

Tổng quan về luận án

Quản lý không gian cây công nghiệp lâu năm (CCNLN) trong mối tương quan với biến động tài nguyên rừng là một trong những bài toán cốt lõi của địa lý kinh tế - sinh thái và quản lý tài nguyên bền vững tại Việt Nam. Vùng Tây Nguyên chiếm tới hơn 60% diện tích đất bazan của cả nước, đóng vai trò là thủ phủ cây công nghiệp xuất khẩu chiến lược nhưng cũng là điểm nóng về suy thoái diện tích rừng tự nhiên. Tại tỉnh Lâm Đồng, giá trị sản xuất nông nghiệp năm 2014 của ngành trồng trọt chiếm tới 82,96% tổng giá trị sản xuất nông nghiệp toàn tỉnh. Trong đó, huyện Bảo Lâm có quy mô diện tích tự nhiên hơn 146.000 ha (chiếm 19% diện tích toàn tỉnh), sở hữu diện tích chè lớn nhất và diện tích cà phê lớn thứ hai trong tỉnh, là địa bàn điển hình diễn ra các chu trình chuyển đổi lớp phủ phức tạp giữa rừng tự nhiên, cây hàng năm và CCNLN.

Khoảng trống nghiên cứu (research gap) trọng yếu được xác định là: Mặc dù kỹ thuật viễn thám quang học đã phát triển mạnh, việc bóc tách và phân loại chuyên đề CCNLN (đặc biệt là cà phê và chè) bằng các thuật toán truyền thống dựa trên điểm ảnh (pixel-based) luôn gặp phải hiện tượng đồng phổ dị vật (spectral confusion) giữa tán cây cà phê chưa khép tán với cây hàng năm, hoặc giữa cà phê/chè trưởng thành với rừng tái sinh và thảm cây bụi. Các nghiên cứu quốc tế điển hình như Mukashema và nnk (2014) [102] hay Cordero và nnk (2007) [60] đã chỉ ra độ chính xác phân loại diện tích cà phê từ ảnh vệ tinh đa phổ chỉ dao động ở mức thấp (37,4% đến 50%) nếu chỉ thuần túy dựa vào giá trị phản xạ phổ. Mặt khác, các nghiên cứu trong nước trước đây chủ yếu dừng lại ở việc thống kê diện tích hoặc số hóa thủ công đơn thuần, thiếu vắng mô hình tích hợp tri thức bối cảnh sinh thái - nông lịch vào cấu trúc phân loại tự động và chưa định lượng hóa được cơ chế động lực chuyển đổi không gian theo chuỗi thời gian.

Luận án tiến sĩ địa lý chuyên ngành Bản đồ, Viễn thám và Hệ thông tin địa lý (Mã số: 62 44 02 14) của nghiên cứu sinh Lê Quang Toan, dưới sự hướng dẫn khoa học của TS. Bùi Quang Thành tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (2019), đã giải quyết triệt để bài toán này qua hai câu hỏi nghiên cứu (Research Questions - RQ) và hai giả thuyết tương ứng:

• RQ1: Làm thế nào để giải quyết tình trạng chồng lấn phổ giữa CCNLN với các trạng thái lớp phủ khác trên ảnh vệ tinh độ phân giải cao nhằm đạt độ tin cậy phân loại tối ưu?
  • Giả thuyết H1: Tích hợp mối quan hệ giữa điều kiện tự nhiên (địa hình, thổ nhưỡng), nông lịch mùa vụ và đặc điểm sinh thái cây trồng với thông tin phổ trong quy trình phân loại định hướng đối tượng (Object-Based Image Analysis - OBIA) sẽ cải thiện đáng kể độ chính xác tổng thể (Overall Accuracy), độ chính xác của người sản xuất (Producer's Accuracy - PA) và người sử dụng (User's Accuracy - UA).
• RQ2: Quy luật không gian - thời gian của quá trình chuyển đổi lớp phủ giữa rừng và CCNLN tại huyện Bảo Lâm diễn ra theo cơ chế nào và chịu sự chi phối định lượng bởi những yếu tố tự nhiên, kinh tế - xã hội nào?
  • Giả thuyết H2: Biến động diện tích CCNLN diễn ra theo chuỗi chuyển đổi có quy luật: Rừng $\rightarrow$ Đất cây hàng năm/Đất trống $\rightarrow$ Cây công nghiệp lâu năm; trong đó xác suất xảy ra chuyển đổi phụ thuộc phi tuyến tính vào các biến không gian (khoảng cách tới đường giao thông, sông suối, trung tâm hành chính, độ dốc, loại đất) thông qua mô hình hồi quy Logistic đa biến.

Khung lý thuyết của công trình được xây dựng vững chắc trên nền tảng Lý thuyết Mô hình hóa Tri thức Không gian (Knowledge-based Spatial Modeling Theory của Burrough, 1998 [54]), Lý thuyết Phân tích Ảnh Định hướng Đối tượng (OBIA Theory của Baatz & Schäpe, 2000 [46]) và Lý thuyết Biến động Sử dụng Đất / Lớp phủ (Land-Use/Cover Change - LUCC Theory của Lambin et al., 2001, 2003 [28, 84]). Phạm vi nghiên cứu bao quát toàn bộ 146.000 ha huyện Bảo Lâm qua chuỗi dữ liệu đa thời gian 12 năm (2004, 2011, 2016) cùng mốc đối chiếu lịch sử 1995, được kiểm chứng bởi 124 điểm khảo sát thực địa GPS chuẩn hóa, tạo nên một công trình có ý nghĩa đột phá cả về học thuật và thực tiễn quản trị tài nguyên lãnh thổ.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước phản ánh quá trình tiến hóa của các phương pháp viễn thám trong giám sát thảm thực vật từ phân tích phổ đơn giản đến các hệ chuyên gia phức hợp. Dòng nghiên cứu thứ nhất tập trung vào việc ứng dụng dữ liệu đa phổ quang học độ phân giải trung bình và cao (Landsat TM/ETM+/8, SPOT-5, Sentinel-2A) để lập bản đồ sinh khối và diện tích cây trồng nhiệt đới. Trong khi các nghiên cứu đối với cây cao su của Suratman và nnk (2004) [117] tại Malaysia hay Zhe Li và nnk (2011) [130] tại Thái Lan đạt độ chính xác khá cao (85% - 87%), việc lập bản đồ cà phê và chè luôn đối mặt với rào cản kỹ thuật nghiêm trọng. Tại Colombia, Cordero và nnk (2007) [60] khi dùng ảnh Landsat TM đơn lẻ chỉ đạt độ chính xác 37,4% cho cây cà phê do độ tán che thấp khiến tín hiệu phổ bị chi phối mạnh bởi nền đất và cây hàng năm. Tại Rwanda, Mukashema và nnk (2014) [102] kết hợp ảnh vệ tinh độ phân giải cao với mạng lưới xác suất Bayesian và bản đồ năng suất nhưng độ chính xác tổng thể cũng chỉ dừng lại ở mức 50%, qua đó đưa ra kết luận mang tính bi quan rằng không thể bóc tách diện tích cà phê chính xác nếu chỉ dựa vào giá trị phản xạ phổ đa kênh.

Dòng nghiên cứu thứ hai liên quan đến động thái biến động lớp phủ và suy thoái rừng tại Đông Nam Á. Nghiên cứu bước ngoặt của Meyfroidt và nnk (2013) [99] trên phạm vi Tây Nguyên đã chỉ rõ vai trò của sự bùng nổ giá nông sản toàn cầu sau giai đoạn mở cửa kinh tế và ảnh hưởng từ việc sụt giảm sản lượng cà phê của Brazil năm 1994 đã kích hoạt làn sóng di dân tự do và mở rộng đất canh tác nông nghiệp lấn sâu vào rừng tự nhiên giai đoạn 2000-2010. Tuy nhiên, nghiên cứu của Meyfroidt chủ yếu sử dụng tư liệu ảnh Landsat độ phân giải 30m ở quy mô cấp vùng, chưa phản ánh được các phân mảnh cảnh quan chi tiết ở cấp độ nông hộ và thửa đất tại quy mô cấp huyện.

Trong y văn tồn tại hai luồng quan điểm đối lập: Một trường phái cho rằng việc tăng cường độ phân giải không gian và độ phân giải phổ của đầu thu vệ tinh (vệ tinh siêu phổ Hyperion, vệ tinh độ phân giải mét như IKONOS, QuickBird) là đủ để giải quyết bài toán phân loại thực phủ. Ngược lại, trường phái thứ hai, được khẳng định bởi Wellens (1997) [123], nhấn mạnh: "Dù sao đi nữa, không có một phương thức 'chuẩn' nào có thể áp dụng được cho mọi vùng trên thế giới, nghĩa là từng phương pháp chỉ áp dụng tốt cho từng hoàn cảnh cụ thể và chỉ có ý nghĩa tham khảo khi áp dụng vào các khu vực khác". Quan điểm này được bổ trợ bởi phát hiện của Biging và nnk (1995) [51] về chuỗi sai số tích lũy từ thu nhận ảnh, tiền xử lý đến thuật toán phân loại.

Công trình của Lê Quang Toan định vị chính xác vào khoảng trống tri thức này: Khắc phục giới hạn của kỹ thuật phân loại điểm ảnh truyền thống bằng cách xây dựng quy trình phân loại định hướng đối tượng tích hợp bối cảnh (ĐHĐT-2). Luận án chứng minh rằng việc lồng ghép thông minh nông lịch (thời điểm thu hoạch cây hàng năm tính chỉ số $NDVI^*$), mô hình số độ cao DEM, bản đồ độ dốc, bản đồ đất và bản đồ quy hoạch 3 loại rừng vào cây quyết định phân cấp (Hierarchical Rule-based System) có thể bẻ gãy ranh giới nhầm lẫn phổ, vượt trội hơn hẳn các nghiên cứu quốc tế của Cordero [60] và Mukashema [102].

Đóng góp lý thuyết và khung phân tích

+-----------------------------------------------------------------------------------+
|               KHUNG PHÂN TÍCH TÍCH HỢP ĐỊNH HƯỚNG ĐỐI TƯỢNG VÀ HỒI QUY             |
+-----------------------------------------------------------------------------------+
|  1. DỮ LIỆU ĐA NGUỒN                                                              |
|     - Ảnh vệ tinh: SPOT-5 (2004, 2011), Sentinel-2A (2016), Landsat TM/8          |
|     - Dữ liệu bối cảnh: DEM, Độ dốc, Bản đồ đất, Nông lịch mùa vụ, Rừng FIPI      |
+-----------------------------------------------------------------------------------+
                                          |
                                          v
|  2. PHÂN PHẠM VI VÀ PHÂN ĐOẠN ĐA ĐỘ PHÂN GIẢI (OBIA Segmentation)                 |
|     - Tối ưu hóa tham số Scale, Shape/Color, Compactness/Smoothness               |
+-----------------------------------------------------------------------------------+
                                          |
                                          v
|  3. CÂY QUY TẮC PHÂN CẤP TÍCH HỢP TRI THỨC (ĐHĐT-2)                               |
|     - Tách Thực vật / Phi thực vật (NDVI*)                                        |
|     - Tách Rừng / Không rừng (Phổ cận hồng ngoại, kết cấu, FIPI)                  |
|     - Tách Cây hàng năm vs CCNLN (Nông lịch mùa vụ, NDVI sau thu hoạch)           |
|     - Bóc tách Cà phê vs Chè (Độ cao DEM, Độ dốc, Thổ nhưỡng Bazan)               |
+-----------------------------------------------------------------------------------+
                                          |
                                          v
|  4. ĐÁNH GIÁ MA TRẬN BIẾN ĐỘNG SAU PHÂN LOẠI (2004 - 2011 - 2016)                 |
|     - Chuỗi biến động: Rừng -> Đất cây hàng năm/Đất trống -> CCNLN                |
+-----------------------------------------------------------------------------------+
                                          |
                                          v
|  5. MÔ HÌNH HÓA ĐỘNG LỰC KHÔNG GIAN (Binary Logistic Regression)                  |
|     - Logit(P) = beta_0 + beta_1(Slope) + beta_2(Dist_Road) + beta_3(Dist_River) + ...|
+-----------------------------------------------------------------------------------+

Đóng góp cho lý thuyết

Luận án đóng góp trực tiếp vào việc mở rộng Lý thuyết Mô hình hóa Tri thức Chuyên gia Địa lý (Knowledge-based Geographic Expert System) của Burrough (1998) [54] và Khung phân tích LUCC của Lambin và nnk (2001, 2003) [28, 84] trong bối cảnh cảnh quan nông - lâm nghiệp nhiệt đới gió mùa:

1. Mệnh đề lý thuyết 1 (Proposition 1 - Sự tương tác bối cảnh - phổ): Độ chính xác định danh đối tượng lớp phủ thực vật nhân tạo không thuần túy là hàm số của độ phân giải phổ hay không gian, mà là hàm tích hợp giữa đặc tính phản xạ phổ bề mặt và các biến ranh giới sinh thái (Ecological Boundary Conditions) gồm độ cao tuyệt đối, độ dốc sườn đồi, nhóm đất phát triển và chu kỳ nông lịch sinh trưởng.
2. Mệnh đề lý thuyết 2 (Proposition 2 - Chuỗi chuyển đổi lớp phủ trung gian): Quá trình chuyển đổi từ rừng tự nhiên sang cây công nghiệp lâu năm không diễn ra trực tiếp tức thời trong một chu kỳ ngắn mà vận động theo cơ chế thế chân tuần tự thông qua giai đoạn trung gian: Đất rừng bị khai phá $\rightarrow$ Trồng cây hàng năm (ngô, sắn, lúa rẫy) hoặc để đất trống trong 1-3 năm đầu $\rightarrow$ Đầu tư vốn chuyển dịch sang cây công nghiệp lâu năm (cà phê, chè). Đây là phát hiện quan trọng bổ sung cho lý thuyết chuyển đổi rừng (Forest Transition Theory).
3. Mệnh đề lý thuyết 3 (Proposition 3 - Tính quy luật không gian phi tuyến): Sự mở rộng không gian của cây công nghiệp lâu năm tuân theo các quy luật kháng trở không gian (Spatial Resistance Laws), trong đó gradient khoảng cách tới mạng lưới hạ tầng giao thông và nguồn nước mặt đóng vai trò là biến số quyết định gia tốc chuyển đổi.

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 chính:

• Lý thuyết Phân tích Ảnh Định hướng Đối tượng (OBIA): Thay thế đơn vị phân loại pixel riêng lẻ bằng "đối tượng ảnh" (Image Object) đồng nhất về phổ, hình dạng (shape), độ nén chặt (compactness) và độ trơn láng (smoothness) thông qua thuật toán phân đoạn đa độ phân giải (Multi-resolution segmentation) của Baatz & Schäpe (2000) [46].
• Hệ thống Phân cấp Quy tắc Chuyên gia (Hierarchical Expert Rule-based System): Cấu trúc cây quyết định logic gồm hai cấp độ phân đoạn và phân loại lồng ghép. Cấp 1 tách nhóm Rừng và Không rừng; Cấp 2 bóc tách chi tiết các lớp đối tượng chuyên đề bằng cách kết hợp ngưỡng phổ vệ tinh với các dữ liệu chuyên đề ngoại sinh (ancillary layers).
• Lý thuyết Hồi quy Không gian Rời rạc (Discrete Spatial Regression): Vận dụng mô hình hồi quy Logistic nhị phân (Binary Logistic Regression) của Cox (1970) [61, 121] để chuyển đổi xác suất biến động nhị phân (có biến động = 1, không biến động = 0) sang hàm Logit liên tục, lượng hóa trọng số tác động của từng nhân tố kinh tế - xã hội và điều kiện tự nhiên.

Điều kiện biên (Boundary Conditions) của mô hình được xác lập rõ: Áp dụng tối ưu cho các hệ sinh thái cao nguyên nhiệt đới có tính phân bậc địa hình mạnh, canh tác nông nghiệp manh mún xen kẽ rừng tự nhiên và có sự phân hóa mùa khô - mùa mưa rõ rệt.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu được định vị theo hệ hình Thực chứng luận không gian (Spatial Positivism) kết hợp Hiện thực luận phản biện (Critical Realism). Nghiên cứu sử dụng thiết kế hỗn hợp đa cấp độ (Multi-level Mixed Design) kết hợp định lượng viễn thám số, phân tích không gian GIS và mô hình kinh tế lượng không gian.

                               CẤU TRÚC THIẾT KẾ ĐA CẤP ĐỘ
+------------------------------------------------------------------------------------+
| Cấp vĩ mô (Macro-level):    | Toàn huyện Bảo Lâm (146.000 ha) - Phân tích cấu trúc |
|                             | cảnh quan và quy hoạch sử dụng đất cấp tỉnh/huyện    |
+-----------------------------+------------------------------------------------------+
| Cấp trung mô (Meso-level):  | 14 đơn vị hành chính xã/thị trấn - Phân tích hồi quy  |
|                             | logistic và thống kê kinh tế - xã hội                |
+-----------------------------+------------------------------------------------------+
| Cấp vi mô (Micro-level):    | Thửa đất/Đối tượng ảnh (Objects) - 124 điểm khảo sát |
|                             | GPS thực địa và điểm ảnh độ phân giải 2.5m - 10m     |
+------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

1. Chiến lược chọn mẫu và dữ liệu đa cảm biến:
   • Ảnh vệ tinh quang học độ phân giải cao: SPOT-5 năm 2004 (kênh đa phổ 10m kết hợp toàn sắc 2.5m), SPOT-5 năm 2011 (kênh đa phổ 10m và toàn sắc 2.5m), Sentinel-2A năm 2016 (4 kênh 10m, 6 kênh 20m).
   • Ảnh vệ tinh đối chiếu mùa vụ và lịch sử: Landsat-5 TM năm 2004, Landsat-5 TM năm 2010, Landsat-8 OLI năm 2016 (độ phân giải 30m, dùng để tính chỉ số thực vật sai phân chuẩn hóa $NDVI$ tại thời điểm sau thu hoạch cây hàng năm - ký hiệu $NDVI^*$), SPOT-4 năm 1995 dùng làm nền đối soát rừng.
   • Bản đồ chuyên đề và dữ liệu phụ trợ: Bản đồ địa hình 1:25.000 (nội suy mô hình số độ cao DEM và bản đồ độ dốc bằng công cụ Spatial Analyst Slope), Bản đồ thổ nhưỡng tỉnh Lâm Đồng tỷ lệ 1:100.000 (MONRE, 2010), Bản đồ hiện trạng 3 loại rừng tỷ lệ 1:100.000 (FIPI, 2010), Bản đồ hiện trạng sử dụng đất năm 2010 huyện Bảo Lâm 1:25.000.
2. Thu thập dữ liệu thực địa và kiểm định tam giác hóa (Triangulation):
   • Thực hiện 2 đợt khảo sát thực địa chuyên sâu vào cuối năm 2011 và đầu năm 2013, thu thập thông tin tọa độ chuẩn xác của 124 điểm mô tả thực địa đại diện cho toàn bộ các kiểu thảm phủ (cà phê kinh doanh, cà phê kiến thiết cơ bản, chè, rừng tự nhiên lá rộng thường xanh, rừng trồng, đất cây hàng năm, đất trống).
   • Tam giác hóa phương pháp và dữ liệu: Kết hợp giữa định vị GPS sai số dưới 5m ngoài thực địa, giải đoán mắt chuyên gia trên bình đồ ảnh vệ tinh độ phân giải 2.5m, và đối chiếu hồ sơ thống kê kiểm kê đất đai hàng năm của Phòng Tài nguyên và Môi trường huyện Bảo Lâm.
3. Quy trình phân đoạn và phân loại ảnh OBIA:
   • Bước 1: Phân đoạn ảnh đa độ phân giải (Multi-resolution segmentation) trên phần mềm chuyên dụng (eCognition/ENVI). Thiết lập các thông số tối ưu: Scale parameter (Thông số tỷ lệ), Shape (Hình dạng = 0.2 - 0.3), Compactness (Độ nén chặt = 0.5), Color (Màu sắc = 0.7 - 0.8).
   • Bước 2: Xây dựng bộ quy tắc phân loại thuần phổ ĐHĐT-1 (chỉ dùng giá trị phổ các kênh SPOT/Sentinel) và bộ quy tắc tích hợp bối cảnh ĐHĐT-2. Trong ĐHĐT-2, áp dụng cây quyết định: $$\text{Nếu } NDVI^* < 0.3 \rightarrow \text{Đất trống / Cây hàng năm đã thu hoạch}$$ $$\text{Nếu } NDVI^* \ge 0.3 \text{ và thuộc vùng đất Bazan, độ cao } 800 - 1000\text{m, độ dốc } < 15^\circ \rightarrow \text{Cây cà phê}$$ $$\text{Nếu } NDVI^* \ge 0.3 \text{ và độ dốc } 15 - 25^\circ \text{ hoặc vùng trũng ẩm } \rightarrow \text{Cây chè}$$ $$\text{Nếu đối tượng nằm trong ranh giới rừng đặc dụng/phòng hộ FIPI } \rightarrow \text{Loại trừ nhầm lẫn CCNLN}$$

Data và phân tích

1. Đánh giá độ chính xác phân loại:
   • Xây dựng ma trận nhầm lẫn (Confusion Matrix / Error Matrix), tính toán Độ chính xác tổng thể (Overall Accuracy - OA), Hệ số Kappa ($\hat{K}$), Độ chính xác của người sản xuất (Producer's Accuracy - PA) và Độ chính xác của người sử dụng (User's Accuracy - UA).
2. Phân tích ma trận biến động sau phân loại (Post-Classification Comparison):
   • Chồng xếp lớp bản đồ phân loại 2004, 2011 và 2016 dạng raster trong môi trường ArcGIS để tính toán ma trận chuyển dịch diện tích giữa các cặp trạng thái: Mất rừng, Suy thoái rừng, Mở rộng cây hàng năm, Mở rộng CCNLN.
3. Mô hình hồi quy Logistic đa biến xác định động lực biến động:
   • Thiết lập phương trình Logit nhị phân theo Cox (1970) [61]: $$\ln\left(\frac{P}{1 - P}\right) = \beta_0 + \beta_1 X_1 + \beta_2 X_2 + \dots + \beta_k X_k$$ Trong đó: $P$ là xác suất xảy ra biến động lớp phủ ($P \in [0, 1]$); $\frac{P}{1-P}$ là tỷ số chênh (Odds Ratio); $\beta_0$ là hằng số hồi quy; $\beta_i$ là các hệ số hồi quy riêng phần; $X_i$ là các biến độc lập không gian gồm: Độ cao ($X_1$), Độ dốc ($X_2$), Khoảng cách tới đường giao thông chính ($X_3$), Khoảng cách tới sông suối ($X_4$), Khoảng cách tới ranh giới rừng ($X_5$), Khoảng cách tới trung tâm xã/UBND ($X_6$), và Mật độ dân số cấp xã ($X_7$).
   • Phân tích thống kê thực hiện trên phần mềm SPSS, đánh giá mức ý nghĩa thống kê qua kiểm định Wald, $p$-value ($p < 0.05$ và $p < 0.01$), độ phù hợp mô hình qua chỉ số $-2\text{ Log Likelihood}$ và hệ số $R^2$ hiệu chỉnh (Cox & Snell $R^2$, Nagelkerke $R^2$).

Phát hiện đột phá và implications

Những phát hiện then chốt

1. Đột phá về kỹ thuật bóc tách phổ thực vật phức tạp: Phương pháp ĐHĐT-2 tích hợp chỉ số nông lịch $NDVI^*$ và dữ liệu địa hình DEM đã giải quyết triệt để sự nhầm lẫn giữa cà phê kiến thiết cơ bản (chưa khép tán) và đất nông nghiệp trồng cây hàng năm. Độ chính xác tổng thể của bản đồ hiện trạng lớp phủ năm 2011 tăng từ 76,5% lên 89,2% (Hệ số Kappa tăng từ 0.71 lên 0.86). Diện tích cà phê và chè được hiệu chỉnh phân loại chính xác hàng ngàn hecta so với các thuật toán truyền thống.
2. Khám phá quy luật động thái chuyển dịch lớp phủ giai đoạn 2004-2016: Kết quả phân tích chuỗi ma trận biến động chỉ ra rằng diện tích CCNLN (chủ yếu là cà phê) tại Bảo Lâm liên tục tăng trưởng mạnh mẽ, từ 58.400 ha (năm 2004) lên 66.800 ha (năm 2011) và đạt trên 71.200 ha (năm 2016). Diện tích này mở rộng chủ yếu thông qua việc chuyển đổi từ đất cây hàng năm và lấn chiếm đất rừng tự nhiên.
3. Lượng hóa chuỗi mất rừng trung gian: Phân tích nguồn gốc đất năm 2004 đối với diện tích CCNLN mới phát sinh năm 2011 và 2016, kết hợp chồng xếp nền rừng lịch sử năm 1995, đã chứng minh: 83,4% diện tích CCNLN mở rộng không trực tiếp xâm lấn vào rừng nguyên sinh trong cùng năm, mà đất rừng ban đầu bị chặt phá để trồng hoa màu ngắn ngày (ngô, sắn), sau 2 đến 4 năm canh tác khi đất đã ổn định và người dân tích lũy đủ vốn, họ mới trồng thay thế bằng cà phê và chè.
4. Phát hiện phản trực giác về động lực giao thông và nguồn nước (Logistic Regression): Mô hình hồi quy Logistic ($p < 0.001$, Nagelkerke $R^2 = 0.432$) chỉ ra mối tương quan nghịch biến cực kỳ chặt chẽ giữa khoảng cách tới đường giao thông và nguy cơ mất rừng:
   • Các khu vực nằm trong bán kính cách trục đường giao thông chính dưới 500m có tỷ số chênh chuyển đổi sang CCNLN cao gấp 4,82 lần ($\text{Odds Ratio} = 4.82, \text{CI } 95%: [3.41 - 6.81]$) so với các khu vực cách xa trên 2.000m.
   • Ngược lại với giả định truyền thống cho rằng địa hình hiểm trở bảo vệ được rừng, tại Bảo Lâm, các sườn đồi có độ dốc từ $8^\circ - 18^\circ$ lại có nguy cơ bị chuyển đổi sang trồng chè và cà phê cao hơn vùng bằng phẳng thung lũng (do vùng thung lũng đã bão hòa và bị ngập úng mùa mưa không thích hợp cho cà phê).

          MÔ THỨC CHUYỂN DỊCH KHÔNG GIAN THEO GRADIENT GIAO THÔNG VÀ ĐỘ DỐC
[Trục giao thông chính] 
       |--- (0 - 500m)   ---> Xác suất chuyển đổi sang Cà phê cực đại (Odds Ratio = 4.82)
       |--- (500 - 1500m) ---> Khu vực mở rộng Đất cây hàng năm (Giai đoạn trung gian)
       |--- (> 2000m)    ---> Rừng tự nhiên còn sót lại / Rừng phòng hộ hiểm trở
[Gradient Độ dốc]
       |--- (0 - 8 độ)   ---> Vùng thung lũng bằng: Đô thị, thổ cư, lúa nước, chè trũng
       |--- (8 - 18 độ)  ---> Sườn đồi thoải Bazan: Điểm nóng mở rộng Cà phê chuyên canh
       |--- (> 25 độ)    ---> Rừng tự nhiên trên núi cao (Tiou Hoan)

Implications đa chiều

• Đóng góp lý thuyết: Khẳng định tính tất yếu của việc tích hợp tri thức địa lý - sinh thái vào các thuật toán viễn thám số, chứng minh rằng xử lý ảnh vệ tinh hiện đại không thể tách rời quy luật địa lý bề mặt và đặc trưng nhân văn lãnh thổ.
• Đổi mới phương pháp luận: Cung cấp bộ quy tắc phân loại ĐHĐT chuẩn hóa có thể tái lập (reproducible) và mở rộng cho toàn bộ 5 tỉnh Tây Nguyên cũng như các vùng sinh thái chuyên canh cây nhiệt đới trên thế giới (vùng cao nguyên Đông Phi, Nam Mỹ, Đông Nam Á).
• Ứng dụng thực tiễn: Cung cấp bộ cơ sở dữ liệu số GIS và tập bản đồ hiện trạng - biến động lớp phủ tỷ lệ 1:25.000 có độ chính xác cao tuyệt đối, phục vụ trực tiếp cho UBND huyện Bảo Lâm và Sở Nông nghiệp & PTNT tỉnh Lâm Đồng trong việc giám sát cơ cấu cây trồng.
• Khuyến nghị chính sách: Cung cấp cơ sở khoa học để chính quyền địa phương ngăn chặn tình trạng "hợp thức hóa đất rừng" thông qua chuỗi chuyển đổi trung gian; thiết lập hành lang giám sát viễn thám tự động đa thời gian cho các dự án chi trả dịch vụ môi trường rừng (PES) và chương trình giảm phát thải khí nhà kính từ mất rừng và suy thoái rừng (REDD+).

Limitations và Future Research

1. Hạn chế về dữ liệu ảnh quang học do mây mù: Do đặc thù khí hậu cao nguyên Bảo Lộc - Di Linh có lượng mưa lớn và độ ẩm cao quanh năm, việc tìm kiếm các cảnh ảnh quang học SPOT-5 và Sentinel-2 hoàn toàn không có mây ở cùng thời điểm mùa khô giữa các năm là thách thức lớn, đòi hỏi phải sử dụng kỹ thuật ghép ảnh đa thời điểm lân cận.
2. Giới hạn trong phân định tuổi cây và cấu trúc tầng tán: Nghiên cứu mới chỉ phân loại được CCNLN ở cấp độ loài chính (cà phê, chè) mà chưa phân định chi tiết được các giai đoạn tuổi sinh trưởng (cây 1-3 năm tuổi, cây kinh doanh thu hoạch, cây già cỗi cần tái canh) do thiếu dữ liệu viễn thám siêu phổ (Hyperspectral) và dữ liệu quét laser LiDAR.
3. Chưa tích hợp dữ liệu kinh tế lượng vi mô hộ gia đình: Mô hình hồi quy Logistic mới chỉ sử dụng các biến số không gian địa lý và mật độ dân số cấp xã, chưa lồng ghép được các biến kinh tế lượng vi mô như biến động giá cà phê thị trường thế giới hàng tháng, thu nhập bình quân hộ nông dân và chính sách tín dụng vi mô.

Chương trình nghiên cứu tương lai (Future Research Agenda) 5-10 năm tới bao gồm:

• Tích hợp chuỗi ảnh radar khẩu độ tổng hợp (SAR) Sentinel-1 (băng C) và ALOS PALSAR-2 (băng L) xuyên mây kết hợp ảnh quang học Sentinel-2 để theo dõi sinh khối cây công nghiệp theo thời gian thực (Real-time monitoring).
• Ứng dụng các thuật toán học sâu (Deep Learning - CNN, U-Net) trên nền tảng điện toán đám mây Google Earth Engine để tự động hóa hoàn toàn quy trình giám sát rừng và CCNLN quy mô quốc gia.
• Mở rộng mô hình kinh tế lượng không gian tích hợp dữ liệu vi mô nông hộ (Agent-Based Modeling - ABM) để mô phỏng và dự báo các kịch bản chuyển đổi sử dụng đất dưới tác động của biến đổi khí hậu toàn cầu.

Tác động và ảnh hưởng

• Tác động học thuật: Luận án mở ra hướng tiếp cận chuẩn mực trong việc kết hợp viễn thám OBIA và mô hình hóa không gian tại Việt Nam; tạo tiền đề cho hàng loạt công bố khoa học trên các tạp chí chuyên ngành uy tín trong nước và quốc tế; dự kiến thu hút trích dẫn cao trong các nghiên cứu về địa lý tài nguyên, lâm nghiệp và biến động môi trường vùng Tây Nguyên.
• Tác động ngành nông nghiệp và lâm nghiệp: Chuyển đổi phương thức quản lý cây công nghiệp truyền thống dựa trên báo cáo hành chính cơ sở sang phương thức quản lý số hóa chính xác dựa trên bản đồ viễn thám và GIS; hỗ trợ tái cơ cấu ngành chè và cà phê Lâm Đồng theo hướng bền vững.
• Tác động chính sách công: Cung cấp luận cứ kỹ thuật chuẩn xác phục vụ việc điều chỉnh quy hoạch sử dụng đất giai đoạn 2020-2030 của tỉnh Lâm Đồng; hỗ trợ thiết lập hệ thống cảnh báo sớm mất rừng phục vụ thực thi Luật Lâm nghiệp và các cam kết quốc tế về bảo tồn đa dạng sinh học.
• Lợi ích xã hội và môi trường: Bảo vệ tính toàn vẹn của lưu vực sông La Ngà và hệ thống thủy điện, hạn chế xói mòn rửa trôi đất bazan trên các sườn dốc, bảo vệ sinh kế lâu dài cho đồng bào dân tộc thiểu số bản địa trước áp lực chuyển đổi đất nông nghiệp.

Đối tượng hưởng lợi

• Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung phương pháp luận hoàn chỉnh, giải quyết bài toán chồng lấn phổ thực vật và quy trình kiểm định mô hình hồi quy không gian chặt chẽ.
• Cơ quan quản lý Nhà nước (Sở Tài nguyên & Môi trường, Sở Nông nghiệp & PTNT Lâm Đồng, UBND huyện Bảo Lâm): Sở hữu công cụ giám sát biến động diện tích đất đai khách quan, minh bạch, giảm thiểu 70% chi phí và thời gian điều tra thực địa thủ công.
• Ban Quản lý Rừng và Doanh nghiệp Nông lâm nghiệp: Xác định chính xác các tọa độ ranh giới xung đột đất đai, tối ưu hóa vùng nguyên liệu chè và cà phê chất lượng cao theo độ cao và thổ nhưỡng phù hợp.
• Các tổ chức quốc tế và dự án phi chính phủ (SNV, UNDP, USAID, IUCN): Khai thác bộ cơ sở dữ liệu tin cậy để đánh giá hiệu quả của các chương trình REDD+, chi trả dịch vụ môi trường rừng và các sáng kiến cảnh quan nông nghiệp bền vữ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?

Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Mô hình hóa Tri thức Không gian (Knowledge-based Spatial Modeling Theory) của Burrough (1998) [54] thông qua việc thiết lập nguyên lý tích hợp "Bộ ba bối cảnh": Nông lịch mùa vụ ($NDVI^$ sau thu hoạch) + Địa hình thích nghi sinh thái (DEM/Độ dốc) + Thổ nhưỡng chuyên biệt (Bazan)* vào cấu trúc phân đoạn ảnh định hướng đối tượng. Luận án đã bác bỏ định kiến kỹ thuật cho rằng không thể phân loại cây cà phê chính xác từ viễn thám quang học, chứng minh rằng tri thức địa lý có thể bù đắp hoàn toàn sự thiếu hụt về tín hiệu phổ.

2. Đổi mới phương pháp luận có gì vượt trội so với các nghiên cứu quốc tế trước đây?

So sánh với nghiên cứu của Cordero và nnk (2007) [60] tại Colombia (chỉ đạt OA 37,4% trên ảnh Landsat) và Mukashema và nnk (2014) [102] tại Rwanda (dùng ảnh phân giải cao kết hợp mạng Bayesian đạt OA ~50%), phương pháp ĐHĐT-2 của luận án đã nâng độ chính xác tổng thể lên 89,2% và độ chính xác phân loại cà phê/chè đạt trên 88%. Sự vượt trội này đạt được nhờ việc sử dụng đối tượng phân đoạn đa độ phân giải kết hợp cây quyết định phân cấp tích hợp dữ liệu chuyên đề FIPI và chỉ số $NDVI^*$ theo mùa.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là cơ chế chuyển đổi đất trung gian kéo dài 2-4 năm: Người dân không chuyển đổi trực tiếp rừng thành vườn cà phê ngay lập tức mà luôn canh tác cây hàng năm ngắn ngày (ngô, sắn) hoặc để đất trống như một bước "đệm" nhằm né tránh sự phát hiện của cơ quan quản lý và tích lũy vốn đầu tư. Đồng thời, các sườn đồi dốc trung bình ($8^\circ - 18^\circ$) có tốc độ mất rừng và chuyển sang CCNLN nhanh hơn các vùng đất bằng thung lũng do đặc tính sinh thái ưa thoát nước của cây cà phê và chè.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án mô tả chi tiết: (1) Bộ tham số phân đoạn ảnh eCognition (Scale: 30-50, Shape: 0.2, Compactness: 0.5); (2) Cây logic bộ quy tắc ĐHĐT-1 và ĐHĐT-2 với ngưỡng giá trị cụ thể cho từng kênh phổ SPOT-5 và Sentinel-2A; (3) Quy trình chuẩn hóa hệ tọa độ VN-2000 / UTM Zone 48N; (4) Ma trận 124 điểm khảo sát GPS thực địa và phương trình hồi quy Logistic chuẩn hóa trong SPSS. Toàn bộ quy trình có thể chuyển giao và áp dụng ngay cho các địa bàn tương đồng.

5. Kế hoạch nghiên cứu 10 năm tới được phác thảo như thế nào?

Kế hoạch nghiên cứu 10 năm tập trung vào 3 hướng chiến lược: (1) Xây dựng hệ thống viễn thám SAR đa thời gian (Sentinel-1/ALOS-2) kết hợp AI để giám sát tự động biến động rừng và CCNLN theo tuần xuyên mây mù; (2) Ứng dụng công nghệ Big Data viễn thám trên Google Earth Engine để mở rộng phạm vi nghiên cứu ra toàn bộ lưu vực sông Đồng Nai và vùng Tây Nguyên; (3) Tích hợp mô hình biến đổi khí hậu RCP 4.5/8.5 để mô phỏng sự dịch chuyển vành đai sinh thái cây cà phê và chè lên các cao độ lớn hơn trong tương lai.

Kết luận

Công trình luận án tiến sĩ của Lê Quang Toan (2019) đại diện cho một bước tiến khoa học xuất sắc, kết hợp nhuần nhuyễn giữa kỹ thuật viễn thám hiện đại, phân tích không gian GIS và các quy luật kinh tế - xã hội học địa lý. Các kết luận cốt lõi bao gồm:

1. Xác lập thành công quy trình phân loại ĐHĐT-2 tích hợp tri thức bối cảnh: Nâng độ chính xác tổng thể phân loại hiện trạng lớp phủ lên 89,2% (Kappa = 0.86), bóc tách chính xác tuyệt đối diện tích cây cà phê và chè huyện Bảo Lâm trên ảnh SPOT-5 và Sentinel-2A.
2. Định lượng hóa toàn diện quy mô và tốc độ chuyển dịch lớp phủ giai đoạn 2004-2016: Ghi nhận sự mở rộng liên tục của diện tích CCNLN từ 58.400 ha lên trên 71.200 ha, đi kèm với sự suy giảm nghiêm trọng của diện tích rừng tự nhiên.
3. Làm sáng tỏ quy luật chuyển đổi rừng theo chuỗi thế chân trung gian: Chứng minh 83,4% diện tích CCNLN mới mở rộng có nguồn gốc sâu xa từ đất rừng bị khai phá thông qua giai đoạn đệm trồng cây hàng năm.
4. Xây dựng mô hình hồi quy Logistic giải thích động lực không gian: Xác định khoảng cách tới trục đường giao thông ($< 500$m, $\text{Odds Ratio} = 4.82$) và sườn dốc $8^\circ - 18^\circ$ là những yếu tố kích hoạt chuyển dịch mạnh nhất.
5. Đề xuất hệ thống giải pháp quản lý lãnh thổ khả thi: Cung cấp công cụ khoa học chuẩn xác phục vụ quy hoạch nông nghiệp bền vững, bảo vệ tài nguyên rừng và thực thi các cam kết giảm phát thải REDD+ tại tỉnh Lâm Đồng và vùng Tây Nguyên.