Description: We report the results of a two-dimensional tomographic inversion of marine seismic refraction data from an array of ocean-bottom seismographs (OBSs), which produced an image of the crustal structure along the axial valley of the ultraslow spreading Mid-Cayman Spreading Center (MCSC). The seismic velocity model shows variations in the thickness and properties of the young oceanic crust that are consistent with the existence of two magmatic-tectonic segments along the 110 km long spreading center. Seismic wave speeds are consistent with exhumed mantle at the boundary between these two segments, but changes in the vertical gradient of seismic velocity suggest that volcanic crust occupies most of the axial valley seafloor along the seismic transect. The two spreading segments both have a low-velocity zone (LVZ) several kilometers beneath the seafloor, which may indicate the presence of shallow melt. However, the northern segment also has low seismic velocities (3 km/s) in a thick upper crustal layer (1.5–2.0 km), which we interpret as an extrusive volcanic section with high porosity and permeability. This segment hosts the Beebe vent field, the deepest known high-temperature black smoker hydrothermal vent system. In contrast, the southern spreading segment has seismic velocities as high as 4.0 km/s near the seafloor. We suggest that the porosity and permeability of the volcanic crust in the southern segment are much lower, thus limiting deep seawater penetration and hydrothermal recharge. This may explain why no hydrothermal vent system has been found in the southern half of the MCSC.

Description: Keypoints This contribution is a reply on a comment submitted by A. Argnani. The alternate interpretation of the wide-angle seismic model is discussed. The Alfeo Fault system is proposed to be the current location of STEP fault. Abstract Andrea Argnani in his comment on Dellong et al., 2020 (Geometry of the deep Calabrian subduction (Central Mediterranean Sea) from wide‐angle seismic data and 3D gravity modeling), proposes an alternate interpretation of the wide-angle seismic velocity models presented by Dellong et al., 2018 and Dellong et al., 2020 and proposes a correction of the literature citations in these paper. In this reply, we discuss in detail all points raised by Andrea Argnani.

Description: Raw-, SEG-Y and other supplementary data of the landside deployment from the amphibious wide-angle seismic experiment ALPHA are presented. The aim of this project was to reveal the crustal and lithospheric structure of the subducting Adriatic plate and the external accretionary wedge in the southern Dinarides. Airgun shots from the RV Meteor were recorded along two profiles across Montenegro and northern Albania.

Description: The Calabrian subduction zone is one of the narrowest arcs on Earth and a key area to understand the geodynamic evolution of the Mediterranean and other marginal seas. Here in the Ionian Sea, the African plate subducts beneath Eurasia. Imaging the boundary between the downgoing slab and the upper plate along the Calabrian subduction zone is important for assessing the potential of the subduction zone to generate mega‐thrust earthquakes and was the main objective of this study. Here we present and analyze the results from a 380 km long, wide‐angle seismic profile spanning the complete subduction zone, from the deep Ionian Basin and the accretionary wedge to NE Sicily, with additional constraints offered by 3‐D Gravity modeling and the analysis of earthquake hypocenters. The velocity model for the wide‐angle seismic profile images thin oceanic crust throughout the basin. The Calabrian backstop extends underneath the accretionary wedge to about 100 km SE of the coast. The seismic model was extended in depth using earthquake hypocenters. The combined results indicate that the slab dip increases abruptly from 2‐3° to 60‐70° over a distance of ≤50 km underneath the Calabrian backstop. This abrupt steepening is likely related to the roll‐back geodynamic evolution of the narrow Calabrian slab which shows great similarity to the shallow and deep geometry of the Gibraltar slab.

Description: Key points:  First insights into the crustal structure of the northeastern Lau Basin, along a 290 km transect at 17°20’S.  Crust in southern Fonualei Rift and Spreading Center was created by extension of arc crust and variable amount of magmatism.  Magmatic underplating is present in some parts of the southern Niuafo’ou Microplate The northeastern Lau Basin is one of the fastest opening and magmatically most active back‐arc regions on Earth. Although the current pattern of plate boundaries and motions in this complex mosaic of microplates is reasonably understood, the internal structure and evolution of the back‐arc crust are not. We present new geophysical data from a 290 km long east‐west oriented transect crossing the Niuafo’ou Microplate (back‐arc), the Fonualei Rift and Spreading Centre (FRSC) and the Tofua Volcanic Arc at 17°20’S. Our P‐wave tomography model and density modelling suggests that past crustal accretion inside the southern FRSC was accommodated by a combination of arc crustal extension and magmatic activity. The absence of magnetic reversals inside the FRSC supports this and suggests that focused seafloor spreading has until now not contributed to crustal accretion. The back‐arc crust constituting the southern Niuafo’ou Microplate reveals a heterogeneous structure comprising several crustal blocks. Some regions of the back‐arc show a crustal structure similar to typical oceanic crust, suggesting they originate from seafloor spreading. Other crustal blocks resemble a structure that is similar to volcanic arc crust or a ‘hydrous’ type of oceanic crust that has been created at a spreading center influenced by slab‐derived water at distances 〈 50 km to the arc. Throughout the back‐arc region we observe a high‐velocity (Vp 7.2‐7.5 km s‐1) lower crust, which is an indication for magmatic underplating, which is likely sustained by elevated upper mantle temperatures in this region.

Description: The dense AlpArray network allows studying seismic wave propagation with high spatial resolution. Here we introduce an array approach to measure arrival angles of teleseismic Rayleigh waves. The approach combines the advantages of phase correlation as in the two-station method with array beamforming to obtain the phase-velocity vector. 20 earthquakes from the first two years of the AlpArray project are selected, and spatial patterns of arrival-angle deviations across the AlpArray are shown in maps, depending on period and earthquake location. The cause of these intriguing spatial patterns is discussed. A simple wave-propagation modelling example using an isolated anomaly and a Gaussian beam solution suggests that much of the complexity can be explained as a result of wave interference after passing a structural anomaly along the wave paths. This indicates that arrival-angle information constitutes useful additional information on the Earth structure, beyond what is currently used in inversions.

Description: New marine geophysical data acquired across the partly ice‐covered northern East Greenland continental margin highlight a complex interaction between tectonic and magmatic events. Breakup‐related lava flows are imaged in reflection seismic data as seaward dipping reflectors (SDRs), which are found to decrease in size both northwards and southwards from a central point at 75° N. We provide evidence that the magnetic anomaly pattern in the shelf area is related to volcanic phases and not to the presence of oceanic crust. The remnant magnetization of the individual lava flows is used to deduce a relative timing of the emplacement of the volcanic wedges. We find that the SDRs have been emplaced over a period of 2‐4 Ma progressively from north to south and from landward to seaward. The new data indicate a major post‐middle Eocene magmatic phase around the landward termination of the West Jan Mayen Fracture Zone. This post‐40 Ma volcanism likely was associated with the progressive separation of the Jan Mayen microcontinent from East Greenland. The break‐up of the Greenland Sea started at several isolated seafloor spreading cells whose location was controlled by rift structures and led to the present‐day segmentation of the margin. The original rift basins were subsequently connected by steady‐state seafloor spreading that propagated southwards, from the Greenland Fracture Zone to the Jan Mayen Fracture Zone. Key Points Polyphase Cenozoic volcanic rifting and consecutive emplacement of breakup‐related lava flows units along the northern East Greenland margin Breakup along restricted margin segments is followed by north to south directed progressive opening of the Greenland Sea Widespread post‐middle Eocene (〈 40 Ma) offshore magmatism, associated with the breakup of the Jan Mayen microcontinent from East Greenland

Description: We present a new 3-D shear-velocity model for the top 30 km of the crust in the wider Vienna Basin region based on surface waves extracted from ambient-noise cross-correlations. We use continuous seismic records of 63 broad-band stations of the AlpArray project to retrieve interstation Green’s functions from ambient-noise cross-correlations in the period range from 5 to 25 s. From these Green’s functions, we measure Rayleigh group traveltimes, utilizing all four components of the cross-correlation tensor, which are associated with Rayleigh waves (ZZ, RR, RZ and ZR), to exploit multiple measurements per station pair. A set of selection criteria is applied to ensure that we use high-quality recordings of fundamental Rayleigh modes. We regionalize the interstation group velocities in a 5 km × 5 km grid with an average path density of ∼20 paths per cell. From the resulting group-velocity maps, we extract local 1-D dispersion curves for each cell and invert all cells independently to retrieve the crustal shear-velocity structure of the study area. The resulting model provides a previously unachieved lateral resolution of seismic velocities in the region of ∼15 km. As major features, we image the Vienna Basin and Little Hungarian Plain as low-velocity anomalies, and the Bohemian Massif with high velocities. The edges of these features are marked with prominent velocity contrasts correlated with faults, such as the Alpine Front and Vienna Basin transfer fault system. The observed structures correlate well with surface geology, gravitational anomalies and the few known crystalline basement depths from boreholes. For depths larger than those reached by boreholes, the new model allows new insight into the complex structure of the Vienna Basin and surrounding areas, including deep low-velocity zones, which we image with previously unachieved detail. This model may be used in the future to interpret the deeper structures and tectonic evolution of the wider Vienna Basin region, evaluate natural resources, model wave propagation and improve earthquake locations, among others.

Description: Crustal structure provides the key to understand the interplay of magmatism and tectonism while oceanic crust is constructed at Mid Ocean Ridges (MOR). At slow spreading rates, magmatic processes dominate central areas of MOR segments, whereas segment ends are highly tectonised. The TAMMAR segment at the Mid-Atlantic Ridge (MAR) between 21°25' N and 22° N is a magmatically active segment. At ~4.5 Ma this segment started to propagate south, causing the termination of the transform fault at 21°40' N. This stopped long-lived detachment faulting and caused the migration of the ridge offset to the south. Here, a segment centre with a high magmatic budget has replaced a transform fault region with limited magma supply. We present results from seismic refraction profiles that mapped the crustal structure across the ridge crest of the TAMMAR segment. Seismic data yield crustal structure changes at the segment centre as a function of melt supply. Seismic Layer 3 underwent profound changes in thickness and became rapidly thicker ~5 Ma. This correlates with the observed “Bull's eye” gravimetric anomaly in that region. Our observations support a temporal change from thick lithosphere with oceanic core complex formation and transform faulting to thin lithosphere with focused mantle upwelling and segment growth. Temporal changes in crustal construction are connected to variations in the underlying mantle. We propose there is a link between the neighbouring segments at a larger scale within the asthenosphere, to form a long, highly magmatically active macro segment, here called the TAMMAR-Kane MacroSegment.

Description: We present the results of P-to-S receiver function analysis to improve the 3D image of the sedimentary layer, the upper crust, and lower crust in the Pannonian Basin area. The Pannonian Basin hosts deep sedimentary depocentres superimposed on a complex basement structure and it is surrounded by mountain belts. We processed waveforms from 221 three-component broadband seismological stations. As a result of the dense station coverage, we were able to achieve so far unprecedented spatial resolution in determining the velocity structure of the crust. We applied a three-fold quality control process; the first two being applied to the observed waveforms and the third to the calculated radial receiver functions. This work is the first comprehensive receiver function study of the entire region. To prepare the inversions, we performed station-wise H-Vp/Vs grid search, as well as Common Conversion Point migration. Our main focus was then the S-wave velocity structure of the area, which we determined by the Neighborhood Algorithm inversion method at each station, where data were sub-divided into back-azimuthal bundles based on similar Ps delay times. The 1D, nonlinear inversions provided the depth of the discontinuities, shear-wave velocities and Vp/Vs ratios of each layer per bundle, and we calculated uncertainty values for each of these parameters. We then developed a 3D interpolation method based on natural neighbor interpolation to obtain the 3D crustal structure from the local inversion results. We present the sedimentary thickness map, the first Conrad depth map and an improved, detailed Moho map, as well as the first upper and lower crustal thickness maps obtained from receiver function analysis. The velocity jump across the Conrad discontinuity is estimated at less than 0.2 km/s over most of the investigated area. We also compare the new Moho map from our approach to simple grid search results and prior knowledge from other techniques. Our Moho depth map presents local variations in the investigated area: the crust-mantle boundary is at 20–26 km beneath the sedimentary basins, while it is situated deeper below the Apuseni Mountains, Transdanubian and North Hungarian Ranges (28–33 km), and it is the deepest beneath the Eastern Alps and the Southern Carpathians (40–45 km). These values reflect well the Neogene evolution of the region, such as crustal thinning of the Pannonian Basin and orogenic thickening in the neighboring mountain belts.