A Bayesian 3-D linear gravity inversion for complex density distributions: application to the Puysegur subduction system

We have developed a linear 3-D gravity inversion method capable of modelling complex geological regions such as subduction margins. Our procedure inverts satellite gravity to determine the best-fitting differential densities of spatially discretized subsurface prisms in a least-squares sense. We use a Bayesian approach to incorporate both data error and prior constraints based on seismic reflection and refraction data. Based on these data, Gaussian priors are applied to the appropriate model parameters as absolute equality constraints. To stabilize the inversion and provide relative equality constraints on the parameters, we utilize a combination of first and second order Tikhonov regularization, which enforces smoothness in the horizontal direction between seismically constrained regions, while allowing for sharper contacts in the vertical. We apply this method to the nascent Puysegur Trench, south of New Zealand, where oceanic lithosphere of the Australian Plate has underthrust Puysegur Ridge and Solander Basin on the Pacific Plate since the Miocene. These models provide insight into the density contrasts, Moho depth, and crustal thickness in the region. The final model has a mean standard deviation on the model parameters of about 17 kg m⁻³, and a mean absolute error on the predicted gravity of about 3.9 mGal, demonstrating the success of this method for even complex density distributions like those present at subduction zones. The posterior density distribution versus seismic velocity is diagnostic of compositional and structural changes and shows a thin sliver of oceanic crust emplaced between the nascent thrust and the strike slip Puysegur Fault. However, the northern end of the Puysegur Ridge, at the Snares Zone, is predominantly buoyant continental crust, despite its subsidence with respect to the rest of the ridge. These features highlight the mechanical changes unfolding during subduction initiation

CAYSEIS - magma-starved oceanic crustal accretion and transform margin formation in the Cayman Trough revealed by seismic and seismological data - Cruise No. M115, April 1 - April 28, 2015 - Kingston (Jamaica) - Pointe-à-Pitre (Guadeloupe)

About 57% of the Earth’s outer surface is oceanic crust and new ocean floor is continuously created along the 55,000-60,000 km long mid-ocean ridge (MOR) system. About 25% of MORs spread at an ultra-slow spreading rate of < 20 mm/yr. Most ultra-slow spreading ridges occur in areas of the world that are difficult to reach, like the Gakkel Ridge in the Arctic Ocean and the Southwest Indian Ridge in the Indian Ocean. It has long been recognized that crustal accretion at ultra-slow spreading rates is fundamentally different from crust generated at faster spreading rates. However, due to the remoteness of ultra-slow ridges, the formation of crust at these magma-starved centres is yet not well understood. During the CAYSEIS cruise we surveyed lithospheric formation at ultra-slow spreading rates at the Mid-Cayman spreading centre (MCSC) in the Caribbean Sea, where oceanic crust is formed at a full rate of ~17 mm/yr. To the northeast and southwest, the MCSC is bound by two major transform faults. Using active-source wide-angle seismic imaging and passive local earthquake monitoring we, studied the balance between magmatic accretion and tectonic stretching (and hence oceanic core complex formation) and the relationship between faulting and hydrothermal activity at ultra-slow spreading rates. In addition, we explored transform margin formation at a unique setting, occurring at the southern terminus of the MCSC. In total, six seismic lines surveyed crust formed at the MCSC, two of these profiles also crossed the Swan Island transform fault. The project was a collaboration between German, British and American scientists

Along-strike structure of the Costa Rican convergent margin from seismic a refraction/reflection survey : evidence for underplating beneath the inner forearc

Author Posting. © American Geophysical Union, 2016. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 17 (2016): 501–520, doi:10.1002/2015GC006029.The convergent margin offshore Costa Rica shows evidence of subsidence due to subduction erosion along the outer forearc and relatively high rates of uplift (∼3–6 mm/yr) along the coast. Recently erupted arc lavas exhibit a low 10Be signal, suggesting that although nearly the entire package of incoming sediments enters the subduction zone, very little of that material is carried directly with the downgoing Cocos plate to the magma generating depths of the mantle wedge. One mechanism that would explain both the low 10Be and the coastal uplift is the underplating of sediments, tectonically eroded material, and seamounts beneath the inner forearc. We present results of a 320 km long, trench-parallel seismic reflection and refraction study of the Costa Rican forearc. The primary observations are (1) margin perpendicular faulting of the basement, (2) thickening of the Cocos plate to the northwest, and (3) two weak bands of reflections in the multichannel seismic (MCS) reflection image with travel times similar to the top of the subducting Cocos plate. The modeled depths to these reflections are consistent with an ∼40 km long, 1–3 km thick region of underplated material ∼15 km beneath some of the highest observed coastal uplift rates in Costa Rica.This work was funded by the U. S. National Science Foundation MARGINS program under grants OCE 0405556, OCE 0405654, and OCE 0625178.2016-08-2

Limited mantle hydration by bending faults at the Middle America Trench

Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(1),(2021): e2020JB020982, https://doi.org/10.1029/2020JB020982.Seismic anisotropy measurements show that upper mantle hydration at the Middle America Trench (MAT) is limited to serpentinization and/or water in fault zones, rather than distributed uniformly. Subduction of hydrated oceanic lithosphere recycles water back into the deep mantle, drives arc volcanism, and affects seismicity at subduction zones. Constraining the extent of upper mantle hydration is an important part of understanding many fundamental processes on Earth. Substantially reduced seismic velocities in tomography suggest that outer rise plate‐bending faults provide a pathway for seawater to rehydrate the slab mantle just prior to subduction. Estimates of outer‐rise hydration based on tomograms vary significantly, with some large enough to imply that, globally, subduction has consumed more than two oceans worth of water during the Phanerozoic. We found that, while the mean upper mantle wavespeed is reduced at the MAT outer rise, the amplitude and orientation of inherited anisotropy are preserved at depths >1 km below the Moho. At shallower depths, relict anisotropy is replaced by slowing in the fault‐normal direction. These observations are incompatible with pervasive hydration but consistent with models of wave propagation through serpentinized fault zones that thin to 1 km below Moho. Confining hydration to fault zones reduces water storage estimates for the MAT upper mantle from ∼3.5 wt% to <0.9 wt% H20. Since the intermediate thermal structure in the ∼24 Myr‐old MAT slab favors serpentinization, limited hydration suggests that fault mechanics are the limiting factor, not temperatures. Subducting mantle may be similarly dry globally.National Science Foundation. Grant Numbers: OCE-0625178, OCE-08410632021-06-1

THE ROLE OF FRICTIONAL STRENGTH

[1] At a subduction zone the amount of friction between the incoming plate and the forearc is an important factor in controlling the dip angle of subduction and the structure of the forearc. In this paper, we investigate the role of the frictional strength of sediments and of the serpentinized peridotite on the evolution of convergent margins. In numerical models, we vary thickness of a serpentinized layer in the mantle wedge (15 to 25 km) and the frictional strength of both the sediments and serpentinized mantle (friction angle 1° to 15°, or static friction coefficient 0.017 to 0.27) to control the amount of frictional coupling between the plates. With plastic strain weakening in the lithosphere, our numerical models can attain stable subduction geometry over millions of years. We find that the frictional strength of the sediments and serpentinized peridotite exerts the largest control on the dip angle of the subduction interface at seismogenic depths. In the case of low sediment and serpentinite friction, the subduction interface has a shallow dip, while the subduction zone develops an accretionary prism, a broad forearc high, a deep forearc basin, and a shallow trench. In the high friction case, the subduction interface is steep, the trench is deeper, and the accretionary prism, forearc high and basin are all absent. The resultant free-air gravity and topographic signature of these subduction zone models are consistent with observations. We believe that the low-friction model produces a geometry and forearc structure similar to that of accretionary margins. Conversely, models with high friction angles in sediments and serpentinite develop characteristics of an erosional convergent margin. We find that the strength of the subduction interface is critical in controlling the amount of coupling at the seismogenic zone and perhaps ultimately the size of the largest earthquakes at subduction zones

A comparison between the transpressional plate boundaries of South Island, New Zealand, and Southern California, USA: the Alpine and San Andreas fault systems

There are clear similarities in structure and tectonics between the Alpine Fault system (AF) of New Zealand’s South Island and the San Andreas Fault system (SAF) of southern California, USA. Both systems are transpressional, with similar right slip and convergence rates, similar onset ages (for the current traces), and similar total offsets. There are also notable differences, including the dips of the faults and their plate-tectonic histories. The crustal structure surrounding the AF and SAF was investigated with active and passive seismic sources along transects known as South Island Geophysical Transect (SIGHT) and Los Angeles Region Seismic Experiment (LARSE), respectively. Along the SIGHT transects, the AF appears to dip moderately southeastward (~50 deg.), toward the Pacific plate (PAC), but along the LARSE transects, the SAF dips vertically to steeply northeastward toward the North American plate (NAM). Away from the LARSE transects, the dip of the SAF changes significantly. In both locations, a midcrustal decollement is observed that connects the plate-boundary fault to thrust faults farther south in the PAC. This decollement allows upper crust to escape collision laterally and vertically, but forces the lower crust to form crustal roots, reaching maximum depths of 44 km (South Island) and 36 km (southern California). In both locations, upper-mantle bodies of high P velocity are observed extending from near the Moho to more than 200-km depth. These bodies appear to be confined to the PAC and to represent oblique downwelling of PAC mantle lithosphere along the plate boundaries

Strike-slip Enables Subduction Initiation beneath a Failed Rift: New Seismic Constraints from Puysegur Margin, New Zealand

Subduction initiation often takes advantage of previously weakened lithosphere and may preferentially nucleate along pre-existing plate boundaries. To evaluate how past tectonic regimes and inherited lithospheric structure might lead to self-sustaining subduction, we present an analysis of the Puysegur Trench, a young subduction zone with a rapidly evolving tectonic history. The Puysegur margin, south of New Zealand, has experienced a transformation from rifting to seafloor spreading to strike-slip, and most recently to incipient subduction, all in the last ~45 million years. Here we present deep-penetrating multichannel reflection (MCS) and ocean-bottom seismometer (OBS) tomographic images to document crustal structures along the margin. Our images reveal that the overriding Pacific Plate beneath the Solander Basin contains stretched continental crust with magmatic intrusions, which formed from Eocene-Oligocene rifting between the Campbell and Challenger plateaus. Rifting was more advanced to the south, yet never proceeded to breakup and seafloor spreading in the Solander Basin as previously thought. Subsequent strike-slip deformation translated continental crust northward causing an oblique collisional zone, with trailing ~10 Myr old oceanic lithosphere. Incipient subduction transpired as oceanic lithosphere from the south forcibly underthrust the continent-collision zone. We suggest that subduction initiation at the Puysegur Trench was assisted by inherited buoyancy contrasts and structural weaknesses that were imprinted into the lithosphere during earlier phases of continental rifting and strike-slip along the plate boundary. The Puysegur margin demonstrates that forced nucleation along a strike-slip boundary is a viable subduction initiation scenario and should be considered throughout Earth's history

Strike-Slip Enables Subduction Initiation Beneath a Failed Rift: New Seismic Constraints From Puysegur Margin, New Zealand

Subduction initiation often takes advantage of previously weakened lithosphere and may preferentially nucleate along pre‐existing plate boundaries. To evaluate how past tectonic regimes and inherited lithospheric structure might lead to self‐sustaining subduction, we present an analysis of the Puysegur Trench, a young subduction zone with a rapidly evolving tectonic history. The Puysegur margin, south of New Zealand, has experienced a transformation from rifting to seafloor spreading to strike‐slip, and most recently to incipient subduction, all in the last ∼45 million years. Here we present deep‐penetrating multichannel reflection and ocean‐bottom seismometer tomographic images to document crustal structures along the margin. Our images reveal that the overriding Pacific Plate beneath the Solander Basin contains stretched continental crust with magmatic intrusions, which formed from Eocene‐Oligocene rifting between the Campbell and Challenger plateaus. Rifting was more advanced to the south, yet never proceeded to breakup and seafloor spreading in the Solander Basin as previously thought. Subsequent strike‐slip deformation translated continental crust northward causing an oblique collisional zone, with trailing ∼10 Myr old oceanic lithosphere. Incipient subduction transpired as oceanic lithosphere from the south forcibly underthrust the continent‐collision zone. We suggest that subduction initiation at the Puysegur Trench was assisted by inherited buoyancy contrasts and structural weaknesses that were imprinted into the lithosphere during earlier phases of continental rifting and strike‐slip along the plate boundary. The Puysegur margin demonstrates that forced nucleation along a strike‐slip boundary is a viable subduction initiation scenario and should be considered throughout Earth's history

A comparison between the transpressional plate boundaries of the South Island, New Zealand, and southern California, USA: the Alpine and San Andreas Fault systems

There are clear similarities in structure and tectonics between the Alpine Fault system (AF) of New Zealand’s South Island and the San Andreas Fault system (SAF) of southern California, USA. Both systems are transpressional, with similar right slip and convergence rates, similar onset ages (for the current traces), and similar total offsets. There are also notable differences, including the dips of the faults and their plate-tectonic histories. The crustal structure surrounding the AF and SAF was investigated with active and passive seismic sources along transects known as South Island Geophysical Transect (SIGHT) and Los Angeles Region Seismic Experiment (LARSE), respectively. Along the SIGHT transects, the AF appears to dip moderately southeastward (~50 deg.), toward the Pacific plate (PAC), but along the LARSE transects, the SAF dips vertically to steeply northeastward toward the North American plate (NAM). Away from the LARSE transects, the dip of the SAF changes significantly. In both locations, a midcrustal decollement is observed that connects the plate-boundary fault to thrust faults farther south in the PAC. This decollement allows upper crust to escape collision laterally and vertically, but forces the lower crust to form crustal roots, reaching maximum depths of 44 km (South Island) and 36 km (southern California). In both locations, upper-mantle bodies of high P velocity are observed extending from near the Moho to more than 200-km depth. These bodies appear to be confined to the PAC and to represent oblique downwelling of PAC mantle lithosphere along the plate boundaries