Analyses of global surface-temperature change have been routinely carried out by several groups, including the NASA Goddard Institute for Space Studies (NASA GISS), the NOAA National Climatic Data Center (NCDC), and a joint effort of the UKMet Office Hadley Centre and the University of East Anglia Climatic Research Unit (HadCRU) . The record of the anomaly of this global temperature with respect to the climatological normal (CLINO) 1961-1990 for 1850 to 2018, the so-called HadCRUT4 data, and the Berkeley record  for the same period as well as the NASA GISS data for 1880 to 2018 are shown in Figure 1. The latter two temperature-anomaly records were originally expressed with respect to the CLINO 1951-1980, but for comparison we related them to the CLINO 1961-1990 as well which has been retained as a standard reference period for long-term climate change assessments . The results of the first GISS analyses of global surface-temperature change were published in 1981 . The first series of the HadCRU observational surface-temperature data set is related to Jones . Since Jones et al.  found for this CLINO 1961-1990 a global mean surface air temperature of ( for the Southern Hemisphere, SH, and for the Northern Hemisphere, NH),we can estimate the global mean temperature for any climatological period (usually, at least, 30 years). The HadCRUT4 records, for instance, provided ( and ) for 1851-1880 and
( and ) for 1881-1910. The Berkeley record provided and
for these periods, respectively. The NASA GISS records yielded
( and ) for 1881-1910. When we considered 1991-2018, the HadCrut4 re-
cord yielded ( and ), the Berkeley record provided
, andthe NASA GISS record provides ( and ). Thus, these results suggest an increase in the globally averaged near-surface temperature during the past 100 years of 0.8˚C, 1.0˚C, and 0.9˚C, respectively.
In 1900, von Bezold  (who became famous due to his five papers on the thermodynamics of the atmosphere and especially by introducing the potential temperature into the literature in 1888 ) estimated the climatological means of solar radiation (in German: “Sonnenstrahlung”), air temperature (“Lufttemperatur”), air pressure (“Luftdruck”), cloudiness (“Bewölkung”), and precipitation (“Niederschlag”) for numerous parallels of latitude, so-called zonal averages (see Figure 2 and Figure 5). His estimates are based on: 1) Meech’s  solar-radiation data, 2) the temperature data of Spitaler  (referred to as vB1901-1) and Batchelder  (referred to as vB1901-2), 3) Ferrel’s  air-pressure data, 4) Murray’s  precipitation data, and 5) Arrhenius’  cloudiness data that are based on Teisserenc de Bort’s chart of cloudiness distribution. All these data are listed in von Hann’s (ennobled by Austrian emperor Franz Joseph in 1910) first of his three volumes of the Handbook of Climatology [15,16]. This epoch-making work on general and regional climatology included data and eyewitness descriptions of weather and climate . The first edition was already published in 1883 . Von Hann [15,16] also reported the results of other climate researchers like Dove , Schoch , Satorius von Waltershausen ,
Figure 1. Time series of the annual temperature anomaly with respect to the climate normal 1961-1990 represented by the HADCRUT4 data for 1850 to 2018 (Morice et al. , but updated). Also shown are the Berkeley record for the same period (Rhode et al. , but updated), and the NASA GISS data for 1880 to 2018 (Hansen et al. , but updated).
Figure 2. Von Bezold’s  table (digitized by Google) of climatological means of solar radiation (“Sonnenstrahlung” in thermal days ), temperature (“Temperatur” in ˚C), air pressure (“Luft- druck” in mm Hg), precipitation (“Niederschlag” in cm), and cloudiness (“Bewölkung” in %) for numerous parallels of latitude (“Geographische Breite” in ˚). Values in parenthesis are uncertain. The normal-temperature-datasets are referred to as vB1901-1 (“Spitaler”) and vB1901-2 (“Bat- chelder”).
Ferrel , and Spitaler . Their results are listed in Table 1. Thus, von Hann [15,16] concluded that the average temperature for the Earth as a whole is about ,where both hemispheres have nearly the same temperatures, but the Southern Hemisphere would probably be slightly cooler because heat is carried across the equator by the strong southeast trade drift. Note that the original datasets of Dove , Ferrel , and Spitaler  used in our study are referred to as Do1852, Fe1877, and Sp1885, respectively (Table 2).
A slight improvement of Dove’s data was derived by Forbes  referred to as Fo1859 (Table 2). Furthermore, Spitaler  used the isothermal charts of Wild  and von Hann . Whereas Batchelder  and Arrhenius [14,26] considered Buchan’s  isothermal charts which are based on the observations of the “Challenger” Expedition. The datasets of Batchelder and Arrhenius are referred to as Ba1894, A1896-1, and A1896-2, respectively (Table 2). Spitaler , Batchelder , and Arrhenius [14,26] considered the 5th and the 10th degree of latitude, respectively. In contrast to this, von Bezold  used an equidistant distribution in ,where is the latitude. Thus, von Bezold’s datasets differ slightly from those of the original ones due to his interpolation procedure.
Von Bezold  also estimated the annual mean temperature for the Earth and confirmed von Hann’s [15,16] value of . Von Bezold’s result is based on the table shown in Figure 2. In his presentation of von Bezold’s  paper, however, von Hann  argued that Batchelder’s  temperature values for the Tropics were undoubtedly too high because the isotherms of the “Challenger” reports were based on too high temperature averages (locally 1˚C - 2˚C). Von Hann  suggested an average temperature for the Earth of . Von Bezold’s  table was eventually adopted by Börnstein , who only considered Spitaler’s  meridional distribution of zonal averages of temperature. Note that the zonal averages of temperature are also called the normal temperatures [11,18].
In 1906, von Bezold  presented slightly improved results for the climatological means of the near-surface air temperature along numerous parallels of latitude (referred to as vB1906, Figure 3). Again, von Bezold estimated the annual mean temperature for the Earth as .
Von Hann  re-examined the climatological mean temperatures along numerous parallels of latitude in the first volume of the third edition of his handbook. He not only considered the normal-tempe- rature-datasets of Spitaler  and Batchelder , but also those of Hopfner . Hopfner followed Dove’s  epoch-making work, but he considered Buchan’s  isothermal charts available for each month of the year and the whole year, respectively. Based on these isothermal charts (already used by Batchelder
Table 1. Average temperature for the Northern Hemisphere (NH), Southern Hemisphere (SH), and the Earth as reported by von Hann ([15,16]).
*) Based on Equation (9).
Table 2. Normal-temperature-datasets for numerous parallels of latitude as published by Dove , Forbes , Ferrel , Spitaler , Batchelder , Arrhenius [14,26], Hopfner , Defant and Obst , von Hann-Süring , and Sellers . Bold numbers denote the Thermal Equator.
*) Re-corrected to mean terrain height above sea level.
 and Arrhenius [14,26]), Hopfner  derived two slightly different normal-temperature-datasets for numerous parallels of latitude: 1) the twelve-monthly means and 2) directly taken from the chart of annual isotherms. These two normal-temperature-datasets are referred to as Ho1906-1 and Ho1906-2 (Table 2). Figure 4 shows von Hann’s normal-temperature dataset (referred to as vH1908). Based on these data, he suggested (see also Lockyer ). Börnstein  also re-examined the climatological mean values of solar radiation, air temperature, air pressure, precipitation, and cloudiness along numerous parallels of latitude (Figure 5). Börnstein’s normal-temperature-dataset, referred to as Bö1913, completely agrees with von Bezold’s  vB1906.
Obviously, the results of well-known climate researchers for the average temperature of the near- surface air temperature published during the second half of the 19th century and the first two decades of the 20th century are notably higher than those derived from the HadCRUT4, Berkeley, and NASA GISS
Figure 3. Von Bezold’s  table (digitized by Google) of mean annual and mean monthly (January and July) temperatures for numerous parallels of latitude (characterized by the sine of the respective angle) of both the Northern Hemisphere (“Nördliche Halbkugel”) and the Southern Hemisphere (“südliche Halbkugel”). Temperatures are in ˚C. Values in parenthesis are uncertain. The normal-temperature-dataset is referred to as vB1906.Note that von Bezold considered ,where is the latitude.
Figure 4. Von Hann’s  table (digitized by Google) of annual and monthly (January, April, July, and October) mean temperatures for numerous parallels of latitude with respect to Spitaler , Batchelder , and Hopfner . Temperatures are in ˚C. Land coverage (“Landbedeckung”) is in %. Values in parenthesis are uncertain. This dataset is referred to as vH1908.
records. The goal our paper is, therefore, to assess the results of von Hann and his fellow climate researchers. Since Defant and Obst , Köppen , von Hann-Süring , and Sellers  also published meridional distributions of the normal temperatures, we used their datasets for additional assessment. These normal-temperature datasets are referred to as De1923, K1936, vHS1939, and Se1965, respectively. Eventually, Haurwitz and Austin  and Blüthgen  adopted vHS1939.
Since three different temperature scales related to de Réaumur (e.g., Dove), Celsius (e.g., von Hann), and Fahrenheit (e.g., Buchan) were used, we used the original sources to check all datasets adopted and converted by others. No serious conversion error was detected.
2. The Globally Averaged Surface Temperature
The average over the Earth’ surface reads [42 - 44]
Here, is an arbitrary variable, is the radius, is the solid angle of the entire planet, and is the differential solid angle, where and are the zenith and azimuthal angles, respectively, in a spherical coordinate frame (Figure 6). Note that ranges from 0 to π, and ranges from 0 to 2π.
Figure 5. Börnstein’s  table (digitized by Google) of climatological means of solar radiation (“Sonnenstrahlung” in thermal days ), temperature (“Temperatur” in ˚C), air pressure (“Luft- druck” in mm Hg), precipitation (“Niederschlag” in mm), and cloudiness (“Bewölkung” in %) for numerous parallels of latitude (“Geographische Breite” in ˚).Values in parenthesis are uncertain. The dataset of the normal temperatures is referred to as Bö1913.
2.1. Spherical Averaging
For a spherical shape of the Earth, we have and, hence, . Usually, is represented by the volumetric mean radius of . Thus, Equation (1) can be written as [10,42,44,46]
Since the average along a parallel of latitude, i.e., the zonal average, is defined by [44,46 - 49]
Equation (2) can be written as
Inserting the averages for both the northern and the southern hemispheres given by 
Figure 6. Mathematical representation of the solid angle. Here, is the differential solid angle, where and are the zenith and azimuthal angles, respectively (adopted from Kasten and Raschke ).
When we set, for instance, ,we obtain the globally averaged near-surface air temperature
Since the angle of the latitude, ,and the zenith angle, ,are related to each other by ,Equation (8) can be written as
2.2. Spheroidal Averaging
The Earth is, of course, not a sphere, but in a first approximation an oblate ellipsoid of revolution, i.e. an oblate spheroid, flattened at the poles and bulging at the Equator due to centrifugal forces. This oblate spheroid serves as a surface of reference for the mathematical reduction of geodetic and cartographic data . Its parameters are the semi-major axis ,the semi-minor axis and the flattening according to the World Geodetic System 1984 (WGS84) . With respect to a Cartesian frame with its origin in the center of mass and its vertical axis congruent with Earth’s rotation axis, we have
The coordinates are given by and . The ellipsoidal radius is
where is the difference between and of the Earth. The ratio ranges from at the Poles to at the Equator (Figure 7). Using Equation (13) leads to
Since the surface of an oblate spheroid is also given by 
Figure 7. The ratio as a function of the zenith angle .
or in the re-arranged form
with ,we have
The equivalent spherical radius is . In the case of a sphere, the rule is exactly fulfilled, where A and V are the surface and the volume, respectively. Since of all surfaces of the same area, the sphere has the greatest volume, we have . Consequently, is marginally larger than ,but this difference is negligible. Thus, we may write
Using again the zonal average defined by Equation (3) leads to
Inserting the averages for both the Northern and Southern Hemispheres defined by
Thus, the globally averaged near-surface air temperature is given by
As in the case of spherical averaging, can be substituted by .
3. Discussion of the Historic Data
3.1. The Problem of the Incomplete Spatial Coverage by Meteorological Stations
Hansen and Lebedeff  already discussed the problem of the incomplete spatial coverage of the Earth’s surface by the network of meteorological stations. They correctly stated:
“The principal limitation of this data set for global or hemispheric analysis is the incomplete spatial coverage, illustrated in Figure 1 [here, Figure 8] for four dates. Although the number and geographical extent of recording stations on land areas increased strongly between 1870 and 1900, there were still large areas in Africa and South America, and all of Antarctica, without coverage in 1900. Substantial station data for Antarctica begins in the 1950s. Large ocean areas remain without fixed meteorological stations at all times.”
However, as expressed by Equations (8) and (9), the uncertainty regarding observations in the south-polar region is of minor importance in determining the globally averaged near-surface temperature. Furthermore, von Hann and his fellow climate researchers not only gathered data from meteorological stations, but also eyewitness descriptions of weather and climate. Over land, the near-surface dry- and wet-bulb temperatures as well as the maximum and minimum temperatures were measured in thermometer boxes (Stevenson’s screens) at a height of 2 m or so above the ground using a standardized thermometer equipment as illustrated in Figure 9. The standardization of land stations goes back to the First Meteorological Congress held in Vienna in September 1873. Such an equipment suitable for sea-going conditions was also used aboard the H.M.S. Challenger during her 1873 to 1876 voyage , where either hourly or two-hourly meteorological observations were performed.
Measurements of the sea-surface temperature have been performed since the mid-18th century; mainly wooden and canvas buckets served to take water probes. Krümmel  already discussed the accuracy of methods to measure the sea-surface and deep-sea temperatures. During the voyage of H.M.S. Challenger, the sea-surface temperature was observed every two hours .
Figure 8. Global distribution of meteorological stations with near-surface air temperature records for 1870, 1900, 1930 and 1960. A circle of 1200-km radius is drawn around each station (adopted from Hansen and Lebedeff ).
Figure 9. Typical thermometer equipment for measuring the dry-bulb (T, “trocken”) and wet-bulb (F, “feucht”) temperatures as well as the maximum (Ma) and minimum (Mi) temperatures in thermometer boxes (Stevenson’s screens) over land (adopted from Börnstein ).
Nonetheless, even today, network density and design affect the regional averages over landscapes. In 2009, PaiMazumber and Mölders  used the results of their WRF simulations performed over Russia for July and December 2005, 2006, and 2007 to create a ‘‘dataset’’ for assessing this impact theoretically. Based on the values at all WRF grid points, they calculated regional averages for various quantities for 2.8˚ × 2.8˚ areas as the ‘‘reference.’’ Then, regional averages determined based on 40 artificial networks and the actual 411 ‘‘sites’’ locations of a real network were compared with the reference regional averages. The 40 networks encompassed 10 networks of 500, 400, 200, or 100 different randomly taken WRF grid points as sites. The results showed that the real network’s site distribution misrepresents the landscape. Such misrepresentation leads to differences in regional averages that show geographical and temporal trends for most quantities: Errors are lower over shores of large lakes than coasts and lowest over flatland followed by low and high mountain ranges; offsets in timing occur during frontal passages when several sites are passed at nearly the same time. Generally, the real network underestimated regional averages of sea-level pressure, wind speed, and precipitation over Russia up to 4.8 hPa (4.8 hPa), 0.7 m∙s−1 (0.5 m∙s−1), and 0.2 mm day−1 and overestimates regional averages of 2 m temperature, downward shortwave radiation, and soil temperature up to 1.9 K (1.4 K), 19 W∙m−2 (14 W∙m−2), and 1.5 K (1.8 K) in July (December) .
3.2. Historical Isothermal Charts
Alexander von Humboldt was the first who drew annual isotherms on a map of the Earth in 1817 . He used the temperature averages of 58 locations available around the world . Based on the course of these isotherms, von Humboldt deduced the simplest laws of the heat distribution on the Earth’s surface. In his blog contribution “The first isothermal world maps”, Mike Klein, Senior Cartographic Specialist in the Geography and Map Division at the Library of Congress, stated:
“Humboldt withheld publishing his idea in the form of a generalized world map while he waited for data from weather stations around the world. Thus, the isotherm concept remained virtually unknown outside of the broader scientific community until 1838, when it was given wider recognition on the map below [here Figure 10]. Titled Alexander von Humboldt’s System der Isotherm-Kurven in Merkator’s Projection, it was published in Germany by Heinrich Berghaus in his Physikalischer Atlas, the first comprehensive physical world atlas.”
Figure 10 shows this isothermal chart published by Heinrich Berghaus in the 1845 Edition of his Physikalischer Atlas. (Note that the correct spelling is Mercator, but not Merkator as used in von Humboldt’s isothermal chart.)
Figure 10. Alexander von Humboldt’s system of the isotherm curves in Mercator’s Projection. It was published in Germany by Heinrich Berghaus in his 1845 Edition of the Physikalischer Atlas (adopted from https://blogs.loc.gov/maps/2018/04/the-first-isothermic-world-maps/). Note that von Humboldt used both the Celsius scale (left side of the chart) and the de Réaumur scale (right side of the chart).
Batchelder  reviewed earlier observations and the use of isothermal charts. By quoting von Hann , he underlined the importance to distinguish between water and land masses:
“Durch den Umstand, dass ein Parallelkreis teils über Land, teils über Wasser verläuft, entstehen klimatische Unterschiede zwischen West und Ost, oder Verschiedenheiten des Klimas nach den Meridianen, welche im solaren Klima nicht vorhanden wären. Neben der ungleichen Erwärmung und Erkaltung von Wasser und Land werden ausserdem durch das Vorhandensein des Landes gewisse konstante Luft- und Meeresströmungen erzeugt, welche ebenfalls eine Verschiedenheit des Klimas unter verschiedenen Meridianen desselben Parallels bedingen.”
The translation from German into English given by Batchelder  reads:
“Owing to the fact that a circle of latitude passes partly over land and partly over water, there arise climatic differences between West and East, or variations of the climate according to the longitudes, that in the solar climate would not exist. Following the dissimilar warming and cooling of water and land there are, moreover, set up, owing to the presence of the land, certain constant air and sea currents which also make possible a difference of climate on different meridians of the same parallel.”
The notion “solar climate” is discussed in Subsection 3.3.
Von Hann  and Spitaler  also discussed Forbes’  formula for determining the average (or normal) temperature of a parallel of latitude derived for the Northern Hemisphere
where n is the fraction of land compared to the circumference of the respective parallel. The coefficients A, B, C, and the exponent m must be derived from observations. Forbes  derived (here expressed with respect to the Celsius scale): A = −10.8˚C, B = 32.9˚C, C = 21.2˚C, and m = 1.25. Forbes also discussed the temperature difference between an aqueous globe (n = 0) and a rocky one (n = 1). Figure 11 shows the meridional distributions of the zonal averages of temperature, ,between 50˚N and 50˚S. Forbes found for the latitude of 45˚ a surface temperature of for both the aqueous globe and the rocky one. For the Equator, he obtained for n = 0 and for n = 1. However, it is unlikely that in both cases the effect of the atmosphere would be the same. With respect to this result, von Hann  stated:
“Between latitude 40˚ and latitude 50˚ is the parallel on which both water and land hemispheres have the same temperature. In higher latitudes, the water hemisphere is warmer than the land hemisphere. Even Dove’s isothermal charts of the northern hemisphere show that the transition takes place between latitudes 40˚ and 45˚, and the graphic representation, above referred to, shows more precisely that this transition occurs at latitude 42.5˚. Beyond this latitude, as far as the pole, a water hemisphere is warmer than a hemisphere wholly covered by land.”
In his paper on von Hann’s contribution to modern climatology, Kahlig  reported that von Hann used a three-term moving average (running mean) procedure with unequal weights (1/4, 1/2, 1/4), called “hanning” (e.g., ), to smooth the raw data of land fraction listed in von Hann’s [15,16] table of the “Mean temperatures of the Parallels of Latitude” in his Handbook of Climatology (pp. 199-200). The raw data and the smoothed ones are listed in columns a and b in von Hann’s table (Figure 4). Thus, von Hann was able to handle Forbes’  formula (Equation (24)) in a sophisticated manner.
Isothermal charts as developed by Dove , Wild , von Hann  and Buchan  played an important role in deriving the normal temperatures for numerous parallels of latitude. Schoch  used Dove’s  isothermal charts, Spitaler  used those of Wild  and von Hann  (eventually published in Berghaus’ Physikalischer Atlas), Batchelder , Arrhenius , and Hopfner  used Buchan’s  isothermal charts. The isothermal charts shown in Figures 12-16 were published by von Hann
Figure 11. Meridional distributions of the zonal averages of temperature, ,between 50˚N and 50˚S for an aqueous globe (n = 0) and a rocky one (n = 1) based on Forbes’  formula (Equation (24)). Dots represent dataset Do1852  (see Table 2).
Figure 12. Von Hann’s  annual isotherms of the Earth published in Berghaus’ Physikalischer Atlasin 1887. Temperatures are in ˚C.
Figure 13. Buchan’s  annual isotherms of the Earth published in the Challenger reports in 1889. Temperatures are in ˚F.
Figure 14. As in Figure 13, but prepared by Bartholomew and Herbertson and edited by Alexander Buchan  (adopted from the ETH-Bibliothek Zürich https://doi.org/10.3931/e-rara-71030). Temperatures are in ˚F.
Figure 15. As in Figure 12, but published by von Hann  in his “Lehrbuch der Meteorologie” (Textbook on Meteorology). Temperatures are in ˚C.
Figure 16. Börnstein’s [31,35] annual isotherms of the Earth (digitized by Google).Temperatures are in ˚C. The dashed lines indicate the 10˚-isotherms of the warmest month.
 in 1887, Buchan [27,60] in 1889 and 1899, von Hann  in 1906, and Börnstein [31,35] in 1906/ 1913, respectively. Dove  used monthly means, seasonal means and annual means determined from the longest so far accessible series of observations at 900 stations. As reported by von Hann , his isothermal chart is based on all information available to him in January 1884. This information also included Wild’s  isothermal charts for the Russian Empire. Buchan’s  isothermal chart not only included the observations of the Challenger Expedition, but also the observations performed at 1620 stations during the fifteen-year period 1870-1884, except in the United States, where the observations started in October 1871, when the Signal Service of the War Department took charge of the Meteorological System of the United States. In his report “Summary of International Meteorological Observations”, Dunwoody  stated, for instance (Figure 17):
Figure 17. Excerpt from Dunwoody’s  “Summary of international meteorological observations”.
DeCourthy Ward  assessed von Hann’s  and Buchan’s  isothermal charts together with those prepared by Bartholomew and Herbertson and edited by Buchan  for the comprehensive work “Atlas of Meteorology: a series of over four hundred maps”. Figure 14 shows the isothermal chart of annual means. It is a reproduction of Buchan’s isothermal chart that appeared in the Challenger reports, but revised by him . Figure 14 also shows Supan’s  temperature zones and climate provinces, and Dunwoody’s  map of annual simultaneous observations. The isothermal charts of von Hann  (Figure 15) and Börnstein [31,35] (Figure 16) slightly differ from each other. Unfortunately, the isotherms were drawn at different intervals. Only the 0˚C - and the 20˚C-isotherms are directly comparable.
More than 85 years later, Levitus  published the Climatological Atlas of the World Ocean. Figure 18 shows his global distribution of the annual mean potential temperature at the sea surface. Figure 19 shows the corresponding figure based on the Ocean Atlas 2018. These figures may be used to assess the distribution of the sea surface in the historical isothermal charts.
Figure 18. Levitus’  distribution of the annual potential temperature at the sea surface in ˚C (but taken from Peixoto and Oort ).
Figure 19. World Ocean Atlas Climatology (decadal average 1955-2017): Distribution of the annual sea-surface temperature in ˚C (one-degree grid) .
Spitaler  described his procedure to determine the normal temperature along the parallels of latitude as follows:
“Die Grundlage für die vorliegende Untersuchung lieferten die neuen Isothermenkarten von Wild und Prof. Hann, welche nach dem sämtlichen bis jetzt vorliegenden Beobachtungsmaterial der Erde gezeichnet wurden. Ich habe für jeden 10. Breitengrad von 5 zu 5 Längengraden, für die dazwischenliegenden Breitengrade aber nur für jeden 10. Längengrad die Temperatur des Jahresmittels bestimmt, sowie die mittleren Temperaturen der beiden extremen Monate Januar and Juli graphisch interpoliert und auf diese Weise einerseits aus je 72, andererseits aus je 36 äquidistanten Temperaturwerten die normale Temperatur der Breitenkreise bestimmt.”
The translation reads:
“The basis for the present study was the new isothermal charts of Wild and Prof. Hann, which were drawn based on all the Earth’s observation material available so far. For each 10th degree of latitude and for 5 to 5 degrees of longitude, but for the intermediate latitudes only for each 10th degree of longitude, I determined the annual-average temperature and interpolated the mean temperatures of the two extreme months of January and July graphically. In doing so, I determined the normal temperature of the parallels of latitude based on 72 and 36 equidistant temperature values, respectively.”
Batchelder  determined the normal temperature of each 10th degree of latitude, from 80˚N to 50˚S, by averaging the observed temperatures at 36 equidistant points or stations on the circle, i.e., at every 10thmeridian. The temperatures at these 36 points were determined, when an isotherm did not happen to fall directly over the point, by interpolation between the nearest isotherms. Referring to Forbes , he stated:
“Although one authority considers that ‘It is by no means an easy matter to deduce the mean temperature of a given parallel correctly from an isothermal chart,’ yet as Buchan constructed his isotherms at intervals of 5˚F, interpolation was usually simple and direct, and the limit of error rarely rose above one degree. The most unsatisfactory cases of interpolation were those succession of isotherms) as in July, 130˚W, 50˚N. Still, even here, 1˚F seems usually the error. A peculiar case occurs in January over the northern Amazon basin; as the observer moves north ward over South America he finds the temperature fall from the central part until above the equator, where occurs an isolated area of +85˚ over the Isthmus of Panama. The neighborhood of this area is very difficult to interpolate for; to the north the nearest isotherm is that of 75˚, a ten-degree skip, to the east and west are no neighboring isotherms, and far away to the south is the isotherm of 80˚ on a falling gradient. In the case of the parallel of 80˚N on the Annual Chart, too, the limit of error between 100˚W and 130˚W may rise as high as 4˚F. On the equator, owing to the slight gradient, the general limit of error is about 2˚F.”
Hopfner  also determined this normal temperature of each 10th degree of latitude, but from 80˚N to 60˚S by averaging the observed temperatures at 36 equidistant points or stations on the respective parallel of latitude. He also discussed the accuracy of the use of isothermal charts in detail.
As argued by Köppen , notable variation exists in the annual mean temperatures of the lower atmosphere, reduced to sea level, for numerous parallels of latitude (cf. Table 3).
As listed in Table 3, the temperature difference between the Equator and 80˚N is 38˚C (maximum) and 44˚C (minimum). Whereas the difference between the Equator and 80˚S is 42˚C (maximum) and 48˚C (minimum) . Thus, when computing the near-surface temperature of the Earth, we must address possible uncertainty related to such variations in annual averages.
The reduction of temperatures to sea level seems to be a problem. Even though the reduction of air pressure to sea level is widely accepted, this is not the case of near-surface air temperature. In his textbook “Lehrbuch der Meteorologie“ published in 1906, von Hann  stated:
Table 3. Annual averages of the temperature of the lower atmosphere, reduced to sea level, for numerous parallels of latitude (adopted from Köppen ).
“Die Isothermen stellen die Temperaturverteilung auf der Erde so dar, wie wenn alle Orte, nach deren Temperaturaufzeichnungen sie entworfen worden sind, im Meeresniveau liegen würden. Die Temperaturmittel müssen daher, bevor man sie auf der Karte einträgt, auf das Meeresniveau reduziert werden, wozu die früher angeführten Erfahrungen über die Wärmeabnahme mit der Höhe benützt werden.”
The translation reads:
“The isotherms represent the temperature distribution on Earth as if all the locations, which contributed to them with their temperature records, were at sea level. The temperature averages must, therefore, be reduced to sea level before entering the chart to which the previously mentioned experience about the heat decrease with height are used.”
Buchan  argued (his temperature notation has been updated):
“The mean temperature falls about 5.5˚C for every 1000 metres of ascent (1˚F per 330 feet or 16˚F per mile); but, in the Northern Hemisphere, between 30˚N and 70˚N, the mean temperature decreases only 0.75˚C (1.35˚F) for every increase of a degree of latitude; i.e. 0.0068˚C for every kilometre nearer the pole (0.0197˚F per mile). The temperature therefore falls 800 times more in the same distance in a vertical direction than in a horizontal one. This may be expressed by saying that the vertical temperature gradient is 800 times steeper than the horizontal one. This contrast can be seen in the sections on Plate 1, where the vertical distances are much exaggerated when compared with the horizontal ones, and hence the sections greatly minimise the contrast.”
Note that von Hann  argued in a similar manner. Based on the CRU dataset, Feulner et al.  recently illustrated the effect of station’s elevation. They argued that the effect of land (and ice sheet) elevation combined with the lapse rate on the temperature distribution is apparent, making Antarctica, Greenland, and the large mountain ranges markedly cooler than their surroundings.
Sellers  did not apply such reduction to sea level in determining his mean annual surface temperatures. He stated:
“Temperatures given in [his] Table 1 are estimated mean annual surface temperatures and are lower than the sea level values often quoted in the meteorological literature (Haurwitz and Austin, 1944 ; Trewartha, 1954 ), especially between 20˚N and 50˚N and south of 60˚S. The average global surface temperature is about 286˚K, compared with 288˚K for the average global sea level temperature.”
We considered Arrhenius’ [14,26] data re-corrected to mean terrain height above sea level and Sellers’  data to illustrate the difference in the global average of the near-surface air temperature caused by zonal averages of temperatures at site elevation and those reduced to sea level.
Sellers considered the physical and climatological data of 10˚-latitude zones. Thus, we assigned his temperature data to the center of the respective latitude zones. Figure 20 reveals that Sellers’zonal temperature means are remarkably lower than those of the other authors. At latitudes beyond 60˚S, they are even lower than those in Köppen’s  distribution of zonal mean minimum temperatures. Beyond 30˚S, Ferrel’s  zonal temperature means slightly exceed Köppen’s.
The normal temperatures determined by the authors considered in this study are not completely independent. As aforementioned, subsequent authors improved them using the rapidly increased number of weather observations in the 2nd half of the 19th century (Figure 8). The differences between Dove’s  and Forbes’  values were still marginal, but Wild , von Hann  and Buchan [27,60] used augmented weather information for editing their isothermal charts.
Graphic solutions have a long tradition in science and technology. Cremona-Maxwell diagrams, for instance, were used in statics with sufficient accuracy long before Konrad Zuse developed the first programmable computer in the late 1930s.
3.3. Solar Climate versus Real Climate
Forbes’  formula (see Equation (24)) of the asymmetric water and land distribution may be considered as the first attempt to explain the difference between the isotherms of the Earth’s real climate and those of its so-called solar climate. According to von Hann [16,18], the solar climate, sometimes also called the mathematical climate, refers to the thought model of an Earth in the absence of its atmosphere. Von Hann  stated
“If the surface of the earth were occupied altogether by land, and if there were no surrounding atmosphere, the condition of our planet would be somewhat similar to that of the moon at the present time. Under these conditions, the distribution of temperature over the earth would depend solely upon the amount of heat received from the sun at any given place, and upon the loss of heat by radiation at that place. As these two factors would necessarily be the same at all points along the same parallel of latitude, the zones of equal temperature would coincide with the parallels of latitude. Even the presence of a vapourless atmosphere would interfere but little with this
Figure 20. Comparison of the meridional distributions of zonal averages of temperature, ,as derived by Dove  (Do1852), Forbes  (Fo1859), Ferrel  (Fe1877), Spitaler  (Sp1885), Batchelder  (Ba1894), Arrhenius  (A1896-1), von Bezold  (vB1901-1), Hopfner  (Ho1906-1), von Hann  (vH1908), von Bezold  (vB1906), Börnstein  (Bö1913), Defant and Obst  (De1923), Köppen , von Hann-Süring  (vHS1939, eventually adopted by Haurwitz and Austin  and Blüthgen ), and Sellers  (Se1965).
distribution of temperature, for only the absolute amounts of heat received at, and radiated from, the surface of the earth would thereby be affected. It is true that convectional currents would be produced under these conditions; but as there would be no reason for the more frequent occurrence of warm or cold air currents along some meridians than others, the distribution of temperature in zones bounded by the parallels of latitude would not thereby be interfered with.”
The solar climate depends on the astronomic aspects like the distance between the Sun’s center and the Earth (Earth-Moon barycenter), the obliquity of Earth’s rotation axis with respect to the normal of the ecliptic plane, the angular velocity of the rotation, and the total solar irradiance (TSI) for 1 AU. Only the absorbed solar radiation and the thermal effect related to the regolith should be considered. With respect to geological time scales, also changes in the precession of Earth’s rotation axis, and long-term variations of the eccentricity, obliquity, and precession of the Perihelion caused by the Sun, Moon and planets of our solar system must be considered [69 - 75].
Figure 21 compares the meridional distributions of the zonal temperature averages of the real climate as given by Spitaler  (Sp1885), Arrhenius  (A1896-1), von Bezold  (vB1906), as well as von Hann-Süring  (vHS1939) and the solar climate predicted by Kramm et al.  for the year 2010. Kramm et al.  used the following local energy budget equation for a thin slab of 2 cm thickness adjacent to the surface of either the Earth’s Moon or the Earth in the absence of its atmosphere:
where is the volume average.Here, t is time, , ,and c, are the temperature, bulk density, and specific heat of this slab of soil, respectively. Furthermore, is the solar radiation absorbed by the slab, where F is the solar irradiance reaching the surface, is the integral albedo inthe solar range, and is the local zenith distance of the Sun’s center. Kramm et al.  calculated by Keihm’s  empirical formula rearranged to
Figure 21. Meridional distributions of (a) the zonal averages of temperature, ,of the real climate, represented by Spitaler  (Sp1885), Arrhenius  (A1896-1), von Bezold  (vB1906), and von Hann-Süring  (vHS1939), the sea surface temperature given by Krümmel , and the solar climate derived by Kramm et al. , as well as (b) the difference, ,between the zonal averages of the real climate, the sea surface temperature only and those of the solar climate.
where is the normal albedo, and and are empirical values. With exception of Keihm’s value for b, all other values are based on observations from the Lunar Reconnaissance Orbiter Diviner Lunar Radiometer Experiment [77,78]. In the case of the Earth, was determined using the rules of spherical trigonometry (e.g., [74,75,79,80])
Here, is the declination of the Sun, is latitude, and h is the hour angle from the local meridian. The solar irradiance, F, is given by
where is the solar emittance [74,75,81,82],  is the radius of the Sun, and r is the actual distance between the Sun’s and Earth’s centers. Inserting the mean distance, (nearly 1 AU), into Equation (28) provides the solar constant for the Earth
Combining formulae (28) and (29) yields
Furthermore, is the emitted radiation according to the power law of Stefan
 and Boltzmann , where is the integral relative emissivity, and is the vertical component of the soil heat flux provided by a multi-layer soil model. The direction of is governed by the difference between the absorbed solar radiation and the emitted infrared radiation. Moreover, the planetary and lunar ephemeris DE430 of the Jet Propulsion Laboratory (JPL) of the California Institute of Technology [86,87] to compute r and . A TSI of 1361 W∙m−2 at 1 AU was taken from Kopp and Lean  and Kopp et al. . According to the TSI reconstruction by Kopp et al. , the TSI only changed by ±1 W∙m−2 during 1850 to 2015. Thus, the solar climate derived by Kramm et al.  may be compared with the real climates represented by the datasets Sp1885 , A1896-1 , vB1906 , and vHS1939 .
Figure 21 also shows Krümmel’s  meridional distribution of the zonal mean sea-surface temperatures (Figure 18 and Figure 19 may be considered to assess Krümmel’s data.). Obviously, between 45˚N and 45˚S, the meridional distributions of normal temperatures mainly follow those of the zonal mean sea surface temperatures. Already Forbes’s formula (Equation (24)) illustrated in Figure 11 documents this fact. Krümmel’s sea-surface data would provide for an aqueous globe. For comparison: Forbes’
formula (Equation (24)) would provide for an aqueous globe and for a
The slab-temperature distribution predicted by Kramm et al.  fulfills the condition of the global radiation balance. Based on the globally averaged solar radiation of and the globally absorbed solar radiation is ,the global albedo in the solar range is about . The globally emitted infrared radiation is . Thus, the radiative imbalance is . It is compensated by the globally averaged soil heat flux. As listed in Table 4, the
global average of this slab-temperature distribution is . Since ,one obtains for the obliquely rotating Earth in the absence of its atmosphere.
Figure 21(a) shows the meridional distributions of zonal averages of the solar irradiance, ,the absorbed solar irradiance, ,and various normal temperatures related to them for the Earth in the absence of its atmosphere. In contrast to Kramm et al. , Defant and Obst  used a solar constant of . Replacing their too high solar constant with the current one would provide a meridional distribution of that substantially agrees with that of Kramm et al. . Based on Equation (4), integration of from to provides for Defant and Obst’s original data and for the corrected data ,respectively. The integration of
from to yields, of course, ,mentioned before.
Defant and Obst  used the power law of Stefan  and Boltzmann  to derive the zonal average of a so-called radiation temperature, ,from . They assumed that the emitted radiation is equal to the completely absorbed solar irradiance. They compared the meridional distribution of this radiation temperature with that of the real zonal temperature means, ,and interpreted the difference as the heat protection due to the atmosphere. The meridional distributions of both
and are illustrated in Figure 21(b). Based on Equation (8), the integration of and from to provides and ,respectively.
This procedure, however, is incorrect from physical and mathematical points of view. Since the local radiation balance is given by 
where is the absorbed solar irradiance, and is the local surface temperature, zonal averaging would provide (see Equation (3))
In accord with the general inequality of Gerlich and Tscheuschner ,
for a non-negative measurable function and a probability measure W, the zonal average of the surface temperature, fulfills the inequality
In the case of completely absorbed solar irradiance, i.e., ,we would have and, hence, . Thus, as illustrated in Figure 22(b), any zonal average of the surface temperature, based on Equation (31) is usually much smaller than the corresponding
radiation temperature . Based on Equation (4), the integration of from to
provided . This global mean surface temperature substantially agrees with result of Gerlich and Tscheuschner  (cf. their Table 12), even though we considered an obliquely rotating Earth in the absence of its atmosphere. In the case ,we would have . Consequently, one must conclude that the solar climate as derived by Defant and Obst  has to be discarded.
As illustrated by Figure 23, for any parallel of latitude, the slab temperature that represents the Earth’s surface temperature in the absence of the atmosphere varies with time due to the Earth’s daily rotation and annual orbit around the Sun. However, this slab temperature is based on a local energy balance equation (25) that also included the soil heat flux density . The meridional distribution of the normal temperature along a parallel of latitude based on a local radiation balance would be the same as proposed
Figure 22. Meridional distributions of (a) the zonal averages of solar radiation, ,the absorbed solar radiation, ,and (b) various normal temperatures, ,for the solar climate, where stands for the radiation temperature of Defant and Obst , for the radiation temperature corrected to , for the radiation temperature related to the absorbed solar radiation, , for the slab temperature according to Kramm et al. , and for the local radiation temperature according to Gerlich and Tscheuschner . Also shown is the meridional distribution of the real normal temperatures, (data adopted from Defant and Obst  and Kramm et al. ).
(a) (b) (c) (d)
Figure 23. Diurnal variation of the slab temperature around (a) Perihelion; (b) Vernal Equinox; (c) Aphelion; and (d) Autumnal Equinox for the Earth in the absence of its atmosphere (data from Kramm et al. ).
by Gerlich and Tscheuschner  because during nighttime the surface temperature derived from the local radiation balance would be zero Kelvin (see Equation (31)). Therefore, as illustrated in Figure 22(b), this meridional distribution of normal temperature is, by far, the lowest one.
Obviously, the existence of the atmosphere, the nonuniformity of the Earth’s surface, the uneven ocean-land distribution, the airflows and ocean currents cause the transition from the solar climate to the real climate . The difference between the solar climate and the real one is usually called the atmospheric effect. Thus, this atmospheric effect would vary with the parallels of latitude.
The atmosphere is a complex thermo-fluid dynamic system with various degrees of freedom. It consists of various layers (e.g., troposphere, stratosphere, mesosphere, thermosphere) separated by conceptual partitions called pauses . It exhibits an impressive amount of detail and huge spatial and temporal variability of its properties . More than 99% of the atmospheric mass of about 5.15 × 1018 kg is below the altitude of 30 km above sea level. The atmospheric response time to an imposed change is of the order of days or weeks owing to its notable compressibility and its low density and specific heat . Consequently, the response time corresponds to weather periods, rather than climate periods of, at least, 30 years. According to Fortak  and Peixoto and Oort , the total energy of the atmosphere is about 1.31 × 1024 J, i.e., less than 0.1% of the total energy of the oceans.
Based on recent observations, the solar insolation at the top of the atmosphere (TOA) corresponds to a globally averaged value of about (100 units). As sketched in Figure 24, only about is absorbed by the entire Earth-atmosphere system (EAS). Compared with of the solar climate of Kramm et al. , the solar radiation absorbed by the EAS is notably lower. The absorbed solar radiation of is approximately balanced, on global average, by the infrared radiation, emitted to space by the entire EAS (Figure 24). Note that the TOA may be interpreted as a height of the intervening atmospheric layer. Above this height, neither solar radiation nor infrared radiation is remarkably affected by atmospheric constituents.
Figure 24. Sketch of the global energy budget of the system Earth-atmosphere (adopted from Mölders and Kramm ). Note that the values of the units slightly vary from author to author. These values just intend to assess the processes relative to each other.
As sketched in Figure 24, a notable portion of the solar insolation penetrating into the atmosphere is absorbed in the ultraviolet and visible ranges as well as in the near infrared range by various gaseous and particulate atmospheric constituents (20 units or 68.1 W∙m−2) and clouds (4 units or 13.6 W∙m−2). Especially the absorption of solar radiation by molecular oxygen (O2) and ozone (O3) heats the atmosphere directly [75,94]. Water vapor (H2O) and O2 are also active in the visible and near infrared range; nitrogen dioxide (NO2) is active in the visible range, too. Furthermore, a considerable portion of the solar radiation is back-scattered by molecules (Rayleigh scattering), cloud and aerosol particles (Lorenz-Mie scattering). A notable amount of solar radiation reaching the Earth’s surface is reflected, either by the soil-vegetation and water systems on land or by the oceans. These processes contribute to a planetary albedo of about 30 units (102.1 W∙m−2), on global average (Figure 24). This means that the remaining 70 units (238.2 W∙m−2) of solar radiation, on global average, feed the EAS with energy. However, only about 46 units (156.5 W∙m−2) are absorbed by water (including ice) and land masses (including vegetation) in the close vicinity of the Earth’s surface (Figure 24). Note that a fraction of the incident solar radiation may penetrate into the water or snow to considerable depths without notable absorption.
Atmospheric motions are stochastic to a certain extent. However, organized patterns like Rossby waves, mountain-induced gravity and inertial-gravity waves, cyclones and anticyclones, jet streams, and circulation pattern of different sizes like Hadley, Ferrel, and polar cells, monsoonal circulation, small-scale land-sea breezes, convection roles, etc. can be identified. Variations of these flow patterns may affect the energy conversion at the Earth’s surface on local and regional scales. Turbulent motion can mainly be observed in the atmospheric boundary layer (ABL) and along jet streams. The exchange of sensible and latent heat between the land or water masses adjacent to the Earth’s surface and the atmosphere is strongly controlled by molecular and turbulent transfer processes within the Prandtl layer (also called the atmospheric surface layer, ASL). Over canopies of tall vegetation, the outer edge of the Prandtl layer may be at 100 m height above ground. The climates of landscapes may be described by the Köppen-Geiger climate classification. An updated version was derived by Peel et al.  in 2007.
During cloud formation huge amounts of latent heat are released that heat the ambient air directly. Clouds of various horizontal and vertical extensions and compositions of hydrometeors mainly occur in the troposphere. They strongly interact with both solar and infrared radiation by absorption, scattering, and emission (only infrared radiation). Clouds also affect the energy conversion at the interface Earth- atmosphere via radiation and by precipitation thereby altering the surface properties of vegetation and soils. Precipitation contributes in small amounts to the conversion of potential energy into kinetic energy which is finally converted into heat.
Solar radiation absorbed by water and soil layers adjacent to the Earth’s surface is converted into heat. Hence, it contributes to warming these layers . These layers also exchange energy with the ABL by flux densities (simply denoted as fluxes) of sensible (5 units or 17 W∙m−2) and latent heat (23 units or 78.3 W∙m−2). These fluxes, on global average, heat the atmosphere from below and cause convective transports of energy and mass into the upper troposphere. There, especially the release of latent heat during phase transition processes contributes to establish atmospheric circulation systems of different spatial and temporal scales [75,93].
As the absorption of solar radiation by atmospheric constituents and the exchange of energy between the soil and/or water layers at the Earth-atmosphere interface by the fluxes of sensible and latent heat already heated the atmosphere (about 177.0 W∙m−2 of the energetically relevant solar radiation, on global average), we have to expect that gaseous atmospheric constituents able to emit and absorb infrared (IR) radiation in finite spectral ranges, will emit energy in the IR range in all directions. The amount of this IR radiation depends on the local temperature of the mixture of these constituents. Therefore, it is indispensable to consider the down-welling IR radiation reaching the Earth’s surface, where most of it is absorbed. The same is true for hydrometeors.
The water and soil layers adjacent to the Earth-atmosphere interface, of course, also emit IR radiation depending on their local temperatures. The net emission in the IR range (emitted radiation minus absorbed down-welling radiation) is about 18 units (61.3 W∙m−2), on global average. A notable portion of this IR net emission is absorbed by atmospheric constituents and emitted in all directions, too. A small fraction propagates through the atmosphere (about 6 units or 20.4 W∙m−2) with marginal extinction by intervening constituents. Such a spectral region is the so-called atmospheric window ranging from 8.3 μm to 12.5 μm (e.g., [74,75,82,97,98]). It only contains the 9.6 μm-band of ozone. Satellite-borne radiometers use the atmospheric window between 10 μm and 12.5 μm to measure radiation up-welling from the Earth’s surface .
Gases like H2O, carbon dioxide (CO2), and O3 and hydrometeors also emit IR radiation to space. As shown in Figure 24, on global average, the net emission by these gases is about 38 units (129.3 W∙m−2), and the emission by clouds to space is 26 units (88.5 W∙m−2).
Besides the solar climate, the oceanic circulation affects the meridional distribution of normal temperatures . As Wunsch pointed out, the mass fluxes in the upper hundred meters of the ocean are primarily wind-driven yielding major features like the Gulf Stream and the Antarctic Circumpolar Current; and secondarily driven by tidal forces.
The role of the thermohaline circulation is assessed somewhat controversy. Wunsch  argued “that the “thermohaline circulation” should be reserved for the separate circulations of heat and salt, and not conflated into one vague circulation with unknown or impossible energetics. No shortcut exists for determining property fluxes from the mass circulation without knowledge of the corresponding property distribution.”On the contrary, Rahmstorf  stated that “the ocean’s density distribution, which determines pressure gradients and thus circulation, is itself affected by currents and mixing. Thermohaline and wind-driven currents therefore interact in non-linear ways and cannot be separated by oceanographic measurements.”
4.1. Global and Hemispheric Temperature Averages
To calculate the global and the hemispheric averages of the near-surface air temperature, we used Equation (8) and the datasets of the historical climatological and annual averages of temperature along numerous parallels of latitude(i.e., the zonal averages of temperature or normal temperatures) published between 1852 (Dove ) and 1913 (Börnstein ). The results are listed in Table 4. For comparison, also the results derived from the meridional distributions of the normal temperatures published by Defant and Obst , von Hann-Süring  (eventually adopted by Haurwitz and Austin  and Blüthgen ) and Sellers  as well as for the solar climate predicted by Kramm et al.  (but not explicitly published) are listed in Table 4. In some cases, we used Equation (23) to calculate the respective spheroidal mean near-surface temperature to assess the accuracy of spherical averaging as compared to spheroidal averaging.
Figure 25 shows results of the integrand ,as requested by Equation (8) for a spherical shape of the Earth, where Köppen’s  maximum and minimum temperatures for numerous parallels of latitude (Table 3) were considered. Using the polygons related to the points of this integrand yields a globally averaged minimum temperature of and a globally averaged maximum temperature of ,respectively. For fitting this small number of points of these polygons, we used the asymmetric double sigmoidal peak function,
because the temperature distributions illustrated by Figure 20 require an asymmetric peak function. The parameters A, B, C, and D are listed in Table 5. For Köppen’s  dataset, this fit function provided and ,respectively.
Table 4. Global and hemispheric averages of the near-surface air temperature based on polygons and Equation (36). Marginal differences in the second decimal place may cause some rounding effects.
*)Re-corrected to mean terrain height above sea level.
Table 5. Parameters provided by fitting the points using Equation (36) for the various datasets used in our study.
*) Re-corrected to mean terrain height above sea level.
Figure 25. The integrand as requested by Equation (8), where Köppen’s  maximum and minimum temperatures for numerous parallels of latitude listed in Table 3 were considered: (a) polygons and (b) fits provided by Equation (36).
With respect to the polygons, Köppen’s  distributions of the zonal averages of minimum and maximum temperatures provided the highest and the lowest global and hemispheric averages, respectively. Beside the latter, only the datasets Do1853  and Fo1859  provided global averages remarkably lower than ,namely and ,respectively. The results derived from the other historical distributions of normal temperatures published between 1877 and 1913 are ranging from (Ba1894 ) to (Fe1877 ), with (vB1906  and Bö1913 ) falling in between.
Obviously, the poor coverage of the Southern Hemisphere by observations during that time indicated by the normal temperatures of Dove  and Forbes  (no zonal temperature averages beyond 40˚S) led to a hemispheric average of . This value is remarkably lower than those derived from the meridional distributions of the historical normal temperatures published by the other authors. The respective averages for the Southern Hemisphere are ranging from (Ho1906-1 ) to
(Fe1877 ). The datasets vB1906 and Bö1913 provided .
On the contrary, the average for the Northern Hemisphere based on Dove  and Forbes  is . Thus, it substantially agrees with the results derived from the datasets of the other authors. Values are ranging from (Ba1894 ) to (vB1906 , Bö1913 ). Our results also confirmed that the average temperature for the Northern Hemisphere slightly exceeds that for the Southern Hemisphere. This is also true in case of the solar climate, but not in the case of Ferrel’s  dataset Fe1877 (Table 4).
For comparison: The datasets of Defant and Obst  and von Hann-Süring  yielded: , ,and for De1923 and , ,and
for vHS1939, respectively. Sellers’  distribution of zonal averages of temperature, which were not reduced to sea level, provided , , . Note that Seller’s distribution prohibits to use Equation (8) in an exact manner to calculate hemispheric averages of temperature because it omitted a zonal average of temperature for the Equator. Thus, we determined the zonal average for the Equator by interpolation. The same was also done for Arrhenius’ datasets A1896-1 and A1896-2. Obviously, the results based on Sellers’ dataset are remarkably lower than , ,and derived from the dataset A1896-2 that represents the normal temperatures of various parallels of latitude, re-corrected by Arrhenius to mean terrain height above sea level. For completeness, the differences between the hemispheric temperature averages expressed by are listed in Table 4 as well.
4.2. Data Fitting
As illustrated in Figure 25, the accuracy of the asymmetric double sigmoidal peak function (Equation (36)) seems to be acceptable, but the global averages of the minimum and maximum temperatures derived from Köppen’s  data (Table 3) notably differ from those provided by the polygons. However, as expected from a mathematical perspective, for higher numbers of points, the difference between polygons and Equation (36) should decrease remarkably, as illustrated in Figure 26 for Spitaler’s  complete dataset Sp1885 (Table 2) and an artificially reduced one based on it. Using Equation (8), the polygon and the Equation (36), provided for the complete dataset (Table 4). Whereas for the reduced dataset, the integration of the polygon yielded and that of Equation (36) provided . As aforementioned, the dataset of von Hann-Süring  (vHS1939) was adopted by Haurwitz and Austin  as well as Blüthgen . However, Haurwitz and Austin  only considered the values at each 10th degree of latitude, while Blüthgen  adopted the complete dataset. The dataset of Haurwitz and Austin  (Figure 27(a)) provided for the polygon and for Equation (36). The higher number of points as shown in Figure 27(b) yields for the polygon and based on Equation (36). The results derived from the reduced datasets are also listed in Table 4, but not further discussed.
Figure 26. Comparison of the integrand (as requested by Equation (8)), derived from Sp1885  using (a) the reduced dataset and (b) the complete one.
Figure 27. As in Figure 26, but for the datasets of (a) Haurwitz and Austin  and (b) von Hann- Süring .
By using Equation (36), we fitted the meridional distributions of the points of the integrand related to the historical distributions of normal temperatures. The parameters of our fitting procedure are listed in Table 5; the calculated global and hemispheric averages are listed in Table 4. Global averages are ranging from (vB1906 , Bö1913 ) to (Fe1877 ). The hemispheric averages are ranging from (Ho1906-1 , vB1906 , Bö1913 ) to (Fe1877 ) for the Southern Hemisphere (i.e., the results reveal that the poor coverage of the Southern Hemisphere as indicated by the normal temperatures of Dove  and Forbes  plays only a minor role) and (vB1901-1 ) to (Ho1906-2 ) for the Northern Hemisphere.
For comparison: Based on Equation (36), the datasets of Defant and Obst  and von Hann-Süring  yielded: , ,and and , ,and ,respectively. Sellers’  distribution of zonal averages of temperature (which are not reduced to sea level) provided , ,and . These results are remarkably lower than , , derived from A1896-2 that repre- sents the normal temperatures of various parallels of latitude, re-corrected by Arrhenius to mean terrain height above sea level.
4.3. Uncertainty Analysis
To estimate the uncertainty of our results, the numerical solution was repeated 50 times to create an ensemble of 50 realizations, where for each solution the zonal averages of the temperature for all parallels of latitude were randomly modified by adding temperature values, ,that are normally distributed within a standard deviation of . No seed was presupposed. This procedure was applied to each of the datasets of the zonal averages of temperature. Figure 28 illustrates examples of these randomly added temperatures. Based on the 16 ensembles of 50 realizations each, we calculated the uncertainty of the global and hemispheric averages and the differences between the hemispheres. The results are listed in Table 6. As expected, the uncertainty depends on the number of zonal averages. The smallest values of uncertainty were obtained for vB1906 (Figure 3) and Bö1913 (Figure 5) followed by vB1901-1 (Figure 2) and vHS1939 (Table 2). The largest uncertainty occurred for Do1852, Fo1859, Fe1877, A1896-1, and A1896-2.
By using Equation (8), the integration of the polygon provided, for instance, for vB1906 and Bö1913 , , ,and (Table 4). Our uncertainty ana-
lysis yielded for these datasets , , ,and
(Table 6). For comparison: vHS1939 yielded ,
, ,and (Table 4). Whereas the respective ensemble
Figure 28. Five instances of randomly added to the climatological mean temperatures of various parallels of latitude listed in Börnstein’s  table (Figure 5).
Table 6. Uncertainty estimates as obtained from the ensemble averaging procedure. Marginal differences in the second decimal place may cause some rounding effects.
*) Re-corrected to mean terrain height above sea level.
provided , , ,and (Table 6).
Presupposing an oblate spheroid, the numerical solution of Equation (23) provided for vB1906 and Bö1913 , , ,and . One of these values is slightly higher than that provided by Equation (8); however, the increase is mainly a rounding effect due to a marginal change in the second decimal place. Based on our uncertainty analysis, we obtained for these datasets , , ,and . Again, any increase is mainly a rounding effect due to a marginal change in the second decimal place. For comparison: vHS1939 provided , , ,and .
The normal temperatures of Kramm et al.  were also used to assess the accuracy of the spherical averaging procedure (Equation (8)) in comparison with spheroidal averaging procedure (Equation (23). The results confirmed that Equation (8) provides results with an acceptable accuracy (±0.04˚C) even in the case of an oblate spheroid.
5. Summary and ConclusionS
Zonal averages of temperature, the so-called normal temperatures, for numerous parallels of latitude published between 1852 and 1913 by Dove , Forbes , Ferrel , Spitaler , Batchelder , Arrhenius [14,26], Hopfner , von Hann , von Bezold [7,29], and Börnstein  served to quantify the global mean near-surface temperature of the terrestrial atmosphere. The investigation showed that only the datasets Do1852 and FO1859 of Dove  and Forbes  provided global averages below ,namely and ,respectively. The results derived from the other historical distributions of normal temperatures published between 1877 and 1913 ranged from (Ba1894 ) to (Fe1877 ), with (vB1906  and Bö1913 ) falling in between.
The poor coverage of the Southern Hemisphere by observations during that time indicated by the normal temperatures of Dove  and Forbes  (no zonal averages of temperature for parallels of latitude beyond 40°S) led to a hemispheric average of . This value is notably below those derived from the other historical datasets investigated. The respective averages calculated for the Southern Hemisphere are ranging from (Ho1906-1 ) to (Fe1877 ). The datasets vB1906 and Bö1913 of von Bezold  and Börstein  provided .
On the contrary, the average for the Northern Hemisphere of Dove  and Forbes  is . It substantially agrees with those calculated from the historical datasets of the other authors. The respective results are ranging from (Batchelder ) to (von Bezold /Börnstein ). Our results confirmed von Hann’s   conclusion that both hemispheres have nearly the same average temperature, but the Southern Hemisphere would probably be slightly cooler than the northern one. This is also true for the solar climate, but not reflected by Ferrel’s  data. Recently, Feulner et al.  confirmed von Hann’s conclusion as well. They found that the annually averaged surface air temperature in the Northern Hemisphere is 1˚C - 2˚C higher than in the Southern Hemisphere.
To estimate the uncertainty of our results, the zonal averages of temperatures for all parallels of latitude were randomly perturbed (without presupposed seed) by adding temperature values, ,that are normally distributed with a standard deviation of . For each of these historical datasets, ensembles of the 50 realizations of perturbed distributions were created. The numerical integrations of these perturbed distributions provided uncertainties in the global averages ranging from ±0.3˚C to ±0.6˚C where the magnitude of uncertainty increases with the decreasing number of normal temperatures available. The global and hemispheric means obtained from the ensembles of perturbed distributions well agreed with those derived from the original unperturbed datasets.
To assess the difference between spherical and spheroidal averaging special attention was paid to the distributions of climatological mean temperatures for numerous parallels of latitude published by von Bezold  and Börnstein’s . Global and spheroidal averaging provided and ,respectively. The corresponding ensembles provided ,and ,respectively. For comparison: The dataset of von Hann-Süring  yielded in both cases . The corresponding ensembles provided for both . This means that spherical averaging is sufficiently accurate.
Compared with our results, the hemispheric averages for the Northern Hemisphere of Dove ,
,Ferrel , ,and Spitaler , ,and for the Southern Hemisphere of Ferrel , ,and Spitaler , and the global averages of Ferrel , ,and Spitaler , ,as reported by von Hann [15,16], are remarkably higher. The same is true for the global average suggested by von Hann [15,16] and von Bezold [7,29]. Von Hann’s  value of (also mentioned by Lockyer ) seems to be sufficiently adequate, but using his numbers we only obtained . The HadCRUT4 records provided ( and ) for 1851-1880 and ( and
) for 1881-1910. The Berkeley record provided and for these periods, respectively. The NASA GISS records yielded ( and
) for 1881-1910. Obviously, these results are notably lower than those calculated from the meridional distributions of historical zonal averages of temperature. Since the HadCrut4 record yielded
( and ), the Berkeley record ,and the NASA GISS
records ( and ) for 1991-2018, the results derived from the
historical data suggest no change in the globally averaged near-surface temperature over the past 100 years.
Our results underline that reviewing the epoch-making literature from the 19th century and the first two decades of the 20th century is indispensable in the assessment of climate change since the end of the Little Ice Age in the first half of the 19th century.
We thank the anonymous reviewers for fruitful comments. We thank Google Books and the Royal College of Physicians of Edinburgh for making the textbooks of von Hann and Börnstein and numerous reports available to us. We also thank ETH-Bibliothek Zürich for making the Atlas of Meteorology: a series of over four hundred maps prepared by John G. Bartholomew and Andrew J. Herbertson and edited by Alexander Buchan available to us as well.
 Rohde, R., Muller, R., Jacobsen, R., Perlmutter, S., Rosenfeld, A., Wurtele, J., Curry, J., Wickham, C. and Mosher, S. (2013) Berkeley Earth Temperature Averaging Process. Geoinformatics & Geostatistics: An Overview, 1, 13.
 Hansen, J., Johnson, D., Lacis, A., Lebedeff, S., Lee, P., Rind, D. and Russell, G. (1981) Climate Impact of Increasing Atmospheric Carbon-Dioxide. Science, 213, 957-966.
 Jones, P.D. (1994) Hemispheric Surface Air Temperature Variations: A Reanalysis and an Update to 1993. Journal of Climate, 7, 1794-1802.
 Jones, P.D., New, M., Parker, D.E., Martin, S. and Rigor, I.G. (1999) Surface Air Temperature and Its Changes over the Past 150 Years. Reviews of Geophysics, 37, 173-199.
 Von Bezold, W. (1901) über klimatologische Mittelwerthe für ganze Breitenkreise. In: Sitzungsberichte der Königlich Preussischen, Verlag der Akademie der Wissenschaften, Berlin, 1330-1343.
 Murray, J. (1887) On the Total Annual Rainfall on the Land of the Globe, and the Relation of Rainfall to the Annual Discharge of Rivers. Scottish Geographical Magazine, 3, 65-77.
 Arrhenius, S. (1896) On the Influence of Carbonic Acid in the Air upon the Temperature of the Ground. Philosophical Magazine and Journal of Science, 41, 237-276.
 Morice, C.P., Kennedy, J.J., Rayner, N.A. and Jones, P.D. (2012) Quantifying Uncertainties in Global and Regional Temperature Change Using an Ensemble of Observational Estimates: The HadCRUT4 Data Set. Journal of Geophysical Research—Atmospheres, 117, D08101.
 Arrhenius, S. (1896) Ueber den Einfluss des atmosphärischen Kohlensäuregehalts auf die Temperatur der Erdoberfläche. In: Bihang till Kongl. Svenska Vetenskaps-Akademiens Handlingar, K. Svenska Vetenskaps-Akademien, Stockholm, 1-102.
 Buchan, A. (1889) Report on Atmospheric Circulation Based on the Observations Made on Board H.M.S. Challenger during the Years 1873-1876, and Other Meteorological Observations. In: Voyage of H.M.S. Challenger, Physics and Chemistry, London, Edinburg, Dublin.
 von Bezold, W. (1906) über Strahlungsnormalen und Mittellinien der Temperatur. In: Pernter, J.M., Von Hann, J. and Hellmann, G., Eds., Meteorologische Zeitschrift: Hann-Band zum vierzigjährigen Redaktionsjubiläum J. Hann’s von Freunden und Kollegen gewidmet, F. Vieweg, 279-287.
 Kramm, G., Dlugi, R. and Zelger, M. (2009) Comments on the “Proof of the Atmospheric Greenhouse Effect” by Arthur P. Smith.
 Kramm, G., Dlugi, R. and Mölders, N. (2017) Using Earth’s Moon as a Testbed for Quantifying the Effect of the Terrestrial Atmosphere. Natural Science, 9, 251-288.
 PaiMazumder, D. and Mölders, N. (2009) Theoretical Assessment of Uncertainty in Regional Averages Due to Network Density and Design. Journal of Applied Meteorology and Climatology, 48, 1643-1666.
 Kahlig, P. (1993) Some Aspects of Julius Von Hann’s Contribution to Modern Climatology. In: Mcbean, G. and Hantel, M., Eds., Interactions between Global Climate Subsystems: The Legacy of Hann, American Geophysical Union, Washington DC, 1-7.
 Bartholomew, J.G., Herbertson, A.J. and Buchan, A. (1899) Atlas of Meteorology: A Series of over Four Hundred Maps Prepared by John G. Bartholomew and Andrew J. Herbertson and Edited by Alexander Buchan. Geographical Institute, Edinburgh.
 Locarnini, R.A., Mishonov, A.V., Baranova, O.K., Boyer, T.P., Zweng, M.M., Garcia, H.E., Reagan, J.R., Seidov, D., Weathers, K.W., Paver, C.R. and Smolyar, I.V. (2019) World Ocean Atlas 2018, Volume 1: Temperature. US Department of Commerce, National Oceanic and Atmospheric Administration, Silver Spring, MD52.
 Feulner, G., Rahmstorf, S., Levermann, A. and Volkwardt, S. (2013) On the Origin of the Surface Air Temperature Difference between the Hemispheres in Earth’s Present-Day Climate. Journal of Climate, 26, 7136-7150.
 Berger, A. (1978) Long-Term Variations of Daily Insolation and Quaternary Climatic Changes. Journal of the Atmospheric Sciences, 35, 2362-2367.
 Keihm, S.J. (1984) Interpretation of the Lunar Microwave Brightness Temperature Spectrum: Feasibility of Orbital Heat Flow Mapping. Icarus, 60, 568-589.
 Paige, D.A., Foote, M.C., Greenhagen, B.T., Schofield, J.T., Calcutt, S., Vasavada, A.R., Preston, D.J., Taylor, F.W., Allen, C.C., Snook, K.J., Jakosky, B.M., Murray, B.C., Soderblom, L.A., Jau, B., Loring, S., Bulharowski, J., Bowles, N.E., Thomas, I.R., Sullivan, M.T., Avis, C., De Jong, E.M., Hartford, W. and McCleese, D.J. (2009) The Lunar Reconnaissance Orbiter Diviner Lunar Radiometer Experiment. Space Science Reviews, 150, 125-160.
 Vasavada, A.R., Bandfield, J.L., Greenhagen, B.T., Hayne, P.O., Siegler, M.A., Williams, J.-P. and Paige, D.A. (2012) Lunar Equatorial Surface Temperatures and Regolith Properties from the Diviner Lunar Radiometer Experiment. Journal of Geophysical Research, 117, E00H18.
 Emilio, M., Kuhn, J.R., Bush, R.I. and Scholl, I.F. (2012) Measuring the Solar Radius from Space during the 2003 and 2006 Mercury Transits. The Astrophysical Journal, 750, 135.
 Boltzmann, L. (1884) Ableitung des Stefan’schen Gesetzes, betreffend die Abhängigkeit der Wärmestrahlung von der Temperatur aus der electromagnetischen Lichttheorie. Wiedemann’s Annalen, 22, 291-294.
 Williams, J.G., Boggs, D.H. and Folkner, W.M. (2013) DE430 Lunar Orbit, Physical Librations and Surface Coordinates. In: JPL Interoffice Memorandum (Internal Document), Jet Propulsion Laboratory, California Institute of Technology, Pasadena, 19.
 Folkner, W.M., Williams, J.G., Boggs, D.H., Park, R.S. and Kuchynka, P. (2014) The Planetary and Lunar Ephemerides DE430 and DE431. In: IPN Progress Report, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, 81.
 Kopp, G., Fehlmann, A., Finsterle, W., Harber, D., Heuerman, K. and Willson, R. (2012) Total Solar Irradiance Data Record Accuracy and Consistency Improvements. Metrologia, 49, S29-S33.
 Kopp, G., Krivova, N., Wu, C.J. and Lean, J. (2016) The Impact of the Revised Sunspot Record on Solar Irradiance Reconstructions. Solar Physics, 291, 2951-2965.
 Gerlich, G. and Tscheuschner, R.D. (2009) Falsication of the Atmospheric CO2 Greenhouse Effects within the Frame of Physics. International Journal of Modern Physics B, 23, 275-364.
 Peel, M.C., Finlayson, B.L. and McMahon, T.A. (2007) Updated World Map of the Koppen-Geiger Climate Classification. Hydrology and Earth System Sciences, 11, 1633-1644.