USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Planets”

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

Models of lithosphere and asthenosphere anisotropic structure of the Yellowstone hot spot from shear wave splitting

Teleseismic shear wave splitting measured at 56 continuous and temporary seismographs deployed in a 500 km by 600 km area around the Yellowstone hot spot indicates that fast anisotropy in the mantle is parallel to the direction of plate motion under most of the array. The average split time from all stations of 0.9 s is typical of continental stations. There is little evidence for plume-induced radial strain, suggesting that any contribution of gravitationally spreading plume material is undetectably small with respect to the plate motion velocity. Two stations within Yellowstone have splitting measurements indicating the apparent fast anisotropy direction (ϕ) is nearly perpendicular to plate motion. These stations are ∼30 km from stations with ϕ parallel to plate motion. The 70° rotation over 30 km suggests a shallow source of anisotropy; however, split times for these stations are more than 2 s. We suggest melt-filled, stress-oriented cracks in the lithosphere are responsible for the anomalous ϕ orientations within Yellowstone. Stations southeast of Yellowstone have measurements of ϕ oriented NNW to WNW at high angles to the plate motion direction. The Archean lithosphere beneath these stations may have significant anisotropy capable of producing the observed splitting.

Yellowstone↗
USGS2005

Waveform tomography of crustal structure in the south San Francisco Bay region

We utilize a scattering-based seismic tomography technique to constrain crustal tructure around the southern San Francisco Bay region (SFBR). This technique is based on coupled traveling wave scattering theory, which has usually been applied to the interpretation of surface waves in large regional-scale studies. Using fully three-dimensional kernels, this technique is here applied to observed P, S, and surface waves of intermediate period (3-4 s dominant period) observed following eight selected regional events. We use a total of 73 seismograms recorded by a U.S. Geological Survey short-period seismic array in the western Santa Clara Valley, the Berkeley Digital Seismic Network, and the Northern California Seismic Network. Modifications of observed waveforms due to scattering from crustal structure include (positive or negative) amplification, delay, and generation of coda waves. The derived crustal structure explains many of the observed signals which cannot be explained with a simple layered structure. There is sufficient sensitivity to both deep and shallow crustal structure that even with the few sources employed in the present study, we obtain shallow velocity structure which is reasonably consistent with previous P wave tomography results. We find a depth-dependent lateral velocity contrast across the San Andreas fault (SAF), with higher velocities southwest of the SAF in the shallow crust and higher velocities northeast of the SAF in the midcrust. The method does not have the resolution to identify very slow sediment velocities in the upper approximately 3 km since the tomographic models are smooth at a vertical scale of about 5 km. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Source mechanism of Vulcanian degassing at Popocatépetl Volcano, Mexico, determined from waveform inversions of very long period signals

The source mechanism of very long period (VLP) signals accompanying volcanic degassing bursts at Popocatépetl is analyzed in the 15–70 s band by minimizing the residual error between data and synthetics calculated for a point source embedded in a homogeneous medium. The waveforms of two eruptions (23 April and 23 May 2000) representative of mild Vulcanian activity are well reproduced by our inversion, which takes into account volcano topography. The source centroid is positioned 1500 m below the western perimeter of the summit crater, and the modeled source is composed of a shallow dipping crack (sill with easterly dip of 10°) intersecting a steeply dipping crack (northeast striking dike dipping 83° northwest), whose surface extension bisects the vent. Both cracks undergo a similar sequence of inflation, deflation, and reinflation, reflecting a cycle of pressurization, depressurization, and repressurization within a time interval of 3–5 min. The largest moment release occurs in the sill, showing a maximum volume change of 500–1000 m 3 , pressure drop of 3–5 MPa, and amplitude of recovered pressure equal to 1.2 times the amplitude of the pressure drop. In contrast, the maximum volume change in the dike is less (200–300 m 3 ), with a corresponding pressure drop of 1–2 MPa and pressure recovery equal to the pressure drop. Accompanying these volumetric sources are single-force components with magnitudes of 10 8 N, consistent with melt advection in response to pressure transients. The source time histories of the volumetric components of the source indicate that significant mass movement starts within the sill and triggers a mass movement response in the dike within a few seconds. Such source behavior is consistent with the opening of a pathway for escape of pent-up gases from slow pressurization of the sill driven by magma crystallization. The opening of this pathway and associated rapid evacuation of volcanic gases induces the pressure drop. Pressure recovery in the magma filling the sill is driven by diffusion of gases from the resulting supersaturated melt into bubbles. Assuming a penny-shaped crack at ambient pressure of 40 MPa, the observed pressure and volume variations can be modeled with the following attributes: crack radius (100 m), crack aperture (5 m), bubble number density (10 10 –10 12 m −3 ), initial bubble radius (10 −6 m), final bubble radius (∼10 −5 m), and net decrease of gas concentration in the melt (0.01 wt %).

Popocatépetl Volcano↗
USGS2005

Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution-dependent limitations

[1] Geophysical imaging has traditionally provided qualitative information about geologic structure; however, there is increasing interest in using petrophysical models to convert tomograms to quantitative estimates of hydrogeologic, mechanical, or geochemical parameters of interest (e.g., permeability, porosity, water content, and salinity). Unfortunately, petrophysical estimation based on tomograms is complicated by limited and variable image resolution, which depends on (1) measurement physics (e.g., electrical conduction or electromagnetic wave propagation), (2) parameterization and regularization, (3) measurement error, and (4) spatial variability. We present a framework to predict how core‐scale relations between geophysical properties and hydrologic parameters are altered by the inversion, which produces smoothly varying pixel‐scale estimates. We refer to this loss of information as “correlation loss.” Our approach upscales the core‐scale relation to the pixel scale using the model resolution matrix from the inversion, random field averaging, and spatial statistics of the geophysical property. Synthetic examples evaluate the utility of radar travel time tomography (RTT) and electrical‐resistivity tomography (ERT) for estimating water content. This work provides (1) a framework to assess tomograms for geologic parameter estimation and (2) insights into the different patterns of correlation loss for ERT and RTT. Whereas ERT generally performs better near boreholes, RTT performs better in the interwell region. Application of petrophysical models to the tomograms in our examples would yield misleading estimates of water content. Although the examples presented illustrate the problem of correlation loss in the context of near‐surface geophysical imaging, our results have clear implications for quantitative analysis of tomograms for diverse geoscience applications.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Spatial and temporal multiyear sea ice distributions in the Arctic: A neural network analysis of SSM/I data, 1988-2001

Arctic multiyear sea ice concentration maps for January 1988-2001 were generated from SSM/I brightness temperatures (19H, 19V, and 37V) using modified multiple layer perceptron neural networks. Learning data for the neural networks were extracted from ice maps derived from Okean and ERS satellite imagery to capitalize on the stability of active radar multiyear ice signatures. Evaluations of three learning algorithms and several topologies indicated that networks constructed with error back propagation learning and 3-20-1 topology produced the most consistent and physically plausible results. Operational neural networks were developed specifically with January learning data, and then used to estimate daily multiyear ice concentrations from daily-averaged SSM/I brightness temperatures during January. Monthly mean maps were produced for analysis by averaging the respective daily estimates. The 14-year series of January multiyear ice distributions revealed dense and persistent cover in the central Arctic surrounded by expansive regions of highly fluctuating interannual cover. Estimates of total multiyear ice area by the neural network were intermediate to those of other passive microwave algorithms, but annual fluctuations and trends were similar among all algorithms. When compared to Radarsat estimates of multiyear ice concentration in the Beaufort and Chukchi Seas (1997-1999), average discrepancies were small (0.9-2.5%) and spatial coherency was reasonable, indicating the neural network's Okean and ERS learning data facilitated passive microwave inversion that emulated backscatter signatures. During 1988-2001, total January multiyear ice area declined at a significant linear rate of -54.3 x 10 3 km 2 /yr -1 (-1.4%/yr -1 ). The most persistent and extensive decline in multiyear ice concentration (-3.3%/yr -1 ) occurred in the southern Beaufort and Chukchi Seas. In autumn 1996, a large multiyear ice recruitment of over 106 km 2 (mostly in the Siberian Arctic) fully replenished the previous 8-year decline in total area, but it was followed by an accelerated and compensatory decline during the subsequent 4 years. Seventy-five percent of the interannual variation in January multiyear sea ice area was explained by linear regression on two atmospheric parameters: the previous inter's (JFM) Arctic Oscillation index as a proxy to melt duration and the previous year's average sea level pressure gradient across the Fram Strait as a proxy to annual ice export. Consecutive year changes (1994-2001) in January multiyear ice volume were significantly correlated with duration of the intervening melt season (R 2 = 0.73, -80.0 km 3 /d -1 ), emphasizing a large thermodynamic influence on the Arctic's mass sea ice balance during summers with anomalous melt durations.

Journal of Geophysical Research↗
USGS2005

Wave spectral energy variability in the northeast Pacific

The dominant characteristics of wave energy variability in the eastern North Pacific are described from NOAA National Data Buoy Center (NDBC) buoy data collected from 1981 to 2003. Ten buoys at distributed locations were selected for comparison based on record duration and data continuity. Long‐period (LP) [ T > 12] s, intermediate‐period [6 ≤ T ≤ 12] s, and short‐period [ T < 6] s wave spectral energy components are considered separately. Empirical orthogonal function (EOF) analyses of monthly wave energy anomalies reveal that all three wave energy components exhibit similar patterns of spatial variability. The dominant mode represents coherent heightened (or diminished) wave energy along the West Coast from Alaska to southern California, as indicated by composites of the 700 hPa height field. The second EOF mode reveals a distinct El Niño‐Southern Oscillation (ENSO)‐associated spatial distribution of wave energy, which occurs when the North Pacific storm track is extended unusually far south or has receded to the north. Monthly means and principal components (PCs) of wave energy levels indicate that the 1997–1998 El Niño winter had the highest basin‐wide wave energy within this record, substantially higher than the 1982–1983 El Niño. An increasing trend in the dominant PC of LP wave energy suggests that storminess has increased in the northeast Pacific since 1980. This trend is emphasized at central eastern North Pacific locations. Patterns of storminess variability are consistent with increasing activity in the central North Pacific as well as the tendency for more extreme waves in the south during El Niño episodes and in the north during La Niña.

Journal of Geophysical Research C: Oceans↗
USGS2005

A comparison of hydrographically and optically derived mixed layer depths

Efforts to understand and model the dynamics of the upper ocean would be significantly advanced given the ability to rapidly determine mixed layer depths (MLDs) over large regions. Remote sensing technologies are an ideal choice for achieving this goal. This study addresses the feasibility of estimating MLDs from optical properties. These properties are strongly influenced by suspended particle concentrations, which generally reach a maximum at pycnoclines. The premise therefore is to use a gradient in beam attenuation at 660 nm (c660) as a proxy for the depth of a particle-scattering layer. Using a global data set collected during World Ocean Circulation Experiment cruises from 1988-1997, six algorithms were employed to compute MLDs from either density or temperature profiles. Given the absence of published optically based MLD algorithms, two new methods were developed that use c660 profiles to estimate the MLD. Intercomparison of the six hydrographically based algorithms revealed some significant disparities among the resulting MLD values. Comparisons between the hydrographical and optical approaches indicated a first-order agreement between the MLDs based on the depths of gradient maxima for density and c660. When comparing various hydrographically based algorithms, other investigators reported that inherent fluctuations of the mixed layer depth limit the accuracy of its determination to 20 m. Using this benchmark, we found a ???70% agreement between the best hydrographical-optical algorithm pairings. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2005

Numerical modeling of an estuary: A comprehensive skill assessment

Numerical simulations of the Hudson River estuary using a terrain-following, three-dimensional model (Regional Ocean Modeling System (ROMS)) are compared with an extensive set of time series and spatially resolved measurements over a 43 day period with large variations in tidal forcing and river discharge. The model is particularly effective at reproducing the observed temporal variations in both the salinity and current structure, including tidal, spring neap, and river discharge-induced variability. Large observed variations in stratification between neap and spring tides are captured qualitatively and quantitatively by the model. The observed structure and variations of the longitudinal salinity gradient are also well reproduced. The most notable discrepancy between the model and the data is in the vertical salinity structure. While the surface-to-bottom salinity difference is well reproduced, the stratification in the model tends to extend all the way to the water surface, whereas the observations indicate a distinct pycnocline and a surface mixed layer. Because the southern boundary coindition is located near the mouth the estuary, the salinity within the domain is particularly sensitive to the specification of salinity at the boundary. A boundary condition for the horizontal salinity gradient, based on the local value of salinity, is developed to incorporate physical processes beyond the open boundary not resolved by the model. Model results are sensitive to the specification of the bottom roughness length and vertical stability functions, insofar as they influence the intensity of vertical mixing. The results only varied slightly between different turbulence closure methods of k-??, k-??, and k-kl. Copyright 2005 by the American Geophysical Union.

New York↗
USGS2005

Nitrogen and carbon flow from rock to water: Regulation through soil biogeochemical processes, Mokelumne River watershed, California, and Grand Valley, Colorado

Soil denitrification is an ecologically important nitrogen removal mechanism that releases to the atmosphere the greenhouse gas N2O, an intermediate product from the reduction of NO3- to N 2. In this study we evaluate the relationship between soil carbon and denitrification potential in watersheds with bedrock acting as a nonpoint source of nitrogen, testing the hypothesis that nitrate leaching to stream water is in part regulated by denitrification. Two sites, one in a Mediterranean climate and the other in an arid climate, were investigated to understand the interplay between carbon and denitrification potential. Both sites included carbonaceous bedrock with relatively high nitrogen concentrations (> 1,000 mg N kg-1) and had low background nitrogen concentrations in surface and groundwater. There was a net accumulation of carbon and nitrogen in soil relative to the corresponding bedrock, with the exception of carbonaceous shale from the arid site. There the concentration of carbon in the soil (15,620 mg C kg-1) was less than the shale parent (22,460 mg C kg-1), consistent with the bedrock being a source of soil carbon. Rates of denitrification potential (0.5-83 ??g N kg-1 hr-1) derived from laboratory incubations appeared to be related to the ratio of dissolved organic carbon and nitrate extracted from soils. These data indicate that microbial processes such as denitrification can help maintain background nitrogen concentrations to tens of ??M N in relatively undisturbed ecosystems when nitrogen inputs from weathering bedrock are accompanied by sufficient organic carbon concentrations to promote microbial nitrogen transformations.

Journal of Geophysical Research F: Earth Surface↗
USGS2005

Critical shear stress for erosion of cohesive soils subjected to temperatures typical of wildfires

[1] Increased erosion is a well-known response after wildfire. To predict and to model erosion on a landscape scale requires knowledge of the critical shear stress for the initiation of motion of soil particles. As this soil property is temperature-dependent, a quantitative relation between critical shear stress and the temperatures to which the soils have been subjected during a wildfire is required. In this study the critical shear stress was measured in a recirculating flume using samples of forest soil exposed to different temperatures (40??-550??C) for 1 hour. Results were obtained for four replicates of soils derived from three different types of parent material (granitic bedrock, sandstone, and volcanic tuffs). In general, the relation between critical shear stress and temperature can be separated into three different temperature ranges (<175??C; 175??C-275??C; >275??C), which are similar to those for water repellency and temperature. The critical shear stress was most variable (1.0-2.0 N m-2) for temperatures <175??C, was a maximum (>2.0 N m-2) between 175?? and 275??C, and was essentially constant (0.5-0.8 N m-2) for temperatures >275??C. The changes in critical shear stress with temperature were found to be essentially independent of soil type and suggest that erosion processes in burned watersheds can be modeled more simply than erosion processes in unburned watersheds. Wildfire reduces the spatial variability of soil erodibility associated with unburned watersheds by eliminating the complex effects of vegetation in protecting soils and by reducing the range of cohesion associated with different types of unburned soils. Our results indicate that modeling the erosional response after a wildfire depends primarily on determining the spatial distribution of the maximum soil temperatures that were reached during the wildfire. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

Influence of a dam on fine-sediment storage in a canyon river

Glen Canyon Dam has caused a fundamental change in the distribution of fine sediment storage in the 99-km reach of the Colorado River in Marble Canyon, Grand Canyon National Park, Arizona. The two major storage sites for fine sediment (i.e., sand and finer material) in this canyon river are lateral recirculation eddies and the main-channel bed. We use a combination of methods, including direct measurement of sediment storage change, measurements of sediment flux, and comparison of the grain size of sediment found in different storage sites relative to the supply and that in transport, in order to evaluate the change in both the volume and location of sediment storage. The analysis shows that the bed of the main channel was an important storage environment for fine sediment in the predam era. In years of large seasonal accumulation, approximately 50% of the fine sediment supplied to the reach from upstream sources was stored on the main-channel bed. In contrast, sediment budgets constructed for two short-duration, high experimental releases from Glen Canyon Dam indicate that approximately 90% of the sediment discharge from the reach during each release was derived from eddy storage, rather than from sandy deposits on the main-channel bed. These results indicate that the majority of the fine sediment in Marble Canyon is now stored in eddies, even though they occupy a small percentage (???17%) of the total river area. Because of a 95% reduction in the supply of fine sediment to Marble Canyon, future high releases without significant input of tributary sediment will potentially erode sediment from long-term eddy storage, resulting in continued degradation in Marble Canyon. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2005

Nearshore wave-induced cyclical flexing of sea cliffs

[1] Evolution of a tectonically active coast is driven by geomorphically destructive energy supplied by ocean waves. Wave energy is episodic and concentrated; sea cliffs are battered by the geomorphic wrecking ball every 4-25 s. We measure the response of sea cliffs to wave assault by sensing the ground motion using near-coastal seismometers. Sea cliffs respond to waves in two distinct styles. High-frequency motion (20 Hz) reflects the natural frequency of the sea cliff as it rings in response to direct wave impact. Low-frequency motion in the 0.1-0.05 Hz (10-20 s) band consistently agrees with the dominant nearshore wave period. Integrating microseismic velocities suggests 50 ??m and 10 ??m displacements in horizontal and vertical directions, respectively. Displacement ellipsoids exhibit simultaneous downward and seaward sea cliff motion with each wave. Video footage corroborates the downward sea cliff flex in response to the imposed water load on the wave cut platform. Gradients in displacement amplitudes documented using multiple seismometers suggest longitudinal and shear strain of the flexing sea cliff on the order of 0.5-4 ?? strains during each wave loading cycle. As this sea cliff flexure occurs approximately 3 million times annually, it has the potential to fatigue the rock through cyclical loading. Local sea cliff retreat rates of 10 cm/yr imply that a given parcel of rock is flexed through roughly 109 cycles of increasing amplitude before exposure to direct wave attack at the cliff face. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2005

Debris-bed friction of hard-bedded glaciers

[1] Field measurements of debris-bed friction on a smooth rock tablet at the bed of Engabreen, a hard-bedded, temperate glacier in northern Norway, indicated that basal ice containing 10% debris by volume exerted local shear traction of up to 500 kPa. The corresponding bulk friction coefficient between the dirty basal ice and the tablet was between 0.05 and 0.08. A model of friction in which nonrotating spherical rock particles are held in frictional contact with the bed by bed-normal ice flow can account for these measurements if the power law exponent for ice flowing past large clasts is 1. A small exponent (n < 2) is likely because stresses in ice are small and flow is transient. Numerical calculations of the bed-normal drag force on a sphere in contact with a flat bed using n = 1 show that this force can reach values several hundred times that on a sphere isolated from the bed, thus drastically increasing frictional resistance. Various estimates of basal friction are obtained from this model. For example, the shear traction at the bed of a glacier sliding at 20 m a-1 with a geothermally induced melt rate of 0.006 m a-1 and an effective pressure of 300 kPa can exceed 100 kPa. Debris-bed friction can therefore be a major component of sliding resistance, contradicting the common assumption that debris-bed friction is negligible. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2005

Generation and verification of theoretical rating curves in the Whitewater River basin, Kansas

[1] A new method for generating stage-discharge relations (rating curves) for geomorphically stable channels is presented and applied to two streams in the Whitewater River basin, Kansas. The approach converts measurements of stage into discharge using a fluid mechanically based model. The model does not use empirical roughness coefficients, such as Manning coefficients, but rather determines channel roughness from field measurements of the (1) channel geometry, (2) the physical roughness of the bed, banks, and floodplain, and (3) the vegetation density on the banks and floodplain. These measurements are used to calculate explicitly the drag on the small-scale topographic features on the boundary, the drag on the stems and branches of woody vegetation, and the friction on the bed, banks, and floodplain. The theoretical rating curves produced by the model for two study reaches, which are near U.S. Geological Survey (USGS) streamflow-gauging stations, are in good agreement with direct measurements of discharge made by the USGS. Our method has the potential of providing accurate estimates of stream flows less expensively than conventional gauging methods. In addition, the method can be used to obtain more accurate discharge estimates than conventional indirect methods for determining discharge, which are based on estimates of Manning's coefficient.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

Effects of tributary debris on the longitudinal profile of the Colorado River in Grand Canyon

The Colorado River in Grand Canyon has long been known as a "rapids-and-pools" river, with the rapids owing their existence primarily to tributary debris flows. The debris flows deposit subaerial debris fans that constrict the channel laterally and, when they enter the river, raise the bed elevation. The rapids are short-wavelength (???0.1 to ???1 km), small-amplitude (??????5 m) convexities in the river's longitudinal profile, arising from the shallow gradient in the upstream pool and the steep gradient through the rapid itself. Analysis of the entire longitudinal profile through Grand Canyon reveals two long-wavelength (???100 km), large-amplitude (15-30 m) river profile convexities: the eastern canyon convexity between river mile (RM) 30 and RM 80 and the western canyon convexity between RM 150 and RM 250. Convexities of intermediate scale are also identified in the longitudinal profile. These longer-wavelength, larger-amplitude convexities have strong spatial correlations with high rates of debris flow occurrence, high densities of Holocene debris fans, the largest debris fans along the river, and alluvial thicknesses of 10 m or more. River profile convexities are unstable and require an active and powerful geologic process to maintain them, in this case the abundant, frequent, and voluminous Holocene debris flow activity in Grand Canyon. At all wavelengths the most likely cause for these river profile convexities is Holocene aggradation of the riverbed beneath them, driven by the coarse particles of tributary debris flows. Large enough debris flows will slow river flow for kilometers upstream, causing it to drop much of its suspended load. Integrated over time and all of the tributary point source contributions, this process will build short-wavelength convexities into long-wavelength convexities. For most if not all of the Holocene the Colorado River has been dissipating most of its energy in the rapids and expending the remainder in transporting fine sediment through Grand Canyon, with little or no regional incision of bedrock.

Journal of Geophysical Research F: Earth Surface↗
USGS2005

Regulation of landslide motion by dilatancy and pore pressure feedback

A new mathematical model clarifies how diverse styles and rates of landslide motion can result from regulation of Coulomb friction by dilation or contraction of water‐saturated basal shear zones. Normalization of the model equations shows that feedback due to coupling between landslide motion, shear zone volume change, and pore pressure change depends on a single dimensionless parameter α, which, in turn, depends on the dilatancy angle ψ and the intrinsic timescales for pore pressure generation and dissipation. If shear zone soil contracts during slope failure, then α < 0, and positive pore pressure feedback and runaway acceleration are inevitable. If the shear zone dilates, then α > 0, and negative feedback permits slow, steady landslide motion to occur while positive pore pressure is supplied by rain infiltration. Steady state slip velocities v 0 obey v 0 = −( K /ψ) p * e , where K is the hydraulic conductivity and p * e is the normalized (dimensionless) negative pore pressure generated by dilation. If rain infiltration and attendant pore pressure growth continue unabated, however, their influence ultimately overwhelms the stabilizing influence of negative p * e . Then, unbounded landslide acceleration occurs, accentuated by an instability that develops if ψ diminishes as landslide motion proceeds. Nonetheless, numerical solutions of the model equations show that slow, nearly steady motion of a clay‐rich landslide may persist for many months as a result of negative pore pressure feedback that regulates basal Coulomb friction. Similarly stabilized motion is less likely to occur in sand‐rich landslides that are characterized by weaker negative feedback.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

The 1923 Kanto earthquake reevaluated using a newly augmented geodetic data set

This study revisits the mechanism of the 1923 Ms = 7.9 Kanto earthquake in Japan. We derive a new source model and use it to assess quantitative and qualitative aspects of the accommodation of plate motion in the Kanto region. We use a new geodetic data set that consists of displacements from leveling and angle changes from triangulation measurements obtained in surveys between 1883 and 1927. Two unique aspects of our analysis are the inclusion of a large number of second-order triangulation measurements and the application of a correction to remove interseismic deformation. The geometry of the fault planes is adopted from a recent seismic reflection study of the Kanto region. We evaluate the minimum complexity necessary in the model to fit the data optimally. Our final uniform-slip elastic dislocation model consists of two adjacent ???20?? dipping low-angle planes accommodating reverse dextral slip of 6.0 in on the larger, eastern plane and 9.5 m on the smaller, western plane with azimuths of 163?? and 121??, respectively. The earthquake was located in the Sagami trough, where the Philippine Sea plate subducts under Honshu. Compared to the highly oblique angle of plate convergence, the coseismic slip on the large fault plane has a more orthogonal orientation to the strike of the plate boundary, suggesting that slip partitioning plays a role in accommodation of plate motion. What other structure is involved in the partitioning is unclear. Uplift records of marine coastal terraces in Sagami Bay document 7500 years of earthquake activity and predict average recurrence intervals of 400 years for events with vertical displacement profiles similar to those of the 1923 earthquake. This means that the average slip deficit per recurrence interval is ???50% of the relative plate convergence. These findings of plate motion partitioning and slip deficit lead us to suggest that instead of a simple recurrence model with characteristic earthquakes, additional mechanisms are necessary to describe the accommodation of deformation in the Kanto region. So far, obvious candidates for these alternative mechanisms have not been discovered. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Coseismic slip distribution of the 1923 Kanto earthquake, Japan

The slip distribution associated with the 1923 M = 7.9 Kanto, Japan, earthquake is reexamined in light of new data and modeling. We utilize a combination of first-order triangulation, second-order triangulation, and leveling data in order to constrain the coseismic deformation. The second-order triangulation data, which have not been utilized in previous studies of 1923 coseismic deformation, are associated with only slightly smaller errors than the first-order triangulation data and expand the available triangulation data set by about a factor of 10. Interpretation of these data in terms of uniform-slip models in a companion study by Nyst et al. shows that a model involving uniform coseismic slip on two distinct rupture planes explains the data very well and matches or exceeds the fit obtained by previous studies, even one which involved distributed slip. Using the geometry of the Nyst et al. two-plane slip model, we perform inversions of the same geodetic data set for distributed slip. Our preferred model of distributed slip on the Philippine Sea plate interface has a moment magnitude of 7.86. We find slip maxima of ???8-9 m beneath Odawara and ???7-8 m beneath the Miura peninsula, with a roughly 2:1 ratio of strike-slip to dip-slip motion, in agreement with a previous study. However, the Miura slip maximum is imaged as a more broadly extended feature in our study, with the high-slip region continuing from the Miura peninsula to the southern Boso peninsula region. The second-order triangulation data provide good evidence for ???3 m right-lateral strike slip on a 35-km-long splay structure occupying the volume between the upper surface of the descending Philippine Sea plate and the southern Boso peninsula. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗