Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research – Earth Surface”

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 217 records · Page 12Linked to original sources
USGS2005

Evolution of melt-vapor surface tension in silicic volcanic systems: Experiments with hydrous melts

We evaluate the melt‐vapor surface tension (σ) of natural, water‐saturated dacite melt at 200 MPa, 950–1055°C, and 4.8–5.7 wt % H 2 O. We experimentally determine the critical supersaturation pressure for bubble nucleation as a function of dissolved water and then solve for σ at those conditions using classical nucleation theory. The solutions obtained give dacite melt‐vapor surface tensions that vary inversely with dissolved water from 0.042 (±0.003) J m −2 at 5.7 wt % H 2 O to 0.060 (±0.007) J m −2 at 5.2 wt % H 2 O to 0.073 (±0.003) J m −2 at 4.8 wt % H 2 O. Combining our dacite results with data from published hydrous haplogranite and high‐silica rhyolite experiments reveals that melt‐vapor surface tension also varies inversely with the concentration of mafic melt components (e.g., CaO, FeO total , MgO). We develop a thermodynamic context for these observations in which melt‐vapor surface tension is represented by a balance of work terms controlled by melt structure. Overall, our results suggest that cooling, crystallization, and vapor exsolution cause systematic changes in σ that should be considered in dynamic modeling of magmatic processes.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Geodetic analysis of reservoir depletion at the Geyser steam field in northern California

Reservoir depletion at the Geysers from 1974 to 1977 is evident in measured changes in gravity, surface strain, and pore pressure drainage. The drainage area increased about 20%, the maximum gravity decrease was about −120 μGal, and the maximum elevation change was about 6 cm during this period. Since the net mass withdrawal is known, it may be combined with the gravity change to estimate a drainage volume. The maximum drainage volume is 25 km 3 . Because the depth of the caprock and extent of the pore pressure drainage is known, this volume limit implies that no significant mass withdrawal occurs below a depth of 4 km. The ratios of surface elevation changes to horizontal contraction coinciding with the drainage area imply an equant drainage geometry, assuming that fluid production produces negative dilatation. Using the same cylindrical geometry used to model the gravity, negative dilatational strain rates of 4 to 5 × 10 −5 /yr were found to produce the observed surface displacements. The likelihood of boiling in the system to produce steam combined with apparent large bulk moduli (from seismic velocities) and small pore pressure declines suggest that most of the reservoir contraction is due to cooling.

Journal of Geophysical Research Solid Earth↗
USGS1997

Martian drainage densities

Drainage densities on Mars range from zero over large areas of volcanic plains to 0.3–0.5 km −1 locally on some volcanoes. These values refer to geologic units, not to drainage basins, as is normal for terrestrial drainage densities. The highest values are close to the lowest terrestrial values derived by similar techniques. Drainage densities were determined for every geologic unit portrayed on the 1:15,000,000 geologic map of Mars. Except for volcanoes the geologic unit with the highest drainage density is the dissected Noachian plains with a drainage density of 0.0074 km −1 . The average drainage density for Noachian units is 0.0032 km −1 , for Hesperian units is 0.00047 km −1 , and for Amazonian units is 0.00007 km −1 , excluding the volcanoes. These values are 2–3 orders of magnitude lower than typical terrestrial densities as determined by similar techniques from Landsat images. The low drainage densities, despite a cumulative record that spans billions of years, indicate that compared with the Earth, the channel-forming processes have been very inefficient or have operated only rarely or that the surface is extremely permeable. The high drainage density on volcanoes is attributed to a local cause, such as hydrothermal activity, rather than to a global cause such as climate change.

Journal of Geophysical Research E: Planets↗
USGS2012

Deep volcanic tremor and magma ascent mechanism under Kilauea, Hawaii

Deep harmonic tremor originating at depths around 40 km under Kilauea was studied using records accumulated since 1962 at the Hawaii Volcano Observatory of the U.S. Geological Survey. The deep source of the tremor was determined by onset times and confirmed by the relative amplitude across the island-wide network of seismometers. The period of tremor was conclusively shown to be determined by the source effect and not by the path or station site effect because the period would change considerably in time but maintained uniformity across the seismic net during the tremor episode. The tremor appeared to be primarily composed of P waves. We interpret the observed period and amplitude in terms of the stationary crack model of Aki et al. (1977) and find that the seismic moment rates for deep tremors are considerably larger than those for shallow-tremors suggesting more vigorous transport for the former. We propose a kinematic source model which may be more appropriate for deep tremor. According to this model, a measurable quantity called ‘reduced displacement’ is directly proportional to the rate of magma flow. A systematic search for deep tremor episodes was made for the period from 1962 through 1979, and the amplitude, period, and duration of the tremor were tabulated. We then constructed a cumulative reduced-displacement plot over the 18-year period. The result shows a generally steady process which does not seem to be significantly affected by major eruptions and large earthquakes near the surface. The total magma flow estimated from the reduced displacement is however, one order of magnitude smaller than that estimated by Swanson (1972). It may be that most channels transport magma aseismically, and only those with strong barriers generate tremor.

Hawaii↗
USGS2007

Hydrothermal fluid flow and deformation in large calderas: Inferences from numerical simulations

[1] Inflation and deflation of large calderas is traditionally interpreted as being induced by volume change of a discrete source embedded in an elastic or viscoelastic half-space, though it has also been suggested that hydrothermal fluids may play a role. To test the latter hypothesis, we carry out numerical simulations of hydrothermal fluid flow and poroelastic deformation in calderas by coupling two numerical codes: (1) TOUGH2 [Pruess et al., 1999], which simulates flow in porous or fractured media, and (2) BIOT2 [Hsieh, 1996], which simulates fluid flow and deformation in a linearly elastic porous medium. In the simulations, high-temperature water (350°C) is injected at variable rates into a cylinder (radius 50 km, height 3–5 km). A sensitivity analysis indicates that small differences in the values of permeability and its anisotropy, the depth and rate of hydrothermal injection, and the values of the shear modulus may lead to significant variations in the magnitude, rate, and geometry of ground surface displacement, or uplift. Some of the simulated uplift rates are similar to observed uplift rates in large calderas, suggesting that the injection of aqueous fluids into the shallow crust may explain some of the deformation observed in calderas.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Geometry of the 1954 Fairview Peak-Dixie Valley earthquake sequence from a joint inversion of leveling and triangulation data

In 1954, four earthquakes greater than M s = 6.0 occurred within a 30-km radius and in a period of 6 months. Elevation and angle changes calculated from repeated leveling and triangulation surveys which span the coseismic period provide constraints on the fault geometries and coseismic slip of the faults which were activated. The quality of the coseismic geodetic data is assessed. Corrections are applied to the leveling data for subsidence due to groundwater withdrawal in the Fallon area, and a rod miscalibration error of 150 ± 30 ppm is isolated in leveling surveys made in 1967. The leveling and triangulation observations are then simultaneously inverted using the single value decomposition (SVD) inversion method to determine fault geometries and coseismic slip. Using SVD, it is possible to determine on which faults slip is resolvable given the data distribution. The faults are found to dip between 50° and 80° and extend to depths of 5 to 14 km. The geodetically derived slip values are generally equal to, or greater than, the maximum observed displacement along the surface scarps. Where slip is resolvable the geodetic data indicates the 1954 sequence contained a significant component of right-lateral slip. This is consistent with the N15°W trending shear zone which geodetic surveys have detected in western Nevada.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Variability of passive gas emissions, seismicity, and deformation during crater lake growth at White Island Volcano, New Zealand, 2002-2006

We report on 4 years of airborne measurements of CO2, SO2, and H2S emission rates during a quiescent period at White Island volcano, New Zealand, beginning in 2003. During this time a significant crater lake emerged, allowing scrubbig processes to be investigated. CO2 emissions varied from a baseline of 250 to >2000 t d-1 and demonstrated clear annual cycling that was consistent with numbers of earthquake detections and annual changes in sea level. The annual variability was found to be most likely related to increases in the strain on the volcano during sea level highs, temporarily causing fractures to reduce in size in the upper conduit. SO2 emissions varied from 0 to >400 t d-1 and were clearly affected by scrubbing processes within the first year of take development. Scrubbing caused increases of SO42- and Cl- in lake waters, and the ratio of carbon to total sulphur suggested that elemental sulphur deposition was also significant in the lake during the first year. Careful measurements of the lake level and chemistry allowed estimates of the rate of H2O(g) and HCl(g) input into the lake and suggested that the molar abundances of major gas species (H2O, CO2, SO2, and HCl) during this quiescent phase were similar to fumarolic ratios observed between earlier eruptive periods. The volume of magma estimated from CO2 emissions (0.0 15-0.04 km3) was validated by Cl- increases in the lake, suggesting that the gas and magma are transported from deep to shallow depths as a closed system and likely become open in the upper conduit region. The absence of surface deformation further leads to a necessity of magma convection to supply and remove magma from the degassing depths. Two models of convection configurations are discussed. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Effect of crustal layering upon dislocation modeling

Slip distribution at depth on a fault may be inferred from the deformation observed on the surface. In inverting the surface deformation data to obtain the slip distribution, the Earth is generally approximated by an elastic half‐space. Slip distributions inferred from a half‐space model may contain artifacts, including zones of reversed slip, due solely to effects of layering in the real Earth. This effect is demonstrated for a vertical strike‐slip fault in an Earth consisting of an elastic layer overlying an elastic half‐space. Slip on the fault is taken to be independent of the along‐strike coordinate (i.e., antiplane strain is assumed). For a given slip distribution in this model the slip distribution on a similar fault in an elastic half‐space is found that produces the identical surface deformation. Comparison of the two slip distributions reveals structure introduced into the half‐space equivalent slip profile by crustal layering. The comparisons suggest that low‐resolution inversion schemes (e.g., single screw dislocation models) are not drastically affected by Earth structure, but attempts at detailed inversion are likely to produce profiles contaminated by artifacts of Earth structure.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Timescales and processes of methane hydrate formation and breakdown, with application to geologic systems

Gas hydrate is an ice-like form of water and low molecular weight gas stable at temperatures of roughly -10ºC to 25ºC and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one-sixth of Earth’s methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean-atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean-atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth’s near-surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some fields studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

Journal of Geophysical Research-Solid Earth↗
USGS1998

Displacement field for an edge dislocation in a layered half-space

The displacement field for an edge dislocation in an Earth model consisting of a layer welded to a half-space of different material is found in the form of a Fourier integral following the method given by Weeks et al. [1968]. There are four elementary solutions to be considered: the dislocation is either in the half-space or the layer and the Burgers vector is either parallel or perpendicular to the layer. A general two-dimensional solution for a dip-slip faulting or dike injection (arbitrary dip) can be constructed from a superposition of these elementary solutions. Surface deformations have been calculated for an edge dislocation located at the interface with Burgers vector inclined 0°, 30°, 60°, and 90° to the interface for the case where the rigidity of the layer is half of that of the half-space and the Poisson ratios are the same. Those displacement fields have been compared to the displacement fields generated by similarly situated edge dislocations in a uniform half-space. The surface displacement field produced by the edge dislocation in the layered half-space is very similar to that produced by an edge dislocation at a different depth in a uniform half-space. In general, a low-modulus (high-modulus) layer causes the half-space equivalent dislocation to appear shallower (deeper) than the actual dislocation in the layered half-space.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗
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↗
USGS2012

Tectonic setting for ophiolite obduction in Oman

The Samail ophiolite is part of an elongate belt in the Middle East that forms an integral part of the Alpine mountain chains that make up the northern boundary of the Arabian-African plate. The Samail ophiolite represents a portion of the Tethyan ocean crust formed at a spreading center of Middle Cretaceous age (Cenomanian). During the Cretaceous spreading of the Tethyan Sea, Gondwana Land continued its dispersal, and the Arabian-African plate drifted northward about 10°. These events combined with the opposite rotation of Eurasia and Africa initiated the closing of the Tethyan during the Late Cretaceous. At the early stages of closure, downwarping of the Arabian continental margin combined with the compressional forces of closure from the Eurasian plate initiated obduction of the Tethyan oceanic crust along preexisting transform faults, and still hot oceanic crust was detached along oblique northeast dipping thrust faults. Amphibolites developed at the base of the detached hot peridotite as it was thrust southward over oceanic volcanic and sedimentary rocks. Plate configurations combined with palinspastic reconstructions show that subduction and attendant large-scale island arc volcanism did not commence until after the Tethyan sea began to close and after the Samail ophiolite was emplaced southward across the Arabian continental margin. The Samail ophiolite nappe now rests upon a melange consisting mainly of pelagic sediments, volcanics, and detached fragments of the basal amphibolites which in turn rest on autochthonous shelf carbonates of the Arabian platform. Laterites and conglomerates with reworked laterites on the eroded upper surface of the ophiolite indicate a period of emergence prior to the deposition of shallow water Maestrichtian carbonates. Following emplacement (Eocene) of the Samail ophiolite, the Tethyan oceanic crust began northward subduction, and active arc volcanism started just north of the present Jaz Murian depression in Iran.

Journal of Geophysical Research Solid Earth↗
USGS2012

Friction of ice

The frictional strength of ice seems to be well below that for all other rocks. Triaxial testing of cylinders of pure water ice containing a 45° inclined sawcut, at temperatures of 77 ≤ T ≤ 115 K and confining pressures 0.1 ≤ P ≤ 250 MPa, reveals the frictional laws τ = 0.20 σ n + 8.3 MPa for P ≥ 10 MPa and τ = 0.55 σ n + 1.0 MPa for P ≤ 5 MPa, where τ and σ n are the shear and normal stresses on the sawcut surface. Friction is independent of T and, over the one order of magnitude tested, of average sliding velocity. The sliding behavior is invariably stick slip. Examination of sliding surfaces shows only minor gouge development. A curious arrangement of densely packed short vertical fractures occurs primarily in samples with anomalously low strength. These results apply to models of near‐surface tectonic activity on Jupiter's three large icy moons, Europa, Callisto, and Ganymede. They indicate, for example, that a global expansion on Ganymede of 3 linear percent will cause extensional movement on preexisting faults at depths to 7±3 km.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Ground surface deformation patterns, magma supply, and magma storage at Okmok volcano, Alaska, from InSAR analysis: 2. Coeruptive deflation, July-August 2008

A hydrovolcanic eruption near Cone D on the floor of Okmok caldera, Alaska, began on 12 July 2008 and continued until late August 2008. The eruption was preceded by inflation of a magma reservoir located beneath the center of the caldera and ∼3 km below sea level (bsl), which began immediately after Okmok's previous eruption in 1997. In this paper we use data from several radar satellites and advanced interferometric synthetic aperture radar (InSAR) techniques to produce a suite of 2008 coeruption deformation maps. Most of the surface deformation that occurred during the eruption is explained by deflation of a Mogi-type source located beneath the center of the caldera and 2–3 km bsl, i.e., essentially the same source that inflated prior to the eruption. During the eruption the reservoir deflated at a rate that decreased exponentially with time with a 1/ e time constant of ∼13 days. We envision a sponge-like network of interconnected fractures and melt bodies that in aggregate constitute a complex magma storage zone beneath Okmok caldera. The rate at which the reservoir deflates during an eruption may be controlled by the diminishing pressure difference between the reservoir and surface. A similar mechanism might explain the tendency for reservoir inflation to slow as an eruption approaches until the pressure difference between a deep magma production zone and the reservoir is great enough to drive an intrusion or eruption along the caldera ring-fracture system.

Alaska↗
USGS2003

Shallow soil CO2 flow along the San Andreas and Calaveras Faults, California

We evaluate a comprehensive soil CO2 survey along the San Andreas fault (SAF) in Parkfield, and the Calaveras fault (CF) in Hollister, California, in the context of spatial and temporal variability, origin, and transport of CO2 in fractured terrain. CO2 efflux was measured within grids with portable instrumentation and continously with meteorological parameters at a fixed station, in both faulted and unfaulted areas. Spatial and temporal variability of surface CO2 effluxes was observed to be higher at faulted SAF and CF sites, relative to comparable background areas. However, ??13C (-23.3 to - 16.4???) and ??14C (75.5 to 94.4???) values of soil CO2 in both faulted and unfaulted areas are indicative of biogenic CO2, even though CO2 effluxes in faulted areas reached values as high as 428 g m-2 d-1. Profiles of soil CO2 concentration as a function of depth were measured at multiple sites within SAF and CF grids and repeatedly at two locations at the SAF grid. Many of these profiles suggest a surprisingly high component of advective CO2 flow. Spectral and correlation analysis of SAF CO2 efflux and meteorological parameter time series indicates that effects of wind speed variations on atmospheric air flow though fractures modulate surface efflux of biogenic CO2. The resulting areal patterns in CO2 effluxes could be erroneously attributed to a deep gas source in the absence of isotopic data, a problem that must be addressed in fault zone soil gas studies.

Journal of Geophysical Research B: Solid Earth↗
USGS2015

Pore-pressure sensitivities to dynamic strains: observations in active tectonic regions

Triggered seismicity arising from dynamic stresses is often explained by the Mohr-Coulomb failure criterion, where elevated pore pressures reduce the effective strength of faults in fluid-saturated rock. The seismic response of a fluid-rock system naturally depends on its hydro-mechanical properties, but accurately assessing how pore-fluid pressure responds to applied stress over large scales in situ remains a challenging task; hence, spatial variations in response are not well understood, especially around active faults. Here I analyze previously unutilized records of dynamic strain and pore-pressure from regional and teleseismic earthquakes at Plate Boundary Observatory (PBO) stations from 2006 through 2012 to investigate variations in response along the Pacific/North American tectonic plate boundary. I find robust scaling-response coefficients between excess pore pressure and dynamic strain at each station that are spatially correlated: around the San Andreas and San Jacinto fault systems, the response is lowest in regions of the crust undergoing the highest rates of secular shear strain. PBO stations in the Parkfield instrument cluster are at comparable distances to the San Andreas fault (SAF), and spatial variations there follow patterns in dextral creep rates along the fault, with the highest response in the actively creeping section, which is consistent with a narrowing zone of strain accumulation seen in geodetic velocity profiles. At stations in the San Juan Bautista (SJB) and Anza instrument clusters, the response depends non-linearly on the inverse fault-perpendicular distance, with the response decreasing towards the fault; the SJB cluster is at the northern transition from creeping-to-locked behavior along the SAF, where creep rates are at moderate to low levels, and the Anza cluster is around the San Jacinto fault, where to date there have been no statistically significant creep rates observed at the surface. These results suggest that the strength of the pore pressure response in fluid-saturated rock near active faults is controlled by shear strain accumulation associated with tectonic loading, which implies a strong feedback between fault strength and permeability: dynamic triggering susceptibilities may vary in space and also in time.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

The growth of geological structures by repeated earthquakes: 2, Field examples of continental dip-slip faults

A strong test of our understanding of the earthquake cycle is the ability to reproduce extant fault-bounded geological structures, such as basins and ranges, which are built by repeated cycles of deformation. Along strike-slip faults, the coseismic and interseismic deformation can be nearly equal in magnitude and opposite in sign, resulting in little permanent deformation except for the fault offset. For dip-slip faults, portions of the crust are lifted and dropped, and so buoyancy forces are exerted. The seismic and interseismic deformations do not balance, and structures grow and become subject to erosion and deposition. We consider three examples for which the structure and fault geometry are well known: the White Wolf reverse fault in California, site of the 1952 Kern County M=7.3 earthquake, the Lost River normal fault in Idaho, site of the 1983 Borah Peak M=7.0 earthquake, and the Cricket Mountain normal fault in Utah, site of Quaternary slip events. Basin stratigraphy and seismic reflection records are used to profile the structure, and coseismic deformation measured by leveling surveys is used to estimate the fault geometry. To reproduce these structures, we add the deformation associated with the earthquake cycle (the coseismic slip and postseismic relaxation) to the flexure caused by the observed sediment load, treating the crust as a thin elastic plate overlying a fluid substrate. The cumulative deformation is principally dependent on the elastic plate thickness, modestly sensitive to the sediment-substrate density difference, and insensitive to the fluid viscosity for the 4- to 8-Ma structures. We deduce a longterm flexural rigidity of 2–15 × 10 19 Nm; this is equivalent to an elastic plate thickness of 2–4 km for a Young's modulus of 2.5 × 10 10 Nm −2 . This value is found where independent estimates of the elastic thickness from the coherence between surface topography and gravity yield values of about 4 km, but where coseismic fault slip extends to a depth of 10–15 km. Thus much of the seismogenic crust must weaken substantially during the life of active faults, causing the fault-bounded basins to narrow over time.

Journal of Geophysical Research Solid Earth↗