Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysical Journal”

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 325 records · Page 18Linked to original sources

Time-dependent earthquake probabilities

We have attempted to provide a careful examination of a class of approaches for estimating the conditional probability of failure of a single large earthquake, particularly approaches that account for static stress perturbations to tectonic loading as in the approaches of Stein et al. (1997) and Hardebeck (2004). We have loading as in the framework based on a simple, generalized rate change formulation and applied it to these two approaches to show how they relate to one another. We also have attempted to show the connection between models of seismicity rate changes applied to (1) populations of independent faults as in background and aftershock seismicity and (2) changes in estimates of the conditional probability of failures of different members of a the notion of failure rate corresponds to successive failures of different members of a population of faults. The latter application requires specification of some probability distribution (density function of PDF) that describes some population of potential recurrence times. This PDF may reflect our imperfect knowledge of when past earthquakes have occurred on a fault (epistemic uncertainty), the true natural variability in failure times, or some combination of both. We suggest two end-member conceptual single-fault models that may explain natural variability in recurrence times and suggest how they might be distinguished observationally. When viewed deterministically, these single-fault patch models differ significantly in their physical attributes, and when faults are immature, they differ in their responses to stress perturbations. Estimates of conditional failure probabilities effectively integrate over a range of possible deterministic fault models, usually with ranges that correspond to mature faults. Thus conditional failure probability estimates usually should not differ significantly for these models. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

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↗

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↗

Frictional strength heterogeneity and surface heat flow: Implications for the strength of the creeping San Andreas fault

Heat flow measurements along much of the San Andreas fault (SAF) constrain the apparent coefficient of friction (??app) of the fault to <0.2, much lower than laboratory-derived friction values for most geologic materials. However, heat flow data are sparse near the creeping section of the SAF, a frictional "asperity" where the fault slips almost exclusively by aseismic creep. We test the hypothesis that the creeping section has a substantially higher or lower ?? app than adjacent sections of the SAF. We use numerical models to explore the effects of faults with spatially and temporally heterogeneous frictional strength on the spatial distribution of surface heat flow. Heat flow from finite length asperities is uniformly lower than predicted by assuming an infinitely long fault. Over geologic time, lateral offset from strike-slip faulting produces heat flow patterns that are asymmetric across the fault and along strike. We explore a range of asperity sizes, slip rates, and displacement histories for comparing predicted spatial patterns of heat flow with existing measurements. Models with ??app ??? 0.1 fit the data best. For most scenarios, heat flow anomalies from a frictional asperity with ??app > 0.2 should be detectable even with the sparse existing observations, implying that ??app for the creeping section is as low as the surrounding SAF. Because the creeping section does not slip in large earthquakes, the mechanism controlling its weakness is not related to dynamic processes resulting from high slip rate earthquake ruptures. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Distribution of recent volcanism and the morphology of seamounts and ridges in the GLIMPSE study area: Implications for the lithospheric cracking hypothesis for the origin of intraplate, non-hot spot volcanic chains

Lithospheric cracking by remotely applied stresses or thermoelastic stresses has been suggested to be the mechanism responsible for the formation of intraplate volcanic ridges in the Pacific that clearly do not form above fixed hot spots. As part of the Gravity Lineations Intraplate Melting Petrology and Seismic Expedition (GLIMPSE) project designed to investigate the origin of these features, we have mapped two volcanic chains that are actively forming to the west of the East Pacific Rise using multibeam echo sounding and side-scan sonar. Side-scan sonar reveals the distribution of rough seafloor corresponding to recent, unsedimented lava flows. In the Hotu Matua volcanic complex, recent flows and volcanic edifices are distributed over a region 450 km long and up to 65 km wide, with an apparent, irregular age progression from older flows in the west to younger in the east. The 550-km-long Southern Cross Seamount/Sojourn Ridge/Brown Ridge chain appears to have been recently active only at its eastern end near the East Pacific Rise. A third region of recent flows is found 120 km north of Southern Cross Seamount in seafloor approximately 9 Myr old. No indication of lithospheric extension in the form of faulting or graben formation paralleling the trend of the volcanic chains is found in the vicinity of recent flows or anywhere else in the study area. Thermoelastic cracking could be a factor in the formation of a few small, very narrow volcanic ridges, but most of the volcanic activity is broadly distributed in wide swaths with no indication of formation along narrow cracks. The Sojourn and Brown chains appear to begin as distributed zones of small seamounts that later develop into segmented ridges, perhaps under the influence of membrane stresses from self-loading. We suggest that the linear volcanic chains are created by moving melting anomalies in the asthenosphere and that lithospheric cracking plays at most a secondary role. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Surface drifter derived circulation in the northern and middle Adriatic Sea: Response to wind regime and season

More than 120 satellite-tracked drifters were deployed in the northern and middle Adriatic (NMA) Sea between September 2002 and November 2003, with the purpose of studying the surface circulation at mesoscale to seasonal scale in relation to wind forcing, river runoff, and bottom topography. Pseudo-Eulerian and Lagrangian statistics were calculated from the low-pass-filtered drifter velocity data between September 2002 and December 2003. The structure of the mean circulation is determined with unprecedented high horizontal resolution by the new data. In particular, mean currents, velocity variance, and kinetic energy levels are shown to be maximal in the Western Adriatic Current (WAC). Separating data into seasons, we found that the mean kinetic energy is maximal in fall, with high values also in winter, while it is significantly weaker in summer. High-resolution Local Area Model Italy winds were used to relate the drifter velocities to the wind fields. The surface currents appear to be significantly influenced by the winds. The mean flow during the northeasterly bora regime shows an intensification of the across-basin recirculating currents. In addition, the WAC is strongly intensified both in intensity and in its offshore lateral extension. In the southeasterly sirocco regime, northward flow without recirculation dominates in the eastern half of the basin, while during northwesterly maestro the WAC is enhanced. Separating the data into low and high Po River discharge rates for low-wind conditions shows that the WAC and the velocity fluctuations in front of the Po delta are stronger for high Po River runoff. Lagrangian covariance, diffusivity, and integral time and space scales are larger in the along-basin direction and are maximal in the southern portion of the WAC. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗

Coseismic source model of the 2003 Mw 6.8 Chengkung earthquake, Taiwan, determined from GPS measurements

A coseismic source model of the 2003 Mw 6.8 Chengkung, Taiwan, earthquake was well determined with 213 GPS stations, providing a unique opportunity to study the characteristics of coseismic displacements of a high-angle buried reverse fault. Horizontal coseismic displacements show fault-normal shortening across the fault trace. Displacements on the hanging wall reveal fault-parallel and fault-normal lengthening. The largest horizontal and vertical GPS displacements reached 153 and 302 mm, respectively, in the middle part of the network. Fault geometry and slip distribution were determined by inverting GPS data using a three-dimensional (3-D) layered-elastic dislocation model. The slip is mainly concentrated within a 44 ?? 14 km slip patch centered at 15 km depth with peak amplitude of 126.6 cm. Results from 3-D forward-elastic model tests indicate that the dome-shaped folding on the hanging wall is reproduced with fault dips greater than 40??. Compared with the rupture area and average slip from slow slip earthquakes and a compilation of finite source models of 18 earthquakes, the Chengkung earthquake generated a larger rupture area and a lower stress drop, suggesting lower than average friction. Hence the Chengkung earthquake seems to be a transitional example between regular and slow slip earthquakes. The coseismic source model of this event indicates that the Chihshang fault is divided into a creeping segment in the north and the locked segment in the south. An average recurrence interval of 50 years for a magnitude 6.8 earthquake was estimated for the southern fault segment. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Northerly surface winds over the eastern North Pacific Ocean in spring and summer

Persistent spring and summer northerly surface winds are the defining climatological feature of the western coast of North America, especially south of the Oregon coast. Northerly surface winds are important for upwelling and a vast array of other biological, oceanic, and atmospheric processes. Intermittence in northerly coastal surface wind is characterized and wind events are quantitatively defined using coastal buoy data south of Cape Mendocino on the northern California coast. The defined wind events are then used as a basis for composites in order to explain the spatial evolution of various atmospheric and oceanic processes. Wind events involve large-scale changes in the three-dimensional atmospheric circulation including the eastern North Pacific subtropical anticyclone and southeast trade winds. Composites of QSCAT satellite scatterometer wind estimates from 1999 to 2005 based on a single coastal buoy indicate that wind events typically last 72-96 h and result in anomalies in surface wind and Ekman pumping that extend over 1000 kin from the west coast of North America. It may be useful to consider ocean circulation and dependent ecosystem dynamics and the distribution of temperature, moisture, and aerosols in the atmospheric boundary layer in the context of wind events defined herein. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

In situ Raman spectroscopic investigation of the structure of subduction-zone fluids

In situ Raman spectra of synthetic subduction-zone fluids (KAlSi3O8-H2O system) were measured to 900?? and 2.3 GPa using a hydrothermal diamond-anvil cell. The structures of aqueous fluid and hydrous melt become closer when conditions approach the second critical endpoint. Almost no three-dimensional network was observed in the supercritical fluid above 2 GPa although a large amount of silicate component is dissolved, suggesting that the physical and chemical properties of these phases change drastically at around the second critical endpoint. Our experimental results indicate that the fluids released from a subducting slab change from aqueous fluid to supercritical fluid with increasing depth under the volcanic arcs. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Constraining earthquake source inversions with GPS data: 1. Resolution-based removal of artifacts

We present a resolution analysis of an inversion of GPS data from the 2004 M w 6.0 Parkfield earthquake. This earthquake was recorded at thirteen 1-Hz GPS receivers, which provides for a truly coseismic data set that can be used to infer the static slip field. We find that the resolution of our inverted slip model is poor at depth and near the edges of the modeled fault plane that are far from GPS receivers. The spatial heterogeneity of the model resolution in the static field inversion leads to artifacts in poorly resolved areas of the fault plane. These artifacts look qualitatively similar to asperities commonly seen in the final slip models of earthquake source inversions, but in this inversion they are caused by a surplus of free parameters. The location of the artifacts depends on the station geometry and the assumed velocity structure. We demonstrate that a nonuniform gridding of model parameters on the fault can remove these artifacts from the inversion. We generate a nonuniform grid with a grid spacing that matches the local resolution length on the fault and show that it outperforms uniform grids, which either generate spurious structure in poorly resolved regions or lose recoverable information in well-resolved areas of the fault. In a synthetic test, the nonuniform grid correctly averages slip in poorly resolved areas of the fault while recovering small-scale structure near the surface. Finally, we present an inversion of the Parkfield GPS data set on the nonuniform grid and analyze the errors in the final model. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Simulation of Na D emission near Europa during eclipse

The Cassini imaging science subsystem observed Europa in eclipse during Cassini's Jupiter flyby. The disk-resolved observations revealed a spatially nonuniform emission in the wavelength range of 200-1050 nm (clear filters). By building on observations and simulations of Europa's Na atmosphere and torus we find that electron-excited Na in Europa's tenuous atmosphere can account for the observed emission if the Na is ejected preferentially from Europa's dark terrain. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Magnetic fabric of sheared till: A strain indicator for evaluating the bed deformation model of glacier flow

Wet-based portions of ice sheets may move primarily by shearing their till beds, resting in high sediment fluxes and the development of subglacial landforms. This model of glacier movement, which requires high bed shear strains, can be tested using till microstructural characteristics that evolve during till deformation. Here we examine the development of magnetic fabric using a ring shear device to defom two Wisconsin-age basal tills to shear strains as high as 70. Hysteresis experiments and the dependence of magnetic susceptibility of these tills on temperature demonstrate that anisotropy of magnetic susceptibility (AMS) develops during shear due to the rotation of primarily magnetite particles that are silt sized or smaller. At moderate shear strains (???6-25), principal axes of maximum magnetic susceptibility develop a strong fabric (S1 eignevalues of 0.83-0.96), without further strengthening at higher strains, During deformation, directions of maximum susceptibility cluster strongly in the direction of shear and plunge 'up-glacier,' consistent with the behavior of pebbles and sand particles studied in earlier experiments. In contrast, the magnitude of AMS does not vary systematically with strain and is small relative to its variability among samples; this is because most magnetite grains are contained as inclusions in larger particles and hence do not align during shear. Although processes other than pervasive bed deformation may result in strong flow parallel fabrics, AMS fabrics provide a rapid and objective means of identifying basal tills that have not been sheared sufficiently to be compatible with the bed deformation model. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

Regional three-dimensional seismic velocity model of the crust and uppermost mantle of northern California

We present a three-dimensional (3D) tomographic model of the P wave velocity (Vp) structure of northern California. We employed a regional-scale double-difference tomography algorithm that incorporates a finite-difference travel time calculator and spatial smoothing constraints. Arrival times from earthquakes and travel times from controlled-source explosions, recorded at network and/or temporary stations, were inverted for Vp on a 3D grid with horizontal node spacing of 10 to 20 km and vertical node spacing of 3 to 8 km. Our model provides an unprecedented, comprehensive view of the regional-scale structure of northern California, putting many previously identified features into a broader regional context and improving the resolution of a number of them and revealing a number of new features, especially in the middle and lower crust, that have never before been reported. Examples of the former include the complex subducting Gorda slab, a steep, deeply penetrating fault beneath the Sacramento River Delta, crustal low-velocity zones beneath Geysers-Clear Lake and Long Valley, and the high-velocity ophiolite body underlying the Great Valley. Examples of the latter include mid-crustal low-velocity zones beneath Mount Shasta and north of Lake Tahoe. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Rupture parameters of the 2003 Zemmouri (M w 6.8), Algeria, earthquake from joint inversion of interferometric synthetic aperture radar, coastal uplift, and GPS

We study the surface deformation associated with the 21 May 2003 (M w = 6.8) Zemmouri (Algeria) earthquake, the strongest seismic event felt in the Algiers region since 1716. The thrust earthquake mechanism and related surface deformation revealed an average 0.50 m coastal uplift along ??55-km-long coastline. We obtain coseismic interferograms using Envisat advanced synthetic aperture radar (ASAR) (IS2) and RADARSAT standard beam (ST4) data from both the ascending and descending orbits of Envisat satellite, whereas the RADARSAT data proved useful only in the descending mode. While the two RADARSAT interferograms cover the earthquake area, Envisat data cover only the western half of the rupture zone. Although the interferometric synthetic aperture radar (InSAR) coherence in the epicenter area is poor, deformation fringes are observed along the coast in different patches. In the Boumerdes area, the maximum coseismic deformation is indicated by the high gradient of fringes visible in all interferograms in agreement with field measurements (tape, differential GPS, leveling, and GPS). To constrain the earthquake rupture parameters, we model the interferograms and uplift measurements using elastic dislocations on triangular fault patches in an elastic and homogeneous half-space. We invert the coseismic slip using first, a planar surface and second, a curved fault, both constructed from triangular elements using Poly3Dinv program that uses a damped least square minimization. The best fit of InSAR, coastal uplift, and GPS data corresponds to a 65-km-long fault rupture dipping 40?? to 50?? SE, located at 8 to 13 km offshore with a change in strike west of Boumerdes from N60??-65?? to N95??-105??. The inferred rupture geometry at depth correlates well with the seismological results and may have critical implications for the seismic hazard assessment of the Algiers region. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

Shoaling of nonlinear internal waves in Massachusetts Bay

The shoaling of the nonlinear internal tide in Massachusetts Bay is studied with a fully nonlinear and nonhydrostatic model. The results are compared with current and temperature observations obtained during the August 1998 Massachusetts Bay Internal Wave Experiment and observations from a shorter experiment which took place in September 2001. The model shows how the approaching nonlinear undular bore interacts strongly with a shoaling bottom, offshore of where KdV theory predicts polarity switching should occur. It is shown that the shoaling process is dominated by nonlinearity, and the model results are interpreted with the aid of a two-layer nonlinear but hydrostatic model. After interacting with the shoaling bottom, the undular bore emerges on the shallow shelf inshore of the 30-m isobath as a nonlinear internal tide with a range of possible shapes, all of which are found in the available observational record. Copyright 2008 by the American Geophysical Union.

Masschusetts↗

Moon meteoritic seismic hum: Steady state prediction

We use three different statistical models describing the frequency of meteoroid impacts on Earth to estimate the seismic background noise due to impacts on the lunar surface. Because of diffraction, seismic events on the Moon are typically characterized by long codas, lasting 1 h or more. We find that the small but frequent impacts generate seismic signals whose codas overlap in time, resulting in a permanent seismic noise that we term the "lunar hum" by analogy with the Earth's continuous seismic background seismic hum. We find that the Apollo era impact detection rates and amplitudes are well explained by a model that parameterizes (1) the net seismic impulse due to the impactor and resulting ejecta and (2) the effects of diffraction and attenuation. The formulation permits the calculation of a composite waveform at any point on the Moon due to simulated impacts at any epicentral distance. The root-mean-square amplitude of this waveform yields a background noise level that is about 100 times lower than the resolution of the Apollo long-period seismometers. At 2 s periods, this noise level is more than 1000 times lower than the low noise model prediction for Earth's microseismic noise. Sufficiently sensitive seismometers will allow the future detection of several impacts per day at body wave frequencies. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

Cretaceous-Eocene (Laramide) landscape development and Oligocene- Pliocene drainage reorganization of transition zone and Colorado Plateau, Arizona

Landscape development of central and northern Arizona can no longer be ascribed mainly to events of Miocene and Pliocene age. New information on the age and distribution of older Cenozoic deposits has led to the recognition of a regional Cretaceous-Paleocene(?) surface of erosion that conforms to major elements of the present topography and to the recognition that a formerly thick deposit of gravel accumulated on this regional surface of erosion. These relations cast new light on the history of evolution of the landscape and indicate a much greater age for the main landscape elements and a more complicated and prolonged history of erosion and deposition than has been previously supposed. The timing of events postulated for development of drainage on the Colorado Plateau can now be compared and partly reconciled with events recognized in the adjacent closely related Mountain Region (Transition Zone) of central Arizona. As a consequence of Late Cretaceous-Paleocene (Laramide) compression, central and northern Arizona underwent at least 1200 m of uplift, documented by paleochannels cut into erosionally truncated Paleozoic strata on the Hualapai Plateau of the southwestern Colorado Plateau. During this time, a highly irregular erosion surface was developed on Proterozoic rocks across the Transition Zone south of the Mogollon Rim, the scarp of the Mogollon Rim was eroded to its present height (600–900 m), and an extensive stripped surface was developed on resistant upper Paleozoic strata north of the rim. Deposition of several hundred meters of Paleocene-Eocene “Rim gravels” derived from highlands south and west of the region followed, covering much of the Cretaceous-Paleocene erosion surface. Nearly complete burial of the rim is suggested by the distribution of remnants of the Rim gravels across the erosional scarps and on high plateau areas north of the rim. A second increment of uplift, apparently occurring in late Eocene time and apparently recorded by a series of fission track cooling ages from the Marble and Grand canyons, is inferred to have been responsible for ending deposition of the Rim gravels, for initiating differential uplift of contemporaneous deposits (Canaan Peak and Claron formations) to their positions in the high plateaus of central Utah, and for causing the drainage reorganization required to explain the extensive removal of Rim gravels from much of the region. A southerly flowing ancestral Verde River related to the drainage reorganization removed much of the older gravel cover from the Transition Zone of central Arizona, resulting in a younger regional erosion surface having 600–900 m of relief, a surface closely approximating the Cretaceous-Paleocene erosion surface. Late Oligocene and early Miocene rocks locally rest unconformably on remnants of Rim gravels in the Transition Zone, indicating that the second episode of regional erosion had been completed by late Oligocene time. North of the Mogollon Rim, a west flowing(?) ancestral Colorado River is inferred to have become established on the Rim gravels, draining the interior parts of the Colorado Plateau and transporting detritus off the plateau. Exhumation of the Mogollon Rim and development of 600–900 m of topographic relief in the Transition Zone by an ancestral Verde River system suggests the potential for a comparable, coeval entrenchment of an ancestral Colorado River in Paleozoic strata north of the Mogollon Rim. Regional extension and volcanic activity ensued in late Oligocene to Pliocene time. The Oligocene erosion surface in the extensional basins of central Arizona became largely concealed by Miocene and Pliocene deposits as the Neogene climate became drier. In late Miocene and Pliocene time, perennial streams appear to have been lacking, transport of detritus appears to have been principally by flash flooding, little or no detritus appears to have been removed from the region, and much of the precipitation presumably moved by groundwater flow through the deposits of aggradation. A coeval episode of aggradation in the Grand Canyon is suggested by deposits that appear to have once choked much of the canyon. If this event parallels the episode of late Miocene and Pliocene aggradation recorded east, south, and west of the Grand Canyon, the Colorado River could have been incised to its present level by late Miocene time. A return to wetter conditions in late Pliocene time presumably was responsible for renewed erosion and reexcavation of older drainages and basins. An understanding of this Tertiary structural, erosional, and depositional history can be important for the geological analysis of geophysical transects across the region.

Journal of Geophysical Research↗

Nonlinear refraction and reflection travel time tomography

We develop a rapid nonlinear travel time tomography method that simultaneously inverts refraction and reflection travel times on a regular velocity grid. For travel time and ray path calculations, we apply a wave front method employing graph theory. The first-arrival refraction travel times are calculated on the basis of cell velocities, and the later refraction and reflection travel times are computed using both cell velocities and given interfaces. We solve a regularized nonlinear inverse problem. A Laplacian operator is applied to regularize the model parameters (cell slownesses and reflector geometry) so that the inverse problem is valid for a continuum. The travel times are also regularized such that we invert travel time curves rather than travel time points. A conjugate gradient method is applied to minimize the nonlinear objective function. After obtaining a solution, we perform nonlinear Monte Carlo inversions for uncertainty analysis and compute the posterior model covariance. In numerical experiments, we demonstrate that combining the first arrival refraction travel times with later reflection travel times can better reconstruct the velocity field as well as the reflector geometry. This combination is particularly important for modeling crustal structures where large velocity variations occur in the upper crust. We apply this approach to model the crustal structure of the California Borderland using ocean bottom seismometer and land data collected during the Los Angeles Region Seismic Experiment along two marine survey lines. Details of our image include a high-velocity zone under the Catalina Ridge, but a smooth gradient zone between. Catalina Ridge and San Clemente Ridge. The Moho depth is about 22 km with lateral variations. Copyright 1998 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗