Accelerating Radiative Transfer Calculations for High‐resolution Atmospheric Models
Howard W. Barker,Jiangnan Li
DOI: https://doi.org/10.1002/qj.3543
2019-01-01
Quarterly Journal of the Royal Meteorological Society
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 ≫ 105). 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 nG GLQ points, and apply 1D RTMs to the nG 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 fRT = N /(nG · M )−1 ∈ [103, 104], MBEs and RMSEs for net short‐wave (SW) surface fluxes FSWNETSFC scale close to (nG · M )−1, RMSEs being typically ∼200× smaller than domain‐average FLWNETSFC . Structure function analyses of FSWNETSFC indicate 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‐average FSWNETSFC . All MBEs are typically <±0.5 W/m2. 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.