Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Applied Earth Sciences”

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 163 records · Page 9Linked to original sources
USGS2015

Slipstream: an early Holocene slump and turbidite record from the frontal ridge of the Cascadia accretionary wedge off western Canada and paleoseismic implications

Slipstream Slump, a well-preserved 3 km wide sedimentary failure from the frontal ridge of the Cascadia accretionary wedge 85 km off Vancouver Island, Canada, was sampled during Canadian Coast Guard Ship (CCGS) John P. Tully cruise 2008007PGC along a transect of five piston cores. Shipboard sediment analysis and physical property logging revealed 12 turbidites interbedded with thick hemipelagic sediments overlying the slumped glacial diamict. Despite the different sedimentary setting, atop the abyssal plain fan, this record is similar in number and age to the sequence of turbidites sampled farther to the south from channel systems along the Cascadia Subduction Zone, with no extra turbidites present in this local record. Given the regional physiographic and tectonic setting, megathrust earthquake shaking is the most likely trigger for both the initial slumping and subsequent turbidity currents, with sediments sourced exclusively from the exposed slump face of the frontal ridge. Planktonic foraminifera picked from the resedimented diamict of the underlying main slump have a disordered cluster of 14 C ages between 12.8 and 14.5 ka BP. For the post-slump stratigraphy, an event-free depth scale is defined by removing the turbidite sediment intervals and using the hemipelagic sediments. Nine 14 C dates from the most foraminifera-rich intervals define a nearly constant hemipelagic sedimentation rate of 0.021 cm/year. The combined age model is defined using only planktonic foraminiferal dates and Bayesian analysis with a Poisson-process sedimentation model. The age model of ongoing hemipelagic sedimentation is strengthened by physical property correlations from Slipstream events to the turbidites for the Barkley Canyon site 40 km south. Additional modelling addressed the possibilities of seabed erosion or loss and basal erosion beneath turbidites. Neither of these approaches achieves a modern seabed age when applying the commonly used regional marine 14 C reservoir age of 800 years (marine reservoir correction Δ R = 400 years). Rather, the top of the core appears to be 400 years in the future. A younger marine reservoir age of 400 years (Δ R = 0 years) brings the top to the present and produces better correlations with the nearby Effingham Inlet paleo-earthquake chronology based only on terrestrial carbon requiring no reservoir correction. The high-resolution dating and facies analysis of Slipstream Slump in this isolated slope basin setting demonstrates that this is also a useful type of sedimentary target for sampling the paleoseismic record in addition to the more studied turbidites from submarine canyon and channel systems. The first 10 turbidites at Slipstream Slump were deposited between 10.8 and 6.6 ka BP, after which the system became sediment starved and only two more turbidites were deposited. The recurrence interval for the inferred frequent early Holocene megathrust earthquakes is 460 ± 140 years, compatible with other estimates of paleoseismic megathrust earthquake occurrence rates along the subduction zone.

Slipstream Slump↗
USGS2008

Interactive visualization to advance earthquake simulation

The geological sciences are challenged to manage and interpret increasing volumes of data as observations and simulations increase in size and complexity. For example, simulations of earthquake-related processes typically generate complex, time-varying data sets in two or more dimensions. To facilitate interpretation and analysis of these data sets, evaluate the underlying models, and to drive future calculations, we have developed methods of interactive visualization with a special focus on using immersive virtual reality (VR) environments to interact with models of Earth's surface and interior. Virtual mapping tools allow virtual "field studies" in inaccessible regions. Interactive tools allow us to manipulate shapes in order to construct models of geological features for geodynamic models, while feature extraction tools support quantitative measurement of structures that emerge from numerical simulation or field observations, thereby enabling us to improve our interpretation of the dynamical processes that drive earthquakes. VR has traditionally been used primarily as a presentation tool, albeit with active navigation through data. Reaping the full intellectual benefits of immersive VR as a tool for scientific analysis requires building on the method's strengths, that is, using both 3D perception and interaction with observed or simulated data. This approach also takes advantage of the specialized skills of geological scientists who are trained to interpret, the often limited, geological and geophysical data available from field observations. ?? Birkhaueser 2008.

Pure and Applied Geophysics↗
USGS2017

Influence of pore pressure change on coseismic volumetric strain

Coseismic strain is fundamentally important for understanding crustal response to changes of stress after earthquakes. The elastic dislocation model has been widely applied to interpreting observed shear deformation caused by earthquakes. The application of the same theory to interpreting volumetric strain, however, has met with difficulty, especially in the far field of earthquakes. Predicted volumetric strain with dislocation model often differs substantially, and sometimes of opposite signs, from observed coseismic volumetric strains. The disagreement suggests that some processes unaccounted for by the dislocation model may occur during earthquakes. Several hypotheses have been suggested, but none have been tested quantitatively. In this paper we first examine published data to highlight the difference between the measured and calculated static coseismic volumetric strains; we then use these data to provide quantitative test of the model that the disagreement may be explained by the change of pore pressure in the shallow crust. The test allows us to conclude that coseismic change of pore pressure may be an important mechanism for coseismic crustal strain and, in the far field, may even be the dominant mechanism. Thus in the interpretation of observed coseismic crustal strain, one needs to account not only for the elastic strain due to fault rupture but also for the strain due to coseismic change of pore pressure.

Earth and Planetary Science Letters↗
USGS2025

Land application of biosolid, livestock, and drilling wastes to US farmland: A potential pathway for the redistribution of contaminants in the environment

In the United States (U.S.), waste byproducts generated from the treatment of municipal waste (biosolids), production of livestock (livestock waste), and drilling of oil and gas wells (drilling waste) are commonly applied to agricultural lands. Although this can be a cost-effective reuse/disposal practice, there is limited research on the potential for contaminant exposures and effects on ecosystems, wildlife, and human health from such land applications. In this study, we conducted extensive chemical, microbial, and toxicity analyses of biosolid, livestock, and drilling wastes just prior to land application on agricultural lands at 34 sites across the U.S. Twenty-two analytical methods were used to determine potential contaminant exposures profiles for 452 organic and 114 inorganic chemicals, nine microbial groups, estrogenicity, and cytotoxicity. Analytical results document unique and substantial chemical, microbial, and toxicity profiles for these land-applied wastes. Of the three waste byproducts, biosolids contained the greatest concentrations of household chemicals, pesticides, pharmaceuticals, per-/polyfluoroalkyl substances, calcium, and phosphorus. Livestock waste contained the greatest concentrations of total and leachable dissolved organic carbon, biogenic hormones, mycotoxins, plant estrogens, total inorganic nitrogen, and potassium. Drilling waste contained the greatest concentrations of BTEX compounds (benzene, toluene, ethylbenzene, and xylenes), polycyclic aromatic hydrocarbons, rare-earth elements, barium, strontium, and uranium–thorium series radioisotopes. Biosolid and livestock wastes had greater culturable heterotrophic bacteria, halophilic bacteria, Escherichia coli ( E. coli ), enterococci, and staphylococci concentrations, and greater microbial diversity than drilling waste. Bioassay analyses indicated that exposure to contaminants in livestock wastes and biosolids could result in estrogenic effects, whereas exposure to contaminants in drilling waste could result in cytotoxic effects. Our study documents that current reuse/disposal practices for biosolid, livestock, and drilling wastes on agricultural lands could provide a potential pathway for the redistribution of unique and complex contaminant mixtures into the environment that have bioactive, endocrine disrupting, and carcinogenic characteristics. Results of this study provide a snapshot of chemical compositions and concentrations that can be used to inform the development of best-management practices to help maximize beneficial reuse of these wastes and minimize risk to the environment and human health.

Environmental Science: Processes & Impacts↗
USGS2011

Flow speed estimated by inverse modeling of sandy sediment deposited by the 29 September 2009 tsunami near Satitoa, east Upolu, Samoa

Sandy deposits from the 29 September 2009 tsunami on the east coast of Upolu, Samoa were investigated to document their characteristics and used to apply an inverse sediment transport model to estimate tsunami flow speed. Sandy deposits 6 to 15 cm thick formed from ~ 25 to ~ 250 m inland. Sedimentary layers in the deposits, that are defined by vertical grain size variation and contacts, are interpreted to have formed during onshore runup of two waves. Deposits at 3 locations (100, 170, and 240 m inland) contained two layers that are predominately normally graded (~ 80%), but contained massive sections (~ 15%) and inversely graded sections (~ 5%) at their bases. About 75% of the total thickness of normally graded intervals exhibits a signature of sediment falling out of suspension at their top. This type of grading, termed suspension grading here, was first recognized in turbidity current deposits and is characterized by the entire distribution shifting finer upwards in a layer as high-settling velocity, coarser material deposits first and low-settling velocity finer material deposits last. The Jaffe and Gelfenbaum (2007) inverse sediment transport model was applied to intervals within layers that exhibited suspension grading to estimate tsunami flow speed and was able to reproduce the general trends of the observed suspension grading. A key unknown input in the modeling is the bottom roughness. For a bottom roughness parameterization using a Manning's n of 0.03 (equivalent to a z 0 ~ 0.006 m for the observed flow depths of 2–3 m) flow speeds calculated for the 2 layers at the 3 locations were 3.8, 3.6, and 3.7 m/s (bottom layer/earlier wave) and 4.4, 4.4, and 4.1 m/s (top layer/later wave) at 100, 170, and 240 m inland, respectively. These estimates are consistent with the ~ 3–8 m/s tsunami flow speed from boulder transport calculations and result in Froude numbers of ~ 0.7–1.0 when maximum measured flow depths are used. Because the inverse model assumes the deposit was formed by sediment falling out of suspension care must be taken to model only intervals of the deposit exhibiting suspension grading. Including intervals deposited by either bedload or suspended load transport convergences result in higher, and sometimes unrealistic, tsunami flow speed estimates.

Satitoa, Upolu↗
USGS2025

Technical note: A low-cost approach to monitoring relative streamflow dynamics in small headwater streams using time lapse imagery and a deep learning model

Despite their ubiquity and importance as freshwater habitat, small headwater streams are under-monitored by existing stream gage networks. To address this gap, we describe a low-cost, non-contact, and low-effort method that enables organizations to monitor relative streamflow dynamics in small headwater streams. The method uses a camera to capture repeat images of the stream from a fixed position. A person then annotates pairs of images, in each case indicating which image has more apparent streamflow or indicating equal flow if no difference is discernible. A deep learning modeling framework called streamflow rank estimation (SRE) is then trained on the annotated image pairs and applied to rank all images from highest to lowest apparent streamflow. From this result a relative hydrograph can be derived. We found that our modeled relative hydrograph dynamics matched the observed hydrograph dynamics well for 11 cameras at 8 streamflow sites in western Massachusetts. Higher performance was observed during the annotation period (median Kendall's Tau rank correlation of 0.75, with a range of 0.6–0.83) than after it (median Kendall's Tau of 0.59, with range 0.34–0.74). We found that annotation performance was generally consistent across the 11 camera sites and 2 individual annotators and was positively correlated with streamflow variability at a site. A scaling simulation determined that model performance improvements were limited after 1000 annotation pairs. Our model's estimates of relative flow, while not equivalent to absolute flow, may still be useful for many applications, such as ecological modeling and calculating event-based hydrological statistics (e.g., the number of out-of-bank floods). We anticipate that this method will be a valuable tool to extend existing stream monitoring networks and provide new insights on dynamic headwater systems.

Massachusetts↗
USGS2024

Rapid simulation of wave runup on morphologically diverse, reef-lined coasts with the BEWARE-2 (Broad-range Estimator of Wave Attack in Reef Environments) meta-process model

Low-lying, tropical, coral-reef-lined coastlines are becoming increasingly vulnerable to wave-driven flooding due to population growth, coral reef degradation, and sea-level rise. Early-warning systems (EWSs) are needed to enable coastal authorities to issue timely alerts and coordinate preparedness and evacuation measures for their coastal communities. At longer timescales, risk management and adaptation planning require robust assessments of future flooding hazard considering uncertainties. However, due to diversity in reef morphologies and complex reef hydrodynamics compared to sandy shorelines, there have been no robust analytical solutions for wave runup to allow for the development of large-scale coastal wave-driven flooding EWSs and risk assessment frameworks for reef-lined coasts. To address the need for fast, robust predictions of runup that account for the natural variability in coral reef morphologies, we constructed the BEWARE-2 (Broad-range Estimator of Wave Attack in Reef Environments) meta-process modeling system. We developed this meta-process model using a training dataset of hydrodynamics and wave runup computed by the XBeach Non-Hydrostatic process-based hydrodynamic model for 440 combinations of water level, wave height, and wave period with 195 representative reef profiles that encompass the natural diversity in real-world fringing coral reef systems. Through this innovation, BEWARE-2 can be applied in a larger range of coastal settings than meta-models that rely on a parametric description of the coral reef geometry. In the validation stage, the BEWARE-2 modeling system produced runup results that had a relative root mean square error of 13% and relative bias of 5% relative to runup simulated by XBeach Non-Hydrostatic for a large range of oceanographic forcing conditions and for diverse reef morphologies (root mean square error and bias 0.63 and 0.26 m, respectively, relative to mean simulated wave runup of 4.85 m). Incorporating parametric modifications in the modeling system to account for variations in reef roughness and beach slope allows for systematic errors (relative bias) in BEWARE-2 predictions to be reduced by a factor of 1.5–6.5 for relatively coarse or smooth reefs and mild or steep beach slopes. This prediction provided by the BEWARE-2 modeling system is faster by 4–5 orders of magnitude than the full, process-based hydrodynamic model and could therefore be integrated into large-scale EWSs for tropical, reef-lined coasts and used for large-scale flood risk assessments.

Natural Hazards and Earth System Sciences↗
USGS2011

Habitat suitability of patch types: a case study of the Yosemite toad

Understanding patch variability is crucial in understanding the spatial population structure of wildlife species, especially for rare or threatened species. We used a well-tested maximum entropy species distribution model (Maxent) to map the Yosemite toad (Anaxyrus (= Bufo) canorus) in the Sierra Nevada mountains of California. Twenty-six environmental variables were included in the model representing climate, topography, land cover type, and disturbance factors (e.g., distances to agricultural lands, fire perimeters, and timber harvest areas) throughout the historic range of the toad. We then took a novel approach to the study of spatially structured populations by applying the species-environmental matching model separately for 49 consistently occupied sites of the Yosemite toad compared to 27 intermittently occupied sites. We found that the distribution of the entire population was highly predictable (AUC = 0.95±0.03 SD), and associated with low slopes, specific vegetation types (wet meadow, alpine-dwarf shrub, montane chaparral, red fir, and subalpine conifer), and warm temperatures. The consistently occupied sites were also associated with these same factors, and they were also highly predictable (AUC = 0.95±0.05 SD). However, the intermittently occupied sites were associated with distance to fire perimeter, a slightly different response to vegetation types, distance to timber harvests, and a much broader set of aspect classes (AUC = 0.90±0.11 SD). We conclude that many studies of species distributions may benefit by modeling spatially structured populations separately. Modeling and monitoring consistently-occupied sites may provide a realistic snapshot of current species-environment relationships, important climatic and topographic patterns associated with species persistence patterns, and an understanding of the plasticity of the species to respond to varying climate regimes across its range. Meanwhile, modeling and monitoring of widely dispersing individuals and intermittently occupied sites may uncover environmental thresholds and human-related threats to population persistence.

Frontiers of Earth Science↗
USGS2003

Lead isotopes in North Pacific deep water - Implications for past changes in input sources and circulation patterns

The sources of non-anthropogenic Pb in seawater have been the subject of debate. Here we present Pb isotope time-series that indicate that the non-anthropogenic Pb budget of the northernmost Pacific Ocean has been governed by ocean circulation and riverine inputs, which in turn have ultimately been controlled by tectonic processes. Despite the fact that the investigated locations are situated within the Asian dust plume, and proximal to extensive arc volcanism, eolian contributions have had little impact. We have obtained the first high-resolution and high-precision Pb isotope time-series of North Pacific deep water from two ferromanganese crusts from the Gulf of Alaska in the NE Pacific Ocean, and from the Detroit Seamount in the NW Pacific Ocean. Both crusts were dated applying 10 Be/9Be ratios and yield continuous time-series for the past 13.5 and 9.6 Myr, respectively. Lead isotopes show a monotonic evolution in 206Pb/204Pb from low values in the Miocene (??? 18.57) to high values at present day (??? 18.84) in both crusts, even though they are separated by more than 3000 km along the Aleutian Arc. The variation exceeds the amplitude found in Equatorial Pacific deep water records by about three-fold. There also is a striking similarity in 207Pb/204Pb and 208Pb/ 204Pb ratios of the two crusts, indicating the existence of a local circulation cell in the sub-polar North Pacific, where efficient lateral mixing has taken place but only limited exchange (in terms of Pb) with deep water from the Equatorial Pacific has occurred. Both crusts display well-defined trends with age in Pb-Pb isotope mixing plots, which require the involvement of at least four distinct Pb sources for North Pacific deep water. The Pb isotope time-series reveal that eolian supplies (volcanic ash and continent-derived loess) have only been of minor importance for the dissolved Pb budget of marginal sites in the deep North Pacific over the past 6 Myr. The two predominant sources have been young volcanic arcs, one located in the northeastern part and one located in the northwestern part of the Pacific margin, from where material has been eroded and delivered to the ocean, most likely via riverine pathways. ?? 2003 Elsevier Science B.V. All rights reserved.

Earth and Planetary Science Letters↗
USGS2014

Hydroclimatic regimes: a distributed water-balance framework for hydrologic assessment, classification, and management

Runoff-based indicators of terrestrial water availability are appropriate for humid regions, but have tended to limit our basic hydrologic understanding of drylands – the dry-subhumid, semiarid, and arid regions which presently cover nearly half of the global land surface. In response, we introduce an indicator framework that gives equal weight to humid and dryland regions, accounting fully for both vertical (precipitation + evapotranspiration) and horizontal (groundwater + surface-water) components of the hydrologic cycle in any given location – as well as fluxes into and out of landscape storage. We apply the framework to a diverse hydroclimatic region (the conterminous USA) using a distributed water-balance model consisting of 53 400 networked landscape hydrologic units. Our model simulations indicate that about 21% of the conterminous USA either generated no runoff or consumed runoff from upgradient sources on a mean-annual basis during the 20th century. Vertical fluxes exceeded horizontal fluxes across 76% of the conterminous area. Long-term-average total water availability (TWA) during the 20th century, defined here as the total influx to a landscape hydrologic unit from precipitation, groundwater, and surface water, varied spatially by about 400 000-fold, a range of variation ~100 times larger than that for mean-annual runoff across the same area. The framework includes but is not limited to classical, runoff-based approaches to water-resource assessment. It also incorporates and reinterprets the green- and blue-water perspective now gaining international acceptance. Implications of the new framework for several areas of contemporary hydrology are explored, and the data requirements of the approach are discussed in relation to the increasing availability of gridded global climate, land-surface, and hydrologic data sets.

Hydrology and Earth System Sciences↗
USGS2016

Towards simplification of hydrologic modeling: Identification of dominant processes

The Precipitation–Runoff Modeling System (PRMS), a distributed-parameter hydrologic model, has been applied to the conterminous US (CONUS). Parameter sensitivity analysis was used to identify: (1) the sensitive input parameters and (2) particular model output variables that could be associated with the dominant hydrologic process(es). Sensitivity values of 35 PRMS calibration parameters were computed using the Fourier amplitude sensitivity test procedure on 110 000 independent hydrologically based spatial modeling units covering the CONUS and then summarized to process (snowmelt, surface runoff, infiltration, soil moisture, evapotranspiration, interflow, baseflow, and runoff) and model performance statistic (mean, coefficient of variation, and autoregressive lag 1). Identified parameters and processes provide insight into model performance at the location of each unit and allow the modeler to identify the most dominant process on the basis of which processes are associated with the most sensitive parameters. The results of this study indicate that: (1) the choice of performance statistic and output variables has a strong influence on parameter sensitivity, (2) the apparent model complexity to the modeler can be reduced by focusing on those processes that are associated with sensitive parameters and disregarding those that are not, (3) different processes require different numbers of parameters for simulation, and (4) some sensitive parameters influence only one hydrologic process, while others may influence many

Hydrology and Earth System Sciences↗
USGS2005

Geologic and geophysical maps of the Las Vegas 30' x 60' quadrangle, Clark and Nye counties, Nevada, and Inyo County, California

Las Vegas and Pahrump are two of the fastest growing cities in the US, and the shortage of water looms as among the greatest future problems for these cities. These new maps of the Las Vegas 30 x 60-minute quadrangle provide a geologic and geophysical framework and fundamental earth science database needed to address societal issues such as ground water supply and contamination, surface flood, landslide, and seismic hazards, and soil properties and their changing impact by and on urbanization. The mountain ranges surrounding Las Vegas and Pahrump consist of Mesozoic, Paleozoic and Proterozoic rocks. A majority of these rocks are Paleozoic carbonate rocks that are part of Nevada's carbonate rock aquifer province. The Spring Mountains represent a major recharge site in the province, where maximum altitude is 3,632 m (Charleston Peak) above sea level. Rocks in the Sheep and Las Vegas Ranges and Spring Mountains contain correlative, northeast-striking, southeast-verging thrust faults that are part of the Cretaceous, Sevier orogenic belt. These thrusts were offset during the Miocene by the Las Vegas Valley shear system (LVVSZ). We conducted new mapping in the Blue Diamond area, highlighting refined work on the Bird Spring thrust, newly studied ancient landslides, and gravity-slide blocks. We conducted new mapping in the Las Vegas Range and mapped previously unrecognized structures such as the Valley thrust and fold belt; recognition of these structures has led to a refined correlation of Mesozoic thrust faults across the LVVSZ. New contributions in the quadrangle also include a greatly refined stratigraphy of Paleozoic bedrock units based on conodont biostragraphy. We collected over 200 conodont samples in the quadrangle and established stratigraphic reference sections used to correlate units across the major Mesozoic thrust faults. Quaternary deposits cover about half of the map area and underlie most of the present urbanized area. Deposits consist of large coalescing alluvial fans that grade downslope to extensive areas of fine-grained sediment indicative of groundwater-discharge during the Pleistocene. In the central areas of Las Vegas and Pahrump valleys, Quaternary fault scarps associated with past ground-water discharge deposits suggest a genetic relationship. In collaboration with NBMG and University of Nevada, a variety of ages of gravelly alluvium are newly mapped using surficial characteristics and soil development, along with reassessment of previously published mapping during compilation. Reconnaissance geochronology (thermoluminescence and U-series) of eolian and authigenic components of surficial and buried soils and spring deposits is applied to test hypotheses of geomorphic and hydrologic response to climate change over the past 100 k.y.). The major structure in the Las Vegas quadrangle is the LVVSZ. Because the LVVSZ is concealed by thick basin-fill deposits of Quaternary and Tertiary age, it was characterized primarily based on geophysics. Likewise, the newly described State line fault system in Pahrump Valley has also been characterized by geophysics, where geophysically inferred structures correlate remarkably with surface structures defined by our new geologic mapping in the Mound Spring and Hidden Hills Ranch 7.5-minute quadrangles.

Scientific Investigations Map↗
USGS1996

Velocity structure of a bottom simulating reflector offshore Peru: Results from full waveform inversion

Much of our knowledge of the worldwide distribution of submarine gas hydrates comes from seismic observations of Bottom Simulating Reflectors (BSRs). Full waveform inversion has proven to be a reliable technique for studying the fine structure of BSRs using the compressional wave velocity. We applied a non-linear full waveform inversion technique to a BSR at a location offshore Peru. We first determined the large-scale features of seismic velocity variations using a statistical inversion technique to maximise coherent energy along travel-time curves. These velocities were used for a starting velocity model for the full waveform inversion, which yielded a detailed velocity/depth model in the vicinity of the BSR. We found that the data are best fit by a model in which the BSR consists of a thin, low-velocity layer. The compressional wave velocity drops from 2.15 km/s down to an average of 1.70 km/s in an 18m thick interval, with a minimum velocity of 1.62 km/s in a 6 m interval. The resulting compressional wave velocity was used to estimate gas content in the sediments. Our results suggest that the low velocity layer is a 6-18 m thick zone containing a few percent of free gas in the pore space. The presence of the BSR coincides with a region of vertical uplift. Therefore, we suggest that gas at this BSR is formed by a dissociation of hydrates at the base of the hydrate stability zone due to uplift and subsequently a decrease in pressure.

Earth and Planetary Science Letters↗
USGS2014

Steady incision of Grand Canyon at the million year timeframe: A case for mantle-driven differential uplift

The Grand Canyon region provides an excellent laboratory to examine the interplay between river incision, magmatism, and the geomorphic and tectonic processes that shape landscapes. Here we apply U-series, Ar–Ar, and cosmogenic burial dating of river terraces to examine spatial variations in incision rates along the 445 km length of the Colorado River through Grand Canyon. We also analyze strath terrace sequences that extend to heights of several hundred meters above the river, and integrate these with speleothem constrained maximum incision rates in several reaches to examine any temporal incision variations at the million-year time frame. This new high-resolution geochronology shows temporally steady long-term incision in any given reach of Grand Canyon but significant variations along its length from 160 m/Ma in the east to 101 m/Ma in the west. Spatial and temporal patterns of incision, and the long timescale of steady incision rule out models where geomorphic controls such as climate oscillations, bedrock strength, sediment load effects, or isostatic response to differential denudation are the first order drivers of canyon incision. The incision pattern is best explained by a model of Neogene and ongoing epeirogenic uplift due to an eastward propagating zone of increased upper mantle buoyancy that we infer from propagation of Neogene basaltic volcanism and a strong lateral gradient in modern upper mantle seismic structure.

Arizona↗
USGS2009

Geometry of the Nojima fault at Nojima-Hirabayashi, Japan - I. A simple damage structure inferred from borehole core permeability

The 1995 Kobe (Hyogo-ken Nanbu) earthquake, M = 7.2, ruptured the Nojima fault in southwest Japan. We have studied core samples taken from two scientific drillholes that crossed the fault zone SW of the epicentral region on Awaji Island. The shallower hole, drilled by the Geological Survey of Japan (GSJ), was started 75 m to the SE of the surface trace of the Nojima fault and crossed the fault at a depth of 624 m. A deeper hole, drilled by the National Research Institute for Earth Science and Disaster Prevention (NIED) was started 302 m to the SE of the fault and crossed fault strands below a depth of 1140 m. We have measured strength and matrix permeability of core samples taken from these two drillholes. We find a strong correlation between permeability and proximity to the fault zone shear axes. The half-width of the high permeability zone (approximately 15 to 25 m) is in good agreement with the fault zone width inferred from trapped seismic wave analysis and other evidence. The fault zone core or shear axis contains clays with permeabilities of approximately 0.1 to 1 microdarcy at 50 MPa effective confining pressure (10 to 30 microdarcy at in situ pressures). Within a few meters of the fault zone core, the rock is highly fractured but has sustained little net shear. Matrix permeability of this zone is approximately 30 to 60 microdarcy at 50 MPa effective confining pressure (300 to 1000 microdarcy at in situ pressures). Outside this damage zone, matrix permeability drops below 0.01 microdarcy. The clay-rich core material has the lowest strength with a coefficient of friction of approximately 0.55. Shear strength increases with distance from the shear axis. These permeability and strength observations reveal a simple fault zone structure with a relatively weak fine-grained core surrounded by a damage zone of fractured rock. In this case, the damage zone will act as a high-permeability conduit for vertical and horizontal flow in the plane of the fault. The fine-grained core region, however, will impede fluid flow across the fault.

Pure and Applied Geophysics↗
USGS2004

Mars chronology: Assessing techniques for quantifying surficial processes

Currently, the absolute chronology of Martian rocks, deposits and events is based mainly on crater counting and remains highly imprecise with epoch boundary uncertainties in excess of 2 billion years. Answers to key questions concerning the comparative origin and evolution of Mars and Earth will not be forthcoming without a rigid Martian chronology, enabling the construction of a time scale comparable to Earth's. Priorities for exploration include calibration of the cratering rate, dating major volcanic and fluvial events and establishing chronology of the polar layered deposits. If extinct and/or extant life is discovered, the chronology of the biosphere will be of paramount importance. Many radiometric and cosmogenic techniques applicable on Earth and the Moon will apply to Mars after certain baselines (e.g. composition of the atmosphere, trace species, chemical and physical characteristics of Martian dust) are established. The high radiation regime may pose a problem for dosimetry-based techniques (e.g. luminescence). The unique isotopic composition of nitrogen in the Martian atmosphere may permit a Mars-specific chronometer for tracing the time-evolution of the atmosphere and of lithic phases with trapped atmospheric gases. Other Mars-specific chronometers include measurement of gas fluxes and accumulation of platinum group elements (PGE) in the regolith. Putting collected samples into geologic context is deemed essential, as is using multiple techniques on multiple samples. If in situ measurements are restricted to a single technique it must be shown to give consistent results on multiple samples, but in all cases, using two or more techniques (e.g. on the same lander) will reduce error. While there is no question that returned samples will yield the best ages, in situ techniques have the potential to be flown on multiple missions providing a larger data set and broader context in which to place the more accurate dates. ?? 2004 Elsevier B.V. All rights reserved.

Earth-Science Reviews↗
USGS1969

Dilemma posed by uranium-series dates on archaeologically significant bones from Valsequillo, Puebla, Mexico

In an attempt to date stone artifacts of Early Man excavated from several sites at the Valsequillo Reservoir, a few kilometers south of Puebla, Mexico, Szabo applied the uranium-series method on bone samples known to be either from the same geologic formation as the sites or in direct association with the artifacts. The geologic context of the bones was studied by Malde, and the archaeological sites were excavated by Irwin-Williams. A date determined for bone associated with an artifact (Caulapan sample M-B-6, see below) agrees with a radiocarbon date for fossil mollusks in the same bed and indicates man's presence more than 20 000 years ago. However, some of these bone dates exceed 200 000 years. Because such dates for man in North America conflict with all prior archaeological evidence here and abroad, we are confronted by a dilemna — either to defend the dates against an onslaught of archaeological thought, or to abandon the uranium method in this application as being so much wasted effort. Faced with these equally undesirable alternatives, and unable to decide where the onus fairly lies (if a choice must be made), we give the uranium-series dates as a possible stimulus for further mutual work in isotopic dating of archaeological material. A sample from the Lindenmeier archaeological site north of Fort Collins and another from a Pleistocene terrace along the Arkansas River, both in Colorado, were also dated.

Puebla↗
USGS2022

Using active source seismology to image the Palos Verdes Fault damage zone as a function of distance, depth, and geology

Fault damage zones provide a window into the non-elastic processes of an earthquake. Geological and seismic tomography methods have been unable to measure damage zones at depth with sufficient spatial sampling to evaluate the relative influence of depth, distance, and lithological variations. Here, we identify and analyze the damage zone of the Palos Verdes Fault offshore southern California using two 3D seismic reflection datasets. We apply a novel algorithm to identify discontinuities attributed to faults and fractures in large seismic volumes and examine the spatial distribution of fault damage in sedimentary rock surrounding the Palos Verdes Fault. Our results show that damage through fracturing is most concentrated around mapped faults and decays exponentially to a distance of ∼2 km, where fracturing reaches a clearly defined and relatively undamaged background for all examined depths and lithologies (450 m to 2.2 km). This decrease in fracturing with distance from the central fault strand exhibits similar functional form to outcrop studies. However, here we extend analysis to distances seldom accessible (∼10 km lateral distance). Separating the data by geologic units we find that the damage decay and background level differs for each unit, with the older and deeper units having higher levels of background fracturing and shallower exponential decays of fracturing with distance from the fault. Surprisingly, these differences in damage decay and background level trade-off result in a consistent damage zone width regardless of lithology or depth. We find that the damage zone has similar decay trends on both sides of the fault. When examining the damage zone at shorter (4 km vs 17 km) along strike distances, the damage zone has a more complex decay trend and at least two strands are resolvable.

California↗