ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN –––––––––––––– LƢƠNG NHƢ NGUYÊN ÁP DỤNG PHƢƠNG PHÁP UMOS ĐỂ CẢI THIỆN DỰ BÁO PHA MƢA TRÊN KHU VỰC VIỆT NAM LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – Năm 2016 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN –––––––––––––––– LƢƠNG NHƢ NGUYÊN ÁP DỤNG PHƢƠNG PHÁP UMOS ĐỂ CẢI THIỆN DỰ BÁO PHA MƢA TRÊN KHU VỰC VIỆT NAM Chuyên ngành : Khí tƣợng Khí hậu học Mã số: 60440222 LUẬN VĂN THẠC SĨ KHOA HỌC NGƢỜI HƢỚNG DẪN KHOA HỌC: TS BÙI MINH TĂNG Hà Nội – Năm 2016 LỜI CAM ĐOAN Tôi xin cam đoan Luận văn “Áp dụng phƣơng pháp UMOS để cải thiện dự báo pha mƣa khu vực Việt Nam” cơng trình nghiên cứu cá nhân thực dƣới hƣớng dẫn khoa học Tiến sỹ Bùi Minh Tăng, khơng chép từ cơng trình nghiên cứu ngƣời khác Số liệu kết Luận văn chƣa đƣợc công bố công trình khoa học khác Các thơng tin thứ cấp sử dụng Luận văn có nguồn gốc rõ ràng, đƣợc trích dẫn đầy đủ, trung thực quy cách Tơi xin hồn tồn chịu trách nhiệm tính xác thực nguyên Luận văn./ Hà Nội, ngày 20 tháng 12 năm 2016 Tác giả Luận văn Lƣơng Nhƣ Nguyên LỜI CẢM ƠN Lời đầu tiên, xin trân trọng cảm ơn Tiến sỹ Bùi Minh Tăng định hƣớng nghiên cứu phƣơng pháp luận cho tơi suốt q trình thực nghiên cứu Luận văn Thạc sỹ Trong trình nghiên cứu học tập Khoa Khí tƣợng, Thủy văn Hải dƣơng học, Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội tơi có hội đƣợc tiếp thu kiến thức chuyên sâu chuyên ngành Khí tƣợng Khí hậu học Qua đó, giúp tơi có đƣợc kiến thức chun mơn nhƣ kinh nghiệm suốt q trình học tập, tạo động lực nghiên cứu khoa học, phục vụ hiệu trình nghiên cứu, thực hồn thiện Luận văn Thạc sỹ tơi Tơi xin trân trọng cảm ơn đồng chí Lãnh đạo, thầy giáo, cô giáo cán Khoa Khí tƣợng, Thủy văn Hải dƣơng học tận tình hƣớng dẫn, giúp đỡ tơi điều kiện suốt trình học tập thực Luận văn Xin chân thành cảm ơn đồng chí Lãnh đạo cán Trung tâm Dự báo khí tƣợng thủy văn Trung ƣơng tạo điều kiện cho tơi tham gia khóa đào tạo Thạc sỹ, đề tài nghiên cứu khoa học Trung tâm Dự báo khí tƣợng thủy văn Trung ƣơng chủ trì cung cấp thông tin, số liệu quan trắc, dự báo, tài liệu tận tình giúp đỡ tơi trình thực Luận văn Trân trọng cảm ơn quan tâm, giúp đỡ bạn bè, đồng nghiệp gia đình ln sát cánh, động viên, giúp đỡ tơi q trình học tập hồn thành Luận văn./ Hà Nội, ngày 20 tháng 12 năm 2016 Tác giả Luận văn Lƣơng Nhƣ Nguyên MỤC LỤC MỤC LỤC DANH MỤC BẢNG DANH MỤC HÌNH VẼ DANH MỤC KÝ HIỆU VIẾT TẮT MỞ ĐẦU CHƢƠNG I: TỔNG QUAN 13 1.1 Tổng quan phƣơng pháp thống kê sau mơ hình 13 1.2 Tình hình nghiên cứu nƣớc ngồi 14 1.3 Tình hình nghiên cứu nƣớc 18 CHƢƠNG II: SỐ LIỆU VÀ PHƢƠNG PHÁP NGHIÊN CỨU 23 2.1 Số liệu 23 2.1.1 Số liệu quan trắc 23 2.1.2 Số liệu dự báo 24 2.2 Phƣơng pháp nghiên cứu 25 2.2.1 Phương pháp UMOS 25 2.2.2 Phương pháp GMOS 31 2.3 Áp dụng phƣơng pháp UMOS cho toán dự báo POP 34 2.3.1 Xác định yếu tố dự báo 34 2.3.2 Xác định nhấn tố dự báo 34 2.3.3 Xây dựng phương trình dự báo POP 37 2.4 Phƣơng pháp đánh giá kết dự báo POP 38 2.4.1 Điểm số Brier 38 2.4.2 Biểu đồ tin cậy 40 2.4.3 Đường đặc trưng hoạt động ROC số ROCA 42 2.5 Xây dựng hệ thống GMOS để đƣa dự báo UMOS từ điểm trạm lƣới phân giải cao 43 CHƢƠNG III: KẾT QUẢ, ĐÁNH GIÁ 46 3.1 Kết xây dựng phƣơng trình dự báo 46 3.2 Kết đánh giá dự báo 49 3.2.1 Đánh giá số kỹ 49 3.2.2 Biểu đồ độ tin cậy 54 3.3 Thử nghiệm kết thị dƣới dạng lƣới GMOS cho dự báo POP 62 KẾT LUẬN VÀ KIẾN NGHỊ 65 TÀI LIỆU THAM KHẢO i PHỤ LỤC 1: Danh sách trạm quan trắc khí tƣợng bề mặt thuộc Việt Nam đƣợc sử dụng Luận văn v PHỤ LỤC 2: Danh sách nhân tố dự báo dự tuyển từ mô hình WRF-ARW ix PHỤ LỤC 3: Minh họa biến đổi phƣơng trình dự báo RV trạm Bắc Quang hạn dự báo 24 với mô hình WRF-ARW vào mùa đơng xiv PHỤC LỤC 4: Minh họa biến đổi phƣơng trình dự báo RV trạm Bắc Quang hạn dự báo 24 với mơ hình WRF-ARW vào mùa hè xv PHỤ LỤC 5: Minh họa biến đổi phƣơng trình dự báo RV trạm Bạch Long Vỹ hạn dự báo 24 với mơ hình WRF-ARW xvii PHỤ LỤC 6: Minh họa biến đổi phƣơng trình dự báo RV trạm Bãi Cháy hạn dự báo 24 với mơ hình WRF-ARW xix PHỤ LỤC 7: Minh họa biến đổi phƣơng trình dự báo RV trạm Buôn Ma Thuột hạn dự báo 24 với mơ hình WRF-ARW xxi PHỤ LỤC 8: Minh họa biến đổi phƣơng trình dự báo RV trạm Cà Mau hạn dự báo 24 với mơ hình WRF-ARW xxiii PHỤ LỤC 9: Minh họa biến đổi phƣơng trình dự báo RV trạm Cần Thơ hạn dự báo 24 với mơ hình WRF-ARW xxv PHỤ LỤC 10: Minh họa biến đổi phƣơng trình dự báo RV trạm Hà Nội hạn dự báo 24 với mơ hình WRF-ARW xxvii PHỤ LỤC 11: Minh họa biến đổi phƣơng trình dự báo RV trạm Vinh hạn dự báo 24 với mơ hình WRF-ARW xxix DANH MỤC BẢNG Bảng 2.1: Số lƣợng trạm quan trắc Đài khu vực 23 Bảng 2.2: Cấu hình chi tiết mơ hình WRF-ARW 25 Bảng 3.1: Tần suất tuyển chọn tần suất tuyển chọn nhân tố phƣơng trình POP hạn 24 dung lƣợng mẫu tăng vào mùa đông mơ hình WRF-ARW 46 Bảng 3.2: Tần suất tuyển chọn tần suất tuyển chọn nhân tố phƣơng trình POP hạn 24 dung lƣợng mẫu tăng vào mùa hè mô hình WRF-ARW 47 Bảng 3.3: Tần suất tuyển chọn tần suất tuyển chọn nhân tố phƣơng trình POP hạn 48 dung lƣợng mẫu tăng vào mùa đơng mơ hình WRF-ARW 48 Bảng 3.4: Tần suất tuyển chọn tần suất tuyển chọn nhân tố phƣơng trình POP hạn 48 dung lƣợng mẫu tăng vào mùa hè mơ hình WRF-ARW 48 Bảng 3.5: Kết tính tốn số BSS cho dự báo POP hạn 24 từ UMOS DMO mơ hình WRF-ARW từ DMO mơ hình GSM trung bình khu vực tồn Việt Nam năm 2014 50 Bảng 3.6: Kết tính tốn số BSS cho dự báo POP hạn 24 từ UMOS DMO mơ hình WRF-ARW từ DMO mơ hình GSM trung bình khu vực toàn Việt Nam năm 2015 51 Bảng 3.7: Kết tính tốn số BSS cho dự báo POP hạn 48 từ UMOS DMO mơ hình WRF-ARW từ DMO mơ hình GSM trung bình khu vực toàn Việt Nam năm 2014 51 Bảng 3.8: Kết tính tốn số BSS cho dự báo POP hạn 48 từ UMOS DMO mơ hình WRF-ARW từ DMO mơ hình GSM trung bình khu vực toàn Việt Nam năm 2015 52 Bảng 3.9: Kết tính toán số ROCA cho dự báo POP hạn 24 từ UMOS DMO mơ hình WRF-ARW từ DMO mơ hình GSM trung bình khu vực toàn Việt Nam năm 2014 52 Bảng 3.10: Kết tính tốn số ROCA cho dự báo POP hạn 24 từ UMOS DMO mơ hình WRF-ARW từ DMO mơ hình GSM trung bình khu vực toàn Việt Nam năm 2015 53 Bảng 3.11: Kết tính tốn số ROCA cho dự báo POP hạn 48 từ UMOS DMO mơ hình WRF-ARW từ DMO mơ hình GSM trung bình khu vực toàn Việt Nam năm 2014 53 Bảng 3.12: Kết tính tốn số ROCA cho dự báo POP hạn 48 từ UMOS DMO mơ hình WRF-ARW từ DMO mơ hình GSM trung bình khu vực toàn Việt Nam năm 2015 54 DANH MỤC HÌNH VẼ Hình 2.1: Miền định hƣớng dự báo mơ hình WRF-ARW 24 Hình 2.2: Biểu đồ tin cậy 41 Hình 2.3: Đƣờng đặc trƣng hoạt động 43 Hình 3.1: Biểu đồ tin cậy cho dự báo POP theo UMOS cho mơ hình WRFARW vào mùa hè năm 2014 cho hạn dự báo 24 khu vực nghiên cứu…………………………………………… ……………………………55 Hình 3.2: Biểu đồ tin cậy cho dự báo POP theo UMOS cho mơ hình WRFARW vào mùa đơng năm 2014 cho hạn dự báo 24 khu vực nghiên cứu 56 Hình 3.3: Biểu đồ tin cậy cho dự báo POP theo UMOS cho mơ hình WRFARW vào mùa hè năm 2015 cho hạn dự báo 24 khu vực nghiên cứu………………………………………………………………………… 57 Hình 3.4: Biểu đồ tin cậy cho dự báo POP theo UMOS cho mơ hình WRFARW vào mùa đơng năm 2015 cho hạn dự báo 24 khu vực nghiên cứu 58 Hình 3.5: Biểu đồ tin cậy cho dự báo POP theo UMOS cho mơ hình WRFARW vào mùa hè năm 2014 cho hạn dự báo 48 khu vực nghiên cứu………………………………………………………………………… 59 Hình 3.6: Biểu đồ tin cậy cho dự báo POP theo UMOS cho mơ hình WRFARW vào mùa đông năm 2014 cho hạn dự báo 48 khu vực nghiên cứu 60 Hình 3.7: Biểu đồ tin cậy cho dự báo POP theo UMOS cho mơ hình WRFARW vào mùa hè năm 2015 cho hạn dự báo 48 khu vực nghiên cứu………………………………………………………………………… 61 Hình 3.8: Biểu đồ tin cậy cho dự báo POP theo UMOS cho mơ hình WRFARWvào mùa đơng năm 2015 cho hạn dự báo 48 khu vực nghiên cứu 62 Hình 3.9: Bản đồ dự báo hạn 24 cho xác suất xảy mƣa từ 00Z ngày 04/08/2015 theo GMOS với dự báo từ UMOS mơ hình WRF-ARW quan trắc tƣơng ứng 63 Hình 3.10: Bản đồ dự báo hạn 48 cho xác suất xảy mƣa từ 00Z ngày 04/08/2015 theo GMOS với dự báo từ UMOS mô hình WRF-ARW quan trắc tƣơng ứng 64 2015090700 516 0.0578 0.628 + 0.121*du87 + -0.046*rain + 0.061*gradttdy500 2015092400 528 0.0492 0.654 + 0.073*du87 + -0.061*rain + -0.034*tdmean xvi PHỤ LỤC Minh họa biến đổi phƣơng trình dự báo RV trạm Bạch Long Vỹ hạn dự báo 24 với mơ hình WRF-ARW 2015053100 0.304 2015060700 0.3001 2015061300 0.2979 2015062100 0.2919 2015062800 0.2925 2015070200 0.2942 2015071200 0.2955 2015071900 0.2997 2015072400 0.2897 2015080200 0.3014 2015080900 0.303 2015081600 0.3024 2015082300 0.3045 2015082800 0.3045 2015090600 0.302 2015091300 0.2878 2015092000 0.2935 2015092700 0.2932 2015100400 0.2789 2015101100 0.288 2015101800 0.2868 0.302 + 0.181*rhmean + 0.093*obs + 0.096*rain + 0.072*u850 + -0.056*graduy500 + -0.041*gradtx500 0.302 + 0.183*rhmean + 0.091*obs + 0.096*rain + 0.072*u850 + -0.056*graduy500 + -0.042*gradtx500 0.303 + 0.182*rhmean + 0.090*obs + 0.095*rain + 0.071*u850 + -0.056*graduy500 + -0.041*gradtx500 0.299 + 0.187*rhmean + 0.089*obs + 0.092*rain + 0.070*u850 + -0.050*graduy500 0.295 + 0.195*rhmean + 0.090*obs + 0.073*rain + 0.073*u850 + -0.048*graduy500 0.296 + 0.193*rhmean + 0.090*obs + 0.073*rain + 0.072*u850 + -0.048*graduy500 0.292 + 0.195*rhmean + 0.090*obs + 0.072*rain + 0.071*u850 + -0.049*graduy500 0.296 + 0.192*rhmean + 0.090*obs + 0.073*rain + 0.073*u850 + -0.055*graduy500 + -0.038*gradtx500 0.293 + 0.193*rhmean + 0.089*obs + 0.075*rain + 0.073*u850 + -0.044*graduy500 0.292 + 0.194*rhmean + 0.090*obs + 0.069*rain + 0.073*u850 + -0.045*graduy500 0.290 + 0.195*rhmean + 0.091*obs + 0.069*rain + 0.072*u850 + -0.044*graduy500 0.289 + 0.196*rhmean + 0.092*obs + 0.068*rain + 0.072*u850 + -0.046*graduy500 0.287 + 0.197*rhmean + 0.092*obs + 0.068*rain + 0.072*u850 + -0.045*graduy500 0.285 + 0.196*rhmean + 0.094*obs + 0.067*rain + 0.071*u850 + -0.046*graduy500 0.284 + 0.197*rhmean + 0.091*obs + -0.072*u850 + 0.065*rain + -0.046*graduy500 0.294 + 0.148*rhmean + 0.109*obs + 0.072*rain + 0.050*u850 + -0.020*graduy500 0.298 + 0.370*rhmean + 0.106*obs + 0.063*rain + 0.044*u850 + 0.238*ttdmean 0.294 + 0.365*rhmean + 0.111*obs + 0.063*rain + 0.035*u850 + 0.239*ttdmean 0.280 + 0.323*rhmean + 0.109*obs + 0.236*ttdmean + 0.069*rh850 + 0.060*rain 0.273 + 0.332*rhmean + 0.115*obs + 0.230*ttdmean + 0.064*rh850 + 0.030*wind700 0.272 + 0.328*rhmean + 0.115*obs + 0.228*ttdmean + xvii 2015102500 0.2989 2015110100 0.2999 2015110800 0.3 2015111200 0.297 2015112200 0.2968 2015112900 0.2954 2015120600 0.2979 2015121100 0.2969 2015122000 0.295 2015122700 0.2953 0.063*rh850 + 0.030*wind700 0.276 + 0.289*rhmean + 0.111*obs + 0.192*ttdmean + 0.064*rh850 + 0.053*wind700 + -0.040*gradttdy500 0.275 + 0.290*rhmean + 0.112*obs + 0.194*ttdmean + 0.064*rh850 + 0.054*wind700 + -0.040*gradttdy500 0.273 + 0.114*obs + 0.288*rhmean + 0.194*ttdmean + 0.064*rh850 + 0.055*wind700 + -0.041*gradttdy500 0.273 + 0.112*obs + 0.286*rhmean + 0.193*ttdmean + 0.065*rh850 + 0.055*wind700 + -0.042*gradttdy500 0.274 + 0.113*obs + 0.288*rhmean + 0.196*ttdmean + 0.066*rh850 + 0.054*wind700 + -0.042*gradttdy500 0.275 + 0.116*obs + 0.283*rhmean + 0.193*ttdmean + 0.066*rh850 + 0.054*wind700 + -0.042*gradttdy500 0.277 + 0.117*obs + 0.286*rhmean + 0.196*ttdmean + 0.067*rh850 + 0.054*wind700 + -0.039*gradttdy500 0.277 + 0.116*obs + 0.287*rhmean + 0.197*ttdmean + 0.067*rh850 + 0.055*wind700 + -0.037*gradttdy500 0.274 + 0.119*obs + 0.283*rhmean + 0.194*ttdmean + 0.065*rh850 + 0.054*wind700 + -0.037*gradttdy500 0.273 + 0.116*obs + 0.287*rhmean + 0.196*ttdmean + 0.065*rh850 + 0.056*wind700 + -0.039*gradttdy500 xviii PHỤ LỤC Minh họa biến đổi phƣơng trình dự báo RV trạm Bãi Cháy hạn dự báo 24 với mơ hình WRF-ARW 2015053100 0.3023 2015060700 0.3026 2015061300 0.3013 2015062100 0.3021 2015062800 0.3061 2015070200 0.3102 2015071200 0.3123 2015071900 0.3107 2015072400 0.309 2015080200 0.3115 2015080900 0.3115 2015081600 0.3092 2015082300 0.3111 2015083000 0.3123 2015090600 0.3118 2015091300 2015092000 2015092700 0.2683 0.2651 0.256 2015100400 0.2593 2015101100 0.2692 2015101800 0.2708 0.502 + 0.210*rhmean + 0.084*rain + -0.066*u700 + 0.055*obs + -0.050*dq87 0.501 + 0.210*rhmean + 0.084*rain + -0.065*u700 + 0.056*obs + -0.050*dq87 0.500 + 0.210*rhmean + 0.084*rain + -0.066*u700 + 0.056*obs + -0.051*dq87 0.500 + 0.211*rhmean + 0.084*rain + -0.067*u700 + 0.055*obs + -0.049*dq87 0.498 + 0.213*rhmean + 0.082*rain + -0.069*u700 + 0.056*obs + -0.050*dq87 0.498 + 0.213*rhmean + 0.082*rain + -0.069*u700 + 0.056*obs + -0.050*dq87 0.497 + 0.213*rhmean + 0.084*rain + -0.068*u700 + 0.055*obs + -0.052*dq87 0.498 + 0.214*rhmean + 0.083*rain + -0.068*u700 + 0.053*obs + -0.052*dq87 0.499 + 0.215*rhmean + 0.084*rain + -0.068*u700 + 0.050*obs + -0.052*dq87 0.498 + 0.215*rhmean + 0.080*rain + -0.068*u700 + 0.051*obs + -0.052*dq87 0.497 + 0.214*rhmean + 0.081*rain + -0.066*u700 + 0.051*obs + -0.053*dq87 0.497 + 0.214*rhmean + 0.082*rain + -0.066*u700 + 0.051*obs + -0.052*dq87 0.496 + 0.215*rhmean + 0.081*rain + -0.068*u700 + 0.051*obs + -0.051*dq87 0.497 + 0.217*rhmean + 0.082*rain + -0.068*u700 + 0.049*obs + -0.051*dq87 0.496 + 0.215*rhmean + 0.085*rain + -0.066*u700 + 0.049*obs + -0.052*dq87 0.466 + 0.118*rhmean + 0.105*rain + 0.092*obs 0.466 + 0.118*rhmean + 0.106*rain + 0.090*obs 0.520 + 0.132*rhmean + 0.118*rain + -0.055*dq87 0.317 + 0.121*obs + -0.153*showalter_index + 0.088*wind700 + -0.035*advttd500 0.322 + 0.122*obs + 0.027*v850 + 0.085*wind700 + 0.130*showalter_index + -0.038*advttd500 0.321 + 0.123*obs + 0.028*v850 + 0.085*wind700 + 0.129*showalter_index + -0.039*advttd500 xix 2015102500 0.2824 2015110100 0.284 2015110800 0.2829 2015111200 0.2831 2015112200 0.2816 2015112900 0.283 2015120600 0.2848 2015121300 0.2851 2015122000 0.2862 2015122700 0.2837 0.303 + 0.120*obs + 0.066*v850 + 0.087*wind700 + 0.078*showalter_index + -0.046*advttd500 0.302 + 0.121*obs + 0.065*v850 + 0.089*wind700 + 0.077*showalter_index + -0.046*advttd500 0.302 + 0.123*obs + 0.065*v850 + 0.088*wind700 + 0.077*showalter_index + -0.046*advttd500 0.303 + 0.123*obs + 0.066*v850 + 0.087*wind700 + 0.077*showalter_index + -0.047*advttd500 0.303 + 0.121*obs + 0.067*v850 + 0.088*wind700 + 0.077*showalter_index + -0.046*advttd500 0.302 + 0.122*obs + 0.064*v850 + 0.089*wind700 + 0.077*showalter_index + -0.047*advttd500 0.303 + 0.123*obs + 0.065*v850 + 0.089*wind700 + 0.078*showalter_index + -0.047*advttd500 0.303 + 0.123*obs + 0.065*v850 + 0.090*wind700 + 0.080*showalter_index + -0.046*advttd500 0.302 + 0.123*obs + 0.065*v850 + 0.089*wind700 + 0.082*showalter_index + -0.046*advttd500 0.301 + 0.120*obs + 0.067*v850 + 0.088*wind700 + 0.080*showalter_index + -0.047*advttd500 xx PHỤ LỤC Minh họa biến đổi phƣơng trình dự báo RV trạm Bn Ma Thuột hạn dự báo 24 với mơ hình WRF-ARW 2015053100 2015060700 2015061300 2015062100 2015062500 2015070200 2015071200 2015071900 2015072400 2015080200 2015080600 2015081600 2015082300 2015083000 2015090600 2015091300 2015092000 2015092700 0.2091 0.2098 0.2069 0.208 0.2104 0.2083 0.209 0.2107 0.2103 0.2127 0.2114 0.2127 0.212 0.2112 0.2067 0.1979 0.1981 0.1918 2015100400 0.3001 2015101100 0.3062 2015101800 0.3074 2015102500 0.3124 2015110100 0.3129 2015110600 0.319 2015111200 0.3169 2015112200 0.3086 2015112900 0.306 2015120600 0.3029 0.095 + 0.295*eptmean + 0.058*rain + 0.064*dq75 + 0.035*obs 0.091 + 0.295*eptmean + 0.057*rain + 0.061*dq75 + 0.036*obs 0.092 + 0.293*eptmean + 0.057*rain + 0.062*dq75 + 0.038*obs 0.092 + 0.295*eptmean + 0.057*rain + 0.061*dq75 + 0.037*obs 0.090 + 0.297*eptmean + 0.056*rain + 0.063*dq75 + 0.037*obs 0.091 + 0.298*eptmean + 0.056*rain + 0.061*dq75 + 0.036*obs 0.092 + 0.300*eptmean + 0.054*rain + 0.060*dq75 + 0.035*obs 0.091 + 0.301*eptmean + 0.054*rain + 0.061*dq75 + 0.034*obs 0.089 + 0.302*eptmean + 0.054*rain + 0.058*dq75 + 0.034*obs 0.084 + 0.305*eptmean + 0.053*rain + 0.058*dq75 + 0.034*obs 0.086 + 0.304*eptmean + 0.054*rain + 0.058*dq75 + 0.034*obs 0.083 + 0.306*eptmean + 0.053*rain + 0.058*dq75 + 0.035*obs 0.081 + 0.310*eptmean + 0.053*rain + 0.060*dq75 + 0.034*obs 0.083 + 0.311*eptmean + 0.052*rain + 0.059*dq75 + 0.034*obs 0.085 + 0.312*eptmean + 0.052*rain + 0.057*dq75 + 0.031*obs 0.129 + 0.238*eptmean + 0.064*rain + 0.045*dq75 + 0.047*obs 0.129 + 0.239*eptmean + 0.064*rain + 0.045*dq75 + 0.046*obs 0.133 + 0.211*eptmean + 0.070*rain + 0.038*dq75 + 0.055*obs 0.173 + 0.322*eptmean + 0.063*rain + 0.052*obs + 0.173*tdmean + -0.021*gradvy850 0.167 + 0.303*eptmean + 0.069*rain + 0.053*obs + 0.164*tdmean + -0.026*gradvy850 0.167 + 0.302*eptmean + 0.070*rain + 0.053*obs + 0.163*tdmean + -0.027*gradvy850 0.151 + 0.235*eptmean + 0.081*rain + 0.051*obs + 0.121*tdmean + -0.031*gradvy850 0.151 + 0.238*eptmean + 0.079*rain + 0.051*obs + 0.121*tdmean + -0.034*gradvy850 0.154 + 0.244*eptmean + 0.079*rain + 0.053*obs + 0.125*tdmean + -0.035*gradvy850 0.154 + 0.246*eptmean + 0.078*rain + 0.051*obs + 0.126*tdmean + -0.036*gradvy850 0.154 + 0.238*eptmean + 0.080*rain + 0.051*obs + 0.120*tdmean + -0.035*gradvy850 0.157 + 0.242*eptmean + 0.078*rain + 0.053*obs + 0.122*tdmean + -0.034*gradvy850 0.157 + 0.242*eptmean + 0.078*rain + 0.051*obs + 0.121*tdmean + -0.034*gradvy850 xxi 2015121300 0.3013 2015122000 0.3026 2015122700 0.3032 0.155 + 0.235*eptmean + 0.079*rain + 0.052*obs + 0.116*tdmean + -0.033*gradvy850 0.156 + 0.239*eptmean + 0.079*rain + 0.052*obs + 0.119*tdmean + -0.033*gradvy850 0.156 + 0.239*eptmean + 0.079*rain + 0.052*obs + 0.119*tdmean + -0.033*gradvy850 xxii PHỤ LỤC Minh họa biến đổi phƣơng trình dự báo RV trạm Cà Mau hạn dự báo 24 với mơ hình WRF-ARW 2015053100 0.2824 2015060700 0.2803 2015061300 0.2772 2015062100 0.2798 2015062800 0.2791 2015070200 0.2797 2015071200 0.2822 2015071900 0.2811 2015072400 0.2743 2015080200 0.2731 2015080800 0.2762 2015081600 0.274 2015082300 0.2759 2015083000 0.2765 2015090600 0.2752 2015091300 0.2768 2015092000 0.2775 2015092700 0.2751 2015100400 0.3482 2015100800 0.3505 0.531 + 0.123*rhmean + 0.066*rain + 0.060*obs + 0.072*umean + 0.042*v500 0.531 + 0.121*rhmean + 0.066*rain + 0.062*obs + 0.072*umean + 0.042*v500 0.527 + 0.116*rhmean + 0.072*rain + 0.063*obs + 0.063*umean + -0.039*advttd700 0.531 + 0.118*rhmean + 0.070*rain + 0.064*obs + 0.064*umean + -0.040*advttd700 0.529 + 0.117*rhmean + 0.071*rain + 0.064*obs + 0.063*umean + -0.039*advttd700 0.529 + 0.117*rhmean + 0.071*rain + 0.062*obs + 0.065*umean + -0.040*advttd700 0.534 + 0.117*rhmean + 0.069*rain + 0.064*obs + 0.069*umean + -0.039*advttd700 0.531 + 0.116*rhmean + 0.070*rain + 0.066*obs + 0.068*umean + -0.039*advttd700 0.531 + 0.110*rhmean + 0.071*rain + 0.068*obs + 0.069*umean + -0.039*advttd700 0.529 + 0.107*rhmean + 0.071*rain + 0.069*obs + 0.072*umean + -0.037*advttd700 0.533 + 0.108*rhmean + 0.070*obs + 0.070*rain + 0.073*umean + -0.037*advttd700 0.533 + 0.108*rhmean + 0.071*obs + 0.069*rain + 0.074*umean + -0.038*advttd700 0.533 + 0.108*rhmean + 0.071*obs + 0.068*rain + 0.076*umean + -0.038*advttd700 0.530 + 0.107*rhmean + 0.074*obs + 0.069*rain + 0.076*umean + -0.037*advttd700 0.529 + 0.106*rhmean + 0.072*obs + 0.071*rain + 0.076*umean + -0.038*advttd700 0.409 + 0.156*qmean + 0.073*obs + 0.079*rain + 0.078*umean + -0.021*advttd700 0.410 + 0.157*qmean + 0.076*obs + 0.078*rain + 0.076*umean + -0.021*advttd700 0.412 + 0.298*qmean + 0.075*umean + 0.083*rain + 0.069*obs + -0.167*kindex 0.299 + 0.080*rain + 0.307*qmean + 0.064*u10m + 0.183*kindex + 0.072*obs 0.253 + 0.097*rain + 0.325*qmean + 0.075*u10m + - xxiii 2015101800 0.361 2015102500 0.3632 2015110100 0.369 2015110800 0.3763 2015111200 0.3819 2015112200 0.3788 2015112900 0.3746 2015120600 0.3648 2015121300 0.3702 2015122000 0.3697 2015122700 0.3591 0.183*kindex + 0.040*adveptmean 0.301 + 0.089*rain + 0.291*qmean + 0.064*u10m + 0.170*kindex + 0.069*obs 0.300 + 0.108*rain + 0.265*qmean + 0.069*u10m + 0.148*kindex + 0.057*obs 0.303 + 0.090*rain + 0.310*qmean + 0.050*u10m + 0.056*obs + -0.199*kindex + 0.053*q2m 0.305 + 0.092*rain + 0.315*qmean + 0.050*u10m + 0.203*kindex + 0.056*obs + 0.053*q2m 0.308 + 0.093*rain + 0.318*qmean + 0.048*u10m + 0.058*obs + -0.207*kindex + 0.055*q2m 0.305 + 0.097*rain + 0.314*qmean + 0.049*u10m + 0.202*kindex + 0.054*q2m + 0.054*obs 0.273 + 0.110*rain + 0.337*qmean + 0.055*u10m + 0.218*kindex + 0.062*q2m + -0.041*gradux500 0.272 + 0.111*rain + 0.338*qmean + 0.054*u10m + 0.215*kindex + 0.058*q2m 0.272 + 0.105*rain + 0.339*qmean + 0.056*u10m + 0.212*kindex + 0.061*q2m + 0.041*adveptmean 0.271 + 0.107*rain + 0.334*qmean + 0.054*u10m + 0.210*kindex + 0.062*q2m + 0.040*adveptmean 0.308 + 0.115*rain + 0.263*qmean + 0.064*u10m + 0.057*obs + -0.143*kindex xxiv PHỤ LỤC Minh họa biến đổi phƣơng trình dự báo RV trạm Cần Thơ hạn dự báo 24 với mơ hình WRF-ARW 2015053100 0.2702 2015060700 0.2729 2015061300 0.2722 2015062100 0.2761 2015062800 0.2715 2015070100 0.2725 2015071200 0.272 2015071900 0.2734 2015072400 0.2758 2015080200 0.2755 2015080900 0.2749 2015081600 0.2745 2015082300 0.2704 2015083000 0.2604 2015090600 0.2598 2015091300 0.2437 2015092000 0.2413 2015092700 0.2363 2015100400 0.327 2015101100 0.3385 0.473 + 0.203*rhmean + 0.065*rain + 0.057*dept75 + 0.046*advt700 0.471 + 0.203*rhmean + 0.067*rain + 0.059*dept75 + 0.047*advt700 0.471 + 0.202*rhmean + 0.068*rain + 0.058*dept75 + 0.047*advt700 0.473 + 0.204*rhmean + 0.067*rain + 0.059*dept75 + 0.048*advt700 0.474 + 0.205*rhmean + 0.064*rain + 0.058*dept75 + 0.046*advt700 0.472 + 0.205*rhmean + 0.064*rain + 0.058*dept75 + 0.047*advt700 0.477 + 0.206*rhmean + 0.062*rain + 0.055*dept75 + 0.047*advt700 0.476 + 0.206*rhmean + 0.063*rain + 0.055*dept75 + 0.050*advt700 0.483 + 0.209*rhmean + 0.059*rain + 0.061*dept75 + 0.052*advt700 + 0.042*v500 0.484 + 0.210*rhmean + 0.059*rain + 0.062*dept75 + 0.052*advt700 + 0.042*v500 0.482 + 0.209*rhmean + 0.060*rain + 0.061*dept75 + 0.053*advt700 + 0.042*v500 0.483 + 0.208*rhmean + 0.062*rain + 0.060*dept75 + 0.051*advt700 + 0.043*v500 0.495 + 0.221*rhmean + 0.059*rain + 0.064*dept75 + 0.050*advt700 + 0.045*vmean 0.478 + 0.200*rhmean + 0.067*rain + 0.054*dept75 + 0.049*advt700 0.476 + 0.200*rhmean + 0.068*rain + 0.054*dept75 + 0.049*advt700 0.406 + 0.157*rhmean + 0.084*rain + 0.056*u850 + -0.121*td700 0.334 + 0.101*rhmean + 0.084*rain + 0.054*u850 + 0.036*advt700 0.319 + 0.089*rhmean + 0.091*rain + 0.058*u850 + 0.030*advt700 0.194 + 0.081*rain + 0.237*eptmean + 0.037*u850 + 0.048*obs + -0.026*vmean + -0.118*td700 0.239 + 0.087*rain + 0.245*eptmean + -0.055*vmean + 0.052*obs + -0.020*ttd2m + -0.113*td700 xxv 2015101800 0.3504 2015102500 0.36 2015110100 0.3589 2015110800 0.3681 2015111200 0.3633 2015112200 0.3622 2015112900 0.3606 2015120600 0.3565 2015121300 0.3543 2015121900 0.3491 2015122700 0.3487 0.240 + 0.088*rain + 0.247*eptmean + -0.055*vmean + 0.051*obs + -0.021*ttd2m + -0.114*td700 0.213 + 0.100*rain + 0.228*eptmean + -0.058*vmean + 0.050*obs + 0.051*u850 + -0.099*td700 0.211 + 0.101*rain + 0.225*eptmean + -0.058*vmean + 0.051*u850 + 0.050*obs + -0.098*td700 0.213 + 0.105*rain + 0.228*eptmean + -0.056*vmean + 0.050*obs + 0.050*u850 + -0.101*td700 0.212 + 0.104*rain + 0.228*eptmean + -0.055*vmean + 0.052*u850 + 0.049*obs + -0.101*td700 0.211 + 0.105*rain + 0.222*eptmean + 0.051*obs + 0.054*vmean + 0.052*u850 + -0.098*td700 0.210 + 0.105*rain + 0.220*eptmean + 0.054*u850 + 0.055*vmean + 0.050*obs + -0.096*td700 0.209 + 0.106*rain + 0.216*eptmean + 0.055*u850 + 0.055*vmean + 0.049*obs + -0.093*td700 0.208 + 0.107*rain + 0.218*eptmean + -0.056*vmean + 0.053*u850 + 0.049*obs + -0.094*td700 0.239 + 0.111*rain + 0.127*eptmean + -0.075*vmean + 0.022*ttd2m + 0.052*obs + 0.034*advvormean 0.238 + 0.110*rain + 0.128*eptmean + -0.076*vmean + 0.022*ttd2m + 0.051*obs + 0.033*advvormean xxvi PHỤ LỤC 10 Minh họa biến đổi phƣơng trình dự báo RV trạm Hà Nội hạn dự báo 24 với mơ hình WRF-ARW 2015053100 0.2161 2015060700 0.2121 2015061300 0.1982 2015062100 0.211 2015062800 0.2127 2015070200 0.2168 2015071200 0.216 2015071900 0.2107 2015072400 0.2094 2015080200 0.2032 2015080900 0.2133 2015081600 0.2145 2015082300 0.2159 2015083000 0.2138 2015090600 0.2148 2015091300 2015092000 0.1872 0.1888 2015092700 0.1805 2015100400 0.2467 2015101100 0.2614 2015101800 0.264 0.491 + 0.162*qmean + -0.107*ttd2m + 0.063*obs + 0.054*wind700 + -0.046*gradttdx500 0.493 + 0.161*qmean + -0.101*ttd2m + 0.066*obs + 0.055*wind700 + -0.044*gradttdx500 0.575 + 0.293*qmean + -0.078*u10m + 0.069*obs + -0.125*td700 0.495 + 0.263*qmean + -0.102*ttd2m + 0.068*obs + 0.062*wind700 + -0.113*td700 0.493 + 0.265*qmean + -0.106*ttd2m + 0.067*obs + 0.061*wind700 + -0.115*td700 0.492 + 0.267*qmean + -0.107*ttd2m + 0.068*obs + 0.062*wind700 + -0.116*td700 0.490 + 0.162*qmean + -0.102*ttd2m + 0.068*obs + 0.053*wind700 + -0.047*gradttdx500 0.491 + 0.162*qmean + -0.102*ttd2m + 0.065*obs + 0.051*wind700 + -0.045*gradttdx500 0.492 + 0.163*qmean + -0.103*ttd2m + 0.062*obs + 0.051*wind700 + -0.045*gradttdx500 0.494 + 0.165*qmean + -0.105*ttd2m + 0.064*obs + 0.053*wind700 0.494 + 0.167*qmean + -0.101*ttd2m + 0.063*obs + 0.052*wind700 + -0.045*gradttdx500 0.570 + 0.300*qmean + -0.078*u10m + 0.067*obs + -0.127*td700 + -0.047*gradttdx500 0.569 + 0.304*qmean + -0.080*u10m + 0.068*obs + -0.131*td700 + -0.047*gradttdx500 0.568 + 0.306*qmean + -0.079*u10m + 0.066*obs + -0.133*td700 + -0.047*gradttdx500 0.571 + 0.312*qmean + -0.081*u10m + 0.064*obs + -0.137*td700 + -0.046*gradttdx500 0.573 + 0.191*qmean + -0.074*u10m + 0.098*obs + -0.063*td700 0.574 + 0.192*qmean + -0.075*u10m + 0.097*obs + -0.063*td700 0.542 + 0.123*qmean + -0.063*u10m + 0.103*obs + 0.064*wind700 0.450 + 0.098*obs + 0.079*q850 + 0.058*wind700 + 0.040*gradty850 + 0.075*qmean + -0.060*u10m 0.354 + 0.106*obs + 0.187*q850 + 0.069*wind850 + -0.092*t2m + -0.039*gradty850 + 0.023*div700 0.343 + 0.107*obs + 0.185*q850 + 0.017*wind700 + 0.036*gradty850 + 0.058*wind850 + -0.085*t2m + 0.025*div700 xxvii 2015102500 0.2776 2015110100 0.2794 2015110800 0.2726 2015111200 0.2732 2015112200 0.271 2015112900 0.2779 2015120600 0.2722 2015121300 0.28 2015122000 0.2711 2015122700 0.2694 0.346 + 0.118*obs + 0.140*q850 + 0.105*wind850 + -0.077*t2m + -0.061*gradty850 + -0.048*u10m 0.347 + 0.119*obs + 0.107*wind850 + 0.140*q850 + -0.079*t2m + -0.061*gradty850 + -0.047*u10m 0.349 + 0.119*obs + 0.139*q850 + 0.104*wind850 + -0.075*t2m + -0.059*gradty850 + -0.047*u10m 0.350 + 0.118*obs + 0.151*q850 + 0.101*wind850 + -0.082*t2m + -0.055*gradty850 + 0.045*div700 0.351 + 0.118*obs + 0.151*q850 + 0.099*wind850 + -0.082*t2m + -0.056*gradty850 + 0.045*div700 0.350 + 0.121*obs + 0.143*q850 + 0.094*wind850 + -0.076*t2m + -0.059*gradty850 + -0.046*u10m + 0.043*div700 0.350 + 0.121*obs + 0.139*q850 + 0.103*wind850 + -0.074*t2m + -0.060*gradty850 + -0.050*u10m 0.351 + 0.118*obs + 0.098*wind850 + 0.143*q850 + -0.079*t2m + -0.063*gradty850 + -0.046*u10m + 0.043*div700 0.349 + 0.118*obs + 0.105*wind850 + 0.139*q850 + -0.074*t2m + -0.062*gradty850 + -0.049*u10m 0.349 + 0.118*obs + 0.107*wind850 + 0.139*q850 + -0.075*t2m + -0.060*gradty850 + -0.048*u10m xxviii PHỤ LỤC 11 Minh họa biến đổi phƣơng trình dự báo RV trạm Vinh hạn dự báo 24 với mơ hình WRF-ARW 2015053100 0.2892 2015060700 0.2916 2015061300 0.2936 2015062100 0.2933 2015062800 0.3023 2015070200 0.3035 2015071200 0.2985 2015071900 0.2915 2015072400 0.2866 2015080200 0.3037 2015080900 0.3088 2015081600 0.3092 2015082300 0.307 2015083000 0.3078 2015090600 0.3074 2015091300 0.2974 2015092000 0.3031 2015092700 0.2893 2015100400 0.338 0.340 + 0.289*rhmean + -0.081*u850 + 0.052*wind500 + 0.072*rain + -0.142*tdmean + 0.053*obs 0.339 + 0.288*rhmean + -0.082*u850 + 0.053*wind500 + 0.071*rain + -0.141*tdmean + 0.054*obs 0.338 + 0.288*rhmean + -0.082*u850 + 0.053*wind500 + 0.071*rain + -0.142*tdmean + 0.054*obs 0.338 + 0.289*rhmean + -0.082*u850 + 0.052*wind500 + 0.071*rain + -0.141*tdmean + 0.053*obs 0.341 + 0.290*rhmean + -0.084*u850 + 0.046*wind500 + 0.057*obs + 0.063*rain + -0.137*tdmean + 0.039*advq500 0.342 + 0.289*rhmean + -0.083*u850 + 0.046*wind500 + 0.057*obs + 0.063*rain + -0.138*tdmean + 0.039*advq500 0.391 + 0.264*rhmean + -0.037*u850 + 0.053*wind500 + 0.079*rain + -0.130*tdmean + 0.050*obs + -0.039*ttd2m 0.336 + 0.289*rhmean + -0.075*u850 + 0.052*wind500 + 0.077*rain + -0.145*tdmean + 0.054*obs 0.333 + 0.270*rhmean + -0.079*u850 + 0.054*wind500 + 0.078*rain + 0.055*obs + -0.130*tdmean 0.393 + 0.272*rhmean + -0.036*u850 + 0.051*wind500 + 0.054*obs + 0.079*rain + -0.139*tdmean + -0.040*ttd2m 0.396 + 0.276*rhmean + -0.035*u850 + 0.050*wind500 + 0.056*obs + 0.081*rain + -0.142*tdmean + -0.040*ttd2m 0.395 + 0.279*rhmean + -0.036*u850 + 0.049*wind500 + 0.054*obs + 0.081*rain + -0.144*tdmean + -0.040*ttd2m 0.365 + 0.299*rhmean + -0.053*u850 + 0.049*wind500 + 0.060*obs + 0.076*rain + -0.187*tdmean + 0.070*rh850 0.364 + 0.303*rhmean + -0.053*u850 + 0.049*wind500 + 0.061*obs + 0.075*rain + -0.189*tdmean + 0.069*rh850 0.363 + 0.307*rhmean + -0.052*u850 + 0.049*wind500 + 0.058*obs + 0.074*rain + -0.193*tdmean + 0.070*rh850 0.449 + 0.104*obs + 0.080*rh850 + 0.076*rain + 0.018*wind500 + 0.172*rhmean + -0.090*tdmean + -0.045*ttd2m + 0.015*advq500 0.449 + 0.104*obs + 0.081*rh850 + 0.074*rain + 0.019*wind500 + 0.172*rhmean + -0.089*tdmean + -0.045*ttd2m + 0.016*advq500 0.436 + 0.118*obs + 0.078*rh850 + 0.069*rain + 0.015*wind500 + 0.133*rhmean + -0.057*tdmean + -0.034*ttd2m + -0.030*u850 0.456 + 0.144*obs + 0.129*rh850 + 0.094*rain + xxix 2015101100 0.3394 2015101800 0.3384 2015102500 0.3452 2015110100 0.3504 2015110800 0.3516 2015111200 0.3523 2015112200 0.3537 2015112900 0.3507 2015120600 0.3478 2015121300 0.3482 2015122000 0.3432 2015122700 0.3404 0.085*gradvy850 + 0.031*gradttdx850 + 0.044*dt87 0.466 + 0.146*obs + 0.123*rh850 + 0.090*rain + 0.092*gradvy850 + 0.042*gradttdx850 + 0.048*dt87 0.464 + 0.148*obs + 0.119*rh850 + 0.092*rain + 0.089*gradvy850 + 0.040*gradttdx850 + 0.046*dt87 0.483 + 0.154*obs + 0.111*rh850 + 0.090*rain + 0.097*gradvy850 + 0.053*gradttdx850 + 0.057*dt87 0.484 + 0.153*obs + 0.110*rh850 + 0.091*rain + 0.097*gradvy850 + 0.051*gradttdx850 + 0.057*dt87 0.485 + 0.155*obs + 0.110*rh850 + 0.091*rain + 0.096*gradvy850 + 0.051*gradttdx850 + 0.058*dt87 0.483 + 0.157*obs + 0.108*rh850 + 0.092*rain + 0.097*gradvy850 + 0.051*gradttdx850 + 0.057*dt87 0.483 + 0.158*obs + 0.108*rh850 + 0.092*rain + 0.096*gradvy850 + 0.051*gradttdx850 + 0.056*dt87 0.482 + 0.162*obs + 0.092*rh850 + 0.091*rain + 0.090*gradvy850 + 0.052*gradttdx850 + -0.052*ttd2m 0.481 + 0.161*obs + 0.092*rh850 + 0.093*rain + 0.088*gradvy850 + 0.052*gradttdx850 + -0.052*ttd2m 0.478 + 0.165*obs + 0.088*rh850 + 0.094*rain + 0.088*gradvy850 + 0.052*gradttdx850 + -0.050*ttd2m 0.474 + 0.169*obs + 0.087*rh850 + 0.092*rain + 0.083*gradvy850 + 0.051*gradttdx850 + -0.053*ttd2m 0.474 + 0.168*obs + 0.088*rh850 + 0.093*rain + 0.080*gradvy850 + 0.051*gradttdx850 + -0.052*ttd2m xxx ... Khu vực Đông Bắc 24 Khu vực Tây Bắc 21 Khu vực Việt Bắc 25 Khu vực Bắc Trung Bộ 20 Khu vực Trung Trung Bộ 15 Khu vực Nam Trung Bộ 12 Khu vực Tây Nguyên 17 Khu vực Nam Bộ 23 23 2.1.2 Số liệu dự báo. .. dự báo POP nhằm nâng cao chất lƣợng dự báo pha mƣa khu vực Việt Nam Trong phạm vi Luận văn áp dụng phƣơng pháp UMOS với sản phẩm dự báo mơ hình khu vực WRF-ARW Các phƣơng trình UMOS đƣợc xây dựng... Phƣơng pháp nghiên cứu 2.2.1 Phƣơng pháp UMOS 2.2.2 Phƣơng pháp GMOS 2.3 Áp dụng phƣơng pháp UMOS cho toán dự báo POP 2.3.1 Xác định yếu tố dự báo 2.3.2 Xác định nhân tố dự báo 2.3.3 Xây dựng phƣơng