Search

Hydration free energies of nonpolarizable monovalent atomic ions in transferable intermolecular potential four point fluctuating charge (TIP4P-FQ) are computed using several commonly employed ion-water force fields including two complete model sets recently developed for use with the simple water model with four sites and Drude polarizability and TIP4P water models. A simulation methodology is presented which incorporates a number of finite-system free energy corrections within the context of constant pressure molecular dynamics simulations employing the Ewald method and periodic boundary conditions. The agreement of the computed free energies and solvation structures with previously reported results for these models in finite droplet systems indicates good transferability of ion force fields from these water models to TIP4Q-FQ even when ion polarizability is neglected. To assess the performance of the ion models in TIP4P-FQ, we compare with consensus values for single-ion hydration free energies arising from recently improved cluster-pair estimates and a reevaluation of commonly cited, experimentally derived single-ion hydration free energies; we couple the observed consistency of these energies with a justification of the cluster-pair approximation in assigning single-ion hydration free energies to advocate the use of these consensus energies as a benchmark set in the parametrization of future ion force fields. [ABSTRACT FROM AUTHOR]

The Canadian Earth System Model version 5 (CanESM5) is a global model developed to simulate historical climate change and variability, to make centennial-scale projections of future climate, and to produce initialized seasonal and decadal predictions. This paper describes the model components and their coupling, as well as various aspects of model development, including tuning, optimization, and a reproducibility strategy. We also document the stability of the model using a long control simulation, quantify the model's ability to reproduce large-scale features of the historical climate, and evaluate the response of the model to external forcing. CanESM5 is comprised of three-dimensional atmosphere (T63 spectral resolution equivalent roughly to 2.8 ∘) and ocean (nominally 1 ∘) general circulation models, a sea-ice model, a land surface scheme, and explicit land and ocean carbon cycle models. The model features relatively coarse resolution and high throughput, which facilitates the production of large ensembles. CanESM5 has a notably higher equilibrium climate sensitivity (5.6 K) than its predecessor, CanESM2 (3.7 K), which we briefly discuss, along with simulated changes over the historical period. CanESM5 simulations contribute to the Coupled Model Intercomparison Project phase 6 (CMIP6) and will be employed for climate science and service applications in Canada. [ABSTRACT FROM AUTHOR]

The Canadian Earth System Model version 5 (CanESM5) is a global model developed to simulate historical climate change and variability, to make centennial scale projections of future climate, and to produce initialized seasonal and decadal predictions. This paper describes the model components and their coupling, as well as various aspects of model development, including tuning, optimization and a reproducibility strategy. We also document the stability of the model using a long control simulation, quantify the model's ability to reproduce large scale features of the historical climate, and evaluate the response of the model to external forcing. CanESM5 is comprised of three dimensional atmosphere (T63 spectral resolution/2.8°) and ocean (nominally 1°) general circulation models, a sea ice model, a land surface scheme, and explicit land and ocean carbon cycle models. The model features relatively coarse resolution and high throughput, which facilitates the production of large ensembles. CanESM5 has a notably higher equilibrium climate sensitivity (5.7 K) than its predecessor CanESM2 (3.8 K), which we briefly discuss, along with simulated changes over the historical period. CanESM5 simulations are contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6), and will be employed for climate science and service applications in Canada. [ABSTRACT FROM AUTHOR]

The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) aims to investigate the spread in simulations of sea-level and ocean climate change in response to CO2 forcing by atmosphere-ocean general circulation models (AOGCMs). It is particularly motivated by the uncertainties in projections of ocean heat uptake, global-mean sealevel rise due to thermal expansion and the geographical patterns of sea-level change due to ocean density and circulation change. FAFMIP has three tier-1 experiments, in which prescribed surface flux perturbations of momentum, heat and freshwater respectively are applied to the ocean in separate AOGCM simulations. All other conditions are as in the pre-industrial control. The prescribed fields are typical of pattern and magnitude of changes in these fluxes projected by AOGCMs for doubled CO2 concentration. Five groups have tested the experimental design with existing AOGCMs. Their results show diversity in the pattern and magnitude of changes, with some common qualitative features. Heat and water flux perturbation cause the dipole in sea-level change in the North Atlantic, while momentum and heat flux perturbation cause the gradient across the Antarctic Circumpolar Current. The Atlantic meridional overturning circulation (AMOC) declines in response to the heat flux perturbation, and there is a strong positive feedback on this effect due to the consequent cooling of sea-surface temperature in the North Atlantic, which enhances the local heat input to the ocean. The momentum and water flux perturbations do not substantially affect the AMOC. Heat is taken up largely as a passive tracer in the Southern Ocean, which is the region of greatest heat input, while the weakening of the AMOC causes redistribution of heat towards lower latitudes. Future analysis of these and other phenomena with the wider range of CMIP6 FAFMIP AOGCMs will benefit from new diagnostics of temperature and salinity tendencies, which will enable investigation of the model spread in behaviour in terms of physical processes as formulated in the models. [ABSTRACT FROM AUTHOR]

The ocean biogeochemistry components of two new versions of the Canadian Earth System Model (CanESM) are presented and compared to observations and other models. CanESM5 employs the same ocean biology model as CanESM2, whereas CanESM5-CanOE (Canadian Ocean Ecosystem model) is a new, more complex model developed for CMIP6, with multiple food chains, flexible phytoplankton elemental ratios, and a prognostic iron cycle. This new model is described in detail and the outputs (distributions of major tracers such as oxygen, dissolved inorganic carbon, and alkalinity, the iron and nitrogen cycles, plankton biomass, and historical trends in CO 2 uptake and export production) compared to CanESM5 and CanESM2, as well as to observations and other CMIP6 models. Both CanESM5 models show gains in skill relative to CanESM2, which are attributed primarily to improvements in ocean circulation. CanESM5-CanOE shows improved skill relative to CanESM5 for most major tracers at most depths. CanESM5-CanOE includes a prognostic iron cycle, and maintains high-nutrient/low-chlorophyll conditions in the expected regions (in CanESM2 and CanESM5, iron limitation is specified as a temporally static "mask"). Surface nitrate concentrations are biased low in the subarctic Pacific and equatorial Pacific, and high in the Southern Ocean, in both CanESM5 and CanESM5-CanOE. Export production in CanESM5-CanOE is among the lowest for CMIP6 models; in CanESM5, it is among the highest, but shows the most rapid decline after about 1980. CanESM5-CanOE shows some ability to simulate aspects of plankton community structure that a single-species model can not (e.g. seasonal dominance of large cells) but is biased towards low concentrations of zooplankton and detritus relative to phytoplankton. Cumulative ocean uptake of anthropogenic carbon dioxide through 2014 is lower in both CanESM5-CanOE (122 PgC) and CanESM5 (132 PgC) than in observation-based estimates (145 PgC) or the model ensemble mean (144 PgC). [ABSTRACT FROM AUTHOR]

The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) aims to investigate the spread in simulations of sea-level and ocean climate change in response to CO2 forcing by atmosphere-ocean general circulation models (AOGCMs). It is particularly motivated by the uncertainties in projections of ocean heat uptake, global-mean sea-level rise due to thermal expansion and the geographical patterns of sea-level change due to ocean density and circulation change. FAFMIP has three tier-1 experiments, in which prescribed surface flux perturbations of momentum, heat and freshwater respectively are applied to the ocean in separate AOGCM simulations. All other conditions are as in the pre-industrial control. The prescribed fields are typical of pattern and magnitude of changes in these fluxes projected by AOGCMs for doubled CO2 concentration. Five groups have tested the experimental design with existing AOGCMs. Their results show diversity in the pattern and magnitude of changes, with some common qualitative features. Heat and water flux perturbation cause the dipole in sea-level change in the North Atlantic, while momentum and heat flux perturbation cause the gradient across the Antarctic Circumpolar Current. The Atlantic Meridional Overturning Circulation (AMOC) declines in response to the heat flux perturbation, and there is a strong positive feedback on this effect due to the consequent cooling of sea surface temperature in the North Atlantic, which enhances the local heat input to the ocean. The momentum and water flux perturbations do not substantially affect the AMOC. Heat is taken up largely as a passive tracer in the Southern Ocean, which is the region of greatest heat input, but elsewhere heat is actively redistributed towards lower latitude. Future analysis of these and other phenomena with the wider range of CMIP6 FAFMIP AOGCMs will benefit from new diagnostics of temperature and salinity tendencies, which will enable investigation of the model spread in behaviour in terms of physical processes as formulated in the models. [ABSTRACT FROM AUTHOR]

Several recent studies have shown that ocean western boundaries are the primary regions of eddy energy dissipation. Globally, the eddy energy sinks have been estimated to integrate to about 0.2 TW. This is a sizable fraction of the tidal energy dissipation in the deep oceanic interior, estimated at about 1.0 TW and contributing to diapycnal mixing. The authors conduct sensitivity experiments with an ocean general circulation model assuming that the eddy energy is scattered into high-wavenumber vertical modes, resulting in energy dissipation and locally enhanced diapycnal mixing. When only the tidal energy dissipation maintains diapycnal mixing, the overturning circulation, and stratification in the deep ocean are too weak. With the addition of the eddy dissipation, the deep-ocean thermal structure becomes closer to that observed and the overturning circulation and stratification in the abyss become stronger. Furthermore, the mixing associated with the eddy dissipation can, on its own, drive a relatively strong overturning. The stratification and overturning in the deep ocean are sensitive to the vertical structure of diapycnal mixing. When most of this energy dissipates within 300 m above the bottom, the abyssal overturning and stratification are too weak. Allowing for the dissipation to penetrate higher in the water column, such as suggested by recent observations, results in stronger stratification and meridional circulation. Zonal circulation is also affected. In particular, the Drake Passage transport becomes closer to its observational estimates with the increase in the vertical scale for turbulence above topography. Consistent with some theoretical models, the Drake Passage transport increases with the increase in the mean upper-ocean diffusivity. [ABSTRACT FROM AUTHOR]

Subduction, water mass transformation, and transport rates in the Indo-Pacific Ocean are diagnosed in a recent version of the Canadian Centre for Climate Modelling and Analysis coupled model. It is found that the subduction across the base of the winter mixed layer is dominated by the lateral transfer, particularly within the relatively dense water classes corresponding to the densest mode and intermediate waters. However, within lighter densities, including those characterizing the lighter varieties of mode waters, the vertical transfer has a strong positive input to the net subduction. The upper-ocean volume transports across 30°°N and 32°°S are largest within the density classes that correspond to mode waters. In the North Pacific, the buoyancy flux converts the near-surface waters mostly to denser water classes, whereas in the Southern Ocean the surface waters are transformed both to lighter and denser water classes, depending on the density. In response to a doubling of CO2, the subduction, transformation, and transport of mode waters in both hemispheres shift to lighter densities but do not change significantly, whereas the subduction of intermediate waters decreases. The area of large winter mixed layer depths decreases, particularly in the Southern Hemisphere. In the low latitudes, the thermocline water flux that enters the tropical Pacific via the western boundary flows generally increases. However, its anomaly has a complex structure, so that integrated estimates can be sensitive to the isopycnal ranges. The upper part of the Equatorial Undercurrent (EUC) strengthens in the warmer climate, whereas its lower part weakens. The anomaly in the EUC closely follows the anomaly in stratification along the equator. The Indonesian Throughflow transport decreases with part of it being redirected eastward. This part joins with the intensified equatorward thermocline flows at the western boundaries and contributes to the EUC anomaly. [ABSTRACT FROM AUTHOR]