Multi-century simulations with the coupled GISS–HYCOM climate model: control experiments

Multi-century climate simulations obtained with the GISS atmospheric general circulation model coupled to the hybrid-isopycnic ocean model HYCOM are described. Greenhouse gas concentrations are held fixed in these experiments to investigate the coupled model’s ability to reproduce the major features of today’s climate with minimal drift. Emphasis is placed on the realism of the oceanic general circulation and its effect on air–sea exchange processes. Several model runs using different closures for turbulent vertical exchange as well as improvements to reduce vertical numerical diffusion are compared with climate observations. As in previous studies, the Southern Ocean emerges as the Achilles Heel of the ocean model; deficiencies in its physical representation had far-reaching consequences in early experiments with the coupled model and have provided the strongest impetus for model improvement. The overarching goal of this work is to add diversity to the pool of ocean models available for climate prediction and thereby reduce biases that may stand in the way of assessing climate prediction uncertainty.

---

Decadal interactions between the western tropical Pacific and the North Atlantic Oscillation

The relationship between interdecadal variations of tropical sea surface temperature (SST) in the last 120 years and circulation anomalies related to the North Atlantic Oscillation (NAO) is investigated in this study. Using an atmospheric general circulation model (AGCM), we confirm observational evidence that variations in the SST gradient in the western tropical Pacific are related to the NAO anomalies on decadal timescale, and may be contributing to the shift towards the positive NAO phase observed in the late 20th century. The role played by the Indian Ocean-NAO teleconnection, advocated in recent studies focused on the last 50 years, is also assessed in the context of the 120-year long record. It is suggested that a positive feedback between the Pacific SST and the hemispheric circulation pattern embedding the decadal NAO signal may act to enhance the internal variability of the coupled ocean–atmosphere system, and justify the stronger teleconnection found in observational data than in SST-forced AGCM experiments.

---

Intraseasonal variability in the equatorial Atlantic-West Africa during March–June

Intraseasonal (30–80 days) variability in the equatorial Atlantic-West African sector during March–June is investigated using various recently-archived satellite measurements and the NCEP/DOE AMIP-II reanalysis daily data. The global connections of regional intraseasonal signals are first examined for the period of 1979–2006 through lag-regression analyses of convection (OLR) and other dynamic components against a regional intraseasonal convective (OLR) index. The eastward-propagating features of convection can readily be seen, accompanied by coherent circulation anomalies, similar to those for the global tropical intraseasonal mode, i.e., the Madden–Julian oscillation (MJO). The regressed TRMM rainfall (3B42) anomalies during the TRMM period (1998–2006) manifest similar propagating features as for the regressed OLR anomalies during 1979–2006. These coherent features hence tend to suggest that the regional intraseasonal convective signals might be mostly a regional response to, or closely associated with the MJO, and probably contribute to the MJO’s global propagation. Atmospheric and surface intraseasonal variability during March–June of 1998–2006 are further examined using the high-quality TRMM Microwave Imager (TMI) sea surface temperature (SST), columnar water vapor, and cloud liquid water, and the QuikSCAT oceanic winds (2000–2006). Enhanced (suppressed) convection or positive (negative) rainfall anomalies approximately cover the entire basin (0°–10°N, 30°W–10°E) during the passage of intraseasonal convective signals, accompanied by anomalous surface westerly (easterly) flow. Furthermore, a unique propagating feature seems to exist within the tropical Atlantic basin. Rainfall anomalies always appear first in the northwestern basin right off the coast of South America, and gradually extend eastward to cover the entire basin. A dipolar structure of rainfall anomalies with cross-equatorial surface wind anomalies can thus be observed during this evolution, similar to the anomaly patterns on the interannual time scale discovered in past studies. Coherent intraseasonal variations and patterns can also be found in other physical components.

Sensitivity of global warming to the pattern of tropical ocean warming

The current generations of climate models are in substantial disagreement as to the projected patterns of sea surface temperatures (SSTs) in the Tropics over the next several decades. We show that the spatial patterns of tropical ocean temperature trends have a strong influence on global mean temperature and precipitation and on global mean radiative forcing. We identify the SST patterns with the greatest influence on the global mean climate and find very different, and often opposing, sensitivities to SST changes in the tropical Indian and West Pacific Oceans. Our work stresses the need to reduce climate model biases in these sensitive regions, as they not only affect the regional climates of the nearby densely populated continents, but also have a disproportionately large effect on the global climate.

---

Simulations of the Madden–Julian oscillation in four pairs of coupled and uncoupled global models

The status of the numerical reproduction of the Madden–Julian Oscillation (MJO) by current global models was assessed through diagnoses of four pairs of coupled and uncoupled simulations. Slow eastward propagation of the MJO, especially in low-level zonal wind, is realistic in all these simulations. However, the simulated MJO suffers from several common problems. The MJO signal in precipitation is generally too weak and often eroded by an unrealistic split of an equatorial maximum of precipitation into a double ITCZ structure over the western Pacific. The MJO signal in low-level zonal wind, on the other hand, is sometimes too strong over the eastern Pacific but too weak over the Indian Ocean. The observed phase relationship between precipitation and low-level zonal wind associated with the MJO in the western Pacific and their coherence in general are not reproduced by the models. The seasonal migration in latitude of MJO activity is missing in most simulations. Air–sea coupling generally strengthens the simulated eastward propagating signal, but its effects on the phase relationship and coherence between precipitation and low-level zonal wind, and on their geographic distributions, seasonal cycles, and interannual variability are inconsistent among the simulations. Such inconsistency cautions generalization of results from MJO simulations using a single model. In comparison to observations, biases in the simulated MJO appear to be related to biases in the background state of mean precipitation, low-level zonal wind, and boundary-layer moisture convergence. This study concludes that, while the realistic simulations of the eastward propagation of the MJO are encouraging, reproducing other fundamental features of the MJO by current global models remains an unmet challenge.

---

Remotely forced variability in the tropical Atlantic Ocean

An ensemble of eight hindcasts has been conducted using an ocean-atmosphere general circulation model fully coupled only within the Atlantic basin, with prescribed observational sea surface temperature (SST) for 1950–1998 in the global ocean outside the Atlantic basin. The purpose of these experiments is to understand the influence of the external SST anomalies on the interannual variability in the tropical Atlantic Ocean. Statistical methods, including empirical orthogonal function analysis with maximized signal-to-noise ratio, have been used to extract the remotely forced Atlantic signals from the ensemble of simulations. It is found that the leading external source on the interannual time scales is the El Niño/Southern Oscillation (ENSO) in the Pacific Ocean. The ENSO signal in the tropical Atlantic shows a distinct progression from season to season. During the boreal winter of a maturing El Niño event, the model shows a major warm center in the southern subtropical Atlantic together with warm anomalies in the northern subtropical Atlantic. The southern subtropical SST anomalies is caused by a weakening of the southeast trade winds, which are partly associated with the influence of an atmospheric wave train generated in the western Pacific Ocean and propagating into the Atlantic basin in the Southern Hemisphere during boreal fall. In the boreal spring, the northern tropical Atlantic Ocean is warmed up by a weakening of the northeast trade winds, which is also associated with a wave train generated in the central tropical Pacific during the winter season of an El Niño event. Apart from the atmospheric planetary waves, these SST anomalies are also related to the sea level pressure (SLP) increase in the eastern tropical Atlantic due to the global adjustment to the maturing El Niño in the tropical Pacific. The tropical SLP anomalies are further enhanced in boreal spring, which induce anomalous easterlies on and to the south of the equator and lead to a dynamical oceanic response that causes cold SST anomalies in the eastern and equatorial Atlantic from boreal spring to summer. Most of these SST anomalies persist into the boreal fall season.

---

Spectrothermography of carbonaceous particles

Spectrothermography, defined as the evaluation of thermograms of carbon evolved from heated aerosol samples, is a technique for evaluating differences in particle characteristics as they relate to emission sources and processes that modify particle evolution. Here we describe the inherent uncertainties and demonstrate the utility of this technique with an evaluation of samples that were collected with eight stage cascade impactors at three sites in Mexico City over a period of 5 months. The study was implemented with statistical analysis based on tests for goodness of fit to separate thermograms with distinctive shapes related to the relative amounts of organic and elemental carbon mass that evolves as a function of temperature. Thermograms with unique shapes were found for particles with aerodynamic diameters of 1–10, 0.56–1, 0.32–0.56 and 0.18–0.32 μm with further differentiation of curves related to the relative amounts of gasoline and diesel fuel that was combusted in the region of the three sites. The common shapes fit 32–42% of samples in each particle size range and indicate that this type of analysis can help distinguish differences in the primary sources of organic and elemental carbon.

The influence of sea surface temperature anomalies on low-frequency variability of the North Atlantic Oscillation

The influence of sea surface temperature anomalies (SSTA) on multi-year persistence of the North Atlantic Oscillation (NAO) during the second half of the twentieth century is investigated using the Center for Ocean-Land-Atmosphere Studies (COLA) Atmospheric GCM (AGCM) with an emphasis on isolating the geographic location of the SSTA that produce this influence. The present study focuses on calculating the atmospheric response to the SSTA averaged over 1988–1995 (1961–1968) corresponding to the observed period of strong persistence of the positive (negative) phase of the decadal NAO. The model response to the global 1988–1995 average SSTA shows a statistically significant large-scale pattern characteristic of the positive phase of the NAO. Forcing with the global 1961–1968 average SSTA generates a NAO of the opposite polarity compared to observations. However, all large-scale features both in the model and observations during this period are weaker in magnitude and less significant compared to 1988–1995. Additional idealized experiments show that over the northern center of the NAO the non-linear component of the forced response appears to be quite important and acts to enhance the positive NAO signal. On the other hand, over the southern center where the model response is the strongest, it is also essentially linear. The 1988–1995 average SSTA restricted to the western tropical Pacific region produce a positive NAO remarkably similar in structure but stronger in magnitude than the model response to the global and tropical Indo-Pacific 1988–1995 forcing. A 200-hPa geopotential height response in these experiments shows a positive anomaly over the southern center of the NAO embedded in the Rossby wave trains propagating from the western tropical Pacific. Indian Ocean SSTA lead to much weaker positive NAO primarily through the effect on its northern center. SST forcing confined to the North Atlantic north of equator does not produce a response statistically different from the control simulation, suggesting that it is not strong enough to significantly affect the phase of the decadal NAO. Inclusion of the South Atlantic north of 45° south does not change this result.

---

Analysis of the projected regional sea-ice changes in the Southern Ocean during the twenty-first century

Using the set of simulations performed with atmosphere-ocean general circulation models (AOGCMs) for the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR4), the projected regional distribution of sea ice for the twenty-first century has been investigated. Averaged over all those model simulations, the current climate is reasonably well reproduced. However, this averaging procedure hides the errors from individual models. Over the twentieth century, the multimodel average simulates a larger sea-ice concentration decrease around the Antarctic Peninsula compared to other regions, which is in qualitative agreement with observations. This is likely related to the positive trend in the Southern Annular Mode (SAM) index over the twentieth century, in both observations and in the multimodel average. Despite the simulated positive future trend in SAM, such a regional feature around the Antarctic Peninsula is absent in the projected sea-ice change for the end of the twenty-first century. The maximum decrease is indeed located over the central Weddell Sea and the Amundsen–Bellingshausen Seas. In most models, changes in the oceanic currents could play a role in the regional distribution of the sea ice, especially in the Ross Sea, where stronger southward currents could be responsible for a smaller sea-ice decrease during the twenty-first century. Finally, changes in the mixed layer depth can be found in some models, inducing locally strong changes in the sea-ice concentration.

Modelling of near-surface ozone over South Asia

Hourly, three-dimensional, fields of tropospheric ozone have been produced for 12 consecutive months on a domain covering South Asia, using the regional Eulerian off-line chemistry transport model MATCH. The results were compared with background observations to investigate diurnal and seasonal variations of near-surface ozone in the region. MATCH reproduced the seasonality of near-surface ozone at most locations in the area. However, the current, and previous, studies indicate that the model consequently overestimate night-time concentrations, while it occasionally underestimates the day-time, near-surface, ozone concentrations. The lowest monthly-mean concentrations of near-surface ozone are typically experienced in June–September, coincident with the rainy season in most areas. The seasonality is not identical across the domain; some locations have a completely different trend. Large areas in Northern India and Nepal show a secondary minimum during the cold winter season (December–January). High concentrations of near-surface ozone are found over the oceans, close to the Indian subcontinent, due to the less efficient dry deposition to water surfaces; over parts of Tibet due to influence of free tropospheric air and little deposition to snow covered surfaces; and along the Gangetic valley due to the large emissions of precursors in this region. Monthly-mean ozone concentrations in the densely populated northern India range from 30–45 ppb(v). The model results were also used to produce maps of AOT40. The results point towards similar levels of AOT40 in India as in Europe: large areas of India show 3-month AOT40 values above 3 ppm(v) hours.

Future changes in cyclone climatology over Europe as inferred from a regional climate simulation

This study analyzes the cyclone climatology in regional climate model simulations of present day (1961–1990) and future (2071–2100, A2 and B2 emission scenarios) european climate conditions. The model domain covers the area from Scandinavia to Northern Africa and from the Eastern Atlantic to Russia at a horizontal grid spacing of 50 km. Compared to present day, in the A2 and B2 scenario conditions the annual average storm track intensity increases over the North-East Atlantic and decreases over Russia and the Eastern Mediterranean region. This overall change pattern is larger in the A2 than in the B2 simulations. However, the cyclone climatology change signal shows a large intermonthly variability and important differences across European regions. The largest changes are found over the North-East Atlantic, where the storm track intensity increases in winter and decreases in summer. A significant reduction of storm track intensity is found during late summer and autumn over the Mediterranean region, and from October to January over Russia. The number of cyclones decreases in future conditions throughout Europe, except over the Central Europe and Mediterranean regions in summer (where it increases). The frequency of intense cyclones and the depth of extreme cyclones increase over the North-East Atlantic, decrease over Russia and show an irregular response over the rest of the domain.

---

Will the tropical land biosphere dominate the climate–carbon cycle feedback during the twenty-first century?

Global warming caused by anthropogenic CO₂ emissions is expected to reduce the capability of the ocean and the land biosphere to take up carbon. This will enlarge the fraction of the CO₂ emissions remaining in the atmosphere, which in turn will reinforce future climate change. Recent model studies agree in the existence of such a positive climate–carbon cycle feedback, but the estimates of its amplitude differ by an order of magnitude, which considerably increases the uncertainty in future climate projections. Therefore we discuss, in how far a particular process or component of the carbon cycle can be identified, that potentially contributes most to the positive feedback. The discussion is based on simulations with a carbon cycle model, which is embedded in the atmosphere/ocean general circulation model ECHAM5/MPI-OM. Two simulations covering the period 1860–2100 are conducted to determine the impact of global warming on the carbon cycle. Forced by historical and future carbon dioxide emissions (following the scenario A2 of the Intergovernmental Panel on Climate Change), they reveal a noticeable positive climate–carbon cycle feedback, which is mainly driven by the tropical land biosphere. The oceans contribute much less to the positive feedback and the temperate/boreal terrestrial biosphere induces a minor negative feedback. The contrasting behavior of the tropical and temperate/boreal land biosphere is mostly attributed to opposite trends in their net primary productivity (NPP) under global warming conditions. As these findings depend on the model employed they are compared with results derived from other climate–carbon cycle models, which participated in the Coupled Climate–Carbon Cycle Model Intercomparison Project (C4MIP).

The relative influence of soil moisture and SST in climate predictability explored within ensembles of AMIP type experiments

Three ensembles of AMIP-type simulations using the Arpege-climat coupled land–atmosphere model have been designed to assess the relative influence of SST and soil moisture (SM) on climate variability and predictability. The study takes advantage of the GSWP2 land surface reanalysis covering the 1986–1995 period. The GSWP2 forcings have been used to derive a global SM climatology that is fully consistent with the model used in this study. One ensemble of ten simulations has been forced by climatological SST and the simulated SM is relaxed toward the GSWP2 reanalysis. Another ensemble has been forced by observed SST and SM is evolving freely. The last ensemble combines the observed SST forcing and the relaxation toward GSWP2. Two complementary aspects of the predictability have been explored, the potential predictability (analysis of variance) and the effective predictability (skill score). An analysis of variance has revealed the effects of the SST and SM boundary forcings on the variability and potential predictability of near-surface temperature, precipitation and surface evaporation. While in the tropics SST anomalies clearly maintain a potentially predictable variability throughout the annual cycle, in the mid-latitudes the SST forced variability is only dominant in winter and SM plays a leading role in summer. In a similar fashion, the annual cycle of the hindcast skill (evaluated as the anomalous correlation coefficient of the three ensemble means with respect to the “observations”) indicates that the SST forcing is the dominant contributor over the tropical continents and in the winter mid-latitudes but that SM is supporting a significant part of the skill in the summer mid-latitudes. Focusing on boreal summer, we have then investigated different aspects of the SM and SST contribution to climate variations in terms of spatial distribution and time-evolution. Our experiments suggest that SM is potentially an additional source of climate predictability. A realistic initialization of SM and a proper representation of the land–atmosphere feedbacks seem necessary to improve state-of-the-art dynamical seasonal predictions, but will be actually efficient only in the areas where SM anomalies are themselves predictable at the monthly to seasonal timescale (since remote effects of SM are probably much more limited than SST teleconnections).     

Cycles and shifts: 1,300 years of multi-decadal temperature variability in the Gulf of Alaska

The Gulf of Alaska (GOA) is highly sensitive to shifts in North Pacific climate variability. Here we present an extended tree-ring record of January–September GOA coastal surface air temperatures using tree-ring width data from coniferous trees growing in the mountain ranges along the GOA. The reconstruction (1514–1999), based on living trees, explains 44% of the temperature variance, although, as the number of chronologies decreases back in time, this value decreases to, and remains around ∼30% before 1840. Verification of the calibrated models is, however, robust. Utilizing sub-fossil wood, we extend the GOA reconstruction back to the early eighth century. The GOA reconstruction correlates significantly (95% CL) with both the Pacific Decadal Oscillation Index (0.53) and North Pacific Index (−0.42) and therefore likely yields important information on past climate variability in the North Pacific region. Intervention analysis on the GOA reconstruction identifies the known twentieth century climate shifts around the 1940s and 1970s, although the mid-1920s shift is only weakly expressed. In the context of the full 1,300 years record, the well studied 1976 shift is not unique. Multi-taper method spectral analysis shows that the spectral properties of the living and sub-fossil data are similar, with both records showing significant (95% CL) spectral peaks at ∼9–11, 13–14 and 18–19 years. Singular spectrum analysis identifies (in order of importance) significant oscillatory modes at 18.7, 50.4, 38.0, 91.8, 24.4, 15.3 and 14.1 years. The amplitude of these modes varies through time. It has been suggested (Minobe in Geophys Res Lett 26:855–858, 1999) that the regime shifts during the twentieth century can be explained by the interaction between pentadecadal (50.4 years) and bidecadal (18.7 years) oscillatory modes. Removal of these two modes of variance from our GOA time series does indeed remove the twentieth century shifts, but many are still identified prior to the twentieth century. Our analysis suggests that climate variability of the GOA is very complex, and that much more work is required to understand the underlying oscillatory behavior that is observed in instrumental and proxy records from the North Pacific region.

---

Past and future changes in climate and hydrological indicators in the US Northeast

To assess the influence of global climate change at the regional scale, we examine past and future changes in key climate, hydrological, and biophysical indicators across the US Northeast (NE). We first consider the extent to which simulations of twentieth century climate from nine atmosphere-ocean general circulation models (AOGCMs) are able to reproduce observed changes in these indicators. We then evaluate projected future trends in primary climate characteristics and indicators of change, including seasonal temperatures, rainfall and drought, snow cover, soil moisture, streamflow, and changes in biometeorological indicators that depend on threshold or accumulated temperatures such as growing season, frost days, and Spring Indices (SI). Changes in indicators for which temperature-related signals have already been observed (seasonal warming patterns, advances in high-spring streamflow, decreases in snow depth, extended growing seasons, earlier bloom dates) are generally reproduced by past model simulations and are projected to continue in the future. Other indicators for which trends have not yet been observed also show projected future changes consistent with a warmer climate (shrinking snow cover, more frequent droughts, and extended low-flow periods in summer). The magnitude of temperature-driven trends in the future are generally projected to be higher under the Special Report on Emission Scenarios (SRES) mid-high (A2) and higher (A1FI) emissions scenarios than under the lower (B1) scenario. These results provide confidence regarding the direction of many regional climate trends, and highlight the fundamental role of future emissions in determining the potential magnitude of changes we can expect over the coming century.

---

Transient simulation of the last glacial inception. Part I: glacial inception as a bifurcation in the climate system

We study the mechanisms of glacial inception by using the Earth system model of intermediate complexity, CLIMBER-2, which encompasses dynamic modules of the atmosphere, ocean, biosphere and ice sheets. Ice-sheet dynamics are described by the three-dimensional polythermal ice-sheet model SICOPOLIS. We have performed transient experiments starting at the Eemiam interglacial, at 126 ky BP (126,000 years before present). The model runs for 26 kyr with time-dependent orbital and CO₂ forcings. The model simulates a rapid expansion of the area covered by inland ice in the Northern Hemisphere, predominantly over Northern America, starting at about 117 kyr BP. During the next 7 kyr, the ice volume grows gradually in the model at a rate which corresponds to a change in sea level of 10 m per millennium. We have shown that the simulated glacial inception represents a bifurcation transition in the climate system from an interglacial to a glacial state caused by the strong snow-albedo feedback. This transition occurs when summer insolation at high latitudes of the Northern Hemisphere drops below a threshold value, which is only slightly lower than modern summer insolation. By performing long-term equilibrium runs, we find that for the present-day orbital parameters at least two different equilibrium states of the climate system exist—the glacial and the interglacial; however, for the low summer insolation corresponding to 115 kyr BP, we find only one, glacial, equilibrium state, while for the high summer insolation corresponding to 126 kyr BP only an interglacial state exists in the model.

---

Internal variability,external forcing and climate trends in multi-decadal AGCM ensembles

An atmospheric general circulation model of intermediate complexity is used to investigate the origin and structure of the climate change in the second half of the twentieth century. The variability of the atmospheric flow is considered as a superposition of an internal part, due to intrinsic dynamical variability, and an external part, due to the variations of the sea surface temperature (SST) forcing. The two components are identified by performing a 50-member ensemble of atmospheric simulations with prescribed, observed SSTs in the period 1949–2002. The large number of realizations allows the estimation of statistics of the atmospheric variability with a high confidence level. The analysis performed focuses on interdecadal and interannual variability of 500 hPa geopotential height in the Northern Hemisphere (NH) during winter. The model reproduces well the structure of the observed trend (defined as the difference in the two 25-year intervals 1977–2001 and 1952–1976), particularly in the Pacific region, and about half of the amplitude of the signal. The trend in 500 hPa height projects mainly onto the second empirical orthogonal function (EOF), both in the observations and in the model ensemble. However, differences between the modelled and the observed variability are found in the pattern of the second EOF in the Atlantic sector. SST changes associated with the El Niño southern oscillation (ENSO) are responsible for about 50% of the signal of the 500 hPa height trend in the Pacific. A second 50-member ensemble is used to evaluate the sensitivity of interdecadal variability to an increase in CO₂ optical depth compatible with observed concentration changes. In this second experiment, the simulated trend includes a statistically significant contribution from the positive phase of the Arctic oscillation (AO). Such a contribution is also found in observations. Furthermore, the additional CO₂ forcing accounts for part of the NH trend in near-surface temperature, and brings the zonal-mean temperature changes in the stratosphere and upper-troposphere closer to observations.

---

The impact of natural and anthropogenic forcings on climate and hydrology since 1550

A climate simulation of an ocean/atmosphere general circulation model driven with natural forcings alone (constant “pre-industrial” land-cover and well-mixed greenhouse gases, changing orbital, solar and volcanic forcing) has been carried out from 1492 to 2000. Another simulation driven with natural and anthropogenic forcings (changes in greenhouse gases, ozone, the direct and first indirect effect of anthropogenic sulphate aerosol and land-cover) from 1750 to 2000 has also been carried out. These simulations suggest that since 1550, in the absence of anthropogenic forcings, climate would have warmed by about 0.1 K. Simulated response is not in equilibrium with the external forcings suggesting that both climate sensitivity and the rate at which the ocean takes up heat determine the magnitude of the response to forcings since 1550. In the simulation with natural forcings climate sensitivity is similar to other simulations of HadCM3 driven with CO₂ alone. Climate sensitivity increases when anthropogenic forcings are included. The natural forcing used in our experiment increases decadal–centennial time-scale and large spatial scale climate variability, relative to internal variability, as diagnosed from a control simulation. Mean conditions in the natural simulation are cooler than in our control simulation reflecting the reduction in forcing. However, over certain regions there is significant warming, relative to control, due to an increase in forest cover. Comparing the simulation driven by anthropogenic and natural forcings with the natural-only simulation suggests that anthropogenic forcings have had a significant impact on, particularly tropical, climate since the early nineteenth century. Thus the entire instrumental temperature record may be “contaminated” by anthropogenic influences. Both the hydrological cycle and cryosphere are also affected by anthropogenic forcings. Changes in tree-cover appear to be responsible for some of the local and hydrological changes as well as an increase in northern hemisphere spring snow cover.

---

The greening of the McGill Paleoclimate Model. Part I: Improved land surface scheme with vegetation dynamics

The formulation of a new land surface scheme (LSS) with vegetation dynamics for coupling to the McGill Paleoclimate Model (MPM) is presented. This LSS has the following notable improvements over the old version: (1) parameterization of deciduous and evergreen trees by using the models climatology and the output of the dynamic global vegetation model, VECODE (Brovkin et al. in Ecological Modelling 101:251–261 (1997), Global Biogeochemical Cycles 16(4):1139, (2002)); (2) parameterization of tree leaf budburst and leaf drop by using the models climatology; (3) parameterization of the seasonal cycle of the grass leaf area index; (4) parameterization of the seasonal cycle of tree leaf area index by using the time-dependent growth of the leaves; (5) calculation of land surface albedo by using vegetation-related parameters, snow depth and the models climatology. The results show considerable improvement of the models simulation of the present-day climate as compared with that simulated in the original physically-based MPM. In particular, the strong seasonality of terrestrial vegetation and the associated land surface albedo variations are in good agreement with several satellite observations of these quantities. The application of this new version of the MPM (the green MPM) to Holocene millennial-scale climate changes is described in a companion paper, Part II.

---

Sensitivity of sea ice to wind-stress and radiative forcing since 1500: a model study of the Little Ice Age and beyond

Three different reconstructed wind-stress fields which take into account variations of the North Atlantic Oscillation, one general circulation model wind-stress field, and three radiative forcings (volcanic activity, insolation changes and greenhouse gas changes) are used with the UVic Earth System Climate Model to simulate the surface air temperature, the sea-ice cover, and the Atlantic meridional overturning circulation (AMOC) since 1500, a period which includes the Little Ice Age (LIA). The simulated Northern Hemisphere surface air temperature, used for model validation, agrees well with several temperature reconstructions. The simulated sea-ice cover in each hemisphere responds quite differently to the forcings. In the Northern Hemisphere, the simulated sea-ice area and volume during the LIA are larger than the present-day area and volume. The wind-driven changes in sea-ice area are about twice as large as those due to thermodynamic (i.e., radiative) forcing. For the sea-ice volume, changes due to wind forcing and thermodynamics are of similar magnitude. Before 1850, the simulations suggest that volcanic activity was mainly responsible for the thermodynamically produced area and volume changes, while after 1900 the slow greenhouse gas increase was the main driver of the sea-ice changes. Changes in insolation have a small effect on the sea ice throughout the integration period. The export of the thicker sea ice during the LIA has no significant effect on the maximum strength of the AMOC. A more important process in altering the maximum strength of the AMOC and the sea-ice thickness is the wind-driven northward ocean heat transport. In the Southern Hemisphere, there are no visible long-term trends in the simulated sea-ice area or volume since 1500. The wind-driven changes are roughly four times larger than those due to radiative forcing. Prior to 1800, all the radiative forcings could have contributed to the thermodynamically driven changes in area and volume. In the 1800s the volcanic forcing was dominant, and during the first part of the 1900s both the insolation changes and the greenhouse gas forcing are responsible for thermodynamically produced changes. Finally, in the latter part of the 1900s the greenhouse gas forcing is the dominant factor in determining the sea-ice changes in the Southern Hemisphere.