{"id":282836,"date":"2023-10-11T14:29:11","date_gmt":"2023-10-11T12:29:11","guid":{"rendered":"https:\/\/climatescience.press\/?p=282836"},"modified":"2023-10-11T14:29:14","modified_gmt":"2023-10-11T12:29:14","slug":"testing-a-constructal-climate-model","status":"publish","type":"post","link":"https:\/\/climatescience.press\/?p=282836","title":{"rendered":"Testing A Constructal Climate Model"},"content":{"rendered":"\n<figure class=\"wp-block-image size-full\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"723\" height=\"475\" data-attachment-id=\"282854\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282854\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/0Global_Climate_Model-1024x673-1.jpg?fit=1024%2C673&amp;ssl=1\" data-orig-size=\"1024,673\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"0Global_Climate_Model-1024&amp;#215;673-1\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/0Global_Climate_Model-1024x673-1.jpg?fit=723%2C475&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/0Global_Climate_Model-1024x673-1.jpg?resize=723%2C475&#038;ssl=1\" alt=\"\" class=\"wp-image-282854\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/0Global_Climate_Model-1024x673-1.jpg?w=1024&amp;ssl=1 1024w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/0Global_Climate_Model-1024x673-1.jpg?resize=300%2C197&amp;ssl=1 300w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/0Global_Climate_Model-1024x673-1.jpg?resize=768%2C505&amp;ssl=1 768w\" sizes=\"auto, (max-width: 723px) 100vw, 723px\" \/><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\"><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong><em>The past climate is being rewritten so fast that we literally don\u2019t know what will happen yesterday \u2026<\/em><\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">From <a href=\"https:\/\/wattsupwiththat.com\/\">Watts Up With That?<\/a><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong><em>Guest Post by Willis Eschenbach<\/em><\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>ABSTRACT<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">A simple constructal model of the operation of the climate system was created by Dr. Adrian Bejan and several others. It posits that the climate system can be modeled very accurately by considering the climate as a giant heat engine turning solar power into mechanical motion. Further, it says that following the constructal law, the heat engine constantly evolves to maximize the heat flow from the tropics to the poles. In this analysis, I examine the inner workings of the model, implement a couple of improvements, and test it against the CERES satellite dataset. Sorry, no spoilers.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>CONSTRUCTAL LAW<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The Constructal Law, formulated by Professor Adrian Bejan in 1996, is a fundamental principle in physics and engineering that describes the natural tendency of all flow systems, whether inanimate or animate, to evolve and organize in a way that maximizes the flow of matter, energy, or information. This law recognizes that patterns and structures in nature, such as river networks, tree branches, and biological organisms, emerge and evolve to enhance their efficiency in the movement of resources. The constructal law explains things like the endlessly meandering nature of rivers seen in the image above. My previous posts on the Constructal Law are&nbsp;<a href=\"https:\/\/wattsupwiththat.com\/2018\/01\/20\/a-new-index-to-williss-posts\/#023\">here<\/a>.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">In essence, the Constructal Law states that the design and development of flow systems, whether the branching of blood vessels in the human body, the structure of transportation networks, or even the layout of technology and information networks, are governed by the imperative to reduce flow resistance and facilitate the transfer of resources.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The Constructal Law, as applied to climate, says that natural climate systems, such as atmospheric and oceanic circulation patterns, evolve and organize in a way that maximizes the efficiency of heat and energy flow on Earth. This principle emphasizes that climate systems, like other flow systems, tend to develop structures and patterns that reduce flow resistance and promote the transfer of heat and energy.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">In a series of three papers, \u201c<a href=\"https:\/\/www.academia.edu\/4061850\/Thermodynamic_optimization_of_global_circulation_and_climate\">Thermodynamic optimization of global circulation and climate<\/a>\u201c, \u201d&nbsp;<a href=\"https:\/\/www.academia.edu\/4061851\/Constructal_theory_of_global_circulation_and_climate\">Constructal theory of global circulation and climate<\/a>\u201c, and \u201c<a href=\"https:\/\/www.academia.edu\/4509041\/Climate_change_in_the_framework_of_the_constructal_law\">Climate change, in the framework of the constructal law<\/a>\u201c, Adrian Bejan and his co-authors show that the climate can be modeled as a heat engine. Following the Constructal Law, this climate heat engine evolves to maximize its mechanical power output. The authors say:<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\">\u201cIn conclusion, the maximization of the mechanical power output is equivalent to the maximization of the heat current from the hot region to the cold region.\u201d<\/p>\n<\/blockquote>\n\n\n\n<p class=\"wp-block-paragraph\">I got to re-reading the final of those three papers the other day, and I realized that I could set up their model on my computer. Let me start with an overview of their model.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"676\" height=\"720\" data-attachment-id=\"282838\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282838\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-273.png?fit=676%2C720&amp;ssl=1\" data-orig-size=\"676,720\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-273\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-273.png?fit=676%2C720&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-273.png?resize=676%2C720&#038;ssl=1\" alt=\"\" class=\"wp-image-282838\" style=\"width:760px;height:auto\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-273.png?w=676&amp;ssl=1 676w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-273.png?resize=282%2C300&amp;ssl=1 282w\" sizes=\"auto, (max-width: 676px) 100vw, 676px\" \/><figcaption class=\"wp-element-caption\"><em>Figure 1. The conceptual model.<\/em><\/figcaption><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\">The top<a>&nbsp;pa<\/a>rt shows the warm (tropical) and cold (poleward) areas of the global climate heat engine. These areas are marked A<sub>H<\/sub>&nbsp;and A<sub>L.&nbsp;<\/sub>for \u201cArea High\u201d and \u201cArea Low\u201d temperatures. They each have a corresponding temperature T<sub>H<\/sub>&nbsp;(temperature high) and T<sub>L<\/sub>&nbsp;(temperature low).<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The lower part of the diagram shows the various heat currents. The far left downward pointing arrow is heat from the sun to the hot zone. The next arrow, pointing up, is heat radiated from the hot zone to space.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Then we have the horizontal arrow \u201cq\u201d, the heat current from the hot zone to the cold zone.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Finally, in the cold zone on the right, we have a downward-pointing solar arrow showing heat from the sun to the cold zone, and an upward-pointing radiation arrow showing heat radiated to space.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">In short, the hot zone gets heat from the sun. Some is radiated back to space. The rest, the flow \u201cq\u201d, is transported to the cold zone. There, the flow \u201cq\u201d gets radiated back to space along with the heat that the cold zone gets from the sun.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">And most important, the Constructal Law says that the system will constantly reorganize itself to maximize the heat flow \u201cq\u201d.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Next, here\u2019s the math of the model, from the third of the papers linked above. Recall from Figure 1 above that \u201cx\u201d is the area fraction, the fraction of the globe occupied by the hot zone.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"720\" height=\"656\" data-attachment-id=\"282839\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282839\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-274.png?fit=720%2C656&amp;ssl=1\" data-orig-size=\"720,656\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-274\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-274.png?fit=720%2C656&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-274.png?resize=720%2C656&#038;ssl=1\" alt=\"\" class=\"wp-image-282839\" style=\"width:760px;height:auto\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-274.png?w=720&amp;ssl=1 720w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-274.png?resize=300%2C273&amp;ssl=1 300w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\">Daunting \u2026 so let me translate for those who like math. For those who don\u2019t, no worries\u2014just skip down to where it says \u201c<strong>THEIR MODEL RESULTS<\/strong>\u201c.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">And for the three folks still reading this section, ignore equation (26) for now. Next, in the above set of equations, rho&nbsp;<a>(<\/a>\u03c1) is the albedo, and gamma (\u03b3) is the \u201cgreenhouse factor\u201d, the fraction of upwelling surface longwave radiation that is absorbed by the atmosphere. And at steady-state, the left-hand side of equations 24 and 25 is zero\u2014there is no change of temperature with time.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">With those as prologue, the first equation (23) describes the hot zone. It says that the hot zone gets heat from the sun. Some is radiated back to space. The rest, the flow \u201cq\u201d, is transported to the cold zone. So \u201cq\u201d is equal to hot zone solar heat input minus hot zone radiation to space. In short, it\u2019s just a mathematical description of the bottom left part of Figure 1 above. Simple<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The second equation (24) describes the cold zone. It says the cold zone gets heat from the sun, plus the flow \u201cq\u201d from the hot zone, and radiates it all to space. So \u201cq\u201d is equal to the cold zone output to space minus the cold zone solar input. This equation is a mathematical description of the bottom right-hand part of Figure 1 above.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The third equation (25) says that the flow \u201cq\u201d is equal to some constant \u201cC\u201d times the 3\/2 power of the difference in temperature between the hot and cold zones.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The final equation (27) specifies that \u201cq\u201d is maximized.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">There are four unknowns in the equations\u2014temperatures of the hot and cold zones \u201cT<sub>H<\/sub>\u201d and \u201cT<sub>L<\/sub>\u201c, the heat flow \u201cq\u201d, and the area fraction \u201cx\u201d. Now, my math-fu is not strong enough to solve those four equations to determine the four unknowns. And unfortunately, the authors of the paper didn\u2019t include the solution. Grrrr.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">However, I\u2019m a determined fellow. After some reflection, I realized that I could use a double optimization process to get the answers.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">I wanted to determine the value of x (the size of the hot zone) which gives the largest value for \u201cq\u201d, the heat flow from the hot zone to the cold zone. But I only had three equations with four unknowns.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">So I divided the problem up by assuming that I knew what \u201cx\u201d was. Using that, I could then use an optimization program to give me the values of T<sub>H<\/sub>, T<sub>L<\/sub>, and q for any given value of x.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">And with that, I could use a second optimization program to give me the value of x that maximized q, the heat flow from the hot zone to the cold zone. See the Appendix below for the R code.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>THEIR MODEL RESULTS<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Here is their report of the first of their calculations. Using their same numbers, I get the same results that they show below.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"720\" height=\"323\" data-attachment-id=\"282841\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282841\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-275.png?fit=720%2C323&amp;ssl=1\" data-orig-size=\"720,323\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-275\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-275.png?fit=720%2C323&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-275.png?resize=720%2C323&#038;ssl=1\" alt=\"\" class=\"wp-image-282841\" style=\"width:760px;height:auto\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-275.png?w=720&amp;ssl=1 720w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-275.png?resize=300%2C135&amp;ssl=1 300w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\">Using their values, I was able to reproduce their results very accurately.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>PROBLEMS WITH THEIR MODEL<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">However, there are a couple of issues with their values. First, as they note, their value for \u201cx\u201d puts the limits of the hot zone at about 57\u00b0N\/S. But that\u2019s not the case in the real world. Here\u2019s the real-world data regarding the heat flow \u201cq\u201d.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"720\" height=\"597\" data-attachment-id=\"282842\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282842\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-276.png?fit=720%2C597&amp;ssl=1\" data-orig-size=\"720,597\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-276\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-276.png?fit=720%2C597&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-276.png?resize=720%2C597&#038;ssl=1\" alt=\"\" class=\"wp-image-282842\" style=\"width:760px;height:auto\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-276.png?w=720&amp;ssl=1 720w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-276.png?resize=300%2C249&amp;ssl=1 300w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><figcaption class=\"wp-element-caption\"><em>Figure 2. How much heat is moved from the tropics to the poles (positive values), and how much heat is absorbed in the polar regions (negative values). The hot zone is the red to yellow part bordered by the black\/white lines. The cold zone is shown in green to blue, outside the black\/white lines.<\/em><\/figcaption><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\">You can see the similarity of this graphic with the model shown in Figure 1 above. However, in the real world, the hot zone fraction \u201cx\u201d is about 0.55 of the total surface. This corresponds with a hot zone extending to about 34\u00b0N\/S. So that was the first problem\u2014the hot zone extends to about 34\u00b0N\/S, not 57\u00b0N\/S.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The second problem is that their equation gives far too cold a result for the cold zone. They say it averages 258.4K, which is -14.75\u00b0C. But in the real world, the cold zone poleward of 57\u00b0N\/S actually has an average temperature of about \u2013 3\u00b0C, far from the minus 14\u00b0C they claim.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>IMPROVING THEIR MODEL<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">So of course, being the eternal tinkerer, I had to see whether I could improve their model. The first thing I noticed was that they are using the same albedo and the same greenhouse factor for both the cold and hot zones. But in the real world, both the albedo and the greenhouse factor are very different for the two areas. As a result, their model is giving inaccurate results<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Using individual albedo and greenhouse factors for the two areas made the model far more accurate. But there was still a problem. The hot temperatures it calculated were too hot and the cold temperatures were too cold to match the real world. Looking at the equations, I realized that this inter-temperature distance is controlled by the constant \u201cC\u201d in Equation (25). This is the \u201cconductance\u201d, a measure of how much heat flow is generated by a given temperature difference between the hot and cold zones. The value they were using for \u201cC\u201d was far too small, which meant it required a much greater temperature difference to get the same flow, resulting in a hot zone that\u2019s too hot and a cold zone that\u2019s too cold.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Once the factor \u201cC\u201d was increased, the results looked very good.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>GROUND-TRUTHING THE MODEL<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">With that model up and running on my computer, I figured that I could test whether in fact, the climate system actually does operate as a gigantic heat engine that is continually evolving to maximize the tropical-polar heat flow. Here was my plan.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The constructal model says that given the albedo and greenhouse factors, for each value of \u201cx\u201d (the area of the hot zone) there will be a preferred temperature for the hot and cold zones. Further, the model says that the average final temperatures will be the ones that maximize \u201cq\u201d, the heat flow from the hot zone to the cold zone. I realized we could test those claims using the CERES data.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">For each year, the average top of atmosphere net radiation CERES data gives us the observed value of \u201cx\u201d in the constructal model. As mentioned above, x is the fraction of the globe that is exporting heat on average. The CERES data also gives us the information needed to calculate rho (<a>\u03c1<\/a>), which is the albedo, and gamma (<a>\u03b3<\/a>) which is the \u201cgreenhouse factor\u201d.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The model says that if we know the albedo \u03c1, the greenhouse factor \u03b3, and the hot zone area x, given those physical constraints the resulting hot and cold temperatures will be the ones that maximize the heat flow \u201cq\u201d from the hot zone to the cold zone.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Here is the performance of the constructal model. Recall that it has only one tuned parameter, C, that regulates how easily the heat flows from the hot zone to the cold zone. I\u2019ll get back in a bit to why I think their value for C (.181) is far too low. In the meantime, these are the actual (blue\/cyan) and modeled (red\/orange) temperatures for the hot and cold zones of the planet.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"720\" height=\"666\" data-attachment-id=\"282844\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282844\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-277.png?fit=720%2C666&amp;ssl=1\" data-orig-size=\"720,666\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-277\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-277.png?fit=720%2C666&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-277.png?resize=720%2C666&#038;ssl=1\" alt=\"\" class=\"wp-image-282844\" style=\"width:760px;height:auto\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-277.png?w=720&amp;ssl=1 720w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-277.png?resize=300%2C278&amp;ssl=1 300w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><figcaption class=\"wp-element-caption\"><em>Figure 3. Modeled and actual temperatures of the world\u2019s hot and cold zones<\/em><\/figcaption><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\">I found this result to be most encouraging. Those model temperatures are calculated based solely on maximizing the heat \u201cq\u201d flowing from the hot zone to the cold zone, subject to the physical constraints of the albedo and the greenhouse factor. And although the conductance C is tuned, all that tunes is the temperature difference between the hot and cold zones. It does not tune the temperatures themselves. There was no guarantee that tuning the conductance would match the absolute temperatures of the hot and cold zones \u2026 but in the event, the match is excellent. I would say that that is very convincing evidence that the constructal model accurately portrays how the climate flow system actually works.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>A SECOND TEST<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">But wait, as they say on TV, there\u2019s more. Here are closeups of the actual and modeled variations in the yearly average temperatures of the hot and cold zones.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"720\" height=\"666\" data-attachment-id=\"282845\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282845\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-278.png?fit=720%2C666&amp;ssl=1\" data-orig-size=\"720,666\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-278\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-278.png?fit=720%2C666&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-278.png?resize=720%2C666&#038;ssl=1\" alt=\"\" class=\"wp-image-282845\" style=\"width:760px;height:auto\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-278.png?w=720&amp;ssl=1 720w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-278.png?resize=300%2C278&amp;ssl=1 300w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><\/figure>\n\n\n\n<figure class=\"wp-block-image size-full\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"720\" height=\"661\" data-attachment-id=\"282847\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282847\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-279.png?fit=720%2C661&amp;ssl=1\" data-orig-size=\"720,661\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-279\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-279.png?fit=720%2C661&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-279.png?resize=720%2C661&#038;ssl=1\" alt=\"\" class=\"wp-image-282847\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-279.png?w=720&amp;ssl=1 720w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-279.png?resize=300%2C275&amp;ssl=1 300w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><figcaption class=\"wp-element-caption\"><em>Figure 4. Modeled and actual annual average temperatures of the world\u2019s hot and cold zones.<\/em><\/figcaption><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\">Not perfect, but not bad either. So not only does the constructal model give good long-term average temperatures. It also does a decent job of replicating the year-by-year variations in temperature.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">And it\u2019s doing all that using nothing more than the hot zone area \u201cx\u201d, the albedo \u201crho\u201d, and the greenhouse factor \u201cgamma\u201d to calculate the temperatures that maximize \u201cq\u201d.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">That\u2019s very clear evidence that in the real world, various physical processes constantly evolve and act to increase the flow of heat from the tropics to the polar regions.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>A FINAL TEST<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Further evidence that the model is an accurate representation of how the climate heat engine really works is visible in both the size and the stability of the area of the hot zone. The model calculates the average of x, the hot zone fraction of the surface, as being 0.564. The actual CERES 22-year average value for x is 0.556. That\u2019s less than a hundredth difference. Once again, the model is accurate.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Regarding stability, remember that x, the hot zone area fraction, is calculated by the model as the hot zone area that maximizes the heat flow \u201cq\u201d. Bear in mind that the hot zone fraction could vary from ~0.1 to ~0.9. And there\u2019s no reason to assume ex-ante that it would remain stable over time.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">However, under the constructal model, since the underlying constraints (annual average albedo and greenhouse fraction) are relatively stable we\u2019d expect the hot zone fraction \u201cx\u201d to be pretty stable as well. In any case, here\u2019s the actual record of the CERES data for \u201cx\u201d, the hot area fraction, along with the constructal model output of the same variable.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"720\" height=\"661\" data-attachment-id=\"282848\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282848\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-280.png?fit=720%2C661&amp;ssl=1\" data-orig-size=\"720,661\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-280\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-280.png?fit=720%2C661&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-280.png?resize=720%2C661&#038;ssl=1\" alt=\"\" class=\"wp-image-282848\" style=\"width:760px;height:auto\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-280.png?w=720&amp;ssl=1 720w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-280.png?resize=300%2C275&amp;ssl=1 300w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><figcaption class=\"wp-element-caption\"><em>Figure 5. The \u201cx\u201d fraction, the amount of the earth\u2019s surface that makes up the hot zone.<\/em><\/figcaption><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\">Clearly, the model is doing an excellent job of representing the real world.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">In Figure 5, as in Figs. 3 and 4 above, it\u2019s important to remember that the output (e.g. the modeled x fraction in Fig. 5 above) is not calculated directly from the input. In Figure 5, for example, the x fraction shown in red is not directly calculated from the albedo and greenhouse fraction figures.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Instead, it is the result of a maximization procedure. The x fraction shown in red in Figure 5 is the value of x that, given the physical constraints of albedo and greenhouse fraction, gives the greatest flow \u201cq\u201d from the hot zone to the cold zone.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>TEMPERATURES<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">For temperatures, I\u2019ve used the CERES surface upwelling longwave data converted using the Stefan-Boltzmann constant and monthly gridded emissivity values. I\u2019ve checked the results and they are extremely similar to both the Berkeley Earth and the HadCRUT datasets. I use it because it is energy-balanced with the rest of the CERES energy flows.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>CONDUCTANCE<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">I mentioned above that I\u2019d explain why I think their value for \u201cC\u201d, the \u201cconductance\u201d, is too low. This conductance is a measure of how much heat flows between the two zones for some given temperature difference between the zones. In their model, they\u2019ve modeled the heat transport via the atmosphere. And they\u2019ve modeled the atmospheric heat transport as being driven by the buoyancy of the warmer, lighter tropical air.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">And that is good as far as it goes. But it leaves out a couple of things. One is a main power source driving the Hadley cell circulation\u2014the perennial line of thunderstorms along the inter-tropical convergence zone (ITCZ). These drive air vertically from the surface up to the upper troposphere, and occasionally even into the stratosphere. These thunderstorms turbocharge the Hadley cell circulation, allowing it to move much more heat polewards than if it were driven solely by the general tropical-extratropical temperature differences as the authors\u2019 analysis assumes. Here\u2019s a map of where the thunderstorms live.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"720\" height=\"593\" data-attachment-id=\"282850\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282850\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-281.png?fit=720%2C593&amp;ssl=1\" data-orig-size=\"720,593\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-281\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-281.png?fit=720%2C593&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-281.png?resize=720%2C593&#038;ssl=1\" alt=\"\" class=\"wp-image-282850\" style=\"width:760px;height:auto\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-281.png?w=720&amp;ssl=1 720w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-281.png?resize=300%2C247&amp;ssl=1 300w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><figcaption class=\"wp-element-caption\"><em>Figure 6. The altitude of the cloud tops, day\/night. High altitude cloud tops are the sign of the tropical thunderstorms driving deep tropical convection. The Inter-Tropical Convergence Zone (ITCZ), where the two atmospheric hemispheres converge, is marked by the band of thunderstorms around the world at 5\u00b0-10\u00b0 north of the equator.<\/em><\/figcaption><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\">The second reason that I think their conductance value is too small is that a large amount of heat is physically moved polewards by the ocean currents. The Agulhas Current in the Indian Ocean and the Gulf Stream in the Atlantic Ocean are constantly transporting warm tropical waters polewards.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">In the Pacific, the El Nino\/La Nina pumping action periodically strips off the warm top layer of vast areas of the tropical Pacific Ocean and moves that warm water first eastwards and then towards both poles.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img data-recalc-dims=\"1\" loading=\"lazy\" decoding=\"async\" width=\"720\" height=\"339\" data-attachment-id=\"282852\" data-permalink=\"https:\/\/climatescience.press\/?attachment_id=282852\" data-orig-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-282.png?fit=720%2C339&amp;ssl=1\" data-orig-size=\"720,339\" data-comments-opened=\"1\" data-image-meta=\"{&quot;aperture&quot;:&quot;0&quot;,&quot;credit&quot;:&quot;&quot;,&quot;camera&quot;:&quot;&quot;,&quot;caption&quot;:&quot;&quot;,&quot;created_timestamp&quot;:&quot;0&quot;,&quot;copyright&quot;:&quot;&quot;,&quot;focal_length&quot;:&quot;0&quot;,&quot;iso&quot;:&quot;0&quot;,&quot;shutter_speed&quot;:&quot;0&quot;,&quot;title&quot;:&quot;&quot;,&quot;orientation&quot;:&quot;0&quot;}\" data-image-title=\"image-282\" data-image-description=\"\" data-image-caption=\"\" data-large-file=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-282.png?fit=720%2C339&amp;ssl=1\" src=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-282.png?resize=720%2C339&#038;ssl=1\" alt=\"\" class=\"wp-image-282852\" style=\"width:760px;height:auto\" srcset=\"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-282.png?w=720&amp;ssl=1 720w, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/image-282.png?resize=300%2C141&amp;ssl=1 300w\" sizes=\"auto, (max-width: 720px) 100vw, 720px\" \/><\/figure>\n\n\n\n<p class=\"wp-block-paragraph\">Because their model doesn\u2019t include either thunderstorms or ocean currents, their estimate of the conductance is an order of magnitude too small.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>CLIMATE SENSITIVITY<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">This constructal model points out some interesting things about climate sensitivity.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">First, sensitivity is a function of changes in rho (albedo) and gamma (greenhouse fraction). But not a direct function. It is the result of physical processes that maximize \u201cq\u201d given the constraints of rho and gamma.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Next, the sensitivity is slightly different depending on whether the changes in albedo and greenhouse fraction are occurring in the hot zone, the cold zone, or both.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Finally, assuming that there is a uniform pole-to-pole increase of 3.7 W\/m2 in downwelling radiation from changes in either albedo or greenhouse fraction, the constructal model shows a temperature increase of ~1.1\u00b0C. (3.7 W\/m2 is the amount of radiation increase predicted to occur from a doubling of CO2.)<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>CONCLUSIONS<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">The CERES data shows that the constructal model of the climate system is very consistent with real-world observations. This model views the climate system as a heat engine that, following the constructal law, constantly acts and evolves to maximize the flow of heat from the warm zone of the planet to the cold zone.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">This simple three-equation constructal climate model, given only information about the earth\u2019s hot zone area and the albedo and greenhouse fractions in the earth\u2019s hot and cold zones, is able to calculate the absolute temperatures of the earth\u2019s hot and cold zones to within a degree or so \u2026 a result that I found quite surprising.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Anyhow, that\u2019s what I did with my weekend. And meanwhile, back in the real world, the past climate is being rewritten so fast that we literally don\u2019t know what will happen yesterday \u2026<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Best to all,<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">w.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>APPENDIX<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">Here is the R code for the optimization programs. Read the linked paper for the full description of their method.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\">First, the inner optimization program that calculates T<sub>H<\/sub>, T<sub>L<\/sub>, and q when given x.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>maxq=function(par2){<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; theansmax=function(par){<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; th=par[1]<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tl=par[2]<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; q=par[3]<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (v1=x*((asin(x)+x*sqrt(1-x^2))\/(2*pi*x))*(1-rhoh)-<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; x*(1-gammah)*th^4-q)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (v2=(1-x)*((pi\/2-asin(x)-x*sqrt(1-x^2))\/(2*pi*(1-x)))*(1-rhoc)-<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (1-x)*(1-gammac)*tl^4+q)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (v3=1.8*(th-tl)^(3\/2)-q)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; sum(v1^2+v2^2+v3^2)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; par=c(.7,.6,.1)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; x=par2<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (par=optim(par,theansmax)$par)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\">&nbsp;&nbsp;<strong>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; par[3]<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>}<\/strong><\/p>\n<\/blockquote>\n\n\n\n<p class=\"wp-block-paragraph\">Next, the outer optimization program that calls the inner program.<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>(bestx=optim(par2,maxq,<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>control=list(fnscale = -1,reltol=1e-10),<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>method=\u201dBrent\u201d,<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>lower=.001,upper=.999)$par)<\/strong><\/p>\n<\/blockquote>\n\n\n\n<p class=\"wp-block-paragraph\">Next, some support functions:<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>surfaream = 5.100656e+14<\/strong>&nbsp;#earth surface area in sq. m.<\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>qtoq= function(q) q*((5.67e-8)*392.8^4*surfaream)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>tunscale=function(tscale) tscale*392.8<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>xtolat=function(x) degrees(asin(x))<\/strong><\/p>\n<\/blockquote>\n\n\n\n<p class=\"wp-block-paragraph\">And to get the final output:<\/p>\n\n\n\n<blockquote class=\"wp-block-quote is-layout-flow wp-block-quote-is-layout-flow\">\n<p class=\"wp-block-paragraph\"><strong>(x=bestx)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>(nupar=optim(par,theansmax)$par)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>q=nupar[3]<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>qtoq(nupar[3])<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>xtolat(x)<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>(th=ktoc(tunscale(nupar[1])))<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><strong>(tl=ktoc(tunscale(nupar[2])))<\/strong><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><\/p>\n\n\n\n<p class=\"wp-block-paragraph\"><\/p>\n<\/blockquote>\n","protected":false},"excerpt":{"rendered":"<p>The past climate is being rewritten so fast that we literally don\u2019t know what will happen yesterday \u2026 From Watts Up With That? Guest Post by Willis Eschenbach ABSTRACT A simple constructal model of the operation of the climate system was created by Dr. Adrian Bejan and several others. It posits that the climate system [&hellip;]<\/p>\n","protected":false},"author":121246920,"featured_media":282855,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_coblocks_attr":"","_coblocks_dimensions":"","_coblocks_responsive_height":"","_coblocks_accordion_ie_support":"","_crdt_document":"","advanced_seo_description":"","jetpack_seo_html_title":"","jetpack_seo_noindex":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_publicize_message":"The past climate is being rewritten so fast that we literally don\u2019t know what will happen yesterday \u2026","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":false,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2},"jetpack_post_was_ever_published":false},"categories":[1],"tags":[691823517,691818056,691821072,691823516,691818267],"class_list":{"0":"post-282836","1":"post","2":"type-post","3":"status-publish","4":"format-standard","5":"has-post-thumbnail","6":"hentry","7":"category-uncategorized","8":"tag-ceres-data","9":"tag-climate-change","10":"tag-climate-model","11":"tag-greenhouse-factor","12":"tag-solar-power","14":"fallback-thumbnail"},"jetpack_publicize_connections":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2023\/10\/0Global_Climate_Model-1024x673-1-1.jpg?fit=1024%2C673&ssl=1","jetpack_likes_enabled":true,"jetpack_sharing_enabled":true,"jetpack_shortlink":"https:\/\/wp.me\/paxLW1-1bzS","jetpack-related-posts":[{"id":337086,"url":"https:\/\/climatescience.press\/?p=337086","url_meta":{"origin":282836,"position":0},"title":"Intro to Climate\u00a0Fallacies","author":"uwe.roland.gross","date":"20\/07\/2024","format":false,"excerpt":"In GCMs, the equilibrium climate sensitivity is an \u2018emergent property\u2019 that is not directly calibrated or tuned.","rel":"","context":"In \"Climate change\"","block_context":{"text":"Climate change","link":"https:\/\/climatescience.press\/?tag=climate-change"},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/07\/2quote-star-trek-scotty.jpg?fit=1200%2C881&ssl=1&resize=350%2C200","width":350,"height":200,"srcset":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/07\/2quote-star-trek-scotty.jpg?fit=1200%2C881&ssl=1&resize=350%2C200 1x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/07\/2quote-star-trek-scotty.jpg?fit=1200%2C881&ssl=1&resize=525%2C300 1.5x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/07\/2quote-star-trek-scotty.jpg?fit=1200%2C881&ssl=1&resize=700%2C400 2x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/07\/2quote-star-trek-scotty.jpg?fit=1200%2C881&ssl=1&resize=1050%2C600 3x"},"classes":[]},{"id":386270,"url":"https:\/\/climatescience.press\/?p=386270","url_meta":{"origin":282836,"position":1},"title":"Toward A New Climate Paradigm","author":"uwe.roland.gross","date":"01\/07\/2025","format":false,"excerpt":"What got me into investigating climate science a quarter-century ago? The astounding stability of the Earth\u2019s global surface mean temperature (GMST). For example, since the year 1900, a century and a quarter, the Earth\u2019s temperature has gone up by 0.46 \u00b1 0.07 percent.","rel":"","context":"In \"Bejan\"","block_context":{"text":"Bejan","link":"https:\/\/climatescience.press\/?tag=bejan"},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/thumb-1920-716585.jpg?fit=1200%2C735&ssl=1&resize=350%2C200","width":350,"height":200,"srcset":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/thumb-1920-716585.jpg?fit=1200%2C735&ssl=1&resize=350%2C200 1x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/thumb-1920-716585.jpg?fit=1200%2C735&ssl=1&resize=525%2C300 1.5x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/thumb-1920-716585.jpg?fit=1200%2C735&ssl=1&resize=700%2C400 2x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/thumb-1920-716585.jpg?fit=1200%2C735&ssl=1&resize=1050%2C600 3x"},"classes":[]},{"id":332388,"url":"https:\/\/climatescience.press\/?p=332388","url_meta":{"origin":282836,"position":2},"title":"Avoiding the Net Zero Trap","author":"uwe.roland.gross","date":"11\/06\/2024","format":false,"excerpt":"Clintel has just put out a new publication, We Can Still Avoid the Net Zero Trap, by Kees de Lange and Guus Berkhout.","rel":"","context":"In \"Climate change\"","block_context":{"text":"Climate change","link":"https:\/\/climatescience.press\/?tag=climate-change"},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/06\/image-257.png?fit=1200%2C783&ssl=1&resize=350%2C200","width":350,"height":200,"srcset":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/06\/image-257.png?fit=1200%2C783&ssl=1&resize=350%2C200 1x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/06\/image-257.png?fit=1200%2C783&ssl=1&resize=525%2C300 1.5x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/06\/image-257.png?fit=1200%2C783&ssl=1&resize=700%2C400 2x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2024\/06\/image-257.png?fit=1200%2C783&ssl=1&resize=1050%2C600 3x"},"classes":[]},{"id":388994,"url":"https:\/\/climatescience.press\/?p=388994","url_meta":{"origin":282836,"position":3},"title":"Learning With The Dumbest Genius Librarian","author":"uwe.roland.gross","date":"16\/07\/2025","format":false,"excerpt":"Since then, I came across a fascinating way to use any AI for research into any topic. Want to learn about something fast? Here\u2019s how. Want to drill down into some recondite corner of a topic? This will do it. Whether you\u2019re looking for the broadest overview or the most\u2026","rel":"","context":"In \"AI (artificial intelligence)\"","block_context":{"text":"AI (artificial intelligence)","link":"https:\/\/climatescience.press\/?tag=ai-artificial-intelligence-2"},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/AQOfkd6UTh06r19sAIITcV0WnBpItcTDIrPR7FcFHsp0SKCOf7TfvK8BqgIVy13UyugSRzSbob6cOGOxR9Lkc6ykoJFx6AER7bSL-01QThGh3agANwnHcsf-jZbjsEhmMKQtLrLcl-PDKqj93eGtCzVIag.jpeg?fit=1200%2C1200&ssl=1&resize=350%2C200","width":350,"height":200,"srcset":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/AQOfkd6UTh06r19sAIITcV0WnBpItcTDIrPR7FcFHsp0SKCOf7TfvK8BqgIVy13UyugSRzSbob6cOGOxR9Lkc6ykoJFx6AER7bSL-01QThGh3agANwnHcsf-jZbjsEhmMKQtLrLcl-PDKqj93eGtCzVIag.jpeg?fit=1200%2C1200&ssl=1&resize=350%2C200 1x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/AQOfkd6UTh06r19sAIITcV0WnBpItcTDIrPR7FcFHsp0SKCOf7TfvK8BqgIVy13UyugSRzSbob6cOGOxR9Lkc6ykoJFx6AER7bSL-01QThGh3agANwnHcsf-jZbjsEhmMKQtLrLcl-PDKqj93eGtCzVIag.jpeg?fit=1200%2C1200&ssl=1&resize=525%2C300 1.5x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/AQOfkd6UTh06r19sAIITcV0WnBpItcTDIrPR7FcFHsp0SKCOf7TfvK8BqgIVy13UyugSRzSbob6cOGOxR9Lkc6ykoJFx6AER7bSL-01QThGh3agANwnHcsf-jZbjsEhmMKQtLrLcl-PDKqj93eGtCzVIag.jpeg?fit=1200%2C1200&ssl=1&resize=700%2C400 2x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/07\/AQOfkd6UTh06r19sAIITcV0WnBpItcTDIrPR7FcFHsp0SKCOf7TfvK8BqgIVy13UyugSRzSbob6cOGOxR9Lkc6ykoJFx6AER7bSL-01QThGh3agANwnHcsf-jZbjsEhmMKQtLrLcl-PDKqj93eGtCzVIag.jpeg?fit=1200%2C1200&ssl=1&resize=1050%2C600 3x"},"classes":[]},{"id":224301,"url":"https:\/\/climatescience.press\/?p=224301","url_meta":{"origin":282836,"position":4},"title":"Climatism is a Logic Fail","author":"uwe.roland.gross","date":"17\/10\/2022","format":false,"excerpt":"Two fallacies ensure meaningless public discussion about climate \u201ccrisis\u201d or \u201cemergency.\u201d H\/T to Terry Oldberg for comments and writings prompting me to post on this topic. One corruption is the numerous times climate claims include fallacies of Equivocation. For instance, \u201cclimate change\u201d can mean all observed events in nature, but\u2026","rel":"","context":"Similar post","block_context":{"text":"Similar post","link":""},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2022\/10\/image-818.png?fit=734%2C652&ssl=1&resize=350%2C200","width":350,"height":200,"srcset":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2022\/10\/image-818.png?fit=734%2C652&ssl=1&resize=350%2C200 1x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2022\/10\/image-818.png?fit=734%2C652&ssl=1&resize=525%2C300 1.5x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2022\/10\/image-818.png?fit=734%2C652&ssl=1&resize=700%2C400 2x"},"classes":[]},{"id":407486,"url":"https:\/\/climatescience.press\/?p=407486","url_meta":{"origin":282836,"position":5},"title":"A Submission to the Australian Senate Investigation into Climate Misinformation","author":"uwe.roland.gross","date":"10\/10\/2025","format":false,"excerpt":"When the Australian government set up a \u201cSenate select committee to investigate misinformation and disinformation in climate science and energy\u201d, I didn\u2019t find out until a few hours before the deadline for submissions from the general public. Things were made a little more complicated by the fact that I was\u2026","rel":"","context":"In \"Australian Government\"","block_context":{"text":"Australian Government","link":"https:\/\/climatescience.press\/?tag=australian-government"},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/10\/0Screenshot-2025-10-10-144810.png?fit=1200%2C788&ssl=1&resize=350%2C200","width":350,"height":200,"srcset":"https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/10\/0Screenshot-2025-10-10-144810.png?fit=1200%2C788&ssl=1&resize=350%2C200 1x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/10\/0Screenshot-2025-10-10-144810.png?fit=1200%2C788&ssl=1&resize=525%2C300 1.5x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/10\/0Screenshot-2025-10-10-144810.png?fit=1200%2C788&ssl=1&resize=700%2C400 2x, https:\/\/i0.wp.com\/climatescience.press\/wp-content\/uploads\/2025\/10\/0Screenshot-2025-10-10-144810.png?fit=1200%2C788&ssl=1&resize=1050%2C600 3x"},"classes":[]}],"_links":{"self":[{"href":"https:\/\/climatescience.press\/index.php?rest_route=\/wp\/v2\/posts\/282836","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/climatescience.press\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/climatescience.press\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/climatescience.press\/index.php?rest_route=\/wp\/v2\/users\/121246920"}],"replies":[{"embeddable":true,"href":"https:\/\/climatescience.press\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=282836"}],"version-history":[{"count":8,"href":"https:\/\/climatescience.press\/index.php?rest_route=\/wp\/v2\/posts\/282836\/revisions"}],"predecessor-version":[{"id":282856,"href":"https:\/\/climatescience.press\/index.php?rest_route=\/wp\/v2\/posts\/282836\/revisions\/282856"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/climatescience.press\/index.php?rest_route=\/wp\/v2\/media\/282855"}],"wp:attachment":[{"href":"https:\/\/climatescience.press\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=282836"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/climatescience.press\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=282836"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/climatescience.press\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=282836"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}