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 361 records · Page 20Linked to original sources
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↗
USGS2002

Satellite-based detection of global urban heat-island temperature influence

This study utilizes a satellite-based methodology to assess the urban heat-island influence during warm season months for over 4400 stations included in the Global Historical Climatology Network of climate stations. The methodology includes local and regional satellite retrievals of an indicator of the presence green photosynthetically active vegetation at and around the stations. The difference in local and regional samples of the normalized difference vegetation index (NDVI) is used to estimate differences in mean air temperature. Stations classified as urban averaged 0.90??C (N. Hemisphere) and 0.92??C (S. Hemisphere) warmer than the surrounding environment on the basis of the NDVI-derived temperature estimates. Additionally, stations classified as rural averaged 0.19??C (N. Hemisphere) and 0.16??C (S. Hemisphere) warmer than the surrounding environment. The NDVI-derived temperature estimates were found to be in reasonable agreement with temperature differences observed between climate stations. The results suggest that satellite-derived data sets can be used to estimate the urban heat-island temperature influence on a global basis and that a more detailed analysis of rural stations and their surrounding environment may be necessary to assure that temperature trends derived from assumed rural environments are not influenced by changes in land use/land cover. Copyright 2002 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2004

Association between atmospheric circulation patterns and firn-ice core records from the Inilchek glacierized area, central Tien Shan, Asia

Glacioclimatological research in the central Tien Shan was performed in the summers of 1998 and 1999 on the South Inilchek Glacier at 5100-5460 m. A 14.36 m firn-ice core and snow samples were collected and used for stratigraphic, isotopic, and chemical analyses. The firn-ice core and snow records were related to snow pit measurements at an event scale and to meteorological data and synoptic indices of atmospheric circulation at annual and seasonal scales. Linear relationships between the seasonal air temperature and seasonal isotopic composition in accumulated precipitation were established. Changes in the ??18O air temperature relationship, in major ion concentration and in the ratios between chemical species, were used to identify different sources of moisture and investigate changes in atmospheric circulation patterns. Precipitation over the central Tien Shan is characterized by the lowest ionic content among the Tien Shan glaciers and indicates its mainly marine origin. In seasons of minimum precipitation, autumn and winter, water vapor was derived from the and and semiarid regions in central Eurasia and contributed annual maximal solute content to snow accumulation in Tien Shan. The lowest content of major ions was observed in spring and summer layers, which represent maximum seasonal accumulation when moisture originates over the Atlantic Ocean and Mediterranean and Black Seas. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
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

Role of transient water pressure in quarrying: A subglacial experiment using acoustic emissions

Probably the most important mechanism of glacial erosion is quarrying: the growth and coalescence of cracks in subglacial bedrock and dislodgement of resultant rock fragments. Although evidence indicates that erosion rates depend on sliding speed, rates of crack growth in bedrock may be enhanced by changing stresses on the bed caused by fluctuating basal water pressure in zones of ice-bed separation. To study quarrying in real time, a granite step, 12 cm high with a crack in its stoss surface, was installed at the bed of Engabreen, Norway. Acoustic emission sensors monitored crack growth events in the step as ice slid over it. Vertical stresses, water pressure, and cavity height in the lee of the step were also measured. Water was pumped to the lee of the step several times over 8 days. Pumping initially caused opening of a leeward cavity, which then closed after pumping was stopped and water pressure decreased. During cavity closure, acoustic emissions emanating mostly from the vicinity of the base of the crack in the step increased dramatically. With repeated pump tests this crack grew with time until the step's lee surface was quarried. Our experiments indicate that fluctuating water pressure caused stress thresholds required for crack growth to be exceeded. Natural basal water pressure fluctuations should also concentrate stresses on rock steps, increasing rates of crack growth. Stress changes on the bed due to water pressure fluctuations will increase in magnitude and duration with cavity size, which may help explain the effect of sliding speed on erosion rates. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
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↗
USGS2007

Glacier changes in southeast Alaska and northwest British Columbia and contribution to sea level rise

The digital elevation model (DEM) from the 2000 Shuttle Radar Topography Mission (SRTM) was differenced from a composite DEM based on air photos dating from 1948 to 1987 to detennine glacier volume changes in southeast Alaska and adjoining Canada. SRTM accuracy was assessed at ??5 in through comparison with airborne laser altimetry and control locations measured with GPS. Glacier surface elevations lowered over 95% of the 14,580 km2 glacier-covered area analyzed, with some glaciers thinning as much as 640 in. A combination of factors have contributed to this wastage, including calving retreats of tidewater and lacustrine glaciers and climate change. Many glaciers in this region are particularly sensitive to climate change, as they have large areas at low elevations. However, several tidewater glaciers that had historically undergone calving retreats are now expanding and appear to be in the advancing stage of the tidewater glacier cycle. The net average rate of ice loss is estimated at 16.7 ?? 4.4 km3/yr, equivalent to a global sea level rise contribution of 0.04 ?? 0.01 mm/yr. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
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↗
USGS2000

Chronological refinement of an ice core record at Upper Fremont Glacier in south central North America

The potential to use ice cores from alpine glaciers in the midlatitudes to reconstruct paleoclimatic records has not been widely recognized. Although excellent paleoclimatic records exist for the polar regions, paleoclimatic ice core records are not common from midlatitude locations. An ice core removed from the Upper Fremont Glacier in Wyoming provides evidence for abrupt climate change during the mid-1800s. Volcanic events (Krakatau and Tambora) identified from electrical conductivity measurements (ECM) and isotopic and chemical data from the Upper Fremont Glacier were reexamined to confirm and refine previous chronological estimates of the ice core. At a depth of 152 m the refined age-depth profile shows good agreement (1736 ± 10 A.D.) with the 14C age date (1729 ± 95 A.D.). The δ 18 O profile of the Upper Fremont Glacier (UFG) ice core indicates a change in climate known as the Little Ice Age (LIA). However, the sampling interval for δ 18 O is sufficiently large (20 cm) such that it is difficult to pinpoint the LIA termination on the basis of δ 18 O data alone. Other research has shown that changes in the δ18O variance are generally coincident with changes in ECM variance. The ECM data set contains over 125,000 data points at a resolution of 1 data point per millimeter of ice core. A 999-point running average of the ECM data set and results from f tests indicates that the variance of the ECM data decreases significantly at about 108 m. At this depth, the age-depth profile predicts an age of 1845 A.D. Results indicate the termination of the LIA was abrupt with a major climatic shift to warmer temperatures around 1845 A.D. and continuing to present day. Prediction limits (error bars) calculated for the profile ages are ±10 years (90% confidence level). Thus a conservative estimate for the time taken to complete the LIA climatic shift to present-day climate is about 10 years, suggesting the LIA termination in alpine regions of central North America may have occurred on a relatively short (decadal) timescale. Copyright 2000 by the American Geophysical Union.

Wyoming↗
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↗