USGS Science⌕ Search

SEARCH · USGS Science

Results for “Proceedings of the National Academy of Sciences”

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

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

1,688 records · Page 33Linked to original sources
USGS2026

Hydrogeologic framework and conceptual groundwater-flow model of the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, 1998–2022

This study was conducted by the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, to update the hydrogeologic framework and conceptual flow model for the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, which together compose the Ogallala aquifer focus area. The study included the construction of a potentiometric surface, and available geologic and hydrologic data were used to evaluate saturated thickness of the aquifer. The water budget for the updated conceptual groundwater-​flow model was based on estimated inflows and outflows for the 1998–2022 study period. Saturated thickness of the Ogallala aquifer averaged 127 and 116 feet for the panhandle and northwest parts, respectively. Groundwater withdrawals from the Ogallala aquifer for 1998–2022 averaged 422,054 and 39,645 acre-feet per year (acre-ft/yr) for the panhandle and northwest parts, respectively. Recharge, the primary inflow, was estimated at 0.63 inch per year for the 1998–2022 study period, with the panhandle part of the Ogallala aquifer receiving 175,068 acre-ft/yr and the northwest part of the Ogallala aquifer receiving 49,376 acre-ft/yr. Additional inflows included irrigation return flows, estimated at 8,111 and 642 acre-ft/yr for the panhandle and northwest parts, respectively, of the Ogallala aquifer. Net lateral groundwater flows, considered to be aquifer outflows, were estimated to account for 31,908 acre-ft/yr for the Ogallala aquifer focus area. Streambed seepage, which was an outflow of 5,535 acre-ft/yr, was only present in the northwest part of the Ogallala aquifer. Vertical leakage and saturated-zone evapotranspiration were considered negligible outflows. These findings provide a revised conceptual groundwater-flow model water budget for the Ogallala aquifer focus area in Oklahoma.

Colorado, Kansas, Oklahoma, Texas↗
USGS1976

Transport and dispersion of fluorescent tracer particles for the dune-bed condition, Atrisco Feeder Canal near Bernalillo, New Mexico

A tracer technique in which mineral particles were coated with fluorescent dyes was used to study the rates of transport and dispersion of sediment particles`of various diameters and specific gravities for a dune-bed condition in an alluvial channel. The experiment was conducted in the Atrisco Feeder Canal near Bernalillo, N. Mex., between May 1 and July 14, 1967. A continuous point source of tracers, approximated by injections at 5- or 10-minute intervals, was maintained for 7 days so that the steady-dilution procedure could be used to calculate the transport rate of bed material. After termination of the injection process, the spatial-integration procedure was used to follow the movement of the tracers downstream and to calculate the transport rate. Samples of the bed material in transport and the accompanying tracers moving along the surface of the dune bed were obtained periodically throughout the study with the "dustpan" sampler especially designed for fluorescent tracer studies. In addition, the spatial distributions of the tracers in the dune bed were determined three times during the study by core sampling.

New Mexico↗
USGS2025

Natural source zone depletion of crude oil in the subsurface: Processes controlling mass losses of individual compounds

At many petroleum hydrocarbon spill sites, residual spilled product forms a long-term source of groundwater contamination. The phrase source zone natural depletion is used to refer to the mass loss rates. Overall mass lost under environmental conditions was analyzed using conservative biomarker concentrations for a 1979 oil spill in northern Minnesota, USA. After 40–41 years, an average of 50% of the mass was lost with values ranging from 22% to 57% depending on location. It is also important to understand the composition changes in the source. To understand controls on the losses of individual compounds, concentrations of volatile hydrocarbons in oil samples were compared with aqueous solubilities, and pore-space oil saturations. The results of the comparison show that losses of the oil compounds were controlled by pore-space oil saturations, solubility, and susceptibility to degradation under methanogenic conditions. Compounds that degrade under methanogenic conditions, including toluene, o -xylene, and n -alkanes are more depleted compared to benzene, ethylbenzene, and m - and p -xylene for which losses are dominated by dissolution. These rates and compound-specific behaviors form a foundation for improved modeling approaches and risk analyses.

Minnesota↗
USGS2025

Horizontal transport of Picture Gorge Basalt magma through the Monument Dike Swarm determined by magnetic fabric

Flood basalts of the mid-Miocene Columbia River Basalt Group (CRBG) cover 210,000 km 2 of Washington, Oregon, and Idaho. The source of CRBG melt is debated; widely spaced feeder dike swarms can be projected toward hypothetical sources near the Oregon-Idaho border. In this study, we use anisotropy of magnetic susceptibility (AMS) to track magma flow in the Monument dike swarm (MDS), the feeder dikes of the Picture Gorge Basalt (PGB). This small formation of the main-phase CRBG eruptions allows us to explore in detail the localized dynamics of a large igneous province feeder system, with implications for the larger CRBG picture. We measured the magnetic fabric of 205 oriented paleomagnetic specimens subsampled from 97 samples collected from 15 dikes of the MDS. Thermal demagnetization and hysteresis loops show that the magnetic minerals are a mixture of single domain and multidomain sized titanomagnetites. At three dikes, the paleodepth of sampling was determined to be shallow (<350 m). Magma flowing through dikes has been shown—in most cases— to acquire an anisotropic magnetic fabric with an AMS ellipsoid minimum axis perpendicular to the wall and maximum axis aligned in the direction of flow. Of 15 dikes, 12 show horizontal flow directions in the plane of the dike. Only one dike displayed imbricated fabrics, showing westward flow away from the Oregon-Idaho border. We conclude that magma flow in the MDS was sub-horizontal from a distal source.

Oregon↗
USGS2026

Shallow differentiation of primitive arc magmas at the Jurassic Emigrant Gap mafic complex, Sierra Nevada, California

The Emigrant Gap composite pluton exposes ultramafic to silicic intrusive rocks that preserve the chemical evolution of primitive mafic arc magmas and their open-system interactions in the upper crust during mid-Jurassic growth of the Sierra Nevada batholith (California). We present field and petrographic observations and mineral and whole-rock chemistry of the ~35-km 2 ultramafic to dioritic Emigrant Gap mafic complex and an adjacent penecontemporaneous ~90-km 2 granodiorite that together make up the composite pluton. In the Emigrant Gap mafic complex, four roughly central masses of dunite, wehrlite, and olivine clinopyroxenite are surrounded by weakly layered gabbronorite and non-layered diorite. The ultramafic rocks are cumulates formed from near-liquidus minerals of primitive arc magmas that accumulated in steep feeder zones with substantial modification by melt–mush reaction as primitive liquids repeatedly transited the mush-filled conduits. The dominant gabbronoritic rocks are the variably accumulative products of more advanced crystallization–differentiation of arc tholeiitic basalts and basaltic andesites. The adjacent granodiorite intrusion originated separately and preserves field and geochemical evidence for assimilation of metasedimentary rocks. Open-system hybridization between the gabbronoritic mushes and the granodioritic magma produced an intervening body of two-pyroxene diorite. We infer that the ultramafic rocks and gabbronorite of the Emigrant Gap mafic complex crystallized from near-primitive arc basaltic to basaltic andesitic magmas at ~0.15–0.3 GPa, with estimated f O 2 of ≥FMQ +1 and dissolved H 2 O concentrations of only ~0.5–2 wt %. Notably, the Emigrant Gap composite pluton is distinct from other Mesozoic plutons in the Sierra Nevada batholith because of (1) its abundance of mafic and ultramafic rocks that crystallized from relatively primitive mafic melts and (2) the low inferred H 2 O concentrations of its parental magmas, indicated by a near absence of igneous amphibole and by the intermediate rather than calcic compositions of plagioclase. A Jurassic regional extension event probably accounts for the formation of relatively dry primitive arc magmas, as well as for their ascending to the upper crust.

California↗
USGS1995

Methods for estimating the magnitude and frequency of peak discharges of rural, unregulated streams in Virginia

Methods are presented for estimating the peak discharges of rural, unregulated streams in Virginia. A Pearson Type III distribution is fitted to the logarithms of the unregulated annual peak-discharge records from 363 stream-gaging stations in Virginia to estimate the peak discharge at these stations for recurrence intervals of 2 to 500 years. Peak-discharge characteristics for 284 unregulated stations are divided into eight regions based on physiographic province, and regressed on basin characteristics, including drainage area, main channel length, main channel slope, mean basin elevation, percentage of forest cover, mean annual precipitation, and maximum rainfall intensity. Regression equations for each region are computed by use of the generalized least-squares method, which accounts for spatial and temporal correlation between nearby gaging stations. This regression technique weights the significance of each station to the regional equation based on the length of records collected at each cation, the correlation between annual peak discharges among the stations, and the standard deviation of the annual peak discharge for each station. Drainage area proved to be the only significant explanatory variable in four regions, while other regions have as many as three significant variables. Standard errors of the regression equations range from 30 to 80 percent. Alternate equations using drainage area only are provided for the five regions with more than one significant explanatory variable. Methods and sample computations are provided to estimate peak discharges at gaged and engaged sites in Virginia for recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, and to adjust the regression estimates for sites on gaged streams where nearby gaging-station records are available.

Virginia↗
USGS2026

Background and geologic model for the 2024 U.S. Geological Survey assessment of undiscovered conventional petroleum resources in the Norphlet Formation, U.S. Gulf Coast

The Upper Jurassic Norphlet Formation is a stratigraphic unit located in the subsurface of the United States Gulf coastal plain and offshore Gulf of America (Gulf of Mexico). The Norphlet consists of clastic sediments derived from the southern Appalachian highlands, and was deposited in a continental dryland setting on the margins of the early Gulf of America rift basin. The formation contains a variety of sedimentary facies, including facies representing updip alluvial fans, wadis, red beds, and a widespread erg, or eolian sand sea. The erg facies form the primary reservoirs from which hydrocarbons sourced from the overlying Smackover Formation are produced. Development of the Norphlet began in earnest in the late 1960s, with new discoveries, including large fields at Flomaton, Hatters Pond, and Mobile Bay, continuing through the 2000s–2010s with the discovery of major oil fields in the offshore Desoto Canyon area. Onshore, Norphlet hydrocarbon traps typically are faulted salt anticlines, although many fields also have a stratigraphic trapping component because of the preservation of dune-scale bedforms. Thermal maturity across the Norphlet play varies significantly with depth, resulting in oil production transitioning downdip to natural gas liquids production, then dry gas production; production depths reach to greater than 21,000 feet. A 2024 U.S. Geological Survey assessment divided the Norphlet into four conventional assessment units defined by the type and frequency of trapping mechanisms, sedimentology, and thermal maturity of each area, with little potential for hydrocarbon development observed west of the Mississippi River. Although petroleum resources in the Norphlet were revised downward, it is probable that significant resources remain to be discovered.

Alabama, Florida, Mississippi↗
USGS2024

A decade of shaking in the Garden City: The dynamics of preparedness, perceptions, and beliefs in Canterbury, New Zealand, and implications for earthquake information

This study explored earthquake preparedness over time - before, during, and 10 years after the Canterbury Earthquake Sequence (CES) in Aotearoa New Zealand (NZ). Surveys of Canterbury residents were conducted in 2009, 2013 and again in 2021, using variables derived from Community Engagement Theory (CET). The surveys measured earthquake perceptions and beliefs, participation and engagement, and preparedness actions. Results were compared across the three samples. Findings indicate that perceptions and beliefs (e.g. risk perception, outcome expectancy beliefs), and types of preparedness actions taken (e.g. collection of survival items, structural preparedness, community and agency relationships), differed over time, depending on people’s experiences before, during, and after the CES. For example, during and after the CES people were more likely to believe that preparing provided a benefit to daily life, but less likely to think it could reduce property damage, perhaps due to people’s experiences of disruption and damage during the earthquakes. An understanding of such dynamics can assist with the provision and timing of risk and preparedness information. This study highlights the importance of providing applicable and actionable preparedness information, that is relevant to people’s experiences, throughout an earthquake sequence. Such information might evolve and change in focus over time depending on risks and needs. Focus could also be given to information that builds peoples beliefs and capacities to undertake preparedness in evolving situations. Understanding preparedness in the context of different experiences and timeframes is useful in helping update models such as the CET, where the dynamics of time might be better incorporated.

Canterbury↗
USGS1976

The water quality of Livingston Reservoir on the Trinity River, Southeastern Texas

The concentrations of dissolved solids, chloride, and sulfate in Livingstone Reservoir on the Trinity River in southeastern Texas usually average less than 250 mg/1 (milligrams per litre), 40 mg/1, and 50 mg/1, respectively. The water is usually hard or moderately hard (61 to 180 mg/1 as calcium carbonate). The concentrations of principal dissolved constituents in the reservoir are usually maximum during summer and fall when evaporation is high and inflow is low. Thermal stratification of the reservoir usually begins in March and persists until September or October. Neither the seasonal variation of dissolved constituents in inflow to the reservoir nor thermal stratification has resulted in significant stratification of the principal dissolved constituents. However, thermal stratification has resulted in significant seasonal and areal variations of dissolved oxygen, which results in higher concentration of dissolved iron, dissolved manganese, total phosphorus, and total inorganic nitrogen.

Texas↗
USGS2024

Numerical model of the groundwater-flow system near the southeastern part of Puget Sound, Washington

Groundwater flow in the active model area (AMA) was simulated using a groundwater-flow model. A steady-state model version of the model simulates equilibrium conditions, and a transient model version simulates monthly variability. The model corresponds to the physical and temporal dimensions of the conceptual model and groundwater budget. The steady-state model version represents average conditions for an 11-year period (January 1, 2005–December 31, 2015), and the transient model represents monthly hydrologic variability within that period. The 13-layer model was constructed using MODFLOW-NWT with a uniformly spaced grid consisting of 416 rows, 433 columns, and cells with a horizontal dimension of 500 feet (ft) on a side. The model was calibrated to measured values of water levels in wells and lakes and estimated base flow for selected streamflow measurement stations, commonly referred to as streamgages. Model calibration was accomplished using a combination of manual and automatic methods, including the Model-Independent Parameter Estimation (PEST) program that adjusted model input parameters with the aim of minimizing the difference between estimated and model-simulated values of hydraulic head and base flow. Model boundary conditions consist of all simulated groundwater inflow to and outflow from the AMA. For example, a stream reach that simulates a gain from or loss to groundwater is a boundary condition that allows water to exit or enter, respectively, the groundwater system. Other boundary conditions include springs, seeps, precipitation recharge, groundwater exchange with lakes and Puget Sound, and groundwater pumping. A comparison of the estimated groundwater budget to that simulated by the steady-state model version indicates that the relative percentages of total inflow or total outflow for six major categories of boundary conditions are similar for the two budgets. The model was used to simulate three suites of scenarios of potential drought and water-use changes. Scenario 1 suite consisted of the steady-state model version that was run with 0, 15, 20, and 25 percent reduction of precipitation recharge to assess the corresponding reductions in base flow with decreasing recharge. The last simulation for the scenario 1 suite consisted of the transient model version simulating 3 years of consecutive seasonal drought, defined by the months of May through September, to assess the corresponding base-flow reductions. Scenario 2 suite consisted of the steady-state model version with all simulated groundwater use removed, compared with a simulation that includes current groundwater use to evaluate changes to potentiometric surfaces and base flows. Scenario 3 suite consisted of a transient model version of the model that simulated pumping increases for four different categories of water-supply wells (compared to no pumping increases) to evaluate resulting reductions in base flow. Although, these scenarios provide examples of model applications and useful insights, many other scenarios could be simulated. A description of how to download the model is described in the body of this report. Uncertainty is associated with most model inputs. Groundwater levels, lake levels, and land-surface altitudes are relatively certain; other model inputs are far less certain, including precipitation recharge, base flow, hydraulic properties, water use, and the three-dimensional structure of subsurface hydrogeologic units. Models are useful not because of high levels of accuracy of all model inputs, but because they combine the best information and estimates available, thereby providing the best predictions available related to physical processes. The model described in this report simulates groundwater flow on a regional scale, which has inherent limitations for simulating hydrologic scenarios at local scales. Model structures and inputs were generalized to be consistent with this regional scale. For example, the actual groundwater system has much greater heterogeneity of hydraulic conductivity than is possible within the model’s degrees of freedom. Variations in hydraulic gradients over distances less than 500 ft cannot be simulated. The distances between model features, such as a pumping well and a stream, must be placed at 500-ft intervals and are co-located if both features are within the same model cell.

Washington↗
USGS2025

Geologic models underpinning the 2024 U.S. Geological Survey assessment of undiscovered oil and gas resources in the Hosston and Travis Peak Formations of the onshore Gulf Coast region, U.S.A.

The Early Cretaceous (Berriasian–Hauterivian) Hosston Formation in Louisiana and eastward is time correlative to the Travis Peak Formation of Texas and southern Arkansas. The formation is a first-order clastic sequence with a regional carbonate transgressive surface as an upper contact. The Hosston and Travis Peak formations contain conventional natural gas and oil accumulations that have been produced for nearly a century. These mature reservoirs contain terrigenous fluvial-deltaic, shore-zone, and paralic deposits across the productive trend; organic-lean mudstone and siltstone lithologies are found outboard of the Lower Cretaceous shelf margin. Producing reservoirs exhibit normal pressure gradients from 0.43 to 0.55 psi/ft (9.7 to 12.4 kpa/m), depths from 4000 to over 20,000 ft (1220 to 6100 m), and temperatures from 150 to 385°F (65 to 196°C). Wells are primarily vertical completions. The number of new field wildcats has been declining since the late 1990s. This paper presents comprehensive geologic models, which include lithofacies maps, structure and isopach maps, burial history models, regional seismic interpretations, and events charts that underpin the recently completed U.S. Geological Survey assessment of undiscovered, technically recoverable hydrocarbons within the Hosston and Travis Peak formations. This study also provides geographic and stratigraphic distributions of Hosston–Travis Peak reservoir properties, including geopressure, reservoir temperature, porosity, permeability, API gravity, and gas-oil ratios. Results indicate estimated undiscovered, technically recoverable mean resources of 28 million barrels of oil and 35.8 trillion cubic ft of gas in conventional and continuous accumulations within the Lower Cretaceous Hosston and Travis Peak formations of the onshore U.S. Gulf Coast region. Quantitative assessment results are detailed in U.S. Geological Survey Fact Sheet 2025–3021 and associated Data Release.

GCAGS Journal↗
USGS2025

Origin and evolution of mafic volcanism associated with 3 m.y. of andesite production at the Goat Rocks volcanic cluster, southern Washington Cascade Range

More than 3 m.y. of mafic volcanism near the Goat Rocks volcanic cluster in the southern Washington Cascade Range, USA, lends insight into the evolution of basalts and the subarc mantle at a long-lived, major arc volcanic locus. We contribute field observations, 40 Ar/ 39 Ar dates, paleomagnetic directions, and bulk rock and mineral compositions to characterize nine mafic units that erupted in association with the Goat Rocks volcanic cluster. The time frame of mafic volcanism, ca. 3.6 Ma to 60 ka, encompasses the lifespan of the central volcanic cluster (3.1 Ma to 115 ka), with a lull from ca. 2.7 Ma to 1.4 Ma. A climactic period of voluminous mafic activity and far-traveled lava flows, including construction of the Hogback Mountain shield volcano, coincided with voluminous andesite eruptions from the central volcanic cluster. The basaltic rocks in the Goat Rocks area are calc-alkaline to barely tholeiitic and have high field strength element depletion relative to large-ion lithophile elements characteristic of calc-alkaline basalts (CAB) of the Cascade volcanic arc. Unlike at neighboring andesitic volcanic centers (Mounts Adams, St. Helens, and Rainier), no other mafic end members such as high-aluminum olivine tholeiite (HAOT) or intraplate-type basalt (IPB) are present at or near the Goat Rocks volcanic cluster, although some of the calc-alkaline basalts in this study have IPB-like affinities. The Goat Rocks mafic units exhibit two main temporal trends in composition: (1) the most primitive basalts erupted earlier, compared to less primitive and more evolved compositions later, and (2) high field strength element concentrations are higher in the younger basalt units relative to the oldest two. In contrast to these temporal trends, the mafic units define two compositional groups that recur through time, a low-Sr and a high-Sr group, each with distinct trace element and Sr and Nd isotope ratios. Although radiogenic isotope ratios are generally aligned with High Cascades CAB and HAOT, some extend toward IPB of Mount Adams and Simcoe Mountains volcanic field. Olivine-dominated crystal fractionation at shallow pressure from a small range of parent magma compositions accounts for much of the variation among the basalts and basaltic andesites. A high-pressure fractionation model is plausible for only one of the youngest basalt units (basalt of Walupt Lake volcano). Mafic recharge and crustal assimilation accounts for the incompatible-element enriched composition of basaltic andesites erupted during construction of the largest andesitic centers, further supporting sustained basalt mass flux and thermal energy driving andesite genesis. We model the most primitive members of the Goat Rocks mafic units as partial melts of successively less depleted mantle in time. Variable degrees of fluxing with fluids and melts from subduction explain the distinction between high-Sr and low-Sr groups. We propose that mantle metasomatism by ancestral subduction and fluid-flux melting is heterogeneously distributed through the local subarc mantle and played a greater role in the genesis of the high-Sr basalt group. The limited range of primitive basalt types around the Goat Rocks volcanic cluster contrasts with the much greater diversity of basalts throughout the southern Washington to northern Oregon Cascade arc. On the other hand, the central volcanic cluster encompasses nearly the entire diversity observed at neighboring composite volcanoes. In the case of the Goat Rocks area at least, and perhaps attributable to the entire region, this means that the genesis of diverse intermediate magmas is independent from and does not require vastly different parental basalt compositions.

California, Oregon, Washington↗
USGS2026

Earthquake catalog for the Fairbanks region of central Alaska, 2014–2024, based on waveform cross-correlation

The Fairbanks region of central Alaska is part of a broad zone of intraplate crustal deformation, situated north of the Denali fault and north of the ongoing collision and flat‐slab subduction of the Yakutat oceanic plateau. Seismicity in the Fairbanks region occurs both in diffuse areas as well as in well‐defined lineaments, such as the left‐lateral Salcha fault, which hosted the 1937 M 8 7.3 earthquake. Starting with the regional seismicity catalog, we perform waveform cross‐correlation, network‐matched filtering, and relative relocation to obtain an enhanced seismicity catalog over the time period 2014–2024. Based on the relocated catalog, we interpret a set of 15 fault segments, including two conjugate faults and two new faults east of the previously documented fault system. Considering the combined seismicity in the Minto and Fairbanks regions, the median depth of seismicity decreases from east (6 km) to west (20 km). Our interpreted faults provide guidance for future tectonic modeling and assessment of seismic hazards in this region.

Alaska↗
USGS2025

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California↗
USGS2023

Eruptive history within the vicinity of Al Madīnah in northern Harrat Rahat, Kingdom of Saudi Arabia

The northernmost part of the Harrat Rahat volcanic field contains early Pleistocene to Holocene mafic eruptive products within the vicinity of the city of Al Madīnah, Kingdom of Saudi Arabia. A detailed geologic investigation into the eruptive history of a 570 square kilometer (km 2 ) area covering Al Madīnah and the surrounding area has yielded 33 mapped Quaternary volcanic units consisting of lava flows, scoria cones, and shield volcanoes. These eruptive products consist of continental, intraplate alkalic and minor transitional basalts, hawaiites, and a single mugearite that were emplaced from at least 1,014±14 thousand years ago (ka) to a single Holocene eruption in 1256 C.E. Lava flows are generally 10 to 15 kilometers (km) long (but can reach 23 km long), 1 to 3 km wide, and at least 10 meters thick. Most of the mapped units erupted episodically between 400 and 340 ka and 180 and 100 ka. Despite small individual volumes (less than 1 cubic kilometers dense rock equivalent), each unit represents eruption of a distinct magma batch that was strongly influenced by clinopyroxene, olivine, and plagioclase fractionation. Some of these units are interpreted to have undergone magma mixing pre- and (or) syneruptively. Integrating eruption ages, geochemistry, and paleomagnetic data yields evidence that some eruptions were temporally and (or) spatially clustered. Aligned scoria cones and elongate vent edifices were constructed atop fissure vent systems that reflect the local stress field, which controls dike ascent through the middle and upper crust.

Professional Paper↗
USGS2025

Integrated stratigraphic and geochemical analysis of organic-rich intervals of the Lewis Shale in the eastern Washakie Basin, Wyoming

Geological studies in the Cretaceous Western Interior Seaway (KWIS) in North America reveal highly variable sedimentological conditions on its western shore caused by rapidly changing sea level and detrital input during the seaway’s closure. Here we examine a 601-foot (183 meters) continuous core through the lower part of the Maastrichtian Lewis Shale in the eastern Washakie Basin, Wyoming, through integrating stratigraphic and geochemical analyses to better understand organic matter deposition and preservation during the final marine transgression within the seaway. The core penetrates eight organic-rich flooding intervals (F1-F8) and a regional condensed section, informally known as the Asquith marker. The lower portions of the core record sediment input from a southern source, likely the Sierra Madre/Park Ranges, while the upper part records sediment input from a northern source, likely the Granite Mountains. This provenance transition is supported by regional stratigraphic analysis and changes in bulk mineral and trace metal composition. The Asquith marker and early flooding surfaces are enriched in oil-prone, marine organic matter deposited under dysoxic to euxinic conditions, whereas younger flooding surfaces show increased terrigenous input and poorer preservation conditions. The Asquith marker is identified as a prime oil-prone source rock. Although younger flooding surfaces also exhibit favorable source-rock properties, their generative potential is reduced due to increased clastic and terrigenous organic matter deposition caused by regression. These results emphasize the importance of sediment source variability, organic matter preservation, and changing redox conditions to provide insight into the sediment provenance and petroleum potential of the Lewis Shale in the context of the final transgression within the KWIS.

Colorado, Wyoming↗
USGS2000

Assessment of potential for natural attenuation of chlorinated ethenes and ethanes in ground water at a petrochemical reclamation site, Harris County, Texas

Redox conditions in the Numerous Sand Channels Zone beneath a petrochemical reclamation site in Harris County, Texas, range from sulfate reducing to methanogenic as indicated by the presence of methane in ground water and the range of molecular hydrogen concentrations. Assessment of the potential for reductive dechlorination using BIOCHLOR as a screening tool indicated conditions favoring anaerobic degradation of chlorinated organic compounds in the Numerous Sand Channels Zone. Evidence supporting reductive dechlorination includes apparently biogenic cis -1,2-dichloroethene; an increased ratio of 1,2-dichloroethane to 1,1,2-trichloroethane downgradient from the assumed contaminant source area; ethene and methane concentrations greater than background concentrations within the area of the contaminant plume; and a positive correlation of the ratio of ethene to vinyl chloride as a function of methane concentrations. The body of evidence presented in this report argues for hydrogenolysis of trichloroethene to cis -1,2-dichloroethene; of 1,1,2-trichloroethane to 1,2-dichloroethane; and of vinyl chloride to ethene within the Numerous Sand Channels Zone. Simulations using BIOCHLOR yielded apparent first-order decay constants for reductive dechlorination in the sequence tetrachloroethene trichloroethene cis -1,2-dichloroethene vinyl chloride ethene within the range of literature values reported for each compound and apparent first-order decay constants for reductive dechlorination in the sequence 1,1,2-trichloroethane 1,2-dichloroethane slightly greater than literature values reported for each compound along the upgradient segment of a simulated ground-water flowpath. Except for vinyl chloride, apparent rates of reductive dechlorination for all simulated species show a marked decrease along the downgradient segment of the simulated ground-water flowpath. Evidence for reductive dechlorination of chlorinated ethenes within the Numerous Sand Channels Zone indicates potential for natural attenuation of chlorinated ethenes. Reductive dechlorination of chlorinated ethanes apparently occurs to a lesser extent, indicating relatively less potential for natural attenuation of chlorinated ethanes. Additional data are needed on the concentrations and distribution of chlorinated ethenes and ethanes in individual fine sand intervals of the Numerous Sand Channels Zone. This information, combined with lower minimum reporting levels for future chloroethane analyses, might enable a more complete and quantitative assessment of the potential for natural attenuation at the site.

Texas↗