Premium
Accelerating radiative transfer calculations for high‐resolution atmospheric models
Author(s) -
Barker Howard W.,
Li Jiangnan
Publication year - 2019
Publication title -
quarterly journal of the royal meteorological society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.744
H-Index - 143
eISSN - 1477-870X
pISSN - 0035-9009
DOI - 10.1002/qj.3543
Subject(s) - radiative transfer , radiative flux , column (typography) , atmospheric radiative transfer codes , meteorology , environmental science , physics , mathematics , geometry , optics , connection (principal bundle)
Abstract This study presents a method for calculating all‐sky, broadband radiative fluxes for high‐resolution atmospheric model domains, such as limited‐area numerical weather prediction models and Cloud System‐Resolving Models, that have many columns (i.e. N ≫ 10 5 ). The aim is to replace the Independent Column Approximation (ICA) which applies 1D radiative transfer models (RTMs) to all N columns; but usually after skipping many dynamical time steps due to the ICA's computational burden. The proposed method, referred to as Partitioned Gauss–Legendre Quadrature (PGLQ), begins by partitioning a domain into M sub‐domains. This study is concerned just with cloudy columns so partitioning was done according to cloud‐top altitudes. Then, for each sub‐domain, sort columns from smallest to largest cloud water path W , assign each column a sequence number, identify n G GLQ points, and apply 1D RTMs to the n G columns with resulting flux profiles assigned to nearby columns in the sorted sequence. Last, reposition columns back to the full domain. Marine boundary‐layer and deep convective cloud fields were used to assess PGLQ. For ratios of RTM calculations between ICA and PGLQ of f RT = N /( n G · M ) −1 ∈ [10 3 , 10 4 ], MBEs and RMSEs for net short‐wave (SW) surface fluxesF SW NET SFCscale close to ( n G · M ) −1 , RMSEs being typically ∼200× smaller than domain‐averageF LW NET SFC. Structure function analyses ofF SW NET SFCindicate that PGLQ errors resemble white noise. Long‐wave (LW) fluxes perform less well due to sorting on W but seem likely to be acceptable, RMSEs being ∼5–10× smaller than domain‐averageF SW NET SFC. All MBEs are typically <±0.5 W/m 2 . Estimates of radiative heating rate profiles are unbiased, and while noisy at small scales, they are portrayed well at scales larger than typical individual clouds. Overall speed‐up of PGLQ relative to ICA, for serial computation as consider here, is expected to be >300×; likely to be comparable to parallel computation. This should allow application of RTMs at every dynamical time step.