WV–PV–SWV prognoz-xatosini diagnostika qilish doirasi, haqiqiy Geo-Kompsat-2A (GK2A) suv bug‘i kuzatuvlarini, Korean Integrated Model (KIM) chiqishlaridan RTTOV yordamida yaratilgan simulyatsiya qilingan suv bug‘i tasvirlarini va modeldan hisoblangan potensial vortisit maydonlarini bir xil dinamik kontekstda taqqoslab, sonli ob-havo prognozi xatolarining joylashuv va intensivlik komponentlarini ajratishga qaratilgan sun’iy yo‘ldoshga asoslangan diagnostika usulidir. Tadqiqotda GK2A ning 6,3, 6,9 va 7,3 µm suv bug‘i kanallari ishlatilgan; KIM atmosfera maydonlari RTTOV 13-versiyasi yordamida sun’iy yo‘ldoshga ekvivalent yorqinlik harorati tasvirlariga aylantirilgan; kuzatilgan va simulyatsiya qilingan tuzilmalar orasidagi gorizontal joylashuv farqi MAPLE asosidagi Variational Echo Tracking yordamida hisoblangan va faza tuzatilgandan keyin qolgan yorqinlik harorati farqi intensivlik/amplituda xatosining ko‘rsatkichi sifatida baholangan.

Usulning miqdoriy tekshiruvi Sharqiy Osiyo ustida 1 Yanvar 2025-1 Yanvar 2026 davridagi KIM +0 dan +18 soatgacha bo‘lgan prognozlaridan foydalanib amalga oshirilgan. Faza tuzatilgandan keyin uchala suv bug‘i kanalida va barcha prognoz muddatlarida o‘rtacha kvadratik ildiz xatosi (RMSE) kamaygan, kuzatuv va simulyatsiya orasidagi fazoviy korrelyatsiya oshgan. Manba Jadval 4 ga ko‘ra RMSE yaxshilanish nisbatlari taxminan %18,59-%22,50; korrelyatsiya yaxshilanish nisbatlari esa %15,99-%27,48 oralig‘ida. Bu natijalar MAPLE orqali chiqarilgan gorizontal siljish vektorlari kuzatilgan va simulyatsiya qilingan suv bug‘i tuzilmalari orasidagi faza xatosining muhim qismini ifodalay olishini ko‘rsatadi.

Doiraning meteorologik ma’nosi faqat “sun’iy yo‘ldosh va model tasviri bir xilmi?” degan savolga qarashdan iborat emas. Simulyatsiya qilingan suv bug‘i tasvirining KIM ning o‘z potensial vortisit tuzilmasi bilan dinamik jihatdan mosligi ham tekshiriladi. Shu tariqa haqiqiy sun’iy yo‘ldosh kuzatuvi bilan model tasviri orasidagi farq ehtimoliy model prognoz xatosimi yoki PV bilan suv bug‘i o‘rtasidagi fizik bog‘lanishning o‘zi zaif bo‘lgan noaniq meteorologik holatmi, yanada ehtiyotkorlik bilan talqin qilinishi mumkin.

Usul ikki Aprel 2021 holatida qo‘llangan. 28 Aprel 2021 da Koreya yarimoroliga ta’sir qilgan Mo‘g‘uliston/shimoli-sharqiy Xitoy manbali botiqlik-vorteks tizimida sun’iy yo‘ldosh kuzatuvlari model prognoz qilganidan sekinroq siljigan va zaifroq rivojlangan tizimni ko‘rsatgan; keyingi 500 hPa tahlil maydonlari ushbu diagnostika bilan mos kelgan. 16-17 Aprel 2021 Bohai ko‘rfazi botiqligi holatida esa kuzatuvlar modelga qaraganda kuchliroq quruq havo kirib kelishi va kuchliroq nam ko‘tarilish/dinamik rivojlanishni ko‘rsatgan; keyinchalik kuzatilgan do‘l va mahalliy kuchli yog‘inning ham prognozdagidan kuchliroq rivojlangani bildirilgan.

Bu natijalar usul har bir meteorologik holatda aniq prognoz-xatosi diagnostikasini beradi degani emas. Mualliflar ayniqsa zaif yoki tartibsiz sinoptik tuzilmalarda PV bilan simulyatsiya qilingan suv bug‘i o‘rtasidagi moslik buzilishi mumkinligini; diabatik isish, chuqur konveksiya, bulut-radiatsiya jarayonlari va namlik noaniqliklari ham WV-PV farqlarini yuzaga keltirishi mumkinligini ta’kidlaydi. Bundan tashqari Type A/B/C tasnifi qisman ekspert talqiniga bog‘liq va doiraning kattaroq, mustaqil, ko‘p-holatli ma’lumotlar to‘plamlarida tizimli diagnostika qobiliyati hali to‘liq tasdiqlanmagan.

Sonli ob-havo prognozidagi “xato” nimani anglatadi?

Sonli ob-havo prognozi modellari atmosferani fizik tenglamalar va sonli usullar bilan ifodalaydi. Masalan, model yuqori troposferadagi quruq havo hududining umumiy shaklini to‘g‘ri prognoz qilsa ham uni haqiqiy joyidan 50 km sharqda ko‘rsatishi mumkin. Bunday holda muammo birinchi navbatda faza yoki joylashuv xatosidir. Boshqa holatda tuzilmaning joylashuvi to‘g‘ri bo‘lishi mumkin, ammo modeldagi quruq/nam signal haqiqatdagidan kuchliroq yoki zaifroq bo‘lishi mumkin. Bu esa amplituda yoki intensivlik xatosi deb qaralishi mumkin.

Bu ikki xato turining bir-biriga aralashishi sun’iy yo‘ldosh-model taqqoslashlarini qiyinlashtiradi. Modeldagi tuzilma haqiqiy tuzilmaga nisbatan gorizontal siljigan bo‘lsa, bir xil koordinatdagi piksel haroratlarini bevosita ayirish faqat intensivlik farqini emas, joylashuv farqini ham bitta son ichida birlashtiradi. Tadqiqotning muhim maqsadi avval gorizontal faza siljishini baholash, simulyatsiya qilingan tasvirni kuzatilgan tuzilma tomon siljitish va shundan keyingina qolgan yorqinlik harorati farqini o‘rganishdir.

WV–PV–SWV Doirasi Nima Uchun Xizmat Qiladi?

WV–PV–SWV doirasi haqiqiy sun’iy yo‘ldosh suv bug‘i tasvirini (WV), modeldan hosil qilingan sun’iy yo‘ldoshga ekvivalent suv bug‘i tasvirini (SWV) va modelning dinamik holatini ifodalovchi potensial vortisit maydonini (PV) birgalikda baholab, prognoz xatosini faqat tasvir farqi sifatida emas, balki atmosfera dinamikasiga mos tarzda talqin qilish imkonini beradi.

WV, SWV va PV nima?

WV (water-vapor imagery), GK2A sun’iy yo‘ldoshining suv bug‘i infraqizil kanallarida haqiqatda o‘lchagan radiatsiyadan hosil qilingan yorqinlik harorati tasviridir. Tadqiqotda ular kuzatuv referensi sifatida ishlatilgan.

SWV (simulated water-vapor imagery), KIM ning harorat, namlik, bulut va sirt holatlaridan foydalanib RTTOV radiatsion uzatish modeli orqali “sun’iy yo‘ldosh bu atmosferani ko‘rsa qanday radiatsiya/yorqinlik haroratini o‘lchardi?” degan savolga berilgan model-ekvivalent tasvirdir.

PV (potential vorticity; potensial vortisit), atmosferaning aylanishi va termodinamik qatlamlanish tuzilmasini birgalikda ifodalovchi dinamik kattalikdir. Tadqiqotda yuqori troposfera va tropopauza atrofidagi botiqlik, tizma, quruq havo kirib kelishi va dinamik rivojlanishni talqin qilishda fizik referens sifatida qo‘llanadi.

GK2A nega uch xil suv bug‘i kanalidan foydalanadi?

Tadqiqot GK2A Advanced Meteorological Imager ning 6,3, 6,9 va 7,3 µm suv bug‘i kanallaridan foydalanadi. Bu kanallar atmosferaning bir xil balandligini teng darajada ko‘rmaydi. Kanal og‘irlik funksiyalari turli bosim qatlamlariga sezgir; shu orqali yuqori troposferadan o‘rta troposferagacha turli qatlamlardagi namlik va oqim tuzilmalari kuzatilishi mumkin.

Manba Shakl 1 ushbu kanallarning normallashtirilgan og‘irlik funksiyalarini bosim o‘qi bo‘ylab sxematik ko‘rsatadi. Og‘irlik funksiyasi bitta bosim sathini o‘lchaydigan ingichka chiziq emas; u atmosfera qatlami bo‘ylab hissaning qanday taqsimlanishini ko‘rsatadi. Shuning uchun MAPLE orqali chiqarilgan gorizontal harakat vektorlari ham bitta aniq balandlikdagi shamol sifatida emas, kanal sezgir bo‘lgan qatlamning o‘rtacha tuzilma harakati sifatida talqin qilinishi kerak.

GK2A/AMI ning kuzatuv xususiyatlari

GEO-KOMPSAT-2A 5 Dekabr 2018 da uchirilgan va 24 Iyul 2019 da operatsion ishni boshlagan. AMI jami 16 kanalga ega. To‘liq disk kuzatuv sikli 10 daqiqa, Koreya yarimoroli hududi uchun kuzatuv sikli 2 daqiqadir. Tadqiqotning suv bug‘i tahlillarida ishlatilgan infraqizil WV kanallarining tabiiy fazoviy aniqligi 2 km.

KIM modelidan qaysi ma’lumotlar ishlatiladi?

Korean Integrated Model (KIM), tadqiqotda model-ekvivalent WV tasvirlari va PV maydonlari yaratiladigan operatsion global sonli ob-havo prognozi tizimidir. Manbada +0 tahlil vaqtidan +18 soat prognoz muddatigacha soatlik model maydonlari ishlatilgan.

KIM, 14 May 2025 gacha 12 km gorizontal aniqlik va 91 vertikal sath bilan ishlagan, shu sanadan keyin 8 km gorizontal aniqlikka o‘tkazilgan. Tadqiqotning bir yillik tekshiruv davri ushbu operatsion o‘zgarishni qamrab oladi. Manba Jadval 4 natijalarini 12 km va 8 km davrlari sifatida alohida ajratmaydi; shu sabab aniqlik o‘zgarishining tekshiruv ko‘rsatkichlariga alohida ta’sirini bu tadqiqotdan chiqarib bo‘lmaydi.

Muhim belgi qaydi: Manba bir xil u va v belgilarini ikki xil kontekstda ishlatadi. KIM ma’lumotlari bo‘limida u va v gorizontal shamol komponentlaridir. MAPLE-VET tenglamalarida esa u(x,y) va v(x,y) kuzatilgan va simulyatsiya qilingan tasvir orasidagi gorizontal tasvir siljishining zonal va meridional komponentlaridir. Bu ikki qo‘llanishni bir-biri bilan adashtirmaslik kerak.

Tadqiqotning usul oqimi

RTTOV nega kerak?

Sun’iy yo‘ldosh kuzatuvi bevosita model namlik maydoni bilan taqqoslansa, ikki xil fizik kattalik taqqoslangan bo‘ladi. GK2A radiatsiyani o‘lchaydi; KIM esa harorat, namlik, bulut va dinamik o‘zgaruvchilarni hisoblaydi. RTTOV model atmosferidan sun’iy yo‘ldosh sensori o‘lchashi mumkin bo‘lgan radiatsiya va yorqinlik haroratini simulyatsiya qilib, kuzatuv bilan modelni bir xil “radiatsiya fazosi”da taqqoslanadigan holga keltiradi.

Tenglama 1: Atmosfera tepasidagi radiatsiya

\[ L_c = \varepsilon_c B_c(T_s)\tau_c(p_s) + \int_{0}^{p_s} B_c(T(p)) \frac{\partial \tau_c(p)}{\partial p}\,dp + L_c^{cld}(x) \tag{1} \]

Formula nimani hisoblaydi? Kanal c uchun atmosfera tepasida sun’iy yo‘ldoshga yetib keladigan jami radiatsiyani ifodalovchi sxematik all-sky RTTOV ifodasidir.

• L_c: c kanalida sun’iy yo‘ldoshga yetib keladigan jami radiatsiya.
• ε_c: sirtning c kanalidagi emissivligi; o‘lchovsiz koeffitsiyent.
• B_c(T): tegishli kanal uchun Planck nurlanish funksiyasi.
• T_s: sirt harorati.
• p_s: sirt bosimi.
• τ_c(p): p bosim sathidan atmosfera tepasigacha kanal o‘tkazuvchanligi.
• T(p): bosimga bog‘liq atmosfera harorati.
• L_c^cld(x): operatsion all-sky konfiguratsiyada bulutlarning sof radiatsiya hissasi.

Birinchi had sirt chiqargan va atmosferadan o‘tib sun’iy yo‘ldoshga yetadigan hissani; integral had atmosfera qatlamlarining o‘z issiqlik nurlanishini; oxirgi had esa bulutlardan keladigan qo‘shimcha nurlanish ta’sirini ifodalaydi. Manba RTTOV all-sky simulyatsiyalarida KIM ning CLW, CIW va CF o‘zgaruvchilaridan foydalanadi; yomg‘ir va qor gidrometeor o‘zgaruvchilari bu qo‘llanmada kirish sifatida ishlatilmagan.

Tenglama 2: Radiatsiyadan yorqinlik haroratiga o‘tish

\[ T_{b,c}^{RTTOV}=B_c^{-1}(L_c) \tag{2} \]

Bu yerda B_c^-1, Planck funksiyasining teskarisidir. Hisoblangan radiatsiya sensor kanalida ayni radiatsiyani hosil qiladigan qora jism haroratiga ekvivalent yorqinlik harorati qiymatiga aylantiriladi. Yorqinlik harorati atmosferadagi bitta fizik nuqtaning termometr harorati emas; u kanal ko‘rgan nurlanishning harorat birliklarida ifodalangan ekvivalentidir.

Tenglama 3: Model holatidan simulyatsiya qilingan yorqinlik harorati

\[ T_{b,c,m}^{RTTOV}(t)=F_c(x_m(t)) \tag{3} \]

F_c, KIM model holat vektori x_m(t) ni kanal c uchun radiatsiya yoki yorqinlik haroratiga aylantiradigan RTTOV oldinga operatoridir. Boshqacha aytganda, model atmosferasi bevosita sun’iy yo‘ldosh tasviriga o‘xshash kuzatuv fazosiga ko‘chiriladi.

Tasvirlar nega 5 km umumiy gridga ko‘chiriladi?

GK2A va KIM ning tabiiy gridlari bir-biridan farq qiladi. MAPLE algoritmi bir xil koordinatlarda ikki tasvir tuzilmasini taqqoslay olishi uchun ikkala maydon ham bir xil geometrik griddan foydalanishi kerak. Tadqiqotda −21,5688° dan 68,7119° gacha kenglik va 45,8222°E dan 204,8646°E gacha uzunliklarni qamrab olgan katta Sharqiy Osiyo hududida umumiy 5 km kuzatuv gridi ishlatilgan va maydonlar bilinear interpolatsiya bilan qayta namunalangan.

MAPLE-VET nima qiladi?

MAPLE ning Variational Echo Tracking yondashuvi kuzatilgan tasvirdagi tuzilma bilan simulyatsiya tasviridagi o‘xshash tuzilma qaysi gorizontal siljishda eng yaxshi ustma-ust tushishini topishga harakat qiladi. Siljish maydoni:

\[ V(x,y)=\left(u(x,y),v(x,y)\right) \]

deb ta’riflanadi. Bu yerdagi u va v ushbu bo‘limda meteorologik shamol emas, tasvir siljishining ikki gorizontal komponentidir.

Tenglama 4: Umumiy MAPLE xarajat funksiyasi

\[ F(V)=F_z(V)+F_V(V) \tag{4} \]

Algoritm faqat ikki tasvirni imkon qadar o‘xshatishga harakat qilmaydi. U ikki maqsadni bir vaqtning o‘zida muvozanatlashtiradi:

• F_z: kuzatilgan va siljitilgan simulyatsiya tasviri orasidagi nomuvofiqlikni kamaytiradi.
• F_V: siljish maydonining fazoda haddan tashqari tartibsiz va shovqinli bo‘lishiga yo‘l qo‘ymaydigan silliqlik/regularization hadidir.

Tenglama 5: Tasvir nomuvofiqligi hadi

\[ F_z(V) = W_z \sum_x\sum_y \left[ Z_o(x,y) - Z_m\left(x-u(x,y),\,y-v(x,y)\right) \right]^2 \,dx\,dy \tag{5} \]

Z_o(x,y) kuzatilgan GK2A maydonini, Z_m simulyatsiya maydonini ifodalaydi. Simulyatsiya maydoni nomzod u,v siljish maydoniga ko‘ra orqaga qarab namunaladi. Ikki tasvir bir-biriga yaxshi mos tushsa kvadrat farq kamayadi. W_z ushbu tasvir saqlanishi/nomuvofiqligi hadining umumiy optimallashtirishdagi og‘irligidir.

Kvadrat olishning matematik natijasi katta nomuvofiqliklar kichik nomuvofiqliklarga qaraganda kuchliroq jazolanishidir. Biroq bu qiymat meteorologik “xato foizi” emas; MAPLE yechishi kerak bo‘lgan optimallashtirish funksiyasining bir qismidir.

Tenglama 6: Silliqlik regularization hadi

\[ F_V(V) = W_V \sum_x\sum_y \left[ \left(\frac{\partial u}{\partial x}\right)^2+ \left(\frac{\partial u}{\partial y}\right)^2+ 2\left(\frac{d^2u}{dx\,dy}\right)^2+ \left(\frac{\partial v}{\partial x}\right)^2+ \left(\frac{\partial v}{\partial y}\right)^2+ 2\left(\frac{d^2v}{dx\,dy}\right)^2 \right] dx\,dy \tag{6} \]

Bu had siljish vektorlarining qo‘shni nuqtalarda ma’nosiz ravishda keskin yo‘nalish va kattalik o‘zgarishlarini ko‘rsatishini jazolaydi. W_V kattalashgani sari yechim fazoviy jihatdan silliqroq bo‘ladi. Buning foydasi shovqindan kelib chiqadigan soxta vektorlarning bostirilishidir; biroq regularization haddan tashqari kuchaytirilsa, haqiqiy kichik miqyosli deformatsiyalar ham ortiqcha silliqlanishi mumkin.

Coarse-to-fine nega ishlatiladi?

Optimallashtirish to‘g‘ridan-to‘g‘ri juda zich griddan boshlansa, lokal minimumlarda qolib ketishi va kichik miqyosli shovqin katta harakatlarni to‘g‘ri ushlashni qiyinlashtirishi mumkin. Shu sabab tadqiqot avval qo‘pol miqyosda katta harakatni baholaydi, keyin natijani tobora zichroq gridlarga boshlang‘ich taxmin sifatida uzatadi.

Qo‘pol bosqichlarda Wz=1,0 va WV=100 ishlatilishi katta miqyosli tasvir harakatlarini ushlashga ko‘proq og‘irlik beradi. 100×100 va 300×300 nozik bosqichlarida Wz=0,5 gacha tushirilgan, WV esa 1000 gacha oshirilgan; shu tariqa kichik miqyosli soxta vektorlar, ma’lumot bo‘shliqlari va chegaralar atrofidagi noreal deformatsiyalarni bostirish kuchaytirilgan. Manba bu operatsion qiymatlar NMSC muhitida takroriy baholashlar natijasida empirik ravishda sozlanganini bildiradi.

Tenglama 7: Fazasi tuzatilgan simulyatsiya tasviri

\[ Z_m^{pc}(x,y) = Z_m\left(x-u(x,y),\,y-v(x,y)\right) \tag{7} \]

Bu amal model tasvirini MAPLE topgan gorizontal siljish miqdorida kuzatilgan tuzilma tomon ko‘chiradi. Maqsad model natijasini meteorologik jihatdan “tuzatish” emas, balki joylashuv xatosini vaqtincha bartaraf etib, qolgan intensivlik farqini alohida tekshirishdir.

Tenglama 8: Faza tuzatilgandan keyingi yorqinlik harorati farqi

\[ \Delta T_b^{pc}(x,y) = T_b^{GK2A}(x,y) - T_{b,m}^{RTTOV,pc}(x,y) \tag{8} \]

ΔT_b^pc, gorizontal joylashuv farqi imkon qadar bartaraf etilgandan keyin kuzatilgan va simulyatsiya qilingan yorqinlik harorati orasidagi farqdir. Shu tariqa bir xil son ichida ham “tuzilma boshqa joyda”, ham “tuzilma boshqa intensivlikda” ta’sirlarini aralashtirish o‘rniga ikkinchi komponentni aniqroq o‘rganish mumkin.

Musbat yoki manfiy ΔT_b^pc ning fizik talqini tasvir tuzilmasining quruq/nam xususiyati va meteorologik muhit bilan birgalikda bajarilishi kerak. Tadqiqot holat tahlillarida bu farqlarni ko‘tarilish va cho‘kish/subsidence jarayonlarining modeldagidan kuchliroq yoki zaifroq bo‘lishi haqidagi talqinlar bilan birga ishlatadi; faqat ΔT belgisi barcha vaziyatlar uchun kontekstdan mustaqil universal tasniflagich sifatida taqdim etilmaydi.

Potensial Vortisit Nega Ishlatiladi?

Potensial vortisit suv bug‘i tasviridagi quruq va nam tuzilmalarni faqat vizual naqshlar sifatida emas, balki yuqori troposfera va tropopauza atrofidagi botiqlik, tizma va dinamik rivojlanish bilan bog‘liq fizik atmosfera tuzilmalari sifatida baholash imkonini beradi; shu sabab u model tasvirining o‘z dinamik maydoni bilan mosligini tekshiruvchi referens qatlam vazifasini bajaradi.

Tenglama 9: Gidrostatik Ertel potensial vortisiti

\[ PV = -g \left[ (\zeta_p+f)\frac{\partial\theta}{\partial p} + \frac{\partial u}{\partial p}\frac{\partial\theta}{\partial y} - \frac{\partial v}{\partial p}\frac{\partial\theta}{\partial x} \right] \tag{9} \]

• g: erkin tushish tezlanishi.
• f: Coriolis parametri.
• ζ_p: izobarik sirt ustidagi nisbiy vortisit.
• θ: potensial harorat.
• p: bosim.
• u, v: ushbu tenglamada KIM modelining gorizontal shamol komponentlari.

Birinchi katta had atmosfera aylanishini potensial haroratning vertikal bosim gradienti bilan bog‘laydi. Qolgan ikki had vertikal shamol kesimi bilan potensial haroratning gorizontal gradientlari orasidagi baroklinik hissani o‘z ichiga oladi.

PV birligi PVU:

\[ 1\;PVU = 10^{-6}\;K\,m^2\,kg^{-1}\,s^{-1} \]

Tadqiqotda oltita PV ko‘rinishi hisoblangan: 400, 500 va 600 hPa izobarik sathlaridagi PV qiymatlari hamda 1,0, 1,5 va 2,0 PVU izosirtlarining bosimlari. Manba 1,5 PVU sirtining tropopauza uchun keng ishlatiladigan dinamik ta’rif ekanini bildiradi.

PV va suv bug‘i nega har doim bir-biriga aynan mos kelmaydi?

Tadqiqotning muhim metodologik ogohlantirishlaridan biri shu. Kuchli va yaxshi tashkil topgan botiqlik yoki tizma tizimlarida PV bilan simulyatsiya qilingan WV tuzilmalari o‘rtasidagi bog‘lanish odatda aniqroq. Ammo zaif tizimlarda, chuqur konveksiya paytida, kuchli yashirin issiqlik ajralishi bo‘lganda, bulut-radiatsiya jarayonlari yoki namlik noaniqligi kuchayganda WV bilan PV o‘rtasida farqlar yuzaga kelishi mumkin.

Shu sabab “sun’iy yo‘ldosh WV model PV dan farq qiladi → model albatta xato” shaklida avtomatik qaror berilmaydi. Avval modelning o‘z SWV tasviri o‘z PV tuzilmasi bilan fizik jihatdan mos yoki mos emasligi tekshiriladi.

Type A, B va C sinflari nimani anglatadi?

Manba Jadval 5 da ishlatilgan “=” va “≠” belgilar bu yerda ikki maydon son jihatdan aynan tengligini emas, balki tahlil qilinayotgan meteorologik tuzilma bo‘yicha dinamik/vizual moslik yoki nomuvofiqlik mavjudligini ko‘rsatuvchi tasnif belgilari.

Type C ning kuchli bo‘lish mantiqi quyidagicha: model o‘z ichida mos; ya’ni model atmosfera dinamikasi hosil qilgan PV tuzilmasi bilan ayni model atmosferasidan RTTOV orqali yaratilgan SWV tasviri bir-birini qo‘llab-quvvatlaydi. Biroq haqiqiy sun’iy yo‘ldosh kuzatuvi bu dinamik rivojlanishdan farq qiladi. Bunday vaziyatda farq faqat model ichidagi SWV-PV nomuvofiqligidan kelib chiqqanini aytish qiyinlashadi va haqiqiy atmosferaning modeldan boshqacha evolyutsiya qilgani haqida kuchliroq diagnostika berish mumkin.

80 km va 2 °C chegaralari universalmikan?

Yo‘q. Tadqiqotda faza siljishi uchun 80 km va yorqinlik harorati farqi uchun 2 °C atrofidagi chegaralar o‘rganilgan holatlar asosida empirik ravishda aniqlangan. Mualliflar ularni kengroq meteorologik sharoitlarda qo‘llashdan oldin qo‘shimcha optimallashtirish va baholash zarurligini ochiq aytadi. Bu qiymatlar universal meteorologik konstantalar emas.

Bir yillik tekshiruv qanday amalga oshirildi?

Faza vektorlari ikki taqqoslash orqali tekshirilgan:

1. GK2A kuzatuvi bilan original RTTOV simulyatsiya tasviri orasidagi yorqinlik harorati farqi.
2. Xuddi shu GK2A kuzatuvi bilan MAPLE siljishi qo‘llangan, fazasi tuzatilgan RTTOV tasviri orasidagi farq.

Agar faza vektori to‘g‘ri bo‘lsa, simulyatsiya tasviri kuzatilgan tuzilmaga yaxshiroq mos tushishi, natijada RMSE kamayishi va fazoviy korrelyatsiya oshishi kerak. Manbada ishlatilgan yaxshilanish nisbati:

\[ \text{İyileştirme oranı} = 100\times \left| \frac{(ii)-(i)}{(i)} \right| \]

shaklida berilgan. Bu yerda (i) tuzatishdan oldingi, (ii) faza tuzatishidan keyingi qiymat. RMSE uchun pastroq qiymat, korrelyatsiya uchun esa yuqoriroq qiymat yaxshilanishni bildiradi.

Jadvalning asosiy xabari yaxshilanish faqat bitta prognoz soatida ko‘rilgani emas. Faza tuzatishi tekshirilgan uch kanalning va +0 dan +18 soatgacha bo‘lgan barcha prognoz muddatlarining barchasida bir xil yo‘nalishdagi ta’sir ko‘rsatgan: RMSE kamaygan, fazoviy korrelyatsiya oshgan. Bu MAPLE ushlagan gorizontal siljishlar tasvir farqlarining haqiqiy bir qismini tushuntirishini qo‘llab-quvvatlaydi.

Eng yuqori RMSE yaxshilanishi %22,50 bilan 6,3 µm kanalida t+06 soatda kuzatilgan. Korrelyatsiya yaxshilanishidagi eng yuqori qiymat %27,48 bilan 7,3 µm kanalida t+18 soatda. Ushbu maksimumlarning turli kanallarda va turli prognoz muddatlarida yuzaga chiqishi kanallarning turli vertikal sezgirligi va prognoz xatolarining vaqt davomida rivojlanishi bilan bog‘lanadi.

Holat 1: 28 Aprel 2021 botiqligi nimani ko‘rsatdi?

Birinchi holat shimoli-sharqiy Xitoy ustidagi cut-off low tizimining quruq havo hududi bilan birga Koreya yarimoroli tomon rivojlanishini o‘rganadi. Tahlilda 6,3 µm WV tasvirlari, 1,5 PVU sirti va zarur hollarda 500 hPa PV maydonlari ishlatilgan.

Region 1 da kuzatilgan WV tuzilmasi modelga nisbatan taxminan 30-50 km shimoli-sharqda; Region 3 da taxminan 20-50 km shimol-shimoli-sharqda; Region 4 da taxminan 40-50 km shimoli-sharqda topilgan. Bu hududlarda fazasi tuzatilgan yorqinlik harorati farqlari taxminan −1 dan −2 °C gacha deb xabar qilingan va kuzatilgan quruq signal modeldagidan zaifroq deb talqin qilingan. Bu hududlar Type C xususiyatini ko‘rsatgan.

Region 2 Type B2 sifatida baholangan. Bu yerdagi faza farqi 30 km dan kichik va mutlaq yorqinlik harorati farqi 2 °C ichida bo‘lgani sababli, manba bu farqlar ma’noli va kuchli prognoz-xatosi talqinini berish uchun juda kichikligini bildirgan.

Umumiy diagnostika haqiqiy atmosferadagi siklon bilan bog‘liq botiqlik model prognoz qilganidan sekinroq aylangan/siljigan va zaifroq rivojlangan degan yo‘nalishda. Keyingi 500 hPa tahlil maydoni model prognoziga nisbatan biroz shimolroqda, iliqroq va zaifroq vortisitga ega tizimni ko‘rsatib, sun’iy yo‘ldoshga asoslangan diagnostikani qo‘llab-quvvatlagan. Manba 29-30 Aprelda ro‘y bergan yog‘inli xavfli ob-havo ham prognoz qilinganidan zaifroq bo‘lganini bildiradi.

Holat 2: 16-17 Aprel 2021 Bohai ko‘rfazi botiqligi

Ikkinchi holat shimoli-sharqiy Xitoydagi cut-off low bilan bog‘liq qisqa to‘lqinli botiqlik Bohai ko‘rfazi orqali Koreya yarimoroli tomon siljigan va kuchli ob-havo hodisalari bilan bog‘liq bo‘lgan 16-17 Aprel 2021 davridir. Hodisa davomida aniq quruq havo kirib kelishi, botiqlik oldida nam frontal bulut tasmasi, do‘l va mahalliy kuchli yog‘in qayd etilgan.

Region 2 da kuzatilgan quruq hudud simulyatsiya qilingan quruq hududdan taxminan 0-20 km torroq bo‘lgan. ΔT_b^pc qiymatlari odatda +1 dan +2 K gacha; quruq-nam chegara yaqinida esa taxminan +1 dan +4 K gacha yetgan. Manba buni kuzatuvlardagi quruq havo kirishi va cho‘kish harakati modeldagidan kuchliroq ekanining ko‘rsatkichi sifatida talqin qiladi.

Region 3 da kuzatilgan nam hududning model tasviriga nisbatan faza farqi taxminan 0-10 km va ΔT_b^pc taxminan −1 K. Bu natija nam ko‘tarilish model prognozidan kuchliroq bo‘lganini ko‘rsatgan.

Har ikki hudud Type C sifatida tasniflangan: modelning SWV tasviri PV tuzilmasi bilan mos qolgan, haqiqiy GK2A WV tasviri esa boshqacha rivojlangan. Keyingi 500 hPa tahlil maydonlari pastroq harorat va kuchliroq vortisit rivojlanishini ko‘rsatib, sun’iy yo‘ldosh diagnostikasidagi “modeldan kuchliroq haqiqiy dinamik rivojlanish” xulosasini qo‘llab-quvvatlagan. Haqiqatda yuz bergan do‘l va mahalliy kuchli yog‘in ham model prognozidan kuchliroq rivojlangan.

Ikki holat bir-biriga qarama-qarshi natija berganda usul nega muhim?

Birinchi holatda model tizimni haddan tashqari kuchli va tez rivojlantirgan, ikkinchi holatda esa model haqiqiy atmosferadagi rivojlanishni yetarlicha kuchli ifodalay olmagan. Ya’ni usul “model doimo juda kuchli” yoki “model doimo zaif” kabi bir yo‘nalishli biasni taklif qilmaydi. Ayni doira turli holatlarda xatoning yo‘nalishini dinamik kontekstga ko‘ra diagnostika qilishga harakat qiladi.

Tadqiqot qo‘llab-quvvatlaydigan xulosalar

• MAPLE faza tuzatishi bir yillik tekshiruvda barcha uch WV kanalida va +0 dan +18 soatgacha bo‘lgan barcha prognoz muddatlarida RMSE ni kamaytirgan va fazoviy korrelyatsiyani oshirgan.
• Kuzatilgan WV bilan model-ekvivalent SWV ni bir xil radiatsiya/yorqinlik harorati fazosida taqqoslash gorizontal siljish va intensivlik xatolarini obyektiv ajratishni qo‘llab-quvvatlaydi.
• PV qatlami suv bug‘i tasvirlaridagi farqlarni yuqori sath atmosfera dinamikasi bilan birga talqin qilish imkonini beradi.
• Type C sharoitlari modelning SWV-PV aloqasi o‘z ichida mos bo‘lgani uchun model prognozi kamchiligi eng aniq diagnostika qilinadigan sinf sifatida ajralib turadi.
• Ikki Aprel 2021 holatida sun’iy yo‘ldoshga asoslangan diagnostikalar keyingi 500 hPa tahlil maydonlarining umumiy rivojlanishi bilan mos kelgan.

Tadqiqot qo‘llab-quvvatlamaydigan umumlashtirishlar

• Doira barcha ob-havo hodisalarida va barcha geografiyalarda bir xil mahorat bilan ishlaydi degan xulosa chiqarib bo‘lmaydi.
• 80 km va 2 °C chegaralari universal optimal qiymatlar emas.
• Har bir WV-PV nomuvofiqligi bevosita NWP model xatosini anglatmaydi.
• Type B sinflarida ishonchli prognoz-xatosi diagnostikasi har doim ham mumkin emas.
• Ikki holat tadqiqoti o‘z-o‘zicha uzoq muddatli operatsion prognoz muvaffaqiyatini isbotlamaydi.
• Tadqiqot doira yakuniy avtomatik tasniflagich ekanini ko‘rsatmaydi; ekspert talqini hali ham zarur.
• Mustaqil kuzatuvlar va keng ko‘p-holatli ma’lumotlar to‘plamlari bilan kompleks diagnostika qobiliyati hali to‘liq baholanmagan.

Tadqiqot Usuli va Natijalari

Tadqiqot dizayni

Tadqiqot kuzatuv sun’iy yo‘ldosh ma’lumotlari, operatsion sonli ob-havo prognozi modeli, oldinga radiatsion uzatish simulyatsiyasi, tasvir moslashtirish/optimallashtirish usuli va sinoptik dinamik tahlilni birlashtirgan hisoblash atmosfera fanlari tadqiqotidir.

Ma’lumot manbalari

• Haqiqiy kuzatuv: GK2A/AMI 6,3, 6,9 va 7,3 µm suv bug‘i yorqinlik haroratlari.
• NWP modeli: Korean Integrated Model (KIM).
• Simulyatsiya qilingan sun’iy yo‘ldosh tasviri: RTTOV 13-versiya, all-sky konfiguratsiya.
• Faza farqi: MAPLE asosidagi Variational Echo Tracking.
• Dinamik referens: Gidrostatik Ertel PV; 400/500/600 hPa va 1,0/1,5/2,0 PVU ko‘rinishlari.

Tekshiruv hududi va vaqt tuzilishi

Miqdoriy faza-tasvir tekshiruvi Sharqiy Osiyo ustida 1 Yanvar 2025 dan 1 Yanvar 2026 gacha amalga oshirilgan. KIM prognozlari t+00 tahlildan t+18 soatgacha soatlik baholangan. PV kompozitlari sinoptik miqyos sabab 00-21 UTC oralig‘ida uch soatlik intervallarda hisoblangan.

KIM konfiguratsiyasi

RTTOV all-sky sharoitlari

RTTOV simulyatsiyalari ham ochiq osmon, ham bulutli sahnalarda ishlatilgan. KIM ning CLW, CIW va CF o‘zgaruvchilari bulut radiativ ta’sirlarini ifodalash uchun qo‘llangan. Manba yomg‘ir va qor gidrometeor o‘zgaruvchilari ushbu infraqizil WV simulyatsiyalariga kiritilmaganini bildiradi. Suyuq bulut zarrachalarining optik xususiyatlari RTTOV ning OPAC asosidagi suv buluti parametrizatsiyasi, muz buluti xususiyatlari esa Baran muz optik parametrizatsiyasi bilan ifodalangan.

Faza tuzatishining miqdoriy natijasi

Bir yillik tekshiruvda uchala WV kanalida faza tuzatilgandan keyin RMSE kamayishi va korrelyatsiya oshishi kuzatilgan. RMSE yaxshilanish nisbatlari taxminan %18,59 dan %22,50 gacha; korrelyatsiya yaxshilanish nisbatlari %15,99 dan %27,48 gacha o‘zgaradi.

Bu natija faza tuzatishi sun’iy yo‘ldosh-model tasvirlaridagi jami xato muammosini to‘liq hal qiladi degani emas. Amal faqat gorizontal siljish komponentini kamaytirishga qaratilgan; faza tuzatilgandan keyin qolgan ΔT_b^pc farqlari modelning intensivlik/termodinamik ifodasidagi ehtimoliy farqlarni o‘rganish uchun saqlanadi.

Holat natijalarini taqqoslash

Cheklovlar

1. PV-SWV mosligi sharoitga bog‘liq. Yaxshi rivojlangan sinoptik botiqlik va tizmalarda bog‘lanish aniqroq, zaif va tartibsiz tizimlarda esa ishonchlilik kamayishi mumkin.

2. Diabatik jarayonlar chalkashlik keltirib chiqarishi mumkin. Yashirin issiqlik ajralishi, chuqur konveksiya, bulut radiativ ta’sirlari va namlik noaniqligi WV-PV farqlarini yuzaga keltirishi mumkin; ularning hammasi bevosita prognoz xatosi emas.

3. Type B sinflari noaniq. Model SWV si model PV si bilan allaqachon nomuvofiq bo‘lsa, haqiqiy sun’iy yo‘ldosh bilan model orasidagi farq manbasini aniq ajratish qiyinroq.

4. Chegaralar empirik. 80 km faza siljishi va 2 °C yorqinlik harorati farqi chegaralari o‘rganilgan holatlarga ko‘ra aniqlangan bo‘lib, boshqa meteorologik sharoitlarda qayta optimallashtirilishi kerak.

5. Ekspert talqini to‘liq yo‘qolmagan. MAPLE va ΔT_b^pc obyektivlikni oshirsa ham, Type A/B/C tasnifining yakuniy bosqichi birlashgan WV-PV-SWV munosabatlarining ekspert talqiniga qisman bog‘liq.

6. Mustaqil keng qamrovli tekshiruv yetishmaydi. Mualliflar tasnif tizimining mustaqil kuzatuvlar, tahlil maydonlari, forecast-error metriklari yoki analysis increment ma’lumotlari bilan keng ko‘p-holatli tizimli mahorat baholashi hali bajarilmaganini ochiq aytadi.

7. Holatlar soni cheklangan. Operatsion qo‘llanish ikki vakillik qiluvchi Aprel 2021 hodisasi orqali ko‘rsatilgan. Turli fasllar, bulut muhitlari va sinoptik rejimlar uchun kengroq ma’lumot bazasi zarur.

Kelajakdagi ish

Mualliflar kattaroq ko‘p-holatli ma’lumotlar to‘plamlarida miqdoriy fazoviy o‘xshashlik metriklari, siljish diagnostikalari va statistik xususiyatlarni ajratib olishdan foydalanib, yanada obyektiv va avtomatik tasnif usullarini ishlab chiqishni rejalashtirmoqda. Shuningdek, uch WV kanalini 1,0-2,0 PVU sirtlari va 400-700 hPa PV maydonlari bilan birgalikda qo‘llaydigan ko‘p qatlamli, yanada keng qamrovli uch o‘lchovli sinoptik tahlillarni rivojlantirish tavsiya etiladi.

Manba va Usul Haqida Izoh

Asl tadqiqot: Diagnosing Numerical Weather Prediction Forecast Errors Using Geo-Kompsat-2A Observed and Simulated Water–Vapor Imagery Within A Potential-Vorticity Dynamical Framework

Mualliflar: Jun-Hyung Heo; Boram Kim; Minsang Kim; In-Chul Shin; Eun-Ha Sohn; KiRyong Kang.

Corresponding author: Jun-Hyung Heo.

ORCID: Springer Nature yozuvi Jun-Hyung Heo uchun 0000-0002-2051-3861 ORCID ni ko‘rsatadi.

Muassasa 1: National Meteorological Satellite Center, Korea Meteorological Administration, Jincheon-Gun 27803, Korea.

Muassasa 2: Korea Ocean Satellite Center, Korea Institute of Ocean Science and Technology, Busan Metropolitan City 49111, Korea.

DOI: 10.1007/s13143-026-00451-w

Jurnal: Asia-Pacific Journal of Atmospheric Sciences.

Jild / maqola raqami: Volume 62, Article 32.

Nashriyot: Springer Nature; jurnal Korean Meteorological Society bilan bog‘liq.

Manba turi: Original Article.

Qabul qilinishi: 6 Mart 2026.

Reviziya: 29 Iyul 2026.

Qabul: 7 Avgust 2026.

Nashr / version of record: 27 Avgust 2026.

Taqriz holati: Jurnalda chop etilgan version-of-record tadqiqot maqolasi; preprint sifatida baholanmagan.

Litsenziya: Creative Commons Attribution 4.0 International (CC BY 4.0). Manba tegishli iqtibos, litsenziya havolasi va o‘zgartirishlar haqida ma’lumot ko‘rsatilishi sharti bilan foydalanish, ulashish, moslashtirish, tarqatish va ko‘paytirishga ruxsat beradi. Uchinchi tomon materiallari uchun alohida kredit satri mavjud bo‘lsa, o‘sha materialning alohida huquqiy holati amal qiladi.

Moliyalashtirish: Tadqiqot Korea Meteorological Administration National Meteorological Satellite Center tomonidan moliyalashtirilgan “Development of Meteorological Satellite Operation and Application Technology” dasturi doirasidagi “Development of Operational Support Technologies for Satellite Utilization in Response to Meteorological Disasters” loyihasi tomonidan qo‘llab-quvvatlangan; grant raqami KMA2025-00121.

Etika: Manba inson ishtirokchilari yoki hayvonlarni o‘z ichiga olgan tadqiqot bo‘lmaganini bildiradi.

Manfaatlar to‘qnashuvi: Mualliflar manfaatlar to‘qnashuvi yo‘qligini bildirgan.

Ma’lumotlar mavjudligi: Manba PDF da alohida “Data Availability” yoki “Code Availability” bayonoti aniqlanmagan. Tadqiqot operatsion GK2A va KIM mahsulotlari, RTTOV va MAPLE asosidagi usullardan foydalanadi; shu sabab Verianla xom ma’lumot yoki kod ommaga ochiq deb taxmin qilmaydi.

Muallif hissalari: Jun-Hyung Heo tadqiqotni loyihalagan. Jun-Hyung Heo, Boram Kim, Minsang Kim va In-Chul Shin eksperimental/hisoblash qo‘llanmalarini bajargan. Jun-Hyung Heo, Minsang Kim va Boram Kim ma’lumotlarni tahlil qilgan. Jun-Hyung Heo, In-Chul Shin va Eun-Ha Sohn resurs va tahlil vositalarini taqdim etgan. Jun-Hyung Heo va Boram Kim maqolaning dastlabki qoralamasini yozgan. Eun-Ha Sohn va KiRyong Kang tadqiqotni nazorat qilgan va loyihani boshqargan. Barcha mualliflar yakuniy maqolani ko‘rib chiqib tasdiqlagan.

Usul xulosasi: GK2A/AMI haqiqiy WV kuzatuvlari → KIM atmosfera maydonlari → RTTOV v13 all-sky simulyatsiya WV → umumiy 5 km grid → MAPLE-VET faza siljish vektorlari → coarse-to-fine yaxshilash → fazasi tuzatilgan SWV → ΔT_b^pc intensivlik farqi → KIM PV maydonlari bilan dinamik moslik tekshiruvi → Type A/B/C prognoz-xatosi diagnostikasi.

Miqdoriy tekshiruv: Sharqiy Osiyo, 1 Yanvar 2025-1 Yanvar 2026; uch GK2A WV kanali; KIM t+00-t+18 soat prognozlari.

Asosiy talqin chegarasi: Faza vektorlarining miqdoriy tekshiruvi kuchli bo‘lsa ham, WV-PV-SWV Type A/B/C doirasining mustaqil ko‘p-holatli diagnostika qobiliyati hali to‘liq tizimli tekshirilmagan. Tadqiqot ikki vakillik qiluvchi yuqori ta’sirli ob-havo hodisasi bilan operatsion qo‘llash mumkinligini ko‘rsatadi.

KIM aniqlik qaydi: Manba KIM 14 May 2025 gacha 12 km, undan keyin 8 km gorizontal aniqlikda ishlaganini bildiradi. Bir yillik tekshiruv davri ikkala rejimni ham qamrab oladi; manba Jadval 4 natijalarini aniqlik davrlariga ko‘ra alohida hisobot qilmaydi.