The Miocene (∼23–5 Ma) is a past warm epoch when global surface temperatures varied between ∼5 and 8°C warmer than today, and CO₂ concentration was ∼400–800 ppm. The narrowing/closing of the tropical ocean gateways and widening of high-latitude gateways throughout the Miocene is likely responsible for the evolution of the ocean's overturning circulation to its modern structure, though the mechanisms remain unclear. Here, we investigate early and middle Miocene ocean circulation in an opportunistic climate model intercomparison (MioMIP1), using 14 simulations with different paleogeography, CO₂, and vegetation. The strength of the Southern Ocean-driven Meridional Overturning Circulation (SOMOC) bottom cell is similar in the Miocene and Pre-Industrial (PI) but dominates the Miocene global MOC due to weaker Northern Hemisphere overturning. The Miocene Atlantic MOC (AMOC) is weaker than PI in all the simulations (by 2–21 Sv), possibly due to its connection with an Arctic that is considerably fresher than today. Deep overturning in the North Pacific (PMOC) is present in three simulations (∼5–10 Sv), of which two have a weaker AMOC, and one has a stronger AMOC (compared to its PMOC). Surface freshwater fluxes control northern overturning such that the basin with the least freshwater gain has stronger overturning. While the orography, which impacts runoff direction (Pacific vs. Atlantic), has an inconsistent impact on northern overturning across simulations, overall, features associated with the early Miocene—such as a lower Tibetan Plateau, the Rocky Mountains, and a deeper Panama Seaway—seem to favor PMOC over AMOC.

The large-scale oceanic gyre circulation in the Pacific regulates its temperature, salinity and nutrient flow, profoundly influencing the biological environment and climate. However, how gyres evolved during the warm Eocene, with its different coastline configurations, remains unknown. Here, we investigate the response of Pacific gyre circulation during the warm early Eocene using eight models from the Deep-Time Model Intercomparison Project (DeepMIP). Our DeepMIP results suggest a northward expansion of the North Pacific subtropical gyre by up to 10 degrees latitude during the Eocene, while maintaining a strength comparable to that of the present day. This simulated poleward expansion of the North Pacific gyre circulation is corroborated by sedimentary evidence, including poleward shifts in the pattern of clay sediments and low sedimentation rate during the Eocene. In the southern Pacific, the super subtropical gyre was much stronger during the Eocene than at present, due to the southward position of Australia, which created a wide-open Indonesian gateway. The poleward shift in the boundary between the subtropical and subpolar gyre in in the North Pacific is attributed to the northward migration of the westerly winds maxima, as confirmed by an analysis of Sverdrup transport. The Sverdrup-balanced upper circulation in the Pacific extends further poleward than in the modern day, mainly due to differences in the position of the continents. Specifically, the circulation corresponds to Sverdrup transport up to ∼53°N in the North Pacific, slightly further north than modern of limit of 50°N, and up to ∼55°S in the South Pacific, further south than in the modern limit of ∼45°S.

The Miocene (23.03-5.33 Ma) is recognized as a period with close to modern-day paleogeography, yet a much warmer climate. With large uncertainties in future hydroclimate projections, Miocene conditions illustrate a potential future analog for the Earth system. A recent opportunistic Miocene Model Intercomparison Project 1 (MioMIP1) focused on synthesizing published Miocene climate simulations and comparing them with available temperature reconstructions. Here, we build on this effort by analyzing the hydrological cycle response to Miocene forcings across early-to-middle (E2MMIO; 20.03-11.6 Ma) and middle-to-late Miocene (M2LMIO; 11.5-5.33 Ma) simulations with CO2 concentrations ranging from 200 to 850 ppm and providing a model-data comparison against available precipitation reconstructions. We find global precipitation increases by similar to 2.1 and 2.3% per degree of warming for E2MMIO and M2LMIO simulations, respectively. Models generally agree on a wetter than modern-day tropics; mid and high-latitude, however, do not agree on the sign of subtropical precipitation changes with warming. Global monsoon analysis suggests most monsoon regions, except the North American Monsoon, experience higher precipitation rates under warmer conditions. Model-data comparison shows that mean annual precipitation is underestimated by the models regardless of CO2 concentration, particularly in the mid- to high-latitudes. This suggests that the models may not be (a) resolving key processes driving the hydrological cycle response to Miocene boundary conditions and/or (b) other boundary conditions or processes not considered here are critical to reproducing Miocene hydroclimate. This study highlights the challenges in modeling and reconstructing the Miocene hydrological cycle and serves as a baseline for future coordinated MioMIP efforts. This study looks at Earth's hydrological cycle during the Miocene (23-5 million years ago). During this period, the Earth's climate was 3-7 degrees C warmer than today, with carbon dioxide (CO2) estimates ranging between 400 and 850 ppm. Understanding how the hydrological cycle responded during warmer climate conditions can give us insight into what might happen as the Earth gets warmer. We analyzed a suite of Miocene paleoclimate simulations with different CO2 concentrations in the atmosphere and compared them against fossil plant data, which gives an estimate of the average annual rainfall during the period. We found that during the Miocene global rainfall increased by about 2.1%-2.3% for each degree of warming. The models agree that the tropics, mid- and high-latitude, became wetter than they are today but have lower agreement on whether subtropical areas got wetter or drier as they warmed. Compared to proxies, models consistently underestimated how much rain fell in a year, especially in the mid- to high-latitude. This illustrates the challenges in reconstructing the Miocene's hydrological cycle and suggests that the models might not fully capture the range of uncertainties associated with changes in the hydrological cycle due to warming or other factors that differentiated the Miocene. A multi-model comparison of the hydrological cycle in early-to-middle and middle-to-late Miocene simulations is conductedModels generally agree on wetter than modern tropics, middle and high latitudes, but not on the sign of subtropical precipitation changesModel-data comparison shows mean annual precipitation is underestimated by the models, particularly in the mid- to high-latitudes

Paleoclimate model simulations provide reference data to help interpret the geological record and offer a unique opportunity to evaluate the performance of current models under diverse boundary conditions. Here, we present a dataset of 35 climate model simulations of the warm early Eocene Climatic Optimum (EECO; ~ 50 million years ago) and corresponding preindustrial reference experiments. To streamline the use of the data, we apply standardised naming conventions and quality checks across eight modelling groups that have carried out coordinated simulations as part of the Deep-Time Model Intercomparison Project (DeepMIP). Gridded model fields can be downloaded from an online repository or accessed through a new web application that provides interactive data exploration. Local model data can be extracted in CSV format or visualised online for streamlined model-data comparisons. Additionally, processing and visualisation code templates may serve as a starting point for advanced analysis. The dataset and online platform aim to simplify accessing and handling complex data, prevent common processing issues, and facilitate the sharing of climate model data across disciplines.

The extent and seasonality of Arctic sea ice during the Last Interglacial (129,000 to 115,000 years before present) is poorly known. Sediment-based reconstructions have suggested extensive ice cover in summer, while climate model outputs indicate year-round conditions in the Arctic Ocean ranging from ice free to fully ice covered. Here we use microfossil records from across the central Arctic Ocean to show that sea-ice extent was substantially reduced and summers were probably ice free. The evidence comes from high abundances of the subpolar planktic foraminifera Turborotalita quinqueloba in five newly analysed cores. The northern occurrence of this species is incompatible with perennial sea ice, which would be associated with a thick, low-salinity surface water. Instead, T. quinqueloba's ecological preference implies largely ice-free surface waters with seasonally elevated levels of primary productivity. In the modern ocean, this species thrives in the Fram Strait-Barents Sea 'Arctic-Atlantic gateway' region, implying that the necessary Atlantic Ocean-sourced water masses shoaled towards the surface during the Last Interglacial. This process reflects the ongoing Atlantification of the Arctic Ocean, currently restricted to the Eurasian Basin. Our results establish the Last Interglacial as a prime analogue for studying a seasonally ice-free Arctic Ocean, expected to occur this century. The warm Last Interglacial led to a seasonally ice-free Arctic Ocean and a transformation to Atlantic conditions, according to planktic foraminifera records from central Arctic Ocean sediment cores.

Earth's hydrological cycle is expected to intensify in response to global warming, with a wet-gets-wetter, dry-gets-drier response anticipated over the ocean. Subtropical regions (similar to 15 degrees-30 degrees N/S) are predicted to become drier, yet proxy evidence from past warm climates suggests these regions may be characterized by wetter conditions. Here we use an integrated data-modeling approach to reconstruct global and zonal-mean rainfall patterns during the early Eocene (similar to 56-48 million years ago). The Deep-Time Model Intercomparison Project (DeepMIP) model ensemble indicates that the mid-(30 degrees-60 degrees N/S) and high-latitudes (>60 degrees N/S) are characterized by a thermodynamically dominated hydrological response to warming and overall wetter conditions. The tropical band (0 degrees-15 degrees N/S) is also characterized by wetter conditions, with several DeepMIP models simulating narrowing of the Inter-Tropical Convergence Zone. However, the latter is not evident from the proxy data. The subtropics are characterized by negative precipitation-evaporation anomalies (i.e., drier conditions) in the DeepMIP models, but there is surprisingly large inter-model variability in mean annual precipitation (MAP). Intriguingly, we find that models with weaker meridional temperature gradients (e.g., CESM, GFDL) are characterized by a reduction in subtropical moisture divergence, leading to an increase in MAP. These model simulations agree more closely with our new proxy-derived precipitation reconstructions and other key climate metrics and imply that the early Eocene was characterized by reduced subtropical moisture divergence. If the meridional temperature gradient was even weaker than suggested by those DeepMIP models, circulation-induced changes may have outcompeted thermodynamic changes, leading to wetter subtropics. This highlights the importance of accurately reconstructing zonal temperature gradients when reconstructing past rainfall patterns. As the world warms, the atmosphere is able to hold more moisture however, this moisture will not fall evenly across the globe. Some regions are expected to become wetter, whereas other regions will become drier. This is the basis of the familiar paradigm wet-gets-wetter, dry-gets-drier and is largely supported by future model projections. However, evidence from the geological record contradicts this hypothesis and suggests that a warmer world could be characterized by wetter (rather than drier) subtropics. Here, we use an integrated data-modeling approach to investigate the hydrological response to warming during an ancient warm interval (the early Eocene, 56-48 million years ago). We show that models with weaker latitudinal temperature gradients are characterized by a reduction in subtropical moisture divergence. However, this was not sufficient to induce subtropical wetting. If the meridional temperature gradient was weaker than suggested by the models, circulation-induced changes may have lead to wetter subtropics. This work shows that the latitudinal temperature gradient is a key factor that influences hydroclimate in the subtropics, especially in past warm climates.

The total meridional heat transport (MHT) is relatively stable across different climates. Nevertheless, the strength of individual processes contributing to the total transport are not stable. Here we investigate the MHT and its main components especially in the atmosphere, in five coupled climate model simulations from the Deep-Time Model Intercomparison Project (DeepMIP). These simulations target the early Eocene climatic optimum, a geological time period with high CO₂ concentrations, analog to the upper range of end-of-century CO₂ projections. Preindustrial and early Eocene simulations, at a range of CO₂ levels are used to quantify the MHT changes in response to both CO₂ and non-CO₂ related forcings. We found that atmospheric poleward heat transport increases with CO₂, while oceanic poleward heat transport decreases. The non-CO₂ boundary conditions cause more MHT toward the South Pole, mainly through an increase in the southward oceanic heat transport. The changes in paleogeography increase the heat transport via transient eddies at the northern mid-latitudes in the Eocene. The Eocene Hadley cells do not transport more heat poleward, but due to the warmer atmosphere, especially the northern cell, circulate more heat in the tropics, than today. The monsoon systems' poleward latent heat transport increases with rising CO₂ concentrations, but this change is counterweighted by the globally smaller Eocene monsoon area. Our results show that the changes in the monsoon systems' latent heat transport is a robust feature of CO₂ warming, which is in line with the currently observed precipitation increase of present day monsoon systems.

A major step in the long-term Cenozoic evolution toward a glacially driven climate occurred at the Eocene–Oligocene transition (EOT), ∼34.44 to 33.65 million years ago (Ma). Evidence for high-latitude cooling and increased latitudinal temperature gradients across the EOT has been found in a range of marine and terrestrial environments. However, the timing and magnitude of temperature change in the North Atlantic remains highly unconstrained. Here, we use two independent organic geochemical palaeothermometers to reconstruct sea surface temperatures (SSTs) from the southern Labrador Sea (Ocean Drilling Program – ODP Site 647) across the EOT. The new SST records, now the most detailed for the North Atlantic through the 1 Myr leading up to the EOT onset, reveal a distinctive cooling step of ∼3 ^∘C (from 27 to 24 ^∘C), between 34.9 and 34.3 Ma, which is ∼500 kyr prior to Antarctic glaciation. This cooling step, when compared visually to other SST records, is asynchronous across Atlantic sites, signifying considerable spatiotemporal variability in regional SST evolution. However, overall, it fits within a phase of general SST cooling recorded across sites in the North Atlantic in the 5 Myr bracketing the EOT.

Such cooling might be unexpected in light of proxy and modelling studies suggesting the start-up of the Atlantic Meridional Overturning Circulation (AMOC) before the EOT, which should warm the North Atlantic. Results of an EOT modelling study (GFDL CM2.1) help reconcile this, finding that a reduction in atmospheric CO₂ from 800 to 400 ppm may be enough to counter the warming from an AMOC start-up, here simulated through Arctic–Atlantic gateway closure. While the model simulations applied here are not yet in full equilibrium, and the experiments are idealised, the results, together with the proxy data, highlight the heterogeneity of basin-scale surface ocean responses to the EOT thermohaline changes, with sharp temperature contrasts expected across the northern North Atlantic as positions of the subtropical and subpolar gyre systems shift. Suggested future work includes increasing spatial coverage and resolution of regional SST proxy records across the North Atlantic to identify likely thermohaline fingerprints of the EOT AMOC start-up, as well as critical analysis of the causes of inter-model responses to help better understand the driving mechanisms.

Based on a 7-member global circulation model ensemble from CMIP6/PMIP4, we compare the regional distribution of Arctic sea ice between a simulation representing the Last Interglacial (LIG) climate, with solar-forced warming, and an idealized future CO2-forced simulation with a similar annual sea-ice volume. The two simulations feature small but robust differences in the Central Arctic and Baffin Bay during summer, and larger differences at the sea-ice margins in the sub-Arctic Atlantic and North Pacific sectors during winter. Our results indicate that, under both forcings, sea ice persists north of Greenland until late summer, suggesting that the assumption that this region is the Last Ice Area is robust and holds for other climate states. However, we show that processes influencing sea-ice distribution in winter, such as Atlantification and sea-ice drift, differ and need to be further investigated.

The Last Interglacial (LIG) period, which had higher summer solar insolation than today, has been suggested as the last time that Arctic summers were ice free. However, the latest suite of Coupled Modelling Intercomparison Project 6 Paleoclimate (CMIP6-PMIP4) simulations of the LIG produce a wide range of Arctic summer minimum sea ice area (SIA) results, with a 30% to 96% reduction from the pre-industrial (PI) period. Sea ice proxies are also currently neither abundant nor consistent enough to determine the most realistic state. Here we estimate LIG minimum SIA indirectly through the use of 21 proxy records for LIG summer surface air temperature (SSAT) and 11 CMIP6-PMIP4 models for the LIG. We use two approaches. First, we use two tests to determine how skilful models are at simulating reconstructed Delta SSAT from proxy records (where Delta refers to LIG-PI). This identifies a positive correlation between model skill and the magnitude of Delta SIA: the most reliable models simulate a larger sea ice reduction. Averaging the two most skilful models yields an average SIA of 1 :3 x 10(6) km(2) for the LIG. This equates to a 4 :5 x 10(6) km(2) or 79% SIA reduction from the PI to the LIG. Second, across the 11 models, the averaged Delta SSAT at the 21 proxy locations and the pan-Arctic average Delta SSAT are inversely correlated with Delta SIA ( r = - 0.86 and 0.79, respectively). In other words, the models show that a larger Arctic warming is associated with a greater sea ice reduction. Using the proxy-record-averaged Delta SSAT of 4 :5 +/- 1 :7K and the relationship between Delta SSAT and Delta SIA suggests an estimated sea ice reduction of 4:2 +/- 1:4 x 10(6) km(2) or about 74% less sea ice than the PI period. The mean proxy-location Delta SSAT is well correlated with the Arctic-wide Delta SSAT north of 60 degrees N (r = D 0:97), and this relationship is used to show that the mean proxy record Delta SSAT is equivalent to an Arctic-wide warming of 3 :7 +/- 1 :5K at the LIG compared to the PI period. Applying this Arctic-wide Delta SSAT and its modelled relationship to Delta SIA, results in a similar estimate of LIG sea ice reduction of 4 :1 +/- 1 :2 x 10(6) km(2). These LIG climatological minimum SIA of 1.3 to 1.5 x 10(6) km(2) are close to the definition of a summer ice-free Arctic, which is a maximum sea ice extent of less than 1 x 10(6) km(2). The results of this study thus suggest that the Arctic likely experienced a mixture of ice-free and near-ice-free summers during the LIG.