ALBERT

Description: 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 degrees S. Over the 1933-2007 period, the CORE-II models consistently simulate spatially variable trends in sea-ice concentration, surface freshwater fluxes, mixed layer depths, and 200-700 in 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

Description: Highlights: • 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. Abstract: 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.

Description: Highlights: • 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.

Description: Highlights: • 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.

Description: Highlights: • We combine high-resolution ocean models with population genetics • Variation in wind-driven ocean currents mediates the collapse of A. anguilla • Female eels are philopatric within the Sargasso Sea, while males maintain gene flow • We present first evidence of the role of ocean currents in shaping species’ evolution Summary: Worldwide, exploited marine fish stocks are under threat of collapse [1]. Although the drivers behind such collapses are diverse, it is becoming evident that failure to consider evolutionary processes in fisheries management can have drastic consequences on a species’ long-term viability [2]. The European eel (Anguilla anguilla; Linnaeus, 1758) is no exception: not only does the steep decline in recruitment observed in the 1980s [ 3 and 4] remain largely unexplained, the punctual detection of genetic structure also raises questions regarding the existence of a single panmictic population [ 5, 6 and 7]. With its extended Transatlantic dispersal, pinpointing the role of ocean dynamics is crucial to understand both the population structure and the widespread decline of this species. Hence, we combined dispersal simulations using a half century of high-resolution ocean model data with population genetics tools. We show that regional atmospherically driven ocean current variations in the Sargasso Sea were the major driver of the onset of the sharp decline in eel recruitment in the beginning of the 1980s. The simulations combined with genotyping of natural coastal eel populations furthermore suggest that unexpected evidence of coastal genetic differentiation is consistent with cryptic female philopatric behavior within the Sargasso Sea. Such results demonstrate the key constraint of the variable oceanic environment on the European eel population.

Description: Highlights: • Global mean sea level simulated in interannual CORE simulations. • Regional sea level patterns simulated in interannual CORE simulations. • Theoretical foundation for analysis of global mean sea level and regional patterns. Abstract: We provide an assessment of sea level simulated in a suite of global ocean-sea ice models using the interannual CORE atmospheric state to determine surface ocean boundary buoyancy and momentum fluxes. These CORE-II simulations are compared amongst themselves as well as to observation-based estimates. We focus on the final 15 years of the simulations (1993–2007), as this is a period where the CORE-II atmospheric state is well sampled, and it allows us to compare sea level related fields to both satellite and in situ analyses. The ensemble mean of the CORE-II simulations broadly agree with various global and regional observation-based analyses during this period, though with the global mean thermosteric sea level rise biased low relative to observation-based analyses. The simulations reveal a positive trend in dynamic sea level in the west Pacific and negative trend in the east, with this trend arising from wind shifts and regional changes in upper 700 m ocean heat content. The models also exhibit a thermosteric sea level rise in the subpolar North Atlantic associated with a transition around 1995/1996 of the North Atlantic Oscillation to its negative phase, and the advection of warm subtropical waters into the subpolar gyre. Sea level trends are predominantly associated with steric trends, with thermosteric effects generally far larger than halosteric effects, except in the Arctic and North Atlantic. There is a general anti-correlation between thermosteric and halosteric effects for much of the World Ocean, associated with density compensated changes.

Description: Highlights: • A joint analysis of deep current meter records in the western North Atlantic. • Intra-seasonal variability dominates the deep boundary current. • Topographic waves near 10d periods trapped over steep topography. • Basin centers are showing longer periods (50d) caused by the eddy field. • Observed variability characteristics compared to high resolution model simulation. Abstract The Deep Western Boundary Current (DWBC) along the western margin of the subpolar North Atlantic is an important component of the deep limb of the Meridional Overturning near its northern origins. A network of moored arrays from Denmark Strait to the tail of the Grand Banks has been installed for almost two decades to observe the boundary currents and transports of North Atlantic Deep Water as part of an internationally coordinated observatory for the Atlantic Meridional Overturning Circulation. The dominant variability in all of the moored velocity time series is in the week-to-month period range. While the temporal characteristics of this variability change only gradually between Denmark Strait and Flemish Cap, a broad band of longer term variability is present farther along the path of the DWBC at the Grand Banks and in the interior basins (Labrador and Irminger Seas). The vigorous intra-seasonal variability may well mask possible interannual to decadal variability that is typically an order of magnitude smaller than the high-frequency fluctuations. Here, the intra-seasonal variability is quantified at key positions along the DWBC path using both, observations and high resolution model data. The results are used to evaluate the model circulation, and in turn the model is used to relate the discrete measurements to the overall pattern of the subpolar circulation. Topographic waves are found to be trapped by the steep topography all around the western basins, the Labrador and Irminger Seas. In the Labrador Sea, the high intra-seasonal variability of the boundary current regime is separated by a region of extremely low variability in narrow recirculation cells from the basin interior. There, the variability is also on intra-seasonal timescales, but at much longer periods around 50 days.

Description: Highlights: • OGCM simulations of the AMOC are highly sensitive to the subarctic freshwater forcing. • Trends in the simulated AMOC are linked to the salinity of the DSOW. • DSOW salinity trends can be traced back to the freshwater transport by the NAC. • The NAC freshwater budget is highly affected by the salinity restoring used in OGCMs. • Modifications in the subarctic precipitation can help to minimize the restoring flux. Global ocean sea-ice models with an atmospheric forcing based on bulk formulations of the air-sea fluxes exhibit spurious trends in key flow indices like the Atlantic Meridional Overturning Circulation (AMOC), constraining their use in investigations of multi-decadal ocean variability. To identify the critical model factors affecting the temporal evolution of the AMOC on time scales of up to 60 years, a series of experiments with both eddy-permitting (0.25°) and non-eddying (0.5°) ocean-ice models has been performed, focusing on the influence of artificial choices for the freshwater forcing, in particular the restoring of sea surface salinity towards climatological values. The atmospheric forcing builds on the proposal for Coordinated Ocean-ice Reference Experiments (CORE), utilizing the refined atmospheric reanalysis products for 1948–2006 compiled by Large and Yeager. Sensitivity experiments with small variations in precipitation (within the observational uncertainty) and sea surface salinity restoring in the subarctic Atlantic produce a wide range of AMOC transports, between upward drifts to more than 22 Sv and nearly-collapsed states with less than 7 Sv, reflecting the excessive role of the salinity feedback in such simulations. In all cases the AMOC is tightly related to the density of the Denmark Strait overflow; changes in that density are governed by the salinity in the Nordic Seas; and in turn, that salinity is strongly affected by the properties of the inflowing North Atlantic water.

Description: The GEOTRACES Intermediate Data Product 2017 (IDP2017) is the second publicly available data product of the international GEOTRACES programme, and contains data measured and quality controlled before the end of 2016. The IDP2017 includes data from the Atlantic, Pacific, Arctic, Southern and Indian oceans, with about twice the data volume of the previous IDP2014. For the first time, the IDP2017 contains data for a large suite of biogeochemical parameters as well as aerosol and rain data characterising atmospheric trace element and isotope (TEI) sources. The TEI data in the IDP2017 are quality controlled by careful assessment of intercalibration results and multi-laboratory data comparisons at crossover stations. The IDP2017 consists of two parts: (1) a compilation of digital data for more than 450 TEIs as well as standard hydrographic parameters, and (2) the eGEOTRACES Electronic Atlas providing an on-line atlas that includes more than 590 section plots and 130 animated 3D scenes. The digital data are provided in several formats, including ASCII, Excel spreadsheet, netCDF, and Ocean Data View collection. Users can download the full data packages or make their own custom selections with a new on-line data extraction service. In addition to the actual data values, the IDP2017 also contains data quality flags and 1-σ data error values where available. Quality flags and error values are useful for data filtering and for statistical analysis. Metadata about data originators, analytical methods and original publications related to the data are linked in an easily accessible way. The eGEOTRACES Electronic Atlas is the visual representation of the IDP2017 as section plots and rotating 3D scenes. The basin-wide 3D scenes combine data from many cruises and provide quick overviews of large-scale tracer distributions. These 3D scenes provide geographical and bathymetric context that is crucial for the interpretation and assessment of tracer plumes near ocean margins or along ridges. The IDP2017 is the result of a truly international effort involving 326 researchers from 22 countries. This publication provides the critical reference for unpublished data, as well as for studies that make use of a large cross-section of data from the IDP2017.

Description: Highlights • Method and application to improve digital soil maps of silt and clay in China • Within the framework of a DSM approach we derived spatial uncertainties. • Spatial uncertainty is based on randomized decision trees. • Model calibration set is refined by purposive sampling in area of high uncertainty. • Method and map refinement is approved using accuracy and uncertainty measures. Digital soil mapping (DSM) products represent estimates of spatially distributed soil properties. These estimations comprise an element of uncertainty that is not evenly distributed over the area covered by DSM. If we quantify the uncertainty spatially explicit, this information can be used to improve the quality of DSM by optimizing the sampling design. This study follows a DSM approach using a Random Forest regression model, legacy soil samples, and terrain covariates to estimate topsoil silt and clay contents in a small catchment of 4.2 km2 in the Three Gorges Reservoir Area, Central China. We aim (i) to introduce a method to derive spatial uncertainty, and (ii) to improve the initial DSM approaches by additional sampling that is guided by the spatial uncertainty. The proposed uncertainty measure is based on multiple realizations of individual and randomized decision tree models. We used the spatial uncertainty of the initial DSM approaches to stratify the study area and thereby to identify potential sampling areas of high uncertainties. Further, we tested how precisely available legacy samples cover the variability of the covariates within each potential sampling area to define the final sampling area and to apply a purposive sampling design. For the final Random Forest model calibration, we combined the legacy sample set with the additional samples. This uncertainty-driven DSM refinement was evaluated by comparing it to a second approach. In this second approach, the additional samples were replaced by a random sample set of the same size, obtained from the entire study area. For the comparative analysis, external, bootstrap-, and cross-validation was applied. The DSM approach using the uncertainty-driven refinement performed best. The averaged spatial uncertainty was reduced by 31% for silt and by 27% for clay compared to the initial DSM approach. Using external validation, the accuracy increased by the same proportions, while showing an overall accuracy of R2 = 0.59 for silt and R2 = 0.56 for clay.