SEARCH · USGS Science

Results for “Near Surface Geophysics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

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

At least 955 records · Page 53Linked to original sources

USGS2007

Questa baseline and pre-mining ground-water quality investigation 15.—Methods of phase II and III well installation and development and results of well logging, hydraulic testing, and water-level measurements in the Red River Valley, New Mexico, 2002-04

In April 2001, the U.S. Geological Survey and the New Mexico Environment Department began a cooperative study to infer the pre-mining ground-water chemistry at the Molycorp molybdenum mine site in the Red River Valley of north- central New Mexico. This report is one in a series of reports that can be used to determine pre-mining ground-water conditions at the mine site. Weathering of hydrothermally altered bedrock in the study area has resulted in steep, highly erosive, and sparsely vegetated scar areas that are clearly visible from the ground and in aerial photographs. Runoff from intense summer rainfall over tributary drainages containing scar areas can transport large quantities of sediment and form debris fans where these tributaries join the Red River. Twenty-nine observation wells were installed in three phases as part of this study in the Red River Valley and tributary drainages. Eight Phase II observation wells were drilled using an air-rotary/hammer rig. Three Phase II and 10 phase III small-diameter wells were installed using a direct-push rig. Lithologic logs were recorded for all eight Phase II drilled wells. Borehole geophysical logging (including natural gamma, induction, and single-detector neutron) was conducted in three Phase II wells. Aquifer tests conducted during 2003 to estimate the hydraulic properties of debris-flow and Red River alluvial deposits in and near Straight Creek included a flow-meter survey, slug tests, and a pumping test. Results of a flow-meter survey in well SC-7A indicated that about 77 percent of the water entered the well from a 10-foot-thick zone near the top of the screened interval and about 23 percent of the water entered the well from a 15-foot-thick zone near the bottom of the screened interval. Slug tests, performed in 11 wells during June 3-5, 2003, indicated that the mean and median estimated hydraulic conductivities for debris-flow deposits were 15.25 and 15.35 feet per day, respectively, for bedrock were 0.12 and 0.08 feet per day, respectively, and for mixed debris flow and Red River alluvium were 73-207 (estimated range) and 80 feet per day. In general, bedrock has the smallest hydraulic conductivity, debris-flow material has the next highest hydraulic conductivity, and mixed debris flow and Red River alluvium has the largest hydraulic conductivity. A pumping test conducted December 3-4, 2003, using well AWWT-1 as the pumped well, and wells AWWT-2, SC-5A, SC-5B, SC-7A, and SC-8A as observation wells, indicated estimated transmissivity of 12,000 to 34,000 feet squared per day and estimated hydraulic conductivity of 230 to 340 feet per day. Water-level measurements in wells SC-6A, SC-7A, SC-8A, and the Hottentot, Hansen, and La Bobita wells show that water levels typically rose rapidly during melting of the winter snowpack in the springtime and then generally declined during the rest of the year. The water-level rise in response to spring snowmelt occurred earlier and was smaller at larger distances from the Red River. Differences between the stage in the Red River and water levels in wells SC-8A and SC-9A, and the absence of water in well SC-9A at the time of well completion, indicate that the Red River has a poor hydraulic connection to the underlying ground-water system and the surface-water system is perched above the ground-water system at this site. Water levels in Phase III wells indicate that the Red River and the shallow ground-water system are connected hydraulically from near wells 4-1D and 4-1S downstream to near wells 2-1 and 2-2 but are poorly connected near the La Bobita well and well 1.

New Mexico↗

USGS2009

Hydrologic Conditions that Influence Streamflow Losses in a Karst Region of the Upper Peace River, Polk County, Florida

The upper Peace River from Bartow to Fort Meade, Florida, is described as a groundwater recharge area, reflecting a reversal from historical groundwater discharge patterns that existed prior to the 1950s. The upper Peace River channel and floodplain are characterized by extensive karst development, with numerous fractures, crevasses, and sinks that have been eroded in the near-surface and underlying carbonate bedrock. With the reversal in groundwater head gradients, river water is lost to the underlying groundwater system through these karst features. An investigation was conducted to evaluate the hydrologic conditions that influence streamflow losses in the karst region of the upper Peace River. The upper Peace River is located in a basin that has been altered substantially by phosphate mining and increases in groundwater use. These alterations have changed groundwater flow patterns and caused streamflow declines through time. Hydrologic factors that have had the greatest influence on streamflow declines in the upper Peace River include the lowering of the potentiometric surfaces of the intermediate aquifer system and Upper Floridan aquifer beneath the riverbed elevation due to below-average rainfall (droughts), increases in groundwater use, and the presence of numerous karst features in the low-water channel and floodplain that enhance the loss of streamflow. Seepage runs conducted along the upper Peace River, from Bartow to Fort Meade, indicate that the greatest streamflow losses occurred along an approximate 2-mile section of the river beginning about 1 mile south of the Peace River at Bartow gaging station. Along the low-water and floodplain channel of this 2-mile section, there are about 10 prominent karst features that influence streamflow losses. Losses from the individual karst features ranged from 0.22 to 16 cubic feet per second based on measurements made between 2002 and 2007. The largest measured flow loss for all the karst features was about 50 cubic feet per second, or about 32 million gallons per day, on June 28, 2002. Streamflow losses varied throughout the year, and were related to seasonal fluctuations in groundwater levels. When groundwater levels were at their lowest level at the end of the dry season (May and June), there was an increased potential for streamflow losses. During this study, the largest streamflow losses occurred at the beginning of the summer rainy season when discharge in the river increased and large volumes of water were needed to replenish unfilled cavities and void spaces in the underlying aquifers. The underlying geology along the upper Peace River and floodplain is highly karstified, and aids in the movement and amount of streamflow that is lost to the groundwater system in this region. Numerous karst features and fractured carbonates and cavernous zones observed in geologic cores and geophysical logs indicate an active, well-connected, groundwater flow system. Aquifer and dye tests conducted along the upper Peace River indicate the presence of cavernous and highly transmissive layers within the floodplain area that can store and transport large volumes of water in underground cavities. A discharge measurement made during this study indicates that the cavernous system associated with Dover Sink can accept over 10 million gallons per day (16 cubic feet per second) of streamflow before the localized aquifer storage volume is replenished and the level of the sink is stabilized.

Scientific Investigations Report↗

USGS2023

Determining seasonal recharge, storage changes, and specific yield using repeat microgravity and water-level measurements in the Mesilla Basin alluvial aquifer, New Mexico, 2016–2018

Increasing water demand and multi-year drought conditions within the Mesilla/Conejos-Médanos Basin near the New Mexico-Texas- Chihuahua border have resulted in diminished surface-water supplies and increased groundwater withdrawals. To better understand recharge to the shallow aquifer, the spatial and temporal groundwater storage changes, and the variability of specific yield (Sy) in the aquifer, seasonal groundwater elevation and repeat microgravity measurements were made during the irrigation release and non-release seasons of 2016, 2017, and 2018 at a network of locations near Las Cruces, New Mexico. The data collected during this investigation were able to capture seasonal change in groundwater elevations and storage from various sources of recharge at multiple sites in the shallow aquifer. Seasonal recharge in the study area was attributed to streamflow, the application and conveyance of irrigation water, and large or sustained precipitation events. However, increasing groundwater gradients in recent decades between piezometers close to the river and those more than a kilometer from the river suggests that recharge from river seepage has become localized at the seasonal scale. Overall, there was a net increase in storage of almost 8.4 cubic hectometers in the study reach between the start and end of the study, largely following the increased surface-water availability and above average precipitation in 2017. Specific yield, estimated by comparing the groundwater-level changes and storage changes at six sites in the study area, ranged from 0.14 (+/− 0.05) to 0.30 (+/− 0.06), which is slightly greater than previously reported estimates (0.10 to 0.25), but still within the error of the estimates. Most of the variability in the estimated storage change, that was not well-correlated with groundwater elevation change, is thought to be from soil moisture in the unsaturated zone. This investigation demonstrates the value of adding repeat microgravity measurements to conventional groundwater monitoring to better understand the sources and extent of recharge as well as the variability of Sy in the aquifer. Continued monitoring, under a variety of available surface water and meteorological conditions, could provide a more comprehensive understanding of the water budget and reduce the specific yield estimation uncertainty. Evaluating water-levels and storage conditions prior to, and following, local recharge events may help managers identify threshold conditions for aquifer storage depletions and recoveries.

New Mexico↗

USGS2000

Geologic mapping of Europa

Galileo data enable the major geological units, structures, and surface features to be identified on Europa. These include five primary units (plains, chaos, band, ridge, and crater materials) and their subunits, along with various tectonic structures such as faults. Plains units are the most widespread. Ridged plains material spans a wide range of geological ages, including the oldest recognizable features on Europa, and appears to represent a style of tectonic resurfacing, rather than cryovolcanism. Smooth plains material typically embays other terrains and units, possibly as a type of fluid emplacement, and is among the youngest material units observed. At global scales, plains are typically mapped as undifferentiated plains material, although in some areas differences can be discerned in the near infrared which might be related to differences in ice grain size. Chaos material is composed of plains and other preexisting materials that have been severely disrupted by inferred internal activity; chaos is characterized by blocks of icy material set in a hummocky matrix. Band material is arrayed in linear, curvilinear, wedge-shaped, or cuspate zones with contrasting albedo and surface textures with respect to the surrounding terrain. Bilateral symmetry observed in some bands and the relationships with the surrounding units suggest that band material forms by the lithosphere fracturing, spreading apart, and infilling with material derived from the subsurface. Ridge material is mapped as a unit on local and some regional maps but shown with symbols at global scales. Ridge material includes single ridges, doublet ridges, and ridge complexes. Ridge materials are considered to represent tectonic processes, possibly accompanied by the extrusion or intrusion of subsurface materials, such as diapirs. The tectonic processes might be related to tidal flexing of the icy lithosphere on diurnal or longer timescales. Crater materials include various interior (smooth central, rough inner, and annular massif) and exterior (continuous ejecta) subunits. Structural features and landforms are shown with conventional symbols. Type localities for the units are identified, along with suggestions for portraying the features on geological maps, including colors and letter abbreviations for material units. Implementing these suggestions by the planetary mapping community would facilitate comparisons of maps for different parts of Europa and contribute to an eventual global synthesis of its complex geology. On the basis of initial mapping results, a stratigraphic sequence is suggested in which ridged plains form the oldest unit on Europa, followed by development of band material and individual ridges. Band materials tend to be somewhat older than ridges, but in many areas the two units formed simultaneously. Similarly, the formation of most chaos follows the development of ridged plains; although chaos is among the youngest materials on Europa, some chaos units might have formed contemporaneously with ridged plains. Smooth plains generally embay all other units and are late-stage in the evolution of the surface. C1 craters are superposed on ridged plains but are crosscut by other materials, including bands and ridges. Most c2 craters postdate all other units, but a few c2 craters are cut by ridge material. C3 craters constitute the youngest recognizable material on Europa. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

USGS1997

Timing of detachment faulting in the Bullfrog Hills and Bare Mountain area, southwest Nevada: Inferences from 40Ar/39Ar, K-Ar, U-Pb, and fission track thermochronology

Crustal extension in the Bullfrog Hills and Bare Mountain area of southwest Nevada is associated with movement along a regional detachment fault. Normal faulting in the upper plate and rapid cooling (denudation) of the lower plate were coeval with Miocene silicic volcanism and with west-northwest transport along the detachment fault. A west-northwest progression of tilting along upper plate normal faults is indicated by ages of the volcanic rocks in relation to angular unconformities. Near the breakaway, tilting in the upper plate occurred between 12.7 and 11.6 Ma, continued less strongly past 10.7 Ma, and was over by 8.2 Ma. Ten to 20 km west of the breakaway, tilting occurred between 10.7 and 10.33 Ma, continued less strongly after 10.33 Ma, and was over by 8.1 Ma. The cooling histories of the lower plate metamorphic rocks were determined by thermochronologic dating methods: K-Ar and 40 Ar/ 39 Ar on muscovite, biotite, and hornblende, 40 Ar/ 39 Ar on K-feldspar, U-Pb on apatite, zircon, and sphene, and fission track on apatite, zircon, and sphene. Lower plate rocks 10 km west of the breakaway cooled slowly from Early Cretaceous lower-amphibolite facies conditions through 350±50° to 300±50°C between 57 and 38 Ma, then cooled rapidly from 205±50° to 120±50°C between 12.6±1.6 and 11.1±1.9 Ma. Lower plate rocks 20 km west of the breakaway cooled slowly from Early Cretaceous upper-amphibolite facies conditions through 500±50°C at 78–67 Ma, passed through 350±50° to 300±50°C between 16.3±0.4 and 10.5±0.3 Ma, then cooled rapidly from 285±50° to 120±50°C between 10.2 and 8.6 Ma. Upper plate tilting and rapid cooling (denudation) of the lower plate occurred simultaneously in the respective areas. The early slow-cooling part of the lower plate thermal histories was probably related to erosion at the Earth's surface, which stripped off about 9 km of material in 50 to 100 m.y. The results indicate an initial fault dip ≥30° and a 12 mm yr −1 west-northwest migration of the locus of rapid tilting in the upper plate.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Reservoir analysis of the Denver earthquakes: A case of induced seismicity

Injection of fluid wastes into the fractured Precambrian crystalline bedrock beneath the Rocky Mountain Arsenal near Denver triggered earthquakes in the 1960's. An analysis, based on the assumption that fluid flow in the fractured reservoir can be approximated by flow in a porous medium, is presented. The configuration and hydrologic properties of the reservoir are determined from two lines of evidence: (1) locations of earthquake hypocenters determined by seismic arrays installed at the Arsenal and (2) observed long-term decline in fluid levels in the injection well. Together these two sets of data indicate that a long, narrow reservoir, aligned in the direction N 60°W, exists. The reservoir is 3.35 km in width, extends 30.5 km to the northwest and infinitely to the southeast, and spans a depth interval from 3.7 to 7.0 km below land surface. It has a transmissivity of 1.08×10-5 m2/s and a storage coefficient of 1.0×10-5. Computed pressure buildup along the length of the reservoir is compared with the spatial distribution of earthquake epicenters. The comparison shows that earthquakes are confined to that part of the reservoir where the pressure buildup exceeds 32 bars. This critical value is interpreted as the pressure buildup above which earthquakes occur. The migration of earthquake epicenters away from the injection well, a phenomenon noted by previous investigators, can be accounted for by the outward propagation of the critical pressure buildup. The analysis is extended to examining the effects of rapid flow in fractures opened by high injection pressure. The results show that the effect is confined to a small region within 1 km of the injection well. The existence of a critical pressure buildup above which earthquakes occur is completely consistent with the theory on the role of fluid pressure in fault movement as presented by Hubbert and Rubey.

Colorado↗

USGS2016

Late Holocene volcanism at Medicine Lake Volcano, northern California Cascades

Late Holocene volcanism at Medicine Lake volcano in the southern Cascades arc exhibited widespread and compositionally diverse magmatism ranging from basalt to rhyolite. Nine well-characterized eruptions have taken place at this very large rear-arc volcano since 5,200 years ago, an eruptive frequency greater than nearly all other Cascade volcanoes. The lavas are widely distributed, scattered over an area of ~300 km 2 across the >2,000-km 2 volcano. The eruptions are radiocarbon dated and the ages are also constrained by paleomagnetic data that provide strong evidence that the volcanic activity occurred in three distinct episodes at ~1 ka, ~3 ka, and ~5 ka. The ~1-ka final episode produced a variety of compositions including west- and north-flank mafic flows interspersed in time with fissure rhyolites erupted tangential to the volcano’s central caldera, including the youngest and most spectacular lava flow at the volcano, the ~950-yr-old compositionally zoned Glass Mountain flow. At ~3 ka, a north-flank basalt eruption was followed by an andesite eruption 27 km farther south that contains quenched basalt inclusions. The ~5-ka episode produced two caldera-focused dacitic eruptions. Quenched magmatic inclusions record evidence of intrusions that did not independently reach the surface. The inclusions are present in five andesitic, dacitic, and rhyolitic host lavas, and were erupted in each of the three episodes. Compositional and mineralogic evidence from mafic lavas and inclusions indicate that both tholeiitic (dry) and calcalkaline (wet) parental magmas were present. Petrologic evidence records the operation of complex, multi-stage processes including fractional crystallization, crustal assimilation, and magma mixing. Experimental evidence suggests that magmas were stored at 3 to 6 km depth prior to eruption, and that both wet and dry parental magmas were involved in generating the more silicic magmas. The broad distribution of eruptive events and the relative accessibility and good exposure of lavas, combined with physical and petrologic evidence for multiple and varied mafic inputs, has created an unusual opportunity to understand the workings of this large magmatic system. A combined total of more than 25 intrusive and extrusive events are indicated for late Holocene time. Plutonic inclusions, some with ages as young as Holocene, were also brought to the surface in five of the eruptions. All eruptions took place along northwest- to northeast-trending alignments of vents, reflecting the overall east-west extensional tectonic environment. The interaction of tectonism and volcanism is a dominant influence at this subduction-related volcano, located where the west edge of the extensional Basin and Range Province impinges on the Cascades arc. Ongoing subsidence focused at the central caldera has been documented along with geophysical evidence for a small magma body. This evidence, combined with the frequency of eruptive and intrusive activity in late Holocene time, an active geothermal system, and intermittent long-period seismic events indicate that the volcano is likely to erupt again.

California↗

USGS1978

The marine geology of the eastern Santa Barbara Channel, with particular emphasis on the ground water basins offshore from the Oxnard Plain, Southern California

Marine geophysical investigations provide new data concerning the stratigraphy, tectonic and sedimentary history, and the ground water geology of the southeastern Santa Barbara Channel region. The offshore stratigraphy identified in seismic reflection profiles includes a succession of Neogene to Quaternary strata. The middle Miocene Conejo volcanics form an acoustical basement and the overlying late Cenozoic sedimentary rocks attain a thickness greater than 2,500 m. These sedimentary deposits fill a structurally controlled, physiographic and depositional depression called the Ventura Basin. Structure consists generally of a gently folded, east-trending Tertiary synclinorium bordered on the north by a regional thrust fault and on the south by a steep asymmetrical anticlinal ridge. Most structures show evidence of north-south compression that occurred during early Pleistocene time. Three well-defined unconformities represent widespread erosion in late Miocene, early to middle Pleistocene, and late Pleistocene time. The boundaries of Miocene, Pliocene, and lower Pleistocene strata continue uninterrupted eastward along the southern part of Santa Barbara basin to Hueneme Canyon, where they turn northeast and can be traced to the coast near Port Hueneme. These limits probably represent the south edge of the Santa Barbara basin during Pliocene and Pleistocene time. Fresh water-bearing materials of the Oxnard plain are unconsolidated Quaternary sediment laid down on more consolidated Tertiary rocks. Offshore, the total fresh water-bearing materials distinguished in the seismic reflection profiles attain a thickness of about 356 m and have an areal extent of over 760 km2. Strata that contain the offshore continuation of the five major on-land aquifers (Grimes Canyon, Fox Canyon, Hueneme, Mugu, and Oxnard aquifers) are identified in the seismic reflection profiles. These strata make up the two offshore ground-water basins, the Mound and Oxnard plain ground-water basins, which are separated by the east-west trending Oak Ridge fault. Possible entrance areas for salt water intrusion into fresh water aquifers are found along the walls of the submarine canyons and along the northern slopes of Santa Barbara and Santa Monica basins. Hueneme and Mugu aquifers are probably exposed locally in all five submarine canyons of the Oxnard offshore area and may also crop out along the upper northern slope of Santa Monica basin. In all of these areas, salt water readily intrudes the aquifers. A salinity-temperature-depth study made in April, 1971, does not indicate any great dilution of surface ocean water by fresh water that could be 'leaking' from the exposed aquifers along the walls of Hueneme Canyon and the landward slope of Santa Barbara Channel. Earthquakes in the vicinity of the Oxnard plain suggest that the region is seismically active. Epicenters are widely dispersed over the region. No distinct trend or alignment of earthquake epicenters occurs near the trace of any of the faults, although many epicenters are scattered around the Oak Ridge zone of deformation in the northern part of the region. The largest magnitude earthquake recorded in the area was a magnitude 5.7 that occurred on February 21, 1973, offshore of Point Mugu, south of the Oxnard plain.

California↗

USGS1981

Surface geology of the Jeptha Knob cryptoexplosion structure, Shelby County, Kentucky

The Jeptha Knob crytoexplosion structure, described by Bucher in 1925, was remapped in 1973 as part of the U.S. Geological Survey and the Kentucky Geological Survey cooperative mapping program. The knob is in the western part of the Blue Grass region. Hilltops in the rolling farmland adjacent to the knob are underlain by the nearly flat-lying Grant Lake and Callaway Creek Limestones of middle Late Ordovician age, and the valleys are cut in interbedded limestone and shale of the Clays Ferry Formation of late Middle and early Late Ordovician age. Precambrian basement is estimated to be 4,000 ft below the surface. The mapped area is 50 miles west of the crest of the Cincinnati arch; the regional dip is westward 16 ft per mile. The 38th parallel lineament is 50 miles to the south. The structure, about 14,000 ft in diameter, consists of a central area 6,300 ft in diameter of uplifted Clays Ferry Formation surrounded by a belt of annular faults that are divided into segments by radial faults. The grass structure of the Clays Ferry Formation is that of a broad dame, but same evidence indicates that, in detail, the beds are complexly folded. The limestone of the Clays Ferry is brecciated and infiltrated by limonite. The brecciation is confined to single beds, and there is no mixing of fragments from different beds. A small plug of the Logana Member of the Lexington Limestone (Middle Ordovician) has been upfaulted at least 700 ft and emplaced within the Clays Ferry. The central uplift is separated by high-angle and, in places, reverse faults from the belt of annular faulting. The concentric faults in the zone of annular faults are extensional, and the general aspect is of collapse and inward movement. Lenses of breccia are present along many of the concentric faults, but not along the radial faults. At least same of the breccia was injected from below. The youngest beds involved in the faulting are in the Bardstown Member of the Drakes Formation of late Late Ordovician age. The faulted and brecciated beds are overlain by nearly horizontal dolomite and shale of Early and Middle Silurian age. The basal 5 ft of the oldest Silurian unit, the Brassfield Formation, contains calcarenite and calcirudite composed, in large part, of locally derived fragments from the Upper Ordovician formations. The Jeptha Knob structure was formed in latest Late Ordovician or earliest Early Silurian time. At the time of formation, the area was either very slightly above or very slightly below sea level; the sediments were already largely indurated. At the onset of Silurian deposition, the area of the central uplift was probably a broad shallow depression not more than about 15 ft deep, possibly surrounded by a rim of Upper Ordovician rocks or rock fragments. The origin of the Jeptha Knob structure cannot be determined from the available data. Shatter cones and coesite, considered by many to be definitive criteria far origin by impact, have not been found. On the other hand, geophysical studies indicate that there is no coincident uplift of the basement, and there is no certain relation of Jeptha Knob to any obvious structural trend.

Professional Paper↗

USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 1. Strong motions

We develop a methodology to perform finite fault source inversions from strong motion data using Green's functions (GFs) calculated for a three-dimensional (3-D) velocity structure. The 3-D GFs are calculated numerically by inserting body forces at each of the strong motion sites and then recording the resulting strains along the target fault surface. Using reciprocity, these GFs can be recombined to represent the ground motion at each site for any (heterogeneous) slip distribution on the fault. The reciprocal formulation significantly reduces the required number of 3-D finite difference computations to at most 3NS, where NS is the number of strong motion sites used in the inversion. Using controlled numerical resolution tests, we have examined the relative importance of accurate GFs for finite fault source inversions which rely on near-source ground motions. These experiments use both 1-D and 3-D GFs in inversions for hypothetical rupture models in order (1) to analyze the ability of the 3-D methodology to resolve trade-offs between complex source phenomena and 3-D path effects, (2) to address the sensitivity of the inversion results to uncertainties in the 3-D velocity structure, and (3) to test the adequacy of the 1-D GF method when propagation effects are known to be three-dimensional. We find that given "data" from a prescribed 3-D Earth structure, the use of well-calibrated 3-D GFs in the inversion provides very good resolution of the assumed slip distribution, thus adequately separating source and 3-D propagation effects. In contrast, using a set of inexact 3-D GFs or a set of hybrid 1-D GFs allows only partial recovery of the slip distribution. These findings suggest that in regions of complex geology the use of well-calibrated 3-D GFs has the potential for increased resolution of the rupture process relative to 1-D GFs. However, realizing this full potential requires that the 3-D velocity model and associated GFs should be carefully validated against the true 3-D Earth structure before performing the inverse problem with actual data.

Journal of Geophysical Research B: Solid Earth↗

USGS2024

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York↗

USGS2010

Outcrops, fossils, geophysical logs, and tectonic interpretations of the Upper Cretaceous Frontier Formation and contiguous strata in the Bighorn Basin, Wyoming and Montana

In the Bighorn Basin of north-central Wyoming and south-central Montana, the Frontier Formation of early Late Cretaceous age consists of siliciclastic, bentonitic, and carbonaceous beds that were deposited in marine, brackish-water, and continental environments. Most lithologic units are laterally discontinuous. The Frontier Formation conformably overlies the Mowry Shale and is conformably overlain by the Cody Shale. Molluscan fossils collected from outcrops of these formations and listed in this report are mainly of marine origin and of Cenomanian, Turonian, and Coniacian ages. The lower and thicker part of the Frontier in the Bighorn Basin is of Cenomanian age and laterally equivalent to the Belle Fourche Member of the Frontier in central Wyoming. Near the west edge of the basin, these basal strata are disconformably overlain by middle Turonian beds that are the age equivalent of the Emigrant Gap Member of the Frontier in central Wyoming. The middle Turonian beds are disconformably overlain by lower Coniacian strata. Cenomanian strata along the south and east margins of the basin are disconformably overlain by upper Turonian beds in the upper part of the Frontier, as well as in the lower part of the Cody; these are, in turn, conformably overlain by lower Coniacian strata. Thicknesses and ages of Cenomanian strata in the Bighorn Basin and adjoining regions are evidence of regional differential erosion and the presence of an uplift during the early Turonian centered in northwestern Wyoming, west of the basin, probably associated with a eustatic event. The truncated Cenomanian strata were buried by lower middle Turonian beds during a marine transgression and possibly during regional subsidence and a eustatic rise. An uplift in the late middle Turonian, centered in north-central Wyoming and possibly associated with a eustatic fall, caused the erosion of lower middle Turonian beds in southern and eastern areas of the basin as well as in an adjoining region of north-central Wyoming. Similarly, in east-central Wyoming and an adjacent area to the south, Cenomanian strata are disconformably overlain by upper middle and lower upper Turonian strata that probably reflect uplift and erosion in that region during the interim period of middle Turonian time. During later subsidence and a marine transgression, upper Turonian deposits buried Cenomanian beds in areas along the south and east margins of the Bighorn Basin and buried lower middle Turonian beds in much of northern Wyoming. Upper Turonian and lower Coniacian strata are apparently conformable in eastern and southern areas of the basin as well as near Riverton, Kaycee, and Casper in central Wyoming. Upper Turonian strata are absent on the west flank of the Bighorn Basin and in outcrops west of the basin, where middle Turonian beds are disconformably overlain by lower Coniacian beds . The conformable upper Turonian and lower Coniacian beds apparently transgressed an eroded middle Turonian surface in the region, but only Coniacian strata overlie middle Turonian beds on the west side of the basin and areas farther west. Coniacian strata onlap the truncated lower middle Turonian surface west of the basin, indicating a region that had higher elevation possibly resulting from tectonic uplift. In east-central Wyoming and an adjoining region to the south, upper middle Turonian and lower upper Turonian strata are disconformably overlain by lower and middle Coniacian beds. That region apparently was uplifted and eroded during the latest Turonian.

Montana, Wyoming↗

USGS1995

An earthquake instability model based on faults containing high fluid-pressure compartments

It has been proposed that large strike-slip faults such as the San Andreas contain water in seal-bounded compartments. Arguments based on heat flow and stress orientation suggest that in most of the compartments, the water pressure is so high that the average shear strength of the fault is less than 20 MPa. We propose a variation of this basic model in which most of the shear stress on the fault is supported by a small number of compartments where the pore pressure is relatively low. As a result, the fault gouge in these compartments is compacted and lithified and has a high undisturbed strength. When one of these locked regions fails, the system made up of the neighboring high and low pressure compartments can become unstable. Material in the high fluid pressure compartments is initially underconsolidated since the low effective confining pressure has retarded compaction. As these compartments are deformed, fluid pressure remains nearly unchanged so that they offer little resistance to shear. The low pore pressure compartments, however, are overconsolidated and dilate as they are sheared. Decompression of the pore fluid in these compartments lowers fluid pressure, increasing effective normal stress and shear strength. While this effect tends to stabilize the fault, it can be shown that this dilatancy hardening can be more than offset by displacement weakening of the fault (i.e., the drop from peak to residual strength). If the surrounding rock mass is sufficiently compliant to produce an instability, slip will propagate along the fault until the shear fracture runs into a low-stress region. Frictional heating and the accompanying increase in fluid pressure that are suggested to occur during shearing of the fault zone will act as additional destabilizers. However, significant heating occurs only after a finite amount of slip and therefore is more likely to contribute to the energetics of rupture propagation than to the initiation of the instability. We present results of a one-dimensional dynamic Burridge-Knopoff-type model to demonstrate various aspects of the fluid-assisted fault instability described above. In the numerical model, the fault is represented by a series of blocks and springs, with fault rheology expressed by static and dynamic friction. In addition, the fault surface of each block has associated with it pore pressure, porosity and permeability. All of these variables are allowed to evolve with time, resulting in a wide range of phenomena related to fluid diffusion, dilatancy, compaction and heating. These phenomena include creep events, diffusion-controlled precursors, triggered earthquakes, foreshocks, aftershocks, and multiple earthquakes. While the simulations have limitations inherent to 1-D fault models, they demonstrate that the fluid compartment model can, in principle, provide the rich assortment of phenomena that have been associated with earthquakes. ?? 1995 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

USGS2012

Basin and range crustal and upper mantle structure, northwest to central Nevada

We present an interpretation of the crustal and uppermost mantle structure of the Basin and Range of northwestern Nevada based on seismic refraction/wide‐angle reflection, near‐vertical reflection, and gravity data. In comparison to most previous estimates, we find that the crust is somewhat thicker (32–36 km versus 22–30 km), and the uppermost mantle velocity is somewhat higher (8.0 km/s versus 7.3–7.9 km/s). Along our transects, the crust is thinnest (32 km) in the Carson Sink‐Buena Vista Valley region and increases by 2–4 km to the west and east, respectively. There is considerable complexity throughout the crust where velocities range from of 2.5 km/s at the surface to 7.4 km/s in the lowermost crust. Variations in velocity and structure of the upper crustal layers reveal apparent basement velocity depressions (areas of lower velocities extending up to 10 km in depth) that underlie some surface ranges as well as the basins. The middle crust rises from about 20 km beneath central Nevada to within 12 km of the surface beneath the area of thinnest crust and is characterized by a modest (∼0.1 km/s) change in velocity and low‐velocity gradients. These midcrustal layers mark the onset of high crustal reflectivity and the apparent limiting depth to which Basin and Range faults can be traced in near‐vertical reflection profiles, suggesting that these midcrustal layers represent the transition between the brittle and ductile zones of the crust. The lower crust is more structurally complex, with layers thickening and thinning in a systematic manner with the upper crustal layers; generally, where there are velocity depressions in the upper crust, the lower crust is thickest and shallowest. The geometry of these lower crustal layers (derived from refraction modeling) coincides with changes in the crustal reflectivity, determined from the Consortium of Continental Reflection Profiling reflection data. The lower crustal layer is unusually high in velocity (7.4 km/s) and is likely the layer identified as mantle in some previous studies. We do not identify the 7.4 km/s layer as mantle because (1) there is an underlying layer with a velocity (8.0 km/s) that is more consistent with the worldwide average velocity for the upper mantle, and (2) the 7.4 km/s layer does not correspond to the “reflection” Moho. Gravity modeling and comparison to existing seismic models show a general consensus in many aspects with respect to crustal structure. This new model forms the basis for speculation on some of the processes associated with rifting of the Basin and Range Province. One such process, lithospheric magmatism, is inferred from the strong attenuation of transmitted seismic waves, which occurs at the same interface at which high‐amplitude, bright spot reflections originate. Unlike previous models, the overall structure and velocity of the crust and uppermost mantle of our new model are similar to other regions worldwide which have undergone high degrees of extension.

Nevada↗

USGS2001

Mars Global Surveyor Thermal Emission Spectrometer experiment: Investigation description and surface science results

The Thermal Emission Spectrometer (TES) investigation on Mars Global Surveyor (MGS) is aimed at determining (1) the composition of surface minerals, rocks, and ices; (2) the temperature and dynamics of the atmosphere; (3) the properties of the atmospheric aerosols and clouds; (4) the nature of the polar regions; and (5) the thermophysical properties of the surface materials. These objectives are met using an infrared (5.8- to 50-μm) interferometric spectrometer, along with broadband thermal (5.1- to 150-μm) and visible/near-IR (0.3- to 2.9-μm) radiometers. The MGS TES instrument weighs 14.47 kg, consumes 10.6 W when operating, and is 23.6×35.5×40.0 cm in size. The TES data are calibrated to a 1-σ precision of 2.5 −6 ×10 −8 W cm −2 sr −1 /cm −1 , 1.6×10 −6 W cm −2 sr −1 , and ∼0.5 K in the spectrometer, visible/near-IR bolometer, and IR bolometer, respectively. These instrument subsections are calibrated to an absolute accuracy of ∼4×10 −8 W cm −2 sr −1 /cm −1 (0.5 K at 280 K), 1–2%, and ∼1–2 K, respectively. Global mapping of surface mineralogy at a spatial resolution of 3 km has shown the following: (1) The mineralogic composition of dark regions varies from basaltic, primarily plagioclase feldspar and clinopyroxene, in the ancient, southern highlands to andesitic, dominated by plagioclase feldspar and volcanic glass, in the younger northern plains. (2) Aqueous mineralization has produced gray, crystalline hematite in limited regions under ambient or hydrothermal conditions; these deposits are interpreted to be in-place sedimentary rock formations and indicate that liquid water was stable near the surface for a long period of time. (3) There is no evidence for large-scale (tens of kilometers) occurrences of moderate-grained (>50-μm) carbonates exposed at the surface at a detection limit of ∼10%. (4) Unweathered volcanic minerals dominate the spectral properties of dark regions, and weathering products, such as clays, have not been observed anywhere above a detection limit of ∼10%; this lack of evidence for chemical weathering indicates a geologic history dominated by a cold, dry climate in which mechanical, rather than chemical, weathering was the significant form of erosion and sediment production. (5) There is no conclusive evidence for sulfate minerals at a detection limit of ∼15%. The polar region has been studied with the following major conclusions: (1) Condensed CO 2 has three distinct end-members, from fine-grained crystals to slab ice. (2) The growth and retreat of the polar caps observed by MGS is virtually the same as observed by Viking 12 Martian years ago. (3) Unique regions have been identified that appear to differ primarily in the grain size of CO 2 ; one south polar region appears to remain as black slab CO 2 ice throughout its sublimation. (4) Regional atmospheric dust is common in localized and regional dust storms around the margin and interior of the southern cap. Analysis of the thermophysical properties of the surface shows that (1) the spatial pattern of albedo has changed since Viking observations, (2) a unique cluster of surface materials with intermediate inertia and albedo occurs that is distinct from the previously identified low-inertia/bright and high-inertia/dark surfaces, and (3) localized patches of high-inertia material have been found in topographic lows and may have been formed by a unique set of aeolian, fluvial, or erosional processes or may be exposed bedrock.

Journal of Geophysical Research E: Planets↗

USGS2012

Strain accumulation rates in the western United States between 1970 and 1978

The rate of dilatation and the rate and direction of shear have been determined from trilateration data for 23 Geodolite networks in the western United States. Sixteen nets are located along the San Andreas fault system between Point Reyes, California, and the United States‐Mexico border. Other locations are across the Garlock fault in California; across Puget Sound near Seattle, Washington; near Hanford in eastern Washington; near Hebgen Lake in Montana; across the Wasatch fault at Ogden, Utah; across the Rio Grande rift at Socorro, New Mexico; and Dixie Valley in Nevada; and at the northern end of Owens Valley on the California‐Nevada border. Implicit in the treatment are the assumptions that the strain was accumulating at a constant rate over the time period (within the interval 1970–1978) and over the local area (usually about 50‐km diameter) covered by the surveys. Of the nets located away from the San Andreas fault, only Ogden and Hebgen show significant strain accumulation. At Ogden the deformation is principally an east‐west compression of 0.23±0.05 μstrain/yr and at Hebgen Lake a northeast‐southwest extension of 0.17±0.03 μstrain/yr. Along the San Andreas fault system the rate of shear is 0.2 to 0.4μ/yr. The direction of shear agrees very well with the surface strike of nearby faults. This agreement is maintained even in regions like the ‘big bend,’ where both the fault strike and the observed shear direction are more westerly than they are elsewhere. Shear strain in northern California appears to be concentrated more closely on the faults, whereas in southern California the strain is a broader, smoother feature. In the San Francisco Bay area the strain data indicate slip at depth on both the San Andreas and the Calaveras faults. In addition to the observed shear the nets in California indicate a negative dilatation (areal decrease) of about 0.2 μstrain/yr. This dilatation is unexplained, but the following sources appear unlikely: (1) systematic survey error, (2) an association with the southern California uplift, (3) an association with the big bend in the San Andreas fault in Southern California, or (4) the result of the superposition of a uniaxial strain on the Pacific‐North American plate boundary shear.

Journal of Geophysical Research B: Solid Earth↗

USGS2026

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Alaska↗

USGS2013

Intermediate‐depth intraslab earthquakes and arc volcanism as physical expressions of crustal and uppermost mantle metamorphism in subducting slabs

We elaborate on the well-known spatial association between axc volcanoes and Wadati Benioff zones and explore in detail their genetic relationships as dual physical expressions of slab metamorphism of the oceanic crust and uppermost mantle. At hypocentral depths less than 200 km intra slab Wadati-Benioff earthquakes tend to occur near the top surfaces of slabs. Subduction of very young lithosphere (age < 15-25 Ma) with high heat flow (> 75 mW/m 2) produces mainly shallow earthquakes and spaxse or absent arc volcanism. Subduction of older crust with normal heat flow (50-65 mW/m 2) produces markedly deeper intraslab earthquakes and generally normal volcanic vigor. Seismological observations show that the low-seismic-velocity gabbroic mineralogy of the crust may persist to depths of as much as 150 km in old, cold lithosphere but only to depths of 50-60 km in young, warm slabs. Metamorphic processes in the crust and shallow upper mantle of subducting slabs and the reactivation of faults originally created at shallow depths in the ocean basins probably control the occurrence of intraslab earthquakes to depths of as much as 350 kin. A conceptual model for this metamorphism incorporates the likely effects of water liberated by dehydration. Such dehydration facilitates both brittle faulting by fault reactivation and promotes the kinetics of the transformation of the anhydrous gabbro component of the crust to eclogite. Finite-element modelling shows that densification to eclogite is expected to produce extensional stresses in transformed crust and a smaller compression in the underlying mantle. This model helps explain why most intermediate-depth intraslab earthquakes occur just below the top surfaces of slabs and why many have focal mechanisms indicating down-dip extension. Young, warm slabs have mostly shallow intraslab earthquakes and sparse axc volcanoes because dehydration and eclogite formation largely cease before such slabs axe in contact with asthenosphere. These processes are evidently delayed by kinetic hindrance and the high pressure stability of hydrous phazes at low temperatures in older, colder subducting crust, and thus earthquake activity and asthenospheric-wedge melting tend to be focused at depths of 100-170 kin. Anomalous behavior correlated with the subduction of island-and seamount chains appears to be associated with anomalou shallow intraplate faulting and with perturbations of slab metamorphism by these chains.

Geophysical Monograph↗