Tổng quan nghiên cứu

Nhiệt độ bề mặt đất đóng vai trò cốt lõi trong việc điều hòa chu trình năng lượng giữa thảm phủ bề mặt và khí quyển, đồng thời là chỉ số sinh thái quan trọng phản ánh trực tiếp mức độ thoái hóa môi trường, hạn hán và biến đổi khí hậu. Tại huyện Cao Phong, tỉnh Hòa Bình, quá trình chuyển dịch cơ cấu kinh tế mạnh mẽ trong giai đoạn 2002 đến 2016 đã làm thay đổi sâu sắc diện mạo cảnh quan; tỷ trọng ngành nông nghiệp giảm từ 72% xuống còn 46%, trong khi tỷ trọng công nghiệp và xây dựng tăng từ 18% lên 28%, thúc đẩy nhanh chóng tốc độ đô thị hóa và bê tông hóa. Thực trạng suy giảm thảm thực vật tự nhiên và gia tăng các bề mặt nhân tạo đã làm xuất hiện các điểm tích tụ nhiệt, tạo ra hiện tượng đảo nhiệt cục bộ, ảnh hưởng trực tiếp đến môi trường sống của hơn 43.000 cư dân trên tổng diện tích tự nhiên 254 km².

Mục tiêu trọng tâm của nghiên cứu là xác định quy luật phân bố không gian và biến động nhiệt độ bề mặt đất của huyện Cao Phong trong năm 2016 thông qua việc khai thác dữ liệu ảnh vệ tinh Landsat 8 OLI/TIRS. Nghiên cứu tập trung so sánh hiệu quả định lượng giữa hai phương pháp trích lọc nhiệt độ tiên tiến, làm rõ mối quan hệ tương quan giữa nhiệt độ bề mặt với các loại hình lớp phủ, và khoanh vùng các khu vực đảo nhiệt. Kết quả nghiên cứu không chỉ thiết lập một cơ sở dữ liệu viễn thám nhiệt có độ tin cậy cao mà còn cung cấp luận cứ khoa học thực tiễn giúp địa phương tối ưu hóa công tác quy hoạch không gian đô thị, quản lý hiệu quả hơn 1.700 ha cây ăn quả có múi chủ lực và thích ứng bền vững với các hiện tượng thời tiết cực đoan.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng các định luật vật lý bức xạ nhiệt kinh điển của vật chất, tiêu biểu là định luật Stefan-Boltzmann xác định tổng năng lượng phát xạ toàn phần tỷ lệ thuận với lũy thừa bậc bốn của nhiệt độ tuyệt đối ($M = \varepsilon \sigma T^4$), định luật dịch chuyển cực trị bức xạ Wien ($\lambda_{max} = A/T$), và định luật phân bố phổ Planck. Trong môi trường tự nhiên, các đối tượng địa lý hoạt động như các vật phát xạ chọn lọc với hệ số phát xạ $\varepsilon$ dao động từ 0 đến 1, khiến nhiệt độ phát xạ bức xạ ghi nhận qua bộ cảm biến viễn thám luôn có sự khác biệt nhất định so với nhiệt độ động lực học bên trong của vật chất.

Hệ thống khái niệm học thuật then chốt trong đề tài bao gồm:

  • Nhiệt độ độ sáng: Nhiệt độ tương đương của vật đen tuyệt đối phát ra cùng một lượng bức xạ tại bước sóng quan trắc, được tính trực tiếp từ giá trị độ chói phổ trên vệ tinh.
  • Độ phát xạ bề mặt: Tỷ số giữa năng lượng bức xạ phát ra từ bề mặt thực tế so với vật đen tuyệt đối ở cùng nhiệt độ và bước sóng.
  • Nhiệt độ bề mặt đất: Giá trị nhiệt độ trung bình bức xạ của bề mặt hỗn hợp bao gồm đất trống, thân cây, tán lá và công trình xây dựng trong một đơn vị không gian điểm ảnh.
  • Chỉ số khác biệt thực vật chuẩn hóa: Tỷ số phản xạ phổ giữa kênh cận hồng ngoại và kênh đỏ nhằm định lượng mật độ sinh khối và sức khỏe của thảm thực vật.
  • Hiện tượng đảo nhiệt: Khu vực cục bộ có nhiệt độ bề mặt và nhiệt độ không khí cao hơn rõ rệt so với các vùng cảnh quan sinh thái tự nhiên bao quanh.

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

Nguồn dữ liệu chính của nghiên cứu là các cảnh ảnh vệ tinh Landsat 8 OLI/TIRS thuộc hệ thống tọa độ Path 127, Row 46, được thu nhận tại hai thời điểm tương phản theo mùa: ngày 10/02/2016 (mùa đông) và ngày 01/06/2016 (mùa hè). Kênh hồng ngoại nhiệt Band 10 (bước sóng 10.60 đến 11.19 $\mu m$) với độ phân giải ban đầu 100m được tái lập mẫu về 30m, kết hợp đồng bộ cùng các kênh ảnh quang học, mô hình số độ cao DEM (độ cao trung bình 399m), bản đồ thổ nhưỡng chi tiết trên diện tích 254 km² và chuỗi số liệu quan trắc khí tượng bề mặt.

Cỡ mẫu nghiên cứu bao quát toàn bộ mạng lưới không gian với hơn 280.000 điểm ảnh phân bố trên toàn lãnh thổ huyện Cao Phong. Phương pháp chọn mẫu phân tầng theo các đơn vị lớp phủ và địa hình được áp dụng để đối soát dữ liệu với các điểm đo nhiệt độ thực tế tại trạm quan trắc địa phương.

Quy trình phân tích thực nghiệm so sánh hai thuật toán trích lọc nhiệt độ bề mặt chuyên sâu:

  1. Phương pháp xác định nhiệt độ bề mặt có tính độ phát xạ từ chỉ số thực vật NDVI: Tiến hành hiệu chỉnh bức xạ từ giá trị số DN sang độ chói phổ $L$, áp dụng hiệu chỉnh khí quyển DOS/ISAC, bóc tách tỷ lệ thực vật $P_v$ theo mô hình của Valor để tính toán độ phát xạ $\varepsilon$ cho từng pixel trước khi đưa vào phương trình hàm sóng nghịch đảo Planck.
  2. Phương pháp chuẩn hóa độ phát xạ NOR: Giả định tồn tại một độ phát xạ cực đại cố định để suy biến nhiệt độ bề mặt lớn nhất qua các kênh nhiệt.

Lý do lựa chọn phương pháp hiệu chỉnh phát xạ dựa trên NDVI là khả năng giải quyết triệt để tính chất không đồng nhất của các bề mặt pha trộn phức tạp tại vùng đồi núi, giúp triệt tiêu sai số do chênh lệch thảm phủ và nâng cao độ chính xác định lượng so với phương pháp chuẩn hóa thông thường.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình trích lọc và đối sánh dữ liệu thực nghiệm đã mang lại những kết quả quan trọng về đặc trưng nhiệt độ bề mặt tại huyện Cao Phong:

  • Phân hóa nhiệt độ theo mùa cực kỳ sâu sắc: Vào mùa hè (ảnh ngày 01/06/2016), nhiệt độ bề mặt dao động mạnh trong khoảng từ 28.5°C đến 42.8°C, với các mức nhiệt độ cực đại tập trung tại vùng đất trống và khu dân cư. Ngược lại, vào mùa đông (ảnh ngày 10/02/2016), nhiệt độ bề mặt toàn huyện hạ thấp rõ rệt, chỉ duy trì trong khoảng từ 14.2°C đến 22.6°C, phản ánh biên độ nhiệt mùa vượt mức 18.0°C.
  • Độ chính xác vượt trội của phương pháp hiệu chỉnh phát xạ từ NDVI: Kết quả đối chiếu với dữ liệu khí tượng mặt đất cho thấy phương pháp xác định nhiệt độ dựa trên độ phát xạ NDVI đạt độ lệch chuẩn rất thấp, sai số tuyệt đối so với nhiệt độ không khí chỉ dao động trong ngưỡng dưới 1.8°C. Trong khi đó, phương pháp chuẩn hóa NOR cho độ lệch nhiệt độ cao hơn, dao động từ 2.5°C đến 3.6°C, chứng minh phương pháp NDVI phản ánh sát thực tế cấu trúc bề mặt hơn.
  • Mối quan hệ tương quan nghịch chặt chẽ giữa NDVI và nhiệt độ bề mặt: Tại các khu vực có thảm thực vật dày đặc với chỉ số NDVI > 0.5 (như rừng tự nhiên xã Yên Thượng, Thung Nai), nhiệt độ bề mặt trung bình thấp hơn từ 4.5°C đến 7.2°C so với các khu vực đất trống có chỉ số NDVI < 0.2.
  • Hình thành các tâm đảo nhiệt rõ nét: Các điểm dị thường nhiệt độ cao tập trung chủ yếu dọc theo hành lang Quốc lộ 6, trung tâm thị trấn Cao Phong và các vùng mặt bằng xây dựng cụm công nghiệp, nơi tỷ lệ che phủ công trình và bề mặt không thấm nước chiếm hơn 65% diện tích bề mặt ô lưới.

Thảo luận kết quả

Sự hình thành các khu vực nhiệt độ cao và đảo nhiệt cục bộ bắt nguồn từ sự kết hợp của nhiều yếu tố tự nhiên và nhân sinh. Về mặt vật lý, vật liệu xây dựng nhân tạo như bê tông, nhựa đường và mái tôn có nhiệt dung riêng thấp và khả năng hấp thụ bức xạ sóng ngắn rất mạnh, đồng thời độ phát xạ bề mặt thấp ($\varepsilon \approx 0.92 - 0.94$) làm tăng nhanh nhiệt bức xạ ban ngày. Ngược lại, hệ thống hơn 108 ha mặt nước hồ sông Đà tại xã Bình Thanh, Thung Nai cùng diện tích rừng trên đất Feralit mùn núi cao đóng vai trò như các bể hấp thụ nhiệt tự nhiên nhờ quá trình bốc thoát hơi nước liên tục, giúp hạ nhiệt độ môi trường xung quanh từ 3.0°C đến 5.0°C.

Về mặt trình bày và trực quan hóa dữ liệu, các kết quả tính toán có thể được biểu diễn một cách trực quan thông qua bản đồ phân vùng đẳng nhiệt đa phổ kết hợp lớp phủ, đồ thị tương quan hồi quy tuyến tính phân tán giữa NDVI và nhiệt độ bề mặt, cùng bảng ma trận đối chiếu nhiệt độ trung bình trên 4 nhóm đất thổ nhưỡng chính (nhóm đất Fv đá vôi 4.495 ha, nhóm đất Fq đá cát 7.953 ha).

Khi so sánh với các công trình nghiên cứu viễn thám nhiệt tại các đô thị khác như Đà Lạt hay Thành phố Hồ Chí Minh, quy luật giảm nhiệt nhờ màng sinh thái thực vật tại Cao Phong hoàn toàn tương đồng. Điều này khẳng định rằng dù ở vùng đồng bằng hay địa hình miền núi Tây Bắc, quá trình chuyển đổi lớp phủ từ sinh thái sang hạ tầng xây dựng luôn là nhân tố hàng đầu chi phối vi khí hậu khu vực.

Đề xuất và khuyến nghị

Dựa trên kết quả phân tích không gian nhiệt độ bề mặt và hiện trạng sử dụng đất, nghiên cứu đưa ra 4 nhóm giải pháp chiến lược nhằm điều hòa vi khí hậu huyện Cao Phong:

  • Quy hoạch không gian xanh và thiết lập hành lang cách ly nhiệt: Ủy ban nhân dân huyện Cao Phong cùng Phòng Kinh tế và Hạ tầng cần ban hành quy định bắt buộc duy trì tỷ lệ mảng xanh tối thiểu 20% trong các dự án mở rộng đô thị và phát triển cụm công nghiệp mới từ nay đến năm 2030, ưu tiên trồng cây bóng mát có tán rộng dọc trục giao thông Quốc lộ 6 để giảm nhiệt độ bề mặt hạ tầng xuống ít nhất 2.5°C.
  • Bảo vệ và phát triển hệ sinh thái nông lâm nghiệp chủ lực: Phòng Nông nghiệp và Phát triển Nông thôn huyện cần đẩy mạnh kế hoạch duy trì ổn định diện tích 1.700 ha cây ăn quả có múi theo tiêu chuẩn VietGAP, kết hợp bảo tồn nghiêm ngặt hơn 108 ha diện tích mặt nước tại vùng hồ Thung Nai và Bình Thanh, nhằm đảm bảo nguồn ẩm điều hòa nhiệt độ cho các mùa khô hạn giai đoạn 2025-2030.
  • Ứng dụng vật liệu xây dựng phản xạ nhiệt và hạ tầng tiêu chuẩn mát: Sở Xây dựng tỉnh Hòa Bình phối hợp cùng các đơn vị chủ đầu tư triển khai sử dụng vật liệu lát vỉa hè thấm nước, sơn phản xạ nhiệt albedo cao ($\ge 0.40$) cho các công trình công cộng và nhà xưởng, hướng tới mục tiêu giảm tích tụ nhiệt tại các tâm đảo nhiệt thị trấn xuống khoảng 1.5°C đến 3.0°C trước năm 2028.
  • Thiết lập hệ thống giám sát viễn thám định kỳ cảnh báo khô hạn: Chi cục Kiểm lâm và Trung tâm Ứng dụng Tiến bộ Khoa học Công nghệ tỉnh Hòa Bình cần định kỳ cập nhật ảnh viễn thám nhiệt 2 lần mỗi năm (vào tháng 2 và tháng 6) để tính toán chỉ số khô hạn TVDI, kịp thời phát hiện các điểm nguy cơ cháy rừng và hỗ trợ điều tiết nguồn nước tưới cho nông nghiệp.

Đối tượng nên tham khảo luận văn

Nội dung và phương pháp tiếp cận của luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng cụ thể:

  • Cơ quan quản lý nhà nước và nhà quy hoạch đô thị nông thôn: Sở Tài nguyên và Môi trường tỉnh Hòa Bình, Ủy ban nhân dân huyện Cao Phong có thể sử dụng bản đồ nhiệt độ bề mặt và các phân vùng đảo nhiệt làm luận cứ khoa học để phân bổ chỉ tiêu sử dụng đất, kiểm soát mật độ xây dựng và quy hoạch hạ tầng sinh thái bền vững.
  • Cán bộ nghiên cứu và chuyên gia nông nghiệp: Các kỹ sư nông nghiệp, cán bộ khuyến nông sử dụng dữ liệu quan hệ nhiệt - thực vật để đánh giá mức độ căng thẳng nhiệt của cây trồng, tối ưu hóa hệ thống thủy lợi và phòng ngừa sâu bệnh trên các vùng chuyên canh cam và mía tím.
  • Giảng viên, học viên cao học và sinh viên chuyên ngành Địa lý, Bản đồ, Viễn thám và GIS: Tiếp cận một khung quy trình xử lý chuẩn hóa ảnh Landsat 8 TIRS, từ bước hiệu chỉnh khí quyển, tính toán phát xạ bề mặt qua NDVI đến trích lọc nhiệt độ thực tế trên địa hình đồi núi phức tạp.
  • Các nhà đầu tư và doanh nghiệp phát triển du lịch sinh thái: Đơn vị quản lý khu du lịch lòng hồ sông Đà, bản Mường Giang Mỗ sử dụng bản đồ vi khí hậu để lựa chọn vị trí xây dựng các khu nghỉ dưỡng sinh thái có điều kiện nhiệt độ mát mẻ tự nhiên.

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

  • Tại sao cần tính toán độ phát xạ bề mặt thay vì sử dụng trực tiếp nhiệt độ độ sáng từ kênh nhiệt vệ tinh? Nhiệt độ độ sáng chỉ phản ánh năng lượng bức xạ của vật đen lý thuyết đo được tại đỉnh khí quyển. Trong thực tế, các lớp phủ tự nhiên có độ phát xạ biến thiên từ 0.77 đến 0.99. Nếu không hiệu chỉnh độ phát xạ, sai số nhiệt độ bề mặt có thể chênh lệch từ 2.0°C đến hơn 5.0°C so với nhiệt độ thực tế.

  • Phương pháp ước tính độ phát xạ từ chỉ số NDVI có ưu điểm vượt trội nào so với thuật toán chuẩn hóa NOR? Phương pháp NDVI bóc tách chính xác tỷ lệ pha trộn giữa đất trống ($\varepsilon = 0.95$) và thực vật thuần nhất ($\varepsilon = 0.98$) trong từng điểm ảnh 30m. Điều này giúp phản ánh chân thực cấu trúc bề mặt phân cắt mạnh tại Cao Phong, giảm sai số kiểm chứng xuống dưới 1.8°C so với mức sai số trên 2.5°C của thuật toán NOR.

  • Yếu tố nào đóng vai trò chính tạo nên hiện tượng đảo nhiệt cục bộ tại huyện Cao Phong? Nguyên nhân chủ yếu là do tốc độ mở rộng không gian xây dựng dọc các trục giao thông chính và cụm công nghiệp, làm thay đổi cơ cấu kinh tế theo hướng tăng tỷ trọng công nghiệp lên 28%. Quá trình này đã thay thế thảm thực vật bằng vật liệu hấp thụ nhiệt mạnh như bê tông và kim loại.

  • Dữ liệu viễn thám nhiệt đóng góp như thế nào cho vùng chuyên canh cam Cao Phong? Dữ liệu nhiệt độ bề mặt LST là biến số đầu vào then chốt để tính toán chỉ số khô hạn TVDI, giúp phát hiện sớm hiện tượng thiếu nước của tán lá trên hơn 1.700 ha cam. Qua đó, chính quyền và người dân có thể chủ động phương án tưới tiêu trong các đợt nhiệt độ vượt ngưỡng 40.0°C.

  • Độ phân giải không gian của ảnh Landsat 8 TIRS có đáp ứng tốt yêu cầu nghiên cứu cấp huyện? Kênh nhiệt Band 10 của Landsat 8 có độ phân giải gốc 100m và được nội suy xử lý về 30m cùng dải phủ 175 km. Quy mô này hoàn toàn phù hợp để phân tích bức tranh nhiệt độ chi tiết trên 254 km² lãnh thổ cấp huyện mà không gặp hiện tượng quá khái quát hóa như ảnh MODIS 1000m.

Kết luận

  • Đã xây dựng thành công quy trình công nghệ chuẩn hóa bóc tách nhiệt độ bề mặt đất LST từ kênh nhiệt Band 10 ảnh vệ tinh Landsat 8 OLI/TIRS tại huyện Cao Phong năm 2016.
  • Chứng minh phương pháp xác định nhiệt độ có hiệu chỉnh độ phát xạ từ chỉ số thực vật NDVI đạt độ tin cậy vượt trội với sai số thực nghiệm thấp dưới 1.8°C so với phương pháp chuẩn hóa NOR.
  • Khẳng định mối quan hệ tương quan nghịch sâu sắc giữa mật độ thảm phủ thực vật và nhiệt độ bề mặt, với mức chênh lệch nhiệt độ làm mát của rừng tự nhiên lên tới 7.2°C so với bề mặt đất trống.
  • Khoanh vùng chính xác các điểm nóng đảo nhiệt cục bộ tại khu vực trung tâm thị trấn Cao Phong và các công trình hạ tầng công nghiệp dọc tuyến Quốc lộ 6.
  • Đề xuất hệ thống 4 nhóm giải pháp đồng bộ kết hợp quy hoạch đô thị xanh, bảo tồn 1.700 ha vùng cây ăn quả và thiết lập mạng lưới quan trắc viễn thám định kỳ phục vụ mục tiêu phát triển bền vững.

Đóng góp khoa học cốt lõi của luận văn là hoàn thiện phương pháp luận viễn thám nhiệt cho vùng bán sơn địa và cung cấp nguồn tư liệu địa không gian chính xác cho công tác quản lý tài nguyên địa phương. Hướng nghiên cứu tiếp theo cần mở rộng chuỗi dữ liệu ảnh nhiệt đa thời gian liên tục từ năm 2016 đến nay nhằm đánh giá xu thế biến động vi khí hậu lâu dài. Các cấp quản lý quy hoạch và cơ quan chuyên môn tại tỉnh Hòa Bình nên khai thác ngay bộ bản đồ phân vùng nhiệt độ này vào kế hoạch hành động ứng phó biến đổi khí hậu và tối ưu hóa cảnh quan sinh thái trên toàn địa bàn.