USGS Science⌕ Search

SEARCH · USGS Science

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 271 records · Page 15Linked to original sources
USGS2009

Volcano-tectonic implications of 3-D velocity structures derived from joint active and passive source tomography of the island of Hawaii

We present a velocity model of the onshore and offshore regions around the southern part of the island of Hawaii, including southern Mauna Kea, southeastern Hualalai, and the active volcanoes of Mauna Loa, and Kilauea, and Loihi seamount. The velocity model was inverted from about 200,000 first-arrival traveltime picks of earthquakes and air gun shots recorded at the Hawaiian Volcano Observatory (HVO). Reconstructed volcanic structures of the island provide us with an improved understanding of the volcano-tectonic evolution of Hawaiian volcanoes and their interactions. The summits and upper rift zones of the active volcanoes are characterized by high-velocity materials, correlated with intrusive magma cumulates. These high-velocity materials often do not extend the full lengths of the rift zones, suggesting that rift zone intrusions may be spatially limited. Seismicity tends to be localized seaward of the most active intrusive bodies. Low-velocity materials beneath parts of the active rift zones of Kilauea and Mauna Loa suggest discontinuous rift zone intrusives, possibly due to the presence of a preexisting volcanic edifice, e.g., along Mauna Loa beneath Kilauea's southwest rift zone, or alternatively, removal of high-velocity materials by large-scale landsliding, e.g., along Mauna Loa's western flank. Both locations also show increased seismicity that may result from edifice interactions or reactivation of buried faults. New high-velocity regions are recognized and suggest the presence of buried, and in some cases, previously unknown rift zones, within the northwest flank of Mauna Loa, and the south flanks of Mauna Loa, Hualalai, and Mauna Kea. Copyright 2009 by the American Geophysical Union.

Hawaii↗
USGS2003

Fault trends on the seaward slope of the Aleutian Trench: Implications for a laterally changing stress field tied to a westward increase in oblique convergence

Normal faults along the seaward trench slope (STS) commonly strike parallel to the trench in response to bending of the oceanic plate into the subduction zone. This is not the circumstance for the Aleutian Trench, where the direction of convergence gradually changes westward, from normal to transform motion. GLORIA side-scan sonar images document that the Aleutian STS is dominated by faults striking oblique to the trench, west of 179??E and east of 172??W. These images also show a pattern of east-west trending seafloor faults that are aligned parallel to the spreading fabric defined by magnetic anomalies. The stress-strain field along the STS is divided into two domains west and east, respectively, of 179??E. Over the western domain, STS faults and nodal planes of earthquakes are oriented oblique (9??-46??) to the trench axis and (69??-90??) to the magnetic fabric. West of 179??E, STS fault strikes change by 36?? from the E-W trend of STS where the trench-parallel slip gets larger than its orthogonal component of convergence. This rotation indicates that horizontal stresses along the western domain of the STS are deflected by the increasing obliquity in convergence. An analytical model supports the idea that strikes of STS faults result from a superposition of stresses associated with the dextral shear couple of the oblique convergence and stresses caused by plate bending. For the eastern domain, most nodal planes of earthquakes strike parallel to the outer rise, indicating bending as the prevailing mechanism causing normal faulting. East of 172??W, STS faults strike parallel to the magnetic fabric but oblique (10??-26??) to the axis of the trench. On the basis of a Coulomb failure criterion the trench-oblique strikes probably result from reactivation of crustal faults generated by spreading. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

Tomographic Rayleigh-wave group velocities in the Central Valley, California centered on the Sacramento/San Joaquin Delta

If shaking from a local or regional earthquake in the San Francisco Bay region were to rupture levees in the Sacramento/San Joaquin Delta then brackish water from San Francisco Bay would contaminate the water in the Delta: the source of fresh water for about half of California. As a prelude to a full shear-wave velocity model that can be used in computer simulations and further seismic hazard analysis, we report on the use of ambient noise tomography to build a fundamental-mode, Rayleigh-wave group velocity model for the region around the Sacramento/San Joaquin Delta in the western Central Valley, California. Recordings from the vertical component of about 31 stations were processed to compute the spatial distribution of Rayleigh wave group velocities. Complex coherency between pairs of stations were stacked over 8 months to more than a year. Dispersion curves were determined from 4 to about 18 seconds. We calculated average group velocities for each period and inverted for deviations from the average for a matrix of cells that covered the study area. Smoothing using the first difference is applied. Cells of the model were about 5.6 km in either dimension. Checkerboard tests of resolution, which is dependent on station density, suggest that the resolving ability of the array is reasonably good within the middle of the array with resolution between 0.2 and 0.4 degrees. Overall, low velocities in the middle of each image reflect the deeper sedimentary syncline in the Central Valley. In detail, the model shows several centers of low velocity that may be associated with gross geologic features such as faulting along the western margin of the Central Valley, oil and gas reservoirs, and large cross cutting features like the Stockton arch. At shorter periods around 5.5s, the model’s western boundary between low and high velocities closely follows regional fault geometry and the edge of a residual isostatic gravity low. In the eastern part of the valley, the boundaries of the low velocity zone and gravity anomaly are better aligned at longer periods (around 10.5s) suggesting that the eastern edge of the gravity low is associated with deeper structure. There is a strong correspondence between a low in gravity near the Kirby Hills fault and low velocities from the ambient noise tomography. At longer periods, higher velocities creep in from the east and narrow the overall dimension defined by the lower velocities. Overall, there is a strong correspondence between the shape and location of low velocities in the Rayleigh wave velocity images, and geological and geophysical features.

California↗
USGS2001

Seismic mapping of shallow fault zones in the San Gabriel Mountains from the Los Angeles Region Seismic Experiment, southern California

During the Los Angeles Region Seismic Experiment (LARSE), a reflection/refraction survey was conducted along a profile (line 1) extending from Seal Beach, California, northeastward to the Mojave Desert and crossing the Los Angeles and San Gabriel Valley basins and San Gabriel Mountains. In most shot gathers from the southern and central San Gabriel Mountains, clear secondary arrivals are seen that merge, or appear to merge, with first arrivals at three locations, including the location of the Vincent thrust fault, an exposed late Mesozoic/early Cenozoic megathrust. These secondary arrivals are interpretable as reflections in the shallow crust (<5 km depth) from a concave-upward interface that projects to the surface in the north near the Vincent thrust fault, is offset in its central part at the San Gabriel fault (an old branch of the San Andreas fault), and terminates in the south at 1 to 2 km depth at the southern mountain front. The velocity structure above and below this interface strongly suggests it is the Vincent thrust fault: intermediate velocities (6.2 km/s), consistent with mylonites overlying the Vincent thrust fault, are observed above it; lower velocities (5.8 km/s), consistent with the Pelona Schist underlying the Vincent thrust fault, are observed below it. Problems arise, however, in attempting to match this reflector to the exposed Vincent thrust fault, which is seen in outcrops east of line 1. The Vincent thrust fault is shallower than the reflector in most places. An unmapped structure (steep fault, monocline, or thrust fault) is required between line 1 and the outcrops that either drops the Vincent thrust fault down to the depths of the reflector or repeats the Vincent thrust fault beneath line 1 in the footwall of another thrust fault. An alternative interpretation of the reflector is a deep greenstone horizon within the Pelona Schist, although this alternative is not favored by the velocity structure. Copyright 2001 by the American Geophysical Union.

California↗
USGS2003

Basin-centered asperities in great subduction zone earthquakes: A link between slip, subsidence, and subduction erosion?

Published areas of high coseismic slip, or asperities, for 29 of the largest Circum-Pacific megathrust earthquakes are compared to forearc structure revealed by satellite free-air gravity, bathymetry, and seismic profiling. On average, 71% of an earthquake's seismic moment and 79% of its asperity area occur beneath the prominent gravity low outlining the deep-sea terrace; 57% of an earthquake's asperity area, on average, occurs beneath the forearc basins that lie within the deep-sea terrace. In SW Japan, slip in the 1923, 1944, 1946, and 1968 earthquakes was largely centered beneath five forearc basins whose landward edge overlies the 350??C isotherm on the plate boundary, the inferred downdip limit of the locked zone. Basin-centered coseismic slip also occurred along the Aleutian, Mexico, Peru, and Chile subduction zones but was ambiguous for the great 1964 Alaska earthquake. Beneath intrabasin structural highs, seismic slip tends to be lower, possibly due to higher temperatures and fluid pressures. Kilometers of late Cenozoic subsidence and crustal thinning above some of the source zones are indicated by seismic profiling and drilling and are thought to be caused by basal subduction erosion. The deep-sea terraces and basins may evolve not just by growth of the outer arc high but also by interseismic subsidence not recovered during earthquakes. Basin-centered asperities could indicate a link between subsidence, subduction erosion, and seismogenesis. Whatever the cause, forearc basins may be useful indicators of long-term seismic moment release. The source zone for Cascadia's 1700 A.D. earthquake contains five large, basin-centered gravity lows that may indicate potential asperities at depth. The gravity gradient marking the inferred downdip limit to large coseismic slip lies offshore, except in northwestern Washington, where the low extends landward beneath the coast. Transverse gravity highs between the basins suggest that the margin is seismically segmented and could produce a variety of large earthquakes. Published in 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Quantifying probabilities of volcanic events: The example of volcanic hazard at Mount Vesuvius

We describe an event tree scheme to quantitatively estimate both long- and short-term volcanic hazard. The procedure is based on a Bayesian approach that produces a probability estimation of any possible event in which we are interested and can make use of all available information including theoretical models, historical and geological data, and monitoring observations. The main steps in the procedure are (1) to estimate an a priori probability distribution based upon theoretical knowledge, (2) to modify that using past data, and (3) to modify it further using current monitoring data. The scheme allows epistemic and aleatoric uncertainties to be dealt with in a formal way, through estimation of probability distributions at each node of the event tree. We then describe an application of the method to the case of Mount Vesuvius. Although the primary intent of the example is to illustrate the methodology, one result of this application merits special mention. The present emergency response plan for Mount Vesuvius is referenced to a maximum expected event (MEE), the largest out of all the possible eruptions within the next few decades. Our calculation suggest that there is a nonnegligible (1-20%) chance that the next eruption could be larger than that stipulated in the present MEE. The methodology allows all assumptions and thresholds to be clearly identified and provides a rational means for their revision if new data or information are obtained. Copyright 2004 by the American Geophysical Union.

Mount Vesuvius↗
USGS2005

Forearc structure beneath southwestern British Columbia: A three-dimensional tomographic velocity model

This paper presents a three-dimensional compressional wave velocity model of the forearc crust and upper mantle and the subducting Juan de Fuca plate beneath southwestern British Columbia and the adjoining straits of Georgia and Juan de Fuca. The velocity model was constructed through joint tomographic inversion of 50,000 first-arrival times from earthquakes and active seismic sources. Wrangellia rocks of the accreted Paleozoic and Mesozoic island arc assemblage underlying southern Vancouver Island in the Cascadia forearc are imaged at some locations with higher than average lower crustal velocities of 6.5-7.2 km/s, similar to observations at other island arc terranes. The mafic Eocene Crescent terrane, thrust landward beneath southern Vancouver Island, exhibits crustal velocities in the range of 6.0-6.7 km/s and is inferred to extend to a depth of more than 20 km. The Cenozoic Olympic Subduction Complex, an accretionary prism thrust beneath the Crescent terrane in the Olympic Peninsula, is imaged as a low-velocity wedge to depths of at least 20 km. Three zones with velocities of 7.0-7.5 km/s, inferred to be mafic and/or ultramafic units, lie above the subducting Juan de Fuca plate at depths of 25-35 km. The forearc upper mantle wedge beneath southeastern Vancouver Island and the Strait of Georgia exhibits low velocities of 7.2-7.5 km/s, inferred to correspond to ???20% serpentinization of mantle peridotites, and consistent with similar observations in other warm subduction zones. Estimated dip of the Juan de Fuca plate beneath southern Vancouver Island is ???11??, 16??, and 27?? at depths of 30, 40, and 50 km, respectively. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2006

Geomorphic control of landscape carbon accumulation

We use the CREEP process-response model to simulate soil organic carbon accumulation in an undisturbed prairie site in Iowa. Our primary objectives are to identify spatial patterns of carbon accumulation, and explore the effect of erosion on basin-scale C accumulation. Our results point to two general findings. First, redistribution of soil carbon by erosion results in a net increase in basin-wide carbon storage relative to a noneroding environment. Landscape-average mean residence times are increased in an eroding landscape owing to the burial/preservation of otherwise labile C. Second, field observations taken along a slope transect may overlook significant intraslope variations in carbon accumulation. Spatial patterns of modeled deep C accumulation are complex. While surface carbon with its relatively short equilibration time is predictable from surface properties, deep carbon is strongly influenced by the landscape's geomorphic and climatic history, resulting in wide spatial variability. Convergence and divergence associated with upland swales and interfluves result in bimodal carbon distributions in upper and mid slopes; variability in carbon storage within modeled mid slopes was as high as simulated differences between erosional shoulders and depositional valley bottoms. The bimodality of mid-slope C variability in the model suggests that a three-dimensional sampling strategy is preferable over the traditional two-dimensional analog or "catena" approach. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗
USGS2008

Geomorphic evidence for former lobate debris aprons at low latitudes on Mars: Indicators of the Martian paleoclimate

Circumferential depressions enclosing mesas and plateaus in the northern Kasei Valles and in the Tartarus Colles regions of Mars are interpreted as indicators of the former extent of lobate debris aprons, thought to be mixtures of ice and elastic particles. These former lobate debris aprons existed about 1 Ga ago and were embayed by lavas or other flow deposits. After the lobate debris aprons had been removed by sublimation and deflation, topographic depressions with a depth of 50 m and a width of several kilometers were left behind between the mesa or plateau scarp and the solidified flow materials. These depressions or moats are located equatorward of ??30?? at significantly lower latitudes than generally observed for occurrences of modem, intact lobate debris aprons. This observation provides evidence that the paleoclimate at that time was different than today, probably due to a higher averaged obliquity of the planet's rotational axis. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2008

Distributed energy balance modeling of South Cascade Glacier, Washington and assessment of model uncertainty

We have developed a physically based, distributed surface energy balance model to simulate glacier mass balance under meteorological and climatological forcing. Here we apply the model to estimate summer ablation on South Cascade Glacier, Washington, for the 2004 and 2005 mass balance seasons. To arrive at optimal mass balance simulations, we investigate and quantify model uncertainty associated with selecting from a range of physical parameter values that are not commonly measured in glaciological mass balance field studies. We optimize the performance of the model by varying values for atmospheric transmissivity, the albedo of surrounding topography, precipitation-elevation lapse rate, surface roughness for turbulent exchange of momentum, and snow albedo aging coefficient. Of these the snow aging parameter and precipitation lapse rates have the greatest influence on the modeled ablation. We examined model sensitivity to varying parameters by performing an additional 103 realizations with parameters randomly chosen over a ??5% range centered about the optimum values. The best fit suite of model parameters yielded a net balance of -1.69??0.38 m water equivalent (WE) for the 2004 water year and -2.10??0.30 m WE up to 11 September 2005. The 2004 result is within 3% of the measured value. These simulations account for 91% and 93% of the variance in measured ablation for the respective years. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2008

Modeling the evolution of channel shape: Balancing computational efficiency with hydraulic fidelity

The cross-sectional shape of a natural river channel controls the capacity of the system to carry water off a landscape, to convey sediment derived from hillslopes, and to erode its bed and banks. Numerical models that describe the response of a landscape to changes in climate or tectonics therefore require formulations that can accommodate evolution of channel cross-sectional geometry. However, fully two-dimensional (2-D) flow models are too computationally expensive to implement in large-scale landscape evolution models, while available simple empirical relationships between width and discharge do not adequately capture the dynamics of channel adjustment. We have developed a simplified 2-D numerical model of channel evolution in a cohesive, detachment-limited substrate subject to steady, unidirectional flow. Erosion is assumed to be proportional to boundary shear stress, which is calculated using an approximation of the flow field in which log-velocity profiles are assumed to apply along vectors that are perpendicular to the local channel bed. Model predictions of the velocity structure, peak boundary shear stress, and equilibrium channel shape compare well with predictions of a more sophisticated but more computationally demanding ray-isovel model. For example, the mean velocities computed by the two models are consistent to within ???3%, and the predicted peak shear stress is consistent to within ???7%. Furthermore, the shear stress distributions predicted by our model compare favorably with available laboratory measurements for prescribed channel shapes. A modification to our simplified code in which the flow includes a high-velocity core allows the model to be extended to estimate shear stress distributions in channels with large width-to-depth ratios. Our model is efficient enough to incorporate into large-scale landscape evolution codes and can be used to examine how channels adjust both cross-sectional shape and slope in response to tectonic and climatic forcing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2009

Global distribution of bedrock exposures on Mars using THEMIS high-resolution thermal inertia

We investigate high thermal inertia surfaces using the Mars Odyssey Thermal Emission Imaging System (THEMIS) nighttime temperature images (100 m/pixel spatial sampling). For this study, we interpret any pixel in a THEMIS image with a thermal inertia over 1200 J m -2 K -1 s -1/2 as "bedrock" which represents either in situ rock exposures or rock-dominated surfaces. Three distinct morphologies, ranked from most to least common, are associated with these high thermal inertia surfaces: (1) valley and crater walls associated with mass wasting and high surface slope angles; (2) floors of craters with diameters >25 km and containing melt or volcanics associated with larger, high-energy impacts; and (3) intercrater surfaces with compositions significantly more mafic than the surrounding regolith. In general, bedrock instances on Mars occur as small exposures (less than several square kilometers) situated in lower-albedo (<0.18), moderate to high thermal inertia (>350 J m -2 K -1 s -1/2 ), and relatively dust-free (dust cover index <0.95) regions; however, there are instances that do not follow these generalizations. Most instances are concentrated in the southern highlands, with very few located at high latitudes (poleward of 45 o N and 58 o S), suggesting enhanced mechanical breakdown probably associated with permafrost. Overall, Mars has very little exposed bedrock with only 960 instances identified from 75 o S to 75 o N with likely <3500 km 2 exposed, representing???1% of the total surface area. These data indicate that Mars has likely undergone large-scale surface processing and reworking, both chemically and mechanically, either destroying or masking a majority of the bedrock exposures on the planet. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2011

Seismically observed seiching in the Panama Canal

A large portion of the seismic noise spectrum is dominated by water wave energy coupled into the solid Earth. Distinct mechanisms of water wave induced ground motions are distinguished by their spectral content. For example, cultural noise is generally <1 s period, microseisms dominate the seismic spectrum from periods of 2 to 20 s, and the Earth's "hum" is in the range of 50 to 600 s. We show that in a large lake in the Panama Canal there is an additional source of long-period noise generated by standing water waves, seiches, induced by disturbances such as passing ships and wind pressure. We compare seismic waveforms to water level records and relate these observations to changes in local tilt and gravity due to an oscillating seiche. The methods and observations discussed in this paper provide a first step toward quantifying the impact of water inundation as recorded by seismometers. This type of quantified understanding of water inundation will help in future estimates of similar phenomena such as the seismic observations of tsunami impact. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Venus cartography

The entire surface of the planet Venus is being mapped at global and regional scales (1:50,000,000 through 1:1,500,000) with synthetic aperture radar (SAR), radar altimeter, and radiometer measurements of physical properties from the Magellan spacecraft. The mapping includes SAR image mosaics, shaded relief maps, and topographic contour overlays made from altimetry data and by radargrammetric methods. Methods used include new techniques of radar image processing that became operational as a result of the Magellan mission. The area being mapped (at a resolution of 75 m/pixel) is roughly equivalent to that of Earth, including seafloors. The mapping is designed to support geologic and geophysical investigations.

Journal of Geophysical Research E: Planets↗
USGS1996

Quantitative investigations of the Missouri gravity low: A possible expression of a large, Late Precambrian batholith intersecting the New Madrid seismic zone

Analysis of gravity and magnetic anomaly data helps characterize the geometry and physical properties of the source of the Missouri gravity low, an important cratonic feature of substantial width (about 125 km) and length (> 600 km). Filtered anomaly maps show that this prominent feature extends NW from the Reelfoot rift to the Midcontinent Rift System. Geologic reasoning and the simultaneous inversion of the gravity and magnetic data lead to an interpretation that the gravity anomaly reflects an upper crustal, 11-km-thick batholith with either near vertical or outward dipping boundaries. Considering the modeled characteristics of the batholith, structural fabric of Missouri, and relations of the batholith with plutons and regions of alteration, a tectonic model for the formation of the batholith is proposed. The model includes a mantle plume that heated the crust during Late Precambrian and melted portions of lower and middle crust, from which the low-density granitic rocks forming the batholith were partly derived. The batholith, called the Missouri batholith, may be currently related to the release of seismic energy in the New Madrid seismic zone (earthquake concentrations occur at the intersection of the Missouri batholith and the New Madrid seismic zone). Three qualitative mechanical models are suggested to explain this relationship with seismicity. Copyright 1996 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Crustal structure of a transform plate boundary: San Francisco Bay and the central California continental margin

Wide-angle seismic data collected during the Bay Area Seismic Imaging Experiment provide new glimpses of the deep structure of the San Francisco Bay Area Block and across the offshore continental margin. San Francisco Bay is underlain by a veneer (<300 m) of sediments, beneath which P wave velocities increase rapidly from 5.2 km/s to 6.0 km/s at 7 km depth, consistent with rocks of the Franciscan subduction assemblage. The base of the Franciscan at-15-18 km depth is marked by a strong wide-angle reflector, beneath which lies an 8- to 10-km-thick lower crust with an average velocity of 6.75??0.15 km/s. The lower crust of the Bay Area Block may be oceanic in origin, but its structure and reflectivity indicate that it has been modified by shearing and/or magmatic intrusion. Wide-angle reflections define two layers within the lower crust, with velocities of 6.4-6.6 km/s and 6.9-7.3 km/s. Prominent subhorizontal reflectivity observed at near-vertical incidence resides principally in the lowermost layer, the top of which corresponds to the "6-s reflector" of Brocher et al. [1994]. Rheological modeling suggests that the lower crust beneath the 6-s reflector is the weakest part of the lithosphere; the horizontal shear zone suggested by Furlong et al. [1989] to link the San Andreas and Hayward/Calaveras fault systems may actually be a broad zone of shear deformation occupying the lowermost crust. A transect across the continental margin from the paleotrench to the Hayward fault shows a deep crustal structure that is more complex than previously realized. Strong lateral variability in seismic velocity and wide-angle reflectivity suggests that crustal composition changes across major transcurrent fault systems. Pacific oceanic crust extends 40-50 km landward of the paleotrench but, contrary to prior models, probably does not continue beneath the Salinian Block, a Cretaceous arc complex that lies west of the San Andreas fault in the Bay Area. The thickness (10 km) and high lower-crustal velocity of Pacific oceanic crust suggest that it was underplated by magmatism associated with the nearby Pioneer seamount. The Salinian Block consists of a 15-km-thick layer of velocity 6.0-6.2 km/s overlying a 5-km-thick, high-velocity (7.0 km/s) lower crust that may be oceanic crust, Cretaceous arc-derived lower crust, or a magmatically underplated layer. The strong structural variability across the margin attests to the activity of strike-slip faulting prior to and during development of the transcurrent Pacific/North American plate boundary around 29 Ma. Copyright 1996 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

Estimating hydrodynamic roughness in a wave-dominated environment with a high-resolution acoustic Doppler profiler

Hydrodynamic roughness is a critical parameter for characterizing bottom drag in boundary layers, and it varies both spatially and temporally due to variation in grain size, bedforms, and saltating sediment. In this paper we investigate temporal variability in hydrodynamic roughness using velocity profiles in the bottom boundary layer measured with a high-resolution acoustic Doppler profiler (PCADP). The data were collected on the ebb-tidal delta off Grays Harbor, Washington, in a mean water depth of 9 m. Significant wave height ranged from 0.5 to 3 m. Bottom roughness has rarely been determined from hydrodynamic measurements under conditions such as these, where energetic waves and medium-to-fine sand produce small bedforms. Friction velocity due to current u*c and apparent bottom roughness z0a were determined from the PCADP burst mean velocity profiles using the law of the wall. Bottom roughness kB was estimated by applying the Grant-Madsen model for wave-current interaction iteratively until the model u*c converged with values determined from the data. The resulting kB values ranged over 3 orders of magnitude (10-1 to 10-4 m) and varied inversely with wave orbital diameter. This range of kB influences predicted bottom shear stress considerably, suggesting that the use of time-varying bottom roughness could significantly improve the accuracy of sediment transport models. Bedform height was estimated from kB and is consistent with both ripple heights predicted by empirical models and bedforms in sonar images collected during the experiment. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2007

Crustal deformation of the Yellowstone-Snake River Plain volcano-tectonic system-Campaign and continuous GPS observations, 1987-2004

The Yellowstone-Snake River Plain tectonomagmatic province resulted from Late Tertiary volcanism in western North America, producing three large, caldera-forming eruptions at the Yellowstone Plateau in the last 2 Myr. To understand the kinematics and geodynamics of this volcanic system, the University of Utah conducted seven GPS campaigns at 140 sites between 1987 and 2003 and installed a network of 15 permanent stations. GPS deployments focused on the Yellowstone caldera, the Hebgen Lake and Teton faults, and the eastern Snake River Plain. The GPS data revealed periods of uplift and subsidence of the Yellowstone caldera at rates up to 15 mm/yr. From 1987 to 1995, the caldera subsided and contracted, implying volume loss. From 1995 to 2000, deformation shifted to inflation and extension northwest of the caldera. From 2000 to 2003, uplift continued to the northwest while caldera subsidence was renewed. The GPS observations also revealed extension across the Hebgen Lake fault and fault-normal contraction across the Teton fault. Deformation rates of the Yellowstone caldera and Hebgen Lake fault were converted to equivalent total moment rates, which exceeded historic seismic moment release and late Quaternary fault slip-derived moment release by an order of magnitude. The Yellowstone caldera deformation trends were superimposed on regional southwest extension of the Yellowstone Plateau at up to 4.3 &plusmn; 0.2 mm/yr, while the eastern Snake River Plain moved southwest as a slower rate at 2.1 &plusmn; 0.2 mm/yr. This southwest extension of the Yellowstone-Snake River Plain system merged into east-west extension of the Basin-Range province. Copyright 2007 by the American Geophysical Union.

Idaho, Oregon, Utah, Wyoming↗