(Luận văn thạc sĩ) các bài toán đa thức

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Mạnh Linh ĐÁNH GIÁ KHẢ NĂNG MÔ PHỎNG LƢỢNG MƢA VÀ NHIỆT ĐỘ CỦA CÁC MƠ HÌNH NHM VÀ REGCM TRÊN KHU VỰC VIỆT NAM LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội - 2014 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Mạnh Linh ĐÁNH GIÁ KHẢ NĂNG MÔ PHỎNG LƢỢNG MƢA VÀ NHIỆT ĐỘ CỦA CÁC MƠ HÌNH NHM VÀ REGCM TRÊN KHU VỰC VIỆT NAM Chuyên ngành: Khí tượng Khí hậu học Mã số: 60 44 0222 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS TSKH KIỀU THỊ XIN Hà Nội - 2014 Lời cảm ơn Tôi xin trân trọng bày tỏ lòng cảm ơn sâu sắc đến PGS TSKH Kiều Thị Xin, Trường Đại học Khoa học Tự nhiên, ĐHQGHN, người hướng dẫn khoa học, GS TS Nguyễn Hữu Dư chủ nhiệm đề tài khoa học “Nghiên cứu xây dựng kịch tượng thời tiết cực đoan trung hạn (2015-2030) cho khu vực Việt Nam – Biển Đơng sử dụng kịch biến đổi khí hậu chương trình Kakushin”, TS Lê Đức giúp đỡ nhiều mặt khoa học hợp tác quốc tế thời gian tơi hồn thành luận văn Tơi xin bày tỏ lịng cảm ơn sâu sắc tới Trung tâm Dự báo Khí tượng Thủy văn Trung ương, phòng Nghiên cứu & Ứng dụng tạo điều kiện thời gian giúp tơi hồn thành luận văn Tơi xin bày tỏ lịng cảm ơn tới thầy giáo, giáo, Khoa Khí tượng Thủy văn & Hải dương học, Khoa sau đại học Trường Đại học Khoa học Tự Nhiên giảng dạy, giúp đỡ cách tận tình chu đáo Trong suốt thời gian thực luận văn, nhận quan tâm giúp đỡ bạn bè, đồng nghiệp, ghi nhớ biết ơn giúp đỡ quý báu Cuối cùng, tơi xin bày tỏ lịng biết ơn sâu sắc tới người thân yêu gia đình tơi, ln nguồn động viên tinh thần q giá giúp tơi hồn thành luận văn Mục lục Mở Đầu _ CHƯƠNG 1: TỔNG QUAN VỀ NGHIÊN CỨU MƠ PHỎNG KHÍ HẬU KHU VỰC BẰNG MƠ HÌNH SỐ _ 1.1.Trên giới 1.2 Tại Việt Nam _ CHƯƠNG 2: SỐ LIỆU VÀ PHƯƠNG PHÁP ĐÁNH GIÁ _ 10 2.1 Số liệu 10 2.2 Mơ hình _ 15 2.2.1 Mơ hình NHRCM 15 2.2.2 Mơ hình RegCM _ 15 2.3 Thiết kế thí nghiệm _ 16 2.4 Phƣơng pháp đánh giá 17 CHƯƠNG 3: KẾT QUẢ ĐÁNH GIÁ _ 20 3.1 Đánh giá kết mô nhiệt độ _ 20 3.2 Đánh giá kết mô mƣa 29 KẾT LUẬN _ 40 Phụ lục _ 45 Danh sách chữ viết tắt AGCM AMeDAS AMV Aphrodite ATOVS CAM CCLM CCSM3 CLM ClWRF CMM5 COSMO COBE-SST CRCM CRIEPI DWD ECMWF ERA GCMs HE-VI ICTP JMA JRA25 MRI NCAR NCDC NCEP NHM NHRCM NICAM NorESM NPD OISST PRECIS TOVS RCMs RegCM REMO SST Atmospheric Global Circulation Model Automated Meteorological Data Acquisition System Atmospheric Motion Vector Asia Precipitation – Highly resolved Observational Data Intergration Towards Evaluations of the Waters Resources Advanced TIROS Operational Vertical Sounder Community Atmospheric Model Cosmo Climate Limited-area Model The Third - Community Climate System Model Community Climate Model Climate Weather Research Forecast Climate Mesoscale Model Consortium for small scale modeling Characteristics of Global Sea Surface Temperature Analysis Data The Third - Generation Canadian Regional Climate Model Central Research Institute of Electric Power Industry Deutscher Wetterdienst European Centre for Medium Range Weather Forecasts European Reanalysis Global Climate Models Horizontally Explicit – Vertically Implicit International Centre for Theoretical Physics Japan Meteorological Agency Japanese 25-year Reanalysis Meteorological Research Institute National Center for Atmospheric Research National Climatic Data Center National Center for Environmental Prediction Nonhydrostatic Regional Model Nonhydrostatic Regional Climate Model Nonhydrostatic ICosahedral Atmospheric Model Norwegian Earth System Model Numerical Prediction Division Optimum Interpolation Sea Surface Temperature Providing Regional Climate for Impacts Studies TIROS Operational Vertical Sounder Regional Climate Models Regional Climate Model Regional Climate Modeling Sea Surface Temperature Danh sách bảng STT Bảng Bảng 2.1 Bảng 2.2 Bảng 2.3 Nội dung Danh sách trạm quan trắc synop Thông tin số liệu tái phân tích Jra25 Các biến sử dụng số liệu Jra25 Trang 11 14 15 Danh sách hình vẽ biểu đồ STT hình Hình 2.1 Hình 2.2 Hình 2.3 Hình 3.1 Hình 3.2 Hình 3.3 Hình 3.4 Hình 3.5 Hình 3.6 Hình 3.7 Hình 3.8 Hình 3.9 Hình 3.10 Hình 3.11 Hình 3.12 Nội dung Phân bố mạng lưới trạm quan trắc (Vùng BB: hình vng xanh dương, BTB: hình vng xanh da trời, NTB: hình vng đen, NB: Hình vng xanh cây) Độ cao địa hình phân bố trạm quan trắc (các trạm quan trắc hình vng màu xanh da trời) Miền tính sử dụng nghiên cứu Kết mơ nhiệt độ khơng khí T2m (0C) trung bình năm thời kỳ 1986-2007 RegCM (a) NHRCM (b) Chênh lệch nhiệt độ khơng khí T2m (0C) trung bình năm thời kỳ 1986-2007 RegCM (a) NHRCM (b) so với số liệu tái phân tích Aphrodite Kết mơ nhiệt độ khơng khí T2m (0C) trung bình tháng mùa đơng DJF thời kỳ 1986-2007 RegCM (a), NHRCM (b) Aphrodite (c) Kết mơ nhiệt độ khơng khí T2m (0C) trung bình tháng mùa hè JJA thời kỳ 1986-2007 RegCM (a), NHRCM (b) Aphrodite (c) Hình 3.5: Biểu đồ tụ điểm nhiệt độ trung bình (0C) (a), tháng JJA (b), tháng DJF (c) thời kỳ 1986-2007 mơ hình RegCM (tam giác màu đỏ), NHRCM (hình vuông màu xanh) với giá trị quan trắc trạm Phân bố nhiệt độ T2m trung bình tháng (0C) giai đoạn (19862007) cho toàn Việt Nam (a), Bắc Bộ (b), Bắc Trung Bộ (c), Nam Trung Bộ (d) Nam Bộ (e) mơ hình RegCM (tam giác màu đỏ), NHRCM (hình vng màu xanh da trời) với giá trị quan trắc trạm (hình trịn xanh cây) Phân bố tần suất T2m trung bình ngày giai đoạn (1986-2007) cho toàn Việt Nam (a), Bắc Bộ (b), Bắc Trung Bộ (c), Nam Trung Bộ (d) Nam Bộ (e) mơ hình RegCM (tam giác màu đỏ), NHRCM (hình vng màu xanh da trời) với giá trị quan trắc trạm (hình trịn xanh cây) Tổng lượng mưa mơ trung bình năm (mm) thời kỳ 1986-2007 RegCM (a), NHRCM (b) Aphrodite (c) Tổng lượng mưa mơ trung bình tháng 6, tháng tháng (mm) thời kỳ 1986-2007 RegCM (a), NHRCM (b) Aphrodite (c) Tổng lượng mưa mô trung bình tháng 9, tháng 10 tháng 11 (mm) thời kỳ 1986-2007 RegCM (a), NHRCM (b) Aphrodite (c) Tỷ lệ mưa mô so với mưa phân tích Aphrodite mơ hình NHRCM (a) RegCM (b) Biểu đồ tụ điểm tổng lượng mưa trung bình năm (mm) (a), tháng JJA (b), tháng SON (c) thời kỳ 1986-2007 Trang 13 13 17 20 21 22 23 25 27 28 30 31 32 33 35 Hình 3.13 Hình 3.14 Hình 4.1 mơ hình RegCM (tam giác màu đỏ), NHRCM (hình vng màu xanh) với giá trị quan trắc trạm Phân bố lượng mưa tháng (mm) trung bình nhiều năm (19862007) cho toàn Việt Nam (a), Bắc Bộ (b), Bắc Trung Bộ (c), Nam Trung Bộ (d) Nam Bộ (e) mơ hình RegCM (tam giác màu đỏ), NHRCM (hình vng màu xanh da trời) với giá trị quan trắc trạm (hình trịn xanh cây) Phân bố tần suất tổng lượng mưa ngày giai đoạn (1986-2007) cho toàn Việt Nam (a), Bắc Bộ (b), Bắc Trung Bộ (c), Nam Trung Bộ (d) Nam Bộ (e) mơ hình RegCM (tam giác màu đỏ), NHRCM (hình vuông màu xanh da trời) với giá trị quan trắc trạm (hình trịn xanh cây) Sai phân theo thời gian thành phần bình lưu trường hợp ns = t /  =7 37 39 56 Mở Đầu Các mơ hình khí hậu tồn cầu mơ tốt q trình quy mơ lớn việc xem xét q trình quy mơ địa phương lại cần thiết Đáp ứng yêu cầu nhiều tổ chức giới phát triển ứng dụng mơ hình khí hậu khu vực với độ phân giải cao, chí vượt qua giới hạn thủy tĩnh (dưới 10km) - phi thủy tĩnh Nhưng để chạy mơ hình khí hậu phân giải cao địi hỏi nhiều chi phí tính tốn, thực tế nước có tài ngun máy tính lớn thực việc Tại Cộng Hòa Liên Bang Đức, CCLM (Cosmo Climate Limited-area Model) phiên dự báo khí hậu mơ hình thời tiết độ phân giải cao phi thủy tĩnh Cosmo, phát triển quan khí tượng quốc gia kết hợp với cộng đồng sử dụng mơ hình quy mơ nhỏ COSMO (Consortium for small scale modeling) sử dụng nghiên cứu khí hậu Bên cạnh đó, ICTP (International Centre for Theoretical Physics) phát triển mơ hình khí hậu khu vực RegCM, ứng dụng nhiều nơi giới Ở Nhật Bản, chạy thử nghiệm mơ hình NICAM (Nonhydrostatic ICosahedral Atmospheric Model) độ phân giải 4km, hay mơ hình NHRCM (Nonhydrostatic Regional Climate model) phiên khí hậu mơ hình dự báo thời tiết NHM (Nonhydrostatic Regional Model) đánh giá mô sát nhiệt độ lượng mưa giai đoạn năm từ 2001 đến 2006 (Hidetaka Sasaki nnk, 2011) với độ phân giải ngang km [17] Theo đó, xu hướng phát triển giới tiến tới chạy mơ hình khí hậu khu vực cao mơ hình phi thủy tĩnh Ở nước ta có nhiều mơ hình khí hậu khu vực sử dụng nghiên cứu ClWRF (phiên khí hậu mơ hình WRF), CMM5 (phiên khí hậu mơ hình MM5), REMO, … đặc biệt mơ hình RegCM sử dụng nhiều để mô yếu tố tượng khí hậu cực đoan (Phan Văn Tân nnk, 2008) [3, 4, 26] Bổ sung vào cơng trình nghiên cứu nước, đề xuất sử dụng phiên khí hậu mơ hình NHM (NHRCM - thuật ngữ sử dụng thay cho NHM phạm vi luận văn) để nghiên cứu khí hậu dựa kết đánh giá khả quan chuyên gia Nhật Bản [17] Nhưng trước đưa vào sử dụng cần phải kiểm chứng khả mơ khí hậu mơ hình NHRCM với mơ hình RegCM, mơ hình dùng nhiều Việt Nam, chuỗi số liệu đủ dài Bởi luận văn tiến hành “Đánh giá khả mô lượng mưa nhiệt độ mơ hình NHRCM RegCM khu vực Việt Nam” thực chuỗi số liệu 22 năm từ 1986 đến 2007 độ phân giải 20km chạy với điều kiện biên điều kiên ban đầu số liệu tái phân tích JRA25 (Japanese 25-year Reanalysis), số liệu tái phân tích nhiệt độ bề mặt nước biển SST Nhật Bản Số liệu lượng mưa nhiệt độ dùng để so sánh lấy từ phân tích Aphrodite(Asia Precipitation – Highly resolved Observational Data Intergration Towards Evaluations of the Waters Resources) độ phân giải 0.25 độ chuỗi quan trắc 58 trạm synop khu vực Việt Nam thời kỳ Luận văn bố cục gồm phần sau: Mở đầu Chƣơng 1: Tổng quan nghiên cứu mơ khí hậu khu vực mơ hình số Chƣơng 2: Số liệu phương pháp đánh giá Chƣơng 3: Kết đánh giá Kết luận Tài liệu tham khảo Trong đó: (38) Các sơ đồ bậc cao sử dụng tốt với sơ đồ sai phân bình lưu trung tâm biến đổi (Kato 1998).Sơ đồ kiểu hạn chế thông lượng hoạt động sơ đồ hiệu chỉnh thơng lượng Nó biến đổi số hạng bình lưu để khử bỏ sai số tính tốn sinh xấp xỉ sai phân hữu hạn để đảm bảo tính đơn điệu d Thiết lập sơ đồ HE-VI Sơ đồ tách thời gian hiển theo phương ngang ẩn theo phương đứng (HEVI) sử dụng phiên nghiệp vụ NHM Thiết lập sơ đồ HE-VI MRI/NPD-NHM tái áp dụng Muroi cộng (1999) có tham khảo phiên HE-VI Ikawa Saito (1991), mà có sóng âm thể ẩn mode cao tần dựa nghiên cứu Klemp Wilhel,son (1978) Các phương trình moment động lượng ngang (13) (14) lấy tích phân tiến theo thời gian m  P U    U  G1/ 2G13 P   1      ADVU  RU   m2  x G1/ z *  (39) m  P V    V  G1/ 2G 23 P         ADVV  RV   m1  y G1/ z *  (40) phương trình dộng lượng theo chiều đứng (15) áp dụng tích phân lùi theo thời gian sau đây: W   W P  g g   P  BUOY   ADVW  RW   1    P 1/  m3G z * m3Cm2 m3 m3Cm2 (41) đó: P         P  P 2 (42) Trong sơ đồ HE-VI MRI/NPD-NHM chọn Saito (2004) lại chọn , phiên , cho bảo toàn khối lượng đề cập đến phần 3f Cần ý, số hạng cuối vế trái phương trình (41) xuất độc 51 lập so với số hạng cuối vế phải đưa vào để đảm bảo tính ổn định xử lý ẩn sóng âm thành phần lực e Tắt dần phân kỳ Để tính tính bất ổn định giảm dần sóng âm sơ đồ tích phân theo thời gian HE-VI, lọc âm sử dụng NHM Mặc dù ý tưởng dựa vào nghiên cứu Skamarock Klemp (1992), thành phần giảm dần tính tổng thơng lượng phân kỳ (DIVT), cơng thức có khác chút so với Skamarock Klemp (1992) Các số hạng có gạch chân bổ sung vào số hạng cịn lại phương trình (13)-(15), áp dụng cho sơ đồ HE-VI Ở hệ số lấy ¼ giá trị cực đảm bảo điều kiện ổn định tuyến tính khuếch tán:  H  0.06 V  0.05  x   y        t  m1   m2   (46) G1/ z   t (47) f Phương trình khí áp Khí áp phương trình (32) tích phân lùi sau:     U  V   mm P   P       Cm2   PFT  m1m2    m3 *  1/2  W  +  G1/2G13 U +G1/2G 23V     -PRC      y  z  m3  x  G      DIF.P (48) DIF.P thành phần tắt dần Rayleigh khí áp cách   xác định U V tương tự phương trình (42) 52    Khử U V , W nhờ phương trình (39)-(42) cách tương tự mô tả (C-3-2) Saito cộng (2001a), ta nhận phương trình khí áp chiều dạng Helmholtz: 2    2 P   g   G1/ *     G1/    P  FP.HE.INV  FP.HE.VAR z *2 z  Cm  C   (1   )  m  (49) FP.HE.INV số hạng bất biến suốt bước thời gian tích phân ngắn số hạng cưỡng  G1/   DIF.P   PFT  PRC    1     Cm2  FP.HE.INV   G1/ 2   g  BUOY- m3  ADVW-RW  + P  *  z  Cm  (50) Và FP.HE.VAR số hạng biến đổi suốt bước thời gian tích phân ngắn số hạng cưỡng bức: FP.HE.VAR    G1/   m W  m1m DIVH  U  ,V y  + 1/32   1     G z *     2G1/  -  P  Cm  (1   )  (51) DIVH biểu diễn phân kỳ ngang DIVH U ,V   U V    *  G13U  G 23V  x y z Điều kiện biên thu với   g   1/2 * m G  z m 3Cm  (52) :    P    ADVW  RW   m   g   BUOY+ 1-  P  Cm   (53) Điều kiện biên giống phương trình (53) Trong trường hợp trượt tự (Free-slip) số hạng phụ bổ sung vào vế phải phương trình    m1m2  G1/ 2G13 U  +G1/ 2G 23V      W + m3   (54) Điều kiện biên (và dưới) (phương trình 53) cần giá trị khí áp biên biên bao gồm số hạng nhiễu động khí áp ( ) vế trái phương trình Đế tính số hạng cần giá trị khí áp kz =1 nz Điều kiện biên xác định kz =1+1/2 nz-1/2 mực W Phương trình khí áp Eliptic chiều (49) đưa lời giải cho khí 53 áp cho mực từ kz=2 tới nz-1 Vì vậy, giá trị khí áp biên ảnh hưởng tới cách giải bên phép ngoại suy đơn giản phá vỡ định luật bảo toàn khối lượng.Đế tránh lỗi này, chúng tơi tính tốn lực nhiễu động mật độ với lựa chọn động khí áp ( (phương trình 19), cho số hạng nhiễu ) vế trái phương trình (53) bù lại số hạng cuối ) việc xác định số hạng trở nên khó khăn vế phải ( g Tách thời gian sóng trọng trường Trong sơ đồ HE-VI (như mục 2.1.1.1d), có sóng âm thể ẩn mode cao tần.Tuy nhiên, mô hình nghiệp vụ, sóng trọng trường thường ảnh hưởng đến hiệu suất tính tốn mạnh Một sơ đồ tách theo thời gian cho sóng trọng trường Saito (2002) áp dụng Để tách sóng trọng trường đánh giá thành phần bình lưu thẳng đứng nhiệt độ vị trạng thái tham chiếu bước thời gian ngắn, phương trình nhiệt động lực học (9) viết sau:        N 2  '    Q    w w u v  DIF.   g z x y  c p  (55) Tuy nhiên, số hạng bình lưu NHM viết dang thơng lượng phương trình (31) thành phần hiệu chỉnh bình lưu áp dụng thực tế đề cập mục 2.2.1.1c Để sử dụng dạng sai phân hữu hạn bình lưu nhiệt độ vị NHM, phương trình (55) viết lại dạng:             d    Q    w w u v   DIF.   dz z x y  c p    d  Q   (w  w)  Adv.  + dz   c p  d  w  w   t dz (56) Ở số hạng cuối vế phải phương trình (56) cho tích phân thăm dò vi vật lý mây số hạng xác định sai phân lùi bậc Vì tích phân thăm dị thực dạng thơng lượng bậc cao 54 nên tính chất bảo tồn dạng thơng lượng bị giới hạn, trong bước thời gian ngắn phù hợp với w bước thời gian dài h Tách thời gian thành phần bình lưu Như đề cập trên, việc tách thời gian sóng trọng trường thực hầu hết trường hợp, nhiên, dẫn đến tình trạng phải đối mặt với tính bất ổn định tính tốn trường hợp có gió mạnh tồn lớp nghịch nhiệt mạnh Để khắc phục nhược điểm Saito (2003) phát triển sơ đồ tách thời gian cho thành phần bình lưu Hsu Sun (2001) đề xuất sơ đồ tách thời gian cho bình lưu mơ hình sơ đồ tiến họ với hy vọng số hạng bình lưu tính tốn bước thời gian ngắn  đảm bảo ổn định tính tốn Tuy nhiên, lựa chọn khơng chấp nhận sơ đồ bình lưu bậc cao với sơ đồ bình lưu biến đổi có chi phí tính tốn lớn làm ảnh hưởng đến hiệu tính tốn sơ đồ HE-VI Trong thực tế, ứng dụng đơn giản kỹ thuật tính thành phần bình lưu cho bước thời gian ngắn không thật hiệu sơ đồ Leapfrog áp dụng làm biến đổi sơ đồ Leapfrog từ sai phân trung tâm sang sai phân tiến theo thời gian Trong sơ đồ sai phân mới, đánh giá đầy đủ số hạng bình lưu bậc cao sơ đồ bình lưu biến đổi tâm bước thời gian Leapfrogvà sau điều chỉnh thành phần bậc thấp (bậc 2) bước thời gian ngắn nửa sau tích phân thời gian Leapfrog minh họa hình 4.1: ADV  ADV  kt   ADVL  kt   ADVL (57) Ở đây, ADV(kt) bình lưu bậc cao với sơ đồ bình lưu biến đổi bước thời gian kt, số hạng cịn lại ADVL(kt) ADVL(kt) thành phần bình lưu bậc thấp kt bước thời gian nhỏ  tương ứng Sự điều chỉnh thực từ (ns-1)/2+1 đến (ns-1), ns tỷ số Để xử lý ẩn sóng âm tách sóng trọng trường, phương trình moment động lượng thẳng đứng viết lại sau: 55 W   W P  g   P  BUOY    m3G1/ z * m3Cm2 m3   ADVW  ADVLM  ADVLW   RW   1    (58) g P mCm2 Hình 4.1: Sai phân theo thời gian thành phần bình lưu trường hợp ns = t /  =7 Việc tách thời gian bình lưu nhiệt độ thếlà cách lựa chọn để tách sóng trọng trường:       Q    Adv.  ADVL  ADVL     DIF.  c p (59)     ADVL  ADVL      t  Ở tính cơng thức sai phân trung tâm bậc dạng thông lượng Trong phiên NHM nghiệp vụ phương trình (59) sử dụng thay cho phương trình (56) để đơn giản tiết kiệm cho tính tốn với độ phân giải 10km Như đề cập phần trước, số hạng cuối vế phải phương trình (59) có tích phân thời gian lựa chọn diễn biến vi vật lý mây i Điều kiện biên xung quanh điều chỉnh thông lượng Ta áp dụng lơpứ đệm xốp tắt dần Rayleigh vào vùng gần biên xung quanh để buộc biến dự báo biến đổi Trong MRI/NPD-NHM lớp tắt dần Rayleigh áp dụng cho moment động lượng, khí áp, nhiệt độ vị, mơ hình gốc NHM tắt dần Rayleigh chỉsử dụng cho độ ẩm riêng Để đáp ứng mối quan hệ xu diễn biến thời gian tổng khối lượng phạm vi tồn miền tính thơng lượng trao đổi qua biên xung quanh áp dụng sơ đồ điều chỉnh thông lượng qua biên xung quanh Lấy tích 56 phân khối phương trình liên tục (29) thu quan hệ sau xu thời gian tổng khối lượng tồn miền mơ hình thong lượng khối qua biên xung quanh:  t    G1/ m2 dxdydz *    Vdxdz  * x Y   Udydz    m1m2 * x 0    Udydz  *  PRC - q w  ' v ' z*  x X    Vdxdz  * y 0 (60) dxdy Số hạng cuối phương trình hiệu số tốc độ mưa trung bình khu vực tổng thơng lượng bề mặt nước (tổng bốc bề mặt) Vì số hạng thường bù trừ (offset) cho nên giá tri hiệu chỉnh thông lượng khối biên xung quanh giảm nhờ bổ sung số hạng khuếch tán ẩm vào phương trình liên tục 2.2.1.2 Các trình vật lý a Vật lý mây Ikawa Saito (1911) áp dụng sơ đồ vi vật lý mây ba pha cho NHM dựa so đồ Lin cộng (1983) Thành phần nhiễu động khí áp tuyến tính hóa để tính nhiễu động nhiệt độ vị MRI/NPD-NHM, gần bị loại bỏ phiên NHM Mặc dù vi vật lý mây đầy đủ dự báo nồng độ ngưng kết mây, tuyết mưa đá, có tỷ số xáo trộn dự báo thức phiên nghiệp vụ Ngoài ra, phiên tiết kiệm phương pháp bulk Yamada (2003a) phát triển, số số hạng khơng quan trọng (chủ yếu cho tạo thành mưa đá) bị loại bỏ Một điều chỉnh bổ sung nghiệp vụ cho phép tốc độ bốc mưa, tuyết, mưa đá thay đổi tới 50% so với tốc độ ban đầu để giảm bốc nhiều tầng thấp khí chạy mơ hình với độ phân giải 10km Bên cạnh giá trị ngưỡng chuyển đổi tự động Kessler từ nước mây sang mưa thay đổi từ đến 0.1g/Kg để tránh xác định thiên thấp lượng mưa nhỏ Sơ đồ hộp Lagarangian (Kaito 1995) không áp dụng tính đến hạt rơi khỏi (fallout of rain) mà mưa đá để đảm báo tính ổn định tính toán nghiệp vụ với = 40s 57 Để chuẩn bị điều kiện ban đầu cho sơ đồ vi vật lý mây sử dụng hệ thống dự báo chu kỳ 6h, biến vi vật lý mây lấy từ dự báo h NHM, nút lưới có độ ẩm tương đối nhỏ 90% giá trị biến gán (Ishida Saito 2005) b Tham số hóa đối lưu Tham số hóa mây tích vấn đề quan trọng để xác định hiệu mơ hình mơ hình hoanf lưu chung khí (GCM) mơ hình dự báo số trị (NWP) Mặc dù MRI/NPD-NHM có sơ đồ tham số hóa đối lưu (sơ đồ Arakawa – Schubert sơ đồ đối lưu ẩm có hiệu chỉnh) cách thiết lập sơ đồ khơng thích hợp chạy mơ hình với độ phân giải ngang cao (10km), sơ đồ ban đầu phát triển cho GCM với giả thiết bước lưới ngang lớn nhiều so với quy mơ mây Trong mơ hình NHM, sơ đồ tham số hóa đối lưu Kain-Fritsch (K-F) (Kain and Fritsch 1993) Yamada áp dụng (2003b) Ngoài thay đổi bổ sung gần mô tả Kain (2004), số sửa đổi Ohmori Yamada (2004) để cải thiện hiệu dự báo mưa lớn với độ phân giải 10km Nhật, nơi mà ẩm khối khí biển khơng ổn định chiếm ưu Trong hàm “khởi động” (“trigger function”) sơ đồ K-F, đưa vào ý tưởng gia tăng nhiệt độ để mô nhiễu động cưỡng địa phương hàm dịng thăng quy mơ lưới Dù sơ đồ gốc gia tăng nhiệt độ coi hàm độ phân giải ngang sơ đồ cho lượng mưa địa hình cao Như vậy, yếu tố gia tăng nhiệt độ khối khí mực ngưng kết (LCL) giảm đi, từ giá trị ban đầu 1K xuống 0.2K.Để giảm bớt xác định thiên cao chuyển động quy mơ lưới, dịng thăng điểm nút mơ hình tính giá trị trung bình điểm xung quanh Vì sơ đồ xác định mưa gốc mà Ogura Cho (1973) sử dụng để tính “đặc điểm dịng thăng”, giáng thủy tạo sinh khu vực có dịng thăng mà khơng tính đến lượng ngưng tụ Với sơ đồ này, kiểu mưa không thật 58 nhạy nhận giống vùng mưa kéo dài có hướng vng góc với đới mưa chính.Các kiểu mưa tương tự mơ đới mưa với áp dụng sơ đồ K-F mô tả cơng trình Knievel (2002) H Kusaka (2004).Để cải thiện dạng mưa vậy, áp dụng sơ đồ dựa ý tưởng Kessler-type chuyển đổi tự động, chất ngưng kết dòng thăng chuyển đổi thành mưa giá trị chúng vượt ngưỡng cho trươc(1.0 g/Kg) Liên quan đến “giả thiết khép kín”, sơ đồ K-F mặc định rằng, đối lưu dùng gần hết lượng CAPE (90%) có trước sơ đồ K-F Tuy nhiên, thử nghiệm dự báo mặc định làm khí ổn định lượng mưa lớn giảm theo thời gian đến hạn dự báo 18h Để ngăn ngừa q bất ổn định khơng muốn khí mơ hình người ta thu nhỏ tốc độ giảm CAPE cột khí đơi giảm cách giảm giá trị mặc định sơ đồ K-F từ 90% xuống 85% c Khuếch tán ẩm Với thay đổi sơ đồ K-F trình bày trên, việc mô mưa lớn cải thiện.Tuy nhiên, cường độ dịng thăng khối khí quy mơ lưới mơ hình thường phát triển với tốc độ tăng phần dư CAPE.Sư thay đổi gây nhiều bất ổn định đối lưu khí mơ hình Để kiểm tra được xốy nút lưới mưa quy mô lưới lớn kèm theo, biện pháp gọi khuyếch tán ẩm kết cục (Targered moisture diffusion) TMD áp dụng NHM (Saito Ishida 2005) Ý tưởng TMD phát triển T.Davies (Met Office) cho mơ hình UM (Davies et al 2005), đó, khuếch tán ngang bậc áp dụng cho nước tồn dịng thăng mạnh để kìm hãm chọn lọc xoáy nút lưới sinh hiệu ứng dương ngược giải phóng ẩn nhiệt ngưng kết gia tốc thăng Trong mơ hình UM quy mơ vừa, 38 mực độ phân giải ngang 0.110 (12-km), TMD áp dụng cho tất điểm lưới cột, nơi có chuyển động dịng thăng vượt q 59 m/s Trong trường hợp mơ hình NHM, khuếch tán ngang hàm e mũ 300s áp dụng cho bốc điểm lưới, nơi có dịng thăng vượt q m/s d Các q trình bề mặt Mơ hình MRI/NPD-NHM có số vấn đề thuộc tính bề mặt: điều kiện không ổn định tồn độ lệch hệ thống (bias) lạnh ẩm Xác định cao dòng ẩn nhiệt bề mặt nguyên nhân gây hai độ lệch Dưới thay đổi thực Kumagai (2004a): Giá trị thông số nhám (hệ số nhám) vô hướng nhiệt ẩm ( zoh ) giảm khoảng 1/7.4 so với hệ số nhám động lượng ( zom ) theo Garatt Francey (1978) Sự thay đổi làm giảm dòng ẩn nhiệt bề mặt Tính tốn hệ số khối (bulk coefficients) cho thông lượng bề mặt đất liền thay đổi từ Sommeria (1976) sang Louis cộng (1982) MSM Sự thay đổi làm giảm xác định cao thông lượng bề mặt trường hợp ổn định mạnh Các hệ số khối bề mặt (bulk coefficients) thong lượng động lượng ( Cdm Cdh ) nhiệt ( ) tính công thức:      am2  2bRiB       dRi B       2 2bRiB  am   za    3bam2 c   z 0m      : RiB     RiB    : Ri  B        am ah   : RiB       3bRiB  dRiB          3bRiB   : Ri  B am ah 1   z a  RiB    3bam ah c  z0    60 (61) (62) số Richardson xác định bởi: RiB  am  gza  a   s  s ua2  va2  z log a z0 m (63)  , ah  log za z0 h (64) Ở b = c = d =  (=0.4) số Kármán, số bên a s biểu giá trị mực thấp mơ hình ( ) bề mặt Theo MSM, biến kháng trở khí khổng (stomatal resistance) đưa vào để biểu diễn biến thiên ngày biến thiên mùa bốc bề mặt Thông lượng nước bề mặt tham số hóa dịng giáng thơng lượng xạ sóng ngắn xuống qv' w'   Cdq ua2  va2 rs  rs ,day  rs , night S 1  S0  rs S :  Qva  Qvs  (65) (66) kháng trở khí khổng W /m2 Vào ban ngày, đạt đến (30s/m vào tháng đến tháng 10 60s/m vào tháng lại) Vào buổi tối bắt đầu + (300m/s) Phép tham số hóa làm giảm thơng lượng ẩn nhiệt làm hạn chế q trình ẩm mực thấp e Các trình lớp biên Độ lệch ẩm bề mặt trình bày gây khuếch tán yếu lớp biên hành tinh (PBL) NHM Để xem xét hiệu ứng phi địa phương, số thay đổi thực Kumagai(2004b) Độ dài xáo trộn phiên gốc NHM xác định bước lưới điển hình độ ổn định khí địa phương công thức Deardorff‟s (1980) thay công thức gần mặt Blackadar (1962) cho độ dài xáo trộn gần tích số Kármán độ cao mặt đáy Tuy 61 nhiên, thay đổi dẫn đến đánh giá thiên thấp độ dài xáo trộngần bề mặt ngăn cản sử phát triển lớp biên đối lưu Trong mơ hình NHM theo Sun Chang (1986), độ dài xáo trộn thẳng đứng phương trình TKE xác định cơng thức đây, có sử dụng đỉnh PBL (  lSC  z   z   1.2    hPBL  max lSC  hPBL  0.2      z   1.2  lz    hPBL  max lSC  hPBL  0.2      s    E1/   s.0.76  N   ): : z  hPBL   s    : hPBL  z  1.2hPBL , N   E1/  s.0.76 N         : hPBL  z  1.2hPBL , N    : z  1.2hPBL , N  : z  1.2hPBL , N  (67) Ở đây, E động rối (TKE), N tần số Brunt-Vaisala      z  z    lSC  z   0.25 1.8hPBL 1  exp  4   0.0003exp    hPBL    hPBL        (68) Đỉnh PBL xác định mực mà nhiệt độ vị ảo vượt giá trị mực thấp Tính vơ hướng rối xem xét, độ dài xáo trộn thẳng đứng xác định riêng theo chiều thẳng đứng chiều ngang với sử dụng bước lưới thẳng đứng ngangtương ứng Số rối Prandtl hệ số rối khuếch tán tính riêng biệt theo hướng thẳng đứng ngang với sử dụng độ dài xáo trộn đứng ngang tương ứng f Cảnh báo lượng rối Mơ hình rối khép kín bậc 1.5 đưa vào mơ hình NHM để xác định hệ sốkhuếch tán Sơ đồ xem TKE biến dự báo tích phân phương trìnhTKE Trong sơ đồ dự báo TKE, điều kiện ban đầu biên xung quanh TKE địi hỏi trì cấu trúc lớp biên Tuy nhiên, điều kiện không 4DVAR Meso RSM cung cấp.Ngoải ra, thời gian tích phân TKE tốn sử dụng sơ đồ bậc cao sơ đồ bình lưu sửa đổi Trong mơ hình NHM, gió bề mặt mạnh (mode) lỗi tính tốn xuất trường hợp chạy 62 với độ phân giải ngang 10km At=40s nếy TKE tính với sử dụng sơ đồ sai phân trung tâm bậc đới giản không dung sơ đồ hiệu chỉnh bình lưu Để khắc phục yếu điểm nêu trên, lựa chọn thực Kumagai Saito (2004), TKE không dự báo mà cảnh báo Giả sử có cân địa phương công thức TKE ta nhận tương quan lực nổi, số hạng tích độ đứt số hạng tiêu hao ma sát nhớt:   u u j gK hz v   K mj  i  *   z i j  x j xi Ở đây, ,  ui Ce  K mz      0 lz  Cmlz   x j (69) hệ số khuếch tán thẳng đứng cho moment động lượng nhiệt tương ứng, sử dụng độ dài xáo trộn thẳng đứng số nghịch đảo rối Prandtl, có cơng thức sau: K mz  Cmlz E1/ (70) K hz  Prz1 Kmz (71) Prz1  1 lz z (72) thu công thức (69) (71): K mz  lj Cm3  lz   Ce  i , j lz  ui u j   xi  x j  ui g 1 v   Prz   z *   x j  (73) lượng rối cảnh báo công thức (70) g Xử lý ẩn khuếch tán thẳng đứng Trong mơ hình MRI/NPD-NHM, khuyếch tán thẳng đứng xử lý hiển, giới hạn hệ số khuếch tán công thức sau: K  K max( z )  0.1  z  t (74) để đáp ứng điều kiện ổn định tuyến tính Khi chọn giá trị lớn, điều kiện (74) làm hạn chế khuyếch tán thẳng đứng gần bề mặt Đơi điều gây sai lệch (bias) ẩm lạnh mơ hình MRI/NPD-NHM Xử lí ẩn (implicit treatment) khuếch tán thẳng đứng ứng dụng đãđược Fujibe 63 cộng (1999) đưa vào MRI-NHM áp dụng lại cho NHM Phương trình dự báo cho  t t   t t 2t   Adv. t   H  K H  H  t t   theo thời gian sau:    t t   Kz  z  z  (75) thời điểm thu cách giải phương trình eliptic sau: Giá trị 2t rời rạc hóa    t t  t t K  2t   Adv. t   H  K H  H  t t     t t  Z   z  z  (76) Ở thử nghiệm đầu tiên, thông lượng bề mặt xử lí thành phần cưỡng điều kiện biên Tuy nhiên, gió mạnh gắn liền với bão dự báo bên sườn dốc núi, xử lí gây bất ổn định tính tốn Trong phiên NHM nghiệp vụ, thơng lượng bề mặt thời điểm tính tốn theo xử lý ẩn hồn tồn, có sử dụng hệ số khối h Bức xạ khí Bức xạ khí thiết lập Sugi cộng (1990) đưa vào MRI/NPD-NHM trình bày chương G-5 Saito et al (2001a) Hàm chuyển đổi xạ sóng dài đượ tính mơ hình ngẫu nhiên Goody (1952), tính đến hấp thụ nước CO2 Bức xạ sóng ngắn tính phương pháp gân hai chiều tham số hóa chiều dài sóng nhỏ 0.9 lớn 0.9 Trong mơ hình phần tử hấp thụ tính đến nước giọt mây Lacis Hansen (1974) đưa giải tham số Phần mây tham số hóa, sử dụng hàm cho độ ẩm tương đối Phần mây tham số hóa hàm cho trước độ ẩm tương đối, mây phân làm ba loại Smagorinsky (1960)và giá trị tới hạn độ ẩm để tính tốn phần mây xác định với tham khảo thống kê Ohno Isa (1984) Trong phiên nghiệp vụ NHM, giới hạn để cảnh báo mây tầng cao thay đổi từ 300hPa MRI/NPD-NHM tới 100hPa đỉnh tầng đối lưu để tránh làm lạnh lựa chọn không thật gần mực 300hPa.Độ dày quang học xạ sóng ngắn mây tầng tham số hóa cho mực áp suất 64 để xem xét thuộc tính xạ Giới hạn thấp để cảnh báo mây tầng thấp thay đổi từ 230m MRI/NPD-NHM tới 500m để cải thiện làm lạnh không thật mực thấp khí biển vào mùa đơng Việc tính tốn q trình xạ thực 15 phút với toán tử làm trơn không gian ngang phân giải 20 km Quá trình làm trơn khác với MSM Ở q trình tính xạ lặp lại tùng 60 phút/lần khơng có làm trơn theo khơng gian 65 ... Dương, Phan Văn Tân, Võ Văn Hịa, 2009: “Đánh giá khả mơ số yếu tố khí hậu cực đoan mơ hình khí hậu khu vực RegCM”, Tạp chí Khí tượng Thủy văn, Hà Nội, (584), tr 15-23 Bùi Hoàng Hải, Phan Văn Tân,... hồn thành luận văn Tơi xin bày tỏ lịng cảm ơn sâu sắc tới Trung tâm Dự báo Khí tượng Thủy văn Trung ương, phòng Nghiên cứu & Ứng dụng tạo điều kiện thời gian giúp tơi hồn thành luận văn Tơi xin... giáo, Khoa Khí tượng Thủy văn & Hải dương học, Khoa sau đại học Trường Đại học Khoa học Tự Nhiên giảng dạy, giúp đỡ cách tận tình chu đáo Trong suốt thời gian thực luận văn, nhận quan tâm giúp