Simulation Model of Benthic Antarctic Assemblages (SIMBAA), link to PC-executable

SIMBAA is a spatially explicit, individual-based simulation model. It was developed to analyse the response of populations of Antarctic benthic species and their diversity to iceberg scouring. This disturbance is causing a high local mortality providing potential space for new colonisation. Traits can be attributed to model species, e.g. in terms of reproduction, dispersal, and life span. Physical disturbances can be designed in space and time, e.g. in terms of size, shape, and frequency. Environmental heterogeneity can be considered by cell-specific capacities to host a certain number of individuals. When grid cells become empty (after a disturbance event or due to natural mortality of of an individual), a lottery decides which individual from which species stored in a pool of candidates (for this cell) will recruit in that cell. After a defined period the individuals become mature and their offspring are dispersed and stored in the pool of candidates. The biological parameters and disturbance regimes decide on how long an individual lives. Temporal development of single populations of species as well as Shannon diversity are depicted in the main window graphically and primary values are listed. Examples for simulations can be loaded and saved as sgf-files. The results are also shown in an additional window in a dimensionless area with 50 x 50 cells, which contain single individuals depicted as circles; their colour indicates the assignment to the self-designed model species and the size represents their age. Dominant species per cell and disturbed areas can also be depicted. Output of simulation runs can be saved as images, which can be assembled to video-clips by standard computer programs (see GIF-examples of which "Demo 1" represents the response of the Antarctic benthos to iceberg scouring and "Demo 2" represents a simulation of a deep-sea benthic habitat).

Water level time series and sediment grain size model

The environment of ebb-tidal deltas between barrier island systems is characterized by a complex morphology with ebb- and flood-dominated channels, shoals and swash bars connecting the ebb-tidal delta platform to the adjacent island. These morphological features reveal characteristic surface sediment grain-size distributions and are subject to a continuous adaptation to the prevailing hydrodynamic forces. The mixed-energy tidal inlet Otzumer Balje between the East Frisian barrier islands of Langeoog and Spiekeroog in the southern North Sea has been chosen here as a model study area for the identification of relevant hydrodynamic drivers of morphology and sedimentology. We compare the effect of high-energy, wave-dominated storm conditions to mid-term, tide-dominated fair-weather conditions on tidal inlet morphology and sedimentology with a process-based numerical model. A multi-fractional approach with five grain-size fractions between 150 and 450 µm allows for the simulation of corresponding surface sediment grain-size distributions. Net sediment fluxes for distinct conditions are identified: during storm conditions, bed load sediment transport is generally onshore directed on the shallower ebb-tidal delta shoals, whereas fine-grained suspended sediment bypasses the tidal inlet by wave-driven currents. During fair weather the sediment transport mainly focuses on the inlet throat and the marginal flood channels. We show how the observed sediment grain-size distribution and the morphological response at mixed-energy tidal inlets are the result of both wave-dominated less frequent storm conditions and mid-term, tide-dominant fair-weather conditions.

Decade Wind and Wave Data for NW Iberia, Cabo Sillerio

Seafloor sediment mobilization on the inner Northwest Iberian continental shelf is caused largely by ocean surface waves. The temporal and spatial variability in the wave height, wave period, and wave direction has a profound effect on local sediment mobilization, leading to distinct sediment mobilization scenarios. Six grain-size specific sediment mobilization scenarios, representing seasonal average and storm conditions, were simulated with a physics-based numerical model. Model inputs included meteorological and oceanographic data in conjunction with seafloor grain-size and the shelf bathymetric data. The results show distinct seasonal variations, most importantly in wave height, leading to sediment mobilization, specifically on the inner shelf shallower than 30 m water depth where up to 49% of the shelf area is mobilized. Medium to severe storm events are modeled to mobilize up to 89% of the shelf area above 150 m water depth. The frequency of each of these seasonal and storm-related sediment mobilization scenarios is addressed using a decade of meteorological and oceanographic data. The temporal and spatial patterns of the modeled sediment mobilization scenarios are discussed in the context of existing geological and environmental processes and conditions to assist scientific, industrial and environmental efforts that are directly affected by sediment mobilization. Examples, where sediment mobilization plays a vital role, include seafloor nutrient advection, recurrent arrival of oil from oil-spill-laden seafloor sediment, and bottom trawling impacts.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Supplementary material: XRF core scanning data, CaCO3 values, foraminiferal carbon and oxygen isotope data, Mg/Ca data, clay mineralogy and helium isotope results

The Paleocene-Eocene Thermal Maximum (PETM; ~56 Ma) is associated with abrupt climate change, carbon cycle perturbation, ocean acidification, as well as biogeographic shifts in marine and terrestrial biota that were largely reversed as the climatic transient waned. We report a clear exception to the behavior of the PETM as a reversing climatic transient in the eastern North Atlantic (Deep-Sea Drilling Project Site 401, Bay of Biscay) where the PETM initiates a greatly prolonged environmental change compared to other places on Earth where records exist. The observed environmental perturbation extended well past the d13C recovery phase and up to 650 kyr after the PETM onset according to our extraterrestrial 3He-based age-model. We observe a strong decoupling of planktic foraminiferal d18O and Mg/Ca values during the PETM d13C recovery phase, which in combination with results from helium isotopes and clay mineralogy, suggests that the PETM triggered a hydrologic change in western Europe that increased freshwater flux and the delivery of weathering products to the eastern North Atlantic. This state change persisted long after the carbon-cycle perturbation had stopped. We hypothesize that either long-lived continental drainage patterns were altered by enhanced hydrological cycling induced by the PETM, or alternatively that the climate system in the hinterland area of Site 401 was forced into a new climate state that was not easily reversed in the aftermath of the PETM.

Magnetic observation of ODP Site 152-919 in the Irminger Basin, North Atlantic

Component natural remanent magnetizations derived from u-channel and 1-qcm discrete samples from ODP Site 919 (Irminger Basin) indicate the existence of four intervals of negative inclinations in the upper Brunhes Chronozone. According to the age model based on planktic oxygen isotope data, these "excursional" intervals occur in sediments deposited during the following time intervals: 32-34 ka, 39-41 ka, 180-188 ka and 205-225 ka. These time intervals correspond to polarity excursions detected elsewhere, known as Mono Lake, Laschamp, Iceland Basin and Pringle Falls. The isotope-based age model is supported by the normalized remanence (paleointensity) record that can be correlated to other calibrated paleointensity records for the 0-500 ka interval, such as that from ODP Site 983. For the intervals associated with the Mono Lake and Laschamp excursions, virtual geomagnetic poles (VGPs) reach equatorial latitudes and mid-southerly latitudes, respectively. For intervals associated with the Iceland Basin and Pringle Falls excursions, repeated excursions of VGPs to high southerly latitudes indicate rapid directional swings rather than a single short-lived polarity reversal. The directional instability associated with polarity excursions is not often recorded, probably due to smoothing of the sedimentary record by the process of detrital remanence (DRM) acquisition.

Extraterrestrial 3He estimation across the Paleocene-Eocene thermal maximum at ODP Site 208-1266

In the deep-sea, the Paleocene-Eocene Thermal Maximum (PETM) is often marked by clay-rich condensed intervals caused by dissolution of carbonate sediments, capped by a carbonate-rich interval. Constraining the duration of both the dissolution and subsequent cap-carbonate intervals is essential to computing marine carbon fluxes and thus testing hypotheses for the origin of this event. To this end, we provide new high-resolution helium isotope records spanning the Paleocene-Eocene boundary at ODP Site 1266 in the South Atlantic. The extraterrestrial 3He, 3HeET, concentrations replicate trends observed at ODP Site 690 by Farley and Eltgroth (2003, doi:10.1016/S0012-821X(03)00017-7). By assuming a constant flux of 3HeET we constrain relative changes in accumulation rates of sediment across the PETM and construct a new age model for the event. In this new chronology the zero carbonate layer represents 35 kyr, some of which reflects clay produced by dissolution of Paleocene (pre-PETM) sediments. Above this layer, carbonate concentrations increase for ~165 kyr and remain higher than in the latest Paleocene until 234 +48/-34 kyr above the base of the clay. The new chronology indicates that minimum d13C values persisted for a maximum of 134 +27/-19 kyr and the inflection point previously chosen to designate the end of the CIE recovery occurs at 217 +44/-31 kyr. This allocation of time differs from that of the cycle-based age model of Röhl et al. (2007, doi:10.1029/2007GC001784) in that it assigns more time to the clay layer followed by a more gradual recovery of carbonate-rich sedimentation. The new model also suggests a longer sustained d13C excursion followed by a more rapid recovery to pre-PETM d13C values. These differences have important implications for constraining the source(s) of carbon and mechanisms for its subsequent sequestration, favoring models that include a sustained release

N2 fixation in eddies of the eastern tropical South Pacific Ocean

Mesoscale eddies play a major role in controlling ocean biogeochemistry. By impacting nutrient availability and water column ventilation, they are of critical importance for oceanic primary production. In the eastern tropical South Pacific Ocean off Peru, where a large and persistent oxygen-deficient zone is present, mesoscale processes have been reported to occur frequently. However, investigations into their biological activity are mostly based on model simulations, and direct measurements of carbon and dinitrogen (N2) fixation are scarce. We examined an open-ocean cyclonic eddy and two anticyclonic mode water eddies: a coastal one and an open-ocean one in the waters off Peru along a section at 16°S in austral summer 2012. Molecular data and bioassay incubations point towards a difference between the active diazotrophic communities present in the cyclonic eddy and the anticyclonic mode water eddies. In the cyclonic eddy, highest rates of N2 fixation were measured in surface waters but no N2 fixation signal was detected at intermediate water depths. In contrast, both anticyclonic mode water eddies showed pronounced maxima in N2 fixation below the euphotic zone as evidenced by rate measurements and geochemical data. N2 fixation and carbon (C) fixation were higher in the young coastal mode water eddy compared to the older offshore mode water eddy. A co-occurrence between N2 fixation and biogenic N2, an indicator for N loss, indicated a link between N loss and N2 fixation in the mode water eddies, which was not observed for the cyclonic eddy. The comparison of two consecutive surveys of the coastal mode water eddy in November 2012 and December 2012 also revealed a reduction in N2 and C fixation at intermediate depths along with a reduction in chlorophyll by half, mirroring an aging effect in this eddy. Our data indicate an important role for anticyclonic mode water eddies in stimulating N2 fixation and thus supplying N offshore.

Abundance, biomass and secondary production of benthic megafauna on the Barents Sea shelf in 2008 and 2009

Megabenthos plays a major role in the overall energy flow on Arctic shelves, but information on megabenthic secondary production on large spatial scales is scarce. Here, we estimated for the first time megabenthic secondary production for the entire Barents Sea shelf by applying a species-based empirical model to an extensive dataset from the joint Norwegian? Russian ecosystem survey. Spatial patterns and relationships were analyzed within a GIS. The environmental drivers behind the observed production pattern were identified by applying an ordinary least squares regression model. Geographically weighted regression (GWR) was used to examine the varying relationship of secondary production and the environment on a shelfwide scale. Significantly higher megabenthic secondary production was found in the northeastern, seasonally ice-covered regions of the Barents Sea than in the permanently ice-free southwest. The environmental parameters that significantly relate to the observed pattern are bottom temperature and salinity, sea ice cover, new primary production, trawling pressure, and bottom current speed. The GWR proved to be a versatile tool for analyzing the regionally varying relationships of benthic secondary production and its environmental drivers (R² = 0.73). The observed pattern indicates tight pelagic? benthic coupling in the realm of the productive marginal ice zone. Ongoing decrease of winter sea ice extent and the associated poleward movement of the seasonal ice edge point towards a distinct decline of benthic secondary production in the northeastern Barents Sea in the future.

Trace elements in benthic foraminifera of sediment core V30-49

Eastern tropical Atlantic benthic foraminiferal Ba/Ca and Cd/Ca data from core V30-949 (3093 m) reveal large inferred changes in nutrient concentrations of deep Atlantic waters during the last 250 kyr. Relative changes in North Atlantic Deep Water contribution to this site are estimated by scaling the V30-49 Cd/Ca record to values of modern end-member water masses; these estimates agree well with the relative structure and timing of circulation changes in the eastern tropical Atlantic reconstructed from a d13C record-based mixing model (Raymo et al., 1997, doi:10.1029/97PA01019). Temporal differences between V30-49 Cd/Ca and Ba/Ca records suggest that the Ba/Ca record reflects changes in circulation with an additional increase in the Ba composition of deep Atlantic water masses during glacial episodes, possibly resulting from increased productivity. Similarity between the d13C and Ba/Ca records suggests that carbon isotopes in the deep glacial Atlantic also reflect productivity increases.