Gas emissions and ROV survey tracks of the Regab pockmark (Northern Congo Fan)

Pockmarks are seafloor depressions commonly associated with fluid escape from the seabed and are believed to contribute noticeably to the transfer of methane into the ocean and ultimately into the atmosphere. They occur in many different areas and geological contexts, and vary greatly in size and shape. Nevertheless, the mechanisms of pockmark growth are still largely unclear. Still, seabed methane emissions contribute to the global carbon budget, and understanding such processes is critical to constrain future quantifications of seabed methane release at local and global scales. The giant Regab pockmark (9°42.6' E, 5°47.8' S), located at 3160 m water depth near the Congo deep-sea channel (offshore southwestern Africa), was investigated with state-of-the-art mapping devices mounted on IFREMER's (French Research Institute for Exploitation of the Sea) remotely operated vehicle (ROV) Victor 6000. ROV-borne micro-bathymetry and backscatter data of the entire structure, a high-resolution photo-mosaic covering 105,000 m2 of the most active area, sidescan mapping of gas emissions, and maps of faunal distribution as well as of carbonate crust occurrence are combined to provide an unprecedented detailed view of a giant pockmark. All data sets suggest that the pockmark is composed of two very distinctive zones in terms of seepage intensity. We postulate that these zones are the surface expression of two fluid flow regimes in the subsurface: focused flow through a fractured medium and diffuse flow through a porous medium. We conclude that the growth of giant pockmarks is controlled by self-sealing processes and lateral spreading of rising fluids. In particular, partial redirection of fluids through fractures in the sediments can drive the pockmark growth in preferential directions.

Effects of ocean acidification and global warming on reef bioerosion-lessons from a clionaid sponge

Coral reefs are under threat, exerted by a number of interacting effects inherent to the present climate change, including ocean acidification and global warming. Bioerosion drives reef degradation by recycling carbonate skeletal material and is an important but understudied factor in this context. Twelve different combinations of pCO2 and temperature were applied to elucidate the consequences of ocean acidification and global warming on the physiological response and bioerosion rates of the zooxanthellate sponge Cliona orientalis-one of the most abundant and effective bioeroders on the Great Barrier Reef, Australia. Our results confirm a significant amplification of the sponges' bioerosion capacity with increasing pCO2, which is expressed by more carbonate being chemically dissolved by etching. The health of the sponges and their photosymbionts was not affected by changes in pCO2, in contrast to temperature, which had significant negative impacts at higher levels. However, we could not conclusively explain the relationship between temperature and bioerosion rates, which were slightly reduced at both colder as well as warmer temperatures than ambient. The present findings on the effects of ocean acidification on chemical bioerosion, however, will have significant implications for predicting future reef carbonate budgets, as sponges often contribute the lion's share of internal bioerosion on coral reefs.

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.

Benthic foraminiferal stable isotope record from ODP Site 138-849 on the East Pacific Rise

Benthic foraminifer and delta13C data from Site 849, on the west flank of the East Pacific Rise (0°11 'N, 110°31'W; 3851 m), give relatively continuous records of deep Pacific Ocean stable isotope variations between 0 and 5 Ma. The mean sample spacing is 4 k.y. Most analyses are from Cibicides wuellerstorfi, but isotopic offsets relative to Uvigerina peregrina appear roughly constant. Because of its location west of the East Pacific Rise, Site 849 yields a suitable record of mean Pacific Ocean delta13C, which approximates a global oceanic signal. The ~100-k.y.-period climate cycle, which is prevalent in delta18O does not dominate the long-term delta13C record. For delta13C, variations in the ~400- and 41-k.y. periods are more important. Phase lags of delta13C relative to ice volume in the 41- and 23-k.y. bands are consistent with delta13C as a measure of organic biomass. A model-calculated exponential response time of 1-2 k.y. is appropriate for carbon stored in soils and shallow sediments responding to glacial-interglacial climate change. Oceanic delta13C leads ice volume slightly in the 100-k.y. band, and this suggests another process such as changes in continental weathering to modulate mean river delta13C at long periods. The delta13C record from Site 849 diverges from that of Site 677 in the Panama Basin mostly because of decay of 13C-depleted organic carbon in the relatively isolated Panama Basin. North Atlantic to Pacific delta13C differences calculated using published data from Sites 607 and 849 reveal variations in Pliocene deep water within the range of those of the late Quaternary. Maximum delta13C contrast between these sites, which presumably reflects maximum influx of high-delta13C northern source water into the deep North Atlantic Ocean, occurred between 1.3 and 2.1 Ma, well after the initiation of Northern Hemisphere glaciation. Export of high-delta13C North Atlantic Deep Water from the Atlantic to the circumpolar Antarctic, as recorded by published delta13C data from Subantarctic Site 704, appears unrelated to the North Atlantic-Pacific delta13C contrast. To account for this observation, we suggest that deep-water formation in the North Atlantic reflects northern source characteristics, whereas export of this water into the circumpolar Antarctic reflects Southern Hemisphere wind forcing. Neither process appears directly linked to ice-volume variations.

Magnetobiostratigraphy and stable isotope record of Paleocene sediments from ODP Hole 122-761B in the subtropical Indian Ocean

Surface and deep water circulation patterns in the eastern Indian Ocean during the Paleocene Epoch are inferred based on an integrated magnetobiostratigraphic and stable isotope investigation of Ocean Drilling Program Hole 761B, drilled on the Wombat Plateau. A combination of magnetostratigraphy, biostratigraphy and isotope stratigraphy demonstrates that numerous deep sea sites that have been considered to show continuous, or nearly continuous sedimentation through the Paleocene are punctuated by a series of hiatuses, some of which exceeding a duration of 1 Myr. Therefore, our study is based on a detailed temporal interpretation of the stratigraphic successions we used for paleoceanographic reconstructions. We compare detailed planktonic and benthic foraminiferal carbon and oxygen isotope records from Hole 761B with several temporally correlative records published from different oceanic provinces in order to distinguish between local and global patterns within the eastern Indian Ocean. Although Site 761 was situated at low latitudes during the Paleocene, its surface waters were predominantly influenced by circulation originating from the Southern Ocean as indicated by inferred cool sea surface temperatures and reduced surface to deep water temperature gradients. We suggest that deep waters in the eastern Indian Ocean were not directly fed by the Southern or Tethys Oceans. Rather, the more negative delta13C composition of the bottom waters recorded by benthic foraminifera implies the presence and/or active contribution of aged deep waters from the Pacific during this time, at least prior to ~60.2 Ma and subsequent to ~59.0 Ma. The Indian continent, Ninetyeast Ridge, Kerguelen Plateau and Broken Ridge may have played a significant role as submarine barriers to deep water circulation during the Paleocene.

Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.

Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo 1.0.5-beta)

Sub-ice shelf circulation and freezing/melting rates in ocean general circulation models depend critically on an accurate and consistent representation of cavity geometry. Existing global or pan-Antarctic data sets have turned out to contain various inconsistencies and inaccuracies. The goal of this work is to compile independent regional fields into a global data set. We use the S-2004 global 1-minute bathymetry as the backbone and add an improved version of the BEDMAP topography for an area that roughly coincides with the Antarctic continental shelf. Locations of the merging line have been carefully adjusted in order to get the best out of each data set. High-resolution gridded data for upper and lower ice surface topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier have been carefully merged into the ambient ice and ocean topographies. Multibeam survey data for bathymetry in the former Larsen B cavity and the southeastern Bellingshausen Sea have been obtained from the data centers of Alfred Wegener Institute (AWI), British Antarctic Survey (BAS) and Lamont-Doherty Earth Observatory (LDEO), gridded, and again carefully merged into the existing bathymetry map. The global 1-minute dataset (RTopo-1 Version 1.0.5) has been split into two netCDF files. The first contains digital maps for global bedrock topography, ice bottom topography, and surface elevation. The second contains the auxiliary maps for data sources and the surface type mask. A regional subset that covers all variables for the region south of 50 deg S is also available in netCDF format. Datasets for the locations of grounding and coast lines are provided in ASCII format.

Osmium isotope composition of ODP Hole 121-758A

This study presents osmium (Os) isotope and elemental data for cleaned planktic foraminifera, authigenic Fe-Mn oxyhydroxides and pelagic carbonate host sediments from ODP site 758 in the southernmost reaches of the Bay of Bengal. The Os in the bulk sediments appears to be dominantly hydrogeneous (sourced by carbonate and Fe-Mn oxyhydroxide), but variations in this particular core are controlled by the presence of volcanic ash. Fe-Mn oxyhydroxide leachates (of the bulk sediments) from Holocene samples also yield an Os isotope composition close to that of seawater, but the record diverges from that of foraminifera at a depth corresponding to the oxic/post-oxic boundary, suggesting diagenetic mobilization of Os at depths below this. Holocene planktic foraminifera, cleaned using oxidative-reductive techniques, also give Os isotope compositions indistinguishable from modern seawater, but the record obtained for the past 150 kyr shows strong covaraitions of 187Os/188Os with both the local and global oxygen isotope record, with less radiogenic Os isotope compositions during glacial intervals. These results indicate that foraminifera provide a robust record of seawater Os isotope compositions, and comparison of the data obtained here with records from the other major oceans demonstrate global changes in 187Os/188Os over this time interval, while the covariation with oxygen isotopes suggest a process controlling the Os isotope composition that is in phase with global climate cycles. Global excursions to relatively unradiogenic 187Os/188Os during glacial intervals are consistent with decreased input of radiogenic continental material, reflecting cooler temperatures and reduced continental runoff. Modelling indicates that the shift to unradiogenic values during glacial intervals could be caused by an ~30% decrease in the global river flux, with an ~5% change in river composition. If the residence time of Os in the oceans is ~5 ka then the post-glacial recovery to present-day seawater values is consistent with a corresponding increase in the river flux of around 30%. However, if the residence time of Os is closer to 40 ka, as is suggested by the global river flux, then this demands either significant changes in both the riverine Os flux and composition of around 40% and 30%, respectively, that closely follow the oxygen isotope record, or else a short-lived post-glacial pulse of weathering some 75% greater than the steady-state flux. In either case, these results clearly indicate that climatic changes affect both the flux and composition of weathered material delivered to the oceans on glacial-interglacial timescales.

Selected age-diagnostic datums of ODP Site 184-1148

The 853 m thick sediment sequence recovered at ODP Site 1148 provides an unprecedented record of tectonic and paleoceanographic evolution in the South China Sea over the past 33 Ma. Litho-, bio-, and chemo-stratigraphic studies helped identify six periods of changes marking the major steps of the South China Sea geohistory. Rapid deposition with sedimentation rates of 60 m/Ma or more characterized the early Oligocene rifting. Several unconformities from the slumped unit between 457 and 495 mcd together erased about 3 Ma late Oligocene record, providing solid evidence of tectonic transition from rifting/slow spreading to rapid spreading in the South China Sea. Slow sedimentation of ~20-30 m/Ma signifies stable seafloor spreading in the early Miocene. Dissolution may have affected the completeness of Miocene-Pleistocene succession with short-term hiatuses beyond current biostratigraphical resolution. Five major dissolution events, D-1 to D-5, characterize the stepwise development of deep water masses in close association to post-Oligocene South China Sea basin transformation. The concurrence of local and global dissolution events in the Miocene and Pliocene suggests climatic forcing as the main mechanism causing deep water circulation changes concomitantly in world oceans and in marginal seas. A return of high sedimentation rate of 60 m/Ma to the late Pliocene and Pleistocene South China Sea was caused by intensified down-slope transport due to frequent sea level fluctuations and exposure of a large shelf area during sea level low-stands. The six paleoceanographic stages, respectively corresponding to rifting (~33-28.5 Ma), changing spreading southward (28.5-23 Ma), stable spreading to end of spreading (23-15 Ma), post-spreading balance (15-9 Ma), further modification and monsoon influence (9-5 Ma), and glacial prevalence (5-0 Ma), had transformed the South China Sea from a series of deep grabens to a rapidly expanding open gulf and finally to a semi-enclosed marginal sea in the past 33 Ma.

1.2 Ma stacked benthic foraminiferal carbon isotope record

Pleistocene stable carbon isotope (d13C) records from surface and deep dwelling foraminifera in all major ocean basins show two distinct long-term carbon isotope fluctuations since 1.00 Ma. The first started around 1.00 Ma and was characterised by a 0.35 per mil decrease in d13C values until 0.90 Ma, followed by an increase of 0.60 per mil lasting until 0.50 Ma. The subsequent fluctuation started with a 0.40 per mil decrease between 0.50 and 0.25 Ma, followed by an increase of 0.30 per mil between 0.25 and 0.10 Ma. Here, we evaluate existing evidence and various hypotheses for these global Pleistocene d13C fluctuations and present an interpretation, where the fluctuations most likely resulted from concomitant changes in the burial fluxes of organic and inorganic carbon due to ventilation changes and/or changes in the production and export ratio. Our model indicates that to satisfy the long-term 'stability' of the Pleistocene lysocline, the ratio between the amounts of change in the organic and inorganic carbon burial fluxes would have to be close to a 1:1 ratio, as deviations from this ratio would lead to sizable variations in the depth of the lysocline. It is then apparent that the mid-Pleistocene climate transition, which, apart from the glacial cycles, represents the most fundamental change in the Pleistocene climate, was likely not associated with a fundamental change in atmospheric pCO2. While recognising that high frequency glacial/interglacial cycles are associated with relatively large (100 ppmv) changes in pCO2, our model scenario (with burial changes close to a 1:1 ratio) produces a maximum long-term variability of only 20 ppmv over the fluctuation between 1.00 and 0.50 Ma.