SEARCH · Geology Reports

Results for “Journal of Geophysical Research-- Solid Earth”

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

SourceResearch connectionFrom yearThrough year

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

At least 577 records · Page 32Linked to original sources

USGS2012

Heat flow in vapor dominated areas of the Yellowstone Plateau volcanic field: implications for the thermal budget of the Yellowstone Caldera

Characterizing the vigor of magmatic activity in Yellowstone requires knowledge of the mechanisms and rates of heat transport between magma and the ground surface. We present results from a heat flow study in two vapor dominated, acid-sulfate thermal areas in the Yellowstone Caldera, the 0.11 km 2 Obsidian Pool Thermal Area (OPTA) and the 0.25 km 2 Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 ± 0.4 MW and 7.5 ± 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 ± 0.3 MW and 1.2 ± 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 ± 1.4 MW. These estimates result in a total heat output of 11.8 ± 1.4 MW and 8.8 ± 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. Extrapolation of the average heat flux from the OPTA (103 ± 2 W·m −2 ) and SPTA (35 ± 3 W·m −2 ) to the ~35 km 2 of vapor dominated areas in Yellowstone yields 3.6 and 1.2 GW, respectively, which is less than the total heat output transported by steam from the Yellowstone Caldera as estimated by the chloride inventory method (4.0 to 8.0 GW).

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Externally triggered renewed bubble nucleation in basaltic magma: the 12 October 2008 eruption at Halema‘uma‘u Overlook vent, Kīlauea, Hawai‘i, USA

From October 2008 until present, dozens of small impulsive explosive eruptions occurred from the Overlook vent on the southeast side of Halema‘uma‘u Crater, at Kīlauea volcano, USA. These eruptions were triggered by rockfalls from the walls of the volcanic vent and conduit onto the top of the lava column. Here we use microtextural observations and data from clasts erupted during the well-characterized 12 October 2008 explosive eruption at Halema‘uma‘u to extend existing models of eruption triggering. We present a potential mechanism for this eruption by combining microtextural observations with existing geophysical and visual data sets. We measure the size and number density of bubbles preserved in juvenile ejecta using 2D images and X-ray microtomography. Our data suggest that accumulations of large bubbles with diameters of >50 μ m to at least millimeters existed at shallow levels within the conduit prior to the 12 October 2008 explosion. Furthermore, a high number density of small bubbles <50 μ m is measured in the clasts, implying very rapid nucleation of bubbles. Visual observations, combined with preexisting geophysical data, suggest that the impact of rockfalls onto the magma free surface induces pressure changes over short timescales that (1) nucleated new additional bubbles in the shallow conduit leading to high number densities of small bubbles and (2) expanded the preexisting bubbles driving upward acceleration. The trigger of eruption and bubble nucleation is thus external to the degassing system.

Hawaii↗

USGS2012

Hydrate morphology: Physical properties of sands with patchy hydrate saturation

The physical properties of gas hydrate-bearing sediments depend on the volume fraction and spatial distribution of the hydrate phase. The host sediment grain size and the state of effective stress determine the hydrate morphology in sediments; this information can be used to significantly constrain estimates of the physical properties of hydrate-bearing sediments, including the coarse-grained sands subjected to high effective stress that are of interest as potential energy resources. Reported data and physical analyses suggest hydrate-bearing sands contain a heterogeneous, patchy hydrate distribution, whereby zones with 100% pore-space hydrate saturation are embedded in hydrate-free sand. Accounting for patchy rather than homogeneous hydrate distribution yields more tightly constrained estimates of physical properties in hydrate-bearing sands and captures observed physical-property dependencies on hydrate saturation. For example, numerical modeling results of sands with patchy saturation agree with experimental observation, showing a transition in stiffness starting near the series bound at low hydrate saturations but moving toward the parallel bound at high hydrate saturations. The hydrate-patch size itself impacts the physical properties of hydrate-bearing sediments; for example, at constant hydrate saturation, we find that conductivity (electrical, hydraulic and thermal) increases as the number of hydrate-saturated patches increases. This increase reflects the larger number of conductive flow paths that exist in specimens with many small hydrate-saturated patches in comparison to specimens in which a few large hydrate saturated patches can block flow over a significant cross-section of the specimen.

Journal of Geophysical Research B: Solid Earth↗

USGS2013

Clustering of GPS velocities in the Mojave Block, southeastern California

We find subdivisions within the Mojave Block using cluster analysis to identify groupings in the velocities observed at GPS stations there. The clusters are represented on a fault map by symbols located at the positions of the GPS stations, each symbol representing the cluster to which the velocity of that GPS station belongs. Fault systems that separate the clusters are readily identified on such a map. The most significant representation as judged by the gap test involves 4 clusters within the Mojave Block. The fault systems bounding the clusters from east to west are 1) the faults defining the eastern boundary of the Northeast Mojave Domain extended southward to connect to the Hector Mine rupture, 2) the Calico-Paradise fault system, 3) the Landers-Blackwater fault system, and 4) the Helendale-Lockhart fault system. This division of the Mojave Block is very similar to that proposed by Meade and Hager. However, no cluster boundary coincides with the Garlock Fault, the northern boundary of the Mojave Block. Rather, the clusters appear to continue without interruption from the Mojave Block north into the southern Walker Lane Belt, similar to the continuity across the Garlock Fault of the shear zone along the Blackwater-Little Lake fault system observed by Peltzer et al. Mapped traces of individual faults in the Mojave Block terminate within the block and do not continue across the Garlock Fault [Dokka and Travis, ].

California↗

USGS2018

Interseismic ground deformation and fault slip rates in the greater San Francisco Bay Area from two decades of space geodetic data

The detailed spatial variations of strain accumulation and creep on major faults in the northern San Francisco Bay Area (North Bay), which are important for seismic potential and evaluation of natural hazards, remain poorly understood. Here we combine interferometric synthetic aperture radar data from the ERS‐1/2 and Envisat satellites between 1992 and 2010 with continuous and campaign GPS data to obtain a high spatial and temporal coverage of ground deformation of the North Bay. The SAR data from both ascending and descending orbits are combined to separate horizontal and vertical components of the deformation. We jointly invert the horizontal component of the mean velocities derived from these data to infer the deep strike‐slip rates on major locked faults. We use the estimated deep rates to simulate the long‐wavelength deformation due to interseismic elastic strain accumulation along these locked faults. After removing the long‐wavelength signal from the InSAR horizontal mean velocity field, we estimate fault‐parallel surface creep rates of up to 2 mm/year along the central section of the Rodgers Creek fault and surface creep rates ranging between 2 and 4 mm/year along the Concord fault. No surface creep is geodetically resolved along the West Napa and Green Valley fault zones. We identified characteristically repeating earthquakes on the Rodgers Creek fault, the West Napa fault, the Green Valley fault, and the Concord fault. Nontectonic deformation in the Geysers geothermal field and in Late Cenozoic basins (Rohnert Park and Sonoma basins) are also observed, likely due to hydrological and sediment‐compaction processes, respectively.

Journal of Geophysical Research B: Solid Earth↗

USGS2019

Heat and mass transport in a vapor-dominated hydrothermal area in Yellowstone National Park, USA: Inferences from magnetic, electrical, electromagnetic, subsurface temperature and diffuse CO2 flux measurements

Vapor‐dominated hydrothermal systems are characterized by localized and elevated heat and gas flux. In these systems, steam and gas ascend from a boiling water reservoir, steam condenses beneath a low‐permeability cap layer, and liquid water descends, driven by gravity (“heat pipe” model). We combine magnetic, electromagnetic, and geoelectrical methods and CO 2 flux and subsurface temperature measurements in the Solfatara Plateau Thermal Area in the Yellowstone Caldera to address several fundamental questions: (1) What are the structural and/or lithological controls on heat and mass transport in vapor‐dominated areas? (2) What is the geometry and size of convecting multiphase thermal plumes? (3) Are thermal plumes associated with subsurface rock alteration and demagnetization? Magnetic and electromagnetic data inversions suggest an asymmetric 50‐ to 100‐m thick basin of glacial deposits with the thickest part adjacent to the margin of a rhyolite flow. The 3‐D electrical conductivity model in the glacial basin reveals a narrow vertical conductor interpreted as a focused multiphase plume, which coincides at the ground surface with the heat and CO 2 flux maxima. The magnetic data suggest that destruction of magnetic minerals due to rock alteration associated with the hydrothermal plume occurs mainly near the ground surface. We propose a model where the buoyant multiphase plume forms in response to decompression, boiling, and phase separation of pressurized thermal groundwater that discharges from the brecciated base of a rhyolite flow into the basin of glacial deposits. Results from multiphase groundwater flow and heat transport numerical simulations corroborate the first‐order characteristics of this model.

Wyoming↗

USGS2019

Lithospheric signature of late Cenozoic extension in electrical resistivity structure of the Rio Grande rift, New Mexico, USA

We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.

New Mexico↗

USGS2018

Induced earthquake families reveal distinctive evolutionary patterns near disposal wells

The timing of events in seismic sequences can provide insights into the physical processes controlling fault slip. In southern Kansas, the rate of earthquakes rose rapidly starting in 2013 following expansion of energy production into the area, demanding the disposal of large volumes of wastewater into deep wells. Seismicity catalogs that are complete to low magnitudes can provide insights into the physical processes that induce seismicity near wastewater disposal. We develop a catalog of over 130,000 earthquakes recorded in southern Kansas from mid‐March 2014 through December 2017 by applying a matched filter algorithm to an original catalog of 5,831 template earthquakes. Detections have nearly identical waveforms to their associated template event and represent slip on nearly co‐located sections of a fault. We select template events with at least 100 associated detections and examine the characteristics of these prolific families of earthquakes. We find that families located close (<10 km) to areas with significant volumes of injected fluids have near‐Poissonian interevent times and the families remain active over longer durations. Families farther from high‐volume injection wells show strong clustering of interevent times and shorter sequence durations. We conclude that increasing pore fluid pressures from nearby disposal of large volumes of wastewater is the primary driver of these long duration episodes, with earthquake‐earthquake interactions driving sequences at greater distance from the wells.

Kansas↗

USGS2018

Quantifying uncertainty in cumulative surface slip along the Cucamonga Fault, a crustal thrust fault in southern California

Studies of historic earthquake ground surface ruptures show that displacements along strike are spatially variable. As a result, latest Quaternary slip rates developed from a spatially restricted set of cumulative displacement measurements may not accurately represent fault velocity. Here we examine the uncertainties associated with slip on the Cucamonga Fault, which is part of a network of faults that have generated damaging historical earthquakes in and around Los Angeles, California. Numerous scarps along its ~25‐km length are well expressed on alluvial fans. We make 310 measurements of vertical separation across the scarps using lidar data. We show that the dispersion of the vertical separations cannot be explained by our best estimates of analytical uncertainties alone. Additional epistemic uncertainties are required. We find that the magnitude of the required epistemic uncertainty is typically larger than analytical uncertainty by a factor of 3 and typically about 22% of the maximum vertical separation. These relationships appear to hold at several spatial scales. We examine three potential sources of epistemic uncertainty and find that none among surface age uncertainty, fault dip, and anthropogenic landscape alteration is likely sufficient to explain the overdispersion of the data, which suggests differences in cumulative strain along the strike of the fault. We calculate a range of dip‐slip rates between 0.4 and 2.6 mm/year. In light of our results, we suggest that future thrust‐fault slip‐rate studies adopt an epistemic uncertainty of 22% of the maximum value for vertical separation measurements, unless there are sufficient data to demonstrate otherwise.

California↗

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

Characteristic earthquake magnitude frequency distributions on faults calculated from consensus data in California

An estimate of the expected earthquake rate at all possible magnitudes is needed for seismic hazard forecasts. Regional earthquake magnitude frequency distributions obey a negative exponential law (Gutenberg‐Richter), but it's unclear if individual faults do. We add three new methods to calculate long‐term California earthquake rupture rates to the existing Uniform California Earthquake Rupture Forecast (UCERF3) efforts to assess method and parameter dependence on magnitude frequency results for individual faults. All solutions show strongly characteristic magnitude‐frequency distributions on the San Andreas and other faults, with higher rates of large earthquakes than would be expected from a Gutenberg‐Richter distribution. This is a necessary outcome that results from fitting high fault slip rates under the overall statewide earthquake rate budget. We find that input data choices can affect the nucleation magnitude‐frequency distribution shape for the San Andreas fault; solutions are closer to a Gutenberg‐Richter distribution if the maximum magnitude allowed for earthquakes that occur away from mapped faults (background events) is raised above the consensus threshold of M =7.6, if the moment rate for background events is reduced, or if the overall maximum magnitude is reduced from M =8.5. We also find that participation magnitude‐frequency distribution shapes can be strongly affected by slip‐rate discontinuities along faults that may be artifacts related to segment boundaries.

California↗

USGS2019

Earthquake stress drop and Arias Intensity

We determine earthquake stress drops directly from the Arias intensity database of NGA-West2. Arias intensity (Arias, 1970) is an engineering measure proportional to the integral of the absolute value of acceleration squared, over the significant duration of the signal. As such, it is closely related to root-mean-square acceleration, and can readily be connected to earthquake stress drop (Hanks and McGuire, 1981). Arias intensity records out to 100 km yield stable stress drops for moderate-to-large magnitude earthquakes, M6.5+; for smaller events ~M4.5 – 6.5, only closer-in records yield stable results. For the 116 events considered, stress drops are about 35% larger for Class 1 mainshocks than for traditional on-fault Class 2 aftershocks, and smaller for those aftershocks close to the main fault plane. Aftershock stress drops show large variability, however, implying that on average they re-rupture weakened patches, but can also rupture intact rock or high-stress asperities. We observe an increase of stress drop with earthquake depth similar to that of other studies but do not find any significant faulting mechanism dependence. The variability of the Arias intensity-based stress drop is lower than that of eGf-based stress drops from Baltay et al. (2010, 2011), and nearly on par with variability seen in ground-motion prediction equations. The Arias intensity stress drop is a novel and promising method to estimate stress drop without the need for path and site corrections, and yields further insight into the connection between source physics and ground-motion.

Journal of Geophysical Research B: Solid Earth↗

USGS2018

Earthquakes induced by hydraulic fracturing are pervasive in Oklahoma

Wastewater disposal is generally accepted to be the primary cause of the increased seismicity rate in Oklahoma within the past decade, but no statewide analysis has investigated the contribution of hydraulic fracturing (HF) to the observed seismicity or the seismic hazard. Utilizing an enhanced seismicity catalog generated with multi‐station template matching from 2010‐2016 and all available hydraulic fracturing information, we identified 274 HF wells that are spatiotemporally correlated with bursts of seismicity. The majority of HF induced seismicity cases occurred in the SCOOP/STACK plays, but we also identified prominent cases in the Arkoma Basin as well as some more complex potential cases along the edge of the Anadarko Platform. For HF treatments where we have access to injection parameters, modeling suggests poroelastic stresses are likely responsible for seismicity, but we cannot rule out direct pore pressure effects as a contributing factor. In all of the 16 regions we identified, ≥75% of the seismicity correlated with reported HF wells. In some regions, >95% of seismicity correlated with HF wells and >50% of the HF wells correlated with seismicity. Overall, we found ~700 HF induced earthquakes with M ≥ 2.0, including 12 events with M 3.0‐3.5. These findings suggest state regulations implemented in 2018 that require operators in the SCOOP/STACK plays to take action if a M > 2 earthquake occurs could have a significant impact on future operations.

Oklahoma↗

USGS2019

Ground-motion residuals, path effects, and crustal properties: A pilot study in southern California

To improve models of ground motion estimation and probabilistic seismic hazard analyses, the engineering seismology field is moving toward developing fully nonergodic ground motion models, models specific for individual source‐to‐site paths. Previous work on this topic has examined systematic variations in ground‐motion along particular paths (from either recorded or simulated earthquake data) and has not included physical properties of the path. We present here a framework to include physical path properties, by seeking correlations between ground motion amplitudes along specific paths and crustal properties, specifically seismic velocity and anelastic attenuation, along that path. Using a large data set of small‐magnitude earthquakes recorded in Southern California, we find a correlation between the gradient of seismic S wave velocity and the path term residual, after accounting for an average geometric spreading and anelastic attenuation, indicating that heterogeneity in crustal velocity primarily controls the path‐specific attenuation. Even in aseismic regions, details of path‐specific ground motion prediction equations can be developed from crustal structure and property data.

California↗

USGS2019

Four major Holocene earthquakes on the Reelfoot fault recorded by sackungen in the New Madrid seismic zone, USA

Three sequences of well-documented, major ~M7+ earthquakes (1811-1812 CE, ~1450 CE, and ~900 CE) in the New Madrid seismic zone, USA, contribute significantly to seismic hazard in the region. However, it is unknown whether this <550 yr recurrence interval has been constant throughout the Holocene given limited geomorphic evidence of prior earthquakes. We extend the record of paleoearthquakes along the Reelfoot fault via investigation of ridge-top gravitational failure features, interpreted as sackungen. The sackungen occur in bluffs along the eastern margin of the Mississippi River floodplain and are concentrated near (<15 km) the southwest-dipping Reelfoot reverse fault. A paleoseismic trench excavated across sackungen at the Paw Paw site exposed four packages of colluvial sediment that postdate 30-11 ka Peoria loess. We interpret the colluvial packages to have been deposited following episodic failure of the sackungen as a result of strong ground motions from the following sequence of earthquakes: event 4, 1640 ± 1730 BCE; event 3, 340 ± 670 CE; event 2, 1430 ± 380 CE; and event 1, 1810 ± 50 CE (2-sigma). Event timing corresponds to previously documented earthquakes and represents the longest archive of paleoearthquakes on the Reelfoot fault. If the trenched sackungen record all major Reelfoot fault earthquakes, our observations in combination with prior investigations indicate a period of quiescence from at least 11 – 4.7 ka, followed by four major seismic events culminating in the 1811-1812 CE sequence. This clustered earthquake recurrence helps place bounds on seismic-hazard and geodynamic models in the New Madrid seismic zone.

Arkansas, Missouri↗

USGS2019

Depth determination of the 2010 El Mayor‐Cucapah earthquake sequence (M ≥ 4.0)

The 2010 M W 7.2 El Mayor‐Cucapah earthquake ruptured a zone of ~120 km in length in northern Baja California. The geographic distribution of this earthquake sequence was well constrained by waveform relocation. The depth distribution, however, was poorly determined as it is near the edge of, or outside, the Southern California Seismic Network. Here we use two complementary methods to constrain the focal depths of moderate‐sized events ( M ≥ 4.0) in this sequence. We first determine the absolute earthquake depth by modeling the regional depth phases at high frequencies (~1 Hz). We mainly focus on Pn and its depth phases pPn and sPn , which arrive early at regional distance and are less contaminated by crustal multiples. To facilitate depth phase identification and to improve signal‐to‐noise ratio, we take advantage of the dense Southern California Seismic Network and use array analysis to align and stack Pn waveforms. For events without clear depth phases, we further determine their relative depths with respect to those with known depths using differential travel times of the Pn , direct P , and direct S phases recorded for event pairs. Focal depths of 93 out of 122 M ≥ 4.0 events are tightly constrained with absolute uncertainty of about 1 km. Aftershocks are clustered in the depth range of 3–10 km, suggesting a relatively shallow seismogenic zone, consistent with high surface heat flow in this region. Most aftershocks are located outside or near the lower terminus of coseismic high‐slip patches of the main shock, which may be governed by residual strains, local stress concentration, or postseismic slip.

Baja California↗

USGS2019

Peak ground displacement saturates exactly when expected: Implications for earthquake early warning

The scaling of rupture properties with magnitude is of critical importance to earthquake early warning (EEW) systems that rely on source characterization using limited snapshots of waveform data. ShakeAlert, a prototype EEW system that is being developed for the western United States, provides real-time estimates of earthquake magnitude based on P-wave peak ground displacements measured at stations triggered by the event. The algorithms used in ShakeAlert assume that the displacement measurements at each station are statistically independent and that there exists a linear and time-independent relation between log peak ground displacement and earthquake magnitude. Here we challenge this basic assumption using a comprehensive database of more than 130,000 vertical component waveforms from M4.5-M9 earthquakes occurring near Japan from 1997 through 2017 and recorded by the K-NET and KiK-net strong-motion networks. By analyzing the time-evolution of P-wave peak ground displacements for these earthquakes, we show that there is a break, or saturation, in the magnitude-displacement scaling that depends on the length of the measurement time window. We demonstrate that the magnitude at which this saturation occurs is well-explained by a simple and non-deterministic model of earthquake rupture growth. We then use the predictions of this saturation model to develop a Bayesian framework for estimating posterior uncertainties in real-time magnitude estimates which incorporates the expected time-dependence of the peak displacement measurements.

Journal of Geophysical Research B: Solid Earth↗

USGS2019

Vertical coseismic offsets from differential high-resolution stereogrammetric DSMs: The 2013 Baluchistan, Pakistan earthquake

The recent proliferation of high-resolution (< 3-m spatial resolution) digital topography datasets opens a spectrum of geodetic applications in differential topography, including the quantification of coseismic vertical displacement fields. Most investigations of coseismic vertical displacements to date rely, in part, on pre- or post-event lidar surveys that are intractable or non-existent in many locales. Stereogrammetric digital surface models (DSMs) derived from high-resolution satellite optical imagery provide a new avenue for the retrieval of spatially-dense vertical coseismic displacements on a global scale. In this study, we generated 2-m resolution pre- and post-seismic DSMs from satellite optical imagery spanning the 2013 Mw7.7 Baluchistan strike-slip earthquake that occurred on the Hoshab fault in southern Pakistan. We applied the Iterative Closest Point algorithm to the DSMs to quantify the coseismic vertical displacement field at a spatial resolution of 10-30 m and to generate 3D coseismic strain tensors. We found that across-fault vertical offsets alternated between uplift and subsidence and varied between ~1-3 m in a non-systematic manner along the Hoshab fault. We show that the pre-existing topography and near-fault geomorphology are variably consistent and inconsistent with the displacement kinematics of the 2013 earthquake, and we argue that these relationships highlight varied slip sense history along the Hoshab fault. Notably, topography along the southern extents of the Hoshab fault requires different surface displacement kinematics than occurred in the 2013 earthquake, suggesting that the Hoshab fault accommodates varying senses of slip (bimodal slip) through time.

Baluchistan Province↗