Search

Plate tectonics describes oceanic transform faults as conservative strike-slip boundaries, where lithosphere is neither created nor destroyed. Therefore, seafloor accreted at ridge-transform intersections should follow a similar subsidence trend with age as lithosphere that forms away from ridge-transform intersections. Yet, recent compilations of high-resolution bathymetry show that the seafloor is significantly deeper along transform faults than at the adjacent fracture zones. We present residual mantle Bouguer anomalies, a proxy for crustal thickness, for 11 transform fault systems across the full range of spreading rates. Our results indicate that the crust is thinner in the transform deformation zone than in either the adjacent fracture zones or the inside corner regions. Consequently, oceanic transform faulting appears not only to thin the transform valley crust but also leads to a secondary phase of magmatic addition at the transition to the passive fracture zones. These observations challenge the concept of transform faults being conservative plate boundaries.

Among the recently evolving transform faults, the Blancotransform fault system (BTFS) in the northwest off the coast of Oregon,is strongly segmented and resembles other segmented oceanic transform systems, such as the Siqueiros in the East Pacific Rise, which developed from a pre-existing transform fault that was subjected to several extensional events due to a known shift in spreading direction. However,plate tectonic reconstructions suggested that the BTFS developed from at least two large ridge offsets rather than a single transform fault, emerging from a series of ridge propagation events after the plate reorientation at ~5 Ma.We investigated the seismotectonic behavior and tectonic evolution of the BTFS using ocean-bottom seismometer data from the Blanco Transform OBS Experiment (2012–2013) as well as high-resolution multibeam bathymetry, aeromagnetic, and gravity datasets. Interestingly, we observed no evidence for the existence of either transform faults or fracture zones around the BTFS before 2 Ma, indicating that there were no pre-existing transform faults before the formation of the BTFS. Therefore, we propose that the BTFS developed from two broad transfer zones rather than pre-existing transform faults. Moreover, local seismicity shows significant along-strike variations, revealing different modes of slip in the eastern and western BTFS.Seismic slip vectors suggest that the western BTFS is an immature transform fault system since its reorganization is still taking place, but the eastern BTFS accommodating the plate motion ismature. The BTFS acts as a natural laboratory to understand the mechanisms controlling the development of oceanic transform plate boundaries., The 28th IUGG General Assembly (IUGG2023) (Berlin 2023)

The Ligurian Sea is located in a complex tectonic area, at the Western Alps – Northern Apennines junction. Little is known about the structure of the basin and the continuation of the Alpine orogen offshore. The Ligurian Sea formed as a back-arc basin related to the retreat of the Apennine subduction zone. During extension, the Corsica-Sardinia block drifted eastwards, rotating counter-clockwise to its current position. The nature of the crust in the basin, whether atypical oceanic crust, thinned continental crust, or even exhumed mantle with thick sediment cover, is still a topic of debate. If oceanic crust was formed, its extent is poorly mapped. Moreover, seismic and geodetic data show that the basin is now under compression. The area is known to have hosted large, tsunamigenic earthquakes in the past, however, its location boarding France, Monaco, and Italy, and climate mean that the coastline is densely populated despite the risk. Understanding the nature of the tectonics in the basin can therefore help assess the seismic hazard of the area. In this project we use P-wave travel time tomography to generate a 3D velocity model of a section of the Gulf of Lion into the Ligurian Sea. This project is part of the DFG Priority Program “Mountain Building Processes in Four Dimensions (4DMB)”. We use data from the LOBSTER-AlpArray OBS as well as from the 2006 SARDINIA Experiment. With FMTOMO we incorporate reflected and refracted rays as well as water multiples which in some cases improves the ray coverage., The 28th IUGG General Assembly (IUGG2023) (Berlin 2023)

The Liguro-Provençal basin was formed as a back-arc basin of the retreating Calabrian-Apennines subduction zone during the Oligocene and Miocene. The resulting rotation of the Corsica-Sardinia block is associated with rifting, shaping the Ligurian Sea. It is still debated whether oceanic or atypical oceanic crust was formed or if the crust is continental and experienced extreme thinning during the opening of the basin. We invert velocity models using an amphibious network of seismic stations, including 22 broadband Ocean Bottom Seismometers (OBS) to investigate the lithospheric structure of the Ligurian sea. The instruments were installed in the Ligurian Sea for eight months between June 2017 and February 2018 as part of the AlpArray seismic network. Because of additional noise sources in the ocean, OBS data are rarely used for ambient noise studies. However, we attentively pre-process the data, including corrections for instrument tilt and seafloor compliance. We took extra care to exclude higher modes of the ambient-noise Rayleigh waves. We calculate daily cross-correlation functions for the LOBSTER array and surrounding land stations. Additionally, we correlate short time windows that include teleseismic earthquakes that allow us to derive surface wave group velocities for longer periods than using ambient noise only. Group velocity maps are obtained by inverting Green’s functions derived from the cross-correlation of ambient noise and teleseismic events, respectively. We then used the resulting 3D group velocity information to calculate 1D depth inversions for S-wave velocities. The shear-wave velocity results show a deepening of the Moho from 12 km at the southwestern basin centre to 20–25 km at the Ligurian coast in the northeast and over 30 km at the Provençal coast. We find no hint on mantle serpentinisation and no evidence for an Alpine slab, at least down to depths of 25 km. However, we see a separation of the southwestern and northeastern Ligurian Basin that coincides with the promoted prolongation of the Alpine front.

The northern margin of the Ligurian Basin shows notable seismicity at the Alpine front, including frequent magnitude 4 events. Seismicity decreases offshore towards the Basin centre and Corsica, revealing a diffuse distribution of low-magnitude earthquakes. We analyse data of the amphibious AlpArray seismic network with focus on the offshore component, the AlpArray ocean bottom seismometer (OBS) network, consisting of 24 broadband OBSs deployed for 8 months, to reveal the seismicity and depth distribution of micro-earthquakes beneath the Ligurian Sea. Two clusters occurred between ∼ 10 km to ∼ 16 km depth below the sea surface, within the lower crust and uppermost mantle. Thrust faulting focal mechanisms indicate compression and an inversion of the Ligurian Basin, which is an abandoned Oligocene–Miocene rift basin. The basin inversion is suggested to be related to the Africa–Europe plate convergence. The locations and focal mechanisms of seismicity suggest reactivation of pre-existing rift-related structures. Slightly different striking directions of presumed rift-related faults in the basin centre compared to faults further east and hence away from the rift basin may reflect the counter-clockwise rotation of the Corsica–Sardinia block. High mantle S-wave velocities and a low Vp/Vs ratio support the hypothesis of strengthening of crust and uppermost mantle during the Oligocene–Miocene rifting-related extension and thinning of continental crust.

The northern margin of the Ligurian Basin shows notable seismicity at the Alpine front, including frequent magnitude 4 events. Seismicity decreases offshore towards the Basin centre and Corsica, revealing a diffuse distribution of low magnitude earthquakes. We analyse data of the amphibious AlpArray seismic network with focus on the offshore component, the AlpArray OBS network, consisting of 24 broadband ocean bottom seismometers deployed for eight months, to reveal the seismicity and depth distribution of micro-earthquakes beneath the Ligurian Sea. Two clusters occurred between ~10 km to ~16 km depth below sea surface, within the lower crust and uppermost mantle. Thrust faulting focal mechanisms indicate compression and an inversion of the Ligurian Basin, which is an abandoned Oligocene rift basin. The Basin inversion is suggested to be related to the Africa-Europe plate convergence. The locations and focal mechanisms of seismicity suggest reactivation of pre-existing rift structures. Slightly different striking directions of faults in the basin centre compared to faults further east and hence away from the abandoned rift may mimic the counter-clockwise rotation of the Corsica-Sardinia block during ~20–16 Ma. The observed cluster events support the hypothesis of strengthening of crust and uppermost mantle during rifting related extension and thinning of continental crust.

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. ISSN:2169-9313 ISSN:0148-0227 ISSN:2169-9356

The region around the town Albstadt, SW Germany, was struck by four damaging earthquakes with magnitudes greater than 5 during the last century. These earthquakes occurred along the Albstadt Shear Zone (ASZ), which is characterized by more or less continuous microseismicity. As there are no visible surface ruptures that may be connected to the fault zone, we study its characteristics by its seismicity distribution and faulting pattern. We use the earthquake data of the state earthquake service of Baden-Württemberg from 2011 to 2018 and complement it with additional phase picks beginning in 2016 at the AlpArray and StressTransfer seismic networks in the vicinity of the ASZ. This extended data set is used to determine new minimum 1-D seismic vp and vs velocity models and corresponding station delay times for earthquake relocation. Fault plane solutions are determined for selected events, and the principal stress directions are derived. The minimum 1-D seismic velocity models have a simple and stable layering with increasing velocity with depth in the upper crust. The corresponding station delay times can be explained well by the lateral depth variation of the crystalline basement. The relocated events align about north–south with most of the seismic activity between the towns of Tübingen and Albstadt, east of the 9∘ E meridian. The events can be separated into several subclusters that indicate a segmentation of the ASZ. The majority of the 25 determined fault plane solutions feature an NNE–SSW strike but NNW–SSE-striking fault planes are also observed. The main fault plane associated with the ASZ dips steeply, and the rake indicates mainly sinistral strike-slip, but we also find minor components of normal and reverse faulting. The determined direction of the maximum horizontal stress of 140–149∘ is in good agreement with prior studies. Down to ca. 7–8 km depth SHmax is bigger than SV; below this depth, SV is the main stress component. The direction of SHmax indicates that the stress field in the area of the ASZ is mainly generated by the regional plate driving forces and the Alpine topography.

The recent development of ambient noise tomography, in combination with the increasing number of permanent seismic stations and dense networks of temporary stations operated during passive seismic experiments, provides a unique opportunity to build the first high-resolution 3-D shear wave velocity (vS) model of the entire crust of the Bohemian Massif (BM). This paper provides a regional-scale model of velocity distribution in the BM crust. The velocity model with a cell size of 22 km is built using a conventional two-step inversion approach from Rayleigh wave group velocity dispersion curves measured at more than 400 stations. The shear velocities within the upper crust of the BM are ∼0.2 km s−1 higher than those in its surroundings. The highest crustal velocities appear in its southern part, the Moldanubian unit. The Cadomian part of the region has a thinner crust, whereas the crust assembled, or tectonically transformed in the Variscan period, is thicker. The sharp Moho discontinuity preserves traces of its dynamic development expressed in remnants of Variscan subductions imprinted in bands of crustal thickening. A significant feature of the presented model is the velocity-drop interface (VDI) modelled in the lower part of the crust. We explain this feature by the anisotropic fabric of the lower crust, which is characterised as vertical transverse isotropy with the low velocity being the symmetry axis. The VDI is often interrupted around the boundaries of the crustal units, usually above locally increased velocities in the lowermost crust. Due to the north-west–south-east shortening of the crust and the late-Variscan strike-slip movements along the north-east–south-west oriented sutures preserved in the BM lithosphere, the anisotropic fabric of the lower crust was partly or fully erased along the boundaries of original microplates. These weakened zones accompanied by a velocity increase above the Moho (which indicate an emplacement of mantle rocks into the lower crust) can represent channels through which portions of subducted and later molten rocks have percolated upwards providing magma to subsequently form granitoid plutons.

We compiled a dataset of continuous recordings from the temporary and permanent seismic networks to compute the high-resolution 3D S-wave velocity model of the Southeastern Alps, the western part of the external Dinarides, and the Friuli and Venetian plains through ambient noise tomography. Part of the dataset is recorded by the SWATH-D temporary network and permanent networks in Italy, Austria, Slovenia and Croatia between October 2017 and July 2018. We computed 4050 vertical component cross-correlations to obtain the empirical Rayleigh wave Green’s functions. The dataset is complemented by adopting 1804 high-quality correlograms from other studies. The fast-marching method for 2D surface wave tomography is applied to the phase velocity dispersion curves in the 2–30 s period band. The resulting local dispersion curves are inverted for 1D S-wave velocity profiles using the non-perturbational and perturbational inversion methods. We assembled the 1D S-wave velocity profiles into a pseudo-3D S-wave velocity model from the surface down to 60 km depth. A range of iso-velocities, representing the crystalline basement depth and the crustal thickness, are determined. We found the average depth over the 2.8–3.0 and 4.1–4.3 km/s iso-velocity ranges to be reasonable representations of the crystalline basement and Moho depths, respectively. The basement depth map shows that the shallower crystalline basement beneath the Schio-Vicenza fault highlights the boundary between the deeper Venetian and Friuli plains to the east and the Po-plain to the west. The estimated Moho depth map displays a thickened crust along the boundary between the Friuli plain and the external Dinarides. It also reveals a N-S narrow corridor of crustal thinning to the east of the junction of Giudicarie and Periadriatic lines, which was not reported by other seismic imaging studies. This corridor of shallower Moho is located beneath the surface outcrop of the Permian magmatic rocks and seems to be connected to the continuation of the Permian magmatism to the deep-seated crust. We compared the shallow crustal velocities and the hypocentral location of the earthquakes in the Southern foothills of the Alps. It revealed that the seismicity mainly occurs in the S-wave velocity range between ∼3.1 and ∼3.6 km/s.

SUMMARYProbing seismic anisotropy of the lithosphere provides valuable clues on the fabric of rocks. We present a 3-D probabilistic model of shear wave velocity and radial anisotropy of the crust and uppermost mantle of Europe, focusing on the mountain belts of the Alps and Apennines. The model is built from Love and Rayleigh dispersion curves in the period range 5–149 s. Data are extracted from seismic ambient noise recorded at 1521 broad-band stations, including the AlpArray network. The dispersion curves are first combined in a linearized least squares inversion to obtain 2-D maps of group velocity at each period. Love and Rayleigh maps are then jointly inverted at depth for shear wave velocity and radial anisotropy using a Bayesian Monte Carlo scheme that accounts for the trade-off between radial anisotropy and horizontal layering. The isotropic part of our model is consistent with previous studies. However, our anisotropy maps differ from previous large scale studies that suggested the presence of significant radial anisotropy everywhere in the European crust and shallow upper mantle. We observe instead that radial anisotropy is mostly localized beneath the Apennines while most of the remaining European crust and shallow upper mantle is isotropic. We attribute this difference to trade-offs between radial anisotropy and thin (hectometric) layering in previous studies based on least-squares inversions and long period data (>30 s). In contrast, our approach involves a massive data set of short period measurements and a Bayesian inversion that accounts for thin layering. The positive radial anisotropy (VSH > VSV) observed in the lower crust of the Apennines cannot result from thin layering. We rather attribute it to ductile horizontal flow in response to the recent and present-day extension in the region.

SUMMARYTo constrain seismic anisotropy under and around the Alps in Europe, we study SKS shear wave splitting from the region densely covered by the AlpArray seismic network. We apply a technique based on measuring the splitting intensity, constraining well both the fast orientation and the splitting delay. Four years of teleseismic earthquake data were processed, from 723 temporary and permanent broad-band stations of the AlpArray deployment including ocean-bottom seismometers, providing a spatial coverage that is unprecedented. The technique is applied automatically (without human intervention), and it thus provides a reproducible image of anisotropic structure in and around the Alpine region. As in earlier studies, we observe a coherent rotation of fast axes in the western part of the Alpine chain, and a region of homogeneous fast orientation in the Central Alps. The spatial variation of splitting delay times is particularly interesting though. On one hand, there is a clear positive correlation with Alpine topography, suggesting that part of the seismic anisotropy (deformation) is caused by the Alpine orogeny. On the other hand, anisotropic strength around the mountain chain shows a distinct contrast between the Western and Eastern Alps. This difference is best explained by the more active mantle flow around the Western Alps. The new observational constraints, especially the splitting delay, provide new information on Alpine geodynamics.

Fracture zones were recognized to be an integral part of the seabed long before plate tectonics was established. Later, plate tectonics linked fracture zones to oceanic transform faults (OTF), suggesting that they are the inactive and hence fossil trace of transforms. Yet, scientists spend little time surveying them in much detail. Recent evidence suggests that the traditional concept of transform faults as being conservative (non-accretionary) plate boundary faults might be wrong. Instead, numerical modelling results suggest that transform faults seem to suffer from extensional tectonics below their strike-slip surface fault zone, and a global compilation of legacy bathymetric data suggest that ridge-transform intersections seem to be settings of magmatic activity, modifying the lithosphere and burying the transform valley before it passes into the fracture zone region. During M170 we tested those hypotheses by collecting a suite of new data from the Oceanographer and Hayes transform faults offsetting the Mid-Atlantic Ridge to the south of the Azores near 35°N and hence conducted a pilot study revealing the state-of-stress derived from micro-earthquakes and bathymetry as well as geological sampling to evaluate magmatic and tectonic processes shaping transform faults. Preliminary analysis of 10-days of seismicity data recorded at hydrophone stations showed 10-15 local earthquakes per day and bathymetric data supports that ridge-transform intersections support indeed a second phase of magmatic accretion.

The 2014 Mw 8.1 Iquique earthquake ruptured the boundary between the subducting Nazca Plate and the overriding South American Plate in the North Chilean subduction zone. The broken segment of the South American subduction zone had likely accumulated elastic strain since an M~9 earthquake in 1877 and what therefore considered a mature seismic gap. The moderate magnitude of the 2014 earthquake and its compact rupture area, which only broke the central part of the seismic gap, did not result in a significant tsunami in the Pacific Ocean. To investigate the seismo-tectonic segmentation of the North Chilean subduction zone in the region of the 2014 Iquique earthquake at the shallow seismic/aseismic transition, we combine two years of local aftershock seismicity observations from ocean bottom seismometers and long- offset seismic reflection data from the rupture area. Our study links short term deformation associated with a single seismic cycle to the permanent deformation history of an erosive convergent margin over millions of years. A high density of aftershocks following the 2014 Iquique earthquake occurred in the up-dip region of the coseismic rupture, where they form a trench parallel band. The events spread from the subducting oceanic plate across the plate boundary and into the overriding continental crust. The band of aftershock seismicity separates a pervasively fractured and likely fluid-filled marine forearc farther seaward from a less deformed section of the forearc farther landward. At the transition, active subduction erosion during the postseismic and possibly coseismic phases of the 2014 Iquique earthquake leads to basal abrasion of the upper plate and associated extensional faulting of the overlying marine forearc. Landward migration of the seismogenic up-dip limit, possibly at similar rates compared to the trench and the volcanic arc, leaves behind a heavily fractured and fluid-filled outermost forearc. This most seaward part of the subduction zone might be too weak to store sufficient elastic strain to nucleate a large megathrust earthquake.

Seafloor geodesy experiments have been expanding considerably in recent years. More and more research teamsaround the globe are conducting projects to monitor the tectonic or volcanic deformation of the seafloor. Theseexperiments are commonly based on limited-duration experiments, but increasingly, permanent observatories arealso being installed. This dynamic development is very encouraging for the establishment of a strong community which arguably willlead to the emergence of a worldwide scientific and technical synergy. However, data andknowledge transferbetween the different groups working on similar topics are still limited at the present time. This can be partlyexplained by the fact that the instruments, infrastructure, and processing software developed are custom-made andthus various file formats are used, although the fundamental observables are most of the time identical. One wayto overcome this limitation is to set up exchange standards in the form of standardized file formats. These fileswould gather and store all the physical quantities observed and will prove useful for the processing simplificationand, in the end, the extraction of the geodetic signal sought. Furthermore, uniformized formats would allow muchmore easily the comparison of software and processing methods between research groups, whether during tests oroperational measurement campaigns. Standardized data will eventually provide a base for the activities of potentialfuture national or international observation services. They would also make it possible to envisage the datadissemination similar to geodetic data recorded on land.