Results of analysis and model simulation of the bed profiles

Bedforms such as dunes and ripples are ubiquitous in rivers and coastal seas, and commonly described as triangular shapes from which height and length are calculated to estimate hydrodynamic and sediment dynamic parameters. Natural bedforms, however, present a far more complicated morphology; the difference between natural bedform shape and the often assumed triangular shape is usually neglected, and how this may affect the flow is unknown. This study investigates the shapes of natural bedforms and how they influence flow and shear stress, based on four datasets extracted from earlier studies on two rivers (the Rio Paraná in Argentina, and the Lower Rhine in The Netherlands). The most commonly occurring morphological elements are a sinusoidal stoss side made of one segment and a lee side made of two segments, a gently sloping upper lee side and a relatively steep (6 to 21°) slip face. A non-hydrostatic numerical model, set up using Delft3D, served to simulate the flow over fixed bedforms with various morphologies derived from the identified morphological elements. Both shear stress and turbulence increase with increasing slip face angle and are only marginally affected by the dimensions and positions of the upper and lower lee side. The average slip face angle determined from the bed profiles is 14°, over which there is no permanent flow separation. Shear stress and turbulence above natural bedforms are higher than above a flat bed but much lower than over the often assumed 30° lee side angle.

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.

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

Stable isotope record and age determination of ODP Site 121-758 in the northeastern Indian Ocean

Numerous studies have shown that delta18O records from benthic and planktonic foraminifera, primarily a proxy of global ice volume variations, reflect Milankovitch periodicities. To study climatic response to orbital forcing at Ocean Drilling Program site 758, we have generated continuous delta18O and delta13C records from a single benthic foraminiferal species Cibicides wuellerstorfi for the last 3.6 m.y. and extended the planktonic foraminiferal isotope records of Farrell and Janecek (1991, doi:10.2973/odp.proc.sr.121.124.1991) (0-2.5 Ma, based on Globigerinoides sacculifer) to 3.6 Ma (Chen, 1994). We then constructed an age model by matching, correlating and tuning the benthic delta18O record to a model simulation of ice volume (Imbrie and Imbrie, 1980, doi:10.1126/science.207.4434.943). The filtered 41- and 23-kyr signals based on the resultant astronomically tuned age model are highly correlated to obliquity (r=0.83) and precession (r=0.75), respectively. Although derived with methodology different from Shackleton et al. (1990) and Hilgen (1991, doi:10.1016/0012-821X(91)90206-W, 1991, doi:10.1016/0012-821X(91)90082-S), our results generally agree with their published astronomical timescales for the time interval from 0 to 3.0 Ma, providing additional support for the newly emerging chronology based on orbital tuning. Slight discrepancies exist in the time interval from 3.0 to 3.6 Ma, suggesting several possibilities, including differences in the approaches of orbital tuning and the relatively low amplitude of delta18O variations in our record. However, even if the discrepancies are due to the relatively low amplitude of the isotope signals in our record at 3.0-3.6 Ma, our resultant timescale as a whole does not adversely affect our evaluation of the paleoclimatology and paleoceanography of the Indian Ocean, such as the evolution of the 100-, 41- and 23-kyr cycles, and variation of global ice volume and deepwater temperature during the past 3.6 m.y.

Field measurements of the atmosphere, ocean, sea ice and sub-ice platelet layer at Atka Bay in 2013

Ice shelves strongly interact with coastal Antarctic sea ice and the associated ecosystem by creating conditions favourable to the formation of a sub-ice platelet layer. The close investigation of this phenomenon and its seasonal evolution remain a challenge due to logistical constraints and a lack of suitable methodology. In this study, we characterize the seasonal cycle of Antarctic fast ice adjacent to the Ekström Ice Shelf in the eastern Weddell Sea. We used a thermistor chain with the additional ability to record the temperature response induced by cyclic heating of resistors embedded in the chain. Vertical sea-ice temperature and heating profiles obtained daily between November 2012 and February 2014 were analyzed to determine sea-ice and snow evolution, and to calculate the basal energy budget. The residual heat flux translated into an ice-volume fraction in the platelet layer of 0.18 ± 0.09, which we reproduced by a independent model simulation and agrees with earlier results. Manual drillings revealed an average annual platelet-layer thickness increase of at least 4m, and an annual maximum thickness of 10m beneath second-year sea ice. The oceanic contribution dominated the total sea-ice production during the study, effectively accounting for up to 70% of second-year sea-ice growth. In summer, an oceanic heat flux of 21 W/m**2 led to a partial thinning of the platelet layer. Our results further show that the active heating method, in contrast to the acoustic sounding approach, is well suited to derive the fast-ice mass balance in regions influenced by ocean/ice-shelf interaction, as it allows sub-diurnal monitoring of the platelet-layer thickness.

Oxygen Utilization and Downward Carbon Flux in an Oxygen-Depleted Eddy in the Eastern Tropical North Atlantic

The occurrence of mesoscale eddies that develop suboxic environments at shallow depth (about 40-100 m) has recently been reported for the eastern tropical North Atlantic (ETNA). Their hydrographic structure suggests that the water mass inside the eddy is well isolated from ambient waters supporting the development of severe near-surface oxygen deficits. So far, hydrographic and biogeochemical characterization of these eddies was limited to a few autonomous surveys, with the use of moorings, under water gliders and profiling floats. In this study we present results from the first dedicated biogeochemical survey of one of these eddies conducted in March 2014 near the Cape Verde Ocean Observatory (CVOO). During the survey the eddy core showed oxygen concentrations as low as 5 µmol kg-1 with a pH of around 7.6 at approximately 100 m depth. Correspondingly, the aragonite saturation level dropped to 1 at the same depth, thereby creating unfavorable conditions for calcifying organisms. To our knowledge, such enhanced acidity within near-surface waters has never been reported before for the open Atlantic Ocean. Vertical distributions of particulate organic matter and dissolved organic matter (POM and DOM), generally showed elevated concentrations in the surface mixed layer (0-70 m), with DOM also accumulating beneath the oxygen minimum. With the use of reference data from the upwelling region where these eddies are formed, the oxygen utilization rate was calculated by determining oxygen consumption through the remineralization of organic matter. Inside the core, we found these rates were almost 1 order of magnitude higher (apparent oxygen utilization rate (aOUR); 0.26 µmol kg-1 day-1) than typical values for the open North Atlantic. Computed downward fluxes for particulate organic carbon (POC), were around 0.19 to 0.23 g C m-2 day-1 at 100 m depth, clearly exceeding fluxes typical for an oligotrophic open-ocean setting. The observations support the view that the oxygen-depleted eddies can be viewed as isolated, westwards propagating upwelling systems of their own, thereby represent re-occurring alien biogeochemical environments in the ETNA.

Turbulent flux and meteorological measurements during ARCTEX-2006 campaign

This dataset present result from the DFG- funded Arctic-Turbulence-Experiment (ARCTEX-2006) performed by the University of Bayreuth on the island of Svalbard, Norway, during the winter/spring transition 2006. From May 5 to May 19, 2006 turbulent flux and meteorological measurements were performed on the monitoring field near Ny-Ålesund, at 78°55'24'' N, 11°55'15'' E Kongsfjord, Svalbard (Spitsbergen), Norway. The ARCTEX-2006 campaign site was located about 200 m southeast of the settlement on flat snow covered tundra, 11 m to 14 m above sea level. The permanent sites used for this study consisted of the 10 m meteorological tower of the Alfred Wegener Institute for Polar- and Marine Research (AWI), the international standardized radiation measurement site of the Baseline Surface Radiation Network (BSRN), the radiosonde launch site and the AWI tethered balloon launch sites. The temporary sites - set up by the University of Bayreuth - were a 6 m meteorological gradient tower, an eddy-flux measurement complex (EF), and a laser-scintillometer section (SLS). A quality assessment and data correction was applied to detect and eliminate specific measurement errors common at a high arctic landscape. In addition, the quality checked sensible heat flux measurements are compared with bulk aerodynamic formulas that are widely used in atmosphere-ocean/land-ice models for polar regions as described in Ebert and Curry (1993, doi:10.1029/93JC00656) and Launiainen and Cheng (1995). These parameterization approaches easily allow estimation of the turbulent surface fluxes from routine meteorological measurements. The data show: - the role of the intermittency of the turbulent atmospheric fluctuation of momentum and scalars, - the existence of a disturbed vertical temperature profile (sharp inversion layer) close to the surface, - the relevance of possible free convection events for the snow or ice melt in the Arctic spring at Svalbard, and - the relevance of meso-scale atmospheric circulation pattern and air-mass advection for the near-surface turbulent heat exchange in the Arctic spring at Svalbard. Recommendations and improvements regarding the interpretation of eddy-flux and laser-scintillometer data as well as the arrangement of the instrumentation under polar distinct exchange conditions and (extreme) weather situations could be derived.

Geochemistry of planktonic foraminifera from DSDP sites in the Atlantic and Pacific Ocean

The lithium content of planktonic foraminiferal calcite has been determined to evaluate temporal variability of seawater Li concentrations over the past 116 m.y. Mean foraminiferal calcite lithium/calcium in each time interval is no more than 16% greater nor 25% less than the mean Li/Ca of all samples. Li/Ca minima are observed for samples from 50-60 m.y. and 80-90 m.y., with Li/Ca about 25% lower than in adjacent time intervals. At no time during the past 40 m.y does mean Li/Ca appear to be higher than that at present. Subject to the limitations imposed by sample coverage and diagenesis, a similar conclusion holds for the past 116 m.y. Coupled with an oceanic mass balance model for Li, these data suggest that: (1) oceanic Li concentrations and, therefore, high-temperature hydrothermal circulation fluxes during the past 40 m.y. (and perhaps the past 100 m.y.) have not been more than perhaps 30-40% greater than at present for intervals any longer than a million years at most, and (2) these fluxes were not a factor of two higher 100 m.y. ago. By inference, variations in oceanic crustal generation rates over these time periods are similarly limited. Decreases in hydrothermal circulation fluxes and crustal generation rates or fluctuations up to 20% in these rates of a few million years duration are not necessarily ruled out by the Li/Ca data. The lack of variability in Li/Ca over time is not unequivocal evidence that hydrothermal fluxes have not varied because the rates of removal processes may be linked to changes in input fluxes.

Calcareous nannofossil datums and control points for the TCMS age model of ODP Site 175-1085

An astronomically calibrated age model for the Pliocene section of Ocean Drilling Program Leg 175 Cape Basin Site 1085 based on magnetic susceptibility data was developed using shipboard biostratigraphic datums. The composite core magnetic susceptibility record was compiled using shipboard correlations between Holes 1085A and 1085B and then tuned to the record of orbital variations in eccentricity to generate an orbitally tuned age model. Magnetic susceptibility apparently records climate variations in the Cape Basin. Strong power spectra values at the 100- and 400-k.y. frequency suggest an orbital control on the beat of Pliocene climate change in the Cape Basin.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Faunal analysis and Mg/Ca records of planktic foraminifera in sediment core SO189/2_039KL

Climate phenomena like the monsoon system, El Niño Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD) are interconnected via various feedback mechanisms and control the climate of the Indian Ocean and its surrounding continents on various timescales. The eastern tropical Indian Ocean is a key area for the interplay of these phenomena and for reconstructing their past changes and forcing mechanisms. Here we present records of upper ocean thermal gradient, thermocline temperatures (TT) and relative abundances of planktic foraminifera in core SO 189-39KL taken off western Sumatra (0°47.400' S, 99°54.510' E) for the last 8 ka that we use as proxies for changes in upper ocean structure. The records suggest a deeper thermocline between 8 ka and ca 3 ka compared to the late Holocene. We find a shoaling of the thermocline after 3 ka, most likely indicating an increased occurrence of upwelling during the late Holocene compared to the mid-Holocene which might represent changes in the IOD-like mean state of the Indian Ocean with a more negative IOD-like mean state during the mid-Holocene and a more positive IOD-like mean state during the past 3 ka. This interpretation is supported by a transient Holocene climate model simulation in which an IOD-like mode is identified that involves an insolation-forced long-term trend of increasing anomalous surface easterlies over the equatorial eastern Indian Ocean.

Static GPS measurements in the vicinity of EDML and an elevation model for Dronning Maud Land

The knowledge of ice sheet surface topography and the location of the ice divides are essential for ice dynamic modeling. An improved digital elevation model (DEM) of Dronning Maud Land (DML), Antarctica, is presented in this paper. It is based on ground-based kinematic GPS profiles, airborne radar altimetry, and data of the airborne radio-echo sounding system, as well as spaceborne laser altimetry from NASA's Ice, Cloud and land Elevation Satellite (ICESat). The accuracy of ICESat ice sheet altimetry data in the area of investigation is discussed. The location of the ice divides is derived from aspect calculation of the topography and is verified with several velocity data derived from repeated static GPS measurements.

Model parameter

Interannual environmental variability in Peru is dominated by the El Niño Southern Oscillation (ENSO). The most dramatic changes are associated with the warm El Niño (EN) phase (opposite the cold La Niña phase), which disrupts the normal coastal upwelling and affects the dynamics of many coastal marine and terrestrial resources. This study presents a trophic model for Sechura Bay, located at the northern extension of the Peruvian upwelling system, where ENSO-induced environmental variability is most extreme. Using an initial steady-state model for the year 1996, we explore the dynamics of the ecosystem through the year 2003 (including the strong EN of 1997/98 and the weaker EN of 2002/03). Based on support from literature, we force biomass of several non-trophically-mediated 'drivers' (e.g. Scallops, Benthic detritivores, Octopus, and Littoral fish) to observe whether the fit between historical and simulated changes (by the trophic model) is improved. The results indicate that the Sechura Bay Ecosystem is a relatively inefficient system from a community energetics point of view, likely due to the periodic perturbations of ENSO. A combination of high system productivity and low trophic level target species of invertebrates (i.e. scallops) and fish (i.e. anchoveta) results in high catches and an efficient fishery. The importance of environmental drivers is suggested, given the relatively small improvements in the fit of the simulation with the addition of trophic drivers on remaining functional groups' dynamics. An additional multivariate regression model is presented for the scallop Argopecten purpuratus, which demonstrates a significant correlation between both spawning stock size and riverine discharge-mediated mortality on catch levels. These results are discussed in the context of the appropriateness of trophodynamic modeling in relatively open systems, and how management strategies may be focused given the highly environmentally influenced marine resources of the region.