Heterogeneous seismic anisotropy beneath Madeira and Canary archipelagos revealed by local and teleseismic shear wave splitting

SUMMARY Mid-plate upward mantle flow is a key component of global mantle convection, but its patterns are poorly constrained. Seismic anisotropy is the most direct way to infer mantle flow as well as melt distribution, yet the convection patterns associated with plume-like mantle upwelling are understudied due to limited seismic data coverage. Here, we investigate seismic anisotropy beneath the Madeira and Canary hotspots using a dense set of shear wave splitting observations and combining teleseismic and local events recorded by three-component broad-band and short-period seismic stations. Using a total of 26 stations in the Madeira archipelago and 43 stations around the Canary Islands, we obtain 655 high-quality measurements that reveal heterogeneous flow patterns. Although local event results are sparse around most islands, we can observe a small average of S-wave splitting times of 0.16 ± 0.01 s, which significantly increase with source depth beneath El Hierro (>20 km) and Tenerife (>38 km) up to 0.58 ± 0.01 and 0.47 ± 0.05 s. This suggests an influence of melt pocket orientation in magma reservoirs developed at uppermost-mantle depths. Likewise, anisotropy increases significantly beneath the islands with shield stage volcanism (up to 9.81 ± 1.78 per cent at El Hierro, western Canaries, against values up to 1.76 ± 0.73 per cent at Lanzarote, eastern Canaries). On average, teleseismic SKS-wave splitting delay times are large (2.19 ± 0.05 s), indicating sublithospheric mantle flow as the primary source for anisotropy in the region. In the Canaries, the western islands show significantly smaller average SKS delay times (1.93 ± 0.07 s) than the eastern ones (2.25 ± 0.11 s), which could be explained by destructive interference above the mantle upwelling. Despite complex patterns of fast polarization directions throughout both regions, some azimuthal pattern across close stations can be observed and related to present-day mantle flow and anisotropy frozen in the lithosphere since before 60 Ma. Additionally, we infer that the current presence of a mantle plume beneath the archipelagos leads to the associated complex, small-scale heterogeneous anisotropy observations

Lateral Variation in Crustal Structure along the Lesser Antilles Arc from Petrology of Crustal Xenoliths and Seismic Receiver Functions

We reconstruct crustal structure along the Lesser Antilles island arc using an inversion approach combining constraints from petrology of magmatic crustal xenoliths and seismic receiver functions. Xenoliths show considerable island-to-island variation in xenolith petrology from plagioclase-free ultramafic lithologies to gabbros and gabbronorites with variable proportions of amphibole, indicative of changing magma differentiation depths. Xenoliths represent predominantly cumulate compositions with equilibration depths in the range 5 to 40 km. We use xenolith mineral modes and compositions to calculate seismic velocities () and density at the estimated equilibration depths. We create a five-layer model of crustal structure for testing against receiver functions (RF) from island seismic stations along the arc. Lowermost layer (5) comprises peridotite with physical characteristics of mantle xenoliths from Grenada. Uppermost layer (1) consists of 5 km of volcaniclastics and sediments, whose physical properties are determined via a grid inversion routine. The three middle layers (2) to (4) comprise igneous arc crust with compositions corresponding to the xenoliths sampled at each island. By inversion we obtain a petrological best-fit for the RF on each island to establish the nature and thicknesses of layers (2) to (4). Along the arc we see variations in the depth and strength of both Moho and mid-crustal discontinuity (MCD) on length-scales of tens of km. Moho depths vary from 25 to 37 km; MCD from 11 and 32 km. The Moho is the dominant discontinuity beneath some islands (St. Kitts, Guadeloupe, Martinique, Grenada), whereas the MCD dominates beneath others (Saba, St. Eustatius). Along-arc variability in MCD depth and strength is consistent with variation in estimated magmatic H2O contents and differentiations depths that, in turn, influence xenolith lithologies. A striking feature is steep, along-arc gradients in similar to those observed at other oceanic arcs. These gradients reflect abrupt changes in rates and processes of magma generation in the underlying crust and mantle. We find no evidence for large, interconnected bodies of partial melt beneath the Lesser Antilles. Instead, the crustal velocity structure is consistent with magma differentiation in vertically-extensive, crystal mush-dominated reservoirs. Along-arc variation in crustal structure may reflect heterogeneous upwelling within the mantle wedge, itself driven by variation in slab-derived H2O fluxes.This work was supported by Natural Environment Research Council grants NE/N001966/1, NE/K004883/1, NE/K014978/1 and NE/K010824/1 and ERC Advanced Grant CRITMAG (Blundy

Heterogeneous seismic anisotropy beneath Madeira and Canary archipelagos revealed by local and teleseismic shear wave splitting

Mid-plate upward mantle flow is a key component of global mantle convection, but its patterns are poorly constrained. Seismic anisotropy is the most direct way to infer mantle flow as well as melt distribution, yet the convection patterns associated with plume-like mantle upwelling are understudied due to limited seismic data coverage. Here, we investigate seismic anisotropy beneath the Madeira and Canary hotspots using a dense set of shear wave splitting observations and combining teleseismic and local events recorded by three-component broad-band and short-period seismic stations. Using a total of 26 stations in the Madeira archipelago and 43 stations around the Canary Islands, we obtain 655 high-quality measurements that reveal heterogeneous flow patterns. Although local event results are sparse around most islands, we can observe a small average of S-wave splitting times of 0.16 ± 0.01 s, which significantly increase with source depth beneath El Hierro (>20 km) and Tenerife (>38 km) up to 0.58 ± 0.01 and 0.47 ± 0.05 s. This suggests an influence of melt pocket orientation in magma reservoirs developed at uppermost-mantle depths. Likewise, anisotropy increases significantly beneath the islands with shield stage volcanism (up to 9.81 ± 1.78 per cent at El Hierro, western Canaries, against values up to 1.76 ± 0.73 per cent at Lanzarote, eastern Canaries). On average, teleseismic SKS-wave splitting delay times are large (2.19 ± 0.05 s), indicating sublithospheric mantle flow as the primary source for anisotropy in the region. In the Canaries, the western islands show significantly smaller average SKS delay times (1.93 ± 0.07 s) than the eastern ones (2.25 ± 0.11 s), which could be explained by destructive interference above the mantle upwelling. Despite complex patterns of fast polarization directions throughout both regions, some azimuthal pattern across close stations can be observed and related to present- day mantle flow and anisotropy frozen in the lithosphere since before 60 Ma. Additionally, we infer that the current presence of a mantle plume beneath the archipelagos leads to the associated complex, small-scale heterogeneous anisotropy observations.info:eu-repo/semantics/publishedVersio

Local seismicity around the Chain Transform Fault at the Mid-Atlantic Ridge from OBS observations

Summary Seismicity along transform faults provides important constraints for our understanding of the factors that control earthquake ruptures. Oceanic transform faults are particularly informative due to their relatively simple structure in comparison to their continental counterparts. The seismicity of several fast-moving transform faults has been investigated by local networks, but as of today there been few studies of transform faults in slow spreading ridges. Here we present the first local seismicity catalogue based on event data recorded by a temporary broadband network of 39 ocean bottom seismometers located around the slow-moving Chain Transform Fault (CTF) along the Mid-Atlantic Ridge (MAR) from March 2016 to March 2017. We locate 972 events in the area by simultaneously inverting for a 1-D velocity model informed by the event P- and S-arrival times. We refine the depths and focal mechanisms of the larger events using deviatoric moment tensor inversion. Most of the earthquakes are located along the CTF (700) and Romanche transform fault (94) and the MAR (155); a smaller number (23) can be observed on the continuing fracture zones or in intraplate locations. The ridge events are characterised by normal faulting and most of the transform events are characterised by strike slip faulting, but with several reverse mechanisms that are likely related to transpressional stresses in the region. CTF events range in magnitude from 1.1 to 5.6 with a magnitude of completeness around 2.3. Along the CTF we calculate a b-value of 0.81 ± 0.09. The event depths are mostly shallower than 15 km below sea level (523), but a small number of high-quality earthquakes (16) are located deeper, with some (8) located deeper than the brittle-ductile transition as predicted by the 600˚C-isotherm from a simple thermal model. The deeper events could be explained by the control of seawater infiltration on the brittle failure limit

Dynamic subsurface changes on El Hierro and La Palma during volcanic unrest revealed by temporal variations in seismic anisotropy patterns

Active hotspot volcanism is the surface expression of ongoing dynamic subsurface changes, such as the generation, transport, and stalling of magmas within the upper mantle and crust. Magmatic influx and migration affects local stress patterns in the crust and lithospheric mantle, which influences seismic anisotropy. A better understanding of those patterns helps improve robustness of models forecasting the likelihood of an eruption and prolonged seismicity, with detailed studies being required to observe the significant variations that can occur on small spatial and temporal scales. Here, we investigate seismic anisotropy before, during and after volcanic eruptions. We use local seismicity around El Hierro and La Palma, the two westernmost islands in the Canaries and sites of the most recent volcanic eruptions in the archipelago. We obtained 215 results in El Hierro during and after the 2011/2012 eruption with five three-component broadband seismic stations and 908 results around the 2021 eruption in La Palma with two three-component broadband stations. On La Palma, the majority of seismicity and splitting results are recorded during the eruption and simultaneous deflation of the island. Seismicity locations do not change significantly and fast shear wave polarisation direction is mostly constant, but some variation can be attributed to changes in the magmatic plumbing system. On El Hierro, the general radial pattern reflects stresses induced by the overall uplift of the island during multiple magma intrusion events. Temporal subsets reveal significant variations in location and depth of the events, as well as significant variations in fast polarisation direction caused by ongoing dynamic changes of under- and overpressurisation. An increase of results starting in 2018 hints towards renewed subsurface activity within deeper parts of the plumbing system, affecting the rate of overall seismicity but not any vertical movement of the island

Variable water input controls evolution of the Lesser Antilles volcanic arc

Oceanic lithosphere carries volatiles, notably water, into the mantle through subduction at convergent plate boundaries. This subducted water exercises control on the production of magma, earthquakes, formation of continental crust and mineral resources. Identifying different potential fluid sources (sediments, crust and mantle lithosphere) and tracing fluids from their release to the surface has proved challenging1. Atlantic subduction zones are a valuable endmember when studying this deep water cycle because hydration in Atlantic lithosphere, produced by slow spreading, is expected to be highly non-uniform2. Here, as part of a multi-disciplinary project in the Lesser Antilles volcanic arc3, we studied boron trace element and isotopic fingerprints of melt inclusions. These reveal that serpentine—that is, hydrated mantle rather than crust or sediments—is a dominant supplier of subducted water to the central arc. This serpentine is most likely to reside in a set of major fracture zones subducted beneath the central arc over approximately the past ten million years. The current dehydration of these fracture zones coincides with the current locations of the highest rates of earthquakes and prominent low shear velocities, whereas the preceding history of dehydration is consistent with the locations of higher volcanic productivity and thicker arc crust. These combined geochemical and geophysical data indicate that the structure and hydration of the subducted plate are directly connected to the evolution of the arc and its associated seismic and volcanic hazards

Local seismicity around the Chain Transform Fault at the Mid-Atlantic Ridge from OBS observations

Seismicity along transform faults provides important constraints for our understanding of the factors that control earthquake ruptures. Oceanic transform faults are particularly useful due to their relatively simple structure in comparison to continental counterparts. The seismicity of several fast-moving transform faults has been investigated by local networks, but as of today there have not been many studies of slower spreading centres. Here we present the first local seismicity catalogue based on event data recorded by a temporary broadband network of 39 ocean bottom seismometers located around the slow-moving Chain Transform Fault (CTF) along the Mid-Atlantic Ridge (MAR) from March 2016 to March 2017. Locations are constrained by simultaneously inverting for a 1-D velocity model informed by the event P- and S-arrival times. Depths and focal mechanisms of the larger events are refined using deviatoric moment tensor inversion. We find a total of 972 events in the area. Most of the seismicity is located at the CTF (700) and Romanche transform fault (94) and the MAR (155); a smaller number (23) can be observed on the continuing fracture zones or in intraplate locations. The ridge events are characterised by normal faulting and most of the transform events are characterised by strike slip faulting, but with several reverse mechanisms that are likely related to transpressional stresses in the region. CTF events range in magnitude from 1.1 to 5.6 with a magnitude of completeness around 2.3. Along the CTF we calculate a b-value of 0.81 ± 0.09. The event depths are mostly shallower than 15 km below sea level (523), but a small number of high-quality earthquakes (16) are located deeper, with some (8) located deeper than the brittle-ductile transition as predicted by the 600˚C-isotherm from a simple thermal model. The deeper events could be explained by the control of seawater infiltration on the brittle failure limit

Gaps, tears and seismic anisotropy around the subducting slabs of the Antilles

Seismic anisotropy in and beneath the subducting slabs of the Antilles is investigated using observations of shear-wave splitting. We use a combination of teleseismic and local events recorded at three-component broadband seismic stations on every major island in the area to map anisotropy in the crust, the mantle wedge and the slab/sub-slab mantle. To date this is the most comprehensive study of anisotropy in this region, involving 52 stations from 8 seismic networks. Local event delay times (0.21 ± 0.12 s) do not increase with depth, indicating a crustal origin in anisotropy and an isotropic mantle wedge. Teleseismic delay times are much larger (1.34 ± 0.47 s), with fast shear-wave polarisations that are predominantly parallel to trend of the arc. These observations can be interpreted three ways: (1) the presence of pre-existing anisotropy in the subducting slab; (2) anisotropy due to sub-slab mantle flow around the eastern margin of the nearly stationary Caribbean plate; (3) some combination of both mechanisms. However, there are two notable variations in the trench-parallel pattern of anisotropy — trench-perpendicular alignment is observed in narrow regions east of Puerto Rico and south of Martinique. These observations support previously proposed ideas of eastward sublithospheric mantle flow through gaps in the slab. Furthermore, the pattern of anisotropy south of Martinique, near Saint Lucia is consistent with a previously proposed location for the boundary between the North and South American plates

The role of the seismically slow Central-East Atlantic anomaly in the genesis of the Canary and Madeira volcanic provinces

The Canary and Madeira provinces in the Central-East Atlantic Ocean are characterized by an irregular spatio-temporal distribution of volcanism along the hotspot tracks, and several alternative scenarios have been suggested to explain it. Here, we combine results from seismic tomography, shear-wave splitting and gravity along with plate reconstruction constraints to investigate the mantle structure and dynamics beneath those provinces. We find that the Central-East Atlantic Anomaly (CEAA), which rises from the core-mantle boundary and stalls in the topmost lower mantle, is the deep source of distinct upper-mantle upwellings beneath the region. The upwellings detach intermittently from the top of the CEAA and appear to be at different evolutionary stages. We argue that the accumulation of plume material in the topmost lower mantle can play a key role in governing the first-order spatio-temporal irregularities in the distribution of hotspot volcanism.info:eu-repo/semantics/publishedVersio