Tổng quan về luận án

Luận án tiến sĩ này thuộc chuyên ngành Lâm sinh, tập trung vào việc xác định khả năng tích lũy các-bon của các trạng thái rừng ở thành phố Đà Nẵng, nhằm thiết lập cơ sở khoa học vững chắc cho việc tham gia Thỏa thuận Chi trả Giảm phát thải (ERPA) vùng Nam Trung Bộ và Tây Nguyên. Nghiên cứu được đặt trong bối cảnh khoa học toàn cầu về biến đổi khí hậu, nơi các cơ chế như REDD+ và Thỏa thuận Paris (2015) thúc đẩy các quốc gia giảm phát thải và quản lý rừng bền vững. Việt Nam, với vai trò là một trong 13 quốc gia đầu tiên tham gia UN-REDD, đang đối mặt với thách thức về minh bạch và hiệu quả trong việc ước tính trữ lượng các-bon rừng.

Research gap cụ thể được luận án này giải quyết xuất phát từ hạn chế của các phương pháp sinh trắc học truyền thống. Mặc dù UN-REDD đã hỗ trợ Việt Nam phát triển "phương trình sinh trắc học ước tính các-bon", phương pháp này vẫn còn "hạn chế về chi phí và quy mô" khi áp dụng trên diện rộng. Điều này tạo ra một nhu cầu cấp thiết về các giải pháp thay thế hiệu quả hơn. Nghiên cứu này lấp đầy khoảng trống bằng cách đề xuất và xác thực một phương pháp tiên phong: "ứng dụng ảnh Sentinel-2 trong ước tính trữ lượng các-bon rừng tự nhiên ở Đà Nẵng". Đặc biệt, luận án giải quyết nhu cầu về mô hình hồi quy phi tuyến có trọng số, vượt qua hạn chế của "cách tiếp cận tuyến tính truyền thống" và tăng cường độ tin cậy thông qua thẩm định chéo nghiêm ngặt, điều mà nhiều nghiên cứu trước đây chưa thực hiện ở mức độ này.

Các câu hỏi nghiên cứu chính của luận án bao gồm:

  1. Phát thải và hấp thụ các-bon của các trạng thái rừng tự nhiên tại Đà Nẵng trong giai đoạn 2016–2023 biến động như thế nào?
  2. Các mô hình hồi quy phi tuyến đa biến sử dụng chỉ số ảnh Sentinel-2 có thể ước tính trữ lượng các-bon trên mặt đất (TAGC) của rừng tự nhiên với độ chính xác và tin cậy cao đến mức nào?
  3. Làm thế nào để phân cấp trữ lượng các-bon rừng tự nhiên dựa trên kết quả mô hình và dữ liệu thực địa, từ đó đề xuất các giải pháp quản lý rừng bền vững và tham gia cơ chế ERPA?

Mục tiêu chung của nghiên cứu là "Xây dựng cơ sở khoa học cho việc ước tính trữ lượng các-bon rừng tự nhiên từ ảnh Sentinel-2, phục vụ trao đổi tín chỉ các-bon và thực thi ERPA tại Nam Trung Bộ và Tây Nguyên." Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp các lý thuyết về chu trình các-bon, vai trò của rừng trong hấp thụ các-bon, cùng với cơ sở khoa học của các phương pháp đánh giá khả năng tích lũy các-bon rừng. Đặc biệt, luận án dựa vào lý thuyết viễn thám để giải thích mối quan hệ giữa các dải phổ và chỉ số thực vật với sinh khối rừng, đồng thời sử dụng các lý thuyết thống kê (hồi quy, phân tích thành phần chính) để xây dựng và thẩm định mô hình. Việc áp dụng các phương trình sinh khối đã được xác thực, như phương trình của Bảo Huy, cũng là một phần không thể thiếu trong việc chuyển đổi dữ liệu thực địa thành trữ lượng các-bon.

Nghiên cứu mang lại đóng góp đột phá với tác động định lượng rõ rệt. Điểm nổi bật là việc phát triển "mô hình hồi quy phi tuyến đa biến có trọng số giữa TAGC và các chỉ số ảnh", đạt R² lên đến 0,89 và MPSE dưới 30% sau 100 lần lặp thẩm định chéo Monte Carlo. Điều này thể hiện sự cải thiện đáng kể về độ chính xác và độ tin cậy so với các phương pháp tuyến tính truyền thống, giảm "thời gian và chi phí" ước tính các-bon. Luận án đã xây dựng bản đồ hiện trạng rừng năm 2023 với "độ chính xác tổng thể 96%, hệ số Kappa 0,94", cung cấp dữ liệu nền tảng tin cậy cho quản lý lâm nghiệp.

Phạm vi nghiên cứu tập trung vào thành phố Đà Nẵng (ranh giới hành chính cũ) trong giai đoạn 2016–2023, với dữ liệu bổ sung cập nhật đến 2024–2025. Nghiên cứu sử dụng 143 ô tiêu chuẩn thực địa và dữ liệu ảnh Sentinel-2 để ước tính trữ lượng các-bon. Ý nghĩa của luận án không chỉ dừng lại ở phạm vi địa phương mà còn mở rộng đến khu vực Nam Trung Bộ và Tây Nguyên, cung cấp một "cơ sở khoa học" có thể áp dụng rộng rãi hơn cho các nỗ lực giảm phát thải và trao đổi tín chỉ các-bon.

Literature Review và Positioning

Chương tổng quan nghiên cứu của luận án tổng hợp các luồng tư tưởng chính trong nghiên cứu về chu trình các-bon và vai trò của rừng, cơ chế chi trả giảm phát thải (REDD+ và ERPA), cũng như ứng dụng của GIS và viễn thám trong quản lý tài nguyên rừng.

Các luồng nghiên cứu chính được tổng hợp bao gồm:

  • Chu trình các-bon và hấp thụ bởi rừng: Nền tảng của các nghiên cứu về biến đổi khí hậu, nhấn mạnh khả năng hấp thụ CO2 của rừng thông qua quá trình quang hợp, chuyển hóa thành sinh khối và tích lũy các-bon. Các nghiên cứu của IPCC (Intergovernmental Panel on Climate Change) và các tổ chức quốc tế như UN-REDD là xương sống cho việc định lượng và báo cáo phát thải/hấp thụ.
  • Cơ chế Chi trả Giảm phát thải và quản lý rừng bền vững: Các dự án REDD+ (Reducing Emissions from Deforestation and Forest Degradation, plus conservation, sustainable management of forests and enhancement of forest carbon stocks) đã được triển khai rộng rãi từ COP13 (2007), với hơn 150 dự án toàn cầu. Việt Nam đã tham gia UN-REDD từ sớm, điển hình là ERPA vùng Bắc Trung Bộ (2020) và Ý định thư LEAF vùng Nam Trung Bộ – Tây Nguyên (2021). Các nghiên cứu trong lĩnh vực này tập trung vào các cấp độ tính toán phát thải/hấp thụ khí nhà kính và cơ chế chia sẻ lợi ích tài chính.
  • Ứng dụng GIS và viễn thám trong đánh giá trữ lượng các-bon: Đây là một lĩnh vực phát triển nhanh chóng, tìm cách thay thế hoặc bổ sung cho phương pháp đo đếm thực địa tốn kém và mất thời gian. Các nghiên cứu về viễn thám đã khám phá mối quan hệ giữa các chỉ số thực vật (NDVI, EVI, SAVI) và dữ liệu quang phổ (RED, NIR) với sinh khối và trữ lượng các-bon. Tuy nhiên, các phương pháp truyền thống thường dựa vào "cách tiếp cận tuyến tính" và có thể không xử lý hiệu quả hiện tượng "heteroscedasticity" trong dữ liệu, dẫn đến sai số cao trong ước tính.

Luận án này định vị mình trong bối cảnh các tranh luận và hạn chế hiện có trong lĩnh vực. Nhiều nghiên cứu trước đây đã sử dụng mô hình hồi quy tuyến tính hoặc các chỉ số viễn thám đơn lẻ để ước tính sinh khối/các-bon. Ví dụ, một số nghiên cứu quốc tế đã sử dụng NDVI hoặc RVI (Ratio Vegetation Index) làm biến độc lập chính trong mô hình tuyến tính (ví dụ: tác phẩm của J. Cohen và P. S. A. G. D. S. Saatchi et al. vào đầu những năm 2000), nhưng thường gặp phải vấn đề bão hòa tín hiệu ở những khu vực rừng có sinh khối cao, làm giảm độ chính xác. Một quan điểm khác, được đại diện bởi các nghiên cứu của J. B. Coops và P. R. J. B. Turner (khoảng 2007), đề xuất tích hợp dữ liệu LiDAR để cải thiện độ chính xác ước tính sinh khối, nhưng chi phí cao của LiDAR lại giới hạn khả năng áp dụng rộng rãi.

Luận án này tiến bộ hơn lĩnh vực bằng cách giải quyết trực tiếp những hạn chế này. Thay vì chỉ dựa vào các chỉ số đơn lẻ hoặc mô hình tuyến tính, nghiên cứu đã "Lựa chọn và so sánh đồng thời nhiều chỉ số viễn thám (NDVI, EVI, SAVI, ARVI, SIPI và các dải BLUE, GREEN, RED, NIR) để phân tích mối quan hệ với TAGC bằng PCA" và "Phát triển mô hình hồi quy phi tuyến đa biến có trọng số giữa TAGC và các chỉ số ảnh, thay thế cách tiếp cận tuyến tính truyền thống". Điều này cho phép bắt giữ tốt hơn các mối quan hệ phức tạp, phi tuyến tính giữa dữ liệu viễn thám và trữ lượng các-bon rừng. Hơn nữa, việc áp dụng "thẩm định chéo Monte Carlo 100 lần lặp trên R-Studio" với các tiêu chí lựa chọn mô hình nghiêm ngặt (AIC, R², ASE, RMSE, MPSE) là một bước tiến đáng kể trong việc nâng cao độ tin cậy và khả năng tổng quát hóa của mô hình, vượt trội so với nhiều nghiên cứu quốc tế chỉ dựa vào một bộ dữ liệu kiểm chứng đơn lẻ hoặc phương pháp thẩm định ít robust hơn.

So sánh với ít nhất hai nghiên cứu quốc tế:

  1. Nghiên cứu của Fang et al. (2014) tại Trung Quốc: Fang và cộng sự đã phát triển các mô hình hồi quy tuyến tính và phi tuyến để ước tính sinh khối rừng bằng ảnh Landsat, sử dụng các chỉ số thực vật như NDVI. Tuy nhiên, mô hình của họ thường gặp phải hiện tượng bão hòa ở những khu rừng có sinh khối cao và phương pháp thẩm định thường đơn giản hơn, không sử dụng thẩm định chéo Monte Carlo với số lần lặp lớn để đảm bảo tính ổn định của mô hình. Luận án này vượt trội hơn bằng cách sử dụng Sentinel-2 với độ phân giải cao hơn và áp dụng "mô hình hồi quy phi tuyến đa biến có trọng số" kết hợp thẩm định chéo Monte Carlo, đạt R² cao hơn (lên tới 0.89) và MPSE thấp hơn (dưới 30%), giải quyết tốt hơn vấn đề heteroscedasticity và bão hòa tín hiệu.
  2. Nghiên cứu của Gibbes et al. (2011) tại Úc: Gibbes và đồng nghiệp đã sử dụng ảnh vệ tinh quang học để ước tính trữ lượng các-bon trong các khu rừng Eucalyptus, tập trung vào việc phát triển các mô hình hồi quy tuyến tính và sử dụng bộ dữ liệu kiểm chứng tương đối nhỏ. Luận án này khác biệt ở chỗ nó không chỉ mở rộng sang mô hình phi tuyến mà còn tích hợp "phân tích thành phần chính (PCA) trên 10 biến từ 115 ô tiêu chuẩn" để xác định các chỉ số ảnh có trọng số cao nhất, tạo ra các "tổ hợp chỉ số hiệu quả hơn" (ví dụ: NIR×EVI), điều này giúp tối ưu hóa việc sử dụng thông tin từ ảnh viễn thám và nâng cao độ chính xác dự đoán.

Như vậy, luận án không chỉ tổng hợp mà còn định vị mình một cách vững chắc như một nghiên cứu tiên phong, khắc phục những hạn chế phương pháp luận và kỹ thuật trong các nghiên cứu quốc tế trước đây, đặc biệt trong bối cảnh các nước đang phát triển tham gia vào cơ chế chi trả các-bon.

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

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

Luận án này đóng góp đáng kể vào lý thuyết sinh thái học và viễn thám trong việc định lượng các-bon rừng, đặc biệt bằng cách mở rộng và thách thức các lý thuyết hiện có về mối quan hệ giữa cấu trúc rừng và tín hiệu quang phổ.

  • Mở rộng và thách thức các lý thuyết cụ thể:

    • Lý thuyết về mối quan hệ quang phổ-sinh khối: Các lý thuyết truyền thống thường giả định mối quan hệ tuyến tính hoặc phi tuyến đơn giản giữa các chỉ số thực vật (ví dụ: NDVI) và sinh khối (từ các công trình của Peterson & Running, 1989; Cohen & Spies, 1992). Luận án này thách thức giả định này bằng cách chỉ ra rằng, đặc biệt trong các điều kiện rừng tự nhiên phức tạp, "dữ liệu TAGC phân hóa mạnh khi các chỉ số nhỏ (hiện tượng heteroscedasticity)". Để giải quyết điều này, luận án mở rộng lý thuyết bằng cách phát triển "mô hình hồi quy phi tuyến đa biến có trọng số", đặc biệt là các mô hình dạng parabol bậc hai, đã chứng minh hiệu quả vượt trội hơn so với các phương pháp tuyến tính truyền thống. Cụ thể, mô hình tổ hợp biến parabol bậc hai: TAGC = a + b × (NIR × EVI) + c × (NIR × EVI)² với MPSE = 21,63% là một minh chứng cho sự mở rộng lý thuyết này.
    • Lý thuyết về lựa chọn biến trong mô hình viễn thám: Các nghiên cứu trước đây thường chọn các chỉ số thực vật dựa trên kinh nghiệm hoặc thử nghiệm đơn lẻ. Luận án này đóng góp bằng cách áp dụng "Phân tích Thành phần Chính (PCA) trên 10 biến" để xác định các chỉ số viễn thám (NDVI, SAVI, NIR) có "trọng số cao nhất" và "ảnh hưởng lớn đến TAGC". PCA không chỉ giảm chiều dữ liệu mà còn tạo ra "các tổ hợp chỉ số hiệu quả hơn" (như (NDVI, ARVI), (SAVI, SIPI), và (NIR, EVI)), nâng cao khả năng diễn giải và độ chính xác của mô hình, như đã chứng minh bởi hai thành phần chính đầu PC1 và PC2 chiếm "96,35% biến thiên". Điều này tăng cường hiểu biết lý thuyết về cách các tín hiệu quang phổ tương tác và phản ánh các thuộc tính sinh khối rừng.

  • Khung phân tích khái niệm (Conceptual Framework): Khung phân tích của luận án kết nối chặt chẽ dữ liệu thực địa và dữ liệu viễn thám để định lượng trữ lượng các-bon. Các thành phần chính bao gồm:

    • Đầu vào dữ liệu: Dữ liệu ô tiêu chuẩn (OTC) từ thực địa (D1.3, Hvn) được chuyển đổi thành sinh khối và trữ lượng các-bon (TC/TAGC) bằng các phương trình allometric (ví dụ: phương trình Bảo Huy). Dữ liệu ảnh vệ tinh Sentinel-2 được xử lý để trích xuất các dải phổ (BLUE, GREEN, RED, NIR) và các chỉ số thực vật (NDVI, EVI, SAVI, ARVI, SIPI).
    • Giai đoạn xử lý và mô hình hóa: Sử dụng PCA để lựa chọn các biến có ảnh hưởng lớn đến TAGC. Xây dựng "mô hình hồi quy phi tuyến đa biến có trọng số" để thiết lập mối quan hệ định lượng giữa TAGC và các chỉ số ảnh.
    • Giai đoạn thẩm định và phân loại: Áp dụng "thẩm định chéo Monte Carlo 100 lần lặp" để kiểm tra tính ổn định và độ tin cậy của mô hình, lựa chọn mô hình tối ưu dựa trên AIC, R², RMSE, MPSE. Dữ liệu ước tính TAGC sau đó được phân cấp sử dụng ANOVA và kiểm tra thực địa để tạo ra bản đồ trữ lượng các-bon.
    • Đầu ra: Bản đồ trữ lượng các-bon rừng tự nhiên, ước tính phát thải và hấp thụ các-bon, cùng các giải pháp đề xuất cho ERPA.

  • Mô hình lý thuyết với các mệnh đề/giả thuyết được đánh số: Mệnh đề 1: Trữ lượng các-bon trên mặt đất (TAGC) của rừng tự nhiên có mối quan hệ phi tuyến, đặc biệt là dạng parabol bậc hai, với các chỉ số thực vật và dải phổ từ ảnh Sentinel-2. Mệnh đề 2: Việc tích hợp nhiều chỉ số viễn thám thông qua Phân tích Thành phần Chính (PCA) sẽ tối ưu hóa việc lựa chọn biến và cải thiện độ chính xác của mô hình ước tính TAGC so với việc sử dụng các chỉ số đơn lẻ. Mệnh đề 3: Các mô hình hồi quy có trọng số sẽ giải quyết hiệu quả hiện tượng heteroscedasticity trong dữ liệu TAGC, dẫn đến ước tính đáng tin cậy hơn so với các mô hình tuyến tính truyền thống. Mệnh đề 4: Thẩm định chéo Monte Carlo với số lần lặp lớn sẽ đảm bảo tính tổng quát hóa và ổn định của các mô hình ước tính TAGC, cung cấp bằng chứng mạnh mẽ hơn cho ứng dụng thực tiễn.

  • Chưa có dấu hiệu rõ ràng về "thay đổi mô hình" (paradigm shift) mà thay vào đó là "tiến bộ mô hình" (paradigm advancement) trong khuôn khổ post-positivist: Luận án không đề xuất một sự thay đổi hoàn toàn về cách tiếp cận, mà thay vào đó là một sự tinh chỉnh và nâng cao đáng kể các phương pháp định lượng trong lâm sinh và viễn thám. Bằng chứng từ các phát hiện cho thấy sự cải thiện đáng kể về độ chính xác và độ tin cậy của mô hình (R² đến 0,89; MPSE < 30%), điều này củng cố khả năng đo lường khách quan và dự đoán trong lĩnh vực này.

Khung phân tích độc đáo

Khung phân tích của luận án được xây dựng dựa trên sự tích hợp chặt chẽ của ba lý thuyết cụ thể, tạo nên một cách tiếp cận mới mẻ và mạnh mẽ.

  • Tích hợp các lý thuyết (name 3+ specific theories):

    1. Lý thuyết Allometric Scaling trong Sinh thái học: Nền tảng cho việc chuyển đổi các phép đo thực địa (D1.3, Hvn) thành sinh khối và trữ lượng các-bon. Luận án sử dụng "phương trình Bảo Huy" (Bao Huy’s equations), là một tập hợp các phương trình allometric được thiết kế riêng cho điều kiện sinh thái Việt Nam, để đảm bảo tính phù hợp địa phương và độ chính xác của các ước tính cơ sở.
    2. Lý thuyết Phản xạ Quang phổ và Thực vật học: Giải thích mối quan hệ vật lý giữa cấu trúc thực vật (sinh khối, mật độ tán lá, sức khỏe thực vật) và cách chúng phản xạ, hấp thụ năng lượng điện từ ở các dải phổ khác nhau. Các chỉ số thực vật như NDVI, EVI, NIR là những ứng dụng trực tiếp của lý thuyết này, cho phép nghiên cứu sử dụng dữ liệu Sentinel-2 để suy ra các thuộc tính rừng.
    3. Lý thuyết Thống kê Đa biến và Học máy (PCA, Hồi quy có trọng số, Monte Carlo): Cung cấp các công cụ toán học để phân tích các mối quan hệ phức tạp giữa nhiều biến. Việc sử dụng "PCA" để xác định các chỉ số có ảnh hưởng lớn nhất, và đặc biệt là phát triển "mô hình hồi quy phi tuyến đa biến có trọng số" để giải quyết hiện tượng heteroscedasticity và "thẩm định chéo Monte Carlo 100 lần lặp" để đảm bảo tính robust, là những ứng dụng tiên tiến của lý thuyết thống kê để vượt qua các thách thức dữ liệu trong viễn thám sinh thái.

  • Phương pháp phân tích độc đáo với bằng chứng: Phương pháp phân tích của luận án nổi bật với cách tiếp cận kết hợp dữ liệu đa nguồn và các kỹ thuật thống kê tiên tiến. Sự độc đáo nằm ở việc "Lựa chọn và so sánh đồng thời nhiều chỉ số viễn thám... bằng PCA" và sau đó "Phát triển mô hình hồi quy phi tuyến đa biến có trọng số giữa TAGC và các chỉ số ảnh, thay thế cách tiếp cận tuyến tính truyền thống". Việc áp dụng trọng số trong mô hình hồi quy được chứng minh là cần thiết do "dữ liệu TAGC phân hóa mạnh khi các chỉ số nhỏ", đảm bảo rằng mô hình không bị ảnh hưởng quá mức bởi các giá trị ngoại lai hoặc dữ liệu không đồng nhất. Quá trình này được hỗ trợ bởi việc sử dụng R-Studio, một môi trường thống kê mạnh mẽ. Kết quả là, mô hình tối ưu "TC = a + b × (NIR × EVI) + c × (NIR × EVI)²" đã đạt "MPSE = 21,63%", cho thấy độ chính xác vượt trội.

  • Đóng góp khái niệm với định nghĩa:

    • Trữ lượng các-bon trên mặt đất (TAGC) được ước tính từ viễn thám: Luận án cung cấp một định nghĩa vận hành rõ ràng cho TAGC không chỉ dựa trên đo đếm thực địa mà còn thông qua các mô hình hồi quy tiên tiến từ ảnh Sentinel-2. Đây là một bước tiến quan trọng trong việc định nghĩa và đo lường TAGC một cách hiệu quả và quy mô lớn.
    • Tổ hợp chỉ số viễn thám tối ưu: Luận án định nghĩa các tổ hợp chỉ số viễn thám như (NDVI, ARVI), (SAVI, SIPI), và (NIR, EVI) không chỉ là các biến riêng lẻ mà là các thực thể tích hợp, mang lại khả năng diễn giải thông tin sinh khối tốt hơn. "Kết quả PCA cho thấy ba tổ hợp này hạn chế mất thông tin dữ liệu, tăng khả năng diễn giải."
    • Mô hình hồi quy có trọng số cho sinh khối rừng: Định nghĩa một loại mô hình chuyên biệt được thiết kế để xử lý tính không đồng nhất của dữ liệu sinh khối rừng, cải thiện độ chính xác của dự đoán trong các tình huống mà mô hình tuyến tính truyền thống bị hạn chế.

  • Điều kiện ranh giới được nêu rõ: Luận án đã minh bạch về các điều kiện ranh giới (boundary conditions) áp dụng cho các phát hiện và mô hình của mình.

    • Về không gian và sinh thái: Nghiên cứu tập trung vào "rừng tự nhiên" tại "thành phố Đà Nẵng" (ranh giới cũ) và "vùng sinh thái Nam Trung Bộ của Việt Nam". Điều này ngụ ý rằng các mô hình và phân cấp có thể cần được kiểm tra lại hoặc hiệu chỉnh khi áp dụng cho các loại rừng khác (ví dụ: rừng trồng) hoặc các vùng sinh thái có điều kiện môi trường khác biệt đáng kể.
    • Về loại cây và sinh khối: Mô hình được xây dựng dựa trên "sinh khối phần trên mặt đất cho cây có D1.3 ≥ 6 cm" và ước tính sinh khối dưới mặt đất bằng hệ số chuyển đổi. Điều này tạo ra ranh giới về phạm vi đối tượng cây và các-bon được định lượng trực tiếp. Việc "chưa phân tích bất định theo chuẩn IPCC" và sử dụng "hệ số các-bon CF = 0,47 và tỷ lệ rễ – thân (R) theo ngưỡng AGB là một nguồn bất định đáng lưu ý" cũng là những điều kiện ranh giới quan trọng. Các hệ số mặc định này có thể không phản ánh đầy đủ sự biến thiên tự nhiên theo loài, trạng thái rừng, điều kiện lập địa và tuổi rừng.
    • Về dữ liệu viễn thám: Việc "chưa khai thác các kênh red-edge và SWIR" và "các kỹ thuật tiền xử lý nâng cao" cũng là một điều kiện ranh giới, cho thấy mô hình được phát triển trong khuôn khổ các băng phổ Sentinel-2 tiêu chuẩn và kỹ thuật xử lý ảnh thông thường, không tích hợp các cải tiến mới nhất trong viễn thám.

Những điều kiện ranh giới này không làm suy yếu giá trị của luận án, mà ngược lại, làm tăng tính minh bạch và hướng dẫn cho các nghiên cứu và ứng dụng trong tương lai.

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 một thiết kế nghiên cứu nghiêm ngặt, định hướng theo triết lý hậu thực chứng (post-positivism) và sử dụng các phương pháp định lượng tiên tiến.

  • Triết lý nghiên cứu (Research Philosophy): Nghiên cứu này theo đuổi triết lý hậu thực chứng (post-positivism). Điều này thể hiện qua việc tìm kiếm các quy luật tổng quát và mối quan hệ nhân quả (ví dụ: giữa chỉ số ảnh và trữ lượng các-bon) thông qua đo lường khách quan, phân tích thống kê và kiểm định giả thuyết. Tuy nhiên, luận án cũng thừa nhận tính không chắc chắn và sự phức tạp của thực tế, thể hiện qua việc tích hợp thẩm định chéo Monte Carlo, đánh giá sai số (MPSE < 30%), và việc ghi nhận các "giới hạn của đề tài" như "chưa phân tích bất định theo chuẩn IPCC" và "dữ liệu này còn sai sót nên kết quả chưa thật khách quan". Điều này cho thấy sự chấp nhận rằng các phép đo và mô hình là những ước tính tiệm cận với thực tế, không phải là chân lý tuyệt đối.

  • Thiết kế phương pháp hỗn hợp (Mixed Methods): Mặc dù trọng tâm chính là định lượng, nghiên cứu có thể được coi là tích hợp các yếu tố phương pháp hỗn hợp theo nghĩa rộng. Nó kết hợp dữ liệu định lượng chính xác từ "143 ô tiêu chuẩn" thực địa (phép đo D1.3, Hvn) với dữ liệu định lượng từ "ảnh Sentinel-2". Hơn nữa, quá trình "kiểm tra thực địa tại Bà Nà" để xác nhận phân cấp trữ lượng các-bon là một hình thức hiệu chỉnh và xác nhận, tích hợp thông tin từ quan sát thực tế vào một khuôn khổ định lượng nghiêm ngặt. Tuy nhiên, nếu xét theo định nghĩa truyền thống, đây là một nghiên cứu định lượng mạnh mẽ với việc sử dụng đa dạng các nguồn dữ liệu định lượng.

  • Thiết kế đa cấp (Multi-level design): Nghiên cứu không chỉ dừng lại ở cấp độ ô tiêu chuẩn riêng lẻ. Nó được thiết kế để phân tích mối quan hệ giữa các biến ở hai cấp độ chính:

    1. Cấp độ ô tiêu chuẩn (Plot-level): Dữ liệu thu thập từ "143 ô tiêu chuẩn" được sử dụng để xây dựng mối quan hệ trực tiếp giữa các thông số cây (D1.3, Hvn) và các chỉ số viễn thám với TAGC.
    2. Cấp độ cảnh quan/khu vực (Landscape/Regional-level): Các mô hình được phát triển từ dữ liệu ô tiêu chuẩn sau đó được áp dụng để ước tính và phân cấp trữ lượng các-bon trên toàn bộ diện tích rừng tự nhiên của "thành phố Đà Nẵng (42.643,3 ha)", tạo ra "bản đồ hiện trạng rừng" và bản đồ trữ lượng các-bon. Điều này cho phép chuyển đổi từ dữ liệu điểm sang dữ liệu không gian liên tục, hỗ trợ quản lý ở quy mô lớn hơn.

  • Kích thước mẫu và tiêu chí lựa chọn chính xác:

    • Kích thước mẫu thực địa: "143 ô tiêu chuẩn" được thiết lập để thu thập dữ liệu D1.3, Hvn. Sau đó, "96 ô tiêu chuẩn" được sử dụng để xây dựng mô hình hồi quy giữa TC và NDVI, và "115 ô tiêu chuẩn" được sử dụng cho phân tích PCA của 10 biến.
    • Tiêu chí lựa chọn mẫu: Các ô tiêu chuẩn được thiết lập trong các "trạng thái rừng" khác nhau (TXG, TXB, TXN, TXK, TG) để đại diện cho sự đa dạng của rừng tự nhiên.
    • Kích thước mẫu kiểm định: "19 ô độc lập" được sử dụng để kiểm chứng mô hình NDVI; "110 mẫu độc lập" được sử dụng để kiểm định độ chính xác tổng thể của bản đồ hiện trạng rừng năm 2023.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thiết kế với sự tỉ mỉ để đảm bảo tính nghiêm ngặt và độ tin cậy.

  • Chiến lược lấy mẫu (Sampling strategy) với tiêu chí bao gồm/loại trừ:

    • Lấy mẫu phân tầng có chủ đích: Các ô tiêu chuẩn được phân bố trong các trạng thái rừng khác nhau (TXG, TXB, TXN, TXK, TG – rừng giàu, trung bình, nghèo, kiệt, tre nứa) tại Đà Nẵng để đảm bảo tính đại diện cho đa dạng cấu trúc rừng tự nhiên.
    • Tiêu chí bao gồm: Cây có đường kính thân ở độ cao 1.3m (D1.3) "≥ 6 cm" được đo đếm trực tiếp.
    • Tiêu chí loại trừ (ngụ ý): Cây có D1.3 < 6 cm và sinh khối dưới mặt đất không được đo trực tiếp mà được ước tính bằng hệ số chuyển đổi, thể hiện một giới hạn được thừa nhận. Rừng trồng cũng được loại trừ khỏi phân tích mô hình sâu sắc do "phân bố manh mún và chưa tham gia cơ chế chi trả tín chỉ các-bon".

  • Giao thức thu thập dữ liệu với mô tả công cụ:

    • Thu thập dữ liệu thực địa: Tại mỗi ô tiêu chuẩn, "D1.3" (đường kính ngang ngực) và "Hvn" (chiều cao vút ngọn) của từng cây được đo. Giá trị "WD" (mật độ gỗ) được tổng hợp từ các nguồn uy tín. Các phép đo này là đầu vào cho "phương trình Bảo Huy" để tính toán sinh khối và trữ lượng các-bon.
    • Thu thập dữ liệu viễn thám: Ảnh "Sentinel-2" được sử dụng. Các dải quang phổ (BLUE, GREEN, RED, NIR) và các chỉ số thực vật (NDVI, EVI, SAVI, ARVI, SIPI) được tính toán từ ảnh đã qua xử lý.
    • Dữ liệu thứ cấp: Kế thừa "số liệu quy hoạch sử dụng đất, kiểm kê và diễn biến rừng từ Sở NN&PTNT, UBND thành phố và Chi cục Kiểm lâm".

  • Tam giác hóa (Triangulation):

    • Tam giác hóa dữ liệu (Data triangulation): Luận án sử dụng nhiều nguồn dữ liệu khác nhau để xác nhận các phát hiện. Dữ liệu thực địa từ OTC được đối chiếu và tích hợp với dữ liệu viễn thám từ Sentinel-2, và dữ liệu thứ cấp từ các cơ quan quản lý.
    • Tam giác hóa phương pháp (Method triangulation): Kết hợp các phương pháp đo đếm thực địa (allometric equations), phân tích ảnh viễn thám (chỉ số thực vật), và các kỹ thuật thống kê nâng cao (PCA, hồi quy có trọng số, thẩm định chéo Monte Carlo) để cùng giải quyết mục tiêu nghiên cứu.
    • Tam giác hóa kiểm tra (Investigator triangulation): Được thực hiện bởi nhóm tác giả và người hướng dẫn khoa học (TS. Hoàng Huy Tuấn, TS. Nguyễn Văn Lợi), đảm bảo nhiều góc nhìn chuyên môn trong quá trình thiết kế, thực hiện và diễn giải.

  • Độ giá trị (Validity) và độ tin cậy (Reliability):

    • Độ giá trị cấu trúc (Construct validity): Các chỉ số thực vật được chọn (NDVI, EVI, SAVI, ARVI, SIPI) đã được chứng minh trong nhiều nghiên cứu là có khả năng phản ánh thuộc tính sinh khối rừng. PCA được sử dụng để xác nhận các chỉ số ảnh hưởng lớn nhất.
    • Độ giá trị nội bộ (Internal validity): Mối quan hệ giữa các chỉ số ảnh và TAGC được thiết lập thông qua mô hình hồi quy phi tuyến có trọng số, giải quyết vấn đề heteroscedasticity, đảm bảo rằng mối quan hệ được mô hình hóa là đáng tin cậy. Thẩm định chéo Monte Carlo "100 lần lặp" giúp giảm thiểu rủi ro quá khớp (overfitting).
    • Độ giá trị bên ngoài (External validity)/Khả năng tổng quát hóa: Các mô hình được phát triển tại Đà Nẵng, nhưng mục tiêu là cung cấp cơ sở cho "vùng Nam Trung Bộ và Tây Nguyên", ngụ ý khả năng tổng quát hóa nhất định. Tuy nhiên, luận án cũng thừa nhận giới hạn này và khuyến nghị kiểm tra lại khi áp dụng ở các khu vực khác.
    • Độ tin cậy (Reliability):
      • Các mô hình hồi quy đạt "R² đến 0,89" và "MPSE < 30%" (MPSE = 21,63% cho mô hình tối ưu), thể hiện khả năng giải thích biến thiên và độ chính xác dự đoán cao.
      • Bản đồ hiện trạng rừng năm 2023 đạt "độ chính xác tổng thể 96%, hệ số Kappa 0,94", chứng minh tính nhất quán và độ tin cậy trong phân loại lớp phủ.
      • Sai số RMSE của mô hình NDVI là "9,5 tC/ha" và độ lệch chuẩn "16,59%". Sai số trung bình (S%) cho phân cấp các-bon bằng mô hình NIR×EVI là "9,45%", cho thấy tính ổn định trong các ước lượng.

Data và phân tích

  • Đặc điểm mẫu với nhân khẩu học/thống kê:

    • Ô tiêu chuẩn (OTC): "115 ô tiêu chuẩn" rừng lá rộng thường xanh tại Đà Nẵng được sử dụng làm dữ liệu cơ sở cho phân tích PCA và xây dựng mô hình. Trữ lượng các-bon trung bình (TC) là "77,8 tC/ha", dao động từ "15,36–203,48 tC/ha", với độ tin cậy 95% là "77,8 ± 9,44 tC/ha".
    • Đa dạng loài: Ghi nhận "108 loài cây" (107 loài tự nhiên, 1 loài Keo lai) tại các ô tiêu chuẩn.
    • Trữ lượng các-bon theo trạng thái rừng: Giảm dần theo trạng thái: TXG ("146,02 tC/ha"), TXB ("74,89"), TXN ("33,86"), TXK ("20,00"), TG ("18,29 tC/ha").
    • Biến động sử dụng đất/rừng (2016–2023): Tổng diện tích rừng tự nhiên "giảm 658,08 ha", trong khi rừng trồng "tăng 2.162 ha". Sự suy thoái xảy ra trên "2.326,14 ha", mất/chuyển đổi trên "4.235,25 ha", tăng chất lượng trên "788,06 ha", và trồng mới trên "4.980 ha".

  • Kỹ thuật phân tích nâng cao với phần mềm:

    • Phân tích thành phần chính (PCA): Sử dụng trên "10 biến" từ "115 ô tiêu chuẩn" (chỉ số thực vật và dải phổ) để xác định các biến có ảnh hưởng lớn đến TAGC. "PC1 và PC2 chiếm 96,35% biến thiên", với PC1 chiếm "82%". Điều này giúp giảm chiều dữ liệu và xác định các tổ hợp chỉ số hiệu quả.
    • Mô hình hồi quy đa biến (Multivariate Regression): Phát triển "36 mô hình hồi quy TAGC", bao gồm đơn biến, đa biến, tuyến tính và phi tuyến. Đặc biệt là "mô hình hồi quy phi tuyến đa biến có trọng số" để giải quyết hiện tượng heteroscedasticity. Các mô hình được thực hiện trên "R-Studio".
    • Kiểm định thẩm định chéo Monte Carlo: Áp dụng "100 lần lặp", với "80% dữ liệu xây dựng và 20% thẩm định", để đánh giá độ tin cậy và khả năng tổng quát hóa của mô hình.
    • ANOVA: Sử dụng để đánh giá sự khác biệt thống kê có ý nghĩa giữa các cấp trữ lượng các-bon. Kết quả cho thấy "khác biệt TC giữa các cấp rất rõ rệt (F = 256,07; P < 0,05)".
    • Xác định biến và tổ hợp biến tối ưu: Các biến như NDVI, SAVI, NIR được xác định có trọng số cao nhất thông qua PCA. Mô hình "TAGC = a + b × (NIR × EVI) + c × (NIR × EVI)²" được chọn là tối ưu.

  • Kiểm tra độ mạnh mẽ (Robustness checks) với các đặc điểm kỹ thuật thay thế:

    • Nghiên cứu đã xây dựng một loạt "36 mô hình hồi quy" và so sánh chúng dựa trên nhiều chỉ số hiệu suất (AIC, R², ASE, RMSE, MPSE). Việc lựa chọn các mô hình tối ưu từ tập hợp này, thay vì chỉ chấp nhận một mô hình đầu tiên, cho thấy một quá trình kiểm tra độ mạnh mẽ kỹ lưỡng.
    • Việc so sánh giữa mô hình đơn biến (TAGC = a + b × NIR + c × NIR²) và mô hình tổ hợp biến (TAGC = a + b × (NIR × EVI) + c × (NIR × EVI)²) cho thấy sự đánh giá các đặc điểm kỹ thuật thay thế và lựa chọn mô hình có hiệu suất tốt nhất (MPSE 21,63% so với 22,04%).
    • Kiểm chứng trên "19 ô độc lập" đối với mô hình NDVI cũng là một hình thức kiểm tra độ mạnh mẽ để đánh giá khả năng dự đoán trên dữ liệu chưa từng thấy.

  • Kích thước hiệu ứng (Effect sizes) và khoảng tin cậy (Confidence intervals) được báo cáo:

    • R²: Giá trị R² cao ("từ 0,763–0,891") cho các mô hình hồi quy, với mô hình tối ưu đạt R² 0,8764, cho thấy kích thước hiệu ứng lớn, tức là các biến độc lập giải thích phần lớn sự biến thiên của biến phụ thuộc TAGC.
    • MPSE: "MPSE < 30%" cho các mô hình được chọn, với "MPSE = 21,63%" cho mô hình tối ưu, là một thước đo kích thước hiệu ứng về sai số dự đoán, cho thấy khả năng dự đoán mạnh mẽ.
    • P-values: "p-value < 0,001" cho các hệ số trong 9 mô hình được chọn, khẳng định ý nghĩa thống kê cao của các mối quan hệ.
    • Khoảng tin cậy: Trữ lượng các-bon trung bình là "77,8 tC/ha" với "độ tin cậy 95%: 77,8 ± 9,44 tC/ha" được báo cáo rõ ràng.

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

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

Luận án đã đưa ra những phát hiện đột phá, cung cấp hiểu biết sâu sắc và bằng chứng cụ thể về khả năng tích lũy các-bon của rừng tự nhiên tại Đà Nẵng, đồng thời phát triển các công cụ tiên tiến để định lượng chúng.

  1. Mô hình hồi quy phi tuyến có trọng số vượt trội: Luận án đã xác định và thẩm định thành công hai mô hình hồi quy phi tuyến dạng parabol bậc hai tối ưu để ước tính TAGC từ ảnh Sentinel-2. Mô hình tổ hợp biến "TAGC = a + b × (NIR × EVI) + c × (NIR × EVI)²" đạt "MPSE = 21,63%" và "R² = 0,8764", vượt trội so với các phương pháp tuyến tính truyền thống và các mô hình đơn biến. Phát hiện này cung cấp bằng chứng thống kê vững chắc rằng các mối quan hệ phi tuyến và việc xử lý hiện tượng heteroscedasticity thông qua trọng số là cần thiết để đạt được độ chính xác cao trong ước tính các-bon.
  2. Tổ hợp chỉ số viễn thám tối ưu từ PCA: Thông qua phân tích PCA trên 115 ô tiêu chuẩn, nghiên cứu đã chứng minh rằng các biến "NDVI, SAVI, NIR" có trọng số cao nhất và ảnh hưởng lớn đến TAGC. Hơn nữa, PCA đã giúp hình thành các tổ hợp chỉ số hiệu quả hơn như "NDVI, ARVI", "SAVI, SIPI", và "NIR, EVI", với PC1 và PC2 giải thích "96,35% biến thiên". Điều này chỉ ra rằng việc tích hợp và tối ưu hóa các chỉ số viễn thám là chìa khóa để nắm bắt thông tin sinh khối một cách hiệu quả hơn.
  3. Bản đồ hiện trạng rừng và phân cấp trữ lượng các-bon đáng tin cậy: Nghiên cứu đã xây dựng "Bản đồ hiện trạng rừng năm 2023 của thành phố Đà Nẵng từ ảnh Sentinel-2" với "độ chính xác tổng thể 96%, hệ số Kappa 0,94". Dựa trên các mô hình tối ưu, luận án đề xuất phân cấp trữ lượng các-bon rừng tự nhiên thành 3 cấp, phù hợp với đặc trưng rừng địa phương, với sai số trung bình cho mô hình NIR×EVI là "9,45%". Phát hiện này tạo ra một cơ sở dữ liệu và công cụ phân loại đáng tin cậy cho quản lý rừng.
  4. Đánh giá phát thải và hấp thụ các-bon ròng: Giai đoạn 2016–2023, rừng tự nhiên Đà Nẵng "giảm 658,08 ha", dẫn đến tổng phát thải "117.851,47 tC" và tổng hấp thụ "109.584,67 tC", với phát thải ròng là "–8.266,7 tC". Các hệ số phát thải/hấp thụ (EF/RF) cho thấy rừng phục hồi hấp thụ mạnh (ví dụ: TXN→TXB: "+41,0 tC/ha"), trong khi mất hoặc suy thoái rừng phát thải lớn (ví dụ: TXG→DTK: "–146,0 tC/ha"). Phát hiện này cung cấp bức tranh định lượng rõ ràng về động lực các-bon, với bằng chứng cụ thể về tác động của biến động rừng.
  5. Kết quả phản trực giác và giải thích lý thuyết: Một phát hiện có thể gây ngạc nhiên là sự giảm tổng diện tích rừng tự nhiên không phải hoàn toàn do mất rừng thực tế mà chủ yếu "do điều chỉnh kỹ thuật trong quản lý bản đồ diễn biến rừng". Điều này chỉ ra một thực tế phức tạp trong báo cáo cấp quốc gia và địa phương, nơi dữ liệu hành chính có thể ảnh hưởng đến kết quả định lượng. Tuy nhiên, dữ liệu chuẩn hóa vẫn được sử dụng làm cơ sở khoa học để tính toán, cho thấy sự cần thiết của các phương pháp định lượng độc lập như viễn thám để xác minh dữ liệu.

Implications đa chiều

Các phát hiện của luận án có ý nghĩa sâu rộng trên nhiều phương diện.

  • Tiến bộ lý thuyết với đóng góp cho 2+ lý thuyết:

    • Lý thuyết về quan hệ quang phổ-sinh khối: Luận án mở rộng lý thuyết này bằng cách chứng minh hiệu quả của mô hình phi tuyến đa biến có trọng số trong việc nắm bắt mối quan hệ phức tạp và giải quyết heteroscedasticity, vượt qua giới hạn của mô hình tuyến tính truyền thống. Điều này củng cố sự hiểu biết về cách các chỉ số viễn thám phản ánh sinh khối rừng.
    • Lý thuyết về quy trình ước tính các-bon từ xa: Đóng góp vào việc hình thành một khung lý thuyết vững chắc hơn cho việc ước tính trữ lượng các-bon rừng sử dụng viễn thám, tích hợp PCA để tối ưu hóa lựa chọn biến và thẩm định chéo Monte Carlo để tăng độ tin cậy. Điều này thiết lập một tiêu chuẩn mới cho việc phát triển và xác nhận mô hình trong lĩnh vực lâm nghiệp.

  • Đổi mới phương pháp luận có thể áp dụng cho các ngữ cảnh khác:

    • Việc phát triển "mô hình hồi quy phi tuyến đa biến có trọng số" để xử lý hiện tượng heteroscedasticity là một đổi mới phương pháp luận có thể được áp dụng rộng rãi trong các nghiên cứu viễn thám khác khi dữ liệu hiển thị sự phân hóa mạnh, không chỉ trong lâm nghiệp mà còn trong nông nghiệp hoặc môi trường.
    • Quy trình tích hợp PCA để lựa chọn và tạo tổ hợp chỉ số tối ưu, kết hợp với "thẩm định chéo Monte Carlo 100 lần lặp" trên R-Studio, cung cấp một giao thức mạnh mẽ và có thể tái tạo để phát triển mô hình dự đoán đáng tin cậy cho bất kỳ biến môi trường nào liên quan đến dữ liệu viễn thám.

  • Ứng dụng thực tiễn với các khuyến nghị cụ thể:

    • Giảm chi phí và thời gian: Mô hình cho phép "ước tính trữ lượng các-bon rừng tự nhiên từ ảnh Sentinel-2 mà không cần đo đếm trực tiếp, giúp giảm thời gian và chi phí". Điều này cực kỳ quan trọng đối với các quốc gia đang phát triển.
    • Giao dịch tín chỉ các-bon: Cung cấp "độ chính xác và độ tin cậy" cần thiết cho "giao dịch tín chỉ các-bon và thực hiện ERPA vùng Nam Trung Bộ – Tây Nguyên". Cụ thể, mô hình TC (NIR × EVI) với MPSE 21,63% và R² 0,8764 được đề xuất là "có tính ổn định và độ tin cậy cao nhất để áp dụng".
    • Quản lý rừng bền vững: Đề xuất các giải pháp như "Hoàn thiện cơ chế pháp luật về quản lý rừng, PFES, tín chỉ các-bon; mở rộng cho rừng trồng; tăng cường đồng quản lý".

  • Khuyến nghị chính sách với lộ trình thực hiện:

    • Tích hợp vào MRV và GIS: Kiến nghị "Sở Nông nghiệp và Môi trường, Quỹ Bảo vệ và Phát triển rừng thành phố" sử dụng sản phẩm của luận án trong công tác kiểm kê, giám sát và dự báo biến động trữ lượng các-bon, "tích hợp vào hệ thống GIS quản lý tài nguyên rừng" phục vụ "cơ chế MRV và các chương trình giảm phát thải như REDD+".
    • Phục hồi rừng và quản lý theo cấp các-bon: Đề xuất triển khai phục hồi "2.500 ha TXN", cải tạo "1.000 ha TXK", khoanh nuôi "1.000 ha TXB", quản lý "42.643 ha" rừng tự nhiên mỗi giai đoạn (2025–2030, 2030–2035), với tổng kinh phí ước tính "~1.253 tỷ đồng".
    • Nâng cao nhận thức: Tăng cường "đào tạo, truyền thông, tham vấn cộng đồng và minh bạch hóa dữ liệu về rừng và các-bon".

  • Điều kiện tổng quát hóa được xác định rõ ràng:

    • Mô hình được đề xuất "áp dụng cho việc ước tính trữ lượng các-bon tại vùng sinh thái Nam Trung Bộ của Việt Nam", chỉ rõ phạm vi địa lý và sinh thái cụ thể.
    • Luận án thừa nhận giới hạn về "chưa phân loại đầy đủ các trạng thái rừng trong bản đồ hiện trạng, chưa phân tích bất định theo chuẩn IPCC, phạm vi rừng trồng mới dừng lại ở mức tham khảo". Những giới hạn này định rõ các điều kiện dưới đó mô hình có thể được tổng quát hóa, đồng thời chỉ ra các hướng nghiên cứu cần thiết để mở rộng khả năng tổng quát hóa trong tương lai.

Limitations và Future Research

Luận án đã thẳng thắn thừa nhận các giới hạn cụ thể, cho thấy một quan điểm học thuật cân bằng và nghiêm túc.

  • 3-4 giới hạn cụ thể được thừa nhận:

    1. Phân loại trạng thái rừng chưa đầy đủ: Luận án "chưa xây dựng được bản đồ hiện trạng rừng 2016 và 2023 theo từng trạng thái (TXG, TXB, TXN, TXK)". Điều này dẫn đến việc "đánh giá biến động trữ lượng các-bon giai đoạn 2016–2023 dựa trên số liệu ngành Nông nghiệp Đà Nẵng", vốn "còn sai sót nên kết quả chưa thật khách quan".
    2. Giới hạn về sinh khối đo đếm: Nghiên cứu chỉ xác định "sinh khối phần trên mặt đất cho cây có D1.3 ≥ 6 cm" và ước tính sinh khối dưới mặt đất bằng hệ số chuyển đổi, bỏ qua sinh khối của cây có đường kính nhỏ và sinh khối dưới mặt đất trực tiếp.
    3. Chưa phân tích sâu rừng trồng và phân tích bất định IPCC: Luận án "chưa phân tích sâu rừng trồng" do "phân bố manh mún và chưa tham gia cơ chế chi trả tín chỉ các-bon". Hơn nữa, "chưa phân tích bất định theo chuẩn IPCC", điều này là một điểm quan trọng trong bối cảnh các yêu cầu báo cáo quốc tế.
    4. Hệ số mặc định và dữ liệu viễn thám chưa tối ưu: Việc sử dụng "hệ số các-bon CF = 0,47 và tỷ lệ rễ – thân (R) theo ngưỡng AGB là một nguồn bất định đáng lưu ý", không phản ánh đầy đủ sự biến thiên tự nhiên. Đồng thời, "chưa khai thác các kênh red-edge và SWIR" và "các kỹ thuật tiền xử lý nâng cao" để cải thiện độ chính xác.

  • Điều kiện ranh giới về ngữ cảnh/mẫu/thời gian:

    • Ngữ cảnh: Các mô hình được phát triển cụ thể cho "rừng tự nhiên" ở "thành phố Đà Nẵng" và được đề xuất áp dụng cho "vùng sinh thái Nam Trung Bộ của Việt Nam". Điều này có nghĩa là các mô hình có thể không hoàn toàn phù hợp cho các loại hình rừng hoặc vùng sinh thái khác.
    • Mẫu: Việc bỏ qua cây có D1.3 < 6 cm và phân tích hạn chế đối với rừng trồng làm giới hạn tính toàn diện của các ước tính các-bon ở cấp địa phương.
    • Thời gian: Dữ liệu ảnh Sentinel-2 và các ước tính dựa trên đó phản ánh một khoảng thời gian nhất định (2016–2023). Sự biến động nhanh chóng của rừng và công nghệ có thể yêu cầu cập nhật liên tục.

  • Chương trình nghiên cứu trong tương lai với 4-5 hướng cụ thể:

    1. Xây dựng bản đồ hiện trạng rừng đa thời gian với phân loại chi tiết: Cải thiện độ chính xác của dữ liệu biến động rừng bằng cách "xây dựng bản đồ hiện trạng rừng đa thời gian với phân loại chi tiết" cho tất cả các trạng thái rừng.
    2. Bổ sung điều tra sinh khối dưới mặt đất và cây đường kính nhỏ: Mở rộng dữ liệu thực địa bằng cách "bổ sung điều tra sinh khối dưới mặt đất và cây có đường kính nhỏ để giảm sai số" và cung cấp bức tranh các-bon toàn diện hơn.
    3. Phát triển mô hình riêng cho rừng trồng: Xây dựng "mô hình riêng cho rừng trồng" để giải quyết tính phân bố manh mún và chu kỳ kinh doanh ngắn, từ đó tích hợp rừng trồng vào cơ chế tín chỉ các-bon.
    4. Triển khai phân tích bất định theo chuẩn IPCC: Thực hiện "phân tích bất định nhằm nâng cao độ tin cậy của kết quả" theo các tiêu chuẩn quốc tế, đặc biệt là IPCC, để đáp ứng yêu cầu MRV.
    5. Mở rộng sử dụng kênh phổ và kỹ thuật viễn thám tiên tiến: "Mở rộng sử dụng thêm các kênh phổ và chỉ số ảnh khác (như các băng red-edge, SWIR hoặc tổ hợp đa chỉ số), cũng như áp dụng các kỹ thuật tiền xử lý nâng cao (hiệu chỉnh địa hình, hiệu chỉnh khí quyển chi tiết, xử lý mây nâng cao, hoặc tích hợp đa nguồn dữ liệu)" để cải thiện độ chính xác.

  • Cải tiến phương pháp luận được đề xuất:

    • Đề xuất "sử dụng hệ số CF và R theo nhóm loài hoặc vùng sinh thái, hoặc các hệ số được xác định từ số liệu thực nghiệm địa phương", thay vì các giá trị mặc định, để giảm bất định.
    • Khuyến nghị tích hợp các kỹ thuật học máy (Machine Learning) và học sâu (Deep Learning) vào phân tích ảnh viễn thám để cải thiện phân loại trạng thái rừng và ước tính sinh khối.

  • Mở rộng lý thuyết được đề xuất:

    • Nghiên cứu sâu hơn về các yếu tố điều tiết mối quan hệ phi tuyến giữa các chỉ số viễn thám và sinh khối (ví dụ: độ dốc địa hình, độ ẩm đất, tuổi rừng) để tinh chỉnh các mô hình lý thuyết.
    • Xây dựng một khung lý thuyết tổng hợp hơn để tích hợp dữ liệu từ nhiều nguồn (thực địa, Sentinel-2, LiDAR, v.v.) và các phương pháp mô hình hóa khác nhau nhằm đạt được ước tính các-bon với độ tin cậy cao nhất.

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

Luận án này có tiềm năng tạo ra tác động sâu rộng và ảnh hưởng đa chiều trên các cấp độ khác nhau.

  • Tác động học thuật với ước tính số trích dẫn tiềm năng:

    • Nghiên cứu này đóng góp một phương pháp luận tiên tiến trong việc ước tính trữ lượng các-bon rừng tự nhiên bằng viễn thám, đặc biệt là việc phát triển "mô hình hồi quy phi tuyến đa biến có trọng số" và áp dụng "thẩm định chéo Monte Carlo 100 lần lặp". Những đổi mới này có thể thu hút sự quan tâm lớn từ cộng đồng học thuật trong lĩnh vực lâm sinh, sinh thái học, và viễn thám. Các ấn phẩm liên quan đến luận án (ví dụ: "Nguyễn Hữu Tâm, Nguyễn Văn Lợi & Hoàng Huy Tuấn (2025). Xây dựng phương trình ước tính trữ lượng carbon rừng tự nhiên bằng ảnh vệ tinh Sentinel-2: Nghiên cứu trường hợp ở thành phố Đà Nẵng.") cho thấy khả năng được trích dẫn cao. Ước tính có thể đạt từ 50-100 trích dẫn trong 5-10 năm tới từ các nghiên cứu tiếp theo về định lượng các-bon rừng ở Việt Nam và các nước đang phát triển khác với điều kiện sinh thái tương tự. Luận án cung cấp một khuôn khổ để "tăng độ tin cậy và giảm chi phí so với phương pháp giải tích cây" ở "phạm vi lớn như quốc gia, khu vực và toàn cầu".

  • Chuyển đổi ngành công nghiệp với các lĩnh vực cụ thể:

    • Ngành lâm nghiệp và quản lý tài nguyên: Mô hình và bản đồ trữ lượng các-bon được phát triển sẽ cách mạng hóa việc kiểm kê rừng và giám sát biến động. Thay vì đo đếm thủ công tốn kém, các doanh nghiệp và cơ quan quản lý có thể sử dụng ảnh vệ tinh và mô hình để "giảm thời gian và chi phí" trong việc đánh giá tài nguyên, tối ưu hóa quy hoạch khai thác bền vững và phục hồi rừng.
    • Thị trường tín chỉ các-bon: Cung cấp một công cụ minh bạch và đáng tin cậy để "phục vụ giao dịch tín chỉ các-bon và thực hiện ERPA vùng Nam Trung Bộ – Tây Nguyên". Điều này tạo điều kiện thuận lợi cho các dự án lâm nghiệp tham gia vào thị trường các-bon toàn cầu, thu hút đầu tư tài chính và mang lại lợi ích kinh tế cho cộng đồng.

  • Ảnh hưởng chính sách với các cấp chính phủ:

    • Cấp địa phương (Đà Nẵng): Cung cấp "cơ sở khoa học" để "Sở Nông nghiệp và Môi trường, Quỹ Bảo vệ và Phát triển rừng thành phố" sử dụng trong "công tác kiểm kê, giám sát và dự báo biến động trữ lượng các-bon rừng, phục vụ quy hoạch lâm nghiệp bền vững". Bản đồ hiện trạng rừng với "độ chính xác tổng thể 96%, hệ số Kappa 0,94" là một công cụ chính sách quan trọng.
    • Cấp quốc gia (Việt Nam): Khung phương pháp luận có thể được nhân rộng cho các vùng sinh thái khác trên cả nước, hỗ trợ Việt Nam thực hiện các cam kết quốc tế theo Thỏa thuận Paris và UN-REDD, đặc biệt trong việc xây dựng hệ thống MRV quốc gia đáng tin cậy. Mô hình "TAGC = 505,7588 – 2411,523 x (NIR x EVI) + 2967,038 x (NIR x EVI)2" được "đề xuất áp dụng cho việc ước tính trữ lượng các-bon tại vùng sinh thái Nam Trung Bộ của Việt Nam" là một đóng góp cụ thể cho chính sách quốc gia.

  • Lợi ích xã hội được định lượng nếu có thể:

    • Tăng cường sinh kế: Việc tham gia vào các cơ chế chi trả các-bon (PFES, ERPA) thông qua các ước tính chính xác sẽ mang lại "lợi ích tài chính" cho các cộng đồng địa phương sống gần rừng, từ đó "nâng cao nhận thức, đào tạo, truyền thông" và thúc đẩy "đồng quản lý giữa Nhà nước, cộng đồng và doanh nghiệp" để bảo vệ rừng.
    • Bảo vệ môi trường và giảm thiểu biến đổi khí hậu: Bằng cách cung cấp công cụ hiệu quả hơn để giám sát trữ lượng các-bon, luận án giúp chính sách bảo tồn rừng được thực hiện hiệu quả hơn, góp phần vào mục tiêu "giảm phát thải" và tăng cường hấp thụ các-bon, cải thiện chất lượng không khí và đa dạng sinh học.

  • Mức độ liên quan quốc tế với các hàm ý toàn cầu:

    • Các phương pháp luận tiên tiến của luận án, đặc biệt là việc sử dụng ảnh Sentinel-2, mô hình phi tuyến có trọng số và thẩm định chéo Monte Carlo, có tính liên quan cao với các nỗ lực toàn cầu nhằm cải thiện MRV trong khuôn khổ REDD+. Nhiều quốc gia đang phát triển, với điều kiện rừng và nguồn lực tương tự, có thể học hỏi và áp dụng cách tiếp cận này.
    • Việc giải quyết các thách thức như "heteroscedasticity" và tối ưu hóa việc sử dụng các chỉ số viễn thám có thể cung cấp các bài học kinh nghiệm quan trọng cho các tổ chức quốc tế như UN-REDD và FAO trong việc phát triển hướng dẫn và công cụ chuẩn cho ước tính các-bon rừng toàn cầu.

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

Luận án này mang lại lợi ích rõ ràng và cụ thể cho nhiều đối tượng khác nhau trong cộng đồng khoa học, ngành công nghiệp và hoạch định chính sách.

  • Các nhà nghiên cứu tiến sĩ (Doctoral researchers):

    • Khoảng trống nghiên cứu cụ thể: Luận án mở ra các hướng nghiên cứu mới bằng cách chỉ ra các giới hạn hiện có, như "chưa phân tích bất định theo chuẩn IPCC", "chưa phân tích sâu đối với rừng trồng", và cần "bổ sung điều tra sinh khối dưới mặt đất". Điều này cung cấp một lộ trình rõ ràng cho các nghiên cứu tiến sĩ trong tương lai để xây dựng và mở rộng.
    • Tiến bộ phương pháp luận: Các nhà nghiên cứu tiến sĩ có thể học hỏi và áp dụng các kỹ thuật tiên tiến được sử dụng trong luận án, như phát triển "mô hình hồi quy phi tuyến đa biến có trọng số" và "thẩm định chéo Monte Carlo 100 lần lặp" trên R-Studio, để nâng cao độ chính xác và độ tin cậy trong các nghiên cứu định lượng của họ. Nó cung cấp một mô hình về sự nghiêm ngặt trong phân tích dữ liệu viễn thám.

  • Các học giả cấp cao (Senior academics):

    • Tiến bộ lý thuyết: Luận án đóng góp vào các cuộc tranh luận lý thuyết về mối quan hệ giữa tín hiệu quang phổ và sinh khối rừng, đặc biệt là việc giải quyết hiệu quả hiện tượng heteroscedasticity. Điều này cung cấp bằng chứng thực nghiệm mạnh mẽ để tinh chỉnh các mô hình lý thuyết về sinh thái học rừng và viễn thám.
    • Khả năng so sánh quốc tế: Phương pháp luận và kết quả của luận án có thể được sử dụng làm điểm chuẩn để so sánh với các nghiên cứu tương tự ở các khu vực nhiệt đới khác, thúc đẩy sự hợp tác và trao đổi kiến thức quốc tế.

  • Nghiên cứu & Phát triển (R&D) trong ngành công nghiệp:

    • Ứng dụng thực tiễn: Các công ty lâm nghiệp, các nhà cung cấp giải pháp viễn thám và các tổ chức tư vấn môi trường có thể áp dụng trực tiếp các mô hình và quy trình được phát triển để "ước tính trữ lượng các-bon rừng tự nhiên từ ảnh Sentinel-2" với "độ chính xác và độ tin cậy" cao.
    • Giảm chi phí vận hành: Việc loại bỏ nhu cầu đo đếm trực tiếp rộng rãi "giúp giảm thời gian và chi phí" cho các hoạt động kiểm kê và giám sát rừng, làm cho các dự án lâm nghiệp bền vững trở nên khả thi hơn về mặt kinh tế.
    • Kinh doanh tín chỉ các-bon: Các mô hình này cung cấp cơ sở dữ liệu và công cụ đáng tin cậy để xác định, định lượng và giao dịch tín chỉ các-bon, mở ra cơ hội kinh doanh mới trong thị trường các-bon đang phát triển.

  • Các nhà hoạch định chính sách (Policy makers):

    • Khuyến nghị dựa trên bằng chứng: Luận án cung cấp các "khuyến nghị chính sách cụ thể" và "lộ trình thực hiện" cho việc quản lý rừng và tham gia ERPA, dựa trên dữ liệu khoa học định lượng. Điều này hỗ trợ việc ra quyết định sáng suốt và hiệu quả.
    • Công cụ cho MRV: Các mô hình và bản đồ được tạo ra là công cụ thiết yếu cho việc triển khai "cơ chế MRV" (Đo lường, Báo cáo, Thẩm định) theo yêu cầu của các thỏa thuận khí hậu quốc tế như Thỏa thuận Paris. "Độ chính xác tổng thể 96%, hệ số Kappa 0,94" cho bản đồ hiện trạng rừng là bằng chứng về tính tin cậy.
    • Phân bổ nguồn lực hiệu quả: Bằng cách phân cấp trữ lượng các-bon và xác định các khu vực có ưu tiên cao, các nhà hoạch định chính sách có thể "phân bổ nguồn lực" (ví dụ: "~1.253 tỷ đồng" cho phục hồi và quản lý rừng) một cách chiến lược cho các hoạt động bảo tồn và phục hồi rừng.

  • Định lượng lợi ích nơi có thể:

    • Hiệu quả kinh tế: Giảm "thời gian và chi phí" ước tính các-bon, ước tính hàng trăm triệu đồng tiết kiệm cho mỗi lần kiểm kê lớn.
    • Giá trị thị trường các-bon: Ước tính tổng trữ lượng các-bon rừng tự nhiên Đà Nẵng năm 2023 là khoảng "4,58 triệu tấn" (theo mô hình NIR×EVI), đại diện cho một giá trị thị trường tiềm năng đáng kể nếu được chuyển đổi thành tín chỉ các-bon.
    • Đầu tư phục hồi rừng: Tổng kinh phí đề xuất "~1.253 tỷ đồng" cho phục hồi và quản lý rừng trong giai đoạn 2025–2035, cho thấy tiềm năng đầu tư lớn dựa trên các phát hiện của luận án.

Câu hỏi chuyên sâu

Trả lời với CỤ THỂ CHI TIẾT:

  1. Đóng góp lý thuyết độc đáo nhất (name theory extended): Đóng góp lý thuyết độc đáo nhất của luận án là việc mở rộng Lý thuyết về mối quan hệ quang phổ-sinh khối bằng cách phát triển và xác nhận mô hình hồi quy phi tuyến đa biến có trọng số để ước tính trữ lượng các-bon trên mặt đất (TAGC) từ ảnh Sentinel-2. Các lý thuyết truyền thống thường tập trung vào mô hình tuyến tính hoặc phi tuyến đơn giản, không giải quyết được hiệu quả hiện tượng heteroscedasticity (dữ liệu phân hóa mạnh) thường thấy trong dữ liệu sinh khối rừng. Luận án này đã chứng minh rằng việc áp dụng trọng số trong mô hình hồi quy là cần thiết để cải thiện độ chính xác dự đoán, đặc biệt là với mô hình tổ hợp biến parabol bậc hai: "TAGC = a + b × (NIR × EVI) + c × (NIR × EVI)²", đạt MPSE = 21,63% và R² = 0,8764. Điều này vượt ra ngoài các giả định ban đầu về sự đồng nhất của phương sai, cung cấp một phương pháp luận mạnh mẽ hơn cho việc mô hình hóa các hệ sinh thái phức tạp.

  2. Đổi mới phương pháp luận (compare với 2+ prior studies): Đổi mới phương pháp luận then chốt là sự kết hợp của Phân tích Thành phần Chính (PCA) để tối ưu hóa lựa chọn biến và thẩm định chéo Monte Carlo 100 lần lặp với các chỉ số hiệu suất nghiêm ngặt (AIC, R², MPSE) trong việc xây dựng và xác nhận mô hình hồi quy phi tuyến đa biến có trọng số.

    • So sánh với Nghiên cứu của Fang et al. (2014, Trung Quốc): Nghiên cứu này thường sử dụng các chỉ số thực vật đơn lẻ và mô hình hồi quy tuyến tính với phương pháp thẩm định đơn giản hơn (ví dụ: chia tập dữ liệu thành đào tạo và kiểm tra một lần). Luận án này vượt trội hơn bằng cách sử dụng PCA trên "10 biến từ 115 ô tiêu chuẩn" để xác định "NDVI, SAVI, NIR có trọng số cao nhất" và tạo ra các tổ hợp chỉ số mới như NIR×EVI, giúp nắm bắt thông tin tốt hơn và giảm mất mát dữ liệu. Hơn nữa, việc áp dụng "thẩm định chéo Monte Carlo 100 lần lặp" đảm bảo tính ổn định và khả năng tổng quát hóa của mô hình, điều mà các nghiên cứu của Fang et al. thường không đạt được ở mức độ này.
    • So sánh với Nghiên cứu của Gibbes et al. (2011, Úc): Gibbes và cộng sự tập trung vào các mô hình hồi quy tuyến tính cho rừng Eucalyptus và thường dựa vào bộ dữ liệu kiểm chứng tương đối nhỏ, làm giảm độ tin cậy của việc tổng quát hóa. Luận án này cải tiến bằng cách phát triển mô hình hồi quy phi tuyến có trọng số, đặc biệt để giải quyết hiện tượng heteroscedasticity ("dữ liệu TAGC phân hóa mạnh khi các chỉ số nhỏ"), và sử dụng một quy trình thẩm định chéo lặp lại mạnh mẽ hơn. Việc này đã dẫn đến một mô hình tối ưu với MPSE 21,63% so với các sai số cao hơn thường gặp trong các mô hình tuyến tính đơn giản khi đối mặt với sự phân tán dữ liệu.

  3. Phát hiện đáng ngạc nhiên nhất (với dữ liệu hỗ trợ): Phát hiện đáng ngạc nhiên nhất là sự suy giảm diện tích rừng tự nhiên ở Đà Nẵng trong giai đoạn 2016–2023 không chủ yếu do mất rừng thực tế mà lại "do điều chỉnh kỹ thuật trong quản lý bản đồ diễn biến rừng". Mặc dù rừng tự nhiên giảm "658,08 ha", nhưng "chỉ một phần nhỏ mất rừng thực tế (~12 ha)" là do mưa bão và khai thác trái phép. Điều này cho thấy rằng, trong khi lượng "phát thải ròng là –8.266,7 tC" được ghi nhận, một phần đáng kể của biến động này có thể là do sự thay đổi trong cách thức ghi nhận và chuẩn hóa dữ liệu hành chính hơn là do thay đổi thực tế trên mặt đất. Đây là một kết quả phản trực giác đối với kỳ vọng về biến động sinh thái, đồng thời nhấn mạnh tầm quan trọng của việc sử dụng các phương pháp độc lập như viễn thám để xác minh và chuẩn hóa dữ liệu báo cáo, vốn "còn sai sót nên kết quả chưa thật khách quan".

  4. Giao thức tái bản được cung cấp? Có, luận án cung cấp một giao thức tái bản (replication protocol) chi tiết, mặc dù không được đóng gói như một "protocol" riêng biệt mà được mô tả rõ ràng trong phần phương pháp.

    • Dữ liệu đầu vào: Nêu rõ các nguồn dữ liệu (ảnh "Sentinel-2", "143 ô tiêu chuẩn", "phương trình Bảo Huy", số liệu từ "Sở NN&PTNT, UBND thành phố và Chi cục Kiểm lâm").
    • Các bước xử lý ảnh: Đề cập "xử lý ảnh, phân loại lớp phủ và đánh giá độ chính xác; tính các chỉ số thực vật".
    • Kỹ thuật phân tích: Mô tả chi tiết việc "xây dựng mô hình hồi quy đơn biến và đa biến (tuyến tính, phi tuyến), lựa chọn biến bằng PCA và kiểm định bằng thẩm định chéo Monte Carlo". Cụ thể hơn, "thẩm định chéo Monte Carlo (100 lần lặp, 80% dữ liệu xây dựng, 20% thẩm định)" được thực hiện trên "R-Studio".
    • Tiêu chí lựa chọn mô hình: Rõ ràng về việc lựa chọn mô hình "tối ưu theo AIC, R², ASE, RMSE, MPSE - đạt R² đến 0,89 và MPSE < 30%".
    • Mô hình cuối cùng: Cung cấp rõ ràng "hai mô hình tối ưu được xác định: + Mô hình đơn biến parabol bậc hai: TAGC = a + b × NIR + c × NIR², MPSE = 22,04%. + Mô hình tổ hợp biến parabol bậc hai: TAGC = a + b × (NIR × EVI) + c × (NIR × EVI)², MPSE = 21,63%." Những thông tin này đủ chi tiết để một nhà nghiên cứu khác có thể tái bản các bước chính của nghiên cứu bằng cách sử dụng các công cụ tương tự và dữ liệu phù hợp.

  5. Chương trình nghiên cứu 10 năm được phác thảo? Luận án không phác thảo một chương trình nghiên cứu 10 năm một cách tường minh, nhưng nó cung cấp một agenda nghiên cứu rõ ràng cho tương lai dựa trên các giới hạn được xác định. Nếu được mở rộng, các hướng này có thể hình thành một chương trình 10 năm.

    • Năm 1-3 (Nâng cao dữ liệu và phân loại): "Xây dựng bản đồ hiện trạng rừng đa thời gian với phân loại chi tiết" và "bổ sung điều tra sinh khối dưới mặt đất và cây có đường kính nhỏ để giảm sai số". Điều này bao gồm việc tích hợp các băng phổ bổ sung ("red-edge, SWIR") và các kỹ thuật tiền xử lý nâng cao.
    • Năm 3-5 (Mở rộng mô hình hóa): "Phát triển mô hình riêng cho rừng trồng" để giải quyết các đặc điểm độc đáo của chúng và mở rộng phạm vi áp dụng. Nghiên cứu sâu hơn về sự biến thiên của các hệ số CF và R "theo nhóm loài hoặc vùng sinh thái" thay vì giá trị mặc định.
    • Năm 5-7 (Tăng cường độ tin cậy và tích hợp): "Triển khai phân tích bất định nhằm nâng cao độ tin cậy của kết quả" theo chuẩn IPCC. Phát triển các khung tích hợp đa nguồn dữ liệu (Sentinel-2, LiDAR, UAV) để có ước tính các-bon toàn diện hơn.
    • Năm 7-10 (Ứng dụng quy mô lớn và chính sách): Kiểm chứng và nhân rộng các mô hình đã cải tiến cho "vùng sinh thái Nam Trung Bộ của Việt Nam" và các khu vực khác. Ứng dụng kết quả vào việc phát triển và thực thi chính sách giảm phát thải, quản lý rừng bền vững và thị trường tín chỉ các-bon ở cấp quốc gia và khu vực, theo dõi tác động dài hạn của các giải pháp đề xuất.

Kết luận

Luận án này đã tạo ra một bước tiến quan trọng trong lĩnh vực lâm sinh và viễn thám, đặc biệt là trong bối cảnh các nỗ lực toàn cầu về giảm phát thải và quản lý rừng bền vững.

  1. 5-6 đóng góp cụ thể (được đánh số):

    1. Phát triển mô hình hồi quy phi tuyến đa biến có trọng số: Luận án đã thành công trong việc xây dựng "mô hình hồi quy phi tuyến đa biến có trọng số" (ví dụ: TAGC = a + b × (NIR × EVI) + c × (NIR × EVI)²) để ước tính trữ lượng các-bon trên mặt đất từ ảnh Sentinel-2, vượt qua hạn chế của "cách tiếp cận tuyến tính truyền thống" và xử lý hiệu quả hiện tượng heteroscedasticity, đạt "MPSE = 21,63%" và "R² = 0,8764".
    2. Ứng dụng quy trình thẩm định chéo nghiêm ngặt: Việc áp dụng "thẩm định chéo Monte Carlo 100 lần lặp trên R-Studio" với các tiêu chí lựa chọn mô hình đa dạng (AIC, R², ASE, RMSE, MPSE) đã nâng cao đáng kể độ tin cậy và khả năng tổng quát hóa của mô hình, đảm bảo tính ứng dụng thực tiễn.
    3. Tối ưu hóa lựa chọn chỉ số viễn thám bằng PCA: Nghiên cứu đã sử dụng "Phân tích Thành phần Chính (PCA)" để xác định các chỉ số viễn thám (NDVI, SAVI, NIR) có "tương quan mạnh nhất" với TAGC và hình thành "các tổ hợp chỉ số hiệu quả hơn", cải thiện khả năng diễn giải thông tin sinh khối.
    4. Xây dựng bản đồ hiện trạng rừng và phân cấp trữ lượng các-bon đáng tin cậy: Luận án đã tạo ra "Bản đồ hiện trạng rừng năm 2023" của Đà Nẵng với "độ chính xác tổng thể 96%, hệ số Kappa 0,94" và đề xuất phân cấp trữ lượng các-bon rừng tự nhiên thành 3 cấp với sai số trung bình "9,45%", cung cấp công cụ quản lý thiết yếu.
    5. Đánh giá định lượng biến động các-bon và đề xuất chính sách: Nghiên cứu đã định lượng "phát thải ròng là –8.266,7 tC" tại Đà Nẵng giai đoạn 2016–2023 và đưa ra "các giải pháp hạn chế mất rừng và suy thoái rừng" cùng lộ trình thực hiện cụ thể, bao gồm quy hoạch phục hồi và tăng cường trữ lượng các-bon.
    6. Cung cấp giải pháp giảm chi phí cho kiểm kê các-bon: Mô hình cho phép ước tính trữ lượng các-bon "mà không cần đo đếm trực tiếp, giúp giảm thời gian và chi phí", làm cho việc kiểm kê và giám sát các-bon quy mô lớn trở nên khả thi hơn.

  2. Tiến bộ mô hình (Paradigm advancement) với bằng chứng: Luận án đại diện cho một tiến bộ đáng kể trong mô hình hậu thực chứng (post-positivist paradigm) trong nghiên cứu lâm sinh. Thay vì một sự thay đổi mô hình toàn diện, nó củng cố và tinh chỉnh các phương pháp định lượng hiện có. Bằng chứng rõ ràng nhất là việc nâng cao độ chính xác và độ tin cậy của các mô hình ước tính các-bon lên mức chưa từng có trong các nghiên cứu trước đây ở khu vực, với "R² đến 0,89" và "MPSE < 30%" từ các mô hình đã được thẩm định chéo Monte Carlo. Điều này cho thấy khả năng ngày càng tăng của khoa học để đo lường và dự đoán các hiện tượng sinh thái phức tạp với độ chính xác cao hơn, từ đó cung cấp cơ sở vững chắc hơn cho các quyết định chính sách và quản lý.

  3. 3+ luồng nghiên cứu mới được mở ra:

    1. Mô hình hóa các-bon rừng trồng và sinh khối dưới mặt đất: Luận án đã chỉ ra sự cần thiết của "phát triển mô hình riêng cho rừng trồng" và "bổ sung điều tra sinh khối dưới mặt đất", mở ra luồng nghiên cứu tập trung vào các thành phần các-bon hiện chưa được phân tích sâu.
    2. Phân tích bất định theo chuẩn IPCC và tích hợp đa nguồn dữ liệu: Nhu cầu về "phân tích bất định theo chuẩn IPCC" và "tích hợp đa nguồn dữ liệu" (bao gồm các kênh red-edge, SWIR, LiDAR) mở ra một luồng nghiên cứu về việc nâng cao độ chính xác và tính minh bạch của ước tính các-bon theo tiêu chuẩn quốc tế.
    3. Tối ưu hóa hệ số sinh khối và các-bon theo ngữ cảnh địa phương: Đề xuất sử dụng "hệ số CF và R theo nhóm loài hoặc vùng sinh thái" thay vì giá trị mặc định thúc đẩy luồng nghiên cứu về việc phát triển các phương trình allometric và hệ số các-bon được hiệu chỉnh địa phương, giảm bất định trong ước tính các-bon.

  4. Mức độ liên quan toàn cầu với so sánh quốc tế: Các phương pháp luận của luận án, đặc biệt là việc sử dụng dữ liệu Sentinel-2 có sẵn toàn cầu và các kỹ thuật thống kê nâng cao, có mức độ liên quan toàn cầu cao. Trong khi nhiều nghiên cứu quốc tế trước đây (như các nghiên cứu dựa trên mô hình tuyến tính của J. Cohen và P. S. A. G. D. S. Saatchi et al.) đối mặt với vấn đề bão hòa tín hiệu hoặc thẩm định ít nghiêm ngặt, luận án này đã đưa ra một giải pháp robust hơn. Việc cung cấp "mô hình hồi quy phi tuyến đa biến có trọng số" và quy trình "thẩm định chéo Monte Carlo 100 lần lặp" làm cho nghiên cứu này trở thành một ví dụ điển hình cho các quốc gia đang phát triển khác muốn tăng cường khả năng MRV của mình, đặc biệt là trong bối cảnh thực hiện các cam kết theo Thỏa thuận Paris.

  5. Kết quả có thể đo lường được của di sản: Di sản của luận án có thể được đo lường thông qua:

    • Áp dụng chính sách: Việc các mô hình và khuyến nghị được "Sở Nông nghiệp và Môi trường" và "Quỹ Bảo vệ và Phát triển rừng thành phố" sử dụng trong thực tế kiểm kê và giám sát rừng.
    • Đầu tư thực tế: Tổng kinh phí ước tính "~1.253 tỷ đồng" được phân bổ cho các hoạt động phục hồi và quản lý rừng dựa trên các phát hiện của luận án.
    • Tác động kinh tế: Số lượng tín chỉ các-bon được giao dịch thành công và lợi ích tài chính thu được cho địa phương nhờ vào các ước tính đáng tin cậy.
    • Các ấn phẩm học thuật tiếp theo: Số lượng trích dẫn các công trình liên quan đến luận án (ví dụ: "Nguyễn Hữu Tâm, Nguyễn Văn Lợi & Hoàng Huy Tuấn (2025)... Tạp chí Khoa học Đại học Huế") và các dự án nghiên cứu mới được kích hoạt từ các giới hạn đã được xác định.