SEARCH · USGS Science

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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15Linked to original sources

USGS2016

Volcano dome dynamics at Mount St. Helens: Deformation and intermittent subsidence monitored by seismicity and camera imagery pixel offsets

The surface deformation field measured at volcanic domes provides insights into the effects of magmatic processes, gravity- and gas-driven processes, and the development and distribution of internal dome structures. Here we study short-term dome deformation associated with earthquakes at Mount St. Helens, recorded by a permanent optical camera and seismic monitoring network. We use Digital Image Correlation (DIC) to compute the displacement field between successive images and compare the results to the occurrence and characteristics of seismic events during a 6 week period of dome growth in 2006. The results reveal that dome growth at Mount St. Helens was repeatedly interrupted by short-term meter-scale downward displacements at the dome surface, which were associated in time with low-frequency, large-magnitude seismic events followed by a tremor-like signal. The tremor was only recorded by the seismic stations closest to the dome. We find a correlation between the magnitudes of the camera-derived displacements and the spectral amplitudes of the associated tremor. We use the DIC results from two cameras and a high-resolution topographic model to derive full 3-D displacement maps, which reveals internal dome structures and the effect of the seismic activity on daily surface velocities. We postulate that the tremor is recording the gravity-driven response of the upper dome due to mechanical collapse or depressurization and fault-controlled slumping. Our results highlight the different scales and structural expressions during growth and disintegration of lava domes and the relationships between seismic and deformation signals.

Washington↗

USGS2019

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

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

Wyoming↗

USGS2020

Shallow basin structure and attenuation are key to predicting long shaking duration in Los Angeles Basin

Ground motions in the Los Angeles Basin during large earthquakes are modulated by earthquake ruptures, path effects into the basin, basin effects, and local site response. We analyzed the direct effect of shallow basin structures on shaking duration at a period of 2–10 s in the Los Angeles region through modeling small magnitude, shallow, and deep earthquake pairs. The source depth modulates the basin response, particularly the shaking duration, and these features are a function of path effect and not site condition. Three-dimensional simulations using the CVM-S4.26.M01 velocity model show good fitting to the initial portion of the waveforms at periods of 5 s and longer but fail to predict the long shaking duration during shallow events, especially at periods less than 5 s. Simulations using CVM-H do not match the timing of the initial arrivals as well as CVM-S4.26.M01, and the strong late arrivals in the CVM-H simulation travel with an apparent velocity slower than observed. A higher-quality factor than traditionally assumed may produce synthetics with longer durations but is unable to accurately match the amplitude and phase. Beamforming analysis using dense array data further reveals the long duration surface waves have the same back azimuth as the direct arrivals and are generated at the basin edges, while the later coda waves are scattered from off-azimuth directions, potentially due to strong, sharp boundaries offshore. Improving the description of these shallow basin structures and attenuation model will enhance our capability to predict long-period ground motions in basins.

California↗

USGS2006

Analysis of gravity and topography in the GLIMPSE study region: Isostatic compensation and uplift of the Sojourn and Hotu Matua Ridge systems

The Gravity Lieations Intraplate Melting Petrologic and Seismic Expedition (GLIMPSE) Experiment investigated the formation of a series of non-hot spot, intraplate volcanic ridges in the South Pacific and their relationship to cross-grain gravity lineaments detected by satellite altimetry. Using shipboard gravity measurements and a simple model of surface loading of a thin elastic plate, we estimate effective elastic thicknesses ranging from ???2 km beneath the Sojourn Ridge to a maximum of 10 km beneath the Southern Cross Seamount. These elastic thicknesses are lower than predicted for the 3-9 Ma seafloor on which the volcanoes lie, perhaps due to reheating and thinning of the plate during emplacement. Anomalously low apparent densities estimated for the Matua and Southern Cross seamounts 2050 and 2250 kg m-3, respectively, probably are artifacts caused by the assumption of only surface loading, ignoring the presence of subsurface loading in the form of underplated crust and/or low-density mantle. Using satellite free-air gravity and shipboard bathymetry, we calculate the age-detrended, residual mantle Bouguer anomaly (rMBA). The rMBA corrects the free-air anomaly for the direct effects of topography, including the thickening of the crust beneath the seamounts and volcanic ridges due to surface loading of the volcanic edifices. There are broad, negative rMBA anomalies along the Sojourn and Brown ridges and the Hotu Matua seamount chain that extend nearly to the East Pacific Rise. These negative rMBA anomalies connect to negative free-air anomalies in the western part of the study area that have been recognized previously as the beginnings of the cross-grain gravity lineaments. Subtracting the topographic effects of surface loading by the ridges and seamounts from the observed topography reveals that the ridges are built on broad bands of anomalously elevated seafloor. This swell topography and the negative rMBA anomalies contradict the predictions of lithospheric cracking models for the origin of gravity lineaments and associated volcanic ridges, favoring models with a dynamic mantle component such as small-scale convection or channelized asthenospheric return flow. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2019

Delayed dynamic triggering of disposal-induced earthquakes observed by a dense array in Northern Oklahoma

Recent increases in earthquake occurrence rates in Oklahoma have been linked to the injection of large volumes of saltwater, a byproduct of oil and gas extraction. Here we present a detailed study of remote earthquake triggering in an area of active injection‐induced seismicity in northern Oklahoma using data from the LArge‐n Seismic Survey in Oklahoma (LASSO) temporary array and nearby permanent broadband seismic stations. We estimate changes in earthquake rates and calculate the Coulomb failure stress changes on potential receiver faults due to passing teleseismic surface waves. A statistically significant increase in seismicity is observed ∼8 hr after the 16 April 2016 M w 7.8 Ecuador earthquake. The Coulomb stress changes associated with the Ecuador earthquake are on the order of ∼1 kPa. Physical mechanisms consistent with the observed dynamic stress threshold include failure driven by activation of aseismic slip or hydrological response of the fault system.

Oklahoma↗

USGS1997

Multiphase groundwater flow near cooling plutons

We investigate groundwater flow near cooling plutons with a computer program that can model multiphase flow, temperatures up to 1200°C, thermal pressurization, and temperature‐dependent rock properties. A series of experiments examines the effects of host‐rock permeability, size and depth of pluton emplacement, single versus multiple intrusions, the influence of a caprock, and the impact of topographically driven groundwater flow. We also reproduce and evaluate some of the pioneering numerical experiments on flow around plutons. Host‐rock permeability is the principal factor influencing fluid circulation and heat transfer in hydrothermal systems. The hottest and most steam‐rich systems develop where permeability is of the order of 10 −15 m 2 . Temperatures and life spans of systems decrease with increasing permeability. Conduction‐dominated systems, in which permeabilities are ≤10 −16 m 2 , persist longer but exhibit relatively modest increases in near‐surface temperatures relative to ambient conditions. Pluton size, emplacement depth, and initial thermal conditions have less influence on hydrothermal circulation patterns but affect the extent of boiling and duration of hydrothermal systems. Topographically driven groundwater flow can significantly alter hydrothermal circulation; however, a low‐permeability caprock effectively decouples the topographically and density‐driven systems and stabilizes the mixing interface between them thereby defining a likely ore‐forming environment.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

The birth and death of transverse aeolian ridges on Mars

Transverse aeolian ridges (TARs) are small bright windblown deposits found throughout the Martian tropics that stand a few meters tall and are spaced a few tens of meters apart. The origin of these features remains mysterious more than 20 years after their discovery on Mars. This paper presents a new hypothesis, that some of the TARs could be indurated dust deposits emplaced millions of years ago during periods of higher axial obliquity. It suggests that these TARs are primary depositional bed forms that accumulated in place from dust carried by the winds in suspension, perhaps in a manner comparable to antidunes on Earth, and were subsequently indurated and eroded to their current states by eons of sandblasting. It points out examples of modern dust drifts and dune‐like features that appear to have been recently formed by dust accumulating directly onto the surface from atmospheric suspension. It shows how these pristine dust deposits could evolve to explain the range of morphologies of the TARs. Finally, it explains how the known properties of many TARs are consistent with this hypothesis, including their composition, thermal behavior, and distribution.

Journal of Geophysical Research E: Planets↗

USGS2016

Fault segmentation: New concepts from the Wasatch Fault Zone, Utah, USA

The question of whether structural segment boundaries along multisegment normal faults such as the Wasatch fault zone (WFZ) act as persistent barriers to rupture is critical to seismic hazard analyses. We synthesized late Holocene paleoseismic data from 20 trench sites along the central WFZ to evaluate earthquake rupture length and fault segmentation. For the youngest (<3 ka) and best-constrained earthquakes, differences in earthquake timing across prominent primary segment boundaries, especially for the most recent earthquakes on the north-central WFZ, are consistent with segment-controlled ruptures. However, broadly constrained earthquake times, dissimilar event times along the segments, the presence of smaller-scale (subsegment) boundaries, and areas of complex faulting permit partial-segment and multisegment (e.g., spillover) ruptures that are shorter (~20–40 km) or longer (~60–100 km) than the primary segment lengths (35–59 km). We report a segmented WFZ model that includes 24 earthquakes since ~7 ka and yields mean estimates of recurrence (1.1–1.3 kyr) and vertical slip rate (1.3–2.0 mm/yr) for the segments. However, additional rupture scenarios that include segment boundary spatial uncertainties, floating earthquakes, and multisegment ruptures are necessary to fully address epistemic uncertainties in rupture length. We compare the central WFZ to paleoseismic and historical surface ruptures in the Basin and Range Province and central Italian Apennines and conclude that displacement profiles have limited value for assessing the persistence of segment boundaries but can aid in interpreting prehistoric spillover ruptures. Our comparison also suggests that the probabilities of shorter and longer ruptures on the WFZ need to be investigated.

Utah↗

USGS2003

Late Pleistocene to Holocene slip rates for the Gurvan Bulag thrust fault (Gobi-Altay, Mongolia) estimated with 10Be dates

We surveyed morphotectonic markers along the central part of the Gurvan Bulag thrust, a fault that ruptured with the Bogd fault during the Gobi-Altay earthquake (1957, M 8.3), to document climatic and tectonic processes along the fault for the late Pleistocene- Holocene period. The markers were dated using 10Be produced in situ. Two major periods of alluviation ended at 131 ?? 20 and 16 ?? 4.8 ka. These appear to be contemporaneous with global climatic changes at the terminations of marine isotope stages (MIS) 6 and 2. The vertical slip rates, determined from offset measurements and surfaces ages, are 0.14 ?? 0.03 mm/yr over the late Pleistocene-Holocene and between 0.44 ?? 0.11 and 1.05 ?? 0.25 mm/yr since the end of the late Pleistocene. The higher of these slip rates for the last ???16 kyr is consistent with paleoseismic investigations along the fault [Prentice et al., 2002], and suggests that, at the end of late Pleistocene, the fault evolved from quiescence to having recurrence intervals of 4.0 ?? 1.2 kyr for surface ruptures with ???4 m vertical offset (similar to that of 1957). The inferred recurrence interval is comparable to that of the Bogd fault (3.7 ?? 1.3 kyr) suggesting that the two faults may have ruptured together also earlier during the last ???16 kyr.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

A magmatic model of Medicine Lake Volcano, California

Medicine Lake volcano is a Pleistocene and Holocene shield volcano of the southern Cascade Range. It is located behind the main Cascade arc in an extensional tectonic setting where high-alumina basalt is the most commonly erupted lava. This basalt is parental to the higher-silica calc-alkaline and tholeiitic lavas that make up the bulk of the shield. The presence of late Holocene, chemically identical rhyolites on opposite sides of the volcano led to hypotheses of a large shallow silicic magma chamber and of a small, deep chamber that fed rhyolites to the surface via cone sheets. Subsequent geophysical work has been unable to identify a large silicic magma body, and instead a small one has apparently been recognized. Some geologic data support the geophysical results. Tectonic control of vent alignments and the dominance of mafic eruptions both in number of events and volume throughout the history of the volcano indicate that no large silicic magma reservoir exists. Instead, a model is proposed that includes numerous dikes, sills, and small magma bodies, most of which are too small to be recognized by present geophysical methods.

Journal of Geophysical Research Solid Earth↗

USGS2012

Tidal triggering of low frequency earthquakes near Parkfield, California: Implications for fault mechanics within the brittle-ductile transition

Studies of nonvolcanic tremor (NVT) have established the significant impact of small stress perturbations on NVT generation. Here we analyze the influence of the solid earth and ocean tides on a catalog of ∼550,000 low frequency earthquakes (LFEs) distributed along a 150 km section of the San Andreas Fault centered at Parkfield. LFE families are identified in the NVT data on the basis of waveform similarity and are thought to represent small, effectively co-located earthquakes occurring on brittle asperities on an otherwise aseismic fault at depths of 16 to 30 km. We calculate the sensitivity of each of these 88 LFE families to the tidally induced right-lateral shear stress (RLSS), fault-normal stress (FNS), and their time derivatives and use the hypocentral locations of each family to map the spatial variability of this sensitivity. LFE occurrence is most strongly modulated by fluctuations in shear stress, with the majority of families demonstrating a correlation with RLSS at the 99% confidence level or above. Producing the observed LFE rate modulation in response to shear stress perturbations requires low effective stress in the LFE source region. There are substantial lateral and vertical variations in tidal shear stress sensitivity, which we interpret to reflect spatial variation in source region properties, such as friction and pore fluid pressure. Additionally, we find that highly episodic, shallow LFE families are generally less correlated with tidal stresses than their deeper, continuously active counterparts. The majority of families have weaker or insignificant correlation with positive (tensile) FNS. Two groups of families demonstrate a stronger correlation with fault-normal tension to the north and with compression to the south of Parkfield. The families that correlate with fault-normal clamping coincide with a releasing right bend in the surface fault trace and the LFE locations, suggesting that the San Andreas remains localized and contiguous down to near the base of the crust. The deep families that have high sensitivity to both shear and tensile normal stress perturbations may be indicative of an increase in effective fault contact area with depth. Synthesizing our observations with those of other LFE-hosting localities will help to develop a comprehensive understanding of transient fault slip below the “seismogenic zone” by providing constraints on parameters in physical models of slow slip and LFEs.

California↗

USGS2003

Evidence for a Battle Mountain-Eureka crustal fault zone, north-central Nevada, and its relation to Neoproterozoic-Early Paleozoic continental breakup

Combined evidence from gravity, radiogenic isotope, and magnetotelluric (MT) data indicates a crustal fault zone that coincides with the northwest-trending Battle Mountain-Eureka (BME) mineral trend in north-central Nevada, USA. The BME crustal fault zone likely originated during Neoproterozoic-Early Paleozoic rifting of the continent and had a large influence on subsequent tectonic events, such as emplacement of allochthons and episodic deformation, magmatism, and mineralization throughout the Phanerozoic. MT models show the fault zone is about 10 km wide, 130-km long, and extends from 1 to 5 km below the surface to deep crustal levels. Isotope data and gravity models imply the fault zone separates crust of fundamentally different character. Geophysical evidence for such a long-lived structure, likely inherited from continental breakup, defies conventional wisdom that structures this old have been destroyed by Cenozoic extensional processes. Moreover, the coincidence with the alignment of mineral deposits supports the assertion by many economic geologists that these alignments are indicators of buried regional structures.

Journal of Geophysical Research B: Solid Earth↗

USGS2016

Tomographic image of a seismically active volcano: Mammoth Mountain, California

High-resolution tomographic P wave, S wave, and V P / V S velocity structure models are derived for Mammoth Mountain, California, using phase data from the Northern California Seismic Network and a temporary deployment of broadband seismometers. An anomalous volume (5.1 × 10 9 to 5.9 × 10 10 m 3 ) of low P and low S wave velocities is imaged beneath Mammoth Mountain, extending from near the surface to a depth of &sim;2 km below sea level. We infer that the reduction in seismic wave velocities is due to the presence of CO 2 distributed in oblate spheroid pores with mean aspect ratio α = 1.6 × 10 −3 to 7.9 × 10 −3 (crack-like pores) and mean gas volume fraction ϕ = 8.1 × 10 −4 to 3.4 × 10 −3 . The pore density parameter κ = 3 ϕ /(4π α ) = na 3 =0.11, where n is the number of pores per cubic meter and a is the mean pore equatorial radius. The total mass of CO 2 is estimated to be 4.6 × 10 9 to 1.9 × 10 11 kg. The local geological structure indicates that the CO 2 contained in the pores is delivered to the surface through fractures controlled by faults and remnant foliation of the bedrock beneath Mammoth Mountain. The total volume of CO 2 contained in the reservoir suggests that given an emission rate of 500 tons day −1 , the reservoir could supply the emission of CO 2 for &sim;25–1040 years before depletion. Continued supply of CO 2 from an underlying magmatic system would significantly prolong the existence of the reservoir.

California↗

USGS2023

Identification of seismo-volcanic regimes at Whakaari/White Island (New Zealand) via systematic tuning of an unsupervised classifier

We present an algorithm based on Self-Organizing Maps (SOM) and k-means clustering to recognize patterns in a continuous 12.5-year tremor time series recorded at Whakaari/White Island volcano, New Zealand (hereafter referred to as Whakaari). The approach is extendable to a variety of volcanic settings through systematic tuning of the classifier. Hyperparameters are evaluated by statistical means, yielding a combination of “ideal” SOM parameters for the given data set. Extending from this, we applied a Kernel Density Estimation approach to automatically detect changes within the observed seismicity. We categorize the Whakaari seismic time series into regimes representing distinct volcano-seismic states during recent unrest episodes at Whakaari (2012/2013, 2016, and 2019). There is a clear separation in classification results between background regimes and those representing elevated levels of unrest. Onset of unrest is detected by the classifier 6 weeks before the August 2012 eruption, and ca. 3.5 months before the December 2019 eruption, respectively. Regime changes are corroborated by changes in commonly monitored tremor proxies as well as with reported volcanic activity. The regimes are hypothesized to represent diverse mechanisms including: system pressurization and depressurization, degassing, and elevated surface activity. Labeling these regimes improves visualization of the 2012/2013 and 2019 unrest and eruptive episodes. The pre-eruptive 2016 unrest showed a contrasting shape and nature of seismic regimes, suggesting differing onset and driving processes. The 2016 episode is proposed to result from rapid destabilization of the shallow hydrothermal system, while rising magmatic gases from new injections of magma better explain the 2012/2013 and 2019 episodes.

Whakaari, White Island↗

USGS2009

On near-source earthquake triggering

[1] When one earthquake triggers others nearby, what connects them? Two processes are observed: static stress change from fault offset and dynamic stress changes from passing seismic waves. In the near-source region ( r ≤ 50 km for M ∼ 5 sources) both processes may be operating, and since both mechanisms are expected to raise earthquake rates, it is difficult to isolate them. We thus compare explosions with earthquakes because only earthquakes cause significant static stress changes. We find that large explosions at the Nevada Test Site do not trigger earthquakes at rates comparable to similar magnitude earthquakes. Surface waves are associated with regional and long-range dynamic triggering, but we note that surface waves with low enough frequency to penetrate to depths where most aftershocks of the 1992 M = 5.7 Little Skull Mountain main shock occurred (∼12 km) would not have developed significant amplitude within a 50-km radius. We therefore focus on the best candidate phases to cause local dynamic triggering, direct waves that pass through observed near-source aftershock clusters. We examine these phases, which arrived at the nearest (200–270 km) broadband station before the surface wave train and could thus be isolated for study. Direct comparison of spectral amplitudes of presurface wave arrivals shows that M ∼ 5 explosions and earthquakes deliver the same peak dynamic stresses into the near-source crust. We conclude that a static stress change model can readily explain observed aftershock patterns, whereas it is difficult to attribute near-source triggering to a dynamic process because of the dearth of aftershocks near large explosions.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Frictional strengths of talc-serpentine and talc-quartz mixtures

Talc is a constituent of faults in a variety of settings, and it may be an effective weakening agent depending on its abundance and distribution within a fault. We conducted frictional strength experiments under hydrothermal conditions to determine the effect of talc on the strengths of synthetic gouges of lizardite and antigorite serpentinites and of quartz. Small amounts of talc weaken serpentinite gouges substantially more than predicted by simple weight averaging. In comparison, mixtures of quartz and talc show a linear trend of strength reduction at talc concentrations 15 wt % and enhanced weakening at higher concentrations. All of the strength data are fit by a modified version of the Reuss mixing law that allows for the dominance of one mineral over the other. The difference in the behavior of serpentinite-talc and quartz-talc mixtures at low talc concentrations is a reflection of their different textures. Lizardite, antigorite, and talc all have platy habits, and displacement within gouges composed of these minerals is localized to narrow shears along which the platy grains have rotated into alignment with the shear surfaces. The shears in the mixed phyllosilicate gouges maximize the proportion of the weaker mineral within them. When mixed with a strong, rounded mineral such as quartz, some minimum concentration of talc is needed to form connected pathways that enhance strength reductions. The typical development of talc by the reaction of Si-rich fluids with serpentinite or dolomite would tend to localize its occurrence in a natural fault and result in enhanced weakening.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Global map of eolian features on Mars

Ten basic categories of eolian features on Mars were identified from a survey of Mariner 9 and Viking orbiter images. The ten features mapped are (1) light streaks (including frost streaks), (2) dark streaks, (3) sand sheets or splotches, (4) barchari dunes, (5) transverse dunes, (6) crescentic dunes, (7) anomalous dunes, (8) yardangs, (9) wind grooves, and (10) deflation pits. The features were mapped in groups, not as individual landforms, and recorded according to their geographic positions and orientations on maps of 1:12.5 million or 1:25 million scale. In the north polar region, light and dark streaks indicate winds from the west and northwest. Frost streaks show northeast and northwest winds. Barchan dunes show southwest and west winds. Transverse dunes show east and west winds. Local erosional features show winds from the northwest or southeast and northeast or southwest. In the middle and low northern latitudes, streaks show northeast winds; yardangs and deflation pits show mostly easterly and northeasterly winds. In the low southern latitudes, light streaks and dunes record northwest and northeast winds; dark streaks record southeast winds. In high southern latitudes, most streaks and dunes record southeast and east winds. In the south polar region, light and frost streaks record southwest and northwest winds, whereas dark streaks and transverse dunes show southeast winds. The patterns recorded for ephemeral features conform to global conditions of strong southern spring and summer (northern fall and winter) wind circulation. Erosional features in bedrock indicate long-term and perhaps ancient wind trends, whereas depositional features may record relatively more recent winds. Wind directions indicated by yardangs, pits, and other erosional features seldom correspond to those shown by streaks. These erosional features indicating winds of different directions than those that formed the streaks may have been carved when the effective regional winds on Mars were different due to relative hemisphere insolation differences, which are related to obliquity variations. Alternatively, many yardangs and pits may be largely controlled by structure and carved along surface structural patterns or features that channeled both ancient and present-day winds. Deflation pits in the mantled terrain may contain the best record of ancient wind trends. These pits seem to have stratigraphically related orientations. Different stratigraphic units composing the mantle contain wind erosion pits that possibly record large-scale wind direction changes because of long-term changes in the obliquity of the Mars.

Journal of Geophysical Research Solid Earth↗

USGS1997

Three-dimensional upper crustal velocity structure beneath San Francisco Peninsula, California

This paper presents new seismic data from, and crustal models of the San Francisco Peninsula. In much of central California the San Andreas fault juxtaposes the Cretaceous granitic Salinian terrane on its west and the Late Mesozoic/Early Tertiary Franciscan Complex on its east. On San Francisco Peninsula, however, the present-day San Andreas fault is completely within a Franciscan terrane, and the Pilarcitos fault, located southwest of the San Andreas, marks the Salinian-Franciscan boundary. This circumstance has evoked two different explanations: either the Pilarcitos is a thrust fault that has pushed Franciscan rocks over Salinian rocks or the Pilarcitos is a transform fault that has accommodated significant right-lateral slip. In an effort to better resolve the subsurface structure of the peninsula faults, we established a temporary network of 31 seismographs arrayed across the San Andreas fault and the subparallel Pilarcitos fault at ∼1–2 km spacings. These instruments were deployed during the first 6 months of 1995 and recorded local earthquakes, air gun sources set off in San Francisco Bay, and explosive sources. Travel times from these sources were used to augment earthquake arrival times recorded by the Northern California Seismic Network and were inverted for three-dimensional velocity structure. Results show lateral velocity changes at depth (∼0.5–7 km) that correlate with downward vertical projections of the surface traces of the San Andreas and Pilarcitos faults. We thus interpret the faults as high-angle to vertical features (constrained to a 70°–110° dip range). From this we conclude that the Pilarcitos fault is probably an important strike-slip fault that accommodated much of the right-lateral plate boundary strain on the peninsula prior to the initiation of the modern-day San Andreas fault in this region sometime after about 3.0 m.y. ago.

Journal of Geophysical Research B: Solid Earth↗