Solar Forcing of Regional Climate Change During the Maunder Minimum

See allHide authors and affiliations

Science  07 Dec 2001:
Vol. 294, Issue 5549, pp. 2149-2152
DOI: 10.1126/science.1064363


We examine the climate response to solar irradiance changes between the late 17th-century Maunder Minimum and the late 18th century. Global average temperature changes are small (about 0.3° to 0.4°C) in both a climate model and empirical reconstructions. However, regional temperature changes are quite large. In the model, these occur primarily through a forced shift toward the low index state of the Arctic Oscillation/North Atlantic Oscillation as solar irradiance decreases. This leads to colder temperatures over the Northern Hemisphere continents, especially in winter (1° to 2°C), in agreement with historical records and proxy data for surface temperatures.

A minimum in solar irradiance, the Maunder Minimum, is thought to have occurred from the mid-1600s to the early 1700s (1–3). Concurrently, surface temperatures appear to have been at or near their lowest values of the past millennium in the Northern Hemisphere (NH) (4–7), and European winter temperatures were reduced by 1° to 1.5°C (8).

We used a version of the Goddard Institute for Space Studies (GISS) general circulation model (GCM), which includes a detailed representation of the stratosphere, to simulate the difference between that period and a century later, when solar output remained relatively high over several decades. The model contains a mixed-layer ocean, allowing sea surface temperatures (SSTs) to respond to radiative forcing, and has been shown to capture observed wintertime solar cycle–induced variations reasonably well (9).

The GCM includes parameterizations of the response of ozone to solar irradiance, temperature, and circulation changes (although transport changes are noninteractive) (9), based on results for preindustrial conditions from our two-dimensional (2D) chemistry model (10). These show more ozone in the upper and lower stratosphere, as expected (11), primarily because of the absence of anthropogenic halogens and a drier stratosphere (12). Ozone's temperature sensitivity was also increased throughout the upper stratosphere, more than doubling at 3 millibars (mbar) relative to the present (13). Solar heating then speeds up chemical destruction more than irradiance enhances photolytic production, so that ozone concentration decreases with solar output. Thus, the reduced Maunder Minimum irradiance leads to increased upper stratospheric ozone, which in turn causes lower stratospheric ozone to decrease. Both cause negative radiative forcing (14). Inclusion of ozone photochemistry therefore enhances radiative forcing from preindustrial solar variations.

Equilibrium simulations were performed for spectrally discriminated irradiances in 1680 and 1780 (1). An alternative reconstruction has a 40% larger long-term solar variation, which would give a larger climate response (15). Initial conditions were taken from 1680 and 1780 in a transient simulation. A spin-up time of about 25 years was nevertheless required for equilibration with the increased preindustrial ozone. Results are therefore based on the last 30 years of 60-year simulations.

Modeled global annual average surface temperatures were 0.34°C cooler in 1680 than in 1780. This is similar to the annual average NH 1680-minus-1780 temperature change of about –0.2° ± 0.2°C from proxy data (6) or an energy balance model (7) (model NH cooling is –0.33° C). Other proxy-based reconstructions for the NH extratropics during summer indicate a –0.3° to –0.4°C cooling (4, 5) (comparable model cooling is –0.3°C).

Modeled surface temperature changes show alternating warm oceans and cold continents at NH mid-latitudes (Fig. 1), with maximum amplitude in winter. To characterize regional changes, we calculated empirical orthogonal functions (EOFs) from the NH extratropical (20° to 90°N) cold season (from November to April) sea level pressure (SLP) time series. The last 40 years of each run were used as an 80-year control.

Figure 1

Surface temperature change (in °C) from 1780 to 1680 in the GCM. NH annual average (left) and November to April NH extratropics (right) are shown. Nearly all points are statistically significant (not shown) because of the large number of model years. Data for all figures are available as Online supplemental material (40).

The leading control EOF (Fig. 2) closely matches the Arctic Oscillation pattern (AO, also called the Northern Hemisphere Annular Mode), defined as the first EOF of 20th-century SLP (16). This pattern contains the North Atlantic Oscillation (NAO), which may be considered a different view of the same phenomenon (17). We define the model's AO pattern as EOF 1, and the AO index as the opposite of the area-weighted change in SLP northward of 60°. Projection of the 1780-to-1680 SLP difference (Fig. 2) onto the first five control EOFs shows that they account for 37.3, 0.0, 3.4, 1.5, and 0.0% of the variance, respectively. Thus, the regional climate variations primarily result from a lower index state of the AO during the Maunder Minimum relative to that of a century later (by –1.1 mbar; similarly, the SLP contrast between Iceland and the Azores decreases by 1.9 mbar). This reduces onshore advection of warm oceanic air, yielding a 1° to 2°C cooling over large areas of North America and Eurasia (Fig. 1).

Figure 2

The leading EOF of NH extratropical November-to-April SLP over the last 40 years of the 1680 and 1780 simulations (left), and the SLP change from 1780 to 1680 (right), both in mbar (mb). The SLP difference is filtered in EOF space by showing the projection onto the first 20 EOFs, which contain 70% of the total difference. This removes some of the high-frequency noise in the model induced by overly large energies at high wavenumbers. Hatched areas indicate statistical significance at the 90% level.

Historical data associate the European Maunder Minimum winter cooling with enhanced northeasterly advection of continental air, consistent with an anomalous negative NAO (18). Reconstructions of large-scale surface temperature patterns in past centuries (19, 20) from networks of diverse proxy information (tree rings, ice cores, corals, and historical records) are particularly useful benchmarks for comparison. Correlations between these reconstructions and solar irradiance reconstructions (1) during the preanthropogenic interval from 1650 to 1850 provide an empirical estimate of the large-scale climate response to solar forcing (21). The instantaneous (zero-lag) frequency-independent regression pattern, primarily indicative of decadal solar-climate relationships, shows little evidence of an AO/NAO signal (Fig. 3, top). This is consistent with the extremely weak modeled AO/NAO response to the 11-year solar cycle using fixed SSTs, equivalent to an instantaneous response (22). However, when filtered to isolate the multidecadal-to-centennial time scales associated with the Maunder Minimum, and lagged to allow for inertia in the ocean's response, the empirical regression (Fig. 3, bottom) shows a clear AO/NAO-type pattern of alternating cold land and warm ocean temperature anomalies. Both paleoclimate reconstructions and the GCM thus indicate in a remarkably consistent manner that solar forcing affects regional scales much more strongly than global or hemispheric scales through forcing of the AO/NAO.

Figure 3

Regression during the period from 1650 to 1850 between reconstructed solar irradiance (1) and annual average surface temperatures (19). The top panel shows the instantaneous regression over all time scales, whereas the bottom panel shows the regression when filtered to only include contributions from time scales longer than 40 years, with a 20-year lag (results are similar for time lags of 10 to 30 years). The correlation is given in °C per –0.32 W/m2 (the solar forcing used in the model) for comparison with the modeled temperature changes shown in Fig. 1. Gray areas indicate no data; hatched areas indicate statistical significance at the 90% level.

Most additional paleoclimate data have insufficient time resolution to distinguish the Maunder Minimum but suggest a shift toward the AO/NAO low-index state during periods of reduced solar forcing. Ocean sediments show a decreased SST of 1° to 2°C in the Sargasso Sea during the 16th to 19th centuries, and warmer Atlantic temperatures north of 44°C, extending to the area off Newfoundland (23), which is consistent with a reduced NAO but not with uniform basin-wide cooling. Other evidence suggests colder North Atlantic temperatures during those centuries, however, including evidence for increased sea ice around Iceland (24), a region of minimal change in the GCM.

Our previous studies have demonstrated how external forcings can excite the AO/NAO in the GISS GCM (22, 25). Briefly, the mechanism works as follows, using a shift toward the high-index AO/NAO as an example: (i) tropical and subtropical SSTs warm, leading to (ii) a warmer tropical and subtropical upper troposphere via moist convective processes. This results in (iii) an increased latitudinal temperature gradient at around 100 to 200 mbar, because these pressures are in the stratosphere at higher latitudes, and so do not feel the surface warming (26). The temperature gradient leads to (iv) enhanced lower stratospheric westerly winds, which (v) refract upward-propagating tropospheric planetary waves equatorward. This causes (vi) increased angular momentum transport to high latitudes and enhanced tropospheric westerlies, and the associated temperature and pressure changes corresponding to a high AO/NAO index. Observations support a planetary wave modulation of the AO/NAO (27,28), and zonal wind and planetary wave propagation changes over recent decades are well reproduced in the model (22).

Reduced irradiance during the Maunder Minimum causes a shift toward the low-index AO/NAO state via this same mechanism. During December to February, the surface in the tropics and subtropics cools by 0.4° to 0.5°C because of reduced incoming radiation and the upper stratospheric ozone increase. Cooling in the tropical and subtropical upper troposphere is even more pronounced (∼0.8°C) because of cloud feedbacks, including an ∼0.5% decrease in high cloud cover induced by ozone through surface effects. A similar response was seen in simulations with a finer resolution version of the GISS GCM (14). This cooling substantially reduces the latitudinal temperature gradient in the tropopause region, decreasing the zonal wind there at ∼40°N. Planetary waves coming up from the surface at mid-latitudes, which are especially abundant during winter, are then deflected toward the equator less than before (equatorward Eliassen-Palm flux is reduced by 0.41 m2/s2, 12° to 35°N, 300 to 100 mbar average), instead propagating up into the stratosphere (increased vertical flux of 6.3 × 10−4m2/s2, 35° to 60°N, 100 to 5 mbar average) (29). This increases the wave-driven stratospheric residual circulation, which warms the polar lower stratosphere (up to 1°C), providing a positive feedback by further weakening the latitudinal temperature gradient. The wave propagation changes imply a reduction in northward angular momentum transport, hence a slowing down of the middle- and high-latitude westerlies and a shift toward the low AO/NAO index. Because the oceans are relatively warm during the winter owing to their large heat capacity, the diminished flow creates a cold-land/warm-ocean surface pattern (Fig. 1).

To quantify the importance of ozone changes, an identical pair of 1680 and 1780 experiments was run, leaving out the ozone response. These simulations show a reduction in the AO index of only 0.5 mbar, and the SLP change projection onto the AO was 15.1%. Thus, there was an analogous but weaker response without interactive ozone. The ozone changes further reduce the tropopause radiative forcing by –0.02 W/m2 (instantaneous, but with adjusted stratospheric temperatures) for 1680 conditions (compared to the direct solar forcing of –0.32 W/m2). Cloud feedbacks caused by ozone magnify its influence somewhat, so that the global annual average surface temperature change increases from –0.28° to –0.34°C, including ozone chemistry. The largest impacts resulted from the enhancement of the tropopause region latitudinal temperature gradient through high cloud changes and subsequent wave-driven feedbacks.

In contrast to the preindustrial simulations, 20th-century simulations (30) show that solar irradiance changes alone warm the surface by ∼0.25°C, but when ozone changes and climate feedbacks are included, the cooling drops to ∼0.19°C. With current atmospheric composition, ozone enhancement from increased irradiance outweighs temperature-induced changes, leading to more middle and upper stratospheric ozone (9), which dampens the net radiative forcing. Ozone's reversal from a positive (preindustrial) to a negative feedback supports results showing that solar forcing has been a relatively minor contributor to late 20th-century surface warming (7, 19, 31). However, ozone feedbacks remain important to the dynamical response associated with decadal and shorter term solar variability, which is initiated by upper stratospheric temperature changes, which are proportional to solar heating and hence to ozone (9).

The AO/NAO response is dependent on the overall forced change in tropical and subtropical SSTs, translation of these changes into upper tropospheric temperature changes by hydrological feedbacks, and the resulting wind and planetary wave refraction changes. These processes vary between models and in particular are sensitive to the inclusion of a realistic stratosphere (22, 25), emphasizing the need for intermodel comparisons of similar experiments (32).

The European Center Hamburg 3 GCM exhibits a similar spatial pattern of response to solar forcing, with warming over the NH continents and cooling over the North Atlantic and North Pacific (33). In that model, however, the surface temperature and AO sensitivity were much less than in our simulation. This is probably due to the lack of interactive ozone and the poorly resolved stratosphere, two important positive feedbacks in our model.

The importance of SST changes to initiate the AO/NAO response is consistent with our fixed-SST solar cycle experiments (9,22), whose response did not project strongly onto the AO/NAO, and with the weak empirical instantaneous climate sensitivity (Fig. 3, top). Present-day observations and modeling also support links between long-term tropical SST changes and the NAO (34, 35).

Surface temperatures in an energy balance model (36) driven by volcanic forcing gave a poor correlation with proxies during the 17th and 18th centuries. Correlations between calculations of solar-induced climate change and temperature proxies were quite good, however, although the impacts were very small. Inclusion of ozone chemical feedback would amplify the response, potentially reproducing observed amplitudes as well as variability.

The GISS model results and empirical reconstructions both suggest that solar-forced regional climate changes during the Maunder Minimum appeared predominantly as a shift toward the low AO/NAO index. Although global average temperature changes were small, modeled regional cooling over the continents during winter was up to five times greater. Changes in ocean circulation were not considered in this model. However, given the sensitivity of the North Atlantic to AO/NAO forcing (37), oceanic changes may well have been triggered as a response to the atmospheric changes (38). Such oceanic changes would themselves further modify the pattern of SST in the North Atlantic (39) and, to a lesser extent, the downstream air temperature anomalies in Europe.

These results provide evidence that relatively small solar forcing may play a significant role in century-scale NH winter climate change. This suggests that colder winter temperatures over the NH continents during portions of the 15th through the 17th centuries (sometimes called the Little Ice Age) and warmer temperatures during the 12th through 14th centuries (the putative Medieval Warm Period) may have been influenced by long-term solar variations.


View Abstract

Navigate This Article