Computation of longwave radiative flux and vertical heating rate with 4A-Flux v1.0 as integral part of the radiative transfer code 4A/OP v1.5
Abstract. Based on advanced spectroscopic databases, line-by-line and layer-by-layer radiative transfer codes numerically solve the radiative transfer equation with a very high accuracy. Taking advantage of its pre-calculated optical depth look-up table, the fast and accurate radiative transfer model Automatized Atmospheric Absorption Atlas OPerational (4A/OP) calculates the transmission and radiance spectra for a user defined layered atmospheric model. Here we present a module, called 4A-Flux, developed and implemented into 4A/OP in order to include the calculation of the clear-sky longwave radiative flux profiles and heating rate profiles at a very high spectral resolution. Calculations are performed under the assumption of local thermodynamic equilibrium, plane-parallel atmosphere and specular reflection on the surface. The computation takes advantage of pre-tabulated exponential integral functions that are used instead of a classic angular quadrature. Furthermore, the sublayer variation of the Planck function is implemented to better represent the emission of layers with a high optical depth. Thanks to the implementation of 4A-Flux, 4A/OP model have participated in the Radiative Forcing Model Intercomparison Project (RFMIP-IRF) along with other state-of-the-art radiative transfer models. 4A/OP hemispheric flux profiles are compared to other models over the 1800 representative atmospheric situations of RFMIP, yielding an Outgoing Longwave Radiation (OLR) mean difference between 4A/OP and other models of −0.148 W .m−2 and a mean standard deviation of 0.218 W .m−2, showing a good agreement between 4A/OP and other models. 4A/OP is applied to the Thermodynamic Initial Guess Retrieval (TIGR) atmospheric database to analyze the response of the OLR and vertical heating rate to several perturbations of temperature or gas concentration. This work shows that 4A/OP with 4A-Flux module can successfully be used to simulate accurate flux and heating rate profiles and provide useful sensitivity studies including sensitivities to minor trace gases such as HFC134a, HCFC22 and CFC113. We also highlight the interest for the modeling community to extend intercomparison between models to comparisons between spectroscopic databases and modelling to improve the confidence in model simulations.
