Using synoptic type analysis to understand New Zealand climate during the Mid-Holocene

Diagnosing the climate of New Zealand from low-resolution General Circulation Models (GCMs) is notoriously difficult due to the interaction of the complex topography and the Southern Hemisphere (SH) mid-latitude westerly winds. Therefore, methods of downscaling synoptic scale model data for New Zealand are useful to help understand past climate. New Zealand also has a wealth of palaeoclimate-proxy data to which the downscaled model output can be compared, and to provide a qualitative method of assessing the capability of GCMs to represent, in this case, the climate 6000 yr ago in the Mid-Holocene. In this paper, a synoptic weather and climate regime classification system using Empirical Orthogonal Function (EOF) analysis of GCM and reanalysis data was used. The climate regimes are associated with surface air temperature and precipitation anomalies over New Zealand. From the analysis in this study, we find at 6000 BP that increased trough activity in summer and autumn led to increased precipitation, with an increased north-south pressure gradient ("zonal events") in winter and spring leading to drier conditions. Opposing effects of increased (decreased) temperature are also seen in spring (autumn) in the South Island, which are associated with the increased zonal (trough) events; however, the circulation induced changes in temperature are likely to have been of secondary importance to the insolation induced changes. Evidence from the palaeoclimate-proxy data suggests that the Mid-Holocene was characterized by increased westerly wind events in New Zealand, which agrees with the preference for trough and zonal regimes in the models.

Climate model benchmarking with glacial and mid-Holocene climates

Past climates provide a test of models’ ability to predict climate change. We present a comprehensive evaluation of state-of-the-art models against Last Glacial Maximum and mid-Holocene climates, using reconstructions of land and ocean climates and simulations from the Palaeoclimate Modelling and Coupled Modelling Intercomparison Projects. Newer models do not perform better than earlier versions despite higher resolution and complexity. Differences in climate sensitivity only weakly account for differences in model performance. In the glacial, models consistently underestimate land cooling (especially in winter) and overestimate ocean surface cooling (especially in the tropics). In the mid-Holocene, models generally underestimate the precipitation increase in the northern monsoon regions, and overestimate summer warming in central Eurasia. Models generally capture large-scale gradients of climate change but have more limited ability to reproduce spatial patterns. Despite these common biases, some models perform better than others.

Evaluation of coupled ocean-atmosphere simulations of northern hemisphere extratropical climates in the mid-Holocene

We have used the BIOME4 biogeography–biochemistry model and comparison with palaeovegetation data to evaluate the response of six ocean–atmosphere general circulation models to mid-Holocene changes in orbital forcing in the mid- to high-latitudes of the northern hemisphere. All the models produce: (a) a northward shift of the northern limit of boreal forest, in response to simulated summer warming in high-latitudes. The northward shift is markedly asymmetric, with larger shifts in Eurasia than in North America; (b) an expansion of xerophytic vegetation in mid-continental North America and Eurasia, in response to increased temperatures during the growing season; (c) a northward expansion of temperate forests in eastern North America, in response to simulated winter warming. The northward shift of the northern limit of boreal forest and the northward expansion of temperate forests in North America are supported by palaeovegetation data. The expansion of xerophytic vegetation in mid-continental North America is consistent with palaeodata, although the extent may be over-estimated. The simulated expansion of xerophytic vegetation in Eurasia is not supported by the data. Analysis of an asynchronous coupling of one model to an equilibrium-vegetation model suggests vegetation feedback exacerbates this mid-continental drying and produces conditions more unlike the observations. Not all features of the simulations are robust: some models produce winter warming over Europe while others produce winter cooling. As a result, some models show a northward shift of temperate forests (consistent with, though less marked than, the expansion shown by data) and others produce a reduction in temperate forests. Elucidation of the cause of such differences is a focus of the current phase of the Palaeoclimate Modelling Intercomparison Project.

Mid-Holocene land-surface conditions in northern Africa and the Arabian peninsula: a data set for the analysis of biogeophysical feedbacks in the climate system

Large changes in the extent of northern subtropical arid regions during the Holocene are attributed to orbitally forced variations in monsoon strength and have been implicated in the regulation of atmospheric trace gas concentrations on millenial timescales. Models that omit biogeophysical feedback, however, are unable to account for the full magnitude of African monsoon amplification and extension during the early to middle Holocene (˜9500–5000 years B.P.). A data set describing land-surface conditions 6000 years B.P. on a 1° × 1° grid across northern Africa and the Arabian Peninsula has been prepared from published maps and other sources of palaeoenvironmental data, with the primary aim of providing a realistic lower boundary condition for atmospheric general circulation model experiments similar to those performed in the Palaeoclimate Modelling Intercomparison Project. The data set includes information on the percentage of each grid cell occupied by specific vegetation types (steppe, savanna, xerophytic woods/scrub, tropical deciduous forest, and tropical montane evergreen forest), open water (lakes), and wetlands, plus information on the flow direction of major drainage channels for use in large-scale palaeohydrological modeling.

Mid-holocene and last glacial maximum climate simulations with the IPSL model-part I: comparing IPSLCM5A to IPSLCM4

The climates of the mid-Holocene (MH), 6,000 years ago, and of the Last Glacial Maximum (LGM), 21,000 years ago, have extensively been simulated, in particular in the framework of the Palaeoclimate Modelling Intercomparion Project. These periods are well documented by paleo-records, which can be used for evaluating model results for climates different from the present one. Here, we present new simulations of the MH and the LGM climates obtained with the IPSL_CM5A model and compare them to our previous results obtained with the IPSL_CM4 model. Compared to IPSL_CM4, IPSL_CM5A includes two new features: the interactive representation of the plant phenology and marine biogeochemistry. But one of the most important differences between these models is the latitudinal resolution and vertical domain of their atmospheric component, which have been improved in IPSL_CM5A and results in a better representation of the mid-latitude jet-streams. The Asian monsoon’s representation is also substantially improved. The global average mean annual temperature simulated for the pre-industrial (PI) period is colder in IPSL_CM5A than in IPSL_CM4 but their climate sensitivity to a CO2 doubling is similar. Here we show that these differences in the simulated PI climate have an impact on the simulated MH and LGM climatic anomalies. The larger cooling response to LGM boundary conditions in IPSL_CM5A appears to be mainly due to differences between the PMIP3 and PMIP2 boundary conditions, as shown by a short wave radiative forcing/feedback analysis based on a simplified perturbation method. It is found that the sensitivity computed from the LGM climate is lower than that computed from 2 × CO2 simulations, confirming previous studies based on different models. For the MH, the Asian monsoon, stronger in the IPSL_CM5A PI simulation, is also more sensitive to the insolation changes. The African monsoon is also further amplified in IPSL_CM5A due to the impact of the interactive phenology. Finally the changes in variability for both models and for MH and LGM are presented taking the example of the El-Niño Southern Oscillation (ENSO), which is very different in the PI simulations. ENSO variability is damped in both model versions at the MH, whereas inconsistent responses are found between the two versions for the LGM. Part 2 of this paper examines whether these differences between IPSL_CM4 and IPSL_CM5A can be distinguished when comparing those results to palaeo-climatic reconstructions and investigates new approaches for model-data comparisons made possible by the inclusion of new components in IPSL_CM5A.

Stacking of discontinuous regional palaeoclimate records: Chironomid-based summer temperatures from the Alpine region

Since multi-site reconstructions are less affected by site-specific climatic effects and artefacts, regional palaeotemperature reconstructions based on a number of sites can provide more robust estimates of centennial- to millennial-scale temperature trends than individual, site-specific records. Furthermore, reconstructions based on multiple records are necessary for developing continuous climate records over time scales longer than covered by individual sequences. Here, we present a procedure for developing such reconstructions based on relatively short (centuries to millennia), discontinuously sampled records as are typically developed when using biotic proxies in lake sediments for temperature reconstruction. The approach includes an altitudinal correction of temperatures, an interpolation of individual records to equal time intervals, a stacking procedure for sections of the interval of interest that have the same records available, as well as a splicing procedure to link the individual stacked records into a continuous reconstruction. Variations in the final, stacked and spliced reconstruction are driven by variations in the individual records, whereas the absolute temperature values are determined by the stacked segment based on the largest number of records. With numerical simulations based on the NGRIP δ18O record, we demonstrate that the interpolation and stacking procedure provides an approximation of a smoothed palaeoclimate record if based on a sufficient number of discontinuously sampled records. Finally, we provide an example of a stacked and spliced palaeotemperature reconstruction 15000–90 calibrated 14C yr BP based on six chironomid records from the northern and central Swiss Alps and eastern France to discuss the potential and limitations of this approach.

Southern Hemisphere high-resolution palaeoclimate records of the last 2000 years

This study presents a comprehensive assessment of high-resolution Southern Hemisphere (SH) paleoarchives covering the last 2000 years. We identified 174 monthly to annually resolved climate proxy (tree ring, coral, ice core, documentary, speleothem and sedimentary) records from the Hemisphere. We assess the interannual and decadal sensitivity of each proxy record to large-scale circulation indices from the Pacific, Indian and Southern Ocean regions over the twentieth century. We then analyse the potential of this newly expanded palaeoclimate network to collectively represent predictands (sea surface temperature, sea level pressure, surface air temperature and precipitation) commonly used in climate reconstructions. The key dynamical centres-of-action of the equatorial Indo-Pacific are well captured by the palaeoclimate network, indicating that there is considerable reconstruction potential in this region, particularly in the post AD 1600 period when a number of long coral records are available. Current spatiotemporal gaps in data coverage and regions where significant potential for future proxy collection exists are discussed. We then highlight the need for new and extended records from key dynamical regions of the Southern Hemisphere. Although large-scale climate field reconstructions for the SH are in their infancy, we report that excellent progress in the development of regional proxies now makes plausible estimates of continental- to hemispheric-scale climate variations possible.

Fossil pollen record of composite sediment core TMD from Tso Moriri, analysis of modern surface pollen samples, mean annual precipitation reconstruction, and digitisation and recalibration of different discussed palaeoclimate proxy records

This paper presents a new fossil pollen record from Tso Moriri (32°54'N, 78°19'E, 4512 m a.s.l.) and seeks to reconstruct changes in mean annual precipitation (MAP) during the last 12,000 years. This high-alpine lake occupies an area of 140 km**2 in a glacial-tectonic valley in the northwestern Himalaya. The region has a cold climate, with a MAP <300 mm, and open vegetation. The hydrology is controlled by the Indian Summer Monsoon (ISM), but winter westerly-associated precipitation also affects the regional water balance. Results indicate that precipitation levels varied significantly during the Holocene. After a rapid increase in MAP, a phase of maximum humidity was reached between ca. 11 to 9.6 cal ka BP, followed by a gradual decline in MAP. This trend parallels the reduction in the Northern Hemisphere summer insolation. Comparison of different palaeoclimate proxy records reveal evidence for a stronger Holocene decrease in precipitation in the northern versus the southern parts of the ISM domain. The long-term trend of ISM weakening is overlaid with several short periods of greater dryness, which are broadly synchronous with the North Atlantic cold spells, suggesting reduced amounts of westerly-associated winter precipitation. Compared to the mid and late Holocene, it appears that westerlies had a greater influence on the western parts of the ISM domain during the early Holocene. During this period, the westerly-associated summer precipitation belt was positioned at Mediterranean latitudes and amplified the ISM-derived precipitation. The Tso Moriri pollen record and moisture reconstructions also suggest that changes in climatic conditions affected the ancient Harappan Civilisation, which flourished in the greater Indus Valley from approximately 5.2 to 3 cal ka BP. The prolonged Holocene trend towards aridity, punctuated by an interval of increased dryness (between ca. 4.5 to 4.3 cal ka BP), may have pushed the Mature Harappan urban settlements (between ca. 4.5 to 3.9 cal ka BP) to develop more efficient agricultural practices to deal with the increasingly acute water shortages. The amplified aridity associated with North Atlantic cooling between ca. 4 to 3.6 and around 3.2 cal ka BP further hindered local agriculture, possibly causing the deurbanisation that occurred from ca. 3.9 cal ka BP and eventual collapse of the Harappan Civilisation between ca. 3.5 to 3 cal ka BP.

Speleothem master chronologies: combined Holocene O-18 and C-13 records from the North Island of New Zealand and their palaeoenvironmental interpretation

The stable isotope records of four stalagmites dated by 19 TIMS uranium series ages are combined to produce master chronologies for delta(18)O and delta(13)C The delta(18)O records display good overall coherence, but considerable variation in detail. Variability in the delta(13)C records is greater, but general trends can still be discerned. This implies that too fine an interpretation of the structure of individual isotopic records can be unreliable. Speleothem delta(18)O values are demonstrated to show a positive relationship with temperature by comparing trends with other proxy records, but also to respond negatively to rainfall amount. Speleothem delta(13)C is considered to be most influenced by rainfall. The postglacial thermal optimum occur-red around 10.8 ka BP, which is similar in timing to Antarctica but up to 2000 years earlier than most Northern Hemisphere sites. Increasingly negative delta(18)O values after 7.5 ka BP indicate that temperatures declined to a late mid-Holocene minimum centred around 3 ka BP, but more positive values followed to mark a warm peak about 750 years ago which coincided with the 'Mediaeval Warm Period' of Europe. Low 5110 values at 325 years BP suggest cooling coincident with the 'Little Ice Age'. A marked feature of the delta(13)C record is an asymmetric periodicity averaging c. 2250 years and amplitude of c. 1.9parts per thousand. It is concluded that this is mainly driven by waterbalance variations with negative swings representing particularly wet intervals. The 5110 record shows a higher-frequency cyclicity with a period of c. 500 years and an amplitude of c. 0.25 parts per thousand. This is most likely to be temperature-driven, but some swings may have been amplified by precipitation.

Late pleistocene to holocene composite speleothem O-18 and C-13 chronologies from south island, new Zealand-did a global younger dryas really exist?

Oxygen and carbon data from eight stalagmites from northwest South Island are combined to produce composite records of delta(18)O and delta(13)C from 23.4 ka to the present. The chronology is anchored by 43 thermal ionization mass spectrometry (TIMS) uranium series ages. Delta O-18 values are interpreted as having a first order positive relationship to temperature, but also to be influenced by precipitation in a complex manner. Delta C-13 is interpreted as responding negatively to increases in atmospheric CO, concentration, biological activity and precipitation amount. Six climatic phases are recognized. After adjustment of 1.2parts per thousand for the ice volume effect, the delta(18)O record between 23 and 18 ka varies around -3.72parts per thousand compared to the Holocene average of -3.17parts per thousand. Late-glacial warming commenced between 18.2 and 17.8 ka and accelerated after 16.7 ka, culminating in a positive excursion between 14.70 and 13.53 ka. This was followed by a significant negative excursion between 13.53 and 11.14 ka of up to 0.55parts per thousand depth that overlapped the Antarctic Cold Reversal (ACR) and spanned the Younger Dryas (YD). Positive delta(18)O excursions at 11.14 ka and 6.91-6.47 ka represent the warmest parts of the Holocene. The mid-Holocene from 6 to 2 ka was marked by negative excursions that coincide with increased glacial activity in the South Island. A short positive excursion from 0.71 to 0.57 ka was slightly later than the Medieval Warm Period of Europe. Delta C-13 values were high until 17.79 ka after which there was an abrupt decrease to 17.19 ka followed by a steady decline to a minimum at 10.97 ka. Then followed a general increase, suggesting a drying trend, to 3.23 ka followed by a further general decline. The abrupt decrease in delta-values after 17.79 ka probably corresponds to an increase in atmospheric CO2 concentration, biological activity and wetness at the end of the Last Glaciation, but the reversal identified in the delta(18)O record from 13.53 to 11.14 ka was not reflected in delta(13)C changes. The lowest delta(13)C values coincided with the early Holocene climatic suboptimum when conditions were relatively wet as well as mild. Major trends in the delta(18)O(c) record are similar to the Northern Hemisphere, but second order detail is often distinctly different. Consequently, at the millennial scale, a more convincing case can be made for asymmetric climatic response between the two hemispheres rather than synchronicity. (C) 2004 Elsevier B.V. All rights reserved.

Age model and palaeoclimate records over the past 22,000 years of sediment core GeoB10053-7

The Australian-Indonesian monsoon is an important component of the climate system in the tropical Indo-Pacific region. However, its past variability, relation with northern and southern high-latitude climate and connection to the other Asian monsoon systems are poorly understood. Here we present high-resolution records of monsoon-controlled austral winter upwelling during the past 22,000 years, based on planktic foraminiferal oxygen isotopes and faunal composition in a sedimentary archive collected offshore southern Java. We show that glacial-interglacial variations in the Australian-Indonesian winter monsoon were in phase with the Indian summer monsoon system, consistent with their modern linkage through cross-equatorial surface winds. Likewise, millennial-scale variability of upwelling shares similar sign and timing with upwelling variability in the Arabian Sea. On the basis of element composition and grain-size distribution as precipitation-sensitive proxies in the same archive, we infer that (austral) summer monsoon rainfall was highest during the Bølling-Allerød period and the past 2,500 years. Our results indicate drier conditions during Heinrich Stadial 1 due to a southward shift of summer rainfall and a relatively weak Hadley cell south of the Equator. We suggest that the Australian-Indonesian summer and winter monsoon variability were closely linked to summer insolation and abrupt climate changes in the northern hemisphere.

Earliest diagenesis in scleractinian coral skeletons : implications for palaeoclimate-sensitive geochemical archives

Live-collected samples of four common reef building coral genera (Acropora, Pocillopora, Goniastrea, Porites) from subtidal and intertidal settings of Heron Reef, Great Barrier Reef, show extensive early marine diagenesis where parts of the coralla less than 3 years old contain abundant macro- and microborings and aragonite, high-Mg calcite, low-Mg calcite, and brucite cements. Many types of cement are associated directly with microendoliths and endobionts that inhabit parts of the corallum recently abandoned by coral polyps. The occurrence of cements that generally do not precipitate in normal shallow seawater (e.g., brucite, low-Mg calcite) highlights the importance of microenvironments in coral diagenesis. Cements precipitated in microenvironments may not reXect ambient seawater chemistry. Hence, geochemical sampling of these cements will contaminate trace-element and stable-isotope inventories used for palaeoclimate and dating analysis. Thus, great care must be taken in vetting samples for both bulk and microanalysis of geochemistry. Visual inspection using scanning electron microscopy may be required for vetting in many cases.