Three time-slice simulations of the Eemian (eem), the mid-Holocene (hol), and a pre-industrial (pri) control run with the coupled atmosphere-ocean-model ECHAM5/MPI-OM

Orbital forcing does not only exert direct insolation effects, but also alters climate indirectly through feedback mechanisms that modify atmosphere and ocean dynamics and meridional heat and moisture transfers. We investigate the regional effects of these changes by detailed analysis of atmosphere and ocean circulation and heat transports in a coupled atmosphere-ocean-sea ice-biosphere general circulation model (ECHAM5/JSBACH/MPI-OM). We perform long term quasi equilibrium simulations under pre-industrial, mid-Holocene (6000 years before present - yBP), and Eemian (125 000 yBP) orbital boundary conditions. Compared to pre-industrial climate, Eemian and Holocene temperatures show generally warmer conditions at higher and cooler conditions at lower latitudes. Changes in sea-ice cover, ocean heat transports, and atmospheric circulation patterns lead to pronounced regional heterogeneity. Over Europe, the warming is most pronounced over the north-eastern part in accordance with recent reconstructions for the Holocene. We attribute this warming to enhanced ocean circulation in the Nordic Seas and enhanced ocean-atmosphere heat flux over the Barents Shelf in conduction with retreat of sea ice and intensified winter storm tracks over northern Europe.

Results of a Finite-Element Sea-Ice Ocean Model (FESOM) set-up to study the interannual to decadal variability in the deep-water formation rates

The characteristics of a global set-up of the Finite-Element Sea-Ice Ocean Model under forcing of the period 1958-2004 are presented. The model set-up is designed to study the variability in the deep-water mass formation areas and was therefore regionally better resolved in the deep-water formation areas in the Labrador Sea, Greenland Sea, Weddell Sea and Ross Sea. The sea-ice model reproduces realistic sea-ice distributions and variabilities in the sea-ice extent of both hemispheres as well as sea-ice transport that compares well with observational data. Based on a comparison between model and ocean weather ship data in the North Atlantic, we observe that the vertical structure is well captured in areas with a high resolution. In our model set-up, we are able to simulate decadal ocean variability including several salinity anomaly events and corresponding fingerprint in the vertical hydrography. The ocean state of the model set-up features pronounced variability in the Atlantic Meridional Overturning Circulation as well as the associated mixed layer depth pattern in the North Atlantic deep-water formation areas.

On the effects of circulation, sediment resuspension and biological incorporation by diatoms in an ocean model of aluminium, link to model data

The distribution of dissolved aluminium in the West Atlantic Ocean shows a mirror image with that of dissolved silicic acid, hinting at intricate interactions between the ocean cycling of Al and Si. The marine biogeochemistry of Al is of interest because of its potential impact on diatom opal remineralisation, hence Si availability. Furthermore, the dissolved Al concentration at the surface ocean has been used as a tracer for dust input, dust being the most important source of the bio-essential trace element iron to the ocean. Previously, the dissolved concentration of Al was simulated reasonably well with only a dust source, and scavenging by adsorption on settling biogenic debris as the only removal process. Here we explore the impacts of (i) a sediment source of Al in the Northern Hemisphere (especially north of ~ 40° N), (ii) the imposed velocity field, and (iii) biological incorporation of Al on the modelled Al distribution in the ocean. The sediment source clearly improves the model results, and using a different velocity field shows the importance of advection on the simulated Al distribution. Biological incorporation appears to be a potentially important removal process. However, conclusive independent data to constrain the Al / Si incorporation ratio by growing diatoms are missing. Therefore, this study does not provide a definitive answer to the question of the relative importance of Al removal by incorporation compared to removal by adsorptive scavenging.

Regional ocean model for the Mediterranean, links to model files

A regional ocean general circulation model of the Mediterranean is used to study the climate of the Last Glacial Maximum. The atmospheric forcing for these simulations has been derived from simulations with an atmospheric general circulation model, which in turn was forced with surface conditions from a coarse resolution earth system model. The model is successful in reproducing the general patterns of reconstructed sea surface temperature anomalies with the strongest cooling in summer in the northwestern Mediterranean and weak cooling in the Levantine, although the model underestimates the extent of the summer cooling in the western Mediterranean. However, there is a strong vertical gradient associated with this pattern of summer cooling, which makes the comparison with reconstructions complicated. The exchange with the Atlantic is decreased to roughly one half of its present value, which can be explained by the shallower Strait of Gibraltar as a consequence of lower global sea level. This reduced exchange causes a strong increase of salinity in the Mediterranean in spite of reduced net evaporation.

Holocene experiment with coupled atmosphere-ocean-model ECHAM5/MPI-OM

Changes in the Earth's orbit lead to changes in the seasonal and meridional distribution of insolation. We quantify the influence of orbitally induced changes on the seasonal temperature cycle in a transient simulation of the last 6000 years - from the mid-Holocene to today - using a coupled atmosphere-ocean general circulation model (ECHAM5/MPI-OM) including a land surface model (JSBACH). The seasonal temperature cycle responds directly to the insolation changes almost everywhere. In the Northern Hemisphere, its amplitude decreases according to an increase in winter insolation and a decrease in summer insolation. In the Southern Hemisphere, the opposite is true. Over the Arctic Ocean, decreasing summer insolation leads to an increase in sea-ice cover. The insulating effect of sea ice between the ocean and the atmosphere leads to decreasing heat flux and favors more "continental" conditions over the Arctic Ocean in winter, resulting in strongly decreasing temperatures. Consequently, there are two competing effects: the direct response to insolation changes and a sea-ice insulation effect. The sea-ice insulation effect is stronger, and thus an increase in the amplitude of the seasonal temperature cycle over the Arctic Ocean occurs. This increase is strongest over the Barents Shelf and influences the temperature response over northern Europe. We compare our modeled seasonal temperatures over Europe to paleo reconstructions. We find better agreements in winter temperatures than in summer temperatures and better agreements in northern Europe than in southern Europe, since the model does not reproduce the southern European Holocene summer cooling inferred from the paleo reconstructions. The temperature reconstructions for northern Europe support the notion of the influence of the sea-ice insulation effect on the evolution of the seasonal temperature cycle.

Simulated Methyl iodide emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by OMIP data

Production pathways of the prominent volatile organic halogen compound methyl iodide (CH3I) are not fully understood. Based on observations, production of CH3I via photochemical degradation of organic material or via phytoplankton production has been proposed. Additional insights could not be gained from correlations between observed biological and environmental variables or from biogeochemical modeling to identify unambiguously the source of methyl iodide. In this study, we aim to address this question of source mechanisms with a three-dimensional global ocean general circulation model including biogeochemistry (MPIOM-HAMOCC (MPIOM - Max Planck Institute Ocean Model HAMOCC - HAMburg Ocean Carbon Cycle model)) by carrying out a series of sensitivity experiments. The simulated fields are compared with a newly available global data set. Simulated distribution patterns and emissions of CH3I differ largely for the two different production pathways. The evaluation of our model results with observations shows that, on the global scale, observed surface concentrations of CH3I can be best explained by the photochemical production pathway. Our results further emphasize that correlations between CH3I and abiotic or biotic factors do not necessarily provide meaningful insights concerning the source of origin. Overall, we find a net global annual CH3I air-sea flux that ranges between 70 and 260 Gg/yr. On the global scale, the ocean acts as a net source of methyl iodide for the atmosphere, though in some regions in boreal winter, fluxes are of the opposite direction (from the atmosphere to the ocean).

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

The CMCC Eddy-permitting Global Ocean Physical Reanalysis (C-GLORS v5, 1980-2014), links to NetCDF files

The CMCC Global Ocean Physical Reanalysis System (C-GLORS) is used to simulate the state of the ocean in the last decades. It consists of a variational data assimilation system (OceanVar), capable of assimilating all in-situ observations along with altimetry data, and a forecast step performed by the ocean model NEMO coupled with the LIM2 sea-ice model. KEY STRENGTHS: - Data are available for a large number of ocean parameters - An extensive validation has been conducted and is freely available - The reanalysis is performed at high resolution (1/4 degree) and spans the last 30 years KEY LIMITATIONS: - Quality may be discontinuos and depend on observation coverage - Uncertainty estimates are simply derived through verification skill scores

FESOM Arctic Ocean sea ice concentration and thickness 1995-2014, links to movies in mp4 format

Northern Hemisphere sea ice from a Finite-Element Sea-Ice Ocean Model (FESOM) 4.5 km resolution simulation carried out by researchers from Alfred Wegener Institute (AWI), Germany. Concentration is shown with color; thickness is shown with shading. A global 1 degree mesh is used, with the "Arctic Ocean" locally refined to 4.5 km. South of CAA and Fram Strait the resolution is not refined in this simulation. The animation indicates that the 4.5 km model resolution helps to represent the small scale sea ice features, although much higher resolution is required to fully resolve the ice leads. The animation is created by Michael Böttinger from DKRZ (https://www.dkrz.de).

Towards convection-resolving, global atmospheric simulations with the Model for Prediction Across Scales (MPAS) v3.1: an extreme scaling experiment, links to sourcecode of model and grides

The Model for Prediction Across Scales (MPAS) is a novel set of Earth system simulation components and consists of an atmospheric model, an ocean model and a land-ice model. Its distinct features are the use of unstructured Voronoi meshes and C-grid discretisation to address shortcomings of global models on regular grids and the use of limited area models nested in a forcing data set, with respect to parallel scalability, numerical accuracy and physical consistency. This concept allows one to include the feedback of regional land use information on weather and climate at local and global scales in a consistent way, which is impossible to achieve with traditional limited area modelling approaches. Here, we present an in-depth evaluation of MPAS with regards to technical aspects of performing model runs and scalability for three medium-size meshes on four different high-performance computing (HPC) sites with different architectures and compilers. We uncover model limitations and identify new aspects for the model optimisation that are introduced by the use of unstructured Voronoi meshes. We further demonstrate the model performance of MPAS in terms of its capability to reproduce the dynamics of the West African monsoon (WAM) and its associated precipitation in a pilot study. Constrained by available computational resources, we compare 11-month runs for two meshes with observations and a reference simulation from the Weather Research and Forecasting (WRF) model. We show that MPAS can reproduce the atmospheric dynamics on global and local scales in this experiment, but identify a precipitation excess for the West African region. Finally, we conduct extreme scaling tests on a global 3?km mesh with more than 65 million horizontal grid cells on up to half a million cores. We discuss necessary modifications of the model code to improve its parallel performance in general and specific to the HPC environment. We confirm good scaling (70?% parallel efficiency or better) of the MPAS model and provide numbers on the computational requirements for experiments with the 3?km mesh. In doing so, we show that global, convection-resolving atmospheric simulations with MPAS are within reach of current and next generations of high-end computing facilities.

Physical oceanography during Maria S. Merian cruise MSM21/2

Flemish Pass, located at the western subpolar margin, is a passage (sill depth 1200 m) that is constrained by the Grand Banks and the underwater plateau Flemish Cap. In addition to the Deep Western Boundary Current (DWBC) pathway offshore of Flemish Cap, Flemish Pass represents another southward transport pathway for two modes of Labrador Sea Water (LSW), the lightest component of North Atlantic Deep Water carried with the DWBC. This pathway avoids potential stirring regions east of Flemish Cap and deflection into the interior North Atlantic. Ship-based velocity measurements between 2009 and 2013 at 47°N in Flemish Pass and in the DWBC east of Flemish Cap revealed a considerable southward transport of Upper LSW through Flemish Pass (15-27%, -1.0 to -1.5 Sv). About 98% of the denser Deep LSW were carried around Flemish Cap as Flemish Pass is too shallow for considerable transport of Deep LSW. Hydrographic time series from ship-based measurements show a significant warming of 0.3°C/decade and a salinification of 0.03/decade of the Upper LSW in Flemish Pass between 1993 and 2013. Almost identical trends were found for the evolution in the Labrador Sea and in the DWBC east of Flemish Cap. This indicates that the long-term hydrographic variability of Upper LSW in Flemish Pass as well as in the DWBC at 47°N is dominated by changes in the Labrador Sea, which are advected southward. Fifty years of numerical ocean model simulations in Flemish Pass suggest that these trends are part of a multidecadal cycle.

Organic geochemistry of sediment cores from the Arbian Sea

Two deep-sea sediment cores from the northeastern and the southeastern Arabian Sea were studied in order to reconstruct the palaeoenvironments of the past glacial cycles. Core 136KL was recovered from the high-productivity area off Pakistan within the modern oxygen-minimum zone (OMZ). By contrast, modern primary productivity at the site of MD900963 close to Maldives is moderate and bottom waters are today well oxygenated. For both cores, we reconstructed the changes in palaeoproductivity using a set of biomarkers (alkenones, dinosterol and brassicasterol); the main result is that primary productivity is enhanced during glacial stages and lowered during interstadials. The proxies associated with productivity show a 23 kyr cyclicity corresponding to the precession-related insolation cycle. Palaeoredox conditions were studied in both cores using a new organic geochemical parameter (C35/C31-n-alkane ratio) developed by analysing surface sediments from a transect across the OMZ off Pakistan. The value of this ratio in core 136KL shows many variations during the last 65 kyr, indicating that the OMZ was not stable during this time: it disappeared completely during Heinrich- and the Younger Dryas events, pointing to a connection between global oceanic circulation and the stability of the OMZ. The C35/C31 ratio determined in sediments of core MD900963 shows that bottom waters remained rather well oxygenated over the last 330 kyr, which is confirmed by comparison with authigenic metal concentrations in the same sediments. A zonally averaged, circulation-biogeochemical ocean model was used to explore how the intermediate Indian Ocean responds to a freshwater flux anomaly at the surface of the North Atlantic. As suggested by the geochemical time series, both the abundance of Southern Ocean Water and the oxygen concentration are significantly increased in response to this freshwater perturbation.