SEARCH · Search USGS

Results for “Earth Surface Dynamics”

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.

262 records · Page 15Linked to original sources

USGS2006

Dating young geomorphic surfaces using age of colonizing Douglas fir in southwestern Washington and northwestern Oregon, USA

Dating of dynamic, young (<500 years) geomorphic landforms, particularly volcanofluvial features, requires higher precision than is possible with radiocarbon dating. Minimum ages of recently created landforms have long been obtained from tree-ring ages of the oldest trees growing on new surfaces. But to estimate the year of landform creation requires that two time corrections be added to tree ages obtained from increment cores: (1) the time interval between stabilization of the new landform surface and germination of the sampled trees (germination lag time or GLT); and (2) the interval between seedling germination and growth to sampling height, if the trees are not cored at ground level. The sum of these two time intervals is the colonization time gap (CTG). Such time corrections have been needed for more precise dating of terraces and floodplains in lowland river valleys in the Cascade Range, where significant eruption-induced lateral shifting and vertical aggradation of channels can occur over years to decades, and where timing of such geomorphic changes can be critical to emergency planning. Earliest colonizing Douglas fir (Pseudotsuga menziesii) were sampled for tree-ring dating at eight sites on lowland (<750 m a.s.l.), recently formed surfaces of known age near three Cascade volcanoes - Mount Rainier, Mount St. Helens and Mount Hood - in southwestern Washington and northwestern Oregon. Increment cores or stem sections were taken at breast height and, where possible, at ground level from the largest, oldest-looking trees at each study site. At least ten trees were sampled at each site unless the total of early colonizers was less. Results indicate that a correction of four years should be used for GLT and 10 years for CTG if the single largest (and presumed oldest) Douglas fir growing on a surface of unknown age is sampled. This approach would have a potential error of up to 20 years. Error can be reduced by sampling the five largest Douglas fir instead of the single largest. A GLT correction of 5 years should be added to the mean ring-count age of the five largest trees growing on the surface being dated, if the trees are cored at ground level. This correction would have an approximate error of ??5 years. If the trees are cored at about 1.4 m above the round surface (breast height), a CTG correction of 11 years should be added to the mean age of the five sampled trees (with an error of about ??7 years).

Earth Surface Processes and Landforms↗

USGS2023

A model integrating satellite-derived shoreline observations for predicting fine-scale shoreline response to waves and sea-level rise across large coastal regions

Satellite-derived shoreline observations combined with dynamic shoreline models enable fine-scale predictions of coastal change across large spatiotemporal scales. Here, we present a satellite-data-assimilated, “littoral-cell”-based, ensemble Kalman-filter shoreline model to predict coastal change and uncertainty due to waves, sea-level rise (SLR), and other natural and anthropogenic processes. We apply the developed ensemble model to the entire California coastline (approximately 1,760 km), much of which is sparsely monitored with traditional survey methods (e.g., Lidar/GPS). Water-level-corrected, satellite-derived shoreline observations (obtained from the CoastSat toolbox) offer a nearly unbiased representation of in situ surveyed shorelines (e.g., mean sea-level elevation contours) at Ocean Beach, San Francisco. We demonstrate that model calibration with satellite observations during a 20-year hindcast period (1995–2015) provides nearly equivalent model forecast accuracy during a validation period (2015–2020) compared to model calibration with monthly in situ observations at Ocean Beach. When comparing model-predicted shoreline positions to satellite-derived observations, the model achieves an accuracy of <10 m RMSE for nearly half of the entire California coastline for the validation period. The calibrated/validated model is then applied for multi-decadal simulations of shoreline change due to projected wave and sea-level conditions, while holding the model parameters fixed. By 2100, the model estimates that 24%–75% of California's beaches may become completely eroded due to SLR scenarios of 1.0–3.0 m, respectively. The satellite-data-assimilated modeling system presented here is generally applicable to a variety of coastal settings around the world owing to the global coverage of satellite imagery.

JGR Earth Surface↗

USGS2010

Geochronology and paleoclimatic implications of the last deglaciation of the Mauna Kea Ice Cap, Hawaii

We present new 3He surface exposure ages on moraines and bedrock near the summit of Mauna Kea, Hawaii, which refine the age of the Mauna Kea Ice Cap during the Local Last Glacial Maximum (LLGM) and identify a subsequent fluctuation of the ice margin. The 3He ages, when combined with those reported previously, indicate that the local ice-cap margin began to retreat from its LLGM extent at 20.5??2.5ka, in agreement with the age of deglaciation determined from LLGM moraines elsewhere in the tropics. The ice-cap margin receded to a position at least 3km upslope for ~4.5-5.0kyr before readvancing nearly to its LLGM extent. The timing of this readvance at ~15.4ka corresponds to a large reduction of the Atlantic meridional overturning circulation (AMOC) following Heinrich Event 1. Subsequent ice-margin retreat began at 14.6??1.9ka, corresponding to a rapid resumption of the AMOC and onset of the B??lling warm interval, with the ice cap melting rapidly to complete deglaciation. Additional 3He ages obtained from a flood deposit date the catastrophic outburst of a moraine-dammed lake roughly coeval with the Younger Dryas cold interval, suggesting a more active hydrological cycle on Mauna Kea at this time. A coupled mass balance and ice dynamics model is used to constrain the climate required to generate ice caps of LLGM and readvance sizes. The depression of the LLGM equilibrium line altitude requires atmospheric cooling of 4.5??1??C, whereas the mass balance modeling indicates an accompanying increase in precipitation of as much as three times that of present. We hypothesize (1) that the LLGM temperature depression was associated with global cooling, (2) that the temperature depression that contributed to the readvance occurred in response to an atmospheric teleconnection to the North Atlantic, and (3) that the precipitation enhancement associated with both events occurred in response to a southward shift in the position of the inter-tropical convergence zone (ITCZ). Such a shift in the ITCZ would have allowed midlatitude cyclones to reach Mauna Kea more frequently which would have increased precipitation at high elevations and caused additional cooling. ?? 2010 Elsevier B.V.

Earth and Planetary Science Letters↗

USGS1997

A dynamic balance between magma supply and eruption rate at Kilauea volcano, Hawaii

The dynamic balance between magma supply and vent output at Kilauea volcano is used to estimate both the volume of magma stored within Kilauea volcano and its magma supply rate. Throughout most of 1991 a linear decline in volume flux from the Kupaianaha vent on Kilauea's east rift zone was associated with a parabolic variation in the elevation of Kilauea's summit as vent output initially exceeded then lagged behind the magma supply to the volcano. The correspondence between summit elevation and tilt established with over 30 years of data provided daily estimates of summit elevation in terms of summit tilt. The minimum in the parabolic variation in summit tilt and elevation (or zero elevation change) occurs when the magma supply to the reservoir from below the volcano equals the magma output from the reservoir to the surface, so that the magma supply rate is given by vent flux on that day. The measurements of vent flux and tilt establish that the magma supply rate to Kilauea volcano on June 19, 1991, was 217,000±10,000 m 3 /d (or 0.079±0.004 km 3 /yr). This is close to the average eruptive rate of 0.08 km 3 /yr between 1958 and 1984. In addition, the predictable response of summit elevation and tilt to each east rift zone eruption near Puu Oo since 1983 shows that summit deformation is also a measure of magma reservoir pressure. Given this, the correlation between the elevation of the Puu Oo lava lake (4 km uprift of Kupaianaha and 18 km from the summit) and summit tilt provides an estimate for magma pressure changes corresponding to summit tilt changes. The ratio of the change in volume to the change in reservoir pressure ( dV / dP ) during vent activity may be determined by dividing the ratio of volume erupted to change in summit tilt ( dV / d tilt ) by the ratio of pressure change to change in summit tilt ( dP / d tilt ). This measure of dV / dP , when combined with laboratory measurements of the bulk modulus of tholeitic melt, provides an estimate of 240±50 km 3 for the volume of Kilauea's magma reservoir. This estimate is much larger than traditional estimates but consistent with seismic tomographic imaging and geophysical modeling of Kilauea's magma system.

Journal of Geophysical Research B: Solid Earth↗

USGS2025

Hysteretic response of suspended-sediment in wildfire affected watersheds of the Pacific Northwest and Southern Rocky Mountains

Wildfires can have a profound impact on hydrosedimentary interactions, or the relationship between sediment and runoff, in forested headwater streams. Quantification of sediment-runoff dynamics at the event scale is integral for understanding source areas and transport of suspended-sediment through a watershed following wildfire. Here we used high-frequency turbidity and stream discharge data, coupled with discrete suspended-sediment measurements, in burned and unburned watersheds in the Southern Rocky Mountains and the western Cascades Range to evaluate the response of fine-grained (clay- and silt-sized particles) suspended-sediment. Hysteresis analysis was conducted on estimated suspended-sediment concentrations (using turbidity as a proxy) and streamflow through measurement of the difference in sediment concentration on the rising and falling limbs of the event hydrograph. All burned watersheds exhibited elevated fine suspended-sediment concentrations relative to concentrations found in pre-fire conditions. Changes to hysteretic response vary and may depend on a watershed's sediment connectivity limitations. Results suggest a watershed's inherent hillslope-to-channel (or lateral) connectivity is the primary factor controlling the relative magnitude of event-driven fine sediment fluxes in watersheds affected by wildfire. While wildfire did promote lateral connectivity through activation of hillslope sources, snowmelt, precipitation characteristics and antecedent conditions were more important drivers of hysteretic response than wildfire. For watersheds influenced by annual snowpack, we identified a predominantly clockwise hysteretic response during snowmelt and counterclockwise events during the late spring and summer months. There were also proportionally more counterclockwise events after wildfire in watersheds with high sediment connectivity. Results suggest contrasting wildfire-related sediment risk potential. Rivers in burned watersheds with high sediment connectivity may pose a higher risk to receiving waterbodies, such as larger tributaries or reservoirs, while rivers with low sediment connectivity may experience long-term sediment-related risk within the watershed above the outlet.

Oregon↗

USGS2007

Geology and Mineral Resources of the East Mojave National Scenic Area, San Bernardino County, California

The rocks of the East Mojave National Scenic Area (EMNSA) record a history of dynamic geologic events that span more than 1,800 million years (m.y.). These geologic events contributed significantly to development of the spectacular vistas and panoramas present in the area today. The oldest rocks underlie much of the northern part of the EMNSA. These rocks were subjected to extreme pressures and temperatures deep in the Earth's crust about 1,700 million years ago (Ma). They were subsequently intruded by granitic magmas from about 1,695 to 1,650 Ma, by additional granitic magmas at about 1,400 Ma and, later, at about 1,100 Ma, by iron-rich magmas that crystallized to form dark igneous rocks termed diabase. Unusual potassium- and magnesium-rich rocks, emplaced at about 1,400 Ma, crop out in a few places within and near the EMNSA. Their distinctive composition results from very small degrees of partial melting of mantle peridotite that was highly enriched in incompatible trace elements. At Mountain Pass, just outside the northeast boundary of the EMNSA, the potassium- and magnesium-rich rocks are accompanied by a rare type of carbonatite, an igneous rock composed of carbonate minerals, that contains high-grade rare earth element mineralization. Subsequent to these igneous-dominated events, sedimentary strata began to be deposited at about 1,000 Ma; mostly sandstone and shale were deposited initially in marine and, less commonly, in continental environments along the west edge of the core of the North American continent. Sedimentation eventually culminated in the widespread deposition of thick marine limestones from about 400 to about 245 Ma. These limestones represent a continental-shelf environment where shallow-water limestone formed to the east and deeper water limestone formed to the west. The end of the formation of these sedimentary deposits probably was caused by uplift of the shelf, which marked the beginning of a long period of tectonic upheaval. At about 170 Ma, widespread emplacement of coarse-grained granitic magmas began again in the region; some of these magmas also erupted as volcanic rocks. Additional episodes of magmatism took place at about 100 Ma and at 75 Ma. Most of the metallic-mineral occurrences in the EMNSA are associated with the igneous rocks that range in age from 170 to 75 Ma. During each of these magmatic events, the previously deposited sedimentary strata were buckled and broken as the entire region, part of a continental-scale fold and thrust belt, underwent crustal shortening and compression. A period of tectonic quiescence characterized the region from about 65 Ma to about 20 Ma. The quiet period ended abruptly with widespread volcanism along the southern and eastern parts of the EMNSA. The major gold deposits in the Castle Mountains are associated with this episode of volcanism. During this volcanic outburst, the crust extended laterally in several areas that border the EMNSA: along the lower Colorado River 65 km to the east, in the Kingston Range 20 km to the north, and in the central Mojave Desert 75 km to the southwest. This extensional deformation is characterized by the superposition of upper-crustal rocks over midcrustal rocks along large flat-lying faults, several of which project beneath rocks now exposed at the surface in the EMNSA. The near-surface rocks of the EMNSA, however, apparently escaped much of this intense extensional deformation. High-angle faults, which cut several of the mountain ranges, possibly have undergone several periods of movement, which date back to approximately 70 to 100 Ma. Some faults are of local importance to the physiographic development of the mountain ranges and intervening basins, and, in places, the faults seem to have localized various kinds of ore bodies and mineral occurrences. Volcanism and extensional deformation waned from 14 to 11 Ma. By approximately 10 Ma, widespread erosion had produced broad erosional dome-shaped mountains in the n

Bulletin↗

USGS2012

Heat flow and energetics of the San Andreas fault zone

Approximately 100 heat flow measurements in the San Andreas fault zone indicate (1) there is no evidence for local factional heating of the main fault trace at any latitude over a 1000-km length from Cape Mendocino to San Bernardino, (2) average heat flow is high (∼2 HFU, ∼80 mW m −2 ) throughout the 550-km segment of the Coast Ranges that encloses the San Andreas fault zone in central California; this broad anomaly falls off rapidly toward the Great Valley to the east, and over a 200-km distance toward the Mendocino Triple Junction to the northwest. As others have pointed out, a local conductive heat flow anomaly would be detectable unless the frictional resistance allocated to heat production on the main trace were ≲100 bars. Frictional work allocated to surface energy of new fractures is probably unimportant, and hydrologic convection is not likely to invalidate the conduction assumption, since the heat discharge by thermal springs near the fault is negligible. Explanations for the low dynamic friction fall into two intergradational classes: those in which the fault is weak all of the time and those in which it is weak only during earthquakes (possibly just large ones). The first class includes faults containing anomalously weak gouge materials and faults containing materials with normal frictional properties under near-lithostatic steady state fluid pressures. In the second class, weakening is caused by the event (for example, a thermally induced increase in fluid pressure, dehydration of clay minerals, or acoustic fluidization). In this class, unlike the first, the average strength and ambient tectonic shear stress may be large, ∼1 kbar, but the stress allocated to elastic radiation (the apparent stress) must be of similar magnitude, an apparent contradiction with seismic estimates. Unless seismic radiation is underestimated for large earthquakes, it is difficult to justify average tectonic stresses on the main trace of the San Andreas fault in excess of ∼200 bars. The development of the broad Coast Range heat flow anomaly southward from Cape Mendocino suggests that heat flow increases by a factor of 2 within 4 m.y. after the passage of the Mendocino Triple Junction. This passage leaves the San Andreas transform fault zone in its wake; the depth of the anomalous sources cannot be much greater than the depth of the seismogenic layer. Some of the anomalous heat may be supplied by conduction from the warmer mantle that must occur south of the Mendocino transform (where there is no subducting slab), and some might be supplied by shear heating in the fault zone. With no contribution from shear heating, extreme mantle upwelling would be required, and asthenosphere conditions should exist today at depths of only ∼20 km in the northernmost Coast Ranges. If there is an appreciable contribution from shear heating, the heat flow constraint implies that the seismogenic layer is partially decoupled at its base and that the basal traction is in the sense that resists right lateral motion on the fault(s). As a result of these basal tractions, the average shearing stress in the seismogenic layer would increase with distance from the main fault, and the seismogenic layer would offer substantial resistance to plate motion even though resistance on the main fault might be negligible. These speculative models have testable consequences.

Journal of Geophysical Research Solid Earth↗

USGS2021

The Boreal-Arctic Wetland and Lake Dataset (BAWLD)

Methane emissions from boreal and arctic wetlands, lakes, and rivers are expected to increase in response to warming and associated permafrost thaw. However, the lack of appropriate land cover datasets for scaling field-measured methane emissions to circumpolar scales has contributed to a large uncertainty for our understanding of present-day and future methane emissions. Here we present the Boreal–Arctic Wetland and Lake Dataset (BAWLD), a land cover dataset based on an expert assessment, extrapolated using random forest modelling from available spatial datasets of climate, topography, soils, permafrost conditions, vegetation, wetlands, and surface water extents and dynamics. In BAWLD, we estimate the fractional coverage of five wetland, seven lake, and three river classes within 0.5 × 0.5∘ grid cells that cover the northern boreal and tundra biomes (17 % of the global land surface). Land cover classes were defined using criteria that ensured distinct methane emissions among classes, as indicated by a co-developed comprehensive dataset of methane flux observations. In BAWLD, wetlands occupied 3.2 × 106 km2 (14 % of domain) with a 95 % confidence interval between 2.8 and 3.8 × 106 km2. Bog, fen, and permafrost bog were the most abundant wetland classes, covering ∼ 28 % each of the total wetland area, while the highest-methane-emitting marsh and tundra wetland classes occupied 5 % and 12 %, respectively. Lakes, defined to include all lentic open-water ecosystems regardless of size, covered 1.4 × 106 km2 (6 % of domain). Low-methane-emitting large lakes (>10 km2) and glacial lakes jointly represented 78 % of the total lake area, while high-emitting peatland and yedoma lakes covered 18 % and 4 %, respectively. Small (<0.1 km2) glacial, peatland, and yedoma lakes combined covered 17 % of the total lake area but contributed disproportionally to the overall spatial uncertainty in lake area with a 95 % confidence interval between 0.15 and 0.38 × 106 km2. Rivers and streams were estimated to cover 0.12 × 106 km2 (0.5 % of domain), of which 8 % was associated with high-methane-emitting headwaters that drain organic-rich landscapes. Distinct combinations of spatially co-occurring wetland and lake classes were identified across the BAWLD domain, allowing for the mapping of “wetscapes” that have characteristic methane emission magnitudes and sensitivities to climate change at regional scales. With BAWLD, we provide a dataset which avoids double-accounting of wetland, lake, and river extents and which includes confidence intervals for each land cover class. As such, BAWLD will be suitable for many hydrological and biogeochemical modelling and upscaling efforts for the northern boreal and arctic region, in particular those aimed at improving assessments of current and future methane emissions. Data are freely available at https://doi.org/10.18739/A2C824F9X (Olefeldt et al., 2021).

Earth System Science Data↗

USGS2016

Relations between rainfall–runoff-induced erosion and aeolian deposition at archaeological sites in a semi-arid dam-controlled river corridor

Process dynamics in fluvial-based dryland environments are highly complex with fluvial, aeolian, and alluvial processes all contributing to landscape change. When anthropogenic activities such as dam-building affect fluvial processes, the complexity in local response can be further increased by flood- and sediment-limiting flows. Understanding these complexities is key to predicting landscape behavior in drylands and has important scientific and management implications, including for studies related to paleoclimatology, landscape ecology evolution, and archaeological site context and preservation. Here we use multi-temporal LiDAR surveys, local weather data, and geomorphological observations to identify trends in site change throughout the 446-km-long semi-arid Colorado River corridor in Grand Canyon, Arizona, USA, where archaeological site degradation related to the effects of upstream dam operation is a concern. Using several site case studies, we show the range of landscape responses that might be expected from concomitant occurrence of dam-controlled fluvial sand bar deposition, aeolian sand transport, and rainfall-induced erosion. Empirical rainfall-erosion threshold analyses coupled with a numerical rainfall–runoff–soil erosion model indicate that infiltration-excess overland flow and gullying govern large-scale (centimeter- to decimeter-scale) landscape changes, but that aeolian deposition can in some cases mitigate gully erosion. Whereas threshold analyses identify the normalized rainfall intensity (defined as the ratio of rainfall intensity to hydraulic conductivity) as the primary factor governing hydrologic-driven erosion, assessment of false positives and false negatives in the dataset highlight topographic slope as the next most important parameter governing site response. Analysis of 4+ years of high resolution (four-minute) weather data and 75+ years of low resolution (daily) climate records indicates that dryland erosion is dependent on short-term, storm-driven rainfall intensity rather than cumulative rainfall, and that erosion can occur outside of wet seasons and even wet years. These results can apply to other similar semi-arid landscapes where process complexity may not be fully understood.

Arizona↗

USGS2012

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗