Search

The Pliocene Yorktown Formation consists of four lithologic units that record three marine transgressive sequences along the U.S. mid-Atlantic margin. These units were deposited during a time when average sea-level and mean global temperatures were ∼25 m and ∼3°C higher than the pre-industrial, respectively, and global atmospheric CO2 concentrations were similar to present. Forty-five samples were collected along the James River near Rushmere, Virginia, and Spring Grove, Virginia, and were analyzed for benthic foraminifera community and sedimentological changes between each member of the formation. These data are useful for developing boundary conditions for shallow, near-shore environments for paleoclimate modeling. Foraminiferal analysis distinguishes six biofacies across the entire formation. The most notable change in grain-size occurs at the conformable boundary between the Rushmere and Morgarts Beach members, where the average percentage of sand decreases from ∼60% in the Rushmere Member to

The presence/absence and abundance of benthic foraminifera in successive discrete beds (Shattuck “zones”) of the Miocene Calvert and Choptank formations, exposed at the Calvert Cliffs, Maryland, USA, allows for investigation of community dynamics over space and time. The stratigraphic distribution of benthic foraminifera is documented and interpreted in the context of sea-level change, sequence stratigraphy, and the previously published distribution of mollusks. Neritic benthic foraminiferal communities of four sea-level cycles over ∼4 million years of the middle Miocene, encompassing the Miocene Climatic Optimum and the succeeding middle Miocene Climate Transition, are dominated by the same abundant species. They differ in the varying abundance of common species that occur throughout most of the studied section and in the different rare species that appear and disappear. Transgressive systems tracts (TSTs) have higher species diversity than highstand systems tracts (HSTs) but much lower density of specimens. In contrast to some previous research, all beds in the studied section are interpreted as being from the inner part of a broad, low gradient shelf and were deposited at water depths of less than ∼50 m. It is suggested that species are recruited from a regional species pool of propagules throughout the duration of TSTs. Recruitment is curtailed during highstands leading to lower diversity in the HSTs.

The Yorktown Formation records paleoclimate conditions along the mid-Atlantic Coastal Plain during the mid-Piacenzian Warm Period (3.264 to 3.025 Ma), a climate interval of the Pliocene in some ways analogous to near future climate projections. To gain insight into potential near future changes, we investigated Yorktown Formation outcrops and cores in southeastern Virginia, refining the stratigraphic framework. We analyzed 485 samples for alkenone-based sea surface temperature (SST) and productivity estimates from the Holland and Dory cores, an outcrop at Morgarts Beach, Virginia, and the lectostratotype of the Yorktown Formation at Rushmere, Virginia, and analyzed planktonic foraminferal assemblage data from the type section. Using the structure of the SST record, we improved the chronology of the Yorktown Formation by establishing the maximum age ranges of the Rushmere (3.3–3.2 Ma) and Morgarts Beach (3.2–3.15 Ma) Members. SST values for these members average ~26 °C, corroborating existing sclerochronological data. Increasing planktonic foraminifer abundance, productivity, and species diversity parallel increasing SST over the MIS M2/M1 transition. These records constitute the greatest temporal concentration of paleoecological estimates within the Yorktown Formation, aiding our understanding of western North Atlantic temperature patterns, seasonality and ocean circulation during this interval. We provide a chronologic framework for future studies analyzing ecological responses to profound climate change.

“Speak to the past and it shall teach thee”. I first read those words on a dedication tablet within the John Carter Brown library at Brown University where I was a graduate student. Little did I know the phrase would accurately describe the next three and a half decades of my career. Paleoclimate data are the language we use to look into the past to understand ourselves and ultimately our future.

Factorisation is widely used in the analysis of numerical simulations. It allows changes in properties of a system to be attributed to changes in multiple variables associated with that system. There are many possible factorisation methods; here we discuss three previously-proposed factorisations that have been applied in the field of climate modelling: the linear factorisation, the Stein and Alpert (1993) factorisation, and the Lunt et al (2012) factorisation. We show that, when more than two variables are being considered, none of these three methods possess all three properties of uniqueness, symmetry, and completeness. Here, we extend each of these factorisations so that they do possess these properties for any number of variables, resulting in three factorisations – the linear-sum factorisation, the shared-interaction factorisation, and the scaled-total factorisation. We show that the linear-sum factorisation and the shared-interaction factorisation reduce to be identical. We present the results of the factorisations in the context of studies that used the previously-proposed factorisations. This reveals that only the linear-sum/shared-interaction factorisation possesses a fourth property – boundedness, and as such we recommend the use of this factorisation in applications for which these properties are desirable.

The Pliocene epoch has great potential to improve our understanding of the long-term climatic and environmental consequences of an atmospheric CO2 concentration near ~ 400 parts per million by volume. Here we present the large-scale features of Pliocene climate as simulated by a new ensemble of climate models of varying complexity and spatial resolution and based on new reconstructions of boundary conditions (the Pliocene Model Intercomparison Project Phase 2; PlioMIP2). As a global annual average, modelled surface air temperatures increase by between 1.4 and 4.7 °C relative to pre-industrial with a multi-model mean value of 2.8 °C. Annual mean total precipitation rates increase by 6 % (range: 2 %–13 %). On average, surface air temperature (SAT) increases are 1.3 °C greater over the land than over the oceans, and there is a clear pattern of polar amplification with warming polewards of 60° N and 60° S exceeding the global mean warming by a factor of 2.4. In the Atlantic and Pacific Oceans, meridional temperature gradients are reduced, while tropical zonal gradients remain largely unchanged. Although there are some modelling constraints, there is a statistically significant relationship between a model's climate response associated with a doubling in CO2 (Equilibrium Climate Sensitivity; ECS) and its simulated Pliocene surface temperature response. The mean ensemble earth system response to doubling of CO2 (including ice sheet feedbacks) is approximately 50 % greater than ECS, consistent with results from the PlioMIP1 ensemble. Proxy-derived estimates of Pliocene sea-surface temperatures are used to assess model estimates of ECS and indicate a range in ECS from 2.5 to 4.3 °C. This result is in general accord with the range in ECS presented by previous IPCC Assessment Reports.

The mid-Piacenzian (Pliocene) climate represents the most geologically recent interval of long-term average warmth, relative to the last million years, sharing similarities with the climate project...

The Arctic cryosphere is changing and making a significant contribution to sea level rise. The Late Pliocene had similar CO2 levels to the present and a warming comparable to model predictions for the end of this century. However, the state of the Arctic cryosphere during the Pliocene remains poorly constrained. For the first time we combine outputs from a climate model with a thermodynamic iceberg model to simulate likely source regions for ice-rafted debris (IRD) found in the Nordic Seas from Marine Isotope Stage M2 to the mid-Piacenzian Warm Period and what this implies about the nature of the Arctic cryosphere at this time. We compare the fraction of melt given by the model scenarios with IRD data from four Ocean Drilling Program sites in the Nordic Seas. Sites 911A, 909C, and 907A show a persistent occurrence of IRD that model results suggest is consistent with permanent ice on Svalbard. Our results indicate that icebergs sourced from the east coast of Greenland do not reach the Nordic Seas sites during the warm Late Pliocene but instead travel south into the North Atlantic. In conclusion, we suggest a continuous occurrence of marine-terminating glaciers on Svalbard and on East Greenland (due to the elevation of the East Greenland Mountains during the Late Pliocene). The study has highlighted the usefulness of coupled climate model-iceberg trajectory modeling for understanding ice sheet behavior when proximal geological records for Pliocene ice presence or absence are absent or are inconclusive.

Ostracoda (bivalved Crustacea) comprise a significant part of the benthic meiofauna in the Pacific-Arctic region, including more than 50 species, many with identifiable ecological tolerances. These species hold potential as useful indicators of past and future ecosystem changes. In this study, we examined benthic ostracodes from nearly 300 surface sediment samples, >34,000 specimens, from three regions—the northern Bering, Chukchi and Beaufort Seas—to establish species’ ecology and distribution. Samples were collected during various sampling programs from 1970 through 2018 on the continental shelves at 20 to ~100m water depth. Ordination analyses using species’ relative frequencies identified six species,Normanicythere leioderma,Sarsicytheridea bradii,Paracyprideis pseudopunctillata,Semicytherura complanata,Schizocythere ikeyai, andMunseyella mananensis, as having diagnostic habitat ranges in bottom water temperatures, salinities, sediment substrates and/or food sources. Species relative abundances and distributions can be used to infer past bottom environmental conditions in sediment archives for paleo-reconstructions and to characterize potential changes in Pacific-Arctic ecosystems in future sampling studies. Statistical analyses further showed ostracode assemblages grouped by the summer water masses influencing the area. Offshore-to-nearshore transects of samples across different water masses showed that complex water mass characteristics, such as bottom temperature, productivity, as well as sediment texture, influenced the relative frequencies of ostracode species over small spatial scales. On the larger biogeographic scale, synoptic ordination analyses showed dominant species—N.leioderma(Bering Sea),P.pseudopunctillata(offshore Chukchi and Beaufort Seas), andS.bradii(all regions)—remained fairly constant over recent decades. However, during 2013–2018, northern Pacific speciesM.mananensisandS.ikeyaiincreased in abundance by small but significant proportions in the Chukchi Sea region compared to earlier years. It is yet unclear if these assemblage changes signify a meiofaunal response to changing water mass properties and if this trend will continue in the future. Our new ecological data on ostracode species and biogeography suggest these hypotheses can be tested with future benthic monitoring efforts.

International audience; Planktonic foraminiferal species identification is central to many paleoceanographic studies, from selecting species for geochemical research to elucidating the biotic dynamics of microfossil communities relevant to physical oceanographic processes and interconnected phenomena such as climate change. However, few resources exist to train students in the difficult task of discerning amongst closely related species, resulting in diverging taxonomic schools that differ in species concepts and boundaries. This problem is exacerbated by the limited number of taxonomic experts. Here we document our initial progress toward removing these confounding and/or rate-limiting factors by generating the first extensive image library of modern planktonic foraminifera, providing digital taxonomic training tools and resources, and automating species-level taxonomic identification of planktonic foraminifera via machine learning using convolution neural networks. Experts identified 34,640 images of modern (extant) planktonic foraminifera to the species level. These images are served as species exemplars through the online portal Endless Forams (endlessforams.org) and a taxonomic training portal hosted on the citizen science platform Zooniverse (zooniverse.org/projects/ahsiang/ endless-forams/). A supervised machine learning classifier was then trained with~27,000 images of these identified planktonic foraminifera. The best-performing model provided the correct species name for an image in the validation set 87.4% of the time and included the correct name in its top three guesses 97.7% of the time. Together, these resources provide a rigorous set of training tools in modern planktonic foraminiferal taxonomy and a means of rapidly generating assemblage data via machine learning in future studies for applications such as paleotemperature reconstruction.

The Piacenzian (Late Pliocene) represents a natural laboratory within which frequency and magnitude of environmental changes during a period of past global warmth can be analyzed, climate models ca...

The mid-Piacenzian is known as a period of relative warmth when compared to the present day. A comprehensive understanding of conditions during the Piacenzian serves as both a conceptual model and a source for boundary conditions as well as means of verification of global climate model experiments. In this paper we present the PRISM4 reconstruction, a paleoenvironmental reconstruction of the mid-Piacenzian ( ∼ 3 Ma) containing data for paleogeography, land and sea ice, sea-surface temperature, vegetation, soils, and lakes. Our retrodicted paleogeography takes into account glacial isostatic adjustments and changes in dynamic topography. Soils and lakes, both significant as land surface features, are introduced to the PRISM reconstruction for the first time. Sea-surface temperature and vegetation reconstructions are unchanged but now have confidence assessments. The PRISM4 reconstruction is being used as boundary condition data for the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2) experiments.

The mid-Pliocene Warm Period (mPWP), analogous to future climate conditions, is considered a test-bed for the predictive capability of climate models. Here, Dowsett et al. review our understanding of the mPWP and discuss recent and future advances in the context of proxy data/model integration., The mid-Pliocene Warm Period (mPWP) offers an opportunity to understand a warmer-than-present world and assess the predictive ability of numerical climate models. Environmental reconstruction and climate modelling are crucial for understanding the mPWP, and the synergy of these two, often disparate, fields has proven essential in confirming features of the past and in turn building confidence in projections of the future. The continual development of methodologies to better facilitate environmental synthesis and data/model comparison is essential, with recent work demonstrating that time-specific (time-slice) syntheses represent the next logical step in exploring climate change during the mPWP and realizing its potential as a test bed for understanding future climate change.

The global distribution of late Pliocene soils and lakes has been reconstructed using a synthesis of geological data. These reconstructions are then used as boundary conditions for the Hadley Centre General Circulation Model (HadCM3) and the BIOME4 mechanistic vegetation model. By combining our novel soil and lake reconstructions with a fully coupled climate model we are able to explore the feedbacks of soils and lakes on the climate of the late Pliocene. Our experiments reveal regionally confined changes of local climate and vegetation in response to the new boundary conditions. The addition of late Pliocene soils has the largest influence on surface air temperatures, with notable increases in Australia, the southern part of northern Africa and in Asia. The inclusion of late Pliocene lakes increases precipitation in central Africa and at the locations of lakes in the Northern Hemisphere. When combined, the feedbacks on climate from late Pliocene lakes and soils improve the data to model fit in western North America and the southern part of northern Africa.

The mid-Pliocene Warm Period (mPWP) bears many similarities to aspects of future global warming as projected by the Intergovernmental Panel on Climate Change (IPCC, 2007). Both marine and terrestrial data point to high-latitude temperature amplification, including large decreases in sea ice and land ice, as well as expansion of warmer climate biomes into higher latitudes. Here we present our most recent simulations of the mid-Pliocene climate using the CMIP5 version of the NASA/GISS Earth System Model (ModelE2-R). We describe the substantial impact associated with a recent correction made in the implementation of the Gent-McWilliams ocean mixing scheme (GM), which has a large effect on the simulation of ocean surface temperatures, particularly in the North Atlantic Ocean. The effect of this correction on the Pliocene climate results would not have been easily determined from examining its impact on the preindustrial runs alone, a useful demonstration of how the consequences of code improvements as seen in modern climate control runs do not necessarily portend the impacts in extreme climates. Both the GM-corrected and GM-uncorrected simulations were contributed to the Pliocene Model Intercomparison Project (PlioMIP) Experiment 2. Many findings presented here corroborate results from other PlioMIP multi-model ensemble papers, but we also emphasise features in the ModelE2-R simulations that are unlike the ensemble means. The corrected version yields results that more closely resemble the ocean core data as well as the PRISM3D reconstructions of the mid-Pliocene, especially the dramatic warming in the North Atlantic and Greenland-Iceland-Norwegian Sea, which in the new simulation appears to be far more realistic than previously found with older versions of the GISS model. Our belief is that continued development of key physical routines in the atmospheric model, along with higher resolution and recent corrections to mixing parameterisations in the ocean model, have led to an Earth System Model that will produce more accurate projections of future climate.

General circulation model (GCM) simulations of the mid-Pliocene Warm Period (mPWP, 3.264 to 3.025 Myr ago) do not reproduce the magnitude of Northern Hemisphere high latitude surface air and sea surface temperature (SAT and SST) warming that proxy data indicate. There is also large uncertainty regarding the state of sea ice cover in the mPWP. Evidence for both perennial and seasonal mPWP Arctic sea ice is found through analyses of marine sediments, whilst in a multi-model ensemble of mPWP climate simulations, half of the ensemble simulated ice-free summer Arctic conditions. Given the strong influence that sea ice exerts on high latitude temperatures, an understanding of the nature of mPWP Arctic sea ice would be highly beneficial. Using the HadCM3 GCM, this paper explores the impact of various combinations of potential mPWP orbital forcing, atmospheric CO 2 concentrations and minimum sea ice albedo on sea ice extent and high latitude warming. The focus is on the Northern Hemisphere, due to availability of proxy data, and the large data–model discrepancies in this region. Changes in orbital forcings are demonstrated to be sufficient to alter the Arctic sea ice simulated by HadCM3 from perennial to seasonal. However, this occurs only when atmospheric CO 2 concentrations exceed 300 ppm. Reduction of the minimum sea ice albedo from 0.5 to 0.2 is also sufficient to simulate seasonal sea ice, with any of the combinations of atmospheric CO 2 and orbital forcing. Compared to a mPWP control simulation, monthly mean increases north of 60°N of up to 4.2 °C (SST) and 9.8 °C (SAT) are simulated. With varying CO 2 , orbit and sea ice albedo values we are able to reproduce proxy temperature records that lean towards modest levels of high latitude warming, but other proxy data showing greater warming remain beyond the reach of our model. This highlights the importance of additional proxy records at high latitudes and ongoing efforts to compare proxy signals between sites.

The Pliocene Model Intercomparison Project (PlioMIP) is a co-ordinated international climate modelling initiative to study and understand climate and environments of the Late Pliocene, as well as their potential relevance in the context of future climate change. PlioMIP examines the consistency of model predictions in simulating Pliocene climate and their ability to reproduce climate signals preserved by geological climate archives. Here we provide a description of the aim and objectives of the next phase of the model intercomparison project (PlioMIP Phase 2), and we present the experimental design and boundary conditions that will be utilized for climate model experiments in Phase 2. Following on from PlioMIP Phase 1, Phase 2 will continue to be a mechanism for sampling structural uncertainty within climate models. However, Phase 1 demonstrated the requirement to better understand boundary condition uncertainties as well as uncertainty in the methodologies used for data–model comparison. Therefore, our strategy for Phase 2 is to utilize state-of-the-art boundary conditions that have emerged over the last 5 years. These include a new palaeogeographic reconstruction, detailing ocean bathymetry and land–ice surface topography. The ice surface topography is built upon the lessons learned from offline ice sheet modelling studies. Land surface cover has been enhanced by recent additions of Pliocene soils and lakes. Atmospheric reconstructions of palaeo-CO2 are emerging on orbital timescales, and these are also incorporated into PlioMIP Phase 2. New records of surface and sea surface temperature change are being produced that will be more temporally consistent with the boundary conditions and forcings used within models. Finally we have designed a suite of prioritized experiments that tackle issues surrounding the basic understanding of the Pliocene and its relevance in the context of future climate change in a discrete way.

High tropical and low polar biodiversity is one of the most fundamental patterns characterising marine ecosystems, and the influence of temperature on such marine latitudinal diversity gradients is increasingly well documented. However, the temporal stability of quantitative relationships among diversity, latitude and temperature is largely unknown. Herein we document marine zooplankton species diversity patterns at four time slices [modern, Last Glacial Maximum (18,000 years ago), last interglacial (120,000 years ago), and Pliocene (~3.3-3.0 million years ago)] and show that, although the diversity-latitude relationship has been dynamic, diversity-temperature relationships are remarkably constant over the past three million years. These results suggest that species diversity is rapidly reorganised as species' ranges respond to temperature change on ecological time scales, and that the ecological impact of future human-induced temperature change may be partly predictable from fossil and paleoclimatological records.

The mid-Pliocene (approximately 3 to 3.3 Ma ago), is a period of sustained global warmth in comparison to the late Quaternary (0 to approximately 1 Ma ago), and has potential to inform predictions of long-term future climate change. However, given that several processes potentially contributed, relatively little is understood about the reasons for the observed warmth, or the associated polar amplification. Here, using a modelling approach and a novel factorisation method, we assess the relative contributions to mid-Pliocene warmth from: elevated CO2, lowered orography, and vegetation and ice sheet changes. The results show that on a global scale, the largest contributor to mid-Pliocene warmth is elevated CO2. However, in terms of polar amplification, changes to ice sheets contribute significantly in the Southern Hemisphere, and orographic changes contribute significantly in the Northern Hemisphere. We also carry out an energy balance analysis which indicates that that on a global scale, surface albedo and atmospheric emmissivity changes dominate over cloud changes. We investigate the sensitivity of our results to uncertainties in the prescribed CO2 and orographic changes, to derive uncertainty ranges for the various contributing processes.

Anthropogenic aerosols over the Chinese Loess Plateau have diminished monsoon precipitation and concomitant soil erosion that plagues the region. Now, a reconstruction documents the differences between historical warming events and the present, highlighting the paradoxical implications of decreasing atmospheric aerosols.

The stability of the Earth's major ice sheets is a critical uncertainty in predictions of future climate and sea level change. One method of investigating the behaviour of the Greenland and the Antarctic ice sheets in a warmer-than-modern climate is to look back at past warm periods of Earth history, for example the Pliocene. This paper presents climate and ice sheet modelling results for the mid-Pliocene warm period (mPWP; 3.3 to 3.0 million years ago), which has been identified as a key interval for understanding warmer-than-modern climates (Jansen et al., 2007). Using boundary conditions supplied by the United States Geological Survey PRISM Group (Pliocene Research, Interpretation and Synoptic Mapping), the Hadley Centre coupled ocean–atmosphere climate model (HadCM3) and the British Antarctic Survey Ice Sheet Model (BASISM), we show large reductions in the Greenland and East Antarctic Ice Sheets (GrIS and EAIS) compared to modern in standard mPWP experiments. We also present the first results illustrating the variability of the ice sheets due to realistic orbital forcing during the mid-Pliocene. While GrIS volumes are lower than modern under even the most extreme (cold) mid-Pliocene orbit (losing at least 35% of its ice mass), the EAIS can both grow and shrink, losing up to 20% or gaining up to 10% of its present-day volume. The changes in ice sheet volume incurred by altering orbital forcing alone means that global sea level can vary by more than 25 m during the mid-Pliocene. However, we have also shown that the response of the ice sheets to mPWP orbital hemispheric forcing can be in anti-phase, whereby the greatest reductions in EAIS volume are concurrent with the smallest reductions of the GrIS. If this anti-phase relationship is in operation throughout the mPWP, then the total eustatic sea level response would be dampened compared to the ice sheet fluctuations that are theoretically possible. This suggests that maximum eustatic sea level rise does not correspond to orbital maxima, but occurs at times where the anti-phasing of Northern and Southern Hemisphere ice sheet retreat is minimised.

The mid-Pliocene warm period (MPWP; ~ 3.3 to 3.0 Ma) is the most recent interval in Earth's history in which global temperatures reached and remained at levels similar to those projected for the near future. The distribution of global warmth, however, was different than today in that the high latitudes warmed more than the tropics. Multiple temperature proxies indicate significant sea surface warming in the North Atlantic and Arctic Oceans during the MPWP, but predictions from a fully coupled ocean–atmosphere model (HadCM3) have so far been unable to fully predict the large scale of sea surface warming in the high latitudes. If climate proxies accurately represent Pliocene conditions, and if no weakness exists in the physics of the model, then model boundary conditions may be in error. Here we alter a single boundary condition (bathymetry) to examine if Pliocene high latitude warming was aided by an increase in poleward heat transport due to changes in the subsidence of North Atlantic Ocean ridges. We find an increase in both Arctic sea surface temperature and deepwater production in model experiments that incorporate a deepened Greenland–Scotland Ridge. These results offer both a mechanism for the warming in the North Atlantic and Arctic Oceans indicated by numerous proxies and an explanation for the apparent disparity between proxy data and model simulations of Pliocene northern North Atlantic and Arctic Ocean conditions. Determining the causes of Pliocene warmth remains critical to fully understanding comparisons of the Pliocene warm period to possible future climate change scenarios.

It is essential to document how well the current generation of climate models performs in simulating past climates to have confidence in their ability to project future conditions. We present the first global, in-depth comparison of Pliocene sea surface temperature (SST) estimates from a coupled ocean–atmosphere climate model experiment and a SST reconstruction based on proxy data. This enables the identification of areas in which both the climate model and the proxy dataset require improvement. In general, the fit between model-produced SST anomalies and those formed from the available data is very good. We focus our discussion on three regions where the data–model anomaly exceeds 2 °C. 1) In the high latitude North Pacific, a systematic model error may result in anomalies that are too cold. Also, the deeper Pliocene thermocline may cause disagreement along the California margin; either the upwelling in the model is too strong or the modeled thermocline is not deep enough. 2) In the North Atlantic, the model predicts cooling in the center of a data-based warming trend that steadily increases with latitude from + 1.5 °C to >+ 6 °C. The discrepancy may arise because the modeled North Atlantic Current is too zonal compared to reality, which is reinforced by the lowering of the altitude of the Pliocene Western Cordillera Mountains. In addition, the model's use of modern bathymetry in the higher latitudes may have led the model to underestimate the northward penetration of warmer surface water into the Arctic. 3) Finally, though the data and model show good general agreement across most of the Southern Ocean, a few locations show offsets due to the modern land–sea mask used in the model. Additional considerations could account for many of the modest data–model anomalies, such as differences between calibration climatologies, the oversimplification of the seasonal cycle, and differences between SST proxies (i.e. seasonality and water depth). New SST estimates from data-sparse and regionally important areas will greatly enhance our ability to judge model performance.