4 Climate System

4.1 EEI - Earth Energy Imbalance

The explanation below from Trenberth & Cheng’s 2022 paper has helped me better understand the extreme weather from the perspective of EEI (Earth’s Energy Imbalance): “About 93% of the extra heat from the EEI ends up in the ocean as increasing ocean heat content (OHC)……On average nearly 3% of the EEI goes into melting ice and another 4% goes into raising temperatures of land and melting permafrost, while less than 1% remains in the atmosphere. The changes in GMST relate to those on land and sea surface temperatures (SSTs) in the ocean, and the increase in GMST is a consequence. The primary reasons for these distributions of excess heat relate to the heat capacity of these climate system components, the specific heat of water versus land versus air, and the masses involved….However, the initial excess energy has profound effects and big impacts along the way to its destination. The radiative energy is mostly absorbed at the surface where it mainly contributes to increased evaporation and moistening of the atmosphere, and is eventually realized as latent heating of the atmosphere during precipitation. This invigorates the hydrological cycle. In regions where it is not raining, it leads to enhanced drying then heating, increasing risk of heat waves and wildfires on land. It increases intensity of droughts. The evaporated moisture, carried by the atmosphere, converges into storms and associated frontal systems, where it increases precipitation rates, increasing risk of flooding on land. The latent heat released may also intensify some weather systems such as hurricanes and convection (thunderstorms). The warmer moister air is more buoyant and, helped by many weather phenomena, it rises and expands in lower pressures and cools adiabatically, then spreads out maybe thousands of km away and warms as it subsides, thus redistributing energy often to higher and drier regions where it can radiate to space. In general, all components of the climate system react to heating by trying to get rid of excess heat in one way or another. The most effective method overall globally is radiative cooling as higher temperatures increase radiation by the fourth power of absolute temperature. Moreover, unequal heating leads to gradients that drive instabilities in the atmosphere (convective, baroclinic) and mixing is pervasive. In the ocean, warmer waters and those with lower salinity are less dense, and thermohaline circulations may develop, although most ocean currents are driven by winds. In ice and land, energy may be redistributed by water flows or conduction, but the latter is a very inefficient process”.

Trenberth (2022) A perspective on climate change from Earth’s energy imbalance (pdf)

Reis (2023) tweet on Trenbech

4.2 Convection

Zhang

Figure: For the first time mesoscale convective systems (MCSs) in both the tropics and midlatitudes and all seasons can be tracked over many years by a new algorithm jointly using satellite observed cloud-top temperature and surface precipitation features at hourly and 10-km resolution globally (top panel). Results show that MCSs account for over 50% of the annual rainfall across the tropics and many regions of the subtropics and midlatitudes (bottom panel).

Mesoscale convective systems (MCSs) are a key component in the Earth’s energy and hydrological cycles. They can grow to hundreds of kilometers in size, last for more than a day, and produce a majority of the annual rainfall in many regions of the world.

Past efforts to develop MCS databases have been limited to the tropics and used methodologies not well tested in the midlatitudes. Feng et al. [2021] developed a new methodology to track MCSs globally using high-resolution satellite observations of both cloud and precipitation. The new method significantly improves the detection of MCSs in the midlatitudes. This new storm tracking database is the first to cover both the tropics and midlatitudes for all seasons.

The study shows that MCSs account for over 50 per cent of the annual rainfall across the tropics and many regions of the subtropics and midlatitudes. Storms over land have more intense convection, while those over oceans produce heavier rainfall and last longer.

This global MCS database supports a broad range of research such as understanding the role of MCSs in global extreme rainfall and circulation, and evaluation of global weather and climate model simulations.

Zhang

Convective storms of mesoscale dimension are a key component in the Earth’s energy and hydrological cycle. Mesoscale storms grow to hundreds of kilometers in size and can last for more than a day, and produce a majority of the annual rainfall in many regions of the world. Past studies of mesoscale storms have been limited to the tropics and used methodologies not well tested in the midlatitudes. Here, we develop a new methodology to track mesoscale storms globally using high-resolution satellite observations of both cloud and precipitation. The satellite-based storm tracking reproduces important storm statistics derived from ground-based radar observations. Our new method significantly improves the detection of mesoscale storms in the midlatitudes. This new storm tracking database is the first to cover both the tropics and midlatitudes for all seasons. Results show that mesoscale convective storms account for over 50% of annual rainfall across the tropics and many regions of the subtropics and midlatitudes. Storms over land have more intense convection, while those over ocean produce heavier rainfall and last longer. This global mesoscale storms tracking database supports a broad range of applications, such as understanding their role in global extreme rainfall and circulation and evaluation of global weather and climate model simulations.

Feng (2021) Mesoscale Convective System (pdf)

4.3 El Nino

Saravan

El Niño, a phenomenon first observed centuries ago by Peruvian fishermen, refers to a sporadic warming of the tropical Pacific Ocean that occurs every 3-7 years. These waters, usually cold and teeming with nutrients, would warm considerably, decimating fish populations. The phenomenon typically occurred around December or January, hence its name, El Niño—Spanish for “the boy child”, a reference to Christ.

El Niño’s impact, propagated by atmospheric waves known as teleconnections, extends to regions far removed from the tropical Pacific. For instance, it intensifies winter rainfall in California while potentially inducing drought in Australia. Recent warming in the Pacific Ocean has prompted the US weather agency to declare an official El Niño event for 2023. When the tropical Pacific warms, it invariably raises the global-average temperature, because the Pacific covers such a large fraction of the globe. This basic arithmetic relationship has fueled a flurry of articles linking El Niño with global warming.

El Niño is a quasi-cyclical natural phenomenon whose physical cause and spatial impacts are distinctly different from human-induced global warming.

When preparing this article, despite a ton of published material on El Niño over several decades, finding graphics that directly compared El Niño and global warming in a relevant way wasn’t easy. So, I wrote code to generate all the figures in this article (except for the schematic). The code is available as a Python notebook on Google Colaboratory

A better set of clickable plots of El Niño impacts

Figure: How El Niño affects weather during December to February.

The trade winds blow from the cool east to the warm west, pushing the tropical Pacific waters westward and creating the cool upwelling region off Peru’s coast. El Niño temporarily disrupts this process and warms the eastern Pacific. Forecasting the specifics of how these east-west temperature gradients will change in the coming decades turns out to be more challenging than predicting changes in global-average temperature. In fact, climate models have yet to reproduce the cooling trend observed in the eastern tropical Pacific over the last 40 years.

Given that the effects of El Niño and La Niña cancel each other in the long-term, it makes little sense to talk about El Niño accelerating or decelerating global warming. However, it is meaningful to wonder whether global warming could amplify El Niño.

Saravan (2023) Who’s afraid of the Big Bad El Niño? (with links to Google Colab Code)

4.4 Earth Energy Imbalance (EEI)

Loeb Abstarct

Earth’s Energy Imbalance (EEI) is a relatively small (presently ∼0.3%) difference between global mean solar radiation absorbed and thermal infrared radiation emitted to space. EEI is set by natural and anthropogenic climate forcings and the climate system’s response to those forcings. It is also influenced by internal variations within the climate system. Most of EEI warms the ocean; the remainder heats the land, melts ice, and warms the atmosphere. We show that independent satellite and in situ observations each yield statistically indistinguishable decadal increases in EEI from mid-2005 to mid-2019 of 0.50±0.47 W m-2 decade-1 (5%-95% confidence interval). This trend is primarily due to an increase in absorbed solar radiation associated with decreased reflection by clouds and sea-ice and a decrease in outgoing longwave radiation (OLR) due to increases in trace gases and water vapor. These changes combined exceed a positive trend in OLR due to increasing global mean temperatures.

Figure: Comparison of overlapping one-year estimates at 6-month intervals of net top-of-the-atmosphere annual energy flux from the CERES EBAF Ed4.1 product (solid red line) and an in situ observational estimate of uptake of energy by Earth climate system (solid blue line). Dashed lines correspond to least squares linear regression fits to the data

Loeb (2021) Satellite and Ocean Data Reveal Marked Increase in Earth’s Heating Rate (pdf) NASA Guardian

Hansen

EEI is the net gain (or loss) of energy by the planet, i.e., the difference between absorbed solar energy and emitted thermal (heat) radiation. EEI is the basic diagnostic informing us where global climate is headed. As long as more energy is coming in than going out, i.e., as long as EEI is positive, Earth will continue to get hotter.

EEI is hard to measure, a small difference between two large quantities (Earth absorbs and emits about 240 W/m2 averaged over the entire planetary surface), but the change of EEI is well-measured from space. Absolute calibration is obtained from the change of heat in the heat reservoirs, mainly the global ocean, over a period of at least a decade, as required to reduce error due to the finite number of places that the ocean is sampled.

EEI varies from year-to-year (Fig. 1), largely because global cloud amount varies with weather and ocean dynamics, but averaged over several years, EEI tells us what is needed to stabilize climate.

Fig: 12-month running-mean of Earth’s energy imbalance, based on CERES satellite data for EEI change normalized to 0.71 W/m2 mean for July 2005 – June 2015 based on in situ data.

10 years ago, EEI was about 0.6 W/m2, averaged over six years (that may not sound like much, but it equals the energy in 400,000 Hiroshima atomic bombs per day, every day). Now, it appears, EEI has approximately doubled, to more than 1 W/m2.

EEI has large natural variability on time scales of a few years and even longer time scales, mainly because clouds are sensitive to fluctuations of weather patterns and ocean dynamics. Thus, we need more time to access the magnitude and significance of the increase of EEI over the past decade.

Accelerated Warming

However, we can say with confidence that EEI increased in the past decade, an increase consistent with independent evidence that the growth rate of human-made climate forcing increased. Increased EEI is the basis for our projection that global warming will accelerate by as much as 50-100% in the few decades following 2010

We are using EEI as a crude, indirect, way to assess changes in global aerosols, and the climate forcing caused by aerosols that occurs principally via effects on clouds. It would be much better if we also had more direct measurement of aerosol.

Hansen (2022) Earth’s Energy Imbalance and Climate Response Time

4.4.1 Energy Imbalance Trend (TEEI)

Raghuraman Abstract

The observed trend in Earth’s energy imbalance (TEEI), a measure of the acceleration of heat uptake by the planet, is a fundamental indicator of perturbations to climate. Satellite observations (2001–2020) reveal a significant positive globally-averaged TEEI of 0.38 ± 0.24 \(Wm^{-2} decade^-1\), but the contributing drivers have yet to be understood. Using climate model simulations, we show that it is exceptionally unlikely (<1% probability) that this trend can be explained by internal variability. Instead, TEEI is achieved only upon accounting for the increase in anthropogenic radiative forcing and the associated climate response. TEEI is driven by a large decrease in reflected solar radiation and a small increase in emitted infrared radiation. This is because recent changes in forcing and feedbacks are additive in the solar spectrum, while being nearly offset by each other in the infrared. We conclude that the satellite record provides clear evidence of a human-influenced climate system.

Raghuraman Memo

Earth’s energy imbalance (EEI) is the difference between the incoming solar radiation (S 0 ), and the reflected shortwave radiation (RSW) plus the outgoing longwave radiation (OLR), at the top of the atmosphere. Thus, EEI is a fundamental measure of the degree to which the Earth’s global climate system is out of balance. A trend in EEI measures the acceleration of heat uptake by the planet and hence is an indicator of perturbations to the coupled atmosphere-ocean-land-ice system. The mean EEI during 2005–2015 was estimated to be 0.71 ± 0.10 \(Wm^{-2}\) . Oceans store 90% of this excess heat. Because of this close relationship between EEI and ocean heating, EEI trends have an important bearing on the warming of the global climate system, sea-level rise, and marine health.

In the contemporary climate system, internal variability (ε), effective radiative forcing (ΔERF), and the radiative response (λΔT s ) change EEI. Thus, an anomaly in EEI can be expressed as the sum of these three terms.

\[ΔEEI = ΔERF + λΔT_s+ ε\]

Climate feedbacks and surface temperature are represented by λ and \(T_s\) respectively.

While we know that TEEI is influenced by internal variability in the climate system, external forcings, and climate feedbacks, the extent of the contribution from each has not been previously determined.

In particular, due to the inherent noise in the Earth system, a single observed 20-year time series of EEI is only one of many possible time series that internal variability could produce, and therefore, it is imperative to quantify the contribution of internal variability (ε) to TEEI.

We focus here on the investigation of the trend in EEI (TEEI). In particular, the contributions to TEEI by the drivers of the OLR and RSW trends, are an unexplored and unquantified area in the explanation of the observed satellite record.

We use Coupled Model Intercomparison Project Phase 6 (CMIP6) experiments and design a hierarchy of climate model experiments with the Geophysical Fluid Dynamics Laboratory Coupled/Atmospheric Model 4.0 (GFDL CM4/AM4)to better understand the contributions of the three components of Eq. (1) to the CERES- observed TEEI, thereby providing an assessment of the relative importance of anthropogenically induced changes versus internal variability changes. We show that anthropogenic forcing and the associated climate response yield the observed positive TEEI.

The observed ΔEEI time series yields a significant positive trend of 0.38 ± 0.24 \(Wm^-2 decade^-1\) There is, however, the potential for systematic errors associated with the observed trend due to instrument drift. The trend is driven by large reductions in RSW, compensated only by a relatively weaker increase in OLR.

The observed trend in EEI is unexplained by internal varia- bility. Our hierarchy of modeling experiments allows us to investigate the possible contributions of trends in ΔERF, \(λΔT_s\) , and ε to TEEI. We compose five sets of estimates of the internal variability (ε) using climate model simulations. For all five sets of estimates, we calculate ε as the ±2σ range of 20-year trends across the realizations.

We find that the CERES-observed TEEI lies outside the range of trends driven by internal variability alone.

ε comprises a range of trends in EEI due to internal variability, the value of ε, ±0.19 \(Wm^{-2} decade^−1\)

The satellite-observed positive EEI trend over the 2001–2020 period is exceptionally unlikely (<1% probability) to be explained by internal variability.

Raghuraman (2021) Anthropogenic forcing and response yield observed positive trend in Earth’s energy imbalance (pdf)

NBCnews

Stability in Earth’s climate hinges on a delicate balance between the amount of energy the planet absorbs from the sun and the amount of energy Earth emits back into space. But that equilibrium has been thrown off in recent years — and the imbalance is growing, according to a paper published Wednesday in the journal Nature Communications.

The Princeton University researchers behind the paper found that there’s a less than 1 percent probability that the changes occurred naturally.

The findings undercut a key argument used by people who do not believe human activity is responsible for the bulk of climate change to explain trends in global warming, demonstrating that the planet’s energy imbalance cannot be explained just by Earth’s own natural variations.

nbcnews

4.5 Hurricanes

Masters

Hurricanes are heat engines that take heat energy out of the ocean and convert it to the kinetic energy of their winds. A hotter ocean will allow hurricanes to grow more powerful, assuming that the other factors that power hurricanes, including low wind shear and a moist atmosphere, are present.

As far back as 1987, MIT hurricane scientist Kerry Emanuel theorized that the wind speeds in hurricanes can be expected to increase about 5% for every one degree Celsius (1.8°F) increase in tropical ocean temperature, assuming that the average wind speed near the surface of the tropical oceans does not change. Computer modeling has found a slightly smaller magnitude (4%) for the increase.

Damage from a hurricane increases exponentially with an increase in winds.

A Category 2 hurricane with 100 mph winds will do 10 times the damage of a Category 1 hurricane with 75 mph winds. This includes damage not only from winds, but also from storm surge, inland flooding, and tornadoes. Bottom line: A 4-5% increase in winds yields about a 40-50% increase in the destructive potential of a hurricane.

Fig: Damage multiplier for hurricane winds compared to a minimal Category 1 hurricane with 75 mph winds. A doubling of the wind speed from 75 mph to 150 mph increases potential damage by a factor of 256. (Image credit: NOAA)

Increased sea surface temperatures of about 1 degree Celsius since preindustrial time in the Gulf of Mexico from human-caused global warming led to a 40-50% increase in Idalia’s destructive power by increasing the hurricane’s winds by at least 4-5%.

Masters (2023) Hotter ocean temperatures from global warming likely increased Idalia’s destructive power by at least 40-50%

GFDL

Summary Statement

Two frequently asked questions on global warming and hurricanes are the following:

What changes in hurricane activity are expected for the late 21st century, given the pronounced global warming scenarios from IPCC models?
Have humans already caused a detectable increase in Atlantic hurricane activity or global tropical cyclone activity?

The IPCC AR6 presents a strong body of scientific evidence that it is unequivocal that humans have caused the earth’s climate to warm, with a likely human contribution of 0.8 to 1.3 degrees Celsius to global mean temperature since the late 1800s. But what does this anthropogenic global warming mean for Atlantic hurricane activity, or global tropical cyclone activity? Here, we address these questions, starting with those conclusions where we have relatively more confidence. The main text of this web page gives more background discussion. “Detectable” change here will refer to a change that is large enough to be clearly distinguishable from the variability due to natural causes. Our main conclusions are: Likelihood Statements

The terminology here for likelihood statements follows these conventions for the assessed likelihood of an outcome or result: Very Likely: > 90%, Likely: > 66% More Likely Than Not (or Better Than Even Odds) > 50%

Sea level rise – which human activity has very likely been the main driver of since at least 1971 according to IPCC AR6 – should be causing higher coastal inundation levels for tropical cyclones that do occur, all else assumed equal.
Tropical cyclone rainfall rates are projected to increase in the future (medium to high confidence) due to anthropogenic warming and accompanying increase in atmospheric moisture content. Modeling studies on average project an increase on the order of 10-15% for rainfall rates averaged within about 100 km of the storm for a 2 degree Celsius global warming scenario.
Tropical cyclone intensities globally are projected to increase (medium to high confidence) on average (by 1 to 10% according to model projections for a 2 degree Celsius global warming). This change would imply an even larger percentage increase in the destructive potential per storm, assuming no reduction in storm size. Storm size responses to anthropogenic warming are uncertain.
The global proportion of tropical cyclones that reach very intense (Category 4 and 5) levels is projected to increase (medium to high confidence) due to anthropogenic warming over the 21st century. There is less confidence in future projections of the global number of Category 4 and 5 storms, since most modeling studies project a decrease (or little change) in the global frequency of all tropical cyclones combined.

For the above tropical cyclone projections, the IPCC AR6 generally concluded there was high confidence as compared to medium-to-high confidence in the WMO assessment. Additional research was published between the reports, which can affect confidence levels. Also, confidence levels for assessment statements can vary between authors within a given report.

In terms of detection and attribution, much less is known about human contributions to hurricane/tropical cyclone activity changes to date, compared to the case for global mean temperature.  Recent findings include:
    In the northwest Pacific basin, observations show a poleward shift in the latitude of maximum intensity of tropical cyclones.  This change is assessed to be detectable (i.e., not explainable by internal variability alone) with medium confidence (IPCC AR6)  or low-to-medium confidence (WMO Task Team report).
     One study finds an increase in the fraction of tropical cyclone intensity estimates of at least Category 3 intensity both globally and in the Atlantic basin, over the past four decades.  These observed changes have not been confidently attributed to anthropogenic forcing.  The global increase was assessed to be detectable (i.e., not explainable by internal variability alone) with medium confidence by IPCC AR6.
    A study of rapid intensification of Atlantic hurricanes finds an observed increase in the probability of rapid intensification (1982-2009) which is highly unusual compared to one climate model’s simulation of internal multidecadal climate variability, and thus is a possibly emerging detectable anthropogenic change.  The increase is consistent in sign with the model’s simulated response to anthropogenic forcing.
    There is increasing evidence from modeling studies at GFDL/NOAA and the UK Met Office/Hadley Centre (UKMO) that the increase in tropical storm frequency in the Atlantic basin since the 1970s has been at least partly driven by decreases in aerosols from human activity and volcanic forcing.  Natural variability or changes in Saharan dust emissions may also have contributed to recent changes.  The recent GFDL and UKMO studies do not imply that the increase in Atlantic tropical storm frequency since the 1970s will continue into the future:   these same models project future decreases in Atlantic tropical storm frequency in response to increasing greenhouse gas concentrations.
    There is evidence for a slowing of tropical cyclone propagation speeds over the continental U.S. over the past century, but these observed changes have not yet been confidently linked to anthropogenic climate change.
    There is no strong evidence of century-scale increasing trends in U.S. landfalling hurricanes or major hurricanes, Similarly for Atlantic basin-wide hurricanes (after adjusting for observing capabilities), there is not strong evidence for an increase since the late 1800s in hurricanes, major hurricanes, or the proportion of hurricanes that reach major hurricane intensity.
In summary, it is premature to conclude with high confidence that human-caused increases in greenhouse gases have caused a change in past Atlantic basin hurricane activity that is outside the range of natural variability, although greenhouse gases are strongly linked to global warming.  Some possible emerging human influences on past tropical cyclone activity were summarized above.  These include, for the Atlantic, recent increases in rapid intensification probability, aerosol-driven changes in hurricane activity, and increases in extreme precipitation in some regions.  At the global scale, increased intensities and fraction of tropical cyclone observations at high intensity are examples, along with a poleward shift of the latitude of maximum tropical cyclone intensity in the Northwest Pacific basin.  Human activities may have already caused other changes in tropical cyclone activity that are not yet clearly apparent due to the small magnitude of these changes compared to estimated natural variability, or due to observational limitations.

GFDL (2022) Global Warming and Hurricanes An Overview of Current Research Results

4.6 Jet Stream

Rousi

Accelerated western European heatwave trends linked to more-persistent double jets over Eurasia

Persistent heat extremes can have severe impacts on ecosystems and societies, including excess mortality, wildfires, and harvest failures. Here we identify Europe as a heatwave hotspot, exhibiting upward trends that are three-to-four times faster compared to the rest of the northern midlatitudes over the past 42 years. This accelerated trend is linked to atmospheric dynamical changes via an increase in the frequency and persistence of double jet stream states over Eurasia. We find that double jet occurrences are particularly important for western European heatwaves, explaining up to 35% of temperature variability. The upward trend in the persistence of double jet events explains almost all of the accelerated heatwave trend in western Europe, and about 30% of it over the extended European region. Those findings provide evidence that in addition to thermodynamical drivers, atmospheric dynamical changes have contributed to the increased rate of European heatwaves, with implications for risk management and potential adaptation strategies.

Rousi (2022) Accelerated western European heatwave trends linked to more-persistent double jets over Eurasia

4.6.1 Quasi-Resonant Amplification (QRA)

Mann Abstract

Persistent episodes of extreme weather in the Northern Hemisphere summer have been associated with high- amplitude quasi-stationary atmospheric Rossby waves, with zonal wave numbers 6 to 8 resulting from the phe- nomenon of quasi-resonant amplification (QRA). A fingerprint for the occurrence of QRA can be defined in terms of the zonally averaged surface temperature field. Examining state-of-the-art [Coupled Model Intercomparison Project Phase 5 (CMIP5)] climate model projections, we find that QRA events are likely to increase by ~50% this century under business-as-usual carbon emissions, but there is considerable variation among climate models. Some predict a near tripling of QRA events by the end of the century, while others predict a potential decrease. Models with amplified Arctic warming yield the most pronounced increase in QRA events. The projections are strongly dependent on assumptions regarding the nature of changes in radiative forcing associated with anthropo- genic aerosols over the next century. One implication of our findings is that a reduction in midlatitude aerosol loading could actually lead to Arctic de-amplification this century, ameliorating potential increases in persistent extreme weather events.

Mann (2018) Projected changes in persistent extreme summer weather events: The role of quasi-resonant amplification (pdf)

4.7 Total Solar Irradiance

Climate.gov: Incoming Sunlight

the Sun likely added at most 0.01 degrees Celsius to the 0.95–1.2 degrees Celsius of global surface warming that’s occurred since pre-industrial times.

Satellite observations through several solar cycles reveal that the difference in total average brightness between solar maxima and minima is very small, on the order of 1 Watt per square meter during strong cycles. On average, the Sun delivers 1,361 Watts of power per square meter at a distance of one astronomical unit. This amount is known as the total solar irradiance. Based on observations and models, experts estimate that the impact of this 11-year variation on global surface temperature is likely around 0.1 degrees Celsius or less.

  • The Sun’s overall brightness varies on timescales from minutes to millennia, and these changes are detectable in the global temperature record.
  • During strong solar cycles, the Sun’s total average brightness varies by up to 1 Watt per square meter; this variation affects global average temperature by 0.1 degrees Celsius or less.
  • Changes in the Sun’s overall brightness since the pre-industrial period have been minimal, likely contributing no more than 0.01 degrees Celsius to the roughly 1 degree of warming that’s occurred over the Industrial period.
  • Projected warming due to increasing greenhouse gas levels in the coming decades will overpower even a very strong Grand Solar Minimum.
  • Rising amounts of atmospheric carbon dioxide have postponed the next Milankovitch-driven ice age by at least tens of thousands of years.

Figure: Yearly total solar irradiance (orange line) from 1610–2020 and the annual global temperature compared to the 20th-century average (red line) from 1880–2020. Since the middle of the 20th-century, solar activity has declined while global temperature increased rapidly. NOAA Climate.gov image, based on solar data from Coddington et al., 2016, and temperature data from NOAA NCEI.

Even if the Sun’s recent quietness—the 11-year cycle minimum in 2011 was the lowest in a century—were to turn into a multi-decade stretch of extremely low activity known as a Grand Solar Minimum, it wouldn’t overpower the amount of global warming projected for the coming century due to increasing greenhouse gas emissions. In fact, as long as atmospheric carbon dioxide remains above 300 parts per million, not even the next ice age, which Milankovitch theory predicts would begin 50,000 years from now, is likely to occur.

Figure: Daily observations of total solar irradiance (orange line) since the start of the satellite era in 1978. Day-to-day, TSI may vary by as much as 0.3 percent, but average differences between maximum and minimum are on the order 0.1 percent, or around 1 Watt per square meter. NOAA Climate.gov image, based on data from LASP Interactive Solar Irradiance Data Center.

Scientists today have close to four decades of overlapping measurements of total solar irradiance and sunspots, which allow them to statistically describe how changes in sunspot numbers relate to variations in total solar irradiance. They’ve used that relationship to model the Sun’s brightness back to the start of the sunspot record in the 1600s.

Figure: Total solar irradiance estimated from sunspot observations since 1610. NOAA Climate.gov image, based on the NOAA TSI Climate Data Record (Coddington et al., 2016).

These historical reconstructions reveal that some solar cycles are more active than others, and that their timing isn’t completely random. Generally, 2-3 relatively strong cycles will be preceded and followed by 2-3 relatively weak ones. These alternating strong and weak epochs tend to be grouped together over approximately 100-year periods, a pattern known as Gleissberg cycles. Over the span of the historical sunspot record, there have been three ~100-year Gleissberg cycles: 1700-1810, 1810-1910, and 1910-2010.

Figure: The amplitude of the 11-year solar cycle (formally called the Schwabe cycle, orange) is modulated by the approximately 100-year Gleissberg cycle (charcoal), in which a number of consecutive cycles of high activity are bracketed by consecutive cycles of lower activity. The start of the Industrial Revolution in the mid-1700s coincided with a Gleissberg maximum. The relatively high activity of the mid 20th-century also coincided with a Gleissberg maximum, while the recent decades coincide with a Gleissberg minimum. NOAA Climate.gov image based on data from Wang and Lean, 2021.

In addition, the record shows that there have been periods when sunspots virtually disappear for several decades. (Other features of the 11-year solar cycle continue to occur, however.) These periods are called Grand Solar Minimums. For example, between 1645-1715, the Sun went through a 70-year quiet period known as the Maunder Minimum. Sunspots disappeared almost completely, and the solar wind was maybe half of its modern velocity. The Maunder Minimum partially overlapped a centuries-long cold spell called the Little Ice Age, which was strongest in the Northern Hemisphere between 1450­-1850.

Figure: A 2000-year temperature history of the Northern Hemisphere outside the tropics shows a warm period that peaked around 1,000 A.D. followed by a multi-century period of cooling: the Little Ice Age. The coldest part of the Little Ice Age overlapped the very low solar activity of the Maunder Minimum, but the cold spell began well before. NOAA Climate.gov graph, based on data from Christiansen and Ljungqvist, 2012.

In contrast, the Sun was unusually active in the twentieth century, a period which solar experts call the Modern Maximum. Starting near the turn of the twentieth century, each solar cycle was increasingly active. This build up was tied to the last Gleissberg Cycle, which peaked during solar cycle 19 in 1957. Solar activity then declined in the second half of the 20th-century. The stretch of high activity drew to a definite close in the first decade of the twenty-first century with solar cycle 23, which had an unusually long and low minimum. Solar cycle 24 went on to have one of the lowest maximums of the last 70 years, and solar cycle 25 is expected to be comparable. Meanwhile, Earth’s surface temperatures continued to rise rapidly.

The increasing solar activity of the first half of the 20th century and the decreasing activity since then have largely canceled each other out in terms of their influence on global temperature.

The modern sunspot record tells us about solar activity over the past four centuries. Indirect evidence for solar activity deeper in the past comes from the presence of cosmogenic isotopes—radioactive atoms that are generated when common isotopes of an element are struck by galactic cosmic rays.

Our solar system is constantly bombarded with galactic cosmic rays, but the Sun’s magnetic field shields us from most of them. When the Sun’s magnetic field is strong, at solar maximum, fewer cosmic rays reach the atmosphere, creating very few cosmogenic isotopes. At solar minimum, when the Sun’s magnetic field is weaker, slightly more cosmic rays reach Earth’s atmosphere, generating more cosmogenic isotopes. The two most common cosmogenic isotopes are carbon-14, which can be found in tree rings, and beryllium-10, which is found in ice cores. Using fluctuations in cosmogenic isotopes, experts have reconstructed solar activity back thousands of years.

Figure: A reconstruction of total solar irradiance over the past 9,400 years based on a combination of carbon-14 isotopes in tree rings and beryllium-10 in ice cores. The record indicates there have been at least 25 Grand Solar Minimums in the Holocene. NOAA Climate.gov image, based on data from Steinhilber et al., 2012.

These paleoclimate reconstructions reveal that the Sun has produced at least 25 grand minimums in the past 9,000 years. Some are short—just two or three decades—and others, like the Maunder Minimum, are five or more decades. They occur every 200 years or so, a period known as the de Vries cycle. Many of them were preceded by a solar cycle with an unusually long and low solar minimum, similar to the minimum of 2008.

Milankovitch cycles and ice ages

The 11-year sunspot cycle and its Gleissberg-cycle modulation cause small changes in the Sun’s actual brightness—how much sunlight the Sun radiates to Earth. Earth’s climate is also affected by how much sunlight reaches us due to changes in our planet’s orbit and position in space relative to the Sun. Called Milankovitch cycles, these predictable orbital patterns have repeat times of tens to hundreds of thousands of years.

For the past million years at least, Milankovitch cycles have coincided with 100,000-year-long ice ages punctuated by short intervals of rapid warming. Although there are pieces of the puzzle experts still don’t understand, the key climate influence seems to be changes in the amount of incoming sunlight, or insolation, reaching the high latitudes of the Northern Hemisphere during the summer. The Northern Hemisphere is key to the ice ages because massive ice sheets can only grow over land, not ocean, and most of Earth’s land area has been concentrated in the Northern Hemisphere for at least tens of millions of years.

Figure: Changes in Northern Hemisphere sunlight are the key to ice ages because massive ice sheets, such as the one that spread across North America and Eurasia roughly 20,000 years ago, can only grow over land, and most land area has been in the Northern Hemisphere for tens of millions of years. (High-resolution without annotations available.) Image by Climate.gov based on data from Science on a Sphere.

The most significant changes in Northern Hemisphere insolation come from three variations in Earth’s orbit:

precession (~26,000 years): the slow rotation or “wobble” in the Earth’s axis of rotation, which changes where in the annual orbital path Northern Hemisphere summer solstice occurs;
obliquity (~41,000 years): how tilted Earth’s axis of rotation is;
eccentricity (~100,000 years): how far Earth’s orbit is from being a perfect circle.

Because these cycles have different lengths, they overlap in complex rhythms, reinforcing one another at some times and offsetting one other at others. Northern Hemisphere summer insolation is maximized when tilt is extreme, eccentricity is extreme, and precession causes Northern Hemisphere summer solstice to occur near perihelion, the place in its orbit when Earth is closest to the Sun. Summer insolation is minimized when tilt is smaller, eccentricity is extreme, and Northern Hemisphere summer solstice occurs near aphelion, when Earth is farthest from the Sun.

Figure: Rows 1-3) Milankovitch cycles over the past million years (tilt, precession, and eccentricity. (Row 4) Northern Hemisphere summer insolation. (Row 5) Global ice volume inferred from oxygen isotopes in sea floor sediments. Light gray column highlights conditions around 220,000 years ago, when overlap among the three orbital cycles brought a peak in Northern Hemisphere insolation, triggering a warming period with low ice sheet volume. NOAA Climate.gov image, based on data from Lisiecki and Raymo, 2005.

How low summer insolation must fall to trigger an ice age depends on how high atmospheric carbon dioxide levels are; the more carbon dioxide, the lower the insolation must be. Looking back over the past million years, the highest carbon dioxide level at the start of any ice age was 300 ppm, and most were far lower.

At pre-industrial carbon dioxide levels, summer insolation at 65° North need only dip 0.75 standard deviations below the mean—about 15 Watts per square meter—for summers to be too cool to melt all the winter snow, a low that Milankovitch cycles predict we will next hit about 50,000 years from now. At 400 parts per million, summer insolation would need to fall twice as much—a low we will next see 125,000 years from now. At carbon dioxide levels above 560 parts per million no Milankovitch variation within the next half million years will be low enough to trigger an ice age.

Climate.gov: Incoming Sunlight

Danish ISAC Study

Conclusions

The Sun - Sunspot number catches only part of solar activity - The Sun’s magnetic field drives short term changes - Use therefore magnetic data for Sun-Climate studies - Solar output predictability is at most years ahead - Due to the complexity and time-varying nature of the Sun solar-climate coupling must be expected to be very complex and probably changing with time

The total solar irradiance (TSI) is in a sense the most obvious solar parameter to single out for study when considering the Sun-climate connection. The power input received by the Earth from the Sun is the driver of Earth’s climate system and the variations in TSI do correlate reasonably well with the variations seen in the climate. However, the brute variations in power received at Earth are not strong enough to explain the variations in climate that are observed. For TSI to be the driver of the climate variations observed some sort of amplification or feedback mechanism must therefore be in place.

The influence of TSI on the climate is multifaceted. Variations in TSI translate to different heating in various layers of the atmosphere (though UV here plays the major role), heat absorption in the oceans with different time scales from those of the atmosphere, differential heating over land and oceans, etc. These primary influences from TSI will in turn alter such things as atmospheric circulation, water vapour content in the atmosphere, and cloud cover. Effects from these secondary phenomena may then feed back into the energy absorbed by the climate system from solar irradiation, for example via changes in the Albedo, possibly enhancing the effects of TSI variations. Due to the geographically heterogeneous nature of solar influence on climate and the rather complicated feedback mechanisms involved in solar forcing of the climate any simple explanation of TSI-climate connections is likely to fail, except in limited cases such as variations in sea surface temperature that seem to be at least partly explainable by simple energy balance considerations.

Danish national Space Center (2007) Influence of Solar Activity Cucles on Earth’s Climate (pdf) Presentation (pdf)

4.8 Spatial Shock

4.8.1 Coastal Flooding

Desmet

Just moving to higher grounds

Nobel Prize winner William Nordhaus has called climate change “the ultimate challenge for economics.”

Economists increasingly have been trying to understand how rising tides and global temperatures will impact resource allocation around the globe, as well as the potential policy tools that can help curb damage to our natural world.

SMU professor Klaus Desmet says that a lot of those analyses are missing a critical factor: migration.

Desmet coauthored a paper in the American Economic JournalEvaluating the Economic Cost of Coastal Flooding (paywall): Macroeconomics that examines how economic output will be affected over the next 200 years as humans move away from coastal areas threatened by rising sea levels. Although losses in vulnerable Southeast Asian cities such as Bangkok and Shanghai will still be very significant, their research shows that overall GDP declines are substantially less than predicted by models that don’t account for spatial shifts in economic activity.

Climate change is to a large extent a spatial shock.

Migration as one of the key responses to climate change.

What we find in our model is that at a global level, flooding decreases real GDP by about 0.1 percent by the year 2200. If we were to completely ignore the dynamic spatial response of the economy, if we were to have everyone stay put in the face of rising seas, then the loss would actually increase to 4.5 percent. So, that difference between 0.1 and 4.5 percent underscores the first-order importance of taking into account moving, migration, and the spatial dynamic response of the economy to rising sea levels.

We develop a high resolution dynamic model of the world economy. This model splits up the world into 64,800 1° by 1° grid cells, which are linked to each other through trade and migration. We feed into this model high-quality projections of both global and local sea level rise over the next 200 years. . . . When we run our model forward, we can then assess what the economic effect will be of having these pieces of land lost for production.

We find large losses in coastal areas of south and east Asia, with countries such as Vietnam, Thailand, and Bangladesh losing up to 10 percent of real GDP in present discounted value terms over the next 200 years. Other areas that will also suffer disproportionately include coastal Northwest Europe, [and] some areas on the US East Coast and Gulf Coast. By contrast, the Pacific Coast of the Americas is much less affected and, in fact, most of the coastal areas of Africa are also a lot less affected.

Desmet