Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Carbon-Nitrogen Coupling Under Three Schemes of Model Representation: A Traceability Analysis

The interaction between terrestrial carbon (C) and nitrogen (N) cycles has been incorporated into more and more land surface models. However, the scheme of C-N coupling differs greatly among models, and how these diverse representations of C-N interactions will affect C-cycle modeling remains unclear. In this study, we explored how the simulated ecosystem C storage capacity in the terrestrial ecosystem (TECO) model varied with three different commonly used schemes of C-N coupling. The three schemes (SM1, SM2, and SM3) have been used in three different coupled C-N models (i.e., TECO-CN, CLM 4.5, and O-CN, respectively). They differ mainly in the stoichiometry of C and N in vegetation and soils, plant N uptake strategies, down regulation of photosynthesis, and the pathways of N import. We incorporated the three C-N coupling schemes into the C-only version of the TECO model and evaluated their impacts on the C cycle with a traceability framework. Our results showed that all three of the C-N schemes caused significant reductions in steady-state C storage capacity compared with the C-only version with magnitudes of −23%, −30%, and −54% for SM1, SM2, and SM3, respectively. This reduced C storage capacity was mainly derived from the combined effects of decreases in net primary productivity (NPP; −29%, −15%, and −45%) and changes in mean C residence time (MRT; 9%, −17%, and −17%) for SM1, SM2, and SM3, respectively. The differences in NPP are mainly attributed to the different assumptions on plant N uptake, plant tissue C : N ratio, down regulation of photosynthesis, and biological N fixation. In comparison, the alternative representations of the plant vs. microbe competition strategy and the plant N uptake, combined with the flexible C : N ratio in vegetation and soils, led to a notable spread in MRT. These results highlight the fact that the diverse assumptions on N processes represented by different C-N coupled models could cause additional uncertainty for land surface models. Understanding their difference can help us improve the capability of models to predict future biogeochemical cycles of terrestrial ecosystems.

Du, Zhenggang↗
NASA NTRS2017

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗
NASA NTRS2023

Convective Impact on the Global Lower Stratospheric Water Vapor Budget

Water vapor in the stratosphere is primarily controlled by temperatures in the tropical upper troposphere and lower stratosphere. However, the direct impact of deep convection on the global lower stratospheric water vapor budget is still an actively debated issue. Two complementary modeling approaches are used to investigate the convective impact in boreal winter and summer. Convective influence is diagnosed by tracing trajectories through convective cloud top altitude fields derived from global rainfall and brightness temperature data. Backward trajectory model simulations coupled with a detailed treatment of cloud microphysical processes indicate that convection moistens the global lower stratosphere by approximately 0.3 ppmv in boreal winter and summer 2010. The diurnal peak in convection is responsible for about half of the total convective moistening during winter and nearly all of the convective moistening during summer. Deep convective clouds overshooting the tropopause have relatively minor effect on global lower stratospheric water vapor. A forward trajectory model coupled with a simplified cloud module is used to estimate the relative magnitude of the interannual variability of the convective impact. Combining the results from the two models, we find that the convective impact on the global lower stratospheric water vapor during 2006-2016 is approximately 0.3 ppmv with year-to-year variations of up to 0.1 ppmv. An important mechanism of convective hydration of the lower stratosphere is via the detrainment of saturated air and ice into the tropical uppermost troposphere and the subsequent upward transport of some of these moist air parcels across relatively warm and subsaturated tropopause.

CONVECTION↗
NASA NTRS1987

A coupled dynamic-thermodynamic model of an ice-ocean system in the marginal ice zone

Thermodynamics are incorporated into a coupled ice-ocean model in order to investigate wind-driven ice-ocean processes in the marginal zone. Upswelling at the ice edge which is generated by the difference in the ice-air and air-water surface stresses is found to give rise to a strong entrainment by drawing the pycnocline closer to the surface. Entrainment is shown to be negligible outside the areas affected by the ice edge upswelling. If cooling at the top is included in the model, the heat and salt exchanges are further enhanced in the upswelling areas. It is noted that new ice formation occurs in the region not affected by ice edge upswelling, and it is suggested that the high-salinity mixed layer regions (with a scale of a few Rossby radii of deformation) will overturn due to cooling, possibly contributing to the formation of deep water.

Hakkinen, Sirpa↗
NASA NTRS2019

The Connected Isotopic Water Cycle in the Community Earth System Model Version 1

Because of the pervasive role of water in the Earth system, the relative abundances of stable isotopologues of water are valuable for understanding atmospheric, oceanic, and biospheric processes, and for interpreting paleoclimate proxy reconstructions. Isotopologues are transported by both largescale and turbulent flows, and the ratio of heavy to light isotopologues changes due to fractionation that can accompany condensation and evaporation processes. Correctly predicting the isotopic distributions requires resolving the relationships between largescale ocean and atmospheric circulation and smallerscale hydrological processes, which can be accomplished within a coupled climate modeling framework. Here we present the water isotopeenabled version of the Community Earth System Model version 1 (iCESM1), which simulates global variations in water isotopic ratios in the atmosphere, land, ocean, and sea ice. In a transient Last Millennium simulation covering the 850-2005 period, iCESM1 correctly captures the latetwentiethcentury structure of δ(exp 18)O and δD over the global oceans, with more limited accuracy over land. The relationship between salinity and seawater δ(exp 18)O is also well represented over the observational period, including interbasin variations. We illustrate the utility of coupled, isotopeenabled simulations using both Last Millennium simulations and freshwater hosing experiments with iCESM1. Closing the isotopic mass balance between all components of the coupled model provides new confidence in the underlying depiction of the water cycle in CESM, while also highlighting areas where the underlying hydrologic balance can be improved. The iCESM1 is poised to be a vital community resource for ongoing model development with both modern and paleoclimate applications.

Brady, E.↗
NASA NTRS2013

Stratospheric Ozone Response in Experiments G3 and G4 of the Geoengineering Model Intercomparison Project (GeoMIP)

Geoengineering with stratospheric sulfate aerosols has been proposed as a means of temporarily cooling the planet, alleviating some of the side effects of anthropogenic CO2 emissions. However, one of the known side effects of stratospheric injections of sulfate aerosols is a decrease in stratospheric ozone. Here we show results from two general circulation models and two coupled chemistry climate models that have simulated stratospheric sulfate aerosol geoengineering as part of the Geoengineering Model Intercomparison Project (GeoMIP). Changes in photolysis rates and upwelling of ozone-poor air in the tropics reduce stratospheric ozone, suppression of the NOx cycle increases stratospheric ozone, and an increase in available surfaces for heterogeneous chemistry modulates reductions in ozone. On average, the models show a factor 20-40 increase of the sulfate aerosol surface area density (SAD) at 50 hPa in the tropics with respect to unperturbed background conditions and a factor 3-10 increase at mid-high latitudes. The net effect for a tropical injection rate of 5 Tg SO2 per year is a decrease in globally averaged ozone by 1.1-2.1 DU in the years 2040-2050 for three models which include heterogeneous chemistry on the sulfate aerosol surfaces. GISS-E2-R, a fully coupled general circulation model, performed simulations with no heterogeneous chemistry and a smaller aerosol size; it showed a decrease in ozone by 9.7 DU. After the year 2050, suppression of the NOx cycle becomes more important than destruction of ozone by ClOx, causing an increase in total stratospheric ozone. Contribution of ozone changes in this experiment to radiative forcing is 0.23 W m-2 in GISS-E2-R and less than 0.1 W m-2 in the other three models. Polar ozone depletion, due to enhanced formation of both sulfate aerosol SAD and polar stratospheric clouds, results in an average 5 percent increase in calculated surface UV-B.

aerosol↗
NASA NTRSDate not supplied

Coupling a Computational Fluid Dynamics (CFD) Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

heat transfer↗
NASA NTRS1992

A numerical model for atomization-spray coupling in liquid rocket thrust chambers

The physical process of atomization is an important consideration in the stable operation of liquid rocket engines. Many spray combustion computational fluid dynamics (CFD) codes do not include an atomization sub-model but assume arbitrary drop size distributions, drop initial locations, and velocities. A method of coupling an atomization model with the spray model in a REFLEQS CFD code is presented. This method is based on a jet-embedding technique in which the equations governing the liquid jet core are solved separately using the surrounding gas phase conditions. The droplet initial conditions are calculated using a stability analysis appropriate for the atomization regime of liquid jet break-up. This novel coupling model is used to analyze the SSME fuel preburner single injector flow. Results of the diffusion flame characteristics in a single injection element will be presented. The effect of relative velocity, mixture ratio and droplet initial conditions will be shown. The predictions of present atomization model is compared with that of the widely used CICM correlation. The results are also compared with the predictions of volume-of-fluid method.

M. G. Giridharan↗
NASA NTRS1971

Apollo cryogenic integrated systems program

The integrated systems program is capable of simulating both nominal and anomalous operation of the Apollo cryogenics storage system (CSS). Two versions of the program exist; one for the Apollo 14 configuration and the other for J Type Mission configurations. The program consists of two mathematical models which are dynamically coupled. A model of the CSS components and lines determines the oxygen and hydrogen flowrate from each storage tank given the tank pressures and temperatures, and the electrical power subsystem and environmental control subsystem flow demands. Temperatures and pressures throughout the components and lines are also determined. A model of the CSS tankage determines the pressure and temperatures in the tanks given the flowrate from each tank and the thermal environment. The model accounts for tank stretch and includes simplified oxygen tank heater and stratification routines. The program is currently operational on the Univac 1108 computer.

Seto, R. K. M.↗
NASA NTRS1988

Abrupt climate change and extinction events

There is a growing body of theoretical and empirical support for the concept of instabilities in the climate system, and indications that abrupt climate change may in some cases contribute to abrupt extinctions. Theoretical indications of instabilities can be found in a broad spectrum of climate models (energy balance models, a thermohaline model of deep-water circulation, atmospheric general circulation models, and coupled ocean-atmosphere models). Abrupt transitions can be of several types and affect the environment in different ways. There is increasing evidence for abrupt climate change in the geologic record and involves both interglacial-glacial scale transitions and the longer-term evolution of climate over the last 100 million years. Records from the Cenozoic clearly show that the long-term trend is characterized by numerous abrupt steps where the system appears to be rapidly moving to a new equilibrium state. The long-term trend probably is due to changes associated with plate tectonic processes, but the abrupt steps most likely reflect instabilities in the climate system as the slowly changing boundary conditions caused the climate to reach some threshold critical point. A more detailed analysis of abrupt steps comes from high-resolution studies of glacial-interglacial fluctuations in the Pleistocene. Comparison of climate transitions with the extinction record indicates that many climate and biotic transitions coincide. The Cretaceous-Tertiary extinction is not a candidate for an extinction event due to instabilities in the climate system. It is quite possible that more detailed comparisons and analysis will indicate some flaws in the climate instability-extinction hypothesis, but at present it appears to be a viable candidate as an alternate mechanism for causing abrupt environmental changes and extinctions.

Crowley, Thomas J.↗
NASA NTRS1992

The Global Ocean Observing System

A Global Ocean Observing System (GOOS) should be established now with international coordination (1) to address issues of global change, (2) to implement operational ENSO forecasts, (3) to provide the data required to apply global ocean circulation models, and (4) to extract the greatest value from the one billion dollar investment over the next ten years in ocean remote sensing by the world's space agencies. The objectives of GOOS will focus on climatic and oceanic predictions, on assessing coastal pollution, and in determining the sustainability of living marine resources and ecosystems. GOOS will be a complete system including satellite observations, in situ observations, numerical modeling of ocean processes, and data exchange and management. A series of practical and economic benefits will be derived from the information generated by GOOS. In addition to the marine science community, these benefits will be realized by the energy industries of the world, and by the world's fisheries. The basic oceanic variables that are required to meet the oceanic and predictability objectives of GOOS include wind velocity over the ocean, sea surface temperature and salinity, oceanic profiles of temperature and salinity, surface current, sea level, the extent and thickness of sea ice, the partial pressure of CO2 in surface waters, and the chlorophyll concentration of surface waters. Ocean circulation models and coupled ocean-atmosphere models can be used to evaluate observing system design, to assimilate diverse data sets from in situ and remotely sensed observations, and ultimately to predict future states of the system. The volume of ocean data will increase enormously over the next decade as new satellite systems are launched and as complementary in situ measuring systems are deployed. These data must be transmitted, quality controlled, exchanged, analyzed, and archived with the best state-of-the-art computational methods.

Kester, Dana↗
NASA NTRS1998

Free-Suspension Residual Flexibility Testing of Space Station Pathfinder: Comparison to Fixed-Base Results

Application of the free-suspension residual flexibility modal test method to the International Space Station Pathfinder structure is described. The Pathfinder, a large structure of the general size and weight of Space Station module elements, was also tested in a large fixed-base fixture to simulate Shuttle Orbiter payload constraints. After correlation of the Pathfinder finite element model to residual flexibility test data, the model was coupled to a fixture model, and constrained modes and frequencies were compared to fixed-base test. modes. The residual flexibility model compared very favorably to results of the fixed-base test. This is the first known direct comparison of free-suspension residual flexibility and fixed-base test results for a large structure. The model correlation approach used by the author for residual flexibility data is presented. Frequency response functions (FRF) for the regions of the structure that interface with the environment (a test fixture or another structure) are shown to be the primary tools for model correlation that distinguish or characterize the residual flexibility approach. A number of critical issues related to use of the structure interface FRF for correlating the model are then identified and discussed, including (1) the requirement of prominent stiffness lines, (2) overcoming problems with measurement noise which makes the antiresonances or minima in the functions difficult to identify, and (3) the use of interface stiffness and lumped mass perturbations to bring the analytical responses into agreement with test data. It is shown that good comparison of analytical-to-experimental FRF is the key to obtaining good agreement of the residual flexibility values.

Tinker, Michael L.↗
NASA NTRS2005

LWS Proposal to Provide Scientific Guidance and Modeling Support for the Ionospheric Mapping Mission

A data assimilation system for specifying the thermospheric density has been developed over the last several years. This system ingests GRACE/CHAMP-type in situ as well as SSULI/SSUSI remote sensing observations while making use of a physical model, the Coupled Thermosphere-Ionosphere Model (CTIM) (Fuller-Rowel1 et al., 1996). The Kalman filter was implemented as the backbone to the data assimilation system, which provides a statistically 'best' estimate as well as an estimate of the error in its state. The system was tested using a simulated thermosphere and observations. CHAMP data were then used to provide the system with a real data source. The results of this study are herein.

Richmond, Arthur D.↗
NASA NTRS2012

Understanding Differences in Upper Stratospheric Ozone Response to Changes in Chlorine and Temperature as Computed Using CCMVal Models

Projections of future ozone levels are made using models that couple a general circulation model with a representation of atmospheric photochemical processes, allowing interactions among photochemical processes, radiation, and dynamics. Such models are known as chemistry and climate models (CCMs). Although developed from common principles and subject to the same boundary conditions, simulated ozone time series vary for projections of changes in ozone depleting substances (ODSs) and greenhouse gases. In the upper stratosphere photochemical processes control ozone level, and ozone increases as ODSs decrease and temperature decreases due to greenhouse gas increase. Simulations agree broadly but there are quantitative differences in the sensitivity of ozone to chlorine and to temperature. We obtain insight into these differences in sensitivity by examining the relationship between the upper stratosphere annual cycle of ozone and temperature as produced by a suite of models. All simulations conform to expectation in that ozone is less sensitive to temperature when chlorine levels are highest because chlorine catalyzed loss is nearly independent of temperature. Differences in sensitivity are traced to differences in simulated temperature, ozone and reactive nitrogen when chlorine levels are close to background. This work shows that differences in the importance of specific processes underlie differences in simulated sensitivity of ozone to composition change. This suggests a) the multi-model mean is not a best estimate of the sensitivity of upper ozone to changes in ODSs and temperature; b) the spread of values is not an appropriate measure of uncertainty.

Douglass, A. R.↗
NASA NTRS2015

CFD Modeling of the Multipurpose Hydrogen Test Bed (MHTB) Self-Pressurization and Spray Bar Mixing Experiments in Normal Gravity: Effect of the Accommodation Coefficient on the Tank Pressure

A CFD model for simulating the self-pressurization of a large scale liquid hydrogen storage tank is utilized in this paper to model the MHTB self-pressurization experiment. The kinetics-based Schrage equation is used to account for the evaporative and condensi ng interfacial mass flows in this model. The effect of the accommodation coefficient for calculating the interfacial mass transfer rate on the tank pressure during tank selfpressurization is studied. The values of the accommodation coefficient which were considered in this study vary from 1.0e-3 to 1.0e-1 for the explicit VOF model and from 1.0e-4 to 1.0e-3 for the implicit VOF model. The ullage pressure evolutions are compared against experimental data. A CFD model for controlling pressure in cryogenic storage tanks by spraying cold liquid into the ullage is also presented. The Euler-Lagrange approach is utilized for tracking the spray droplets and for modeling the interaction between the droplets and the continuous phase (ullage). The spray model is coupled with the VOF model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. Droplet-ullage heat and mass transfer are modeled. The flow, temperature, and interfacial mass flux, as well as droplets trajectories, size distribution and temperatures predicted by the model are presented. The ul lage pressure and vapor temperature evolutions are compared with experimental data obtained from the MHTB spray bar mixing experiment. The effect of the accommodation coefficient for calculating the interfacial and droplet mass transfer rates on the tank pressure during mixing of the vapor using spray is studied. The values used for the accommodation coefficient at the interface vary from 1.0e-5 to 1.0e-2. The droplet accommodation coefficient values vary from 2.0e-6 to 1.0e-4.

cryogenic fluid storage↗
NASA NTRS2016

Modeling Droplet Heat and Mass Transfer during Spray Bar Pressure Control of the Multipurpose Hydrogen Test Bed (MHTB) Tank in Normal Gravity

A CFD model for simulating pressure control in cryogenic storage tanks through the injection of a subcooled liquid into the ullage is presented and applied to the 1g MHTB spray bar cooling experiments. An Eulerian-Lagrangian approach is utilized to track the spray droplets and capture the interaction between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. A new model for calculating the droplet-ullage heat and mass transfer is developed. In this model, a droplet is allowed to warm up to the saturation temperature corresponding to the ullage vapor pressure, after which it evaporates while remaining at the saturation temperature. The droplet model is validated against the results of the MHTB spray-bar cooling experiments with 50% and 90% tank fill ratios. The predictions of the present T-sat based model are compared with those of a previously developed kinetic-based droplet mass transfer model. The predictions of the two models regarding the evolving tank pressure and temperature distributions, as well as the droplets' trajectories and temperatures, are examined and compared in detail. Finally, the ullage pressure and local vapor and liquid temperature evolutions are validated against the corresponding data provided by the MHTB spray bar mixing experiment.

mass transfer↗
NASA NTRS2021

Mid-Pliocene Atlantic Meridional Overturning Circulation simulated in PlioMIP2

In the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), coupled climate models have been used to simulate an interglacial climate during the mid-Piacenzian warm period (mPWP; 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), poleward ocean heat transport and sea surface warming in the Atlantic simulated with these models. In PlioMIP2, all models simulate an intensified mid-Pliocene AMOC. However, there is no consistent response in the simulated Atlantic ocean heat transport nor in the depth of the Atlantic overturning cell. The models show a large spread in the simulated AMOC maximum, the Atlantic ocean heat transport and the surface warming in the North Atlantic. Although a few models simulate a surface warming of ~8–12°C in the North Atlantic, similar to the reconstruction from Pliocene Research, Interpretation and Synoptic Mapping (PRISM) version 4, most models appear to underestimate this warming. The large model spread and model–data discrepancies in the PlioMIP2 ensemble do not support the hypothesis that an intensification of the AMOC, together with an increase in northward ocean heat transport, is the dominant mechanism for the mid-Pliocene warm climate over the North Atlantic.

Pliocene Model Intercomparison Project Phase 2 (Pl↗