USGS Science⌕ Search

SEARCH · USGS Science

Results for “Modern Geology”

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 451 records · Page 25Linked to original sources
USGS2025

Hydrologic budgets and water availability of six bedrock aquifers in the Black Hills area, South Dakota and Wyoming, 1931–2022

Population growth and recurring droughts in the Black Hills region raised interest in water resources and future availability. The Black Hills hydrology study (BHHS) was initiated in the early 1990s to address questions regarding water resources. Since completion of the BHHS in the early 2000s, the population of the Black Hills region increased by about 39 percent, which has renewed interest in water demand and availability in the Black Hills. The U.S. Geological Survey, in cooperation with the Western Dakota Regional Water System, completed a study to update hydrologic budgets from the BHHS for six of the most used aquifers in the Black Hills. Water availability was determined by comparing results from hydrologic budgets to modern well withdrawals (2003–22) and water rights information. Key updates to the BHHS budgets included adding available data from 1999 to 2022 and determining hydrologic budgets for six aquifers in nine smaller areas (called “subareas”). Inflows for the hydrologic budget included recharge from precipitation and streamflow losses to aquifers. Total mean annual recharge for the six aquifers in the study area was estimated at 278,900 acre-feet, with 205,100 acre-feet from precipitation recharge and 73,800 acre-feet from streamflow recharge. Mean annual precipitation recharge for the Madison and Minnelusa aquifers together accounted for 76 percent of the total mean annual precipitation recharge, with the Madison aquifer contributing 57,000 acre-feet and the Minnelusa aquifer contributing 98,100 acre-feet. Outflow components estimated for the hydrologic budget include artesian springflow and well withdrawals. Total mean annual artesian springflow in the study area was estimated as 166,100 acre-feet for the combined Madison and Minnelusa aquifers. Mean total annual well withdrawals for 2003–22 in the study area were about 50,000 acre-feet. No increased well withdrawal patterns corresponding to population increases were observed between 2003 and 2022. Water availability was determined by comparing total annual appropriations and mean and maximum annual well withdrawals for 2003–22 to mean annual recharge for 1931–2022 for each aquifer in subareas 1–9. Modern well withdrawals (mean and maximum for 2003–22) exceeded mean annual recharge for only the Deadwood and Inyan Kara aquifers in subareas 9 and 4, respectively. Additionally, total annual appropriations did not exceed mean annual recharge in most subareas, except most notably in subarea 4 (Rapid City area) where appropriations exceeded recharge for the Madison, Minnelusa, and Inyan Kara aquifers. Total annual appropriations also exceeded mean annual recharge for the Inyan Kara aquifer in subareas 3 and 5. In addition to recharge, water availability includes the water stored in pore spaces of aquifer materials. Estimates of total volume of recoverable water in storage were updated as part of this study to include the portion of aquifers in Wyoming, which were omitted during the BHHS. In total, the estimated total amount of recoverable water in storage in the study area was 356.9 million acre-feet for six major aquifers in the Black Hills area of South Dakota and Wyoming.

South Dakota, Wyoming↗
USGS2004

The Evolving Landscape of the Columbia River Gorge: Lewis and Clark and Cataclysms on the Columbia

TAVELERS RETRACING LEWIS AND CLARKE JOURNEY to the Pacific over the past two hundred years have witnessed tre mendous change to the Columbia River Gorge and its pri mary feature, the Columbia River. Dams, reservoirs, timber harvest, altered fisheries, transportation infrastructure, and growth and shrinkage of communities have transformed the river and valley.1 This radically different geography of human use and habitation is commonly contrasted with the sometimes romantic view of a prior time provided both by early nineteenth-century chroniclers and present-day critics of the modern condition ? an ecotopia of plentiful and perpetual resources sustaining a stable culture from time immemorial. Reality is more com plicated. Certainly the human-caused changes to the Columbia River and the gorge since Lewis and Clark have been profound; but the geologic his tory of immense floods, landslides, and volcanic eruptions that occurred before their journey had equally, if not more, acute effects on landscapes and societies of the gorge. In many ways, the Lewis and Clark Expedi tion can be viewed as a hinge point for the Columbia River, the changes engineered to the river and its valley in the two hundred years since their visit mirrored by tremendous changes geologically engendered in the thousands of years before.

Oregon Historical Society Quarterly↗
USGS2012

Modification of wave-cut and faulting-controlled landforms

From a casual observation that the form of degraded fault scarps resembles the error function, this investigation proceeds through an elementary diffusion equation representation of landform evolution to the application of the resulting equations to the modern topography of scarplike landforms. The morphologic observations can be analyzed either in the form of one or more cross-strike elevation profiles or in the form of the slope-offset plot, a point plot of maximum scarp slope versus scarp offset. Working with either or both of these data representations for nine geologic structures, which range in age from 3 to 400 ka B.P. and in offset from 1 to 50 m, we apply analytical solutions for the vertical initial value scarp, the vertical continuous offset scarp, and the finite slope, initial value scarp. The model calculations are intrinsically ambiguous, yielding as the final answer only the product κ t (in the case of the initial value problem) or the product κ A −1 (in the case of the repeated faulting problem); here t is the age of a single scarp-forming event, 2 A is the vertical slip rate, and κ is the “mass diffusivity.” A single profile across three sea cliffs along the Santa Cruz, California, coast is analyzed as three separate initial value problems. A reasonably constrained age for the sea cliff standing above the Highway 1 platform returns κ = 11 GKG (1 GKG = 1 m 2 /ka). With this κ, we can date the two older sea cliffs. In fact, we do the converse: age estimates for these two older sea cliffs based on a uniform rate of uplift both yield the same κ as for the lower sea cliff. We treat a single profile of the Raymond fault in Pasadena/San Marino in terms of the repeated faulting problem; for it the uplift rate of R. Crook and others yields κ = 16 GKG. The very substantial preexisting offset across the Raymond fault must have been buried/leveled some 230 ka B.P., when the modern topography began to form. Our analysis of the Lake Bonneville shoreline scarps reveals a dependence of κ t on 2a, suggestive of nonlinear modification processes. This appearance is treated with the finite slope initial value scarp model to determine κ=1.1 GKG for the Lake Bonneville shoreline scarps. The suggestion of M. N. Machette that approximately 100,000-year-old, meter-high scarps are “unobservable” in weakly consolidated alluvial terranes of the Basin and Range and Rio Grande Rift Valley provinces can be formulated as κ ≳ 1 GKG. The coincidence between this inequality and the Lake Bonneville shoreline κ is striking, and it suggests that the value of κ = 1 GKG may be generally applicable, as a good first approximation, to the modification of alluvial terranes within the semiarid regions of the western United States. The Lake Bonneville shoreline κ is the basis for dating four sets of fault scarps in west-central Utah. The Drum Mountains fault scarps can be modeled in several different circumstances, but the most likely interpretation is that these fault scarps formed as the result of a single episode of normal faulting 3.6 to 5.7 ka B.P. The younger age is associated with quite low initial slope angles (25°). The other three sets of fault scarps show no evidence for finite initial value slopes. Fault scarps along the eastern base of the Fish Springs Range are very young, 3 ka B.P. We estimate the age of fault scarps along the western flank of the Oquirrh Mountains to be 32 ka B.P., which meets the weak geologic constraint that they be older than the Lake Bonneville shoreline. Fault scarps along the northeastern margin of the Sheeprock Mountains are even older, 53 ka B.P. An intriguing consequence of our single-event analysis of these scarps is that an 11.5-m offset occurred in a single earthquake.

Journal of Geophysical Research Solid Earth↗
USGS1965

Late quaternary geologic history of the lower Chippewa Valley, Wisconsin

The lower Chippewa Valley in west-central Wisconsin extends 65 miles from the Cary terminal moraine in Chippewa County to the Mississippi River Valley. The Chippewa Valley and its tributaries were filled with a valley train of sand and gravel during the maximum stand of the Cary ice, and entrenchment of this deposit has formed the Wissota terrace, a prominent geomorphic feature that can be traced the length of the valley. Several lower terraces in the valley indicate progressive downcutting of the Wissota terrace sediments. Erosion and deposition in the Mississippi Valley are closely linked to the post-Cary history of the lower Chippewa Valley, for these factors controlled the outlet level of the Chippewa River. This outlet was substantially lower than at present throughout much of post-Cary Pleistocene and early Recent time. The modern Chippewa River has built a delta into the Mississippi Valley. The Chippewa River is aggrading the lower part of its valley, a meandering river is slowly eroding the central part; stream erosion in the upper part is restricted by sills of hard bedrock.

Wisconsin↗
USGS2013

Mass fractionation of noble gases in synthetic methane hydrate: Implications for naturally occurring gas hydrate dissociation

As a consequence of contemporary or longer term (since 15 ka) climate warming, gas hydrates in some settings may presently be dissociating and releasing methane and other gases to the ocean-atmosphere system. A key challenge in assessing the impact of dissociating gas hydrates on global atmospheric methane is the lack of a technique able to distinguish between methane recently released from gas hydrates and methane emitted from leaky thermogenic reservoirs, shallow sediments (some newly thawed), coal beds, and other sources. Carbon and deuterium stable isotopic fractionation during methane formation provides a first-order constraint on the processes (microbial or thermogenic) of methane generation. However, because gas hydrate formation and dissociation do not cause significant isotopic fractionation, a stable isotope-based hydrate-source determination is not possible. Here, we investigate patterns of mass-dependent noble gas fractionation within the gas hydrate lattice to fingerprint methane released from gas hydrates. Starting with synthetic gas hydrate formed under laboratory conditions, we document complex noble gas fractionation patterns in the gases liberated during dissociation and explore the effects of aging and storage (e.g., in liquid nitrogen), as well as sampling and preservation procedures. The laboratory results confirm a unique noble gas fractionation pattern for gas hydrates, one that shows promise in evaluating modern natural gas seeps for a signature associated with gas hydrate dissociation.

Chemical Geology↗
USGS1992

Diagenetic formation of bedded chert: Evidence from chemistry of the chert-shale couplet

Theories concerning the formation of bedded chert traditionally have emphasized either depositional or diagenetic processes. Major and rare earth element data from Franciscan assemblage (Mesozoic) and Claremont Formation (Miocene) bedded chert sequences, along with physical observations such as the presence of rare and highly corroded radiolarians in shale interbeds, are most consistent with a dominantly diagenetic origin of chert-shale couplets and are incompatible with many depositional theories. Chemical distributions between Franciscan and Claremont bedded chert=shale closely match chemical fractionations recorded by Monterey Formation and Deep Sea Drilling Project-sampled cherts formed by diagenetic SiO 2 dissolution, transport, and reprecipitation, suggesting that diagenetic migration of SiO 2 from proto-shale to proto-chert is also largely responsible for chert-shale couplets. Identical Ce anomalies (Ce/Ce*) found in immediately adjacent chert-shale layers indicate that turbidites or other transport mechanisms are not responsible for the alternating beds. Neither the chemistry of the chert-shale couplet nor the overall stratigraphy of the sequences is consistent with couplet formation being caused by productivity fluctuations. Chemical mass balance calculations reconstructing the total bulk sediment composition suggest that modern siliceous sequences do not contain enough labile biogenic SiO 2 to form entire stratigraphies of bedded chert.

Geology↗
USGS2010

Snowmelt hydrograph interpretation: Revealing watershed scale hydrologic characteristics of the Yellowstone volcanic plateau

Snowmelt hydrograph analysis and groundwater age dates of cool water springs on the Yellowstone volcanic plateau provide evidence of high volumes of groundwater circulation in watersheds comprised of quaternary Yellowstone volcanics. Ratios of maximum to minimum mean daily discharge and average recession indices are calculated for watersheds within and surrounding the Yellowstone volcanic plateau. A model for snowmelt recession is used to separate groundwater discharge from overland runoff, and compare groundwater systems. Hydrograph signal interpretation is corroborated with chlorofluorocarbon (CFC) and tritium concentrations in cool water springs on the Yellowstone volcanic plateau. Hydrograph parameters show a spatial pattern correlated with watershed geology. Watersheds comprised dominantly of quaternary Yellowstone volcanics are characterized by slow streamflow recession, low maximum to minimum flow ratios. Cool springs sampled within the Park contain CFC's and tritium and have apparent CFC age dates that range from about 50 years to modern. Watersheds comprised of quaternary Yellowstone volcanics have a large volume of active groundwater circulation. A large, advecting groundwater field would be the dominant mechanism for mass and energy transport in the shallow crust of the Yellowstone volcanic plateau, and thus control the Yellowstone hydrothermal system. ?? 2009 Elsevier B.V.

Journal of Hydrology↗
USGS2013

Geologic effects on groundwater salinity and discharge into an estuary

Submarine groundwater discharge (SGD) can be an important pathway for transport of nutrients and contaminants to estuaries. A better understanding of the geologic and hydrologic controls on these fluxes is critical for their estimation and management. We examined geologic features, porewater salinity, and SGD rates and patterns at an estuarine study site. Seismic data showed the existence of paleovalleys infilled with estuarine mud and peat that extend hundreds of meters offshore. A low-salinity groundwater plume beneath this low-permeability fill was mapped with continuous resistivity profiling. Extensive direct SGD measurements with seepage meters ( n = 551) showed fresh groundwater discharge patterns that correlated well with shallow porewater salinity and the hydrogeophysical framework. Small-scale variability in fresh and saline discharge indicates influence of meter-scale geologic heterogeneity, while site-scale discharge patterns are evidence of the influence of the paleovalley feature. Beneath the paleovalley fill, fresh groundwater flows offshore and mixes with saltwater before discharging along paleovalley flanks. On the adjacent drowned interfluve where low-permeability fill is absent, fresh groundwater discharge is focused at the shoreline. Shallow saltwater exchange was greatest across sandy sediments and where fresh SGD was low. The geologic control of groundwater flowpaths and discharge salinity demonstrated in this work are likely to affect geochemical reactions and the chemical loads delivered by SGD to coastal surface waters. Because similar processes are likely to exist in other estuaries where drowned paleovalleys commonly cross modern shorelines, the existence and implications of complex hydrogeology are important considerations for studies of groundwater fluxes and related management decisions.

Delaware↗
USGS2006

Attention turns to naturally occurring methane seepage

Methane is the most abundant organic compound in the Earth's atmosphere. As a powerful greenhouse gas, it has implications for global climate change. Sources of methane to the atmosphere are varied. Depending on the source, methane can contain either modern or ancient carbon. Methane exiting from swamps and wetlands contains modern carbon, whereas methane leaking from petroleum reservoirs contains ancient carbon. The total annual source of methane to the atmosphere has been constrained to about 540 teragrams (Tg) per year “ Cicerone and Oremland , 1988”. Notably absent from any identified sources is the contribution of geologically sourced methane from naturally occurring seepage.

Eos, Earth and Space Science News↗
USGS2016

Spatial and temporal assessment of back-barrier erosion on Cumberland Island National Seashore, Georgia, 2011–2013

Much research has been conducted to better understand erosion and accretion processes for the seaward zones of coastal barrier islands; however, at Cumberland Island National Seashore, Georgia, the greater management concern is the effect that erosion is having on the resources of the island’s western shoreline, or the back barrier. Catastrophic slumping and regular rates of erosion greater than 1 meter per year threaten important habitat, historical and pre-historical resources, and modern infrastructure on the island. Prior research has helped National Park Service (NPS) staff identify the most severe and vulnerable areas, but in order to develop effective management actions, information is needed on what forces and conditions cause erosion. To this end, the U.S. Geological Survey, in cooperation with the NPS, conducted two longitudinal surveys, one each at the beginning and end of the approximately year-long monitoring period from late 2011 to early 2013, along five selected segments of the back barrier of the Cumberland Island National Seashore. Monitoring stations were constructed at four of these locations that had previously been identified as erosional hotspots. The magnitude of erosion at each location was quantified to determine the relative influence of causative agents. Results indicate that erosion is, in general, highly variable within and among these segments of the Cumberland Island National Seashore’s back barrier. Observed erosion ranged from a maximum of 2.5 meters of bluff-line retreat to some areas that exhibited no net erosion over the 1-year study period. In terms of timing of erosion, three of the four sites were primarily affected by punctuated erosional events that were coincident with above-average high tides and elevated wind speeds. The fourth site exhibited steady, low-magnitude retreat throughout the study period. While it is difficult to precisely subscribe certain amounts of erosion to specific agents, this study provides insight into the mode of erosion among sites and the interaction among factors that set up conditions that may be leading to punctuated events. Estimates of sea-level rise were incorporated into the results of this study to project conditions that could be in place by the end of the 21st century. When using the erosion rates observed in this study to extrapolate future shoreline position, results indicate an average retreat (across all monitored locations) of 15 meters by 2050 and approximately 37 meters by 2100.

Georgia↗
USGS2007

Surficial Geologic Map of the West Franklin Quadrangle, Vanderburgh and Posey Counties, Indiana, and Henderson County, Kentucky

The valley of the Ohio River is filled with alluvium and outwash (unit Qal), which total 33-39 m thick under the land surface in the southeast part of the West Franklin quadrangle in Indiana, and 30.5-35 m thick under Diamond Island in the southwest corner of the quadrangle. The deposits are chiefly fine- to medium-grained, lithic quartzose sand, interbedded by lenses of clay, clayey silt, silt, coarse sand, granules, and gravel. Although grain size of the river alluvium varies widely, in general it fines upward-being gravelly sand to sandy gravel in the lower part, mainly sand in the middle part, and silty and clayey in the upper part (Holocene). The middle and lower parts probably accumulated during the Wisconsin Episode (late Pleistocene). The sandy middle part contains interbeds of clay, silt, and minor gravel. At the base is highly consolidated mud (silt and clay), sand, and gravel 2-10 m thick. This unit may be valley train that predates the Wisconsin Episode. Creek alluvium (unit Qa) is silt, clayey silt, and subordinate intercalated fine sand, granules, and pebbles; the coarser grains are generally concentrated in the basal 1-2 m of the deposit. Lenses and beds of clay are present locally. Fossil wood collected from an auger hole in the alluvial deposits of Little Creek, at depths of 10.6 m and 6.4 m, were dated 16,650?50 and 11,120?40 radiocarbon years, respectively. Probable lacustrine terrace silt and clay (Qlt), so-called slackwater-lake or backwater deposits, form deposits 12-22 m thick in the lowest reaches of tributary creeks to the Ohio River. The surfaces of the lacustrine deposits are terraces a few meters higher than the modern creek flood plains. Covering the bedrock upland is loess (Ql) 3-7.5 m thick, deposited about 18,000-12,000 years before present. Most surficial deposits in the quadrangle are probably no older than about 35,000 yrs. Lithologic logs, shear-wave velocities, and other cone penetrometer data are used to interpret depositional environments and geologic history of the alluvium and lacustrine deposits.

Scientific Investigations Map↗
USGS2019

Groundwater inflow toward a preheated volcanic conduit: Application to the 2018 eruption at Kīlauea Volcano, Hawai’i

The many successes in volcano forecasting over the past several decades owe mainly to pattern recognition, both in monitoring data and the geologic record. During the early stages of the 2018 Kīlauea eruption, the conceptual model of Stearns (1925), based on the explosive 1924 Kīlauea eruption, was highly influential. This model postulates that explosions are triggered by liquid-water inflow into a recently vacated magma conduit. Modern quantitative modeling approaches, supplemented by hydrogeologic data unavailable in 1925, yield a more nuanced view. Results demonstrate that liquid-water inflow would likely be delayed by months to years, owing to the inability of liquid water to transit a zone of very hot rock surrounding the conduit. The exercise demonstrates the use of physically based modeling to supplement traditional volcano-forecasting approaches during an ongoing event.

Hawaii↗
USGS2014

Holocene earthquakes and right-lateral slip on the left-lateral Darrington-Devils Mountain fault zone, northern Puget Sound, Washington

Sources of seismic hazard in the Puget Sound region of northwestern Washington include deep earthquakes associated with the Cascadia subduction zone, and shallow earthquakes associated with some of the numerous crustal (upper-plate) faults that crisscross the region. Our paleoseismic investigations on one of the more prominent crustal faults, the Darrington–Devils Mountain fault zone, included trenching of fault scarps developed on latest Pleistocene glacial sediments and analysis of cores from an adjacent wetland near Lake Creek, 14 km southeast of Mount Vernon, Washington. Trench excavations revealed evidence of a single earthquake, radiocarbon dated to ca. 2 ka, but extensive burrowing and root mixing of sediments within 50–100 cm of the ground surface may have destroyed evidence of other earthquakes. Cores in a small wetland adjacent to our trench site provided stratigraphic evidence (formation of a laterally extensive, prograding wedge of hillslope colluvium) of an earthquake ca. 2 ka, which we interpret to be the same earthquake documented in the trenches. A similar colluvial wedge lower in the wetland section provides possible evidence for a second earthquake dated to ca. 8 ka. Three-dimensional trenching techniques revealed evidence for 2.2 ± 1.1 m of right-lateral offset of a glacial outwash channel margin, and 45–70 cm of north-side-up vertical separation across the fault zone. These offsets indicate a net slip vector of 2.3 ± 1.1 m, plunging 14° west on a 286°-striking, 90°-dipping fault plane. The dominant right-lateral sense of slip is supported by the presence of numerous Riedel R shears preserved in two of our trenches, and probable right-lateral offset of a distinctive bedrock fault zone in a third trench. Holocene north-side-up, right-lateral oblique slip is opposite the south-side-up, left-lateral oblique sense of slip inferred from geologic mapping of Eocene and older rocks along the fault zone. The cause of this slip reversal is unknown but may be related to clockwise rotation of the Darrington–Devils Mountain fault zone into a position more favorable to right-lateral slip in the modern N-S compressional stress field.

Washington↗
USGS2012

History of surface displacements at the Yellowstone Caldera, Wyoming, from leveling surveys and InSAR observations, 1923-2008

Modern geodetic studies of the Yellowstone caldera, Wyoming, and its extraordinary tectonic, magmatic, and hydrothermal systems date from an initial leveling survey done throughout Yellowstone National Park in 1923 by the U.S. Coast and Geodetic Survey. A repeat park-wide survey by the U.S. Geological Survey (USGS) and the University of Utah during 1975-77 revealed that the central part of the caldera floor had risen more than 700 mm since 1923, at an average rate of 14±1 mm/yr. From 1983 to 2007, the USGS conducted 15 smaller surveys of a single level line that crosses the northeast part of the caldera, including the area where the greatest uplift had occurred from 1923 to 1975-77. The 1983 and 1984 surveys showed that uplift had continued at an average rate of 22±1 mm/yr since 1975-77, but no additional uplift occurred during 1984-85 (-2±5 mm/yr), and during 1985-95 the area subsided at an average rate of 19±1 mm/yr. The change from uplift to subsidence was accompanied by an earthquake swarm, the largest ever recorded in the Yellowstone area (as of March 2012), starting in October 1985 and located near the northwest rim of the caldera. Interferometric synthetic aperture radar (InSAR) images showed that the area of greatest subsidence migrated from the northeast part of the caldera (including the Sour Creek resurgent dome) during 1992-93 to the southwest part (including the Mallard Lake resurgent dome) during 1993-95. Thereafter, uplift resumed in the northeast part of the caldera during 1995-96, while subsidence continued in the southwest part. The onset of uplift migrated southwestward, and by mid-1997, uplift was occurring throughout the entire caldera (essentially rim to rim, including both domes). Consistent with these InSAR observations, leveling surveys indicated 24±3 mm of uplift in the northeast part of the caldera during 1995-98. The beginning of uplift was coincident with or followed shortly after an earthquake swarm near the north caldera rim during June-July 1995 - the strongest swarm since 1985. Rather than a single deformation source as inferred from leveling surveys, the InSAR images revealed two distinct sources - one beneath each resurgent dome on the caldera floor. Subsequently, repeated GPS surveys (sometimes referred to as "campaign" surveys to distinguish them from continuous GPS observations) and InSAR images revealed a third deformation source beneath the north caldera rim. The north-rim source started to inflate in or about 1995, resulting in as much as 80 mm of surface uplift by 2000. Meanwhile, motion of the caldera floor changed from uplift to subsidence during 1997-8. The north rim area rose, while the entire caldera floor (including both domes) subsided until 2002, when both motions paused. Uplift in the northeast part of the caldera resumed in mid-2004 at a historically unprecedented rate of as much as 70 mm/yr, while the north rim area subsided at a lesser rate. Resurveys of the level line across the northeast part of the caldera in 2005 and 2007 indicated the greatest average uplift rate since the initial survey in 1923-53±3 mm/yr. Data from a nearby continuous GPS (CGPS) station showed that the uplift rate slowed to 40-50 mm/yr during 2007-8 and to near zero by September 2009. Following an intense earthquake swarm during January-February 2010, this one near the northwest caldera rim and the strongest since the 1985 swarm in the same general area, CGPS stations recorded the onset of subsidence throughout the entire caldera. Any viable model for the cause(s) of ground deformation at Yellowstone should account for (1) three distinct deformation sources and their association with both resurgent domes and the north caldera rim; (2) interplay among these sources, as suggested by the timing of major changes in deformation mode; (3) migration of the area of greatest subsidence or uplift from the northeast part of the caldera to the southwest part during 1992-95 and 1995-97, respectively; (4) repeated cycles of uplift and subsidence and sudden changes from uplift to subsidence or vice versa; (5) spatial and temporal relationships between changes in deformation mode and strong earthquake swarms; and (6) lateral dimensions of all three deforming areas that indicate source depths in the range of 5 to 15 km. We prefer a conceptual model in which surface displacements at Yellowstone are caused primarily by variations in the flux of basaltic magma into the crust beneath the caldera. Specifically, we envision a magmatic conduit system beneath the northeast part of the caldera that supplies basalt from a mantle source to an accumulation zone at 5-10 km depth, perhaps at a rheological boundary within a crystallizing rhyolite body remnant from past eruptions. Increases in the magma flux favor uplift of the caldera and decreases favor subsidence. A delicate equilibrium exists among the mass and heat flux from basaltic intrusions, heat and volatile loss from the crystallizing rhyolite body, and the overlying hydrothermal system. In the absence of basalt input, steady subsidence occurs mainly as a result of fluid loss from crystallizing rhyolite. At times when a self-sealing zone in the deep hydrothermal system prevents the escape of magmatic fluid, the resulting pressure increase contributes to surface uplift within the caldera; such episodes end when the seal ruptures during an earthquake swarm. To account for the north rim deformation source, we propose that magma or fluid exsolved from magma episodically escapes the caldera system at the three-way structural intersection of (1) the northern caldera boundary, (2) an active seismic belt to the north-northwest that is associated with the Hebgen Lake fault zone, and (3) the Norris - Mammoth corridor - a zone of faults, volcanic vents, and thermal activity that strikes north from the north rim of the caldera near Norris Geyser Basin to Mammoth Hot Springs near the northern boundary of Yellowstone National Park. Increased fluid flux out of the caldera by way of this intersection favors subsidence of the north rim area, and decreased flux favors uplift. This model does not account for poroelastic and thermoelastic effects, nonelastic rheology, or heat and mass transport in the hot and wet subcaldera crust. Such effects almost surely play a role in caldera deformation and are an important topic of ongoing research.

Wyoming↗
USGS2014

Petrographic maturity parameters of a Devonian shale maturation series, Appalachian Basin, USA. ICCP Thermal Indices Working Group interlaboratory exercise

This paper presents results of an interlaboratory exercise on organic matter optical maturity parameters using a natural maturation series comprised by three Devonian shale samples (Huron Member, Ohio Shale) from the Appalachian Basin, USA. This work was conducted by the Thermal Indices Working Group of the International Committee for Coal and Organic Petrology (ICCP) Commission II (Geological Applications of Organic Petrology). This study aimed to compare: 1. maturation predicted by different types of petrographic parameters (vitrinite reflectance and spectral fluorescence of telalginite), 2. reproducibility of the results for these maturation parameters obtained by different laboratories, and 3. improvements in the spectral fluorescence measurement obtained using modern detection systems in comparison with the results from historical round robin exercises. Mean random vitrinite reflectance measurements presented the highest level of reproducibility (group standard deviation 0.05) for low maturity and reproducibility diminished with increasing maturation (group standard deviation 0.12). Corrected fluorescence spectra, provided by 14 participants, showed a fair to good correspondence. Standard deviation of the mean values for spectral parameters was lowest for the low maturity sample but was also fairly low for higher maturity samples. A significant improvement in the reproducibility of corrected spectral fluorescence curves was obtained in the current exercise compared to a previous investigation of Toarcian organic matter spectra in a maturation series from the Paris Basin. This improvement is demonstrated by lower values of standard deviation and is interpreted to reflect better performance of newer photo-optical measuring systems. Fluorescence parameters measured here are in good agreement with vitrinite reflectance values for the least mature shale but indicate higher maturity than shown by vitrinite reflectance for the two more mature shales. This red shift in λ max beyond 0.65% vitrinite reflectance was also observed in studies of Devonian shale in other basins, suggesting that the accepted correlation for these two petrographic thermal maturity parameters needs to be re-evaluated. A good linear correlation between λ max and T max for this maturation series was observed and λ max 600 nm corresponds to T max of 440 °C. Nevertheless if a larger set of Devonian samples is included, the correlation is polynomial with a jump in λ max ranging from 540 to 570 nm. Up to 440 °C of T max , the λ max , mostly, reaches up to 500 nm; beyond a T max of 440 °C, λ max is in the range of 580–600 nm. This relationship places the “red shift” when the onset of the oil window is reached at T max of 440 °C. Moreover, the correlation between HI and λ max (r 2 = 0.70) shows a striking inflection and decrease in HI above a λ max of 600 nm, coincident with the approximate onset of hydrocarbon generation in these rocks.

International Journal of Coal Geology↗
USGS2021

National-scale reservoir thermal energy storage pre-assessment for the United States

The U.S. Geological Survey is performing a pre-assessment of the cooling potential for reservoir thermal energy storage (RTES) in five generalized geologic regions (Basin and Range, Coastal Plains, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Reservoir models are developed for the metropolitan areas of eight cities (Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon) so that computed metrics can be compared to evaluate RTES potential across diverse climates, geologic settings, and physiography. Permeable, semi-confined/confined units that underlie more-utilized aquifers and contain low-quality groundwater are selected for each city. Energy storage metrics are computed for the anticipated total thickness of stratigraphy for which RTES might be feasible, including estimated required well spacing, thermal storage capacity, and thermal recovery efficiency over time. Falta et al. (2016) showed that for a modern 25,000 square-foot (2,323 square-meter), two-story office building, cooling needs exceed heating demand for almost every region of the country. We therefore use Falta et al.’s cooling demand for each city as the representative RTES stress condition for metric computation, allowing comparisons across regions. Results indicate that favorable RTES conditions exist in each region, particularly in the Illinois Basin, Coastal Plains, and Basin and Range. Thermal recovery efficiencies are very high in all regions and increase over time. The thermal storage capacity metric is most informative in the pre-assessment and underscores the importance of mapping reservoir thicknesses and porosities to permit detailed mapping of thermal storage capacity per unit area as a key RTES resource classification standard. This assessment provides a basic understanding of the RTES potential in several metropolitan areas and geologic regions throughout the United States and will aid further evaluation of national RTES efficacy.

Conference Paper↗
USGS2018

Titanium mineral resources in heavy-mineral sands in the Atlantic coastal plain of the southeastern United States

This study examined titanium distribution in the Atlantic Coastal Plain of the southeastern United States; the titanium is found in heavy-mineral sands that include the minerals ilmenite (Fe 2+ TiO 3 ), rutile (TiO 2 ), or leucoxene (an alteration product of ilmenite). Deposits of heavy-mineral sands in ancient and modern coastal plains are a significant feedstock source for the titanium dioxide pigments industry. Currently, two heavy-mineral sands mining and processing operations are active in the southeast United States producing concentrates of ilmenite-leucoxene, rutile, and zircon. The results of this study indicate the potential for similar deposits in many areas of the Atlantic Coastal Plain. This study used the titanium analyses of 3,457 stream sediment samples that were analyzed as part of the U.S. Geological Survey’s National Geochemical Survey program. This data set was analyzed by an integrated spatial modeling technique known as Bayesian hierarchical modeling to map the regional-scale, spatial distribution of titanium concentrations. In particular, clusters of anomalous concentrations of titanium occur: (1) along the Fall Zone, from Virginia to Alabama, where metamorphic and igneous rocks of the Piedmont region contact younger sediments of the Coastal Plain; (2) a paleovalley near the South Carolina and North Carolina border; (3) the upper and middle Atlantic Coastal Plain of North Carolina; (4) the majority of the Atlantic Coastal Plain of Virginia; and (5) barrier islands and stretches of the modern shoreline from South Carolina to northeast Florida. The areas mapped by this study could help mining companies delimit areas for exploration.

Scientific Investigations Report↗
USGS2003

High resolution seismic imaging of faults beneath Limón Bay, northern Panama Canal, Republic of Panama

High-resolution seismic reflection profiles from Limo??n Bay, Republic of Panama, were acquired as part of a seismic hazard investigation of the northern Panama Canal region. The seismic profiles image gently west and northwest dipping strata of upper Miocene Gatu??n Formation, unconformably overlain by a thin (<20 m) sequence of Holocene muds. Numerous faults, which have northeast trends where they can be correlated between seismic profiles, break the upper Miocene strata. Some of the faults have normal displacement, but on many faults, the amount and type of displacement cannot be determined. The age of displacement is constrained to be Late Miocene or younger, and regional geologic considerations suggest Pliocene movement. The faults may be part of a more extensive set of north- to northeast-trending faults and fractures in the canal region of central Panama. Low topography and the faults in the canal area may be the result of the modern regional stress field, bending of the Isthmus of Panama, shearing in eastern Panama, or minor deformation of the Panama Block above the Caribbean subduction zone. For seismic hazard analysis of the northern canal area, these faults led us to include a source zone of shallow faults proximal to northern canal facilities. ?? 2003 Elsevier B.V. All rights reserved.

Tectonophysics↗