Dự báo đông bắc bộ

T¹p chÝ Khoa häc ®hqghn, KHTN & CN, T.xxII, Sè 1PT., 2006 Dù b¸o d«ng ë b¾c bé TrÇn T©n TiÕn, NguyÔn ThÕ Vinh, §Æng ThÞ An Khoa KhÝ t−îng-Thuû v¨n vµ H¶i d−¬ng häc Tr−êng §¹i häc Khoa häc Tù nhiªn, §HQG Hµ Néi 1. Më ®Çu Mäi ho¹t ®éng trong ®êi sèng x· héi cña con ng−êi ®Òu trùc tiÕp chÞu ¶nh h−ëng cña c¸c ®iÒu kiÖn thêi tiÕt, nhÊt lµ c¸c hiÖn t−îng thêi tiÕt nguy hiÓm nh− b·o, m−a lò, lèc, tè vµ d«ng. Nh÷ng hiÖn t−îng thêi tiÕt nguy hiÓm trªn ®Æc biÖt lµ d«ng ¶nh h−ëng lín ®Õn hiÖu qu¶ cña c¸c ngµnh kinh tÕ nh− hµng kh«ng, hµng h¶i, du lÞch, b−u chÝnh viÔn th«ng. V× vËy viÖc dù b¸o d«ng lµ mét vÊn ®Ò cÊp b¸ch ®Æt ra hiÖn nay. Cã rÊt nhiÒu ph−¬ng ph¸p ®Ó dù b¸o d«ng song theo −¬c tÝnh cho thÊy dù b¸o b»ng m« h×nh chØ ®óng ®−îc kho¶ng 10% cßn l¹i do kinh nghiÖm ®¹t 30% vµ theo b¶n n¨ng lµ 60%. C¸c ph−¬ng ph¸p dù b¸o d«ng trªn thÕ giíi chñ yÕu dùa vµo c¸c chØ sè ®Æc tr−ng cho ®é bÊt æn ®Þnh cña khÝ quyÓn. C¸c chØ sè th−êng dïng lµ LI ( Liftid index), CAPE (Convective Available Potention Eneregy ), ®é ®øt cña giã, giã trªn c¸c mÆt ®¼ng ¸p, nhiÖt ®é, ®é Èm v.v…Trong bµi b¸o nµy chóng t«i sö dông kÕt qu¶ dù b¶o tr−êng ®Þa thÕ vÞ , nhiÖt ®é, ®iÓm s−¬ng cña m« h×nh bÊt thñy tÜnh ETA ®Ó tÝnh to¸n c¸c chØ sè bÊt æn ®Þnh vµ x©y dùng ph−¬ng tr×nh dù b¸o d«ng cho tõng tr¹m thuéc B¾c bé. 2. M« h×nh ETA M« h×nh khu vùc h¹n chÕ ETA do tr−êng ®¹i häc Belgrade vµ ViÖn KhÝ t−îng Thñy v¨n Federal – Belgrade cïng víi Trung t©m KhÝ t−îng Quèc gia Mü x©y dùng. M« h×nh liªn tôc ®−îc thay ®æi th«ng qua c¸c héi th¶o hµng n¨m cña c¸c chuyªn gia vÒ m« h×nh trªn thÕ giíi ®ang sö dông ETA. M« h×nh ETA hiÖn nay ®−îc NCEP c¶i tiÕn trë thµnh mét trong c¸c m« h×nh sè trÞ ch¹y nghiÖp vô dù b¸o ng¾n h¹n t¹i Mü. Trªn thÕ giíi cã nhiÒu n−íc sö dông nh− Nam T−, Hy L¹p, Italy, Nam Mü… Phiªn b¶n míi nhÊt cña m« h×nh ETA kh«ng thuû tÜnh ®−îc hoµn thiÖn vµ sö dông t¹i §¹i häc tæng hîp Hy L¹p ®øng ®Çu lµ gi¸o s− G.Kaloss. HiÖn nay nhãm nghiªn cøu dù b¸o khÝ t−îng tr−êng §¹i häc Khoa häc tù nhiªn - §HQGHN víi sù gióp ®ì cña gi¸o s− G. Kaloss ®· ¸p dông thµnh c«ng m« h×nh nµy vµo dù b¸o thêi tiÕt ë ViÖt Nam. 2.1. HÖ ph−¬ng tr×nh thñy nhiÖt ®éng lùc häc cña m« h×nh ETA Tr−íc hÕt ta xÐt hÖ täa ®é th¼ng ®øng sigma (σ) µ ππσ )( t−= (1) trong ®ã, π lµ ¸p suÊt thñy tÜnh, µ lµ hiÖu sè gi÷a ¸p suÊt thñy tÜnh bÒ mÆt vµ ®Ønh m« h×nh, tøc lµ: ts ππµ −= . (2) 56 Dù b¸o d«ng ë B¾c Bé 57 ë ®©y, πs vµ πt lµ ¸p suÊt thñy tÜnh t¹i bÒ mÆt vµ ®Ønh cña m« h×nh. Víi hÖ täa ®é nµy th× hÖ ph−¬ng tr×nh nhiÖt ®éng lùc bÊt thñy tÜnh cã d¹ng sau (Janjic et al., 2001) [27]: 1) Ph−¬ng tr×nh cho xu thÕ khÝ ¸p: σµµ σ ′∫∇−=∂ ∂ • d t )( 1 0 v (3) 2) Ph−¬ng tr×nh tr¹ng th¸i: RTp =α (4) 3) Ph−¬ng tr×nh cña ®é cao ®Þa thÕ vÞ: ∫+= 1 σ s dσp RTµΦΦ (5) 4) Ph−¬ng tr×nh cña c¸c thµnh phÇn giã ngang: vkv ×+∇−Φ∇+−= rfp dt d σσ αε )1( (6) 5) Ph−¬ng tr×nh cña xu thÕ nhiÖt ®é (b¶o toµn n¨ng l−îng): pp p c Q tt p c dp c TT= t T +∂ ∂+−∂ ∂+ ′∫ ∇+−∇+∂ ∂−∇−∂ ∂ ••• ])1([ ])()1([ 0 πεα σµεασσ σ σσσ vvv & (7) επ +=∂ ∂ 1p (8) 6) Ph−¬ng tr×nh cña tèc ®é th¼ng ®øng: )(11 σ σ tgdt d g w σ ∂ ∂+∇+∂ ∂== • ΦΦΦΦ &v (9) )(11 σ wσw t w gdt dw g σ ∂ ∂+∇+∂ ∂== • &vε (10) 7) Ph−¬ng tr×nh ®èi víi tû sè hçn hîp q, ph−¬ng tr×nh dù b¸o cã d¹ng: S dt dq = (11) ë ®©y, S lµ nguån Èm, v lµ vector giã ngang, p lµ ¸p suÊt kh«ng thñy tÜnh, R lµ h»ng sè khÝ cña kh«ng khÝ kh«, T lµ nhiÖt ®é, Φ lµ ®Þa thÕ vÞ vµ Φs lµ ®Þa thÕ vÞ t¹i bÒ mÆt ®Êt. HÖ ph−¬ng tr×nh (1)-(11) lµ hÖ ph−¬ng tr×nh thuû nhiÖt ®éng lùc häc cho khÝ quyÓn kh«ng thuû tÜnh. Tr−êng hîp ε = 0, m« h×nh trë vÒ tr−êng hîp thuû tÜnh. Trong khÝ quyÓn rèi, th× c¸c ph−¬ng tr×nh chuyÓn ®éng, nhËp nhiÖt, Èm cã thªm c¸c thµnh phÇn m« t¶ x¸o trén rèi th¼ng ®øng. TrÇn T©n TiÕn, NguyÔn ThÕ Vinh, §Æng ThÞ An 58 Theo ®¸nh gi¸ cña NCEP th× m« h×nh ETA cho kÕt qu¶ dù b¸o kh¸ tèt ë miÒn nhiÖt ®íi, ®Æc biÖt lµ khu vùc ®åi nói. Víi to¹ ®é ®Þa h×nh η, m« h×nh cã thÓ n¾m b¾t ®−îc dßng chuyÓn ®éng c−ìng bøc ®Þa h×nh t¹o ®iÒu kiÖn thuËn lîi cho ®èi l−u ph¸t triÓn. V× vËy viÖc sö dông kÕt qu¶ cña m« h×nh bÊt thuû tÜnh ETA tÝnh to¸n c¸c chØ sè bÊt æn ®Þnh ®Ó dù b¸o d«ng lµ rÊt kh¶ quan. 3. TÝnh to¸n c¸c chØ sè ®Ó dù b¸o d«ng theo m« h×nh ETA Chóng t«i tÝnh to¸n 3 chØ sè ®Ó dù b¸o d«ng theo m« h×nh ETA lµ: chØ sè Boyd, chØ sè TT, chØ sè KI. 3.1. ChØ sè Boyd: C«ng thøc to¸n häc cña chØ sè Boyd: ( ) 2001.0 7001000700 −−−= TZZBoyd , trong ®ã: Z700 lµ ®é cao ®Þa thÕ vÞ ë mùc 700mb (m2/s2) Z1000 lµ ®é cao ®Þa thÕ vÞ ë mùc 1000mb (m2/s2) T700 lµ nhiÖt ®é t¹i mùc 700mb (K0) ChØ sè Boyd m« t¶ profile nhiÖt ®é th¼ng ®øng gi÷a mùc 1000mb vµ 700mb. 3.2. ChØ sè KI: C«ng thøc to¸n häc cña chØ sè KI: ( ) TTTT ddKI 700850500850 ++−= , trong ®ã: T850 lµ nhiÖt ®é t¹i mùc 850mb (K0) T500 lµ nhiÖt ®é t¹i mùc 500mb (K0) Td850 lµ nhiÖt ®é ®iÓm s−¬ng t¹i mùc 850mb (K0) Td700 lµ nhiÖt ®é ®iÓm s−¬ng t¹i mùc 700mb (K0) 3.3. ChØ sè TT: C«ng thøc to¸n häc cña chØ sè TT: TTTT dTT 500850500850 −+−= , trong ®ã: T850 lµ nhiÖt ®é t¹i mùc 850mb (K0) T500 lµ nhiÖt ®é t¹i mùc 500mb (K0) Td850 lµ nhiÖt ®é ®iÓm s−¬ng t¹i mùc 850mb (K0) Sè liÖu cho m« h×nh ETA ®−îc cËp nhËt tõ m« h×nh AVN hµng ngµy tõ 10h30 ®Õn 11h30. Tr−êng sè liÖu cña m« h×nh AVN cã ®é ph©n gi¶i ngang 10×10 kinh vÜ (t−¬ng Dù b¸o d«ng ë B¾c Bé 59 ®−¬ng 111km) trªn 26 mùc theo chiÒu th¼ng ®øng nh÷ng tÖp sè liÖu nµy ®−îc cung cÊp miÔn phÝ trªn Internet ë c¸c thêi ®iÓm 00Z, 06Z, 12Z vµ 18Z. Nh−ng tr−êng sè liÖu cña m« h×nh AVN ch−a phï hîp víi tr−êng sè liÖu cña m« h×nh ETA nªn ®−îc néi suy vÒ tr−êng sè liÖu phï hîp víi m« h×nh ETA b»ng mét ch−¬ng tr×nh Fortran. Qu¸ tr×nh ch¹y vµ hiÓn thÞ kÕt qu¶ cña m« h×nh lµ 1h30 phót, kÕt qu¶ dù b¸o c¸c tr−êng khÝ t−îng nh− khÝ ¸p, nhiÖt ®é, giã, ®é Èm ®−îc ghi ra tõng tiÕng mét. §Ó tÝnh ®−îc c¸c chØ sè d«ng Boyd, TT vµ KI ë c¸c thêi ®iÓm 1,2,3... giê chóng t«i cÇn lÊy ra c¸c tr−êng sè liÖu cña c¸c yÕu tè: Z1000, Z700, T850, Td850, T700, Td700, T500 t¹i thêi ®iÓm ®ã. M« h×nh ETA ch¹y dù b¸o cho 48 giê chóng t«i lÊy ra 336 tr−êng sè liÖu mçi tr−êng gåm 12285 ®iÓm (105 x117 víi 7 yÕu tè tõng giê mét. MiÒn tÝnh tõ kinh ®é tõ 96 ®Õn 1250 vµ vÜ ®é tõ 10 ®Õn 250, ®é ph©n gi¶i ngang lµ 0.25 x 0.25. §· lËp tr×nh tÝnh to¸n c¸c chØ sè Boyd, KI, TT cho tÊt c¶ c¸c ®iÓm (tøc 12285 ®iÓm) dùa vµo 7 yÕu tè: Z1000, Z700, T850, Td850, T700, Td700, T500 vµ ®−îc tr−êng chØ sè d«ng. Chóng t«i thu ®−îc chØ sè d«ng t¹i c¸c thêi ®iÓm c¸ch nhau 1 giê mét theo kÕt qu¶ dù b¸o 48 giê cña m« h×nh ETA. TiÕn hµnh néi suy chØ sè d«ng vÒ tr¹m ®Ó nhËn ®−îc c¸c chØ sè ë mçi giê t¹i c¸c tr¹m quan tr¾c .C¸c tr¹m lµ c¸c s©n bay cã quan tr¨c d«ng : Néi Bµi, C¸t Bi, §µ N½ng, §iÖn Biªn, Vinh, Phó Bµi, Nµ S¶n.Víi chuçi sè liÖu trªn ®· sö dông ph−¬ng ph¸p ph©n líp ®Ó t×m ph−¬ng tr×nh dù b¸o. 4. X©y dùng ph−¬ng tr×nh dù b¸o d«ng cho khu vùc B¾c Bé 4.1. Nguån sè liÖu: Chóng t«i xö lý kÕt qu¶ cña m« h×nh ETA vµ tÝnh chØ sè d«ng cho tõng ®iÓm tr¹m vµo tÊt c¶ c¸c giê 00h, 01h, 02h,..., 24h,…, 48h. Nh− vËy chóng t«i sÏ thu ®−îc chØ sè d«ng tõng giê mét trong hai ngµy t¹i 7 tr¹m. T¹i 7 tr¹m nµy th× cø 30 phót mét lÇn quan tr¾c d«ng. Chóng t«i xö lý kÕt qu¶ quan tr¾c theo quy ®Þnh sau: Tæng hîp 2 lÇn ph¸t b¸o liªn tiÕp l¹i tÝnh cho kÕt qu¶ quan tr¾c d«ng mét giê (ph¸t b¸o mét lÇn t¹i thêi ®iÓm trïng víi thêi ®iÓm tÝnh c¸c chØ sè vµ ph¸t b¸o lÇn liªn tiÕp cã nghÜa lµ 30 phót sau). B¶ng 1: Minh ho¹ sù quy ®Þnh hiÖn t−îng d«ng. Ph¸t b¸o t giê (t=1, 2…24) Ph¸t b¸o ( t+30 phót ) HiÖn t−îng d«ng (1giê) Cã d«ng Cã d«ng Cã d«ng Cã d«ng Kh«ng cã d«ng Cã d«ng Kh«ng cã d«ng Cã d«ng Cã d«ng Kh«ng cã d«ng Kh«ng cã d«ng Kh«ng cã d«ng Sau khi quy ®Þnh nh− trªn, chóng t«i lËp tr×nh ch−¬ng tr×nh ®äc sè liÖu ph¸t b¸o 7 tr¹m trªn, kÕt hîp víi sè liÖu tÝnh to¸n c¸c chØ sè Boyd, chØ sè KI vµ chØ sè TT ta lËp ®−îc chuçi sè liÖu ®Ó dù b¸o d«ng ë c¸c tr¹m. TrÇn T©n TiÕn, NguyÔn ThÕ Vinh, §Æng ThÞ An 60 M« h×nh ETA ch¹y tõ th¸ng 2 ®Õn th¸ng 5 chóng t«i thu ®−îc mét tËp sè liÖu nh− d· m« t¶. Sö dông ch−¬ng tr×nh ph©n ®Ó tÝnh to¸n víi chuçi sè liÖu ®· cã ®· t×m ®−îc hµm ph©n líp dù b¸o d«ng. 4. 2. Ph−¬ng ph¸p ®¸nh gi¸ kÕt qu¶ Ng−êi ta cã nhiÒu ph−¬ng ph¸p ®Ó ®¸nh gi¸ kÕt qu¶, nh−ng trong bµi b¸o vÒ dù b¸o d«ng nµy chóng t«i sö dông ph−¬ng ph¸p ®¸nh gi¸ kÕt qu¶ dù b¸o pha (cã hay kh«ng cã d«ng). §é chÝnh x¸c cña ph−¬ng ph¸p dù b¸o ®−îc x¸c ®Þnh theo c«ng thøc sau: N U NN 2211 += , trong ®ã: N: Tæng sè lÇn dù b¸o. N11: Sè lÇn dù b¸o ®óng pha cã d«ng. N22: Sè lÇn dù b¸o ®óng pha kh«ng cã d«ng. VÒ nguyªn t¾c, chän ph−¬ng ph¸p dù b¸o nµo cã U lín nhÊt. Thùc tÕ chØ sö dông U ®Ó ®¸nh gi¸ th× kh«ng ®ñ, v× vËy cÇn ph¶i ®¸nh gi¸ theo mét tiªu chuÈn n÷a gäi tiªu chuÈn ®é tin cËy H: U UUH 0 0 1− −= , trong ®ã: U0 lµ ®é chÝnh x¸c toµn phÇn cña dù b¸o ngÉu nhiªn. Víi U0 = 0,5. H > 0,2 th× kÕt luËn chØ tiªu nµy sö dông cho kÕt qu¶ tèt. 4.3. §¸nh gi¸ kÕt qu¶ Chóng t«i x¸c ®Þnh: - YÕu tè dù b¸o lµ hiÖn t−îng d«ng cã xuÊt hiÖn hay kh«ng. - Nh©n tè dù b¸o lµ ba chØ sè Boyd, KI vµ TT. Chóng t«i lÇn l−ît ®æ c¸c nh©n tè dù b¸o vµ ph©n líp thu ®−îc c¸c kÕt qu¶ nh− sau: 4.3.1. Sö dông mét chØ sè lµm nh©n tè dù b¸o: B¶ng 2: KÕt qu¶ ph©n líp cña chØ sè Boyd: U(%) H (®é tin cËy) Ng−ìng 63.985 (%) 0.3197 -176.004 B¶ng 3: KÕt qña ph©n líp cña chØ sè TT: U (%) H (®é tin cËy) Ng−ìng 59.032 (%) 0.18 43.236 Dù b¸o d«ng ë B¾c Bé 61 B¶ng 4: KÕt qu¶ ph©n líp cña chØ sè KI: U (%) H (®é tin cËy) Ng−ìng 55.040 (%) 0.1008 582.8 Sau khi ®æ vµo ch−¬ng tr×nh ph©n líp ta thu ®−îc ng−ìng cña c¸c chØ sè TT (43.236), chØ sè KI (582.8) vµ chØ sè Boyd (-176.004), trong ®ã hai chØ sè KI vµ TT dù b¸o kh«ng tèt b»ng chØ sè Boyd v× chØ sè Boyd dù b¸o ®óng 65.985% víi ®é tin cËy lµ 0.3197 trong khi ®ã hai chØ sè KI vµ TT chØ dù b¸o ®óng 55.040% vµ 59.032% víi ®é tin cËy nhá h¬n 0.2. V× vËy ta cã thÓ kÕt luËn ®−îc r»ng chØ sè Boyd cã thÓ ®−îc sö dông ®Ó dù b¸o d«ng. 4.3.2. Sö dông hai chØ sè lµm hai nh©n tè dù b¸o: B¶ng 5: KÕt qu¶ ph©n líp cña chØ sè Boyd vµ TT: U(%) H (§é tin cËy) Hµm ph©n líp 65.985% 0.3197 I=14.744TT + 6.199Boyd + 823.13 B¶ng 6: KÕt qu¶ ph©n líp cña chØ sè Boyd vµ KI: U(%) H (§é tin cËy) Hµm ph©n líp 69.051% 0.381 I=133.97KI + 6.2127Boyd – 76984.256 B¶ng 7: KÕt qu¶ ph©n líp cña chØ sè KI vµ TT: U(%) H (§é tin cËy) Hµm ph©n líp 58.905% 0.1781 I=22.094TT + 137.567KI – 81129.4 KÕt qu¶ ph©n líp thu ®−îc cho thÊy viÖc sö dông hai chØ sè lµm hai nh©n tè dù b¸o cho kÕt qu¶ kh¶ quan h¬n vµ kh¶ n¨ng dù b¸o chÝnh x¸c h¬n viÖc sö dông mét chØ sè lµm nh©n tè dù b¸o. KÕt qu¶ dù b¸o kÕt hîp hai chØ sè Boyd vµ KI cã thÓ ®óng tíi 69.051% víi hµm ph©n líp I = 133.97KI + 6.2127Boyd – 76984.256 vµ ®é tin cËy 0.381. Tæ hîp hai chØ sè Boyd vµ TT còng cho kÕt qu¶ tèt h¬n ®óng tíi 65.985% víi ®é tin cËy 0.3197 vµ còng t×m ®−îc hµm ph©n líp I = 14.744TT + 6.199Boyd + 823.13. ViÖc tæ hîp hai chØ sè KI vµ TT cho kÕt qu¶ kh«ng tèt chØ dù b¸o ®óng 58,905% víi ®é tin cËy 0.1781 Nh− vËy ta cã thÓ t×m ®−îc hai ph−¬ng tr×nh ph©n líp: • I = 133.97KI + 6.2127Boyd – 76984.256 • I = 14.744TT + 6.199Boyd + 823.13 Cã thÓ dù b¸o d«ng víi ®é tin cËy lín h¬n 0.2 TrÇn T©n TiÕn, NguyÔn ThÕ Vinh, §Æng ThÞ An 62 4.3.3. Sö dông ba nh©n tè lµm ba nh©n tè dù b¸o: B¶ng 9: KÕt qu¶ ph©n líp cña ba chØ sè Boyd, KI vµ TT: U(%) H (§é tin cËy) Hµm ph©n líp 69.270% 0.3854 I = 22.231TT + 139.27KI + 11.864Boyd - 80521.2 KÕt qu¶ ph©n líp cho thÊy sö dông c¶ ba chØ sè ®Ó dù b¸o cho kÕt qu¶ tèt nhÊt víi dù b¸o ®óng lµ 69.270% vµ ®é tin cËy 0.385 Hµm ph©n líp trªn ®· kiÓm nghiÖm trªn chuçi sè liÖu ®éc lËp cña 5 tr¹m trong 6 ngµy 06, 07, 14, 15, 25, 26 cña th¸ng 6 n¨m 2004 gåm 720 tr−êng hîp. KÕt qu¶ dù b¸o U = 68.841% víi ®é tin cËy H = 0.3768.Nh− vËy ph−¬ng tr×nh t×m ®−îc cã thÓ ®−a vµo thö nghiÖm cho dù b¸o nghiÖp vô ë ViÖt Nam. §©y lµ kÕt qu¶ nghiªn cøu cña ®Ò tµi CB 733104 .C«ng tr×nh hoµn thµnh víi sù hç trî cña ch−¬ng tr×nh nghiªn cøu c¬ b¶n n¨m 2005 . Tµi liÖu tham kh¶o 1. TrÇn C«ng Minh, KhÝ t−îng Syn«p phÇn NhiÖt ®íi, §¹i häc Quèc gia Hµ Néi, Hµ Néi, 2001. 2. TrÇn T©n TiÕn, §èi l−u khÝ quyÓn, §¹i häc Quèc gia Hµ Néi, Hµ Néi, 2002. 3. TrÇn T©n TiÕn, NguyÔn §¨ng QuÕ, Xö lý sè liÖu khÝ t−îng vµ dù b¸o thêi tiÕt b»ng ph−¬ng ph¸p thèng kª vËt lý, §¹i häc Quèc gia Hµ Néi, Hµ Néi, 2002. 4. Haklander A.J., Delden A.V. 5. Thunderstorm predictors and their forecast skill for the Netherlands 6. Atmospheric Research 67-68(2003) 273-299 7. Thunderstorm forecasting guide 8. VNU. JOURNAL OF SCIENCE, Nat., Sci., & Tech., T.xXII, n01AP., 2006 THUNDERSTORM FORECAST IN THE NORTH OF VIETNAM Tran Tan Tien, Nguyen The Vinh, Dang Thi An Department of Hydro-Meteorology & Oceanography College of Science, VNU Using prediction fields of geopotential, temperature, moistue given by ETA every hour, some unstable indexes of the atmosphere at a model grid point are calculated. The indexes include Boyd, KI and TT. These indexes are interpolated into the positions of weather stations, where thunderstorms are observed every 30 minute. The series of these indexes and thunderstorms are used to establish discriminative functions. The functions are tested using independent and dependent data. The results of the tests show that discriminative function with 3 variables may be used for thunderstorm forecast in the North of Vietnam in meteorological operations. Dù b¸o d«ng ë B¾c Bé 63