•We focus on ACC and Southern Ocean MOC during 1958–2007 in 17 CORE-II forced models.•Most CORE-II simulations are close to eddy saturation.•Most CORE-II simulations are far from showing signs of ...eddy compensation.•Constant in time or space k results in poor representation of mesoscale eddy effects.•MOC has larger sensitivity than ACC transport even in eddy saturated state.
In the framework of the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II), we present an analysis of the representation of the Antarctic Circumpolar Current (ACC) and Southern Ocean meridional overturning circulation (MOC) in a suite of seventeen global ocean–sea ice models. We focus on the mean, variability and trends of both the ACC and MOC over the 1958–2007 period, and discuss their relationship with the surface forcing. We aim to quantify the degree of eddy saturation and eddy compensation in the models participating in CORE-II, and compare our results with available observations, previous fine-resolution numerical studies and theoretical constraints. Most models show weak ACC transport sensitivity to changes in forcing during the past five decades, and they can be considered to be in an eddy saturated regime. Larger contrasts arise when considering MOC trends, with a majority of models exhibiting significant strengthening of the MOC during the late 20th and early 21st century. Only a few models show a relatively small sensitivity to forcing changes, responding with an intensified eddy-induced circulation that provides some degree of eddy compensation, while still showing considerable decadal trends. Both ACC and MOC interannual variabilities are largely controlled by the Southern Annular Mode (SAM). Based on these results, models are clustered into two groups. Models with constant or two-dimensional (horizontal) specification of the eddy-induced advection coefficient κ show larger ocean interior decadal trends, larger ACC transport decadal trends and no eddy compensation in the MOC. Eddy-permitting models or models with a three-dimensional time varying κ show smaller changes in isopycnal slopes and associated ACC trends, and partial eddy compensation. As previously argued, a constant in time or space κ is responsible for a poor representation of mesoscale eddy effects and cannot properly simulate the sensitivity of the ACC and MOC to changing surface forcing. Evidence is given for a larger sensitivity of the MOC as compared to the ACC transport, even when approaching eddy saturation. Future process studies designed for disentangling the role of momentum and buoyancy forcing in driving the ACC and MOC are proposed.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
•Evaluation of 15 models participating in Phase II of the Coordinated Ocean-ice Reference Experiments (CORE-II).•Focus on the representation of upper ocean water mass properties and sea-ice processes ...in the Southern Ocean.•Both the mean states and trends over the 1988–2007 period are presented.•The model solutions agree in their sea ice concentration, surface buoyancy fluxes, and summer mixed layer depths.•Significant differences in model representation exist for the winter mixed layer depths, and the trends in sea ice extent.
We characterise the representation of the Southern Ocean water mass structure and sea ice within a suite of 15 global ocean-ice models run with the Coordinated Ocean-ice Reference Experiment Phase II (CORE-II) protocol. The main focus is the representation of the present (1988–2007) mode and intermediate waters, thus framing an analysis of winter and summer mixed layer depths; temperature, salinity, and potential vorticity structure; and temporal variability of sea ice distributions. We also consider the interannual variability over the same 20 year period. Comparisons are made between models as well as to observation-based analyses where available.
The CORE-II models exhibit several biases relative to Southern Ocean observations, including an underestimation of the model mean mixed layer depths of mode and intermediate water masses in March (associated with greater ocean surface heat gain), and an overestimation in September (associated with greater high latitude ocean heat loss and a more northward winter sea-ice extent). In addition, the models have cold and fresh/warm and salty water column biases centred near 50°S. Over the 1988–2007 period, the CORE-II models consistently simulate spatially variable trends in sea-ice concentration, surface freshwater fluxes, mixed layer depths, and 200–700 m ocean heat content. In particular, sea-ice coverage around most of the Antarctic continental shelf is reduced, leading to a cooling and freshening of the near surface waters. The shoaling of the mixed layer is associated with increased surface buoyancy gain, except in the Pacific where sea ice is also influential. The models are in disagreement, despite the common CORE-II atmospheric state, in their spatial pattern of the 20-year trends in the mixed layer depth and sea-ice.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
GISS‐E2.1: Configurations and Climatology Kelley, Maxwell; Schmidt, Gavin A; Nazarenko, Larissa ...
Journal of advances in modeling earth systems,
August 2020, Volume:
12, Issue:
8
Journal Article
Peer reviewed
Open access
This paper describes the GISS‐E2.1 contribution to the Coupled Model Intercomparison Project, Phase 6 (CMIP6). This model version differs from the predecessor model (GISS‐E2) chiefly due to ...parameterization improvements to the atmospheric and ocean model components, while keeping atmospheric resolution the same. Model skill when compared to modern era climatologies is significantly higher than in previous versions. Additionally, updates in forcings have a material impact on the results. In particular, there have been specific improvements in representations of modes of variability (such as the Madden‐Julian Oscillation and other modes in the Pacific) and significant improvements in the simulation of the climate of the Southern Oceans, including sea ice. The effective climate sensitivity to 2xCO2 is slightly higher than previously at 2.7‐‐3.1°C (depending on version), and is a result of lower CO2 radiative forcing and stronger positive feedbacks.
Full text
Available for:
DOBA, FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, SIK, UILJ, UKNU, UL, UM, UPUK
•Mesoscale eddies modulate the strength of global overturning cells in ocean models.•The strength and depth of the Atlantic Meridional Overturning Circulation (AMOC) cell co-vary.•The AMOC is ...sensitive to the representation of eddies in the Southern Ocean.
The dependence of the depth and strength of the ocean’s global meridional overturning cells (MOC) on the specification of mesoscale eddy diffusivity (K) is explored in two ocean models. The GISS and MIT ocean models are driven by the same prescribed forcing fields, configured in similar ways, spun up to equilibrium for a range of K’s and the resulting MOCs mapped and documented. Scaling laws implicit in modern theories of the MOC are used to rationalize the results. In all calculations the K used in the computation of eddy-induced circulation and that used in the representation of eddy stirring along neutral surfaces, is set to the same value but is changed across experiments. We are able to connect changes in the strength and depth of the Atlantic MOC, the southern ocean upwelling MOC, and the deep cell emanating from Antarctica, to changes in K.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK, ZRSKP
We present different simulations of primary atomization using an adaptive Volume-of-Fluid method based on octree meshes. The use of accurate numerical schemes for mesh adaptation, Volume-of-Fluid ...advection and balanced force surface tension calculation implemented in Gerris, the code used to perform the simulations included in this work, has made possible to carry out accurate simulations with characteristic scales spreading over several orders of magnitude. The code is validated by comparisons with the temporal linear theory for moderate density and viscosity ratios, which basically corresponds to atomization processes in high pressure chambers. In order to show the potential of the code in different scenarios related to atomization, preliminary results are shown in relation with the study of the two-dimensional and 3D temporal and spatial problem, the influence of the injector and the vortex generated inside the chamber, and the effect of swirling at high Reynolds numbers.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
•Phase II of the Coordinated Ocean-ice Reference Experiments (CORE-II) is introduced.•Solutions from CORE-II simulations from eighteen participating models are presented.•Mean states in the North ...Atlantic with a focus on AMOC are examined.•The North Atlantic solutions differ substantially among the models.•Many factors, including parameterization choices, contribute to these differences.
Simulation characteristics from eighteen global ocean–sea-ice coupled models are presented with a focus on the mean Atlantic meridional overturning circulation (AMOC) and other related fields in the North Atlantic. These experiments use inter-annually varying atmospheric forcing data sets for the 60-year period from 1948 to 2007 and are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The protocol for conducting such CORE-II experiments is summarized. Despite using the same atmospheric forcing, the solutions show significant differences. As most models also differ from available observations, biases in the Labrador Sea region in upper-ocean potential temperature and salinity distributions, mixed layer depths, and sea-ice cover are identified as contributors to differences in AMOC. These differences in the solutions do not suggest an obvious grouping of the models based on their ocean model lineage, their vertical coordinate representations, or surface salinity restoring strengths. Thus, the solution differences among the models are attributed primarily to use of different subgrid scale parameterizations and parameter choices as well as to differences in vertical and horizontal grid resolutions in the ocean models. Use of a wide variety of sea-ice models with diverse snow and sea-ice albedo treatments also contributes to these differences. Based on the diagnostics considered, the majority of the models appear suitable for use in studies involving the North Atlantic, but some models require dedicated development effort.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
•Inter-annual to decadal variability in AMOC from CORE-II simulations is presented.•AMOC variability shows three stages, with maximum transports in mid- to late-1990s.•North Atlantic temporal ...variability features are in good agreement among simulations.•Such agreements suggest variability is dictated by the atmospheric data sets.•Simulations differ in spatial structures of variability due to ocean dynamics.
Simulated inter-annual to decadal variability and trends in the North Atlantic for the 1958–2007 period from twenty global ocean – sea-ice coupled models are presented. These simulations are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The study is Part II of our companion paper (Danabasoglu et al., 2014) which documented the mean states in the North Atlantic from the same models. A major focus of the present study is the representation of Atlantic meridional overturning circulation (AMOC) variability in the participating models. Relationships between AMOC variability and those of some other related variables, such as subpolar mixed layer depths, the North Atlantic Oscillation (NAO), and the Labrador Sea upper-ocean hydrographic properties, are also investigated. In general, AMOC variability shows three distinct stages. During the first stage that lasts until the mid- to late-1970s, AMOC is relatively steady, remaining lower than its long-term (1958–2007) mean. Thereafter, AMOC intensifies with maximum transports achieved in the mid- to late-1990s. This enhancement is then followed by a weakening trend until the end of our integration period. This sequence of low frequency AMOC variability is consistent with previous studies. Regarding strengthening of AMOC between about the mid-1970s and the mid-1990s, our results support a previously identified variability mechanism where AMOC intensification is connected to increased deep water formation in the subpolar North Atlantic, driven by NAO-related surface fluxes. The simulations tend to show general agreement in their temporal representations of, for example, AMOC, sea surface temperature (SST), and subpolar mixed layer depth variabilities. In particular, the observed variability of the North Atlantic SSTs is captured well by all models. These findings indicate that simulated variability and trends are primarily dictated by the atmospheric datasets which include the influence of ocean dynamics from nature superimposed onto anthropogenic effects. Despite these general agreements, there are many differences among the model solutions, particularly in the spatial structures of variability patterns. For example, the location of the maximum AMOC variability differs among the models between Northern and Southern Hemispheres.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
Simulations of the CMIP6 historical period 1850‐‐2014, characterized by the emergence of anthropogenic climate drivers like greenhouse gases, are presented for different configurations of the NASA ...Goddard Institute for Space Studies (GISS) Earth System ModelE2.1. The GISS‐E2.1 ensembles are more sensitive to forcing than their CMIP5 predecessors (GISS‐E2), but warm less during recent decades due to reduced total forcing. This forcing reduction is attributed to an increase in longwave opacity in pre‐industrial simulations, resulting in an atmosphere less sensitive to further increases in opacity that result from rising greenhouse gas concentrations. This demonstrates the importance of the base climatology to forcing and forced climate trends. Most model versions match observed temperature trends since 1979 from the ocean to the stratosphere. The choice of ocean model is important to the transient climate response, as found in CMIP5 GISS‐E2: the model that more efficiently exports heat to the deep ocean shows a smaller rise in tropospheric temperature. Model sea level rise over the historical period is traced to excessive drawdown of aquifers to meet irrigation demand with a smaller contribution from thermal expansion. This shows how fully coupled models can provide indirect observational constraints upon forcing, in this case, constraining irrigation rates with observed sea level changes. The overall agreement of GISS‐E2.1 with observed trends is familiar from evaluation of its predecessors, as is the conclusion that these trends are almost entirely anthropogenic in origin.
Full text
Available for:
DOBA, FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, SIK, UILJ, UKNU, UL, UM, UPUK
This paper presents the response to anthropogenic forcing in the GISS-E2.1 climate models for the 21st century Shared Socioeconomic Pathways emission scenarios within the Coupled Model ...Intercomparison Project Phase 6 (CMIP6). The experiments were performed using an updated and improved version of the NASA Goddard Institute for Space Studies (GISS) coupled general circulation model that includes two different versions for atmospheric composition: A non-interactive version (NINT) with prescribed composition and a tuned aerosol indirect effect and the One-Moment Aerosol model (OMA) version with fully interactive aerosols which includes a parameterized first indirect aerosol effect on clouds. The effective climate sensitivities are 3.0°C and 2.9°C for the NINT and OMA models, respectively. Each atmospheric version is coupled to two different ocean general circulation models: The GISS ocean model (E2.1-G) and HYCOM (E2.1-H). We describe the global mean responses for all future scenarios and spatial patterns of change for surface air temperature and precipitation for four of the marker scenarios: SSP1-2.6, SSP2-4.5, SSP4-6.0, and SSP5-8.5. By 2100, global mean warming ranges from 1.5°C to 5.2°C relative to 1,850–1,880 mean temperature. Two high-mitigation scenarios SSP1-1.9 and SSP1-2.6 limit the surface warming to below 2°C by the end of the 21st century, except for the NINT E2.1-H model that simulates 2.2°C of surface warming. For the high emission scenario SSP5-8.5, the range is 4.6–5.2°C at 2100. Due to about 15% larger effective climate sensitivity and stronger transient climate response in both NINT and OMA CMIP6 models compared to CMIP5 versions, there is a stronger warming by 2100 in the SSP emission scenarios than in the comparable Representative Concentration Pathway (RCP) scenarios in CMIP5. Changes in sea ice area are highly correlated to global mean surface air temperature anomalies and show steep declines in both hemispheres, with the largest sea ice area decreases occurring during September in the Northern Hemisphere in both E2.1-G (−1.21 × 106 km2/°C) and E2.1-H models (−0.94 × 106 km2/°C). Both coupled models project decreases in the Atlantic overturning stream function by 2100. The largest decrease of 56%–65% in the 21st century overturning stream function is produced in the warmest scenario SSP5-8.5 in the E2.1-G model, comparable to the reduction in the corresponding CMIP5 GISS-E2 RCP8.5 simulation. Both low-end scenarios SSP1-1.9 and SSP1-2.6 also simulate substantial reductions of the overturning (9%–37%) with slow recovery of about 10% by the end of the 21st century (relative to the maximum decrease at the middle of the 21st century).
Full text
Available for:
DOBA, FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, SIK, UILJ, UKNU, UL, UM, UPUK
Large-eddy simulation (LES) is conducted of a three-dimensional temporal mixing layer whose lower stream is initially laden with liquid drops which may evaporate during the simulation. The gas-phase ...equations are written in an Eulerian frame for two perfect gas species (carrier gas and vapour emanating from the drops), while the liquid-phase equations are written in a Lagrangian frame. The effect of drop evaporation on the gas phase is considered through mass, species, momentum and energy source terms. The drop evolution is modelled using physical drops, or using computational drops to represent the physical drops. Simulations are performed using various LES models previously assessed on a database obtained from direct numerical simulations (DNS). These LES models are for: (i) the subgrid-scale (SGS) fluxes and (ii) the filtered source terms (FSTs) based on computational drops. The LES, which are compared to filtered-and-coarsened (FC) DNS results at the coarser LES grid, are conducted with 64 times fewer grid points than the DNS, and up to 64 times fewer computational than physical drops. It is found that both constant-coefficient and dynamic Smagorinsky SGS-flux models, though numerically stable, are overly dissipative and damp generated small-resolved-scale (SRS) turbulent structures. Although the global growth and mixing predictions of LES using Smagorinsky models are in good agreement with the FC-DNS, the spatial distributions of the drops differ significantly. In contrast, the constant-coefficient scale-similarity model and the dynamic gradient model perform well in predicting most flow features, with the latter model having the advantage of not requiring a priori calibration of the model coefficient. The ability of the dynamic models to determine the model coefficient during LES is found to be essential since the constant-coefficient gradient model, although more accurate than the Smagorinsky model, is not consistently numerically stable despite using DNS-calibrated coefficients. With accurate SGS-flux models, namely scale-similarity and dynamic gradient, the FST model allows up to a 32-fold reduction in computational drops compared to the number of physical drops, without degradation of accuracy; a 64-fold reduction leads to a slight decrease in accuracy.