Bản đồ trường nhiệt độ trung bình tháng (hình 3.1) dự báo cho tháng 3 đến tháng 8 có hạn dự báo lần lượt từ 1 đến 6 tháng với thời điểm dự báo ban đầu là ngày 13 tháng 02 năm 2012. Kết quả dự báo cho thấy trường nhiệt độ trung bình đã được mô phỏng khá tốt. Vào tháng 1, nhiệt độ không khí thấp ở phía Nam Trung Quốc do khối không khí lạnh lục địa hoạt động. Nhiệt độ tăng dần khi di chuyển xuống phía Nam. Vào các tháng tiếp theo, khi sự hoạt động của khối không khí lạnh lục địa suy yếu, nhiệt độ không khí đã tăng dần. Ngoài ra nhiệt độ không khí thấp hơn ở những nơi có địa hình cao cũng được mô phỏng khá hợp lý.
41
42
(a) (b)
(c) (d)
Hình 3.2 Nhiệt độ trung bình tháng (a), hệ số tương quan (b), sai số ME (c) và sai số quân phương (d)
Hình 3.2 (a) biểu diễn giá trị nhiệt độ trung bình tháng với trục x là tháng dự báo (từ tháng 1 đến tháng 10), hạn dự báo từ 0 đếm 6 tháng được biểu diễn theo trục y. Thang bên dưới là giá trị quan trắc từng tháng tương ứng, thang giá trị biểu diễn nhiệt độ đơn vị là độ C. Giá trị tại mỗi ô thể hiện giá trị trung bình cho một tháng cụ thể với 1 hạn dự báo xác định. Khi so sánh giá trị trung bình tháng kết quả mô hình và quan trắc chúng ta có thể thấy kết quả mô hình thấp hơn quan trắc ở hầu hết các tháng và hầu hết các hạn dự báo. Khi xem xét với các hạn dự báo, với các thời điểm dự báo ban đầu khác nhau, chúng ta có thể thấy kết quả dự báo nhiệt độ khá đồng nhất với các đích
43
dự báo, điều này cho thấy, mô hình cho kết quả dự báo khí hậu khá ổn định, ko bị ảnh hưởng nhiều bởi điều kiện thời gian ban đầu khác nhau. Cần phải lưu ý rằng, kết quả giá trị dự báo tại một tháng được so sánh với các hạn dự báo khác nhau cùng một cột giá trị.
Hình 3.2 (b) biểu diễn giá trị tương quan của chuỗi giá trị quan trắc và mô hình, trục x là tháng dự báo (từ tháng 1 đến tháng 10), trục y là hạn dự báo (từ 0 đến 6 tháng). Thang màu biểu diễn giá trị tương quan tương ứng (từ 0.0 đến 1.0) với mỗi tháng cụ thể và 1 hạn dự báo xác định. Khi xem xét giá trị tương quan giữa kết quả mô hình và quan trắc, chúng ta có thể thấy, đối với trường nhiệt độ, mô hình cho kết quả tương quan khá cao, hầu hết các dự báo cho kết quả tương quan từ 0,4 đến 0,8. Đối với dự báo cho tháng 6, tháng 7 và tháng 10 cho tương quan cao hơn cả. Ngoài ra, với hạn dự báo càng ngắn cho giá trị tương quan càng cao, điều này cho thấy, đối với trường nhiệt độ, hạn dự báo càng ngắn cho kết quả dự báo càng tốt hơn.
Hình 3.2 (c) biểu diễn tương tự với giá trị sai số trung bình ME. Kết quả cho thấy mô hình cho dự báo giá trị nhiệt độ thiên thấp ở tất cả các tháng và tất cả các hạn dự báo, giá trị sai số chủ yếu trong khoảng từ -1 đến -3.50C. Trong các tháng dự báo, các tháng mùa hè tháng 7, tháng 8, tháng 9 là cho sai số dự báo nhỏ hơn cả. Các tháng chuyển mùa, tháng 4, tháng 5 và tháng 10 cho sai số lớn hơn cả. Hạn dự báo càng dài cho sai số dự báo càng lớn. Với các thời điểm dự báo khác nhau không cho giá trị sai số đáng kể.
Hình 3.2 (d) biểu diễn tương tự với giá trị sai số quân phương (RMSE). Ở trường nhiệt độ, RMSE chủ yếu biến động trong khoảng từ 3.50C đến 5.50C. Giá trị sai số quân phương lớn do sai số dự báo giá trị từng ngày của mô hình sai khác với giá trị hàng ngày của quan trắc. Trong dự báo hạn dài tối đa đến 6 tháng, rất khó để dự báo chính xác giá trị nhiệt độ từng ngày cụ thể. Kết quả cũng cho thấy dự báo cho các tháng 6, 7, 8 với hạn dự báo từ 1 đến 3 tháng cho kết quả sai số quân phương là thấp nhất, Với thời điểm dự báo ban đầu là tháng 1, cho khoảng sai số quân phương là lớn nhất.
Từ những đánh giá cho trường nhiệt độ bề mặt cho ta thấy một số kết luận: - Đối với trường nhiệt độ bề mặt, dự báo cho các tháng mùa hè (từ tháng 6 đến tháng 9) với hạn dự báo ngắn (khoảng từ 1 đến 3 tháng) cho kết quả tốt nhất. Với các thời điểm dự báo khác nhau hầu như ít cho kết quả khác nhau với cùng đích dự báo.
44
- Với trường nhiệt độ, có thể thấy rõ ràng là mô hình cho dự báo thiên âm, với các hạn dự báo khác nhau đều cho sai số khá đồng nhất, ngoài ra chuỗi số liệu được phân tích theo giá trị từng ngày, vì mục đích của bài toán là dự báo hạn mùa, hàm phân bố các giá trị trong tháng của dự báo mới là đặc trưng quan trọng nhất, nó cho phép đánh giá dự báo có nắm bắt được phổ giá trị hay không, đó là quyết định có dự báo được các khoảng giá trị cực trị hay không. Sau đây chúng tôi sẽ phân tích hàm phân bố tần suất (hình 3.3) và phân bố đồng thời (hình 3.4) đối với trường nhiệt độ trung bình.
(a) (b) (c)
(d) (e)
Hình 3.3 Đồ thị phân bố tần suất trường nhiệt độ trung bình cho tháng 8 và tháng 4
Hình 3.3 biểu diễn phân bố tần suất (phân bố một chiều) trường nhiệt độ trung bình cho tháng 8 và tháng 4 với các hạn dự báo khác nhau. Trục x là giá trị nhiệt độ (đơn vị độ C), trục y là tần suất (đơn vị %). Đồ thị màu đỏ biểu diễn phân bố quan trắc, đồ thị màu xanh biểu diễn phân bố mô hình. Dự báo tháng 8 với hạn dự báo 1 tháng (a), 3 tháng (b), 6 tháng (c) và dự báo tháng 4 với hạn dự báo 1 tháng (d) và 3 tháng (e).
Dựa vào đồ thị phân bố với dự báo trường nhiệt độ tháng 8 với các hạn dự báo 1 tháng, 3 tháng, 6 tháng cho thấy phân bố nhiệt độ tháng 8 chủ yếu khoảng 24 đến 260
45
kỳ vọng khoảng 250C trong khi giá trị quan trắc có ngưỡng giá trị lớn hơn, khoảng 26 đến 300C, trung bình khoảng 280C. Đối với dạng phân bố của đồ thị dự báo có dạng chuẩn, đối xứng hơn, trong khi phân bố hàm quan trắc có dạng lệch phải. Giá trị quan trắc phân bố trong khoảng từ 18 đến 320C, trong khi dự báo giá trị trong khoảng 18 đến 300C. Với các hạn dự báo khác nhau hầu như có sự đồng nhất về dạng phân bố. Từ dạng phân bố của mô hình và quan trắc có thể thấy, đối với dự báo tháng 8, mô hình dự báo thiên thấp so với quan trắc từ 2-40C có tính hệ thống.
Phân bố tần suất giá trị nhiệt độ tháng 4 có sự khác biệt so với các tháng 6,7,8 khi mà phân bố giá trị quan trắc trải rộng trong khoảng từ 16 đến 340C và chủ yếu phân bố trong khoảng 24 đến 300C, trong khi phân bố giá trị quan trắc trong tháng 8 chủ yếu trong khoảng 28 đến 300C. Sự khác biệt trong phân bố tần suất tháng 4 và tháng 8 cho phép khảo sát về khả năng dự báo với các phân bố tần suất khác nhau.
Sự phân bố tần suất giữa 3 mùa hè, tháng 7, tháng 8, tháng 9 (xem phụ lục) với các hạn dự báo khác nhau cũng được cho thấy nền nhiệt độ dự báo đều thấp hơn quan trắc, giá trị kỳ vọng của dự báo đều thấp hơn quan trắc khoảng 2 đến 40C. Về dạng phân bố nhiệt độ dự báo tháng 7, 8, 9 hầu như tuân theo dạng phân bố chuẩn, trong khi phân bố quan trắc thường có dạng chuẩn (tháng 8, 9) hoặc lệch phải (tháng 7). Giá trị dự báo hầu hết rơi vào khoảng 18 đến 300C, trong khi giá trị quan trắc phổ biến trong khoảng 20-320C. Ngoài ra, khi xem xét độ nhọn của phân bố, chúng ta có thể thấy dự báo nắm bắt khá tốt. Phân bố tần suất cao nhất khoảng 35% đối với tháng 7, 38% đối với tháng 8 và khoảng 40% đối với tháng 9. Từ các kết quả trên cũng cho thấy mô hình cho giá trị dự báo nắm bắt khá tốt phân bố giá trị quan trắc, từ đó củng cố thêm nhận định sai số có tính hệ thống của mô hình.
46
(d) (e)
Hình 3.4: Phân bố đồng thời giá trị nhiệt độ tháng 8 và tháng 4
Tiếp theo, trong hình 3.4 chúng tôi sẽ xem xét phân bố đồng thời hay phân bố tần suất 2 chiều của trường nhiệt độ trung bình. Trục x là giá trị mô hình, trục y là giá trị quan trắc, giá trị tại mỗi ô là xác suất (%) có điều kiện với ngưỡng giá trị tương ứng. Phân bố đồng thời nhiệt độ tháng 4 và tháng 8 cho thấy các dự báo luôn cho giá trị thiên âm tuy nhiên các giá trị dự báo và quan trắc tập trung xung quanh đường đẳng trị. Đường tần suất cực đại luôn song song với đường đẳng trị cho thấy sai số của mô hình là đồng nhất và mang tính hệ thống. Cực đại tần suất cũng cho thấy với khoảng nhiệt độ từ 24 đến 280C đối với tháng 4 và 26 đến 300C đối với tháng 8 mô hình nắm bắt phân bố khá tốt. So với tháng 4, các giá trị nhiệt độ trung bình tháng 8 có mức độ tập trung cao hơn, do nhiệt độ vào mùa hè đồng nhất giữa các trạm ở miền Bắc và các trạm ở miền Nam, trong khi vào tháng 4 có sự chênh lệch nhiệt độ rõ rệt hơn. Ngoài ra giá trị tần suất cực đại trong tháng 8 cũng cao hơn so với tháng 4, cho thấy độ tin cậy dự báo tháng 8 là tốt hơn. Với hạn dự báo khác nhau hầu như không có sự khác biệt. Với hạn dự báo càng dài, độ tán của dự báo càng rộng cho thấy độ tin cậy của dự báo giảm khi hạn dự báo tăng.
Phân bố đồng thời với các tháng 7, 8, 9 (xem phụ lục) cũng cho thấy với dự báo các tháng khác nhau, mô hình đều cho dự báo khá đồng nhất, đường cực đại tần suất gần như song song với đường đẳng trị, nhiệt độ dự báo. Ở cả 3 tháng dự báo chúng ta đều có thể nhận thấy với hạn dự báo càng dài, độ tán của giá trị mô hình càng rộng, điều này có nghĩa thông tin dự báo tản mạn hơn và độ tin cậy thấp hơn.
b) Lượng mưa tháng
Bản đồ tổng lượng mưa tháng được biểu diễn trong hình 3.5 dự báo cho tháng 3 đến tháng 8 với hạn dự báo lần lượt từ 1 đến 6 tháng với thời điểm điểm dự báo ngày 13 tháng 02 năm 2012.
47
48
(a) (b)
Hình 3.6 : Lượng mưa trung bình tháng (a), hệ số tương quan (b), sai số ME (c) và sai số quân phương (d)
Trong dự báo thời tiết cũng như khí hậu, trường mưa là một yếu tố rất khó dự báo, nhất là dự báo bằng mô hình. Giá trị dự báo mưa cho thấy, mô hình cũng đã dự báo được những tháng mưa có mưa và những tháng không mưa, tuy nhiên giá trị tương quan giữa dự báo và quan trắc là không lớn. Với tháng 4 là tháng chuyển mùa, bắt đầu từ mùa khô sang mùa mưa, hầu như mô hình không có khả năng dự báo. Tương tự với các dự báo hạn trên 3 tháng với đích dự báo là các tháng mùa hè cũng cho tương quan rất nhỏ.
49
Với dự báo mưa hầu hết các dự báo có xu hướng thiên âm với các tháng có lượng mưa nhiều và không rõ rệt xu hướng đối với các tháng có lượng mưa nhỏ. RMSE khá đồng nhất với các tháng, điều này cho thấy dự báo mưa hạn mùa phụ thuộc chủ yếu và đích dự báo, đối với dự báo mưa RMSE lớn vào các tháng mùa hè và nhỏ vào các tháng đầu năm, điều này không có nghĩa là dự báo các tháng mùa hè thì xấu hơn các tháng đầu năm. Nguyên nhân chủ yếu là do lượng mưa nhiều hơn vào các tháng mùa hè, còn vào các tháng đầu năm lượng mưa nhỏ.
(a) (b) (c)
(d) (e)
Hình 3.7 Đồ thị phân bố tần suất trường mưa cho tháng 8 và tháng 4
Phân bố tần suất có thể thấy mô hình chưa dự báo tốt trường mưa, khi dạng phân bố quan trắc và dự báo có sự khác biệt. Mô hình dự báo không mưa (lượng mưa bằng 0mm/ngày) với tần suất khoảng 15% trong khi quan trắc có tần suất trên 50%. Trong khi quan trắc lượng mưa trải đều trong khoảng từ 3 đến khoảng 30mm/ngày với tần suất thấp thì trong dự báo, hầu hết giá trị trong khoảng từ 0 đến 10mm/ngày với tần suất cao hơn. Với các hạn dự báo khác nhau đều cho kết quả dự báo khá đồng nhất.
Đối với tháng 4, dự báo giá trị chủ yếu trong khoảng giá trị từ 1 đến 5 mm/ngày với tần suất lớn, trong khi quan trắc giá trị có độ tán rộng, với tần suất khá nhỏ, sự
50
chênh lệch tần suất giữa quan trắc và dự báo là khá rõ rệt. So với tháng 8 thì dự báo phân bố tần suất tháng 4 kém hơn.
Kết quả cho thấy dự báo sự phân bố 3 tháng 7, 8, 9 (xem phụ lục) không có sự khác biệt rõ nét, dự báo lượng mưa chủ yếu ở ngưỡng thấp với tần suất lớn, trong khi giá trị quan trắc có độ tán lớn hơn, trải rộng và tần suất nhỏ hơn dự báo. Dự báo không mưa ở tháng 9 là tốt nhất so với tháng 7 và tháng 8. Dạng phân bố của hầu hết các dự báo đều có dạng lệch phải. Sự khác biệt với các hạn dự báo khác nhau không được thể hiện rõ.