A first-order global model of late Cenozoic climatic change II. Further analysis based on a simplification of CO2 dynamics

A model developed recently for the long-term variations of global ice mass, carbon dioxide, and mean ocean temperature through the late Cenozoic is simplified by hypothesizing a new equation for the CO₂ variations containing one less adjustable parameter, but retaining the essential physical content of the previous equation (including nonlinearity and the potential for instability). By assuming plausible time constants for the glacial ice mass and global mean ocean temperature, and setting the values of six adjustable parameters (rate constants), a solution for the last 5 My is obtained displaying many of the features observed over this period, including the transition to the near-100 ky major ice-age oscillations of the late Pleistocene. In obtaining this solution it is also assumed that variations in tectonic forcing lead to a reduction of the equilibrium CO₂ concentration (perhaps due to increased weathering of rapidly uplifted mountain ranges over this period). As a consequence of this CO₂ reduction, the model dynamical system can bifurcate to a free oscillatory ice-age regime that is under the pacemaker influence of earthorbital (Milankovitch) forcing. Expanded discussions are given of the surface temperature variations accompanying the evolution of ice, CO₂, and ocean temperature, and of the bifurcation properties of the model from both mathematical and physical viewpoints.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dumenil     

The initiation of ice sheet growth,Milankovitch solar radiation variations,and the 100 ky ice age cycle

Much work has gone into deciphering the causes of the large scale glacial/interglacial variations in the climate system over the last 900 000 years. While variations on the 41 thousand year (ky) and 23 ky time scales seem to be linearly linked to the variations in the distribution of solar radiation at the top of the atmosphere, Milankovitch solar radiation variations, the causes of the dominant 100 ky cycle in the geologic record are still unknown. One of the aspects of this cycle that is not well understood is how large scale ice sheet growth is initiated. Here we describe the mechanisms by which large scale ice sheet growth may have been initiated by the changes in the seasonal and latitudinal distribution of solar radiation over the past 160 ky. This is done through the use of a coupled energy balance climate-thermodynamic sea ice model that includes a hydrologic cycle which computes precipitation, and a land surface energy balance which determines the net accumulation of snow and ice. Results indicate that the initiation of ice sheet growth is possible during times of extremely low summer solstice solar radiation as a result of a large decrease in ablation during the critical melt season.     

Impacts of shortwave radiation forcing on ENSO: a study with a coupled tropical ocean-atmosphere model

 We describe a coupled tropical ocean-atmosphere model that represents a new class of models that fill the gap between anomaly coupled models and fully coupled general circulation models. Both the atmosphere and ocean are described by two and half layer primitive equation models, which emphasize the physical processes in the oceanic mixed layer and atmospheric boundary layer. Ocean and atmosphere are coupled through both momentum and heat flux exchanges without explicit flux correction. The coupled model, driven by solar radiation, reproduces a realistic annual cycle and El Nino-Southern Oscillation (ENSO). In the presence of annual mean shortwave radiation forcing, the model exhibits an intrinsic mode of ENSO. The oscillation period depends on the mean forcing that determines the coupled mean state. A perpetual April (October) mean forcing prolongs (shortens) the oscillation period through weakening (enhancing) the mean upwelling and mean vertical temperature gradients. The annual cycle of the solar forcing is shown to have fundamental impacts on the behavior of ENSO cycles through establishing a coupled annual cycle that interacts with the ENSO mode. Due to the annual cycle solar forcing, the single spectral peak of the intrinsic ENSO mode becomes a double peak with a quasi-biennial and a low-frequency (4–5 years) component; the evolution of ENSO becomes phase-locked to the annual cycle; and the amplitude and frequency of ENSO become variable on an interdecadal time scale due to interactions of the mean state and the two ENSO components. The western Pacific monsoon (the annual shortwave radiation forcing in the western Pacific) is primarily responsible for the generation of the two ENSO components. The annual march of the eastern Pacific ITCZ tends to lock ENSO phases to the annual cycle. The model's deficiencies, limitations, and future work are also discussed. Received: 15 June 1999 / Accepted: 11 December 1999     

Three basic problems of paleoclimatic modeling: a personal perspective and review

The development of a theory of the evolution of the climate of the earth over millions of years can be subdivided into three fundamental, nested, problems: Firstly, to establish by equilibrium climate models (e.g., general circulation models) the diagnostic relations, valid at any time, between the fast-response climate variables (i.e., the weather statistics) and both the prescribed external radiative forcing and the prescribed distribution of the slow-response variables (e.g., the ice sheets and shelves, the deep ocean state, and the atmospheric CO₂ concentration). Secondly, to construct, by an essentially inductive process, a model of the time-dependent evolution of the slow-response climatic variables over time scales longer than the damping times of these variables but shorter than the time scale of ultra-slow tectonic and astronomical changes in the boundary conditions (e.g., altered geography and elevation of the continents, slow outgassing and weathering and and solar radiative output). Thirdly, to determine the nature of these ultra-slow processes and their effects on the evolution of the equilibrium state of the climatic system about which the previously mentioned time-dependent variations occur. In this review we discuss the basis for this resolution, and give a broad overview of the contributions that have been made thus far in each area, emphasizing the work of the Yale climate group.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil     

The UVic earth system climate model: Model description,climatology, and applications to past,present and future climates

Abstract

A new earth system climate model of intermediate complexity has been developed and its climatology compared to observations. The UVic Earth System Climate Model consists of a three‐dimensional ocean general circulation model coupled to a thermodynamic/dynamic sea‐ice model, an energy‐moisture balance atmospheric model with dynamical feedbacks, and a thermomechanical land‐ice model. In order to keep the model computationally efficient a reduced complexity atmosphere model is used. Atmospheric heat and freshwater transports are parametrized through Fickian diffusion, and precipitation is assumed to occur when the relative humidity is greater than 85%. Moisture transport can also be accomplished through advection if desired. Precipitation over land is assumed to return instantaneously to the ocean via one of 33 observed river drainage basins. Ice and snow albedo feedbacks are included in the coupled model by locally increasing the prescribed latitudinal profile of the planetary albedo. The atmospheric model includes a parametrization of water vapour/planetary longwave feedbacks, although the radiative forcing associated with changes in atmospheric CO₂ is prescribed as a modification of the planetary longwave radiative flux. A specified lapse rate is used to reduce the surface temperature over land where there is topography. The model uses prescribed present‐day winds in its climatology, although a dynamical wind feedback is included which exploits a latitudinally‐varying empirical relationship between atmospheric surface temperature and density. The ocean component of the coupled model is based on the Geophysical Fluid Dynamics Laboratory (GFDL) Modular Ocean Model 2.2, with a global resolution of 3.6° (zonal) by 1.8° (meridional) and 19 vertical levels, and includes an option for brine‐rejection parametrization. The sea‐ice component incorporates an elastic‐viscous‐plastic rheology to represent sea‐ice dynamics and various options for the representation of sea‐ice thermodynamics and thickness distribution. The systematic comparison of the coupled model with observations reveals good agreement, especially when moisture transport is accomplished through advection.

Global warming simulations conducted using the model to explore the role of moisture advection reveal a climate sensitivity of 3.0°C for a doubling of CO₂, in line with other more comprehensive coupled models. Moisture advection, together with the wind feedback, leads to a transient simulation in which the meridional overturning in the North Atlantic initially weakens, but is eventually re‐established to its initial strength once the radiative forcing is held fixed, as found in many coupled atmosphere General Circulation Models (GCMs). This is in contrast to experiments in which moisture transport is accomplished through diffusion whereby the overturning is reestablished to a strength that is greater than its initial condition.

When applied to the climate of the Last Glacial Maximum (LGM), the model obtains tropical cooling (30°N‐30°S), relative to the present, of about 2.1°C over the ocean and 3.6°C over the land. These are generally cooler than CLIMAP estimates, but not as cool as some other reconstructions. This moderate cooling is consistent with alkenone reconstructions and a low to medium climate sensitivity to perturbations in radiative forcing. An amplification of the cooling occurs in the North Atlantic due to the weakening of North Atlantic Deep Water formation. Concurrent with this weakening is a shallowing of, and a more northward penetration of, Antarctic Bottom Water.

Climate models are usually evaluated by spinning them up under perpetual present‐day forcing and comparing the model results with present‐day observations. Implicit in this approach is the assumption that the present‐day observations are in equilibrium with the present‐day radiative forcing. The comparison of a long transient integration (starting at 6 KBP), forced by changing radiative forcing (solar, CO₂, orbital), with an equilibrium integration reveals substantial differences. Relative to the climatology from the present‐day equilibrium integration, the global mean surface air and sea surface temperatures (SSTs) are 0.74°C and 0.55°C colder, respectively. Deep ocean temperatures are substantially cooler and southern hemisphere sea‐ice cover is 22% greater, although the North Atlantic conveyor remains remarkably stable in all cases. The differences are due to the long timescale memory of the deep ocean to climatic conditions which prevailed throughout the late Holocene. It is also demonstrated that a global warming simulation that starts from an equilibrium present‐day climate (cold start) underestimates the global temperature increase at 2100 by 13% when compared to a transient simulation, under historical solar, CO₂ and orbital forcing, that is also extended out to 2100. This is larger (13% compared to 9.8%) than the difference from an analogous transient experiment which does not include historical changes in solar forcing. These results suggest that those groups that do not account for solar forcing changes over the twentieth century may slightly underestimate (～3% in our model) the projected warming by the year 2100.     

Phase Modulation Effect of the Rubincam Insolation Variations on Climate Change

Summary For astronomical seasons, Rubincam insolation deviations at latitude 65° N varied from 218.50 Wm⁻² to 225.75 Wm⁻² (3%). The periodicity of the insolation cycles varied from 36.7 Kyr to 44.7 Kyr (20%) due to phase shift. Phase shift of insolation variations can induce asymmetry of the insolation cycles, permitting rapid melting and prolonged glaciation of ice sheets to occur. For instance, an abnormal decrease of the insolation frequency during the longer period of glacial interval would prolong glaciation into deep ice age. In this study, we apply Rubincam’s insolation equations to investigate the phase shift effect of insolation variations on climate change. Using complex transforms of the changing insolation, we have detected a phase modulation signal in the insolation variations. As a result, an especially new and interesting series of the phase-related insolation pulsation is established. The phase modulated insolation is then introduced as a forcing function into energy balance climate models. Results of model computations shed new insights into the spectrum of the paleoclimatic proxy-data. It is shown that phase modulation of the insolation may provide an appropriate and complete external forcing mechanism to which the climate system would respond. The 100 Kyr cycle of the frequency modulation of the Rubincam’s insolation variations does seem adequate to change the climate. Received July 16, 1997 Revised May 18, 1998     

Past millennial solar forcing magnitude

A set of global climate model simulations for the last thousand years developed by the Max Planck Institute is compared with paleoclimate proxy data and instrumental data, focusing on surface temperatures for land areas between 30° and 75°N. The proxy data are obtained from six previously published Northern Hemispheric-scale temperature reconstructions, here re-calibrated for consistency, which are compared with the simulations utilizing a newly developed statistical framework for ranking several competing simulations by means of their statistical distance against past climate variations. The climate model simulations are driven by either “low” or “high” solar forcing amplitudes (0.1 and 0.25 % smaller total solar irradiance in the Maunder Minimum period compared to the present) in addition to several other known climate forcings of importance. Our results indicate that the high solar forcing amplitude results in a poorer match with the hemispheric-scale temperature reconstructions and lends stronger statistical support for the low-amplitude solar forcing. However, results are likely conditional upon the sensitivity of the climate model used and strongly dependent on the choice of temperature reconstruction, hence a greater consensus is needed regarding the reconstruction of past temperatures as this currently provides a great source of uncertainty.     

The importance of precessional signals in the tropical climate

Past research on the climate response to orbital forcing has emphasized the glacial-interglacial variations in global ice volume, global-mean temperature, and the global hydrologic cycle. This emphasis may be inappropriate in the tropics, where the response to precessional forcing is likely to be somewhat independent of the glacial-interglacial variations, particularly in variables relating to the hydrologic cycle. To illustrate this point, we use an atmospheric general circulation model coupled to a slab ocean model, performing experiments that quantify the tropical climates response to (1) opposite phases of precessional forcing, and (2) Last Glacial Maximum boundary conditions. While the glacially-forced tropical temperature changes are typically more than an order of magnitude larger than those arising from precessional forcing, the hydrologic signals stemming from the two forcings are comparable in magnitude. The mechanisms behind these signals are investigated and shown to be quite distinct for the precessional and glacial forcing. Because of strong dynamical linkages in the tropics, the model results illustrate the impossibility of predicting the local hydrologic response to external forcing without understanding the response at much larger spatial scales. Examples from the paleoclimate record are presented as additional evidence for the importance of precessional signals in past variations of the tropical climate.     

Late Pliocene to Pleistocene sensitivity of the Greenland Ice Sheet in response to external forcing and internal feedbacks

The timing and nature of ice sheet variations on Greenland over the last ～5 million years remain largely uncertain. Here, we use a coupled climate-vegetation-ice sheet model to determine the climatic sensitivity of Greenland to combined sets of external forcings and internal feedbacks operating on glacial-interglacial timescales. In particular, we assess the role of atmospheric pCO₂, orbital forcing, and vegetation dynamics in modifying thresholds for the onset of glaciation in late Pliocene and Pleistocene. The response of circum-Arctic vegetation to declining levels of pCO₂ (from 400 to 200 ppmv) and decreasing summer insolation includes a shift from boreal forest to tundra biomes, with implications for the surface energy balance. The expansion of tundra amplifies summer surface cooling and heat loss from the ground, leading to an expanded summer snow cover over Greenland. Atmospheric and land surface fields respond to forcing most prominently in late spring-summer and are more sensitive at lower Pleistocene-like levels of pCO₂. We find cold boreal summer orbits produce favorable conditions for ice sheet growth, however simulated ice sheet extents are highly dependent on both background pCO₂ levels and land-surface characteristics. As a result, late Pliocene ice sheet configurations on Greenland differ considerably from late Pleistocene, with smaller ice caps on high elevations of southern and eastern Greenland, even when orbital forcing is favorable for ice sheet growth.     

Radiative Forcing and Climate Response Due to Black Carbon in Snow and Ice

The radiative forcing and climate response due to black carbon(BC) in snow and/or ice were investigated by integrating observed effects of BC on snow/ice albedo into an atmospheric general circulation model(BCC AGCM2.0.1) developed by the National Climate Center(NCC) of the China Meteorological Administration(CMA).The results show that the global annual mean surface radiative forcing due to BC in snow/ice is +0.042 W m 2,with maximum forcing found over the Tibetan Plateau and regional mean forcing exceeding +2.8 W m 2.The global annual mean surface temperature increased 0.071 C due to BC in snow/ice.Positive surface radiative forcing was clearly shown in winter and spring and increased the surface temperature of snow/ice in the Northern Hemisphere.The surface temperatures of snow-covered areas of Eurasia and North America in winter(spring) increased by 0.83 C(0.6 C) and 0.83 C(0.46 C),respectively.Snowmelt rates also increased greatly,leading to earlier snowmelt and peak runoff times.With the rise of surface temperatures in the Arctic,more water vapor could be released into the atmosphere,allowing easier cloud formation,which could lead to higher thermal emittance in the Arctic.However,the total cloud forcing could decrease due to increasing cloud cover,which will offset some of the positive feedback mechanism of the clouds.     

A nonlinear impulse response model of the coupled carbon cycle-climate system (NICCS)

 Impulse-response-function (IRF) models are designed for applications requiring a large number of climate change simulations, such as multi-scenario climate impact studies or cost-benefit integrated-assessment studies. The models apply linear response theory to reproduce the characteristics of the climate response to external forcing computed with sophisticated state-of-the-art climate models like general circulation models of the physical ocean-atmosphere system and three-dimensional oceanic-plus-terrestrial carbon cycle models. Although highly computer efficient, IRF models are nonetheless capable of reproducing the full set of climate-change information generated by the complex models against which they are calibrated. While limited in principle to the linear response regime (less than about 3 ^∘C global-mean temperature change), the applicability of the IRF model presented has been extended into the nonlinear domain through explicit treatment of the climate system's dominant nonlinearities: CO₂ chemistry in ocean water, CO₂ fertilization of land biota, and sublinear radiative forcing. The resultant nonlinear impulse-response model of the coupled carbon cycle-climate system (NICCS) computes the temporal evolution of spatial patterns of climate change for four climate variables of particular relevance for climate impact studies: near-surface temperature, cloud cover, precipitation, and sea level. The space-time response characteristics of the model are derived from an EOF analysis of a transient 850-year greenhouse warming simulation with the Hamburg atmosphere-ocean general circulation model ECHAM3-LSG and a similar response experiment with the Hamburg carbon cycle model HAMOCC. The model is applied to two long-term CO₂ emission scenarios, demonstrating that the use of all currently estimated fossil fuel resources would carry the Earth's climate far beyond the range of climate change for which reliable quantitative predictions are possible today, and that even a freezing of emissions to present-day levels would cause a major global warming in the long term. Received: 28 January 2000 / Accepted: 9 March 2001     

Global warming and the growth of ice sheets

Recent research has suggested that warmer conditions, that may result from increased levels of CO₂ in the atmosphere, may induce the growth of the Northern Hemisphere ice sheets (Miller and de Vernal 1992) through the impact of warmer temperature on the water carrying capacity of the atmosphere and thus on precipitation. In this study we examine this possibility by using a coupled energy balance climate-thermo-dynamic sea ice model. Results indicate that if summer ice albedo is high enough, and there is some mechanism for initially maintaining ice through the summer season, then it may be possible to have ice sheet growth under the conditions of CO₂ induced warming.This paper was presented at the Second International Conference on Modelling of Global Climate Variability, held in Hamburg 7–11 September 1992 under the auspices of the Max Planck Institute for Meteorology. Guest Editor for these papers is L. Dümenil     

Design and performance of a new cloud microphysics scheme developed for the ECHAM general circulation model

A new cloud microphysics scheme including a prognostic treatment of cloud ice (PCI) is developed to yield a more physically based representation of the components of the atmospheric moisture budget in the general circulation model ECHAM. The new approach considers cloud water and cloud ice as separate prognostic variables. The precipitation formation scheme for warm clouds distinguishes between maritime and continental clouds by considering the cloud droplet number concentration, in addition to the liquid water content. Based on several observational data sets, the cloud droplet number concentration is derived from the sulfate aerosol mass concentration as given from the sulfur cycle simulated by ECHAM. Results obtained with the new scheme are compared to satellite observations and in situ measurements of cloud physical and radiative properties. In general, the standard model ECHAM4 and also PCI capture the overall features, and the simulated results usually lie within the range of observed uncertainty. As compared to ECHAM4, only slight improvements are achieved with the new scheme. For example, the overestimated liquid water path and total cloud cover over convectively active regions are reduced in PCI. On the other hand, some shortcomings of the standard model such as underestimated shortwave cloud forcing over the extratropical oceans of the respective summer hemisphere are more pronounced in PCI.This paper was presented at the Third International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 4–9 Sept. 1995 under the auspices of the Max Planck Institute for Meteorology, Hamburg. Editor for these papers is L. Dümenil.     

Natural forcing of climate during the last millennium: fingerprint of solar variability

The variability of the climate during the last millennium is partly forced by changes in total solar irradiance (TSI). Nevertheless, the amplitude of these TSI changes is very small so that recent reconstruction data suggest that low frequency variations in the North Atlantic Oscillation (NAO) and in the thermohaline circulation may have amplified, in the North Atlantic sector and mostly in winter, the radiative changes due to TSI variations. In this study we use a state-of-the-art climate model to simulate the last millennium. We find that modelled variations of surface temperature in the Northern Hemisphere are coherent with existing reconstructions. Moreover, in the model, the low frequency variability of this mean hemispheric temperature is found to be correlated at 0.74 with the solar forcing for the period 1001?C1860. Then, we focus on the regional climatic fingerprint of solar forcing in winter and find a significant relationship between the low frequency TSI forcing and the NAO with a time lag of more than 40?years for the response of the NAO. Such a lag is larger than the around 20-year lag suggested in other studies. We argue that this lag is due, in the model, to a northward shift of the tropical atmospheric convection in the Pacific Ocean, which is maximum more than four decades after the solar forcing increase. This shift then forces a positive NAO through an atmospheric wave connection related to the jet-stream wave guide. The shift of the tropical convection is due to the persistence of anomalous warm SST forcing the anomalous precipitation, associated with the advection of warm SST by the North Pacific subtropical gyre in a few decades. Finally, we analyse the response of the Atlantic meridional overturning circulation to solar forcing and find that the former is weakened when the latter increases. Changes in wind stress, notably due to the NAO, modify the barotropic streamfunction in the Atlantic 50?years after solar variations. This implies a wind-driven modification of the oceanic circulation in the Atlantic sector in response to changes in solar forcing, in addition to the variations of the thermohaline circulation.     

Potential Role of Solar Variability as an Agent for Climate Change

Numerical experiments have been carried out with a two-dimensional sector averaged global climate model in order to assess the potential impact of solar variability on the Earth's surface temperature from 1700 to 1992. This was done by investigating the model response to the variations in solar radiation caused by the changes in the Earth's orbital elements, as well as by the changes intrinsic to the Sun. In the absence of a full physical theory able to explain the origin of the observed total solar irradiance variations, three different total solar irradiance reconstructions have been used. A total solar irradiance change due to the photospheric effects incorporated in the Willson and Hudson (1988) parameterization, and the newly reconstructed solar total irradiance variations from the solar models of Hoyt and Schatten (1993) and Lean et al. (1995). Our results indicate that while the influence of the orbital forcing on the annual and global mean surface temperature is negligible at the century time scale, the monthly mean response to this forcing can be quite different from one month to another. The modelled global warming due to the three investigated total solar irradiance reconstructions is insufficient to reproduce the observed 20th century warming. Nevertheless, our simulated surface temperature response to the changes in the Sun's radiant energy output suggests that the Gleissberg cycle (88 years) solar forcing should not be neglected in explaining the century-scale climate variations. Finally, spectral analysis seems to point out that the 10- to 12-year oscillations found in the recorded Northern Hemisphere temperature variations from 1700 to 1992 could be unrelated to the solar forcing. Such a result could indicate that the eleven-year period which is frequently found in climate data might be related to oscillations in the atmosphere or oceans, internal to the climate system.     

Long-term patterns of solar irradiance forcing in model experiments and proxy based surface temperature reconstructions

Comparisons are made of long-term empirical and model-estimated patterns of solar irradiance forcing during a 200-year period (1650-1850), which precedes any apparent anthropogenic influence on climate. This interval encompasses a considerable range (approximately 4 W/m²) of estimated variation in solar output, including the "Maunder" and "Dalton" Minima of solar irradiance, and an intervening interval of relatively high values of irradiance, but does not encroach into the industrial era wherein it is difficult to separate solar and anthropogenic influences. Particular emphasis is placed on comparing empirical and modeled patterns of forced surface temperature variation. The empirical patterns bear a greater similarity to the pattern of forced response of a coupled ocean-atmosphere general circulation model (AOGCM) than with an independent model simulation result using an ocean with specified heat transport, both in terms of the spatial pattern of response and implied global mean sensitivity to forcing. Heightened sensitivity in the western Pacific warm pool apparent in the empirical response pattern, is not observed in the forced response of the coupled model. It is possible that this pattern is the result of feedback processes not currently reproduced in course-resolution coupled models. The greatest empirical response is found at the multidecadal-to-century (> 40 year period) time scale, for which the forcing is dominated by the roughly 90-year Gleissberg Cycle of irradiance. This indicates a global-mean sensitivity (approximately 0.3 K/W/m²), which is close to the coupled model result (approximately 0.4 K/W/m²). At decadal time scales (8-25 year period), for which the forcing is dominated by the 11-year and 22-year period solar cycles), the temperature sensitivity is moderately reduced, and its spatial pattern of response is dominated by an apparent resonance with known decadal modes of climate variability.     

Solar irradiance during the last 1200 years based on cosmogenic nuclides

Based on a quantitative study of the common fluctuations of ¹⁴C and ¹⁰Be production rates, we have derived a time series of the solar magnetic variability over the last 1200 years. This record is converted into irradiance variations by linear scaling based on previous studies of sun‐like stars and of the sun's behavior over the last few centuries. The new solar irradiance record exhibits low values during the well‐known solar minima centered at about 1900, 1810 (Dalton) and 1690 ad (Maunder). Further back in time, a rather long period between 1450 and 1750 ad is characterized by low irradiance values. A shorter period is centered at about 1200 ad , with irradiance slightly higher or similar to present day values. It is tempting to correlate these periods with the so‐called "little ice age" and "medieval warm period", respectively. An accurate quantification of the climatic impact of this new irradiance record requires the use of coupled atmosphere–ocean general circulation models (GCMs). Nevertheless, our record is already compatible with a global cooling of about 0.5‐1°C during the "little ice age", and with a general cooling trend during the past millenium followed by global warming during the 20th century (Mann et al., 1999).     

东亚地区小冰期气候的模拟

本文使用ECHO-G全球气候模式对1550～1850年的小冰期气候进行了300个模式年的模拟,着重分析了东亚地区小冰期的温度变化特征,并与目前所得到的小冰期气候重建结果进行了对比。结果表明,在考虑了太阳辐射、火山活动、CO2和CH4等主要气候影响驱动因子的条件下,较好地模拟出了东亚地区的小冰期气候特征,并与其它手段的气候重建结果相吻合,显示太阳活动和火山活动是小冰期气候形成的主要原因。     

Central antarctic climate response to the solar cycle

Antarctic “Vostok” station works most closely to the center of the ice cap among permanent year-around stations. Climate conditions are exclusively stable: low precipitation level, cloudiness and wind velocity. These conditions can be considered as an ideal model laboratory to study the surface temperature response on solar irradiance variability during 11-year cycle of solar activity. Here we solve an inverse heat conductivity problem: calculate the boundary heat flux density (HFD) from known evolution of temperature. Using meteorological temperature record during (1958–2011) we calculated the HFD variation about 0.2–0.3 W/m² in phase with solar activity cycle. This HFD variation is derived from 0.5 to 1 °C temperature variation and shows relatively high climate sensitivity per 0.1 % of solar radiation change. This effect can be due to the polar amplification phenomenon, which predicts a similar response 0.3–0.8 °C/0.1 % (Gal-Chen and Schneider in Tellus 28:108–121, 1975). The solar forcing (TSI) is disturbed by volcanic forcing (VF), so that their linear combination TSI + 0.5VF empirically provides higher correlation with HFD (r = 0.63 ± 0.22) than TSI (r = 0.50 ± 0.24) and VF (r = 0.41 ± 0.25) separately. TSI shows higher wavelet coherence and phase agreement with HFD than VF.     

An Antarctic view of Beryllium-10 and solar activity for the past millennium

Beryllium-10 in ice provides a valuable proxy of solar activity. However, complex production pathways, atmospheric transport, and deposition processes impede its quantitative interpretation. Here, we examine the influence of deposition processes on two Be-10 ice core records from Central Antarctica (South Pole and Dome Fuji stations), covering the last millennium. We try to quantify how Be-10 variations in ice relate to variations in Be-10 production, and the bias associated to this relationship. An independent bias estimation is provided by comparing atmospheric radiocarbon variations reconstructed from tree rings and deduced from Be-10 variations. Both techniques suggest an uncertainty of the order of 10% in Be-10 production. This uncertainty estimate does not account for the geographical origin of Be-10, which remains a major issue. Because both Be-10 records are so similar, we propose to average them as a means to decrease the unshared (non solar) variability. This average record provides a new reconstruction of solar modulation parameter ?? and total solar irradiance over the last ~1,300?years. The lowest solar activity is found during the so-called Sp?rer Minimum (around AD 1450). The highest activities are found during the 8th century and over the last decades: as shown in previous studies, our results suggest that the recent solar activity is not exceptionally high for the last millennium.