Results of PMIP2 coupled simulations of the mid-Holocene and Last Glacial Maximum, Part 1: experiments and large-scale features

A set of coupled ocean-atmosphere simulations using state of the art climate models is now available for the Last Glacial Maximum and the Mid-Holocene through the second phase of the Paleoclimate Modeling Intercomparison Project (PMIP2). This study presents the large-scale features of the simulated climates and compares the new model results to those of the atmospheric models from the first phase of the PMIP, for which sea surface temperature was prescribed or computed using simple slab ocean formulations. We consider the large-scale features of the climate change, pointing out some of the major differences between the different sets of experiments. We show in particular that systematic differences between PMIP1 and PMIP2 simulations are due to the interactive ocean, such as the amplification of the African monsoon at the Mid-Holocene or the change in precipitation in mid-latitudes at the LGM. Also the PMIP2 simulations are in general in better agreement with data than PMIP1 simulations.

Technical note: a simple procedure for removing temporal discontinuities in ERA-Interim upper stratospheric temperatures for use in nudged chemistry-climate model simulations

This note describes a simple procedure for removing unphysical temporal discontinuities in ERA-Interim upper stratospheric global mean temperatures in March 1985 and August 1998 that have arisen due to changes in satellite radiance data used in the assimilation. The derived temperature adjustments (offsets) are suitable for use in stratosphere-resolving chemistry-climate models that are nudged (relaxed) to ERA-Interim winds and temperatures. Simulations using a nudged version of the Canadian Middle Atmosphere Model (CMAM) show that the inclusion of the temperature adjustments produces temperature time series that are devoid of the large jumps in 1985 and 1998. Due to its strong temperature dependence, the simulated upper stratospheric ozone is also shown to vary smoothly in time, unlike in a nudged simulation without the adjustments where abrupt changes in ozone occur at the times of the temperature jumps. While the adjustments to the ERA-Interim temperatures remove significant artefacts in the nudged CMAM simulation, spurious transient effects that arise due to water vapour and persist for about 5 yr after the 1979 switch to ERA-Interim data are identified, underlining the need for caution when analysing trends in runs nudged to reanalyses.

How well can we simulate past climates? Evaluating the models using global palaeoenvironmental datasets.

Global syntheses of palaeoenvironmental data are required to test climate models under conditions different from the present. Data sets for this purpose contain data from spatially extensive networks of sites. The data are either directly comparable to model output or readily interpretable in terms of modelled climate variables. Data sets must contain sufficient documentation to distinguish between raw (primary) and interpreted (secondary, tertiary) data, to evaluate the assumptions involved in interpretation of the data, to exercise quality control, and to select data appropriate for specific goals. Four data bases for the Late Quaternary, documenting changes in lake levels since 30 kyr BP (the Global Lake Status Data Base), vegetation distribution at 18 kyr and 6 kyr BP (BIOME 6000), aeolian accumulation rates during the last glacial-interglacial cycle (DIRTMAP), and tropical terrestrial climates at the Last Glacial Maximum (the LGM Tropical Terrestrial Data Synthesis) are summarised. Each has been used to evaluate simulations of Last Glacial Maximum (LGM: 21 calendar kyr BP) and/or mid-Holocene (6 cal. kyr BP) environments. Comparisons have demonstrated that changes in radiative forcing and orography due to orbital and ice-sheet variations explain the first-order, broad-scale (in space and time) features of global climate change since the LGM. However, atmospheric models forced by 6 cal. kyr BP orbital changes with unchanged surface conditions fail to capture quantitative aspects of the observed climate, including the greatly increased magnitude and northward shift of the African monsoon during the early to mid-Holocene. Similarly, comparisons with palaeoenvironmental datasets show that atmospheric models have underestimated the magnitude of cooling and drying of much of the land surface at the LGM. The inclusion of feedbacks due to changes in ocean- and land-surface conditions at both times, and atmospheric dust loading at the LGM, appears to be required in order to produce a better simulation of these past climates. The development of Earth system models incorporating the dynamic interactions among ocean, atmosphere, and vegetation is therefore mandated by Quaternary science results as well as climatological principles. For greatest scientific benefit, this development must be paralleled by continued advances in palaeodata analysis and synthesis, which in turn will help to define questions that call for new focused data collection efforts.

Hurricanes and climate: the U.S. CLIVAR working group on hurricanes

While a quantitative climate theory of tropical cyclone formation remains elusive, considerable progress has been made recently in our ability to simulate tropical cyclone climatologies and understand the relationship between climate and tropical cyclone formation. Climate models are now able to simulate a realistic rate of global tropical cyclone formation, although simulation of the Atlantic tropical cyclone climatology remains challenging unless horizontal resolutions finer than 50 km are employed. This article summarizes published research from the idealized experiments of the Hurricane Working Group of U.S. CLIVAR (CLImate VARiability and predictability of the ocean-atmosphere system). This work, combined with results from other model simulations, has strengthened relationships between tropical cyclone formation rates and climate variables such as mid-tropospheric vertical velocity, with decreased climatological vertical velocities leading to decreased tropical cyclone formation. Systematic differences are shown between experiments in which only sea surface temperature is increased versus experiments where only atmospheric carbon dioxide is increased, with the carbon dioxide experiments more likely to demonstrate the decrease in tropical cyclone numbers previously shown to be a common response of climate models in a warmer climate. Experiments where the two effects are combined also show decreases in numbers, but these tend to be less for models that demonstrate a strong tropical cyclone response to increased sea surface temperatures. Further experiments are proposed that may improve our understanding of the relationship between climate and tropical cyclone formation, including experiments with two-way interaction between the ocean and the atmosphere and variations in atmospheric aerosols.

A lagged response to the 11 year solar cycle in observed winter Atlantic/European weather patterns

The surface response to 11 year solar cycle variations is investigated by analyzing the long-term mean sea level pressure and sea surface temperature observations for the period 1870–2010. The analysis reveals a statistically significant 11 year solar signal over Europe, and the North Atlantic provided that the data are lagged by a few years. The delayed signal resembles the positive phase of the North Atlantic Oscillation (NAO) following a solar maximum. The corresponding sea surface temperature response is consistent with this. A similar analysis is performed on long-term climate simulations from a coupled ocean-atmosphere version of the Hadley Centre model that has an extended upper lid so that influences of solar variability via the stratosphere are well resolved. The model reproduces the positive NAO signal over the Atlantic/European sector, but the lag of the surface response is not well reproduced. Possible mechanisms for the lagged nature of the observed response are discussed.

Respostas dinâmicas em sistemas distribuídos e decomposição forçada da superfície livre para um modelo acoplado Oceano-Atmosfera

O objetivo deste trabalho é a introdução e desenvolvimento de uma metodologia analítico-simbólica para a obtenção de respostas dinâmicas e forçadas (soluções homogêneas e não homogêneas) de sistemas distribuídos, em domínios ilimitados e limitados, através do uso da base dinâmica gerada a partir da resposta impulso. Em domínios limitados, a resposta impulso foi formulada pelo método espectral. Foram considerados sistemas com condições de contorno homogêneas e não homogêneas. Para sistemas de natureza estável, a resposta forçada é decomposta na soma de uma resposta particular e de uma resposta livre induzida pelos valores iniciais da resposta particular. As respostas particulares, para entradas oscilatórias no tempo, foram calculadas com o uso da fun»c~ao de Green espacial. A teoria é desenvolvida de maneira geral permitindo que diferentes sis- temas evolutivos de ordem arbitrária possam ser tratados sistematicamente de uma forma compacta e simples. Realizou-se simulações simbólicas para a obtenção de respostas dinâmicas e respostas for»cadas com equações do tipo parabólico e hiperbólico em 1D,2D e 3D. O cálculo das respostas forçadas foi realizado com a determinação das respostas livres transientes em termos dos valores iniciais das respostas permanentes. Foi simulada a decomposição da resposta forçada da superfície livre de um modelo acoplado oceano-atmosfera bidimensional, através da resolução de uma equação de Klein-Gordon 2D com termo não-homogêneo de natureza dinâmica, devido a tensão de cisalhamento na superfície do oceano pela ação do vento.

Sazonalidade da precipitação para a Amazônia usando o modelo REGCM3: avaliando apenas a forçante do Atlântico Equatorial

O presente estudo tenta trazer um novo ponto de vista sobre a importância do Oceano Atlântico Equatorial no padrão sazonal, durante o período de verão e de outono da região Amazônica. Desta maneira, contribui de uma forma diferente para os estudos na área de modelagem climática e variabilidade pluviométrica na região. Para isto, foi utilizado o modelo RegCM3, com esquema de convecção Grell, aplicando a técnica donwscaling e utilizando como condição inicial os dados de reanalise do NCEP/NCAR. Após a obtenção das simulações foi aplicado o método estatístico Viés para avaliar o quanto o modelo consegue reproduzir a precipitação sazonal. Em uma primeira análise, constata-se que o modelo é sensível a entrada de SF ou ZCAS, bem como, aproxima-se da dinâmica dos trópicos e extratropicos. Os resultados sugerem que a baixa resolução e a parametrização Grell são os principais fatores para a inadequação da qualidade das simulações. Apesar de o esquema Grell ser indicado para regiões com intensa convecção e fortes movimentos verticais, o mesmo necessita ser ajustado com as características físicas da região. Este estudo contribuiu para o aprimoramento dos modelos climáticos regionais para a região amazônica, considerando a contribuição dos processos oceano-atmosfera, durante o período de verão e outono do hemisfério sul.

Accumulation of semivolatile organic compounds in Antarctic vegetation: A case study of polybrominated diphenyl ethers

Antarctic plant communities are dominated by lichens and mosses which accumulate semivolatile organic compounds (SOCs) such as polybrominated diphenyl ethers (PBDEs) directly from the atmosphere. Differences in the levels of PBDEs observed in lichens and mosses collected at King George Island in the austral summers 2004-05 and 2005-06 are probably explained by environmental and/or plant parameters. Contamination of lichens showed a positive correlation with local precipitation, suggesting that wet deposition processes are a major mechanism controlling the uptake of most PBDE congeners. These findings are in agreement with physical-chemical data supporting that tetra- through hepta-BDEs in the Antarctic atmosphere are basically bound to aerosols. Conversely, accumulation of PBDEs in mosses appears to be controlled by other environmental factors and/or plant-specific characteristics. Model simulations demonstrated that an ocean-atmosphere coupling may have played a role in the long-range transport of less volatile SOCs such as PBDEs to Antarctica. According to simulations, the atmosphere is the most important transport medium for PBDEs while the surface ocean serves as a temporary storage compartment, boosting the deposition/volatilization ""hopping"" effect similarly to vegetation on continents. (C) 2011 Elsevier B.V. All rights reserved.

Reproducibility of South American precipitation due to subtropical South Atlantic SSTs

This work investigates the eproducibility of precipitation simulated with an atmospheric general circulation model (AGCM) forced by subtropical South Atlantic sea surface temperature (SST) anomalies. This represents an important test of the model prior to investigating the impact of SSTs on regional climate. A five-member ensemble run was performed using the National Center for Atmospheric Research (NCAR) Community Climate Model, version 3 (CCM3). The CCM3 was forced by observed monthly SST over the South Atlantic from 20 to 60 S. The SST dataset used is from the Hadley Centre covering the period of September 1949-October 2001; this covers more than 50 yr of simulation. A statistical technique is used to determine the reproducibility in the CCM3 runs and to assess potential predictability in precipitation. Empirical orthogonal function analysis is used to reconstruct the ensemble using the most reproducible forced modes in order to separate the atmospheric response to local SST forcing from its internal variability. Results for reproducibility show a seasonal dependence, with higher values during austral autumn and spring. The spatial distribution of reproducibility shows that the tropical atmosphere is dominated by the underlying SSTs while variations in the subtropical-extratropical regions are primarily driven by internal variability. As such, changes in the South Atlantic convergence zone (SACZ) region are mainly dominated by internal atmospheric variability while the ITCZ has greater external dependence, making it more predictable. The reproducibility distribution reveals increased values after the reconstruction of the ensemble.

Causes and consequences of past and projected Scandinavian summer temperatures, 500-2100 AD

Tree rings dominate millennium-long temperature reconstructions and many records originate from Scandinavia, an area for which the relative roles of external forcing and internal variation on climatic changes are, however, not yet fully understood. Here we compile 1,179 series of maximum latewood density measurements from 25 conifer sites in northern Scandinavia, establish a suite of 36 subset chronologies, and analyse their climate signal. A new reconstruction for the 1483–2006 period correlates at 0.80 with June–August temperatures back to 1860. Summer cooling during the early 17th century and peak warming in the 1930s translate into a decadal amplitude of 2.9°C, which agrees with existing Scandinavian tree-ring proxies. Climate model simulations reveal similar amounts of mid to low frequency variability, suggesting that internal ocean-atmosphere feedbacks likely influenced Scandinavian temperatures more than external forcing. Projected 21st century warming under the SRES A2 scenario would, however, exceed the reconstructed temperature envelope of the past 1,500 years.

Sr/Ca ratios of coccoliths from the Paleocene-Eocene Thermal Maximum

During the Paleocene-Eocene Thermal Maximum (PETM), rapid release of isotopically light C to the ocean-atmosphere system elevated the greenhouse effect and warmed temperatures by 5-7 °C for 105 yr. The response of the planktic ecosystems and productivity to the dramatic climate changes of the PETM may represent a significant feedback to the carbon cycle changes, but has been difficult to document. We examine Sr/Ca ratios in calcareous nannofossils in sediments spanning the PETM in three open ocean sites as a new approach to examine productivity and ecological shifts in calcifying plankton. The large heterogeneity in Sr/Ca among different nannofossil genera indicates that nannofossil Sr/Ca reflects primary productivity-driven geochemical signals and not diagenetic overprinting. Elevated Sr/Ca ratios in several genera and constant ratios in other genera suggest increased overall productivity in the Atlantic sector of the Southern Ocean during the PETM. Dominant nannofossil genera in tropical Atlantic and Pacific sites show Sr/Ca variations during the PETM which are comparable to background variability prior to the PETM. Despite acidification of the ocean there was not a productivity crisis among calcifying phytoplankton. We use the Pandora ocean box model to explore possible mechanisms for PETM productivity change. If independent proxy evidence for more stratified conditions in the Southern Ocean during the PETM is robust, then maintenance of stable or increased productivity there likely reflects increased nutrient inventories of the ocean. Increased nutrient inventories could have resulted from climatically enhanced weathering and would have important implications for burial rates of organic carbon and stabilization of climate and the carbon cycle.

Modern and fossil corals from Palmyra and Christmas islands

The relationship between decadal to centennial changes in ocean circulation and climate is difficult to discern using the sparse and discontinuous instrumental record of climate and, as such, represents a large uncertainty in coupled ocean-atmosphere general circulation models. We present new modern and fossil coral radiocarbon (D14C) records from Palmyra (6°N, 162°W) and Christmas (2°N, 157°W) islands to constrain central tropical Pacific ocean circulation changes during the last millennium. Seasonally to annually resolved coral D14C measurements from the 10th, 12th-17th, and 20th centuries do not contain significant interannual to decadal-scale variations, despite large changes in coral d18O on these timescales. A centennial-scale increase in coral radiocarbon from the Medieval Climate Anomaly (~900-1200 AD) to the Little Ice Age (~1500-1800) can be largely explained by changes in the atmospheric D14C, as determined with a box model of Palmyra mixed layer D14C. However, large 12th century depletions in Palmyra coral D14C may reflect as much as a 100% increase in upwelling rates and/or a significant decrease in the D14C of higher-latitude source waters reaching the equatorial Pacific during this time. SEM photos reveal evidence for minor dissolution and addition of secondary aragonite in the fossil corals, but our results suggest that coral D14C is only compromised after moderate to severe diagenesis for these relatively young fossil corals.

Results of modelled Arctic river discharge trends since 7 ka BP

The Arctic hydrological cycle throughout the Holocene is analyzed based on the results of transient simulations with the coupled atmosphere-ocean circulation model ECHO-G. The results suggest a ~ 2 % increase of mid-Holocene to preindustrial Arctic river discharges for the Eurasian continent. However, rivers of the North America Arctic realm show a moderate runoff decline of approximately 4 to 5 % for the same period. The total river discharge into the Arctic Ocean has remained at an approximately constant preindustrial level since the mid Holocene. The positive discharge trend within Eurasia is caused by a more rapid decrease in local net evaporation compared to a smaller decline in advected moisture and hence precipitation. This effect is neither recognized within the North American Arctic domain nor in the far eastern part of the Eurasian Arctic realm. A detailed comparison of these model findings with a variety of proxy studies is conducted. The collected proxy records show trends of continental surface temperatures and precipitation rates that are consistent with the simulations. A continuation of the transient Holocene runs for the 19th and 20th century with increased greenhouse gases indicates an increase of the total river influx into the Arctic Ocean of up to 7.6 %. The Eurasian river discharges increase by 7.5 %, the North American discharges by up to 8.4 %. The most rapid increases have been detected since the beginning of the 20th century. These results are corroborated by the observed rising of Arctic river discharges during the last century which is attributed to anthropogenic warming. The acceleration of the Arctic hydrological cycle in the 20th century is without precedence in the Holocene.

Age models and stable isotope ratios from the Mid-Miocene

Mid-Miocene pelagic sedimentary sections can be correlated using intermediate and high resolution oxygen and carbon isotopic records of benthic foraminifera. Precision of a few tens of thousands of years is readily achievable at sites with high sedimentation rates, for example, Deep Sea Drilling Project sites 289 and 574. The mid-Miocene carbon isotope records are characterized by an interval of high d13C values between 17 and 13.5 Ma (the Monterey Excursion of Vincent and Berger 1985) upon which are superimposed a series of periodic or quasi-periodic fluctuations in d13C values. These fluctuations have a period of approximately 440 kyr, suggestive of the 413 kyr cycle predicted by Milankovitch theory. Vincent and Berger proposed that the Monterey Excursion was the result of increased organic carbon burial in continental margins sediments. The increased d13C values (called 13C maxima) superimposed on the generally high mid-Miocene signal coincide with increases in d18O values suggesting that periods of cooling and/or ice buildup were associated with exceptionally rapid burial of organic carbon and lowered atmospheric CO2 levels. It is likely that during the Monterey Excursion the ocean/atmosphere system became progressively more sensitive to small changes in insolation, ultimately leading to major cooling of deep water and expansion of continental ice. We have assigned an absolute chronology, based on biostratigraphic and magneto-biostratigraphic datum levels, to the isotope stratigraphy and have used that chronology to correlate unconformities, seismic reflectors, carbonate minima, and dissolution intervals. Intervals of sediment containing 13C maxima are usually better preserved than the overlying and underlying sediments, indicating that the d13C values of TCO2 in deep water and the corrosiveness of seawater are inversely correlated. This again suggests that the 13C maxima were associated with rapid burial of organic carbon and reduced levels of atmospheric CO2. The absolute chronology we have assigned to the isotopic record indicates that the major mid-Miocene deepwater cooling/ice volume expansion took 2 m.y. and was not abrupt as had been reported previously. The cooling appears abrupt at many sites because the interval is characterized by a number of dissolution intervals. The cooling was not monotonic, and the 2 m.y. interval included an episode of especially rapid cooling as well as a brief return to warmer conditions before the final phase of the cooling period. The increase in d18O values of benthic foraminifera between 14.9 and 12.9 Ma was greatest at deeper water sites and at sites closest to Antarctica. The data suggest that the d18O value of seawater increased by no more than about 1.1 per mil during this interval and that the remainder of the change in benthic d18O values resulted from cooling in Antarctic regions of deepwater formation. Equatorial planktonic foraminifera from sites 237 and 289 exhibit a series of 0.4 per mil steplike increases in d13C values. Only one of these increases in planktonic d13C is correlated with any of the features in the mid-Miocene benthic carbon isotope record.

Stable isotope ratios of foraminifera from sediment cores off northeastern Brazil

The western South Atlantic boundary currents represent a sensitive system within the global thermohaline circulation (THC). We investigated the impact of deglacial THC changes on the western tropical Atlantic studied in six high resolution sediment cores from the upper continental slope of Brazil. The stratigraphy of the cores is mainly based on 14C AMS dating of monospecific foraminiferal samples. Changes in the upper layer tropical ocean during the deglaciation are inferred from stable oxygen isotope measurements on planktic and benthic foraminifera. Variations in the delta18O residuals are assumed to be mainly temperature related. During the Oldest and Younger Dryas cooling periods, two major deglacial THC disturbances are reported from North Atlantic sediment cores. Concomitant to the repeated THC slowdown, we observe an upper layer warming in the tropical ocean. A reduced northward heat export from the tropical areas during these periods (weak North Brazil Current) is additionally reflected by low meridional gradients in the stable oxygen records. This generally agrees with results from coupled ocean atmosphere models.