USGS Science⌕ Search

SEARCH · USGS Science

Results for “Offshore”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,459 records · Page 81Linked to original sources
USGS2023

Volcanic and tectonic sources of seismicity near the Tanaga Volcanic Cluster, Alaska

Tanaga Island in the Central Aleutian Islands includes four stratovolcanoes: Sajaka, Tanaga, and East Tanaga in the northwest, and Takawangha in the central part of the island. Of these volcanoes, only Tanaga has a confirmed record of historical eruptive activity. We use double-difference methods to relocate Tanaga Island earthquakes from the period 2003–2017 to gain insight into volcanic and tectonic activity in the region. High precision relative relocations show subsurface structures in the Tanaga region related to volcanism and crustal tectonics. In 2005, a swarm of ∼600 volcano-tectonic events located below the NW portion of the island culminated with an episode of volcanic tremor. Although there was no known eruption associated with this swarm, we suggest that this activity is associated with fluids beneath Takawangha volcano. Overall, relatively little seismicity appears to be associated with volcanism: of 5,680 earthquakes relocated in this study, only ∼700 took place within 10 km of the volcanoes themselves. Other regions appear to experience primarily, if not exclusively, tectonic seismicity. Two regions on Tanaga Island became seismically active following a M6.6 tectonic earthquake east of the island on 2 May 2008, suggesting the triggering of tectonic earthquakes by the M6.6. Seismic activity recorded offshore and below the southern part of the island is interpreted as due to the clockwise rotation of the Delarof Block in the forearc, and bookshelf faulting along the volcanic arc. This suggests a complex pattern of earthquake hypocenters governed by both volcanic and tectonic processes surrounding Tanaga and Takawangha volcanoes.

Alaska↗
USGS2024

Modeled coastal-ocean pathways of land-sourced contaminants in the aftermath of Hurricane Florence

Extreme precipitation during Hurricane Florence, which made landfall in North Carolina in September 2018, led to breaches of hog waste lagoons, coal ash pits, and wastewater facilities. In the weeks following the storm, freshwater discharge carried pollutants, sediment, organic matter, and debris to the coastal ocean, contributing to beach closures, algae blooms, hypoxia, and other ecosystem impacts. Here, the ocean pathways of land-sourced contaminants following Hurricane Florence are investigated using the Regional Ocean Modeling System (ROMS) with a river point source with fixed water properties from a hydrologic model (WRF-Hydro) of the Cape Fear River Basin, North Carolina's largest watershed. Patterns of contaminant transport in the coastal ocean are quantified with a finite duration tracer release based on observed flooding of agricultural and industrial facilities. A suite of synthetic events also was simulated to investigate the sensitivity of the river plume transport pathways to river discharge and wind direction. The simulated Hurricane Florence discharge event led to westward (downcoast) transport of contaminants in a coastal current, along with intermittent storage and release of material in an offshore (bulge) or eastward (upcoast) region near the river mouth, modulated by alternating upwelling and downwelling winds. The river plume patterns led to a delayed onset and long duration of contaminants affecting beaches 100 km to the west, days to weeks after the storm. Maps of the onset and duration of hypothetical water quality hazards for a range of weather conditions may provide guidance to managers on the timing of swimming/shellfishing advisories and water quality sampling.

Journal of Geophysical Research: Oceans↗
USGS2023

A large sediment accretion wave along a northern California littoral cell

The northern California littoral cell of the Klamath River, which is a mixed rocky and sandy system with significant shoreline curvature, was investigated by examining ∼40 yr of satellite-derived shoreline positions and historical records. We find that an accretion wave of sediment was initiated near the Klamath River mouth in the late 1980s and translated downcoast over the subsequent decades. The wave passed rapidly (∼2,500 m/yr) through a rocky coastal reach with more oblique wave directions and slowly through a sandy reach (∼200 m/yr) where wave crests approach at more normal angles. Within the sandy reach, the accretion wave extended over 200 m offshore, was ∼10 km long, incorporated 20 ± 6 million m 3 of sediment, and averaged 1.3 ± 0.4 million m 3 /yr of longshore sediment transport over a 20-yr interval. Diffusion of the accretion wave was observed, but the diffusivity coefficient ( ε obs ∼0.01 m 2 /s) was lower than values predicted by theory, which we attribute to net sediment transport convergence in the study area caused by the curvature of the shoreline. Examining historical records, we find that increased sediment discharge in the Klamath River occurred during the 20th century from industrial-scale logging and climatic extremes. Thus, we hypothesize that increased river sediment discharge introduced new sediment to the littoral cell that initiated the observed accretion wave. These hypotheses can be tested with stratigraphic and mineralogic investigations of the broad study area beach that has formed during the past 150 years.

California↗
USGS2024

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗
USGS2023

Reconnaissance basement geology and tectonics of North Zealandia

New rock dredge samples supply key information to establish the tectonic and geological framework of the northern two-thirds of the 95% submerged Zealandia continent. The R/V Investigator voyage IN2016T01 to the Fairway Ridge, Coral Sea, obtained poorly sorted poly-lithologic pebbly to cobbly sandstones, well sorted fine grained sandstones, mudstones, bioclastic limestones, and basaltic lavas. Post-cruise analytical work comprised petrography, whole rock geochemical and Sr and Nd isotopic analyses, and U-Pb zircon, Rb-Sr, and Ar-Ar geochronology. A Fairway Ridge cobbly sandstone has a ∼95 Ma (early Late Cretaceous) depositional age; two biotite granite cobbles are 111 ± 1 and 128 ± 1 Ma in age, and some volcanic pebbles are also likely Early Cretaceous. Fairway Ridge basalts have intraplate alkaline chemistry and are of Late Eocene age (∼40–36 Ma). By analogy with South Zealandia, we interpret strong positive continental magnetic anomalies of North Zealandia to mainly result from Late Cretaceous to Cenozoic intraplate basalts, many of them rift-related lavas. A new basement geological map of North Zealandia shows the position of the Mesozoic Gondwana magmatic arc axis (Median Batholith) and other major geological units. This study completes onland and offshore reconnaissance geological mapping of the entire 5 Mkm 2 Zealandia continent.

Tectonics↗
USGS2025

Slow rupture, long rise times, and multi-fault geometry: The 2020 M6.4 southwestern Puerto Rico mainshock

The M 6.4 mainshock of the southwestern Puerto Rico seismic sequence on 7 January 2020, was one of the most impactful modern earthquakes in the northeastern Caribbean. Due to its offshore location and complex aftershock distribution, its source kinematics remain poorly constrained. This active sequence illuminated a complex set of previously unrecognized structures that indicate multiple causative faults may have slipped during its rupture. Here, we utilize seismic and geodetic observations to enhance model resolution, estimate the finite slip of the mainshock, and test a multi-segment, geologically realistic fault geometry. Our refined model finds a lower rupture velocity and longer rise times than typical for an event of this magnitude. This indicates a slow-evolving rupture process that resembles characteristics of a tsunami earthquake. Although this normal/strike-slip faulting event was not tsunamigenic, these qualities, if pervasive for this region, have important implications for future seismic monitoring and hazards in southwestern Puerto Rico.

Puerto Rico↗
USGS2025

A high-resolution 3-D P-wave velocity structure of the south-central Cascadia subduction zone from wide-angle shore-crossing seismic refraction data

This study addresses a significant gap in understanding the features of the south-central Cascadia subduction zone, a region characterized by complex geologic, tectonic, and seismic transitions both offshore and onshore. Unlike other segments along this margin, this area lacks a 3-D velocity model to delineate its structural and geological features on a fine scale. To address this void, we developed a high-resolution 3-D P-wave velocity model using active source seismic data from ship-borne seismic shots recorded on temporary and permanent onshore seismic stations and ocean-bottom seismometers. Our model shows velocity variations across the region with distinct velocity-depth profiles for the Siletz, Franciscan, and Klamath terranes in the overlying plate. We identified seaward dipping high-velocity static backstops associated with the Siletz and Klamath terranes, situated near the shoreline and further inland, respectively. Regions of reduced crustal velocity are associated with crustal faults. Moreover, there is significant along-strike depth variation in the subducting slab, which is about 4 km deeper near the thick, dense Siletz terrane and becomes shallower near the predominantly less-dense Franciscan terrane. This highlights a sudden tectonic and geologic transition at the southern boundary of the Siletz terrane. Our velocity model also indicates slightly increased hydration, though still minimal, in both the oceanic crust and the upper mantle of the subducting plate compared to other parts of the margin.

British Columbia, California, Oregon, Washington↗
USGS2026

Spatial heterogeneity of salt marsh vulnerability to sea-level rise: Dual controls of hydrological setting and salinity regime

Salt marsh vulnerability to sea-level rise (SLR) is typically assessed using point measurements of vertical accretion, neglecting three-dimensionality of geomorphic evolution and spatial variability. Recent studies suggest links between vertical and horizontal vulnerability, with differences between oligohaline and polyhaline marshes, yet these relationships remain untested in estuary-marsh systems. Here we combine geospatial analysis with hydrodynamic modeling to evaluate how unvegetated/vegetated marsh ratio (UVVR), a metric of marsh degradation, relates to elevation across hydrological regions and salinity regimes in the Albemarle-Pamlico Estuarine System, the largest lagoonal estuary in U.S. We show that at given normalized elevation, UVVR decreases across hydrological regions and salinity regimes from offshore to inland. UVVR-elevation relationship varies systematically with both hydrological setting and salinity regime, with hydrology exerting stronger influence. These findings challenge the assumption of a universal marsh deterioration trajectory and underscore the need to account for spatial heterogeneity when predicting responses to SLR.

North Carolina↗
USGS2026

Ice sheet dynamics drive pronounced changes in the subsurface freshwater-saltwater interface

Saltwater is migrating into freshwater aquifers globally with water quality and biogeochemical implications, yet saltwater intrusion in glaciated regions is sparsely investigated. Field observations suggest that groundwater head in glaciated systems is influenced by ice sheet forcings and provides evidence that seawater infiltrated into offshore aquifers during past deglaciation events. To understand links between ice sheet dynamics, groundwater head, and saltwater intrusion, we use numerical models to explore the effects of deglaciation on nearshore head and salinity distributions. We find that ice sheet thinning diminishes groundwater head, and the resulting shift in subsurface pressure gradients drives rapid landward movement of the subsurface freshwater-saltwater interface up to 4.0 km or 1.3 m per m ice sheet loss. Results highlight an overlooked saltwater intrusion mechanism that aligns with field observations and affects glaciated coastlines undergoing ice sheet retreat, underscoring the need to consider this mechanism in studies of contemporary coastal water quality.

Geophysical Research Letters↗
USGS2025

A glimpse into the future of tectonic tremor monitoring

Tectonic tremor is a weak, long-duration seismic signal often observed in subduction zones and on some other plate-bounding faults. Because of tremor's characteristically low amplitude (and low signal-to-noise) and lack of clear phase arrivals, detecting and locating tremor usually requires techniques distinct from those applied to typical earthquakes. Major advances in detection and understanding of tremor have derived in the past from a powerful combination of new data and new analysis techniques. In a recent study, Sagae et al. (2025, https://doi.org/10.1029/2025jb031348 ) exploit that combination again, developing a new machine-learning based workflow and applying it to the S-net cabled seismic network in the Japan trench offshore northern Honshu. Their approach, although complex, succeeds in detecting several times more tremor activity than earlier studies, resulting in new insights and providing a blueprint for similar approaches that could be applied elsewhere. As real-time earthquake monitoring adopts similar tools, it may present an opportunity to bring tremor monitoring into operational workflows. In turn, this could solidify tremor monitoring as a component of future operational earthquake forecasting.

JGR Solid Earth↗
USGS2012

Imaging the midcontinent rift beneath Lake Superior using large aperture seismic data

We present a detailed velocity model across the 1.1 billion year old Midcontinent Rift System (MRS) in central Lake Superior. The model was derived primarily from onshore-offshore large-aperture seismic and gravity data. High velocities obtained within a highly reflective half-graben that was imaged on coincident seismic reflection data demonstrate the dominantly mafic composition of the graben fill and constrain its total thickness to be at least 30km. Strong wide-angle reflections are observed from the lower crust and Moho, indicating that the crust is thickest (55–60km) beneath the axis of the graben. The total crustal thickness decreases rapidly to about 40 km beneath the south shore of the lake and decreases more gradually to the north. Above the Moho is a high-velocity lower crust interpreted to result from syn-rift basaltic intrusion into and/or underplating beneath the Archean lower crust. The lower crust is thickest beneath the axis of the main rift half-graben. A second region of thick lower crust is found approximately 100km north of the axis of the rift beneath a smaller half graben that is interpreted to reflect an earlier stage of rifting. The crustal model presented here resembles recent models of some passive continental margins and is in marked contrast to many models of both active and extinct Phanerozoic continental rift zones. It demonstrates that the Moho is a dynamic feature, since the pre-rift Moho is probably within or above the high-velocity lower crust, whereas the post-rift Moho is defined as the base of this layer. In the absence of major tectonic activity, however, the Moho is very stable, since the large, abrupt variations in crustal thickness beneath the MRS have been preserved for at least a billion years.

Lake Superior↗
USGS2012

On plate tectonics and the geologic evolution of southwestern North America

Very rapid subduction of the Farallon plate under southwestern North America between 60 and 40 Ma was accompanied by a relatively low volume of magmatism throughout the southwestern United States and northern Mexico. Between 40 and 20 Ma, when subduction slowed significantly and in one area may have even stopped, magmatism became widespread and voluminous from Nevada and Utah to central Mexico. This correlation of rapid subduction with a relatively low volume of magmatism can be explained by the observation that subduction-related andesitic arc volcanism, often formed in a Laramide-style compressional regime, is relatively low volume compared to continental volcanism. The shallow roots of arc volcanic systems are clearly exposed in the porphyry copper deposits found in currently active arcs and common throughout southwestern North America between 60 and 50 Ma. By 43 Ma, worldwide plate motions changed, the Pacific plate began moving away from North America, and subduction of the Farallon plate slowed. By around 36 Ma, the easternmost part of the East Pacific Rise, which was located between the Pioneer and Murray fracture zones, approached the trench and the young, hot, buoyant lithosphere appears to have clogged part of the subduction zone. Uplift on land became widespread. Voluminous continental magmatism formed the Sierra Madre Occidental (SMO) of Mexico, one of the largest batholiths in the world, as well as volcanic centers now exposed in the San Juan Mountains of Colorado and the Rio Grande Rift of New Mexico. Vectors of motion of the Pacific plate relative to the North American plate determined by Stock and Molnar (1988) are consistent with formation of a transtensional environment along the plate boundary sufficient to create a 100- to 200-km-wide void just landward of the old volcanic arc. While the SMO batholith was forming within this void, the Monterey and Arguello microplates just offshore to the west were broken off from the Farallon plate and rotated so that the East Pacific Rise in this immediate area became nearly perpendicular to the trench and perpendicular to the vector of motion of the Pacific plate relative to North America. Formation of the SMO batholith was followed between 24 and 20 Ma by a major increase in the rate of subduction of the Guadalupe plate, a fragment of the former Farallon plate, and by increasing mylonitization, extension, and uplift in the metamorphic core complexes that extend northwestward through southern Arizona from the northern end of the SMO batholith. The plate margin underwent another major change between 12.5 and 10 Ma when subduction again stopped, strike slip faulting became dominant along the coast, the Basin and Range Province opened, and numerous tectonostratigraphic terranes in southern California underwent large rotations. By 3 Ma a large, new terrane had been severed from North America immediately west of the SMO batholith as the Gulf of California opened. These observations can be explained by a model for the weakening and ultimate falling apart of the uppermost part of the subducted oceanic plate in the 20–30 m.y. after the end of rapid subduction. As the plate falls apart, not only is compressional stress relieved, but significant backslip along the old subduction zone is also possible, perhaps bringing blueschists rapidly upward from 20- to 30-km depths.

Journal of Geophysical Research↗
USGS2012

Strain accumulation in western Washington

The Juan de Fuca plate is subducted beneath the North American plate off the coast of Washington at a rate of about 40 mm/yr N68°E. The average principal strain rates (extension reckoned positive) measured in northwestern Washington are as follows: Olympic peninsula 25 km south of Port Angeles from 1982 through 1990, and and near Seattle from 1972 through 1985, and . Both strain measurements are consistent with uniaxial contraction in the direction of plate convergence. Uplift rates inferred from tide gage recordings are about 4 mm/yr on the Pacific coast and near 0 mm/yr farther inland near Seattle. These deformation rates are consistent with a model of the Cascadia subduction zone in which the plate interface beneath the continental slope and outer continental shelf is locked but free to slip farther landward. The limited downdip extent of the locked segment of the plate interface is consistent with a shallow depth (∼20 km) of the isotherm (∼450°C) that defines the brittle‐ductile transition. Small thrust events diagnostic of seismic subduction should then occur only offshore and at shallow depths. The principal strain rates measured from 1972 through 1983 in the back arc region near Richland, Washington, are and .

Washington↗
USGS2010

Gas hydrates—Geological perspective and global change

Natural gas hydrates occur worldwide in polar regions, normally associated with onshore and offshore permafrost, and in sediment of outer continental and insular margins. The total amount of methane in gas hydrates likely exceeds 10 19 g of methane carbon. Three aspects of gas hydrates are important: their fossil fuel resource potential, their role as a submarine geohazard, and their effects on global climate change. Because gas hydrates represent a large amount of methane within 2000 m of the Earth's surface, they are considered to be an unconventional, unproven source of fossil fuel. Because gas hydrates are metastable, changes of pressure and temperature affect their stability. Destabilized gas hydrates beneath the seafloor lead to geologic hazards such as submarine slumps and slides, examples of which are found worldwide. Destabilized gas hydrates may also affect climate through the release of methane, a “greenhouse” gas, which may enhance global warming and be a factor in global climate change.

Reviews of Geophysics↗
USGS2008

Aeromagnetic evidence for a buried Early Cretaceous magmatic arc, northeast Japan

Positive aeromagnetic anomalies, recent drilling, and models constructed from these data delineate the plutonic roots of the Early Cretaceous Kitakami magmatic arc in northeast Japan. Buried plutons, mostly offshore, produce belts of positive magnetic anomalies. These anomalies and magnetotelluric data suggest that the plutons form a batholith 70–120 km wide, nearly 800 km long, and 10–15 km thick. The batholith may mark the location of the main Kitakami arc. Most of the exposed Kitakami plutons are 2–20 km in diameter; some are 3 km thick. The small plutons line up along NW trending faults; some may have been satellite vents that tapped into the magma supply of the main arc. The batholithic roots of the main arc now compose almost half of the modern Japan forearc basement. Steep magnetic gradients, offset anomalies, and basin stratigraphy portray extensive faulting of the Kitakami batholith during oblique subduction in the Late Cretaceous and rifting in the Miocene. The eastern boundary of the Kitakami batholith lies between 90 and 140 km west of the modern trench, much closer than the 300-km distance between the active arc and trench. The Early Cretaceous forearc basin and accretionary prism may underlie the modern forearc basin east of the batholith, but clear evidence is lacking. Much of the Early Cretaceous margin, including most of its forearc therefore is missing. How the material was removed is unknown: it could have been strike-slip faulted, eroded by subduction-related processes, or both.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Pulling the rug out from under California: Seismic images of the Mendocino Triple Junction region

The active tectonic regime of northwestern California changes abruptly from transform motion to subduction at the Mendocino Triple Junction. Northward migration of the triple junction has been a major factor in the tectonic history of the continental margin of California since the Oligocene and continues at present. Understanding the effects of triple junction migration on the structure of the crust and upper mantle in this region is therefore necessary for reconstructing the geologic evolution of the continental margin of California and accurately assessing seismic hazards associated with the San Andreas fault system and the Cascadia subduction zone. In 1993 and 1994 a network of large-aperture seismic profiles was collected to image the crustal and upper-mantle structure beneath northern California and the adjacent continental margin. The data include approximately 650 km of onshore seismic refraction/reflection data, 2000 km of off-shore multichannel seismic (MCS) reflection data, and simultaneous onshore and offshore recording of the MCS airgun source to yield large-aperture data. Scientists from more than 12 institutions were involved in data acquisition.

California↗
USGS1996

The 1992 M=7 Cape Mendocino, California, earthquake: Coseismic deformation at the south end of the Cascadia megathrust

We invert geodetic measurements of coseismic surface displacements to determine a dislocation model for the April 25, 1992, M = 7 Cape Mendocino, California, earthquake. The orientation of the model slip vector, which nearly parallels North America-Juan de Fuca relative plate convergence, and the location and orientation of the model fault relative to the offshore Cascadia megathrust, suggest that the 1992 Cape Mendocino earthquake is the first well-recorded event to relieve strain associated with the Cascadia subduction zone. We use data from three geodetic techniques: (1) the horizontal and vertical displacements of 13 monuments surveyed with the Global Positioning System, corrected for observed horizontal interseismic strain accumulation, (2) 88 section-elevation differences between leveling monuments, and (3) the uplift of 12 coastal sites observed from the die-off of intertidal marine organisms. Maximum observed displacements are 0.4 m of horizontal movement and 1.5 m of uplift along the coast. We use Monte Carlo techniques to estimate an optimal uniform slip rectangular fault geometry and its uncertainties. The optimal model using all the data resolves 4.9 m of slip on a 14 by 15 km fault that dips 28° SE. The fault extends from 1.5 to 8.7 km in depth and the main-shock hypocenter is close to the downdip projection of the fault. The shallowly dipping fault plane is consistent with the observed aftershock locations, and the estimated geodetic moment is 3.1 × 10 19 N m, 70% of the seismic moment. Other models that exclude leveling data collected in 1935 and 1942 are more consistent with seismological estimates of the fault geometry. If the earthquake is characteristic for this segment, the estimated horizontal slip vector compared with plate convergence rates suggests a recurrence interval of 140 years, with a 95% confidence range of 100–670 years. The coseismic uplift occurred in a region that also has high Quaternary uplift rates determined from marine terrace studies. If repeated ruptures of this southernmost segment of the Cascadia megathrust are responsible for the Quaternary uplift, a comparison of the coseismic uplift with coastal uplift rates suggests a recurrence interval of 200–400 years. Thus comparing horizontal and vertical coseismic to long-term deformation suggests a recurrence interval of about 100–300 years for M = 7 events at the south end of the Cascadia megathrust.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

The July 12, 1993, Hokkaido-Nansei-Oki, Japan, earthquake: Coseismic slip pattern from strong-motion and teleseismic recordings

We employ a finite fault inversion scheme to infer the distribution of coseismic slip for the July 12, 1993 Hokkaido-Nansei-Oki earthquake using strong ground motions recorded by the Japan Meteorological Agency within 400 km of the epicenter and vertical P waveforms recorded by the Global Digital Seismograph Network at teleseismic distances. The assumed fault geometry is based on the location of the aftershock zone and comprises two fault segments with different orientations: a northern segment striking at N20°E with a 30° dip to the west and a southern segment with an N20°W strike. For the southern segment we use both westerly and easterly dip directions to test thrust orientations previously proposed for this portion of the fault. The variance reduction is greater using a shallow west dipping segment, suggesting that the direction of dip did not change as the rupture propagated south from the hypocenter. This indicates that the earthquake resulted from the shallow underthrusting of Hokkaido beneath the Sea of Japan. Static vertical movements predicted by the corresponding distribution of fault slip are consistent with the general pattern of surface deformation observed following the earthquake. Fault rupture in the northern segment accounts for about 60% of the total P wave seismic moment of 3.4 × 10 20 N m and includes a large circular slip zone (4-m peak) near the earthquake hypocenter at depths between 10 and 25 km. Slip in the southern segment is also predominantly shallower than 25 km, but the maximum coseismic displacements (2.0–2.5 m) are observed at a depth of about 5 km. This significant shallow slip in the southern portion of the rupture zone may have been responsible for the large tsunami that devastated the small offshore island of Okushiri. Localized shallow faulting near the island, however, may require a steep westerly dip to reconcile the measured values of ground subsidence.

Journal of Geophysical Research B: Solid Earth↗