Sunlight drives the Earth's weather, climate, chemistry, and biosphere. Recent efforts to improve solar heating codes in climate models focused on more accurate treatment of the absorption spectrum or fractional clouds. A mostly forgotten assumption in climate models is that of a flat Earth atmosphere. Spherical atmospheres intercept 2.5 W⋅m-2 more sunlight and heat the climate by an additional 1.5 W⋅m-2 globally. Such a systematic shift, being comparable to the radiative forcing change from preindustrial to present, is likely to produce a discernible climate shift that would alter a model's skill in simulating current climate. Regional heating errors, particularly at high latitudes, are several times larger. Unlike flat atmospheres, constituents in a spherical atmosphere, such as clouds and aerosols, alter the total amount of energy received by the Earth. To calculate the net cooling of aerosols in a spherical framework, one must count the increases in both incident and reflected sunlight, thus reducing the aerosol effect by 10 to 14% relative to using just the increase in reflected. Simple fixes to the current flat Earth climate models can correct much of this oversight, although some inconsistencies will remain.

In calculating solar radiation, climate models make many simplifications, in part to reduce computational cost and enable climate modeling, and in part from lack of understanding of critical atmospheric information. Whether known errors or unknown errors, the community's concern is how these could impact the modeled climate. The simplifications are well known and most have published studies evaluating them, but with individual studies it is difficult to compare. Here, we collect a wide range of such simplifications in either radiative transfer modeling or atmospheric conditions and assess potential errors within a consistent framework on climate-relevant scales. We build benchmarking capability around a solar heating code (Solar-J) that doubles as a photolysis code for chemistry and can be readily adapted to consider other errors and uncertainties. The broad classes here include: use of broad wavelength bands to integrate over spectral features; scattering approximations that alter phase function and optical depths for clouds and gases; uncertainty in ice-cloud optics; treatment of fractional cloud cover including overlap; and variability of ocean surface albedo. We geographically map the errors in W m−2 using a full climate re-creation for January 2015 from a weather forecasting model. For many approximations assessed here, mean errors are ∼2 W m−2 with greater latitudinal biases and are likely to affect a model's ability to match the current climate state. Combining this work with previous studies, we make priority recommendations for fixing these simplifications based on both the magnitude of error and the ease or computational cost of the fix.