Tổng quan nghiên cứu

Quá trình đô thị hóa nhanh chóng và biến đổi khí hậu toàn cầu đang tác động sâu sắc đến môi trường sinh thái của các đô thị cao nguyên. Tại thành phố Đà Lạt, số liệu quan trắc từ trạm khí tượng cho thấy nhiệt độ không khí trung bình năm đã tăng khoảng 0,6°C trong giai đoạn 2003–2012, với biên độ chênh lệch nhiệt độ ngày và đêm dao động mạnh từ 12°C đến 15°C. Tuy nhiên, việc đo đạc nhiệt độ bằng các trạm khí tượng mặt đất chỉ phản ánh điều kiện vi khí hậu cục bộ tại từng điểm riêng lẻ, không thể bao quát toàn bộ không gian lãnh thổ rộng 394,64 km² gồm 12 phường và 4 xã.

Vấn đề nghiên cứu trọng tâm là sự gia tăng nhiệt độ bề mặt đất do suy giảm thảm thực vật tự nhiên và sự mở rộng của các bề mặt nhân tạo như bê tông, nhà kính nông nghiệp. Mục tiêu cụ thể của đề tài là hệ thống hóa cơ sở khoa học, ứng dụng công nghệ viễn thám đa phổ, hồng ngoại nhiệt và hệ thống thông tin địa lý (GIS) để chiết xuất thông tin lớp phủ, tính toán nhiệt độ bề mặt qua 4 mốc thời gian (1996, 2003, 2010, 2014). Phạm vi không gian nghiên cứu giới hạn trong ranh giới hành chính thành phố Đà Lạt (tọa độ từ 11°48’36" đến 12°01'07" vĩ độ Bắc và 108°19'23" đến 108°36'27" kinh độ Đông) ở độ cao trung bình 1.520 m. Nghiên cứu mang ý nghĩa thực tiễn lớn khi định lượng được mức tăng nhiệt độ bề mặt trung bình 6,8°C sau 18 năm, cung cấp cơ sở khoa học để chính quyền địa phương quy hoạch mảng xanh và quản lý tài nguyên bền vững.

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 dựa trên định luật bức xạ nhiệt Stefan-Boltzmann về mối liên hệ giữa năng lượng phát xạ và nhiệt độ nhiệt động lực học, kết hợp với định luật Planck và lý thuyết cân bằng năng lượng bề mặt đất - khí quyển. Mô hình nghiên cứu tích hợp công nghệ Viễn thám và GIS để phân tích mối tương quan giữa biến động lớp phủ mặt đất và sự thay đổi nhiệt độ bề mặt.

Khung lý thuyết của đề tài xoay quanh 5 khái niệm cốt lõi:

  • Nhiệt độ bề mặt động lực (Land Surface Temperature - Ts): Tham số vật lý đo lường trạng thái cân bằng bức xạ thực tế của lớp vật chất trên cùng của bề mặt đất.
  • Nhiệt độ sáng (Brightness Temperature - Tb): Giá trị nhiệt độ bức xạ tương đương của vật đen tuyệt đối thu nhận bởi cảm biến nhiệt vệ tinh.
  • Độ phát xạ bề mặt (Surface Emissivity - e): Tỷ số giữa năng lượng bức xạ phát ra từ một bề mặt thực tế so với vật đen ở cùng nhiệt độ, được ước tính qua chỉ số thực vật.
  • Chỉ số thực vật chuẩn hóa (NDVI): Chỉ số không gian tính toán từ kênh đỏ và cận hồng ngoại nhằm định lượng sinh khối và mật độ xanh của thảm thực vật.
  • Đảo nhiệt đô thị bề mặt (Surface Urban Heat Island - SUHI): Hiện tượng nền nhiệt độ bề mặt tại các khu vực tập trung công trình xây dựng cao hơn rõ rệt so với vùng phủ xanh và mặt nước lân cận.

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

Cỡ mẫu nghiên cứu bao gồm 4 cảnh ảnh vệ tinh Landsat đa thời gian: Landsat 5 TM (ngày 01/3/1996), Landsat 7 ETM+ (ngày 29/3/2003 và 08/3/2010), Landsat 8 OLI/TIRS (ngày 16/3/2014). Toàn bộ dữ liệu ảnh có độ phân giải không gian 30 m đối với các kênh quang học và từ 60 m đến 120 m đối với kênh nhiệt, kết hợp cùng dữ liệu ranh giới GIS của 16 đơn vị hành chính và bản đồ hiện trạng sử dụng đất.

Phương pháp chọn mẫu là chọn mẫu có chủ đích vào cùng thời điểm tháng 3 (mùa khô) nhằm đảm bảo độ che phủ mây dưới 10%, loại bỏ sự sai khác do góc chiếu mặt trời và độ ẩm theo mùa. Đối với phân loại ảnh, mẫu huấn luyện được chọn lọc đại diện cho 4 lớp phủ chính: đất đô thị, đất trống, thực phủ và mặt nước.

Lý do lựa chọn phương pháp phân tích:

  • Thuật toán phân loại có giám định theo xác suất cực đại (Maximum Likelihood Classifier - MLC) được chọn vì có cơ sở toán học vững chắc, đánh giá xác suất phân bố chuẩn của các pixel để gán nhãn lớp phủ chính xác nhất. Kết quả phân loại được lọc mịn bằng thuật toán Majority và kiểm định bằng ma trận sai số.
  • Phương pháp ước tính nhiệt độ bề mặt sử dụng thuật toán NDVI của Valor và Caselles để tính độ phát xạ bề mặt (e), hiệu chỉnh chính xác năng lượng bức xạ phổ cho khu vực địa hình cao nguyên phức tạp.
  • Phương pháp hồi quy tuyến tính đa biến với kỹ thuật Bayesian Model Averaging (BMA) và hồi quy không gian (Spatial Regression) được lựa chọn nhằm giải quyết hiện tượng tự tương quan không gian của phần dư, cho phép mô hình hóa trực quan mức độ tác động của từng loại lớp phủ đến sự tăng nhiệt độ bề mặt tại từng tọa độ địa lý.

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

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

  • Thứ nhất, nhiệt độ bề mặt trung bình của thành phố Đà Lạt trong giai đoạn 1996–2014 đã gia tăng 6,8°C. Xu hướng tăng nhiệt diễn ra liên tục qua các mốc thời gian, với các khu vực nhiệt độ cao mở rộng rõ rệt từ vùng lõi trung tâm ra các khu vực vùng ven.
  • Thứ hai, cơ cấu lớp phủ mặt đất có sự chuyển dịch sâu sắc: diện tích đất đô thị và đất trống tăng mạnh, trong khi diện tích thực phủ (rừng thông và cây bụi) sụt giảm đáng kể. Hàng nghìn hecta đất tự nhiên đã bị chuyển đổi sang đất xây dựng và đất sản xuất nông nghiệp thâm canh.
  • Thứ ba, phân tích hồi quy tuyến tính đa biến xác nhận mối tương quan chặt chẽ giữa lớp phủ và nhiệt độ bề mặt: diện tích đất đô thị và đất trống có tương quan dương với nhiệt độ bề mặt, trong khi diện tích thực phủ và mặt nước có tương quan âm. Đất trống là yếu tố có tác động làm gia tăng nhiệt độ bề mặt mạnh nhất, tiếp theo là đất đô thị; ngược lại, thực phủ và mặt nước đóng vai trò hạ nhiệt hiệu quả.
  • Thứ tư, phân tích thống kê không gian bằng tâm trung bình trọng số và ellipse độ lệch chuẩn chỉ ra hướng dịch chuyển của các điểm nóng nhiệt độ bề mặt theo trục Đông Bắc – Tây Nam, trùng khớp hoàn toàn với hành lang phát triển hạ tầng và vùng canh tác nông nghiệp công nghệ cao của thành phố.

Thảo luận kết quả

Nguyên nhân chính của sự gia tăng nhiệt độ bề mặt tại Đà Lạt là do bề mặt tự nhiên có khả năng bốc thoát hơi nước và điều hòa nhiệt cao (rừng thông, hồ nước) bị thay thế bởi các bề mặt nhân tạo có nhiệt dung riêng thấp và khả năng hấp thụ nhiệt lớn (bê tông, mái tôn, màng nilon nhà kính nông nghiệp). Đất trống bị phơi lộ trực tiếp dưới bức xạ mặt trời mùa khô cũng hấp thụ một lượng nhiệt năng khổng lồ, tạo thành các trung tâm nhiệt cục bộ.

Về mặt trực quan, dữ liệu nghiên cứu được trình bày chi tiết thông qua hệ thống bản đồ phân cấp nhiệt độ bề mặt theo 5 mức nhiệt độ, bản đồ ma trận chuyển dịch lớp phủ giai đoạn 1996–2003, 2003–2010, 2010–2014, cùng các biểu đồ phân tích phần dư chuẩn hóa và biểu đồ hộp (box-plot) thể hiện phân bố nhiệt độ theo từng loại vật liệu. Bảng ma trận biến động diện tích 2 chiều đã lượng hóa chính xác tỷ lệ chuyển đổi từ thảm thực vật sang đất phi nông nghiệp qua từng thời kỳ.

Kết quả này hoàn toàn tương thích với các nghiên cứu tương tự tại Thành phố Hồ Chí Minh (nhiệt độ bề mặt tăng bình quân 0,26°C/năm trong giai đoạn 1989–2006) hay vùng châu thổ sông Châu Giang. Tuy nhiên, tại Đà Lạt, do yếu tố địa hình đồi núi dốc và mật độ nhà kính nông nghiệp dày đặc, tác động làm tăng nhiệt của đất trống và bề mặt nông nghiệp nhân tạo diễn ra mạnh mẽ hơn cả khu vực đô thị thuần túy, làm suy giảm vi khí hậu mát mẻ truyền thống của cao nguyên Lâm Viên.

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

  • Quy hoạch và kiểm soát mật độ xây dựng nhà kính: Ban hành quy chuẩn quản lý diện tích nhà kính tại 12 phường nội ô; bắt buộc áp dụng tỷ lệ mảng xanh tối thiểu 40% trong các dự án chỉnh trang đô thị mới trước năm 2030. Đơn vị chủ trì thực hiện là Sở Xây dựng Lâm Đồng phối hợp cùng UBND TP. Đà Lạt.
  • Tăng cường bảo vệ và phục hồi thảm thực vật rừng: Triển khai trồng mới và phục hồi 1.500 hecta rừng phòng hộ, rừng thông cảnh quan đầu nguồn quanh các hồ Dankia, Tuyền Lâm và các sườn đồi trước năm 2028 nhằm nâng chỉ số NDVI bình quân của thành phố lên mức trên 0,55. Đơn vị thực hiện là Chi cục Kiểm lâm tỉnh Lâm Đồng và các ban quản lý rừng.
  • Bảo tồn và mở rộng không gian mặt nước: Nạo vét, bảo vệ nguồn nước và mở rộng diện tích hành lang cảnh quan mặt nước xung quanh Hồ Xuân Hương, Hồ Than Thở và hệ thống hồ chứa thêm 15% vào năm 2027 để tăng khả năng bốc thoát hơi nước làm mát tự nhiên. Đơn vị thực hiện là Sở Tài nguyên và Môi trường tỉnh Lâm Đồng.
  • Thiết lập hệ thống giám sát viễn thám định kỳ: Xây dựng cơ sở dữ liệu viễn thám tự động cập nhật biến động lớp phủ và nhiệt độ bề mặt định kỳ 6 tháng một lần bằng ảnh vệ tinh Landsat và Sentinel; tích hợp cảnh báo sớm các điểm nóng nhiệt độ bề mặt vượt ngưỡng 32°C vào Trung tâm Giám sát Điều hành Thông minh (IOC) tỉnh Lâm Đồng giai đoạn 2026–2030. Đơn vị thực hiện là Sở Thông tin và Truyền thông tỉnh Lâm Đồng.

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

  • Cơ quan quản lý nhà nước và nhà quy hoạch đô thị: Lãnh đạo Sở Tài nguyên và Môi trường, Sở Xây dựng tỉnh Lâm Đồng và UBND TP. Đà Lạt có thể sử dụng dữ liệu chuỗi 18 năm và mô hình hồi quy không gian làm cơ sở xây dựng đồ án quy hoạch tổng thể thành phố Đà Lạt tầm nhìn 2045, nhằm cân đối tỷ lệ cây xanh và kiểm soát hiệu ứng đảo nhiệt.
  • Chuyên gia, nhà nghiên cứu và học viên ngành Viễn thám - GIS: Các nhà khoa học trong lĩnh vực bản đồ, trắc địa và quản lý tài nguyên có thể tham khảo quy trình tiền xử lý, thuật toán tính độ phát xạ bề mặt theo NDVI trên địa hình cao nguyên và phương pháp mô hình hóa hồi quy không gian kết hợp BMA.
  • Các tổ chức bảo tồn môi trường và biến đổi khí hậu: Các tổ chức phi chính phủ và viện nghiên cứu sinh thái có nguồn tư liệu định lượng xác thực về mức tăng nhiệt 6,8°C để làm bằng chứng khoa học cho các dự án bảo vệ đa dạng sinh học và vận động chính sách phục hồi rừng thông.
  • Doanh nghiệp đầu tư du lịch, nghỉ dưỡng và nông nghiệp: Các nhà đầu tư tại 16 phường, xã của Đà Lạt có thể khai thác bản đồ phân bố nhiệt độ để đánh giá vi khí hậu khu vực, từ đó thiết kế các công trình kiến trúc xanh giảm bức xạ nhiệt và định hướng mô hình canh tác nông nghiệp thân thiện với cảnh quan.

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

  • Tại sao nghiên cứu lại chọn ảnh vệ tinh vào tháng 3 thay vì các tháng khác trong năm? Tháng 3 là thời kỳ đỉnh điểm mùa khô tại Đà Lạt với tỷ lệ mây che phủ dưới 10%, giúp giảm thiểu tối đa hiện tượng tán xạ và hấp thụ của khí quyển. Việc cố định thời điểm tháng 3 trong chuỗi 4 mốc 1996, 2003, 2010, 2014 đảm bảo sự đồng nhất về góc chiếu mặt trời, độ ẩm bề mặt và độ chính xác khi so sánh chuỗi biến động đa thời gian.

  • Độ phát xạ bề mặt được tính toán như thế nào để phù hợp với địa hình đồi núi cao nguyên? Nghiên cứu sử dụng phương pháp chỉ số thực vật NDVI của Valor và Caselles, kết hợp công thức thực nghiệm của Van de Griend và Owe để xác định ngưỡng phát xạ riêng cho đất trống và thảm thực vật. Phương pháp này tính toán chính xác đặc trưng độ nhám của bề mặt hỗn hợp ở độ cao 1.520 m thay vì gán giá trị độ phát xạ cố định bằng 1.

  • Thuật toán phân loại ảnh có giám định nào được áp dụng và mức độ tin cậy đạt bao nhiêu? Thuật toán phân loại gần đúng nhất (Maximum Likelihood Classifier) được sử dụng kết hợp bộ lọc Majority để loại bỏ nhiễu điểm ảnh. Kết quả kiểm tra ma trận sai số cho thấy độ chính xác toàn cục của các bản đồ lớp phủ qua 4 thời điểm đều đạt trên 85% và hệ số Kappa vượt mức 0,82, khẳng định độ tin cậy khoa học rất cao.

  • Nhiệt độ bề mặt từ ảnh viễn thám khác gì so với nhiệt độ không khí tại trạm quan trắc? Nhiệt độ không khí đo tại lều khí tượng ở độ cao 2 m chỉ phản ánh điều kiện vi khí hậu cục bộ tại điểm đặt trạm. Trong khi đó, nhiệt độ bề mặt (Ts) thu nhận từ kênh hồng ngoại nhiệt vệ tinh phản ánh trực tiếp năng lượng bức xạ phát ra từ vật liệu bề mặt trên toàn bộ diện tích 394,64 km², cho phép giám sát diện rộng và liên tục.

  • Mô hình hồi quy không gian mang lại ưu điểm gì so với mô hình hồi quy truyền thống? Hồi quy truyền thống (OLS) chỉ đưa ra một hệ số toàn cục duy nhất và dễ bị sai số khi phần dư có hiện tượng tự tương quan không gian. Mô hình hồi quy không gian liên kết trực tiếp dữ liệu thống kê với tọa độ địa lý, phản ánh rõ nét mức độ tác động mang tính cục bộ của từng loại lớp phủ lên nhiệt độ bề mặt tại từng đơn vị hành chính.

Kết luận

  • Ứng dụng thành công công nghệ tích hợp Viễn thám và GIS để giám sát biến động lớp phủ và nhiệt độ bề mặt thành phố Đà Lạt qua chuỗi thời gian 18 năm (1996–2014).
  • Ghi nhận nhiệt độ bề mặt trung bình toàn thành phố tăng 6,8°C trong giai đoạn nghiên cứu, tập trung cao nhất tại các khu vực đô thị hóa và vùng nông nghiệp nhà kính.
  • Chứng minh diện tích đất trống và đất đô thị có tương quan dương mạnh nhất với nhiệt độ bề mặt, trong khi thực phủ và mặt nước đóng vai trò chủ đạo trong việc giảm nhiệt cảnh quan.
  • Xác định trục biến đổi không gian nhiệt độ bề mặt theo hướng Đông Bắc – Tây Nam bằng các công cụ phân tích thống kê không gian hiện đại.
  • Lộ trình giai đoạn 2026–2030 cần tiếp tục mở rộng chuỗi dữ liệu với vệ tinh Landsat 9 và Sentinel-2, tích hợp trí tuệ nhân tạo để dự báo kịch bản đảo nhiệt đô thị phục vụ phát triển Đà Lạt thành đô thị xanh bền vững.