USGS Science⌕ Search

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.

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

At least 595 records · Page 33Linked to original sources
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain 2. Numerical predictions and experimental tests

Numerical solutions of the equations describing flow of variably fluidized Coulomb mixtures predict key features of dry granular avalanches and water-saturated debris flows measured in physical experiments. These features include time-dependent speeds, depths, and widths of flows as well as the geometry of resulting deposits. Threedimensional (3-D) boundary surfaces strongly influence flow dynamics because transverse shearing and cross-stream momentum transport occur where topography obstructs or redirects motion. Consequent energy dissipation can cause local deceleration and deposition, even on steep slopes. Velocities of surge fronts and other discontinuities that develop as flows cross 3-D terrain are predicted accurately by using a Riemann solution algorithm. The algorithm employs a gravity wave speed that accounts for different intensities of lateral stress transfer in regions of extending and compressing flow and in regions with different degrees of fluidization. Field observations and experiments indicate that flows in which fluid plays a significant role typically have high-friction margins with weaker interiors partly fluidized by pore pressure. Interaction of the strong perimeter and weak interior produces relatively steep-sided, flat-topped deposits. To simulate these effects, we compute pore pressure distributions using an advection-diffusion model with enhanced diffusivity near flow margins. Although challenges remain in evaluating pore pressure distributions in diverse geophysical flows, Riemann solutions of the depthaveraged 3-D Coulomb mixture equations provide a powerful tool for interpreting and predicting flow behavior. They provide a means of modeling debris flows, rock avalanches, pyroclastic flows, and related phenomena without invoking and calibrating Theological parameters that have questionable physical significance.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Quantifying canal leakage rates using a mass-balance approach and heat-based hydraulic conductivity estimates in selected irrigation canals, western Nebraska, 2007 through 2009

The water supply in areas of the North Platte River Basin in the Nebraska Panhandle has been designated as fully appropriated or overappropriated by the Nebraska Department of Natural Resources (NDNR). Enacted legislation (Legislative Bill 962) requires the North Platte Natural Resources District (NPNRD) and the NDNR to develop an Integrated Management Plan (IMP) to balance groundwater and surface-water supply and demand in the NPNRD. A clear understanding of the groundwater and surface-water systems is critical for the development of a successful IMP. The primary source of groundwater recharge in parts of the NPNRD is from irrigation canal leakage. Because canal leakage constitutes a large part of the hydrologic budget, spatially distributing canal leakage to the groundwater system is important to any management strategy. Surface geophysical data collected along selected reaches of irrigation canals has allowed for the spatial distribution of leakage on a relative basis; however, the actual magnitude of leakage remains poorly defined. To address this need, the U.S. Geological Survey, in cooperation with the NPNRD, established streamflow-gaging stations at upstream and downstream ends from two selected canal reaches to allow a mass-balance approach to be used to calculate daily leakage rates. Water-level and sediment temperature data were collected and simulated at three temperature monitoring sites to allow the use of heat as a tracer to estimate the hydraulic conductivity of canal bed sediment. Canal-leakage rates were estimated by applying Darcy's Law to modeled vertical hydraulic conductivity and either the estimated or measured hydraulic gradient. This approach will improve the understanding of the spatial and temporal variability of canal leakage in varying geologic settings identified in capacitively coupled resistivity surveys. The high-leakage potential study reach of the Tri-State Canal had two streamflow-gaging stations and two temperature monitoring sites along its length. Calculated leakage rates from the mass-balance approach varied from year to year and were generally dependent on local climatic conditions, and the timing and magnitude of the initial seasonal diversion into the Tri-State Canal. Leakage rates ranged from 0.98 meter per day (m/d) on June 22, 2007, to about to 0 m/d during July 2009. Drier conditions generally resulted in higher leakage rates because of reduced flow from Spottedtail Creek, lower groundwater levels near Spottedtail Creek, and no unmeasured flow entering the reach. Of the three years studied (2007-09), 2007 was the driest, and therefore had the highest canal leakage rates. The moderately low leakage potential study reach of Interstate Canal had two streamflow-gaging stations and one temperature monitoring site along its length. Excluding the leakage calculations from early May 2007, leakage rates ranged from 0.08 to 0.7 m/d. Less variability in leakage from year to year indicates that climatic conditions may have less of an effect for Interstate Canal compared to Tri-State Canal. This may be because Interstate Canal was cut into the northern edge of the North Platte alluvial valley and consequently the canal bed is well above the local groundwater table resulting in a constant (1 meter per meter [m/m]) hydraulic gradient. Interstate Canal also does not receive any captured flow that can vary substantially year to year. Two temperature monitoring sites were installed within the high-leakage potential reach of Tri-State Canal. Site TCTEMP1 was established in 2007 where the water table was well below the canal bed surface. The vertical hydraulic conductivity of the poorly sorted sand and gravel beneath site TCTEMP1 was estimated using a calibrated one-dimensional VS2DH model. Using a trial-and-error approach, the best-fit vertical hydraulic conductivity for the site TCTEMP1 model domain was 1.1 m/d. Site TCTEMP2 was established at the mouth of Spottedtail Creek where a shallow

Scientific Investigations Report↗
USGS2005

Modeling effects of bank friction and woody bank vegetation on channel flow and boundary shear stress in the Rio Puerco, New Mexico

[1] We have applied a physically based model for steady, horizontally uniform flow to calculate reach-averaged velocity and boundary shear-stress distributions in a natural stream with woody vegetation on the channel banks. The model calculates explicitly the form drag on woody plant stems and includes the effects of vegetation on the boundary shear stress, velocity, and turbulence fields. Average channel shapes, bed gradients, and shrub characteristics were determined for four long, internally similar reaches covering an 81 -km segment of the lower Rio Puerco, in north-central New Mexico. Hydraulic geometries for each of three flow events were determined from reach-average channel shapes, slopes, and measured silt-line heights above the bed. Model results show that friction on the lateral boundaries reduced the boundary shear stress in the center of the Rio Puerco channel by as much as 20%. In reaches with moderate to dense bank shrubs, perimeter-averaged boundary shear stresses were reduced by almost 40% during near bankfull flows. Model-calculated discharges for all three flow events indicate there was a loss of about 40% of the upstream discharge through the 81-km river segment. Infiltration rates estimated from discharge losses during the quasi-steady peak flows and areas of the bed and banks suggest rates of flow loss into the silty sand of the upper banks are about 8 times greater than rates of loss into the clayey silt-covered bed and lower banks.

Journal of Geophysical Research F: Earth Surface↗
USGS1996

The July 12, 1993, Hokkaido-Nansei-Oki, Japan, earthquake: Coseismic slip pattern from strong-motion and teleseismic recordings

We employ a finite fault inversion scheme to infer the distribution of coseismic slip for the July 12, 1993 Hokkaido-Nansei-Oki earthquake using strong ground motions recorded by the Japan Meteorological Agency within 400 km of the epicenter and vertical P waveforms recorded by the Global Digital Seismograph Network at teleseismic distances. The assumed fault geometry is based on the location of the aftershock zone and comprises two fault segments with different orientations: a northern segment striking at N20°E with a 30° dip to the west and a southern segment with an N20°W strike. For the southern segment we use both westerly and easterly dip directions to test thrust orientations previously proposed for this portion of the fault. The variance reduction is greater using a shallow west dipping segment, suggesting that the direction of dip did not change as the rupture propagated south from the hypocenter. This indicates that the earthquake resulted from the shallow underthrusting of Hokkaido beneath the Sea of Japan. Static vertical movements predicted by the corresponding distribution of fault slip are consistent with the general pattern of surface deformation observed following the earthquake. Fault rupture in the northern segment accounts for about 60% of the total P wave seismic moment of 3.4 × 10 20 N m and includes a large circular slip zone (4-m peak) near the earthquake hypocenter at depths between 10 and 25 km. Slip in the southern segment is also predominantly shallower than 25 km, but the maximum coseismic displacements (2.0–2.5 m) are observed at a depth of about 5 km. This significant shallow slip in the southern portion of the rupture zone may have been responsible for the large tsunami that devastated the small offshore island of Okushiri. Localized shallow faulting near the island, however, may require a steep westerly dip to reconcile the measured values of ground subsidence.

Journal of Geophysical Research B: Solid Earth↗
USGS1999

Geologic map of the Galaxias quadrangle (MTM 35217) of Mars

The Galaxias region (MTM 35217) is one of a series of 1:500,000-scale science study areas on Mars sponsored by NASA's Planetary Geology and Geophysics Program. Situated near the northern limit of lava flows associated with Elysium Mons, this region includes a mixture of volcanic and nonvolcanic terrains. The region is also of interest for the fluvial systems that originate along the distal margins of the Elysium lava flows. Resolution of Viking Orbiter images used to prepare the base map ranges from 40 to 160 m/pixel. High-resolution frames (40 to 80 m/pixel) are found in the southeastern part of the map area and along the north edge of the quadrangle, but over half the quadrangle is included in medium-resolution frames (150 m/pixel). Two 8 m/pixel, very high resolution scenes are available (see fig. 1). Interpretation is complicated by variable resolution and sun angles that vary from east to west illumination on different images. Mapping methods and principles are adapted from those developed for lunar photogeologic mapping by Shoemaker and Hackman (1962), refined by Wilhelms (1972), and successfully applied by many workers to a variety of planetary surfaces. Mapping units are distinguished by topography and texture and are ranked by relative age on the basis of superposition and transection relations. Material units are assigned to time-stratigraphic systems defined by Scott and Carr (1978) and Tanaka (1986). This area is included within earlier maps that used Mariner 9 images at 1:5,000,000 scale (Elston, 1979) and globally at 1:25,000,000 scale (Scott and Carr, 1978). Regional maps based on the much higher resolutions of Viking Orbiter allowed more detailed discrimination of materials by Greeley and Guest (1987) at 1:15,000,000 scale and Tanaka and others (1992) at 1:5,000,000 scale. Some map units on this 1:500,000-scale map correspond to, or are partially equivalent to, units on the larger scale maps of Greeley and Guest (1987) and Tanaka and others (1992). Established terminology is used where feasible, but the scale of this map requires that some new units be introduced and that some previous terminology be redefined. Photogeologic methods are limited; therefore, more than one geologic explanation is given for some material units that do not readily lend themselves to an unequivocal interpretation.

IMAP↗
USGS2021

Hydrogeologic framework, water levels, and selected contaminant concentrations at Valmont TCE Superfund Site, Luzerne County, Pennsylvania, 2020

The Valmont TCE Superfund Site, Luzerne County, Pennsylvania is underlain by fractured and folded sandstones and shales of the Pottsville and Mauch Chunk Formations, which form a fractured-rock aquifer recharged locally by precipitation. Industrial activities at the former Chromatex Plant resulted in trichloroethene (TCE) contamination of groundwater at and near the facility, which was identified in 1987 and led to listing as a Superfund site by the U.S. Environmental Protection Agency (EPA) in 1989. To address the problem of TCE concentrations in nearby residential wells that exceed the maximum contaminant level (MCL) of 5 micrograms per liter (μg/L), alternate water supplies were provided. A 2015 review of initial characterization and subsequent remediation by the EPA identified the need for an updated understanding of the complex hydrogeology and the conceptual site model. Additional contaminants present in groundwater at the site include some other volatile organic compounds (VOCs) and per- and polyfluoroalkyl substances (PFAS), predominantly consisting of perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS) present in concentrations that exceeded the EPA Health Advisory (HA) level of 5 nanograms per liter (ng/L) for combined PFOA and PFOS. In response to a request from the EPA in 2019, the U.S. Geological Survey (USGS) prepared cross sections and maps to provide more information about the hydrogeologic framework at and near the site and assist in improving the conceptual site model using water level and contaminant data collected by the EPA in 2020. The cross sections present lithologic correlations from available geophysical logs collected in wells from 2002 to 2014; they show alternating intervals of relatively elevated and reduced natural gamma activity that correspond to changes in lithology, with water-bearing zones and well screens commonly located at lithologic contacts, sometimes near thin coal seams. Water-bearing zones commonly are associated with fractures at or near lithologic contacts but also may be associated with fractures at or near apparent faulting. Recent (March 2020) water-level data shown on cross sections and maps indicate large downward vertical gradients and apparent radial gradients laterally to the northeast, northwest, and southwest that generally following topography. Recent (February to March 2020) data for TCE groundwater concentration shown on cross sections and maps indicate the highest TCE concentrations (greater than 3,000 μg/L and as much as 75,000 μg/L) and combined PFOA and PFOS concentrations (greater than 1,000 ng/L and up to at least 2,350 ng/L) are from shallow (less than 60 feet [ft] below land surface [bls]) and intermediate depth (60 to 100 ft bls) wells near the center of the former Chromatex Plant. TCE and PFAS (as combined PFOA and PFOS) contamination is present at greater depths, as much as 304 ft bls, as evidenced by samples collected from one well (a reconstructed former production well) near the plant, that contained concentrations of about 240 μg/L and 508 ng/L, respectively. The 2020 data also indicate that TCE and PFAS concentrations which exceed drinking-water MCL or HA levels are present in groundwater depths of less than 200 ft in an area that extends predominantly in a northeast direction from the former Chromatex Plant, and is apparently influenced by hydraulic gradients, lithology, and geologic structure.

Pennsylvania↗
USGS1988

Volumetric strain in relation to particle displacements for body and surface waves in a general viscoelastic half-space

Dilatational earth strain, associated with the radiation fields for several hundred local, regional, and teleseismic earthquakes, has been recorded over an extended bandwidth and dynamic range at four borehole sites near the San Andreas fault, CA. The general theory of linear viscoelasticity is applied to account for anelasticity of the near-surface materials and to provide a mathematical basis for interpretation of seismic radiation fields as detected simultaneously by co-located volumetric strain meters and seismometers. The general theory is applied to describe volumetric strain and displacement for general (homogeneous or inhomogeneous) P and S waves in an anelastic whole space. Solutions to the free-surface reflection problems for incident general P and S-I waves are used to evaluate the effect of the free surface on observations from co-located sensors. Corresponding expressions are derived for a Rayleigh-type surface wave on a linear viscoelastic half-space. The theory predicts a number of anelastic wave field characteristics that can be inferred from observation of volumetric strains and displacement fields as detected by co-located sensors that cannot be inferred from either sensor alone. -from Author

Geophysical Journal↗
USGS2021

The biophysical role of water and ice within permafrost nearing collapse: Insights from novel geophysical observations

The impact of permafrost thaw on hydrologic, thermal, and biotic processes remains uncertain, in part due to limitations in subsurface measurement capabilities. To better understand subsurface processes in thermokarst environments, we collocated geophysical and biogeochemical instruments along a thaw gradient between forested permafrost and collapse-scar bogs at the Alaska Peatland Experiment (APEX) site near Fairbanks, Alaska. Ambient seismic noise monitoring provided continuous high-temporal resolution measurements of water and ice saturation changes. Maps of seismic velocity change identified areas of large summertime velocity reductions nearest the youngest bog, indicating potential thaw and expansion at the bog margin. These results corresponded well with complementary borehole nuclear magnetic resonance measurements of unfrozen water content with depth, which showed permafrost soils nearest the bog edges contained the largest amount of unfrozen water along the study transect, up to 25% by volume. In situ measurements of methane within permafrost soils revealed high concentrations at these bog-edge locations, up to 30% soil gas. Supra-permafrost talik zones were observed at the bog margins, indicating talik formation and perennial liquid water may drive lateral bog expansion and enhanced permafrost carbon losses preceding thaw. Comparison of seismic monitoring with wintertime surface carbon dioxide fluxes revealed differential responses depending on time and proximity to the bogs, capturing the controlling influence of subsurface water and ice on microbial activity and surficial emissions. This study demonstrates a multidisciplinary approach for gaining new understanding of how subsurface physical properties influence greenhouse gas production, emissions, and thermokarst development.

Alaska↗
USGS2001

Mineralogy of the last lunar basalts: Results from Clementine

The last major phase of lunar volcanism produced extensive high-titanium mare deposits on the western nearside which remain unsampled by landing missions. The visible and near-infrared reflectance properties of these basalts are examined using Clementine multispectral images to better constrain their mineralogy. A much stronger 1 μm ferrous absorption was observed for the western high-titanium basalts than within earlier maria, suggesting that these last major mare eruptions also may have been the most iron-rich. These western basalts also have a distinctly long-wavelength, 1 μm ferrous absorption which was found to be similar for both surface soils and materials excavated from depth, supporting the interpretation of abundant olivine within these deposits. Spectral variation along flows within the Imbrium basin also suggests variations in ilmenite content along previously mapped lava flows as well as increasing olivine content within subsequent eruptions.

Journal of Geophysical Research E: Planets↗
USGS1981

Significant results from using earth observation satellites for mineral and energy resource exploration

A large number of Earth-observation satellites orbit our world several times each day, providing new information about the land and sea surfaces and the overlying thin layer of atmosphere that makes our planet unique. Meteorological satellites have had the longest history of experimental use and most are now considered operational. The geologic information collected by the Landsat, Polar Orbiting Geophysical Observatory (POGO), Magsat, Heat Capacity Mapping Mission (HCMM) and Seasat land and ocean observation systems is being thoroughly tested, and some of these systems are now approaching operational use. Landsat multispectral images provide views of large areas of the Earth under uniform lighting conditions and can be obtained at a variety of scales and formats. Not only do the Landsat data provide highly useful images showing surficial materials and structures such as folds and faults, but also measurements and computer-derived ratios of the brightness of different rock types, alteration zones, and mineral associations. These data have led to the finding of a variety of new ore deposits. In addition, the combination of Landsat digital data and aeromagnetic data has extended the use of Landsat as an exploration tool which can be used to readily relate surface features to subsurface anomalies. Magsat data, now being collected, are helping refine information on major crustal anomalies that were first recognized during the analysis of POGO data. The more nearly circular orbit, lower altitude, and increased sophistication of its vector magnetometer enable Magsat to provide more precise information than POGO. Information of this type is required to develop crustal models. Although Magsat is designed to operate for only 4–8 months, the number of orbits that it should be able to make will be sufficient to accomplish its mission and to record a major magnetic storm expected in 1980. HCMM is a two-band visible to near-IR (0.55–1.1 μm) and thermal infrared (10.2–12.5 μm) system designed to measure reflected solar energy, determine the heat capacity of rocks and to monitor soil moisture, thermal effluents, plant canopy temperatures and snow cover. Launched in April 1978, it is in sun-synchronous, circular orbit at an altitude of 620 km. It is a relatively low-resolution system with an instantaneous field of view (IFOV) of 500–600 m and a swath width of 716 km. However, the system is designed to detect objects in the range of 260°–340° K with a sensitivity (NEδT) of 0.4°K at 280°. Recording the thermal radiation of urban heat islands and high thermal inertia of quartzite strata in the Appalachian region are two examples of its land applications. Launched in June 1978, Seasat operated for only 100 days, but successfully acquired much information over both sea and land. The collection of synthetic aperture radar (SAR) imagery and radar altimetry was particularly important to geologists. Although there are difficulties in processing and distributing these data in a timely manner, initial evaluations indicate that the radar imagery supplements Landsat data by increasing the spectral range and offering a different look angle. The radar altimeter provides accurate profiles over narrow strips of land (1 km wide) and has demonstrated usefulness in measuring icecap surfaces (Greenland, Iceland, and Antarctica). The Salar of Uyuni in southern Bolivia served as a calibration site for the altimeter and has enabled investigators to develop a land-based smoothing algorithm that is believed to increase the accuracy of the system to 10 cm. Data from the altimeter are currently being used to measure subsidence resulting from ground water withdrawal in the Phoenix-Tucson area.

Advances in Space Research↗
USGS2011

Collection, processing, and interpretation of ground-penetrating radar data to determine sediment thickness at selected locations in Deep Creek Lake, Garrett County, Maryland, 2007

The U.S. Geological Survey collected geophysical data in Deep Creek Lake in Garrett County, Maryland, between September 17 through October 4, 2007 to assist the Maryland Department of Natural Resources to better manage resources of the Lake. The objectives of the geophysical surveys were to provide estimates of sediment thickness in shallow areas around the Lake and to test the usefulness of three geophysical methods in this setting. Ground-penetrating radar (GPR), continuous seismic-reflection profiling (CSP), and continuous resistivity profiling (CRP) were attempted. Nearly 90 miles of GPR radar data and over 70 miles of CSP data were collected throughout the study area. During field deployment and testing, CRP was determined not to be practical and was not used on a large scale. Sediment accumulation generally could be observed in the radar profiles in the shallow coves. In some seismic profiles, a thin layer of sediment could be observed at the water bottom. The radar profiles appeared to be better than the seismic profiles for the determination of sediment thickness. Although only selected data profiles were processed, all data were archived for future interpretation. This investigation focused on selected regions of the study area, particularly in the coves where sediment accumulations were presumed to be thickest. GPR was the most useful tool for interpreting sediment thickness, especially in these shallow coves. The radar profiles were interpreted for two surfaces of interest-the water bottom, which was defined as the "2007 horizon," and the interface between Lake sediments and the original Lake bottom, which was defined as the "1925 horizon"-corresponding to the year the Lake was impounded. The ground-penetrating radar data were interpreted on the basis of characteristics of the reflectors. The sediments that had accumulated in the impounded Lake were characterized by laminated, parallel reflections, whereas the subsurface below the original Lake bottom was characterized by more discontinuous and chaotic reflections, often with diffractions indicating cobbles or boulders. The reflectors were picked manually along the water bottom and along the interface between the Lake sediments and the pre-Lake sediments. A simple graphic approach was used to convert traveltimes to depth through water and depth through saturated sediments using velocities of the soundwaves through the water and the saturated sediments. Nineteen cross sections were processed and interpreted in 9 coves around Deep Creek Lake, and the difference between the 2007 horizon and the 1925 horizon was examined. In most areas, GPR data indicate a layer of sediment between 1 and 7 feet thick. When multiple cross sections from a single cove were compared, the cross sections indicated that sediment thickness decreased toward the center of the Lake.

Maryland↗
USGS1984

Evaluation of the ground-water resources of coastal Georgia: preliminary report of the data available as of July 1983

A compilation of ground-water data that have been collected for nearly 100 years in the coastal area of Georgia is presented in this report. The compilation of pertinent data indicates what information is available for use in the evaluation of the ground-water resources of the 13 counties of coastal Georgia. Also included in this report is a fairly complete discussion of previous and ongoing investigations and monitoring networks, and an extensive list of references. Maps at 1:24,000 and 1:1,000,000 scales contain well locations and identifiers for all wells in the Ground Water Site Inventory (GWSI) data base of the National Water Data Storage and retrieval System (WATSTORE). Tabular summaries of selected site information from GWSI, including well identifiers and names, latitude-longitude location, depth of well, altitude of land surface, and use of water are presented. Water-use data from the National Water Use Data System, and water use for irrigation from the University of Georgia, Department of Agriculture survey, also are tabulated. Also included are pertinent information on geophysical surveys and data obtained, and proposed project activities, particularly test-monitor well drilling. The data in this report were collected and compiled as part of the cooperative activities between the U.S. Geological Survey and other agencies.

Georgia↗
USGS2022

Can lava flow like water? Assessing applications of critical flow theory to channelized basaltic lava flows

Flowing lava and water have dramatically different physical properties but can form similar hydraulic structures, including undular hydraulic jumps, or standing wave trains. In water flows, undular hydraulic jumps are evidence of critical flow (Froude number ∼1) and open-channel hydraulic theory provides a powerful tool for estimating flow depth and velocity. Monitoring these parameters in an active lava channel is inherently challenging, but essential for calculating lava discharge (effusion rate), a primary control on the rate of flow front advance and ultimate flow runout distance. We analyze undular hydraulic jumps in both water and lava flows to assess the conditions under which they form and, by extension, the potential use of critical flow theory to estimate, in real time, lava flow velocity, depth, and discharge. Experimental data for water flows show that these structures mark the transition from supercritical to subcritical flow. Undular hydraulic jumps in the near-vent lava channel of the 2018 lower East Rift Zone eruption of Kīlauea, Hawaiʻi also reflect critical flow conditions; their wavelengths scale with flow depth and velocity, consistent with hydraulic theory. Calculated lava effusion rates are similar to estimates made using more traditional approaches (Jeffreys', 1925, https://doi.org/10.1080/14786442508634662 , equation based on lava viscosity, density, and channel slope) and with lava volumes derived from topographic-change mapping. From this we conclude that critical flow phenomena show great potential to track flow dynamics and inform hazard assessment for a wide range of geophysical fluids.

Journal of Geophysical Research - Earth Surface↗
USGS2012

Permafrost, heat flow, and the geothermal regime at Prudhoe Bay, Alaska

Temperature measurements through permafrost in the oil field at Prudhoe Bay, Alaska, combined with laboratory measurements of the thermal conductivity of drill cuttings permit an evaluation of in situ thermal properties and an understanding of the general factors that control the geothermal regime. A sharp contrast in temperature gradient at ∼600 m represents a contrast in thermal conductivity caused by the downward change from interstitial ice to interstitial water at the base of permafrost under near steady state conditions. Interpretation of the gradient contrast in terms of a simple model for the conductivity of an aggregate yields the mean ice content (∼39%), and thermal conductivities for the frozen and thawed sections (8.1 and 4.7 mcal/cm s °C, respectively). These results yield a heat flow of ∼1.3 HFU, which is similar to other values on the Alaskan Arctic Coast; the anomalously deep permafrost is a result of the anomalously high conductivity of the siliceous ice-rich sediments. Curvature in the upper 160 m of the temperature profiles represents a warming of ∼1.8°C of the mean surface temperature and a net accumulation of 5–6 kcal/cm 2 by the solid earth surface during the last 100 years or so. Rising sea level and thawing of ice-rich sea cliffs probably caused the shoreline to retreat tens of kilometers in the last 20,000 years, inundating a portion of the continental shelf that is presently the target of intensive oil exploration. A simple conduction model suggests that this recently inundated region is underlain by near-melting ice-rich permafrost to depths of 300–500 m; its presence is important to seismic interpretations in oil exploration and to engineering considerations in oil production. With confirmation of the permafrost configuration by offshore drilling, heat conduction models can yield reliable new information on the chronology of arctic shorelines.

Journal of Geophysical Research Solid Earth↗
USGS1989

Isostatic residual gravity and crustal geology of the United States

A new isostatic residual gravity map of the conterminous United States presents continent-wide gravity data in a form that can be readily used, with geologic information and other geophysical data, in studies of the composition and structure of the continental crust. This map was produced from the gridded gravity data on which the recently released Gravity Anomaly Map of the United States is based. About 1 million onland and 0.8 million offshore gravity observations interpolated to a 4- by 4-km grid serve as the basis for both maps. The Airy-Heiskanen model of isostatic compensation of topography applied to topographic and bathymetric data averaged over 5- by 5-min compartments was used to remove, to first order, the large, long-wavelength Bouguer gravity anomalies caused by deep density distributions that support topographic loads. The parameters used in the Airy-Heiskanen model were topographic density, 2.67 g/cm 3 ; sea-level crustal thickness, 30 km; and density contrast across the base of the model crust, 0.35 g/cm 3 . Many of the conspicuous short-wavelength anomalies (widths less than several hundred kilometers) on the isostatic residual gravity map correlate with mapped or near-surface geologic features, and primarily reflect shallow-density distributions rather than any departures from isostatic equilibrium. In general, gravity highs occur over (1) mafic igneous bodies emplaced in rift or magmatic arc settings or as isolated intrusions controlled by structures; (2) accreted slices of mafic oceanic, island-arc, or transitional crust; and (3) uplifted crystalline basement. Gravity lows are found over (1) thick bodies of felsic intrusive or extrusive rocks; (2) sedimentary deposits in extensional, convergent, or transform settings; and (3) depressed crystalline basement. Anomalies with widths as much as 1,000 km or more also appear to reflect crustal properties in many cases—several broad gravity highs are associated with crust having a high average seismic wave velocity, and comparable broad gravity lows occur over areas of low average seismic velocity. Alternative ways of viewing the isostatic residual gravity data provide additional information about density distributions in the crust. The first-vertical derivative map accentuates gravity anomalies over shallow, abrupt density changes at the expense of those resulting from deeper or more gradual density transitions. The maximum horizontal gradient map contains information about the locations of pronounced density boundaries. Two-dimensional spectral analysis of the gravity data provides a quantitative means for identifying dominant fabrics in the gravity field and for distinguishing various terranes from each other. Neither Bouguer nor isostatic residual gravity anomalies are particularly well suited for practical modeling of deep structure in conjunction with deep seismic information. However, a scheme in which the entire Earth outside the area of interest is approximated by laterally homogeneous layers and isostatically compensated topography, and in which the area of interest is modeled using the seismic constraints applied in a two-and-one-half-dimensional geometry, holds promise for exploiting useful features of both the Bouguer and isostatic residual gravity anomalies.

conterminous United States↗
USGS2019

Shear-wave seismic reflection studies of unconsolidated sediments in the near surface

We have successfully applied of SH-wave seismic reflection methods to two different near-surface problems targeting unconsolidated sediments. At the former Fort Ord, where the water table is approximately 30m deep, we imaged aeolian and marine aquifer and aquitard stratigraphy to a depth of approximately 80m. We identified reflections from sand/clay and sand/silt interfaces and we mapped these interfaces along our transects. At an aggregate study site in Indiana, where the water table is at a depth of 1to2m, we imaged stratigraphy in alluvial sand and gravel, and observe a strong reflection from the 32-m-deep bedrock surface. In both cases, we exploited the high resolution potential of SH waves, their insensitivity to water content, and the possibility of reducing Love wave contamination by working along a roadway. We accomplished our results using only sledgehammer sources and simple data processing flows.

California, Indiana↗
USGS2012

Airborne electromagnetic imaging of discontinuous permafrost

The evolution of permafrost in cold regions is inextricably connected to hydrogeologic processes, climate, and ecosystems. Permafrost thawing has been linked to changes in wetland and lake areas, alteration of the groundwater contribution to streamflow, carbon release, and increased fire frequency. But detailed knowledge about the dynamic state of permafrost in relation to surface and groundwater systems remains an enigma. Here, we present the results of a pioneering ∼1,800 line-kilometer airborne electromagnetic survey that shows sediments deposited over the past ∼4 million years and the configuration of permafrost to depths of ∼100 meters in the Yukon Flats area near Fort Yukon, Alaska. The Yukon Flats is near the boundary between continuous permafrost to the north and discontinuous permafrost to the south, making it an important location for examining permafrost dynamics. Our results not only provide a detailed snapshot of the present-day configuration of permafrost, but they also expose previously unseen details about potential surface – groundwater connections and the thermal legacy of surface water features that has been recorded in the permafrost over the past ∼1,000 years. This work will be a critical baseline for future permafrost studies aimed at exploring the connections between hydrogeologic, climatic, and ecological processes, and has significant implications for the stewardship of Arctic environments.

Alaska↗
USGS2026

Dynamic coupling between faulting, rifting and magmatism during 2021-2025 unrest on Reykjanes Peninsula, Iceland

Interactions among faulting, earthquakes, and eruptions are fundamental to plate tectonics and hazard forecasting yet rarely observed along mid-ocean ridges. On Iceland's Reykjanes Peninsula, seismotectonic–volcanic unrest resumed after nearly 800-year hiatus, providing an opportunity to observe these interactions during 2021–2025 activity. By integrating high-resolution seismicity, focal mechanisms, satellite geodesy, surface deformation, and eruption data, we document ∼4 m of total extension accommodated through 14 rifting episodes. The largest, in 2023, involved graben reactivation and diking, with seismic swarms and earthquake faulting that matched the surface ruptures, where strike-slip faulting preceded normal-faulting earthquakes and extension. The accrued extension was released by extension fractures triggered by magma accumulation. Long-term observations show no correlation between erupted magma volume, seismicity, and crustal extension. This highlights dynamic relation between rifting, faulting, and magmatism in transtensional settings and their implications for hazard assessment.

Reykjanes Peninsula↗