Observation Operator for Visible and Near-Infrared Satellite Reflectances

Abstract

Abstract Operational numerical weather prediction systems currently only assimilate infrared and microwave satellite observations, whereas visible and near-infrared reflectances that comprise information on atmospheric clouds are not exploited. One of the reasons for that is the absence of computationally efficient observation operators. To remedy this issue in anticipation of the future regional Kilometer-Scale Ensemble Data Assimilation (KENDA) system of Deutscher Wetterdienst, we have developed a version that is fast enough for investigating the assimilation of cloudy reflectances in a case study approach. The operator solves the radiative transfer equation to simulate visible and near-infrared channels of satellite instruments based on the one-dimensional (1D) discrete ordinate method. As input, model output of the operational limited-area Consortium for Small-Scale Modeling (COSMO) model of Deutscher Wetterdienst is used. Assumptions concerning subgrid-scale processes, calculation of in-cloud values of liquid water content, ice water content, and cloud microphysics are summarized, and the accuracy of the 1D simulation is estimated through comparison with three-dimensional (3D) Monte Carlo solver results. In addition, the effects of a parallax correction and horizontal smoothing are quantified. The relative difference between the 1D simulation in “independent column approximation” and the 3D calculation is typically less than 9% between 0600 and 1500 UTC, computed from four scenes during one day (with local noon at 1115 UTC). The parallax-corrected version reduces the deviation to less than 6% for reflectance observations with a central wavelength of 810 nm. Horizontal averaging can further reduce the error of the 1D simulation. In all cases, the bias is less than 1% for the model domain.