Tree-ring-reconstructed summer temperatures from northwestern North America during the last nine centuries

Northwestern North America has one of the highest rates of recent temperature increase in the world, but the putative “divergence problem” in dendroclimatology potentially limits the ability of tree-ring proxy data at high latitudes to provide long-term context for current anthropogenic change. Here, summer temperatures are reconstructed from a Picea glauca maximum latewood density (MXD) chronology that shows a stable relationship to regional temperatures and spans most of the last millennium at the Firth River in northeastern Alaska. The warmest epoch in the last nine centuries is estimated to have occurred during the late twentieth century, with average temperatures over the last 30 yr of the reconstruction developed for this study [1973–2002 in the Common Era (CE)] approximately 1.3° ± 0.4°C warmer than the long-term preindustrial mean (1100–1850 CE), a change associated with rapid increases in greenhouse gases. Prior to the late twentieth century, multidecadal temperature fluctuations covary broadly with changes in natural radiative forcing. The findings presented here emphasize that tree-ring proxies can provide reliable indicators of temperature variability even in a rapidly warming climate.

Tracing the origin of Arctic driftwood

Arctic environments, where surface temperatures increase and sea ice cover and permafrost depth decrease, are very sensitive to even slight climatic variations. Placing recent environmental change of the high-northern latitudes in a long-term context is, however, complicated by too short meteorological observations and too few proxy records. Driftwood may represent a unique cross-disciplinary archive at the interface of marine and terrestrial processes. Here, we introduce 1445 driftwood remains from coastal East Greenland and Svalbard. Macroscopy and microscopy were applied for wood anatomical classification; a multi-species subset was used for detecting fungi; and information on boreal vegetation patterns, circumpolar river systems, and ocean current dynamics was reviewed and evaluated. Four conifer (Pinus, Larix, Picea, and Abies) and three deciduous (Populus, Salix, and Betula) genera were differentiated. Species-specific identification also separated Pinus sylvestris and Pinus sibirica, which account for ~40% of all driftwood and predominantly originate from western and central Siberia. Larch and spruce from Siberia or North America represents ~26% and ~18% of all materials, respectively. Fungal colonization caused different levels of driftwood staining and/or decay. Our results demonstrate the importance of combining wood anatomical knowledge with insight on boreal forest composition for successfully tracing the origin of Arctic driftwood. To ultimately reconstruct spatiotemporal variations in ocean currents, and to better quantify postglacial uplift rates, we recommend consideration of dendrochronologically dated material from many more circumpolar sites.

Northern Hemisphere temperature reconstruction during the last millennium using multiple annual proxies

Previous studies have either exclusively used annual tree-ring data or have combined tree-ring series with other, lower temporal resolution proxy series. Both approaches can lead to significant uncertainties, as tree-rings may underestimate the amplitude of past temperature variations, and the validity of non-annual records cannot be clearly assessed. In this study, we assembled 45 published Northern Hemisphere (NH) temperature proxy records covering the past millennium, each of which satisfied 3 essential criteria: the series must be of annual resolution, span at least a thousand years, and represent an explicit temperature signal. Suitable climate archives included ice cores, varved lake sediments, tree-rings and speleothems. We reconstructed the average annual land temperature series for the NH over the last millennium by applying 3 different reconstruction techniques: (1) principal components (PC) plus second-order autoregressive model (AR2), (2) composite plus scale (CPS) and (3) regularized errors-in-variables approach (EIV). Our reconstruction is in excellent agreement with 6 climate model simulations (including the first 5 models derived from the fifth phase of the Coupled Model Intercomparison Project (CMIP5) and an earth system model of intermediate complexity (LOVECLIM), showing similar temperatures at multi-decadal timescales; however, all simulations appear to underestimate the temperature during the Medieval Warm Period (MWP). A comparison with other NH reconstructions shows that our results are consistent with earlier studies. These results indicate that well-validated annual proxy series should be used to minimize proxy-based artifacts, and that these proxy series contain sufficient information to reconstruct the low-frequency climate variability over the past millennium.

A case-control study of bilateral compared with unilateral breast cancer

Approximately 10 to 15% of breast cancer patients develop a primary cancer in the contralateral breast. This study examined differences between women with unilateral compared with bilateral primary breast cancer. It focused on hormonal factors and family history, and evaluated the prevalences of invasive lobular histology and the replication error phenotype in the tumors. ^ Cases (n = 82) were patients at M.D. Anderson Cancer Center (MDACC) in Houston, Texas diagnosed with primary breast cancer in each breast between 1985 and 1994 inclusive. Controls (n = 82) were MDACC patients with primary cancer in a single breast diagnosed during the same interval, individually matched to cases. Data were obtained by in-person and/or telephone interview with the patient and/or proxy. Replication error phenotype was determined from archival tissue. ^ Diagnosis of breast, but not ovarian, cancer in a female first-degree relative (FFDR) was a strong risk factor for bilateral cancers. Cases had a significantly 3-fold higher excess of familial breast cancer than did controls (cases: O/E = 2.65, 95% CI = 1.85–3.69; controls: 0.86, 0.46–1.47; homogeneity: p = 0.00). Risk did not vary with menopausal status of the patient, but was greatest if a relative was diagnosed before age 45 (O/E = 38.9; 95% CI = 21.7–64.1). By implication, young first-degree relatives of patients with bilateral breast cancer are at very high risk of breast cancer themselves. Cases also had significantly fewer siblings than did controls. ^ Earlier menarche, and parity in the absence of lactation, were associated with bilateral cancers; age at menopause and parity with lactation were not. A history of alcohol consumption, particularly if heavy, carried a 3.4-fold risk (p = 0.03). The data suggested a slightly different pattern in risk factors according to menopausal status and interval between cancers. ^ Replication error phenotype was available for 59 probands. It was associated with bilateral cancers (particularly if diagnosed within one year of each other), increased age (p = 0.02) and negative nodal status. Invasive lobular histology was associated with bilateral disease but numbers were small. ^ These data suggest bilateral breast cancer arises in the context of a combination of familial and hormonal factors, and alcohol consumption. The relative importance of each factor may vary by age of the patient. ^

Chronostratigraphy of the 600,000 year old continental record of Lake Van

Lake Van sediment cores from the Ahlat Ridge and Northern Basin drill sites of the ICDP project PALEOVAN contain a wealth of information about past environmental processes. The sedimentary sequence was dated using climatostratigraphic alignment, varve chronology, tephrostratigraphy, argon-argon single-crystal dating, radiocarbon dating, magnetostratigraphy, and cosmogenic nuclides. Based on the lithostratigraphic framework, the different age constraints are compiled and a robust and precise chronology of the 600,000 year-old Lake Van record is constructed. Proxy records of total organic carbon content and sediment color, together with the calcium/potassium-ratios and arboreal pollen percentages of the 174-meter-long Ahlat Ridge record, mimic the Greenland isotope stratotype (NGRIP). Therefore, the proxy records are systematically aligned to the onsets of interstadials reflected in the NGRIP or synthesized Greenland ice-core stratigraphy. The chronology is constructed using 27 age control points derived from visual synchronization with the GICC05 timescale, an absolutely-dated speleothem record (e.g., Hulu, Sanbao, Linzhu cave) and the Epica Dome C timescale. In addition, the uppermost part of the sequence is complemented with four ages from Holocene varve chronology and two calibrated radiocarbon ages. Furthermore, nine argon-argon ages and a comparison of the relative paleointensity record of the magnetic field with reference curve PISO-1500 confirm the accuracy of the age model. Also the identification of the Laschamp event via measurements of 10Be in the sediment confirms the presented age model. The chronology of the Ahlat Ridge record is transferred to the 79-meter-long event-corrected composite record from the Northern Basin and supplemented by additional radiocarbon dating on organic marco-remains. The basal age of the Northern Basin record is estimated at ~90 ka. The variations of the time series of total organic carbon content, the Ca/K ratio, and the arboreal pollen percentages illustrate that the presented chronology and paleoclimate data are suited for reconstructions and modeling of the Quaternary and Pleistocene climate evolution in the Near East at millennial timescales. Furthermore, the chronology of the last 250 kyr can be used to test other dating techniques.

The potential of stomata analysis in conifers to estimate presence of conifer trees: examples from the Alps

To estimate whether or not a plant taxon found in the fossil record was locally present may be difficult if only pollen is analyzed. Plant macrofossils, in contrast, provide a clear indication of a taxon's local presence, although in some lake sediments or peats, macrofossils may be rare or degraded. For conifers, the stomata found on pollen slides are derived from needles and thus provide a valuable proxy for local presence and they can be identified to genus level. From previously published studies, a transect across the Alps based on 13 sites is presented. For basal samples in sandy silt above the till with high pollen values of Pinus, for example, we may distinguish pine pollen from distant sources (samples with no stomata), from reworked pollen (samples with stomata present). The first apparent local presence of most conifer genera based on stomata often but not always occurs together with the phase of rapid pollen increase (rational limit). An exception is Larix, with its annual deposition of needles and heavy poorly dispersed pollen, for it often shows the first stomata earlier, at the empirical pollen limit. The decline and potential local extinction of a conifer can sometimes be shown in the stomata record. The decline may have been caused by climatic change, competition, or human impact. In situations where conifers form the timberline, the stomata record may indicate timberline fluctuations. In the discussion of immigration or migration of taxa we advocate the use of the cautious term ``apparent local presence'' to include some uncertainties. Absence of a taxon is impossible to prove.

A 2600-year record of past polycyclic aromatic hydrocarbons (PAHs) deposition at Holzmaar (Eifel, Germany)

Polycyclic aromatic hydrocarbons (PAHs) are a proxy for climate- and human-related historical fire activity which has rarely been used beyond 1800 AD. We explored the concentration and composition patterns of PAHs together with other proxies (charcoal, C, N, S, δ13C, δ15N, and δ34S) in a sediment core of Holzmaar as indicators of variations in climate and anthropogenic activity over the past 2600 years. The concentrations of pyrogenic PAHs remained low (< 500 ng g− 1) from the pre-Roman Iron Age (600 BC) until the first significant increases to ca. 1000–1800 ng g− 1 between 1700 and 1750 AD related to regional iron production. The highest increases in pyrogenic PAH concentrations occurred with industrialization peaking in the 1960s. PAH concentrations in most recent sediments decreased to pre-industrial levels because of emission control measures and the switch from coal to oil and gas as major fuel sources. Fluxes of PAHs (mg km− 2 yr− 1) increased in the order Migration Period and Early Middle Ages < pre-Roman Iron Ages < Roman era < High Middle Ages and Renaissance < pre-industrial modern period < industrial modern period. The changes of PAHs fluxes in pre-industrial times parallel known changes in local, regional and continental anthropogenic activity and climatic variations or their interactions across these historical periods. Perylene, a mainly biologically produced compound, was the dominant PAH in pre-industrial times. The Migration Period and the Early Middle Ages witnessed the most profound and simultaneous changes to the sedimentary organic matter properties in the past 2600 years including the lowest PAH fluxes because of reduced human activity and more negative δ13C and δ15N values probably in response to colder and wetter conditions.

Inter-hemispheric temperature variability over the past millennium

The Earth’s climate system is driven by a complex interplay of internal chaotic dynamics and natural and anthropogenic external forcing. Recent instrumental data have shown a remarkable degree of asynchronicity between Northern Hemisphere and Southern Hemisphere temperature fluctuations, thereby questioning the relative importance of internal versus external drivers of past as well as future climate variability1, 2, 3. However, large-scale temperature reconstructions for the past millennium have focused on the Northern Hemisphere4, 5, limiting empirical assessments of inter-hemispheric variability on multi-decadal to centennial timescales. Here, we introduce a new millennial ensemble reconstruction of annually resolved temperature variations for the Southern Hemisphere based on an unprecedented network of terrestrial and oceanic palaeoclimate proxy records. In conjunction with an independent Northern Hemisphere temperature reconstruction ensemble5, this record reveals an extended cold period (1594–1677) in both hemispheres but no globally coherent warm phase during the pre-industrial (1000–1850) era. The current (post-1974) warm phase is the only period of the past millennium where both hemispheres are likely to have experienced contemporaneous warm extremes. Our analysis of inter-hemispheric temperature variability in an ensemble of climate model simulations for the past millennium suggests that models tend to overemphasize Northern Hemisphere–Southern Hemisphere synchronicity by underestimating the role of internal ocean–atmosphere dynamics, particularly in the ocean-dominated Southern Hemisphere. Our results imply that climate system predictability on decadal to century timescales may be lower than expected based on assessments of external climate forcing and Northern Hemisphere temperature variations5, 6 alone.

Vegetation responses to rapid warming and to minor climatic fluctuations during the Late-Glacial Interstadial (GI-1) at Gerzensee (Switzerland)

High-resolution pollen analyses made on the same samples on which the ratios of oxygen isotopes were measured that provided the time scale and a temperature proxy after correlation to NorthGRIP. (1) A primary succession: The vegetation responded to the rapid rise of temperatures around 14,685 yr BP, with a primary succession on a decadal to centennial time scale. The succession between ca 15,600 and 13,000 yr BP included: (1.1.) The replacement of shrub-tundra by woodland of Juniperus and tree birch (around 14,665 yr BP) (1.2.) The response of Juniperus pollen to the shift in oxygen isotopes in less than 20 yr, (1.3.) A sequence of population increases of Hippophaë rhamnoides (ca 14,600 yr BP), Salix spp. (ca 14,600 yr BP), Betula trees (ca.14,480 yr BP), Populus cf. tremula (ca. 14,300 yr BP), and Pinus cf. sylvestris (ca. 13,830 yr BP). (2) Biological processes: Plants responded to the rapid increase of summer temperatures on all organisational levels: (2.1) Individuals may have produced more pollen (e.g. Juniperus); (2.2) Populations increased or decreased (e.g. Juniperus, Betula, later Pinus), and (2.3) Populations changed their biogeographical range and may show migrational lags. (2.4) Plant communities changed in their composition because the species pools changed through immigration and (local) extinction. Some plant communities may have been without modern analogue.These mechanisms require increasing amounts of time. (2.5) Processes on the level of ecosystems, with species interactions, may involve various time scales. Besides competition and facilitation, nitrogen fixation is discussed. (3) The minor fluctuations of temperature during the Late-Glacial Interstadial, which are recorded in δ18O, resulted in only very minor changes in pollen during the Aegelsee Oscillation (Older Dryas biozone, GI-1d) and the Gerzensee Oscillation (GI-1b). (4) Biodiversity: The afforestation at the onset of Bølling coincided with a gradual increase of taxonomic diversity up to the time of the major Pinus expansion.

Solar Forcing of the Polar Atmosphere

We present highly resolved, annually dated, calibrated proxies for atmospheric circulation from several Antarctic ice cores (ITASE (International Trans-Antarctic Scientific Expedition), Siple Dome, Law Dome) that reveal decadal-scale associations with a South Pole ice-core Be-10 proxy for solar variability over the last 600 years and annual-scale associations with solar variability since AD 1720. We show that increased (decreased) solar irradiance is associated with increased (decreased) zonal wind strength near the edge of the Antarctic polar vortex. The association is particularly strong in the Indian and Pacific Oceans and as such may contribute to understanding climate forcing that controls drought in Australia and other Southern Hemisphere climate events. We also include evidence suggestive of solar forcing of atmospheric circulation near the edge of the Arctic polar vortex based on ice-core records from Mount Logan, Yukon Territory, Canada, and both central and south Greenland as enticement for future investigations. Our identification of solar forcing of the polar atmosphere and its impact on lower latitudes offers a mechanism for better understanding modern climate variability and potentially the initiation of abrupt climate-change events that operate on decadal and faster scales.

Experimental Determination of Salinity, Temperature, Growth, and Metabolic Effects on Shell Isotope Chemistry of Mytilus edulis Collected from Maine and Greenland

To study the effects of temperature, salinity, and life processes (growth rates, size, metabolic effects, and physiological/ genetic effects) on newly precipitated bivalve carbonate, we quantified shell isotopic chemistry of adult and juvenile animals of the intertidal bivalve Mytilus edulis (Blue mussel) collected alive from western Greenland and the central Gulf of Maine and cultured them under controlled conditions. Data for juvenile and adult M. edulis bivalves cultured in this study, and previously by Wanamaker et al. (2006), yielded statistically identical paleotemperature relationships. On the basis of these experiments we have developed a species-specific paleotemperature equation for the bivalve M. edulis [T degrees C = 16.28 (+/- 0.10) -4.57 (+/- 0.15) {delta(18)O(c) VPBD - delta(18)O(w) VSMOW} + 0.06 (+/- 0.06) {delta(18)O(c) VPBD - delta(18)O(w) VSMOW}(2); r(2) = 0.99; N = 323; p < 0.0001]. Compared to the Kim and O'Neil (1997) inorganic calcite equation, M. edulis deposits its shell in isotope equilibrium (delta(18)O(calcite)) with ambient water. Carbon isotopes (delta(13)C(calcite)) from sampled shells were substantially more negative than predicted values, indicating an uptake of metabolic carbon into shell carbonate, and delta(13)C(calcite) disequilibrium increased with increasing salinity. Sampled shells of M. edulis showed no significant trends in delta(18)O(calcite) based on size, cultured growth rates, or geographic collection location, suggesting that vital effects do not affect delta(18)O(calcite) in M. edulis. The broad modern and paleogeographic distribution of this bivalve, its abundance during the Holocene, and the lack of an intraspecies physiologic isotope effect demonstrated here make it an ideal nearshore paleoceanographic proxy throughout much of the North Atlantic Ocean.

Atmospheric Soluble Dust Records from a Tibetan Ice Core: Possible Climate Proxies and Teleconnection With the Pacific Decadal Oscillation

In autumn 2005, a joint expedition between the University of Maine and the Institute of Tibetan Plateau Research recovered three ice cores from Guoqu Glacier (33 degrees 34'37.80 '' N, 91 degrees 10'35.3 '' E, 5720 m above sea level) on the northern side of Mt. Geladaindong, central Tibetan Plateau. Isotopes ( delta(18)O), major soluble ions (Na(+), K(+), Mg(2+), Ca(2+), Cl(-), NO(3)(-), SO(4)(2-)), and radionuclide (beta-activity) measurements from one of the cores revealed a 70-year record (1935-2005). Statistical analysis of major ion time series suggests that atmospheric soluble dust species dominate the chemical signature and that background dust levels conceal marine ion species deposition. The soluble dust time series have interspecies relations and common structure (empirical orthogonal function (EOF) 1), suggesting a similar soluble dust source or transport route. Annual and seasonal correlations between the EOF 1 time series and National Centers for Environmental Prediction/National Center for Atmospheric Research reanalysis climate variables (1948-2004) suggest that the Mt. Geladaindong ice core record provides a proxy for local and regional surface pressure. An approximately threefold decrease of soluble dust concentrations in the middle to late 1970s, accompanied by regional increases in pressure and temperature and decreases in wind velocity, coincides with the major 1976-1977 shift of the Pacific Decadal Oscillation (PDO) from a negative to a positive state. This is the first ice core evidence of a potential teleconnection between central Asian atmospheric soluble dust loading and the PDO. Analysis of temporally longer ice cores from Mt. Geladaindong may enhance understanding of the relationship between the PDO and central Asian atmospheric circulation and subsequent atmospheric soluble dust loading.

State of the Antarctic and Southern Ocean Climate System

This paper reviews developments in our understanding of the state of the Antarctic and Southern Ocean climate and its relation to the global climate system over the last few millennia. Climate over this and earlier periods has not been stable, as evidenced by the occurrence of abrupt changes in atmospheric circulation and temperature recorded in Antarctic ice core proxies for past climate. Two of the most prominent abrupt climate change events are characterized by intensification of the circumpolar westerlies (also known as the Southern Annular Mode) between similar to 6000 and 5000 years ago and since 1200-1000 years ago. Following the last of these is a period of major trans-Antarctic reorganization of atmospheric circulation and temperature between A. D. 1700 and 1850. The two earlier Antarctic abrupt climate change events appear linked to but predate by several centuries even more abrupt climate change in the North Atlantic, and the end of the more recent event is coincident with reorganization of atmospheric circulation in the North Pacific. Improved understanding of such events and of the associations between abrupt climate change events recorded in both hemispheres is critical to predicting the impact and timing of future abrupt climate change events potentially forced by anthropogenic changes in greenhouse gases and aerosols. Special attention is given to the climate of the past 200 years, which was recorded by a network of recently available shallow firn cores, and to that of the past 50 years, which was monitored by the continuous instrumental record. Significant regional climate changes have taken place in the Antarctic during the past 50 years. Atmospheric temperatures have increased markedly over the Antarctic Peninsula, linked to nearby ocean warming and intensification of the circumpolar westerlies. Glaciers are retreating on the peninsula, in Patagonia, on the sub-Antarctic islands, and in West Antarctica adjacent to the peninsula. The penetration of marine air masses has become more pronounced over parts of West Antarctica. Above the surface, the Antarctic troposphere has warmed during winter while the stratosphere has cooled year-round. The upper kilometer of the circumpolar Southern Ocean has warmed, Antarctic Bottom Water across a wide sector off East Antarctica has freshened, and the densest bottom water in the Weddell Sea has warmed. In contrast to these regional climate changes, over most of Antarctica, near-surface temperature and snowfall have not increased significantly during at least the past 50 years, and proxy data suggest that the atmospheric circulation over the interior has remained in a similar state for at least the past 200 years. Furthermore, the total sea ice cover around Antarctica has exhibited no significant overall change since reliable satellite monitoring began in the late 1970s, despite large but compensating regional changes. The inhomogeneity of Antarctic climate in space and time implies that recent Antarctic climate changes are due on the one hand to a combination of strong multidecadal variability and anthropogenic effects and, as demonstrated by the paleoclimate record, on the other hand to multidecadal to millennial scale and longer natural variability forced through changes in orbital insolation, greenhouse gases, solar variability, ice dynamics, and aerosols. Model projections suggest that over the 21st century the Antarctic interior will warm by 3.4 degrees +/- 1 degrees C, and sea ice extent will decrease by similar to 30%. Ice sheet models are not yet adequate enough to answer pressing questins about the effect of projected warming on mass balance and sea level. Considering the potentially major impacts of a warming climate on Antarctica, vigorous efforts are needed to better understand all aspects of the highly coupled Antarctic climate system as well as its influence on the Earth's climate and oceans.

LGM Summer Climate on the Southern Margin of the Laurentide Ice Sheet: Wet Or Dry?

Regional climate simulations are conducted using the Polar fifth-generation Pennsylvania State University (PSU)-NCAR Mesoscale Model (MM5) with a 60-km horizontal resolution domain over North America to explore the summer climate of the Last Glacial Maximum (LGM: 21 000 calendar years ago), when much of the continent was covered by the Laurentide Ice Sheet (LIS). Output from a tailored NCAR Community Climate Model version 3 (CCM3) simulation of the LGM climate is used to provide the initial and lateral boundary conditions for Polar MM5. LGM boundary conditions include continental ice sheets, appropriate orbital forcing, reduced CO2 concentration, paleovegetation, modified sea surface temperatures, and lowered sea level. The simulated LGM summer climate is characterized by a pronounced low-level thermal gradient along the southern margin of the LIS resulting from the juxtaposition of the cold ice sheet and adjacent warm ice-free land surface. This sharp thermal gradient anchors the midtropospheric jet stream and facilitates the development of synoptic cyclones that track over the ice sheet, some of which produce copious liquid precipitation along and south of the LIS terminus. Precipitation on the southern margin is orographically enhanced as moist southerly low-level flow (resembling a contemporary, Great Plains low-level jet configuration) in advance of the cyclone is drawn up the ice sheet slope. Composites of wet and dry periods on the LIS southern margin illustrate two distinctly different atmospheric flow regimes. Given the episodic nature of the summer rain events, it may be possible to reconcile the model depiction of wet conditions on the LIS southern margin during the LGM summer with the widely accepted interpretation of aridity across the Great Plains based on geological proxy evidence.

Ice Layers as an Indicator of Summer Warmth and Atmospheric Blocking in Alaska

Samples were collected from a snow pit and shallow urn core near Kahiltna Pass (2970 m a.s.l.), Denali National Park, Alaska, USA, in May 2008. The record spans autumn 2003 to spring 2008 and reveals clusters of ice layers interpreted as summertime intervals of above-freezing temperatures. High correlation coefficients (0.75-1.00) between annual ice-layer thickness and regional summertime station temperatures for 4 years (n=4) indicate ice-layer thickness is a good proxy for mean and extreme summertime temperatures across Alaska, at least over the short period of record. A Rex-block (aka high-over-low) pattern, a downstream trough over Hudson Bay, Canada, and an upstream trough over eastern Siberia occurred during the three melting events that lasted at least 2 weeks. About half of all shorter melting events were associated with a cut-off low traversing the Gulf of Alaska. We hypothesize that a surface-to-bedrock core extracted from this location would provide a high-quality record of summer temperature and atmospheric blocking variability for the last several hundred years.