ALBERT

Description: Highlights • We review the knowledge on modern high-latitude planktic foraminifers. • Subpolar species currently invade higher latitudes. • Climate change affects phenology, seawater pH, and carbon turnover. • Modern planktic foraminifers are briefly discussed for their paleoceanographic significance. Abstract Planktic foraminifers can be sensitive indicators of the changing environment including both the Arctic Ocean and Southern Ocean. Due to variability in their ecology, biology, test characteristics, and fossil preservation in marine sediments, they serve as valuable archives in paleoceanography and climate geochemistry over the geologic time scale. Foraminifers are sensitive to, and can therefore provide proxy data on ambient water temperature, salinity, carbonate chemistry, and trophic conditions through shifts in assemblage (species) composition and the shell chemistry of individual specimens. Production and dissolution of the calcareous shell, as well as growth and remineralization of the cytoplasm, affect the carbonate counter pump and to a lesser extent the soft-tissue pump, at varying regional and temporal scales. Diversity of planktic foraminifers in polar waters is low in comparison to lower latitudes and is limited to three native species: Neogloboquadrina pachyderma, Turborotalita quinqueloba, and Globigerina bulloides, of which N. pachyderma is best adapted to polar conditions in the surface ocean. Neogloboquadrina pachyderma hibernates in brine channels in the lower layers of the Antarctic sea ice, a strategy that is presently undescribed in the Arctic. In open Antarctic and Arctic surface waters T. quinqueloba and G. bulloides increase in abundance at lower polar to subpolar latitudes and Globigerinita uvula, Turborotalita humilis, Globigerinita glutinata, Globorotalia inflata, and Globorotalia crassaformis complement the assemblages. Over the past two to three decades there has been a marked increase in the abundance of Orcadia riedeli and G. uvula in the subpolar and polar Indian Ocean, as well as in the northern North Atlantic. This paper presents a review of the knowledge of polar and subpolar planktic foraminifers. Particular emphasis is placed on the response of foraminifers to modern warming and ocean acidification at high latitudes and the implications for data interpretation in paleoceanography and paleoclimate research.

Description: Highlights • Deep-sea mineral exploration and exploitation licenses have been issued recently. • Mining will modify the abiotic and biotic environment. • At directly mined sites, species are removed and cannot resist disturbance. • Recovery is highly variable in distinct ecosystems and among benthic taxa. • Community changes may persist over geological time-scales at directly mined sites. Abstract With increasing demand for mineral resources, extraction of polymetallic sulphides at hydrothermal vents, cobalt-rich ferromanganese crusts at seamounts, and polymetallic nodules on abyssal plains may be imminent. Here, we shortly introduce ecosystem characteristics of mining areas, report on recent mining developments, and identify potential stress and disturbances created by mining. We analyze species’ potential resistance to future mining and perform meta-analyses on population density and diversity recovery after disturbances most similar to mining: volcanic eruptions at vents, fisheries on seamounts, and experiments that mimic nodule mining on abyssal plains. We report wide variation in recovery rates among taxa, size, and mobility of fauna. While densities and diversities of some taxa can recover to or even exceed pre-disturbance levels, community composition remains affected after decades. The loss of hard substrata or alteration of substrata composition may cause substantial community shifts that persist over geological timescales at mined sites.

Description: As a major source for atmospheric CO2, the Peruvian upwelling region exhibits strong variability in surface fCO2 on short spatial and temporal scales. Understanding the physical processes driving the strong variability is of fundamental importance for constraining the effect of marine emissions from upwelling regions on the global CO2 budget. In this study, a frontal decay on length scales of 𝒪(10 km) was observed off the Peruvian coast following a pronounced decrease in down-frontal (equatorward) wind speed with a time lag of 9 h. Simultaneously, the sea-to-air flux of CO2 on the inshore (cold) side of the front dropped from up to 80 to 10 mmol m−2 day−1, while the offshore (warm) side of the front was constantly outgassing at a rate of 10–20 mmol m−2 day−1. Based on repeated ship transects the decay of the front was observed to occur in two phases. The first phase was characterized by a development of coherent surface temperature anomalies which gained in amplitude over 6–9 h. The second phase was characterized by a disappearance of the surface temperature front within 6 h. Submesoscale mixed-layer instabilities were present but seem too slow to completely remove the temperature gradient in this short time period. Dynamics such as a pressure-driven gravity current appear to be a likely mechanism behind the evolution of the front.

Description: Highlights • New multiscale seismic data of the Carboneras Fault Zone (CFZ). • The tectonic architecture and depth geometry of the Carboneras Fault is examined. • We characterise fault segments and sub-segments to estimate their seismic potential. • The basement plays a key role in the actual configuration of the fault. • We explore CFZ terminations to know how strain is transferred to nearby structures. Abstract In the SE Iberian Margin, which hosts the convergent boundary between the European and African Plates, Quaternary faulting activity is dominated by a large left-lateral strike–slip system referred to as the Eastern Betic Shear Zone. This active fault system runs along more than 450 km and it is characterised by low to moderate magnitude shallow earthquakes, although large historical events have also occurred. The Carboneras Fault is the longest structure of the Eastern Betic Shear Zone, and its southern termination extends further into the Alboran Sea. Previously acquired high-resolution data (i.e. swath-bathymetry, TOBI sidescan sonar and sub-bottom profiler) show that the offshore Carboneras Fault is a NE–SW-trending upwarped zone of deformation with a length of 90 km long and a width of 0.5 to 2 km, which shows geomorphic features typically found in subaerial strike–slip faults, such as deflected drainage, pressure ridges and “en echelon” folds. However, the neotectonic, depth architecture, and Neogene evolution of Carboneras Fault offshore are still poorly known. In this work we present a multiscale seismic imaging of the Carboneras Fault (i.e. TOPAS, high-resolution multichannel-seismic reflection, and deep penetration multichannel-seismic reflection) carried out during three successive marine cruises, from 2006 to 2010. The new dataset allowed us to define a total of seven seismostratigraphic units (from Tortonian to Late Quaternary) above the basement, to characterise the tectonic architecture and structural segmentation of the Carboneras Fault, and to estimate its maximum seismic potential. We finally discuss the role of the basement in the present-day tectonic evolution of the Carboneras Fault, and explore the northern and southern terminations of the fault and how the strain is transferred to nearby structures.

Description: High-resolution acoustic and seismic data acquired 100 km offshore Cape São Vicente, image with unprecedented detail one of the largest active reverse faults of the SW Iberian Margin, the Horseshoe Fault (HF). The HF region is an area seismogenically active, source of the largest magnitude instrumental and historical earthquake (Mw〉6) occurred in the SW Iberian Margin. The HF corresponds to a N40 trending, 110 km long, and NW-verging active thrust that affects the whole sedimentary sequence and reaches up to the seafloor, generating a relief of more than 1 km. The along-strike structural variability as well as fault trend suggests that the HF is composed by three main sub-segments: North (N25), Central (N50) and South (N45). Swath-bathymetry, TOBI sidescan sonar backscatter and parametric echosounder TOPAS profiles reveal the surface morphology of the HF block, characterized by several, steep (20º) small scarps located on the hangingwall, and a succession of mass transport deposits (i.e. turbidites) on its footwall, located in the Horseshoe Abyssal Plain. A succession of pre-stack depth-migrated multichannel seismic reflection profiles across the HF and neighboring areas allowed us to constrain their seismo-stratigraphy, structural geometry, tectono-sedimentary evolution from Upper Jurassic to present-day, and to calculate their fault parameters. Finally, on the basis of segment length, surface fault area and seismogenic depth we evaluated the seismic potential of the HF, which in the worst-case scenario may generate an earthquake of magnitude Mw 7.8 ± 0.1. Thus, considering the tectonic behavior and near-shore location, the HF should be recognized in seismic and tsunami hazard assessment models of Western Europe and North Africa.

Description: An important tool for deep-sea temperature reconstruction is Mg/Ca paleothermometry applied to benthic foraminifera. Foraminifera of the genus Melonis appear to be promising candidates for temperature reconstructions due to their wide geographical and bathymetric distribution, and their infaunal habitat, which was suggested to reduce secondary effects from carbonate ion saturation (Δ[CO3 2−]). Here, we make substantial advances to previous calibration efforts and present new multi-lab Mg/Ca data for Melonis barleeanum and Melonis pompilioides from more than one hundred core top samples spanning in situ bottom temperatures from −1 to 16 °C, coupled with morphometric analyses of the foraminifer tests. Both species and their morphotypes seem to have a similar response of Mg/Ca to growth temperature. Compilation of new and previously published data reveals a linear dependence of temperature on Mg/Ca, with a best fit of Mg/Ca (mmol/mol) = 0.113 ± 0.005 ∗ BWT (°C) + 0.792 ± 0.036 (r2 = 0.81; n = 120; 1σ SD). Salinity, bottom water Δ[CO3 2−], and varying morphotypes have no apparent effect on the Mg/Ca-temperature relationship, but pore water Δ[CO3 2−] might have had an influence on some of the samples from the tropical Atlantic.

Description: Complex multifault earthquake ruptures involving secondary faults emphasize the necessity to characterize their seismogenic potential better and study their relationship with major faults to improve the seismic hazard assessment of a region. High-resolution geophysical data were interpreted to make a detailed characterization of the Averroes Fault and the North Averroes Faults, which are poorly known secondary right-lateral strike-slip faults located in the central part of the Alboran Sea (western Mediterranean). These faults appear to have evolved since the Pliocene as part of a distributed dextral strike-slip shear zone in response to local strain engendered by the diverging movement of the Carboneras Fault to the north, and the Yusuf and Alboran Ridge faults to the south. In addition, the architecture of these faults suggests that the Averroes Fault may eventually link with the Yusuf fault, thus leading to a higher seismogenic potential. Therefore, these secondary faults represent a hitherto unrecognized seismogenic hazard since they could produce earthquakes up to moment magnitude (Mw) 7.6. Our results highlight the importance of the role played by secondary faults in a specific kinematic framework. Their reciprocal linkage and their mechanical relationship with the main faults could lead to future complex fault ruptures. This information could improve fault source and earthquake models used in seismic and tsunami hazard assessment in this and similar regions.

Description: The ocean is a substantial source of nitrous oxide (N2O) to the atmosphere, but little is known about how this flux might change in the future. Here, we investigate the potential evolution of marine N2O emissions in the 21st century in response to anthropogenic climate change using the global ocean biogeochemical model NEMO-PISCES. Assuming nitrification as the dominant N2O formation pathway, we implemented two different parameterizations of N2O production which differ primarily under low-oxygen (O2) conditions. When forced with output from a climate model simulation run under the business-as-usual high-CO2 concentration scenario (RCP8.5), our simulations suggest a decrease of 4 to 12 % in N2O emissions from 2005 to 2100, i.e., a reduction from 4.03/3.71 to 3.54/3.56 TgN yr−1 depending on the parameterization. The emissions decrease strongly in the western basins of the Pacific and Atlantic oceans, while they tend to increase above the oxygen minimum zones (OMZs), i.e., in the eastern tropical Pacific and in the northern Indian Ocean. The reduction in N2O emissions is caused on the one hand by weakened nitrification as a consequence of reduced primary and export production, and on the other hand by stronger vertical stratification, which reduces the transport of N2O from the ocean interior to the ocean surface. The higher emissions over the OMZ are linked to an expansion of these zones under global warming, which leads to increased N2O production, associated primarily with denitrification. While there are many uncertainties in the relative contribution and changes in the N2O production pathways, the increasing storage seems unequivocal and determines largely the decrease in N2O emissions in the future. From the perspective of a global climate system, the averaged feedback strength associated with the projected decrease in oceanic N2O emissions amounts to around −0.009 W m−2 K−1, which is comparable to the potential increase from terrestrial N2O sources. However, the assessment for a potential balance between the terrestrial and marine feedbacks calls for an improved representation of N2O production terms in fully coupled next-generation Earth system models.

Description: IODP Expedition 340 successfully drilled a series of sites offshore Montserrat, Martinique and Dominica in the Lesser Antilles from March to April 2012. These are among the few drill sites gathered around volcanic islands, and the first scientific drilling of large and likely tsunamigenic volcanic island-arc landslide deposits. These cores provide evidence and tests of previous hypotheses for the composition and origin of those deposits. Sites U1394, U1399, and U1400 that penetrated landslide deposits recovered exclusively seafloor-sediment, comprising mainly turbidites and hemipelagic deposits, and lacked debris avalanche deposits. This supports the concepts that i/ volcanic debris avalanches tend to stop at the slope break, and ii/ widespread and voluminous failures of pre-existing low-gradient seafloor sediment can be triggered by initial emplacement of material from the volcano. Offshore Martinique (U1399 and 1400), the landslide deposits comprised blocks of parallel strata that were tilted or micro-faulted, sometimes separated by intervals of homogenized sediment (intense shearing), while Site U1394 offshore Montserrat penetrated a flat-lying block of intact strata. The most likely mechanism for generating these large-scale seafloor-sediment failures appears to be propagation of a decollement from proximal areas loaded and incised by a volcanic debris avalanche. These results have implications for the magnitude of tsunami generation. Under some conditions, volcanic island landslide deposits comprised of mainly seafloor sediment will tend to form smaller magnitude tsunamis than equivalent volumes of subaerial block-rich mass flows rapidly entering water. Expedition 340 also successfully drilled sites to access the undisturbed record of eruption fallout layers intercalated with marine sediment which provide an outstanding high-resolution dataset to analyze eruption and landslides cycles, improve understanding of magmatic evolution as well as offshore sedimentation processes. This article is protected by copyright. All rights reserved.

Description: The surface of the solid Earth is effectively stress free in its subaerial portions, and hydrostatic beneath the oceans. Unfortunately, this type of boundary condition is difficult to treat computationally, and for computational convenience, numerical models have often used simpler approximations that do not involve a normal stress-loaded, shear-stress free top surface that is free to move. Viscous flow models with a computational free surface typically confront stability problems when the time step is bigger than the viscous relaxation time. The small time step required for stability (〈2. Kyr) makes this type of model computationally intensive, so there remains a need to develop strategies that mitigate the stability problem by making larger (at least ~10 Kyr) time steps stable and accurate. Here we present a new free-surface stabilization algorithm for finite element codes which solves the stability problem by adding to the Stokes formulation an intrinsic penalization term equivalent to a portion of the future load at the surface nodes. Our algorithm is straightforward to implement and can be used with both Eulerian or Lagrangian grids. It includes α and β parameters to respectively control both the vertical and the horizontal slope-dependent penalization terms, and uses Uzawa-like iterations to solve the resulting system at a cost comparable to a non-stress free surface formulation. Four tests were carried out in order to study the accuracy and the stability of the algorithm: (1) a decaying first-order sinusoidal topography test, (2) a decaying high-order sinusoidal topography test, (3) a Rayleigh-Taylor instability test, and (4) a steep-slope test. For these tests, we investigate which α and β parameters give the best results in terms of both accuracy and stability. We also compare the accuracy and the stability of our algorithm with a similar implicit approach recently developed by Kaus et al. (2010). We find that our algorithm is slightly more accurate and stable for steep slopes, and also conclude that, for longer time steps, the optimal α controlling factor for both approaches is ~2/3, instead of the 1/2 Crank-Nicolson parameter inferred from a linearized accuracy analysis. This more-implicit value coincides with the velocity factor for a Galerkin time discretization applied to our penalization term using linear shape functions in time.