ALBERT

Description: Continental hyperextension during magma-poor rifting at the Deep Galicia Margin is characterised by a complex pattern of faulting, thin continental fault blocks, and the serpentinisation, with local exhumation, of mantle peridotites along the S-reflector, interpreted as a detachment surface. In order to understand fully the evolution of these features, it is important to image seismically the structure and to model the velocity structure to the greatest resolution possible. Travel-time tomography models have revealed the long-wavelength velocity structure of this hyperextended domain, but are often insufficient to match accurately the short-wavelength structure observed in reflection seismic imaging. Here we demonstrate the application of two-dimensional (2D) time-domain acoustic full-waveform inversion to deep water seismic data collected at the Deep Galicia Margin, in order to attain a high resolution velocity model of continental hyperextension. We have used several quality assurance procedures to assess the velocity model, including comparison of the observed and modelled waveforms, checkerboard tests, testing of parameter and inversion strategy, and comparison with the migrated reflection image. Our final model exhibits an increase in the resolution of subsurface velocities, with particular improvement observed in the westernmost continental fault blocks, with a clear rotation of the velocity field to match steeply dipping reflectors. Across the S-reflector there is a sharpening in the velocity contrast, with lower velocities beneath S indicative of preferential mantle serpentinisation. This study supports the hypothesis that normal faulting acts to hydrate the upper mantle peridotite, observed as a systematic decrease in seismic velocities, consistent with increased serpentinisation. Our results confirm the feasibility of applying the full-waveform inversion method to sparse, deep water crustal datasets.

Description: Ocean acidification is an escalating environmental issue and associated changes in the ocean carbonate system have implications for many calcifying organisms. The present study followed the growth of Sepia officinalis from early-stage embryos, through hatching, to 7-week-old juveniles. Responses of cuttlefish to elevated pCO(2) (hypercapnia) were investigated to test the impacts of near-future and extreme ocean acidification conditions on growth, developmental time, oxygen consumption, and yolk utilization as proxies for individual fitness. We further examined gross morphological characteristics of the internal calcareous cuttlebone to determine whether embryonically secreted shell lamellae are impacted by environmental hypercapnia. Embryonic growth was reduced and hatching delayed under elevated pCO(2), both at environmentally relevant levels (0.14 kPa pCO(2) similar to predicted ocean conditions in 2100) and extreme conditions (0.40 kPa pCO(2)). Comparing various metrics from control and intermediate treatments generally showed no significant difference in experimental measurements. Yet, results from the high pCO(2) treatment showed significant changes compared with controls and revealed a consistent general trend across the three treatment levels. The proportion of animal mass contributed by the cuttlebone increased in both elevated pCO(2) treatments. Gross cuttlebone morphology was affected under such conditions and cuttlebones of hypercapnic individuals were proportionally shorter. Embryonic shell morphology was maintained consistently in all treatments, despite compounding hypercapnia in the perivitelline fluid; however, post-hatching, hypercapnic animals developed denser cuttlebone laminae in shorter cuttlebones. Juvenile cuttlefish in acidified environments thus experience lower growth and yet increased calcification of their internal shell. The results of this study support recent findings that early cuttlefish life stages are more vulnerable towards hypercapnia than juveniles and adults, which may have negative repercussions on the biological fitness of cuttlefish hatchlings in future oceans.

Description: The Mozambique Ridge, a prominent basement high in the southwestern Indian Ocean, consists of four major geomorphological segments associated with numerous phases of volcanic activity in the Lower Cretaceous. The nature and origin of the Mozambique Ridge have been intensely debated with one hypothesis suggesting a Large Igneous Province origin. High-resolution seismic reflection data reveal a large number of extrusion centres with a random distribution throughout the southern Mozambique Ridge and the nearby Transkei Rise. Intrabasement reflections emerge from the extrusion centres and are interpreted to represent massive lava flow sequences. Such lava flow sequences are characteristic of eruptions leading to the formation of continental and oceanic flood basalt provinces, hence supporting a Large Igneous Province origin of the Mozambique Ridge. We observe evidence for widespread post-sedimentary magmatic activity that we correlate with a southward propagation of the East African Rift System. Based on our volumetric analysis of the southern Mozambique Ridge we infer a rapid sequential emplacement between ∼131 and ∼125 Ma, which is similar to the short formation periods of other Large Igneous Provinces like the Agulhas Plateau.

Description: New marine geophysical data recorded across the Tonga-Kermadec subduction zone are used to image deformation and seismic velocity structures of the forearc and Pacific Plate where the Louisville Ridge seamount chain subducts. Due to the obliquity of the Louisville Ridge to the trench and the fast 128 mm yr−1 south–southwest migration of the ridge-trench collision zone, post-, current and pre-seamount subduction deformation can be investigated between 23°S and 28°S. We combine our interpretations from the collision zone with previous results from the post- and pre-collision zones to define the along-arc variation in deformation due to seamount subduction. In the pre-collision zone the lower-trench slope is steep, the mid-trench slope has ∼3-km-thick stratified sediments and gravitational collapse of the trench slope is associated with basal erosion by subducting horst and graben structures on the Pacific Plate. This collapse indicates that tectonic erosion is a normal process affecting this generally sediment starved subduction system. In the collision zone the trench-slope decreases compared to the north and south, and rotation of the forearc is manifest as a steep plate boundary fault and arcward dipping sediment in a 12-km-wide, ∼2-km-deep mid-slope basin. A ∼3 km step increase in depth of the middle and lower crustal isovelocity contours below the basin indicates the extent of crustal deformation on the trench slope. At the leading edge of the overriding plate, upper crustal P-wave velocities are ∼4.0 km s−1 and indicate the trench fill material is of seamount origin. Osbourn Seamount on the outer rise has extensional faulting on its western slope and mass wasting of the seamount provides the low Vp material to the trench. In the post-collision zone to the north, the trench slope is smooth, the trench is deep, and the crystalline crust thins at the leading edge of the overriding plate where Vp is low, ∼5.5 km s−1. These characteristics are attributed to a greater degree of extensional collapse of the forearc in the wake of seamount subduction. The northern end of a seismic gap lies at the transition from the smooth lower-trench slope of the post-collision zone, to the block faulted and elevated lower-trench slope in the collision zone, suggesting a causative link between the collapse of the forearc and seismogenesis. Along the forearc, the transient effects of a north-to-south progression of ridge subduction are preserved in the geomorphology, whereas longer-term effects may be recorded in the ∼80 km offset in trench strike at the collision zone itself.

Description: Gelatinous zooplankton (GZ) comprise a taxonomically and functionally diverse group of marine organisms which includes ctenophores, cnidarians and pelagic tunicates, sharing a soft, mostly transparent body texture, a high body water content and a lack of exoskeleton. They range in size from less than a millimetre to nearly 2 m for the cnidarian jellyfish Nemopilema nomurai, and comprise some of the fastest growing metazoans on Earth (Hopcroft et al., 1998), sometimes surpassing crustacean zooplankton in their contribution to secondary production (i.e. in subtropical waters; Jaspers et al., 2009). They feed on a wide range of prey sizes, with predator–prey ratios comparable in some cases to those of baleen whales and krill (Deibel and Lee, 1992), and with prey removal rates which are similar to those of their non-gelatinous competitors (Acuña et al., 2011). In spite of early work pointing to gelatinous zooplankton as a trophic dead end (Verity and Smetacek, 1996), evidence is rapidly accumulating which shows that they may potentially channel energy from the picoplankton-sized, microbial loop organisms up to the higher trophic levels, including fish (Llopiz et al., 2010). However, this pathway is still largely neglected in most food web investigations even though it is now becoming clear that GZ represent a major fraction of the diet of several commercially important fish species such as bluefin tuna (Thunnus thynnus) (Cardona et al., 2012).

Description: This paper applies nonlinear Bayesian inversion to marine controlled source electromagnetic (CSEM) data collected near two sites of the Integrated Ocean Drilling Program (IODP) Expedition 311 on the northern Cascadia Margin to investigate subseafloor resistivity structure related to gas hydrate deposits and cold vents. The Cascadia margin, off the west coast of Vancouver Island, Canada, has a large accretionary prism where sediments are under pressure due to convergent plate boundary tectonics. Gas hydrate deposits and cold vent structures have previously been investigated by various geophysical methods and seabed drilling. Here, we invert time-domain CSEM data collected at Sites U1328 and U1329 of IODP Expedition 311 using Bayesian methods to derive subsurface resistivity model parameters and uncertainties. The Bayesian information criterion is applied to determine the amount of structure (number of layers in a depth-dependent model) that can be resolved by the data. The parameter space is sampled with the Metropolis–Hastings algorithm in principal-component space, utilizing parallel tempering to ensure wider and efficient sampling and convergence. Nonlinear inversion allows analysis of uncertain acquisition parameters such as time delays between receiver and transmitter clocks as well as input electrical current amplitude. Marginalizing over these instrument parameters in the inversion accounts for their contribution to the geophysical model uncertainties. One-dimensional inversion of time-domain CSEM data collected at measurement sites along a survey line allows interpretation of the subsurface resistivity structure. The data sets can be generally explained by models with 1 to 3 layers. Inversion results at U1329, at the landward edge of the gas hydrate stability zone, indicate a sediment unconformity as well as potential cold vents which were previously unknown. The resistivities generally increase upslope due to sediment erosion along the slope. Inversion results at U1328 on the middle slope suggest several vent systems close to Bullseye vent in agreement with ongoing interdisciplinary observations.

Description: The ocean load in glacial isostatic adjustment (GIA) modelling is represented by the so-called sea level equation (SLE). The SLE describes the mass redistribution of water between ice sheets and oceans on a deforming Earth. Despite various teams independently investigating GIA, there has been no systematic intercomparison among the numerical solvers of the SLE through which the methods may be validated. The goal of this paper is to present a series of synthetic examples designed for testing and comparing the numerical implementations of the SLE in GIA modelling. The 10 numerical codes tested combine various temporal and spatial parametrizations. The time-domain or Laplace-domain discretizations are used to solve the SLE through time, while spherical harmonics, finite differences or finite elements parametrize the GIA-related field variables spatially. The surface ice-water load and solid Earth’s topography are represented spatially either on an equiangular grid, a Gauss–Legendre or an equiarea grid with icosahedron-shaped spherical pixels. Comparisons are made in a series of five benchmark examples with an increasing degree of complexity. Due to the complexity of the SLE, there is no analytical solution to it. The accuracy of the numerical implementations is therefore assessed by the differences of the individual solutions with respect to a reference solution. While the benchmark study does not result in GIA predictions for a realistic loading scenario, we establish a set of agreed-upon results that can be extended in the future by including more complex case studies, such as solutions with realistic loading scenarios, the rotational feedback in the linear-momentum equation, and by considering a 3-D viscosity structure of the Earth’s mantle. The test computations performed so far show very good agreement between the individual results and their ability to capture the main features of sea-surface variation and the surface vertical displacement. The differences found can often be attributed to the different approximations inherent in the various algorithms. This shows the accuracy that can be expected from different implementations of the SLE, which helps to assess differences noted in the literature between predictions for realistic loading cases.