Mục tiêu của nghiên cứu này là phân tích xu thế biến đổi lượng mưa của tỉnh Ninh Thuận trong bối cảnh Biến đổi khí hậu. Để thực hiện nghiên cứu này, tác giả đã sử dụng dữ liệu về lượng mưa tại các trạm đo mưa với thời gian tối thiểu là 24 năm và tối đa là 36 năm. Phương pháp chính được sử dụng trong nghiên cứu này là thống kê phi tham số (non-parametric method), cụ thể của phương pháp này là phân tích Mann-Kendall và phân tích độ dốc dựa vào phân tích Theil-Sen.
Tạp chí Phát triển Khoa học Cơng nghệ – Khoa học Tự nhiên, 5(1):957-968 Bài nghiên cứu Open Access Full Text Article Nghiên cứu xu biến đổi lượng mưa Ninh Thuận bối cảnh biến đổi khí hậu phương pháp phi tham số Nguyễn Hoàng Tuấn* , Trương Thanh Cảnh TÓM TẮT Use your smartphone to scan this QR code and download this article Mục tiêu nghiên cứu phân tích xu biến đổi lượng mưa tỉnh Ninh Thuận bối cảnh Biến đổi khí hậu Để thực nghiên cứu này, tác giả sử dụng liệu lượng mưa trạm đo mưa với thời gian tối thiểu 24 năm tối đa 36 năm Phương pháp sử dụng nghiên cứu thống kê phi tham số (non-parametric method), cụ thể phương pháp phân tích Mann-Kendall phân tích độ dốc dựa vào phân tích Theil-Sen Nghiên cứu tiến hành hỗ trợ phần mền ProULC 5.1 Makesen 1.0 Kết nghiên cứu đề tài giai đoạn tính lượng mưa năm địa bàn tỉnh Ninh Thuận có xu hướng tăng nhiều từ tháng đến tháng 11 giảm vào tháng Cùng với đó, thơng qua số xu hướng tăng giảm lượng mưa, nghiên cứu dự báo xu hướng tăng lượng mưa cho khu vực Kết dự báo cho thấy đến năm 2035 lượng mưa tăng trung bình 7,7% đến năm 2050 lượng mưa tăng 13,8% so với giai đoạn dự báo Kết nghiên cứu phản ánh thực trạng biến đổi lượng mưa bối cảnh biến đổi khí hậu với trạm có ý nghĩa thống kê (p Sgn(x j − xi ) = i f (x j − xi ) = −1 i f (x − x ) < j i 958 (2) [n(n − 1)(2n + 5)− 18 ( )( ) g ∑ p=1 t p t p − 2t p + ] (3) Trong đó, t p coi số lượng mối quan hệ với mẫu p, g số nhóm có giá trị liệu giống VAR(S) độ dốc ước lượng Sen sử dụng để ước tính xu hướng liệu chuỗi thời gian; Nếu giá trị S tích cực cho biết xu hướng tích cực giá trị âm cho biết xu hướng tiêu cực chuỗi liệu thời gian 29 Để tuân theo luật phân phối chuẩn trung bình 0, phương sai 1, số MK - Z tính theo cơng thức: S−1 , S>0 (4) Z = 0, S = S+1 Z= , S0 xu hướng xác định tăng; Z0) cho thấy thời điểm thuộc trạm lượng mưa luôn tăng Tại Quán Thẻ số Z tháng 10 2,455 1,999 độ lệch chuẩn 57,850 65,540; Sông Pha số Z tháng 1, 7, 8, 10 1,653; 1,985; 2,172; 2,359 1,938, độ lệch chuẩn 38,730 (tháng 1) 48,820 (tháng 7-10); Tân Mỹ số Z độ lệch chuẩn 2,657 67,740; Chỉ số Z Ba tháp tháng 1, 2,329; 1,756; 1,768 độ lệch chuẩn xác định 56,680 ; 57,530 61,670 ; cuối Cà Ná, số Z tháng 1,707, tháng 2,783 tháng 1,768 Độ lệch chuẩn tháng 53,310; 56,050 47,970 Còn số khác kết phân tích MK mô tả đầy đủ chi tiết Bảng Thông qua độ dốc Theil-Sen cho thấy 959 Tạp chí Phát triển Khoa học Cơng nghệ – Khoa học Tự nhiên, 5(1):957-968 Hình 1: Biểu đồ thể xu hướng biển đổi lượng mưa trạm Sông Phaa a 960 Nguồn : Kết nghiên cứu Trạm Kết Nhị Hà N S SD Z pValue Sen’s Slope Quán N Thẻ S SD Z pValue Sen’s Slope Sông N Pha S SD Z pValue Tháng Mùa 32 32 32 32 32 32 32 32 32 10 32 11 32 12 32 Xuân Hạ 32 32 Thu 32 Đông Mưa 32 28 Khô 10 32 46,66 0,664 0,253 -9 41,77 -0,192 0,424 48 55,72 0,844 0,199 -12 60,90 -0,181 0,428 -8 61,65 -0,114 0,455 -37 61,64 -0,584 0,280 54 61,65 0,86 0,195 61,64 0,065 0,474 -64 61,65 -1,022 0,153 47 61,64 0,746 0,228 61,63 0,081 0,468 19 61,29 0,294 0,384 28 61,65 0,438 0,331 -1 61,64 0,000 0,500 60 61,58 0,958 0,169 1,645 61,65 0,146 0,442 0,000 0,000 0,000 0,000 -0,067 -0,607 0,668 0,012 -2,037 1,452 0,081 0,000 0,848 -,317 0,000 0,603 2,018 0,300 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 32 143 57,85 2,455 0,007 -45 55,39 -0,794 0,213 21 62,86 0,318 0,375 73 62,86 1,145 0,126 16 64,54 0,232 0,408 31 64,53 0,465 0,321 48 64,54 0,728 0,233 64,54 0,047 0,481 -40 64,54 -0,604 0,273 130 64,54 1,999 0,023 86 64,54 1,317 0,094 57 64,53 0,868 0,193 39 64,53 0,589 0,278 22 64,54 0,325 0,372 122 64,54 1,875 0,0304 124 1,645 64,54 61,66 1,906 0,746 0,0283 0,228 0,000 0,000 0,000 0,018 0,427 0,367 0,439 0,025 -1,119 3,316 2,354 0,9 0,692 0,37 6,278 2,252 3,641 4,638 25 25 25 25 25 25 25 25 25 25 24 24 25 25 25 25 1,645 38,73 1,653 0,049 1,645 37,74 1,457 0,073 -1,645 42,3 -1,277 0,101 1,645 42,81 0,047 0,481 1,645 42,82 0,584 0,28 1,645 42,81 0,701 0,242 1,645 42,82 1,985 0,024 1,645 42,82 2,172 0,015 1,645 42,82 2,359 0,009 1,645 42,82 1,938 0,026 1,645 40,32 1,017 0,155 1,645 40,2 0,398 0,345 4,0 42,82 0,070 0,472 122 42,82 2,826 0,002 106 42,82 2,452 0,007 58 1,645 42,76 42,82 1,333 2,452 0,0913 0,0071 -16 61,65 -,243 0,404 1,645 61,63 0,438 0,331 47 106 136 1,645 61,67 2,189 0,0143 98 1,645 42,82 2,265 0,0117 961 Continued on next page Tạp chí Phát triển Khoa học Công nghệ – Khoa học Tự nhiên, 5(1):957-968 Bảng 2: Kết kiểm định Mann – Kendall Sen’s Slope lượng mưa tháng trạm quan trắc Ninh Thuận Ba Tháp Cà Ná 0,000 -,086 0,021 2,752 2,515 Table continued 5,628 5,861 9,194 11,970 2,954 0,321 0,737 17,250 29,260 2,817 30,750 21,160 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 82 60 S SD Z pValue Sen’s Slope N 1,645 67,74 2,657 0,004 1,645 52 1,461 0,072 1,645 69,97 0,186 0,426 1,645 76,17 0,683 0,247 1,645 76,46 0,955 0,17 -1,645 76,46 -0,929 0,177 1,645 76,46 0,536 0,296 1,645 76,45 0,445 0,328 -1,645 76,46 -0,667 0,252 -1,645 76,46 -0,379 0,352 1,645 76,46 0,955 0,17 1,645 76,4 1,152 0,125 46 76,46 0,589 0,278 -7 76,45 -,0785 0,469 40 76,46 0,51 0,305 165 76,45 2,145 0,016 1,645 76,46 1,059 0,145 1,645 76,46 0,772 0,22 0,000 0,000 0,000 0,205 0,977 -0,912 0,322 0,340 -0,653 -0,675 1,958 0,658 1,447 -,133 1,919 1,698 4,835 2,831 20 22 21 29 30 29 30 31 33 32 29 27 32 32 32 32 S SD Z pValue Sen’s Slope N 133 56,68 2,329 0,01 69 50,45 1,348 0,089 102 57,53 1,756 0,039 22 60,3 0,348 0,364 47 61,66 0,746 0,228 110 61,67 1,768 0,037 66 61,67 1,054 0,146 40 61,67 0,632 0,264 -32 61,67 -0,503 0,308 61,66 0,065 0,474 74 61,67 1,184 0,118 30 60,92 0,476 0,317 62 61,67 0,989 0,161 136 61,67 2,189 0,014 60 61,67 0,957 0,169 91 1,645 61,64 61,66 1,46 0,746 0,0721 0,228 0,000 0,000 0,000 0,000 0,594 1,524 0,739 0,323 -1,015 0,121 2,988 0,000 1,294 2,384 4,768 2,198 3,641 4,638 20 22 21 29 30 29 30 31 33 32 29 27 33 33 33 33 174 S SD Z pValue Sen’s Slope 37 30,09 1,196 0,116 33,6 0,119 0,453 -2 35,5 -0,03 0,488 21 53,3 0,375 0,354 55 56,05 0,963 0,168 92 53,31 1,707 0,044 157 56,05 2,783 0,003 13 58,84 0,204 0,419 46 64,54 0,697 0,243 110 61,67 1,768 0,039 42 53,31 0,769 0,221 -3 47,97 -,042 0,483 121 64,51 1,86 0,031 202 64,54 3,114 0,001 166 64,54 2,557 0,005 135 64,2 2,087 0,018 1,645 64,54 2,681 0,004 0,100 0,000 0,000 0,066 0,780 1,957 2,604 0,213 1,344 2,540 1,184 -,088 Ghi chú: S: Giá trị S; SD: Độ lệch chuẩn; Z: Giá trị chuẩn S, p-Value: Mức ý nghĩa; Sen’s Slope: Độ dốc Theil-Sen Nguồn: Kết nghiên cứu 2,625 6,485 47 136 1,645 61,67 2,189 0,0143 232 1,645 64,54 3,579 1,72E04 9,739 1,997 11,690 10,700 Tạp chí Phát triển Khoa học Công nghệ – Khoa học Tự nhiên, 5(1):957-968 962 Tân Mỹ Sen’s 0,000 Slope N 37 Tạp chí Phát triển Khoa học Công nghệ – Khoa học Tự nhiên, 5(1):957-968 Xu biến đổi lượng mưa theo mùa Kết kiểm định theo mùa (Bảng 2) cho thấy, lượng mưa theo mùa Ninh Thuận có xu hướng tăng, tăng nhiều vào tháng mùa Hạ mùa Thu, tăng chậm vào tháng mùa Xuân Đơng Cụ thể vào mùa Xn có trạm Cà Ná giá trị S 121, độ lệch chuẩn 64,51 số Z dương với Z 1,86 ; vào mùa Hạ, có ba trạm Sơng Pha, Ba Tháp Cà Ná có hệ số Z>0 tương ứng 2,826; 2,189 3,114 trị số S 122; 136 202, độ lệch chuẩn trạm 42,82; 61,67 64,54; giai đoạn mùa Thu, kết kiểm định cho thấy có Qn Thẻ, Sơng Pha Cà Ná, với p-value 0,030; 0,007 0,005, hệ số Z dương với Z 1,875; 2,452 2,557, hệ số S 122; 106 166 độ lệch chuẩn trạm 64,54; 42,82 65,54; cuối mùa Đơng, có trạm có hệ số p-value