Geology topics

Nathaniel C. Miller

Publications and source records attributed to Nathaniel C. Miller.

At least 19 recordsLinked to original sources

USGS2026

Improving offshore 3D splay fault geometries and slip histories using seismic data reprocessing and structural modeling

The goal of this project as written in the CRESCENT seed grant proposal was as follows: 1) reprocess selected profiles along strike from 45° to 48°N from the CASIE21 crustal-scale seismic data to obtain higher-resolution and higher-quality imaging of the uppermost 1-2 km of the accretionary wedge; 2) convert high-resolution USGS sparker seismic data from the time to depth domain to constrain near-surface fault geometries; 3) use kinematic modeling in the MOVE software to derive individual fault slip rates and per-event-displacements; and 4) work with the CFM group to create updated 3D models of identified faults based on the new data sources. At this stage, we have begun the reprocessing of the CASIE21 seismic reflection data (Carbotte et al., 2023). The time-migration reprocessing of the CASIE21 dataset has improved imaging of the near-surface structure by incorporating usable frequencies up to ~220 Hz, in contrast to the ~50 Hz maximum usable content in the currently available pre-stack depth migrated (PSDM) profiles (Fig. 1C). An example of the image quality gained through reprocessing is shown in Fig. 1, where the dominant wavelength of each reflector is reduced to ~10 m (Fig. 1B) from ~20 m in the existing PSDM product (Fig. 1A). The reprocessing work is still currently in progress. While the imaging in the shallow section has been much improved, reprocessing below the first multiple is still needed to create a seamless image from the plate boundary to the surface.

Final Report↗

USGS2025

Fiber-optic sensing for earthquake hazards research, monitoring and early warning

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

California, Oregon, Washington↗

USGS2025

Slab tearing and segmented subduction termination driven by transform tectonics

Subduction may terminate when a mid-ocean ridge approaches a trench, introducing buoyant lithosphere that resists subduction, leading to slab detachment and plate boundary reconfiguration. Yet, the spatial and temporal dynamics of slab tearing remain enigmatic due to a lack of modern examples. Here, we integrate new seismic images with regional seismicity to investigate an actively fragmenting subduction system at northern Cascadia’s ridge-trench-fault triple junction where subduction termination is imminent. Our analyses reveal a broad shear-zone initiated at ~4 Ma by exploiting ridge-parallel fabrics of nascent oceanic lithosphere and progressively localized into a mature trench-perpendicular transform boundary. This process severed an oceanic microplate and enabled its diminished subduction relative to adjacent subducting lithosphere. Downdip, we image trench-parallel slab tears offset by the transform, suggesting lateral tear propagation was intersected by the transform boundary, facilitating efficient decoupling of the microplate while allowing adjacent subduction to continue. We propose a 4D model where transform boundaries drive laterally diachronous slab fragmentation and subduction termination.

Vancouver, Washington↗

USGS2024

Late-Quaternary surface displacements on accretionary wedge splay faults in the Cascadia Subduction Zone: Implications for megathrust rupture

Because splay faults branch at a steep dip angle from the plate-boundary décollement in an accretionary wedge, their coseismic displacement can potentially result in larger tsunamis with distinct characteristics compared to megathrust-only fault ruptures, posing an enhanced hazard to coastal communities. Elsewhere, there is evidence of coseismic slip on splay faults during many of the largest subduction zone earthquakes, but our understanding of potentially active splay faults and their hazards at the Cascadia subduction zone remains limited. To identify the most recently active splay faults at Cascadia, we conduct stratigraphic and structural interpretations of near-surface deformation in the outer accretionary wedge for the ~400 km along-strike length of the landward vergence zone. We analyze recently acquired high-frequency sparker seismic data and crustal-scale multi-channel seismic data to examine the record of deformation in shallow slope basins and the upper ~1 km of the surrounding accreted sediments and to investigate linkages to deeper décollement structure. We present a new fault map for widest, most completely locked portion of Cascadia from 45 to 48°N latitude, which documents the distribution of faults that show clear evidence of recent late Quaternary activity. We find widespread evidence for active splay faulting up to 30 km landward of the deformation front, in what we define as the active domain, and diminished fault activity landward outside of this zone. The abundance of surface-deforming splay faults in the active outer wedge domain suggests Cascadia megathrust events may commonly host distributed shallow rupture on multiple splay faults located within 30 km of the deformation front.

British Columbia, California, Oregon, Washington↗

USGS2024

Methane seeps on the U.S. Atlantic margin: An updated inventory and interpretative framework

Since the discovery of >570 methane flares on the northern U.S. Atlantic margin between Cape Hatteras and Georges Bank in the last decade, the acquisition of thousands of kilometers of additional water column imaging data has provided greater coverage at water depths between the outer continental shelf and the lower continental slope. The additional high-resolution data reveal >1400 gas flares, but the removal of probable duplicates from the combined database of new flares and those recognized in 2014 yields ∼1139 unique sites. Most of these sites occur in clusters of 5 or more seeps, leaving about 275 unique locations (including 47 clusters) for seepage along the margin. As a function of depth, seep distribution is heavily skewed toward the upper continental slope at water depths shallower than 400 m on the southern New England margin and ∼ 550 m in the Mid-Atlantic Bight, with additional seeps clustered at ∼1100 m and just deeper than ∼1400 m in both sectors. Despite little ongoing tectonic deformation or active faulting on this passive margin, a variety of processes driven from below the seafloor (e.g., migration of fluids along faults or through permeable strata, seepage above diapirs or other pre-existing structures) and from above (e.g., erosion, sapping, unroofing) contribute to the development of seeps in different settings along the margin. In addition, the prevalence of seeps on promontories overlooking shelf-breaking canyons may be directly related to the three-dimensional nature of the hydrate stability zone in these locations. As a function of depth, the parts of the slope at the contemporary landward limit of gas hydrate stability are devoid of seeps, and the upper slope zones with the most concentrated seepage were not within the gas hydrate stability zone even during the Last Glacial Maximum. Thus, if the large number of upper slope seeps is at least partially sourced in gas hydrate degradation, the gas emitted at these seeps must have migrated there from greater depths on the continental slope.

Marine Geology↗

USGS2024

A scaling relationship for the width of secondary deformation around strike-slip faults

Simple mechanical arguments suggest that slip along interlocked, rough faults, damages surrounding rocks. The same arguments require that the scale of secondary damage is proportional to the size of geometric irregularities along the main fault. This relationship could apply at all scales, but has, so far, been difficult to observe at the 10s to 100 s of km scales of large, natural faults, often because large-scale deformation is distributed across wide, complex plate-boundary fault systems, like the San Andreas Fault. The geometry and geology of another large-scale plate-boundary strike slip fault—the Queen Charlotte Fault (QCF)—is, in contrast, especially simple. Here, we show that observations of secondary deformation are well-aligned with predictions of stress variations caused by geometric irregularities along the QCF, suggesting a geometric relationship between primary fault geometry and secondary deformation. The analytic stress solution reveals that the highest stresses and highest likelihood of failure are confined to a zone of influence (ZOI) with a width quantified by ZOI = λ / 2 π "> ZOI=λ/2π , where λ is the wavelength of geometric variations along the main fault. This simple model is consistent with ∼100-km-scale observations along the QCF and can theoretically be used to predict the width of secondary deformation at all scales.

Queen Charlotte Fault↗

USGS2024

Subducting plate structure and megathrust morphology from deep seismic imaging linked to earthquake rupture segmentation at Cascadia

The origin of rupture segmentation along subduction zone megathrusts and linkages to the structural evolution of the subduction zone are poorly understood. Here, regional-scale seismic imaging of the Cascadia margin is used to characterize the megathrust spanning ~900 km from Vancouver Island to the California border, across the seismogenic zone to a few tens of kilometers from the coast. Discrete domains in lower plate geometry and sediment underthrusting are identified, not evident in prior regional plate models, which align with changes in lithology and structure of the upper plate and interpreted paleo-rupture patches. Strike-slip faults in the lower plate associated with oblique subduction mark boundaries between regions of distinct lower plate geometry. Their formation may be linked to changes in upper plate structure across long-lived upper plate faults. The Juan de Fuca plate is fragmenting within the seismogenic zone at Cascadia as the young plate bends beneath the heterogeneous upper plate resulting in structural domains that coincide with paleo-rupture segmentation.

Oregon, Washington↗

USGS2023

Diving deeper into seep distribution along the Cascadia Convergent Margin, USA

Previous margin-wide studies of methane seep distribution along the Cascadia Subduction Zone indicate peaks in seep density within the landward limit of the of gas hydrate stability zone (GHSZ; ≤500 m depth), suggesting a link between current ocean warming, acceleration of hydrate dissociated, and methane emissions. This inferred connection, however, may not account for regional geologic and/or structural complexities driving methane seepage. Expanding upon an existing seep database by adding new seeps data, we conducted statistical and spatial analyses to determine margin-wide distribution trends and offer a tectonic framework for understanding the tendency toward non-normality and spatial clustering. We then highlight the role of local-scale drivers of seep formation in addition to the first-order tectonic framework, using systematic geologic/geomorphic characterization of seep emission sites in southern Cascadia and case studies using meta-attribute analysis of seismic reflection data.. Seep distribution along the margin is non-random, but instead of clustering along the 500-m isobath, regions of high seep density occur in canyons and topographic highs. New findings from this study conclude that co-location of the outer arc high (OAH) and the landward limit of the GHSZ may explain high concentrations of seeps where deformation is the greatest and hydrates are unstable. Detailed analysis of the spatial relationships between seep sites and geologic-geomorphic features in southern Cascadia reveal a link between seeps and anticlines, with 52% of the seeps found in association with anticlines, 36% found at faults, 16% associated with canyons, and 11% at seafloor failure scarps. Given that a majority of anticlines are located along or seaward of the OAH in the actively deforming outer wedge, we suggest that the location of the OAH is a primary structural control on seep distribution. This scenario is supported by neural network analysis of multichannel seismic data revealing zones of probable fluid migration along vertical pipes, faults, and chimneys in the vicinity of active seep sites on anticlines. Determining linkages between seeps and submarine tectonic geomorphology is a crucial first step for understanding and forecasting the distribution of methane seepage, but also a necessity for evaluating causal relationships between ocean warming and gas-hydrate stability.

California, Oregon, Washington↗

USGS2023

Episodic intraplate magmatism fed by a long-lived melt channel of distal plume origin

In the past decade, marine geophysical observations have led to the discovery of thin channels at the base of oceanic plates with anomalous physical properties that indicate the presence of low-degree partial melts. However, mantle melts are buoyant and should migrate toward the surface. We show abundant observations of widespread intraplate magmatism on the Cocos Plate where a thin partial melt channel was imaged at the lithosphere-asthenosphere boundary. We combine existing geophysical, geochemical, and seafloor drilling results with seismic reflection data and radiometric dating of drill cores to constrain the origin, distribution, and timing of this magmatism. Our synthesis indicates that the sublithospheric channel is a regionally extensive (>100,000 km 2 ) and long-lived feature that originated from the Galápagos Plume more than 20 Ma ago, supplying melt for multiple magmatic events and persisting today. Plume-fed melt channels may be widespread and long-lived sources for intraplate magmatism and mantle metasomatism.

Science Advances↗

USGS2022

Systematic mapping of the ocean-continent transform plate boundary of the Queen Charlotte fault system, southeastern Alaska and western British Columbia—A preliminary bathymetric terrain model

In 2015, U.S. Geological Survey scientists in collaboration with scientists from other institutions began a study of the Queen Charlotte fault—the first systematic study of the fault in more than three decades. The primary goal of the study was to gain a better understanding of the earthquake, tsunami, and underwater-landslide hazards throughout southeastern Alaska, as well as gather data to develop geologic models that can be applied to similar plate boundaries around the globe, such as the San Andreas fault system in southern California, the Alpine fault in New Zealand, and the North Anatolian fault in Turkey. A bathymetric terrain model was compiled from six different multibeam surveys of the previously unmapped Queen Charlotte fault offshore of southeastern Alaska and Haida Gwaii archipelago.

Alaska, British Columbia↗

USGS2022

Neural net detection of seismic features related to gas hydrates and free gas accumulations on the northern U.S. Atlantic margin

Bottom-simulating reflections (BSRs) that sometimes mark the base of the gas hydrate stability zone in marine sediments are often identified based on the reverse polarity reflections that cut across stratigraphic layering in seismic amplitude data. On the northern U.S. Atlantic margin (USAM) between Cape Hatteras and Hudson Canyon, legacy seismic data have revealed pronounced BSRs south of the deepwater extension of Hudson Canyon and more subtle ones from offshore Delaware south to Cape Hatteras, where the reflections sometimes follow stratigraphic layering. Using high-resolution seismic data acquired during the 2018 Mid-Atlantic Resource Imaging Experiment and a supervised neural net, we identify seismic features associated with gas hydrates and/or the top of gas between Hudson Canyon and Cape Hatteras. Using seismic attributes especially sensitive to the presence of gas, we train a neural network algorithm on seismic data from an area with strong BSRs and then apply the model to the rest of the data set. The results indicate that gas hydrate and/or shallow free gas are significantly more widespread on the northern part of the USAM than previously known. Seismic indicators of gas extend landward from the 2000 m isobath to the upper continental slope in sectors with (offshore Virginia) and, to a lesser extent, without (offshore New Jersey) pervasive upper slope methane seeps. Higher sand content and intermediate sediment thickness, factors related to the container size and gas charge in a petroleum systems framework, are associated with more robust gas indicators.

U.S. Atlantic margin↗

USGS2022

U.S. Atlantic margin gas hydrates

The minimum distribution of gas hydrates on the U.S. Atlantic margin is from offshore South Carolina northward to the longitude of Shallop Canyon on the southern New England margin. Few wells have logged or sampled the gas hydrate zone on this margin, meaning that the presence of gas hydrates is inferred primarily based on seismic data that reveal bottom simulating reflections, mostly at water depths greater than 2000 m. The highest hydrate saturations most likely exist in sandy sediments of the Whale Prospect offshore New Jersey, New York, and the western part of Cape Cod, an area characterized by strong bottom simulating reflections. Such reflections are also imaged on the well-studied Blake Ridge, where fine-grained sediments host lower hydrate saturations that have been constrained by drilling. Within the section of the margin stretching from south of Cape Hatteras to nearly Hudson Canyon, the diagnostic seismic reflections are hard to discern, making inferences about gas hydrate distributions more uncertain. The recognition of as-yet unmapped bottom simulating reflections or top of gas features seaward of the 2000 m bathymetric contour (e.g., Cape Fear Slide, Currituck slide, beneath deepwater gas seeps) within the Mid-Atlantic Bight expands the area of probable gas hydrates on this margin.

Delaware, Maryland, New Jersey, New York, North Ca↗

USGS2021

Limited mantle hydration by bending faults at the Middle America Trench

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 <100‐m in width at depths >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% H 2 0. 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.

Journal of Geophysical Research↗

USGS2020

Along-margin variations in breakup volcanism at the Eastern North American Margin

We model the magnetic signature of rift-related volcanism to understand the distribution and volumeofmagmatic activity that occurred during the breakup of Pangaea and early Atlantic opening at the Eastern North American Margin (ENAM).Along-strike variations in the amplitude and character of the prominent East Coast Magnetic Anomaly (ECMA) suggest that the emplacement of the volcanic layers producing this anomaly similarly varied along the margin. We use three-dimensional magnetic forward modeling constrained by seismic interpretationsto identify along-margin variations in volcanic thickness and width that can explain the observed amplitude and character of the ECMA. Our model results suggest that the ECMA is produced by a combination of both first-order (~600-1000 km)and second-order (~50-31100 km) magmatic segmentation. The first-order magmatic segmentation could have resulted from preexisting variations in crustal thickness and rheology developed during the tectonic amalgamation of Pangaea. The second-order magmatic segmentation developed during continental breakup and likely influenced the segmentation and transform fault spacing of the initial, and modern, Mid-Atlantic Ridge. These variations in magmatism showhow extension and thermal weakening was distributed at the ENAM during continental breakup and how this breakup magmatism was related to both previous and subsequent Wilson Cycle stages.

East Coast of United States, Atlantic Ocean↗

USGS2020

Plate boundary localization, slip-rates and rupture segmentation of the Queen Charlotte Fault based on submarine tectonic geomorphology

Linking fault behavior over many earthquake cycles to individual earthquake behavior is a primary goal in tectonic geomorphology, particularly across an entire plate boundary. Here, we examine the 1150-km-long, right-lateral Queen Charlotte-Fairweather fault system using comprehensive multibeam bathymetry data acquired along the Queen Charlotte Fault (QCF) offshore southeastern Alaska and western British Columbia. Fine-scale analysis of tectonic geomorphology allowed us to identify and reconstruct 184 strike-slip piercing points over a 630 km stretch of the QCF. Age constraints from glacial recession and offshore sedimentation patterns yield a consistent slip-rate of ∼50–57 mm/yr since ∼17–12 ka, the fastest rate for a continent-ocean strike-slip fault on Earth. These slip-rates equal or exceed estimates of Pacific-North America (PA-NA) relative motion from global plate reconstructions, indicating that PA-NA motion is highly localized. The QCF cuts the seafloor along a narrow and unusually straight trace for its entire length and multiple fault traces are observed only at local step-overs. The geometry and behavior of the QCF over many earthquake cycles is simple and typical of mature faults with relatively homogeneous stress fields. Since the QCF is the primary PA-NA plate boundary, we used the trace of the QCF to define the small circle path for relative plate motion and computed the associated Euler pole. Predicted along-strike obliquity variations based on the new pole agree with observed tectonic geomorphology and suggest that previous global plate reconstructions overestimated the degree of oblique convergence along the QCF. We also find that subtle, long-wavelength (75–150 km) bends and discrete step-overs appear to define the endpoints of M>7 earthquakes, suggesting that obliquity and resultant fault geometry may control rupture segmentation and asperity development. Lastly, the agreement between predicted obliquity and tectonic geomorphology along the entire length of QCF compelled a reevaluation of regional tectonic models. In the north, the eastern Yakatat Terrane appears to be translating northwest with the Pacific plate, and slip transferred from the QCF to the Fairweather Fault results in ∼20 mm/yr of convergence along the southern St. Elias mountains. In the south, we predict a reduced rate of convergence along the QCF west of Haida Gwaii (∼5–6 mm/yr of shortening, on average) relative to previous studies. Our results support a model for transpression and strike-slip partitioning along the edge of a hot and weak Pacific Plate, leading to crustal thickening and growth of the Queen Charlotte Terrace to the west of Haida Gwaii.

Alaska, British Columbia↗

USGS2020

Mysterious tsunami in the Caribbean Sea following the 2010 Haiti earthquake possibly generated by dynamically triggered early aftershocks

Dynamically triggered offshore aftershocks, caused by passing seismic waves from main shocks located on land, are currently not considered in tsunami warnings. The M7.0 2010 Haiti earthquake epicenter was located on land 27 km north of the Caribbean Sea and its focal mechanism was oblique strike-slip. Nevertheless, a tsunami recorded on a Caribbean Deep-Ocean Assessment and Reporting of Tsunami (DART) buoy and a tide gauge produced runup heights of 1–3 m along Haiti southeast coast. Earthquake finite-fault model inversions of the DART waveform suggest that a reverse fault doublet with magnitudes of M6.8 and M6.5 located 85 km southwest of the epicenter may have excited the tsunami. This doublet collocates with dynamically triggered aftershocks, derived from back-projection analysis, that occurred 20-60 s after the main shock of the Haiti earthquake. The aftershocks are within a region of maximum dynamic strain predicted by the main shock, on a possibly tectonically active submarine ridge southwest of Haiti's Southern Peninsula. The agreement between the tsunami finite-fault source models and the seismic and tectonic evidence suggests that earthquakes on land, even strike-slip faults, can generate tsunamis by dynamically triggering offshore aftershocks.

Earth and Planetary Science Letters↗

USGS2019

Azimuthal seismic anisotropy of 70 Ma Pacific‐plate upper mantle

Plate formation and evolution processes are predicted to generate upper mantle seismic anisotropy and negative vertical velocity gradients in oceanic lithosphere. However, predictions for upper mantle seismic velocity structure do not fully agree with the results of seismic experiments. The strength of anisotropy observed in the upper mantle varies widely. Further, many refraction studies observe a fast direction of anisotropy rotated several degrees with respect to the paleospreading direction, suggesting that upper mantle anisotropy records processes other than 2D corner flow and plate‐driven shear near mid‐ocean ridges. We measure 6.0 ± 0.3% anisotropy at the Moho in 70 Ma lithosphere in the central Pacific with a fast direction parallel to paleospreading, consistent with mineral alignment by 2D mantle flow near a mid‐ocean ridge. We also find an increase in the strength of anisotropy with depth, with vertical velocity gradients estimated at 0.02 km/s/km in the fast direction and 0 km/s/km in the slow direction. The increase in anisotropy with depth can be explained by mechanisms for producing anisotropy other than intrinsic effects from mineral fabric, such as aligned cracks or other structures. This measurement of seismic anisotropy and gradients reflects the effects of both plate formation and evolution processes on seismic velocity structure in mature oceanic lithosphere, and can serve as a reference for future studies to investigate the processes involved in lithospheric formation and evolution.

Journal of Geophysical Research B: Solid Earth↗

USGS2018

Practical approaches to maximizing the resolution of sparker seismic reflection data

Sparkers are a type of sound source widely used by the marine seismic community to provide high-resolution imagery of the shallow sub-bottom (i.e., < 1000 m). Although sparkers are relatively simple, inexpensive, and high-frequency (100–2500 Hz) sources, they have several potential pitfalls due to their complicated and unpredictable signature. In this study we quantify the source characteristics of several sparker systems and develop a suite of simple processing approaches for both single channel and multi-channel sparker data. In all cases, the results show improved vertical resolution and reflection coherency. Correcting for small static variations in multi-channel seismic (MCS) data is a critical first step to preserve the broad frequency content during stacking, and to reduce the shot-to-shot variability of outgoing and incoming signals. Application of predictive deconvolution to static-corrected, post-stack traces suppresses short-path multiples and restores the latent high-resolution reflection patterns. However, if shot-to-shot source signatures are recorded directly, pre-stack deterministic deconvolution followed by post-stack predictive deconvolution produces the most robust results. Processing sparker data without broadband techniques results in less confident or completely missed interpretations when compared to the broadband equivalent. If processed correctly, marine sparker data can provide exceptional sub-bottom imagery that rivals other more repeatable marine seismic sources (e.g., high-frequency air-guns).

Marine Geophysical Research↗