NGHIÊN CỨU ĐÁNH GIÁ BIẾN ĐỘNG VÀ DỰ TÍNH HẠN KHÍ TƯỢNG THEO CHỈ SỐ ẨM DƯỚI TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU TẠI TỈNH NINH THUẬN - BÌNH THUẬN Đặng Quốc Khánh(1), Dương Văn Khảm(2), Dương Hải Yến(2), Nguyễn Văn Sơn(2) (1) Tổng Cục Khí tượng Thủy văn (2) Viện Khoa học Khí tượng Thủy văn Biến đổi khí hậu Ngày nhận bài: 25/4/2022; ngày chuyển phản biện: 26/4/2022; ngày chấp nhận đăng: 20/5/2022 Tóm tắt: Ninh Thuận - Bình Thuận tỉnh có điều kiện khí hậu khơ hạn Việt Nam Đây bất lợi lớn thiên nhiên phát triển nơng nghiệp nói riêng phát triển kinh tế nói chung tỉnh Trên sở chuỗi số liệu khí tượng thủy văn kịch biến đổi khí hậu (BĐKH), áp dụng phương pháp quan trắc hạn khí tượng thơng qua số ẩm (Moist index-MI) mơ hình thống kê, báo nghiên cứu tính toán biến động hạn hán, mức độ khắc nghiệt hạn khí tượng khả xảy hạn hán tương lai tỉnh Ninh Thuận - Bình Thuận Theo tính tốn trạm Cam Ranh, Phan Thiết, Phan Rang, số MI đại đa số năm nhỏ 0,4 (mức độ hạn nghiêm trọng) Tần suất xuất cấp độ từ hạn nhẹ đến hạn nghiêm trọng vào mùa khô chiếm tới 57,1% đến 92,9% tùy trạm Đặc biệt mùa mưa, trạm Hàm Tân Phan Thiết, hạn nhẹ chiếm đến gần 60% số năm nghiên cứu Theo kịch BĐKH, năm tới, thời gian xuất khơ hạn tỉnh khơng có nhiều biến động, nhiên mức độ khơ hạn có xu tăng lên cường độ tần suất Vì vậy, địa phương cần chủ động việc quy hoạch phát triển kinh tế - xã hội, có biện pháp thích ứng với hạn hán đặc biệt bối cảnh BĐKH nhằm hạn chế thấp thiệt hại thiên tai nói chung hạn hán nói riêng gây địa bàn tỉnh khu vực nghiên cứu Từ khóa: Chỉ số ẩm (MI), biến đổi khí hậu, hạn hán, Ninh Thuận, Bình Thuận Giới thiệu chung Hạn hán phân loại gồm có: Hạn khí tượng (thiếu hụt lượng mưa cán cân mưa - bốc hơi), hạn thủy văn (dòng chảy sông suối giảm rõ rệt, mực nước tầng chứa nước đất hạ thấp), hạn nông nghiệp (thiếu hụt nước mưa dẫn tới cân lượng nước thực tế nhu cầu nước trồng), hạn kinh tế - xã hội (thiếu hụt nguồn nước cấp cho hoạt động kinh tế - xã hội) [6, 14, 15] Việc đánh giá tổng hợp trạng, nguyên nhân, diễn biến xu loại hạn dựa số hạn ngưỡng hạn Hiện nay, nhiều số/hệ số hạn khác phát triển ứng dụng nước giới như: Chỉ số ẩm Ivanov (1948), Chỉ số khô Budyko Liên hệ tác giả: Đặng Quốc Khánh Email: khanhdangkhtc@gmail.com 36 TẠP CHÍ KHOA HỌC BIẾN ĐỔI KHÍ HẬU Số 22 - Tháng 6/2022 (1950), Chỉ số khơ Penman, Chỉ số gió mùa GMI, Chỉ số mưa chuẩn hóa SPI, Chỉ số Sazonov, Chỉ số Koloskov (1925), Hệ số khô, Hệ số cạn, Chỉ số Palmer (PDSI), Chỉ số độ ẩm trồng (CMI), Chỉ số cấp nước mặt (SWSI), Chỉ số RDI (Reclamation Drought Index), số hạn viễn thám VTCI, VCI LSWI [2, 3, 6, 9, 11, 12, 16, 17, 18, 20] Ở Việt Nam, có nhiều cơng trình nghiên cứu ứng dụng số hạn khác phục vụ việc đánh giá trạng, biến đổi, giám sát, cảnh báo dự báo hạn hán GS.TS Nguyễn Đức Ngữ Nguyễn Trọng Hiệu sử dụng số khô hạn K nghiên cứu hạn [6] PGS.TS Nguyễn Văn Thắng (2007, 2014) sử dụng số SPI, K, KBDI, PDSI nghiên cứu đánh giá, giám sát, cảnh báo dự báo hạn hán Việt Nam [9] GS.TS Trần Thục cs sử dụng số K, SPI, tỷ chuẩn lượng mưa (PN), thiếu hụt lượng mưa (D) số hạn thực tế (EDI), Penman để đánh giá xây dựng đồ hạn hán thiếu nước sinh hoạt khu vực Tây Nguyên Nam Bộ [13] Chỉ số hạn tích lũy sử dụng nhằm đánh giá xu biến đổi hạn hán khứ tương lai [10] PGS TS Dương Văn Khảm cs sử dụng số viễn thám để xây dựng đồ hạn hán cho Việt Nam [2] Thực tế, việc ứng dụng số ẩm (MI) để xác định điều kiện khơ hạn dự tính biến đổi tương lai theo kịch nhiều tác giả quan tâm [7, 8] Theo tác giả, nghiên cứu điều kiện khô hạn, MI phản ánh thiếu hụt nước mưa so với bốc thoát Sự thiếu hụt nước mưa khoảng thời gian coi điều kiện khơ hạn khí tượng, tính tốn thơng qua số MI Tuy nhiên, ngưỡng số MI cho xảy khô hạn tùy thuộc vào điều kiện khí hậu khu vực xác định [7, 8] Xuất phát từ ý tưởng đó, nghiên cứu thực đánh giá biến đổi số ẩm (MI) theo kịch BĐKH Các kết nghiên cứu nhằm góp phần cung cấp thơng tin quan trọng hạn hán phục vụ đánh giá tác động, tổn thương ứng phó với hạn hán tương lai Số liệu phương pháp nghiên cứu 2.1 Giới thiệu khu vực nghiên cứu Ninh Thuận - Bình Thuận tỉnh ven biển thuộc vùng Duyên hải Nam Trung Bộ (Hình 1), song thiên nhiên không thật ưu đãi cho người dân đây: Khơ hạn nắng gió nhắc đến biểu trưng khí hậu khắc nghiệt bất lợi lớn thiên nhiên phát triển nơng nghiệp nói riêng phát triển kinh tế - xã hội nói chung tỉnh Do có mùa khơ kéo dài - tháng, nên năm thời gian hạn, với mức độ khác Lượng mưa trung bình từ 1.000 - 1.400 mm có biến đổi mạnh vùng (từ 600 - 2.400 mm) Hệ thống sông suối ven biển Ninh Thuận - Bình Thuận chảy trực tiếp Biển Đông, ngắn dốc, mùa lũ nước lên nhanh xuống nhanh, mùa kiệt số sông nhỏ cạn nước Tuy nhiều năm qua tỉnh Ninh Thuận Bình Thuận đầu tư thủy lợi lớn, song đến số nơi cịn thiếu cơng trình hồ chứa chủ động tạo nguồn nên vào mùa khơ cịn tình trạng hạn hán, thiếu nước, chí thiếu nước gây hạn hán nghiêm trọng Ninh Thuận có lượng mưa trung bình năm khoảng 1.000 mm, song phân bố không Lượng mưa năm tập trung vào tháng, từ tháng - 12 Hạ lưu lưu vực sông Cái thuộc tỉnh xem vùng khô hạn nước Cộng thêm vào đấy, biến động mưa năm lại cao Những năm khô hạn, lượng mưa 60 70% trung bình Trong chuỗi số liệu gần 80 năm qua Phan Rang, có số lần xuất năm hạn (năm 1982 lượng mưa đạt 449 mm) Mưa xảy - tháng, lại mùa khô kéo dài - tháng, nên hạn hán nghiêm trọng nghiêm trọng Mùa khô 2015 - 2016, Bình Thuận trải qua đợt hạn hán khốc liệt vòng 10 năm qua Lượng mưa mùa mưa 2015 thiếu hụt 20 - 30% so với TBNN kết thúc sớm Cuối mùa mưa 2015, hồ thủy lợi - thủy điện tỉnh tích khoảng 60 - 70% dung tích thiết kế Đến đầu mùa khơ 2016, tổng dung tích hồ chứa tỉnh 111 tỷ m3 (đạt 48% dung tích thiết kế), hai hồ thủy điện mức 50% dung tích thiết kế, thấp kỳ đến 132 triệu m3 Đến mùa khô 2016, hồ thủy điện Đại Ninh, nguồn bổ sung nước quan trọng cho hệ thống thủy lợi Phan Rí - Phan Thiết (bao gồm huyện Hàm Thuận Bắc Bắc Bình) cịn chưa đầy 29% dung tích thiết kế Lượng nước tích trữ địa bàn tồn tỉnh cịn khoảng 87 triệu m3 (40% dung tích thiết kế) Theo Sở Nông nghiệp Phát triển nông thôn (NN&PTNT) tỉnh Ninh Thuận, sản xuất nông nghiệp, thiếu nước tưới nên diện tích phải dừng sản xuất vụ Đông - Xuân 2019 - 2020 7.874 Diện tích lúa bị thiệt hại hạn hán gây địa bàn tỉnh 378,7 Nhiều loại trồng lâu năm có nguy chết, giảm suất sản lượng thiếu nước tưới Đàn gia súc có nguy thiếu thức ăn, nước uống, phát sinh dịch bệnh… Bên cạnh đó, thiên tai hạn hán cịn tác động tới tình hình cháy rừng địa bàn tỉnh Tính đến ngày 19/3/2020, tồn tỉnh xuất 72 điểm cháy rừng, với tổng diện tích rừng bị cháy 45,7 TẠP CHÍ KHOA HỌC BIẾN ĐỔI KHÍ HẬU Số 22 - Tháng 6/2022 37 Đối với tỉnh Bình Thuận, khơng đảm bảo nguồn nước tưới nên từ đầu vụ Đông - Xuân 2019 - 2020, toàn tỉnh gieo trồng 32.000 trồng loại, giảm 20.000 so với kế hoạch Trong đó, diện tích lúa hoa màu đạt 12.500 ha, giảm đến 62% so với kế hoạch Đặc biệt, huyện canh tác nơng nghiệp trọng điểm tỉnh Bắc Bình Hàm Thuận Bắc, tình hình khơ hạn ảnh hưởng nghiêm trọng đến đời sống người dân, đặc biệt cung cấp nước sinh hoạt tưới rau màu [4, 5] Hình Khu vực nghiên cứu 2.2 Số liệu nghiên cứu (1) Số liệu quan trắc trạm: Nghiên cứu sử dụng số liệu khí tượng ngày từ năm 1993 đến năm 2020, bao gồm yếu tố: Nhiệt độ trung bình, nhiệt độ tối thấp, nhiệt độ tối cao, độ ẩm khơng khí, lượng mưa, số nắng, tốc độ gió khu vực nghiên cứu vùng phụ cận, gồm trạm: Phan Rang, Cam Ranh, Phan Thiết Hàm Tân (Nguồn: Trung tâm thông tin liệu Khí tượng Thủy văn) (2) Số liệu mơ phỏng: Kịch biến đổi khí hậu sử dụng nghiên cứu bao gồm kịch RCP4.5 RCP8.5 (3) Bản đồ: Địa hình, hành giao thơng thủy hệ tỉnh Ninh Thuận - Bình Thuận tỷ lệ 1/50.000 2.3 Phương pháp nghiên cứu - Sử dụng số ẩm (MI) để đánh giá mức độ hạn Trong báo này, nghiên cứu lựa chọn số ẩm (MI) Tổ chức Nông Lương Thế giới giới thiệu: Chỉ số ẩm (MI) định nghĩa tỷ số lượng mưa (X) với lượng bốc thoát tiềm (PET) [5, 6] Phân cấp mức độ hạn theo số MI xác định Bảng (1) Bảng Phân cấp mức độ hạn theo số hạn (MI) [14] 38 Chỉ số MI Cấp hạn MI < 0,4 Nghiêm trọng 0,4