Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Computers and Fluids”

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

Simulation of the effects of different inflows on hydrologic conditions in Lake Houston with a three-dimensional hydrodynamic model, Houston, Texas, 2009–10

Lake Houston, an important water resource for the Houston, Texas, area, receives inflows from seven major tributaries that compose the San Jacinto River Basin upstream from the reservoir. The effects of different inflows from the watersheds drained by these tributaries on the residence time of water in Lake Houston and closely associated physical and chemical properties including lake elevation, salinity, and water temperature are not well known. Accordingly, the U.S. Geological Survey (USGS), in cooperation with the City of Houston, developed a three-dimensional hydrodynamic model of Lake Houston as a tool for evaluating the effects of different inflows on residence time of water in the lake and associated physical and chemical properties. The Environmental Fluid Dynamics Code (EFDC), a grid-based, surface-water modeling package for simulating three-dimensional circulation, mass transport, sediments, and biogeochemical processes, was used to develop the model of Lake Houston. The Lake Houston EFDC model was developed and calibrated by using 2009 data and verified by using 2010 data. Three statistics (mean error, root mean square error, and the Nash-Sutcliffe model efficiency coefficient) were used to evaluate how well the Lake Houston EFDC model simulated lake elevation, salinity, and water temperature. The residence time of water in reservoirs is associated with various physical and chemical properties (including lake elevation, salinity, and water temperature). Simulated and measured lake-elevation values were compared at USGS reservoir station 08072000 Lake Houston near Sheldon, Tex. The accuracy of simulated salinity and water temperature values was assessed by using the salinity (computed from measured specific conductance) and water temperature at two USGS monitoring stations: 295826095082200 Lake Houston south Union Pacific Railroad Bridge near Houston, Tex., and 295554095093401 Lake Houston at mouth of Jack’s Ditch near Houston, Tex. Specific conductance and water temperature were measured at as many as four different depths at each of the two monitoring stations during 2009 and then used for assessing the accuracy of simulated values of salinity and water temperature during 2010. The performance evaluation statistics indicate that the model performed satisfactorily. The calibrated model was used to simulate two possible inflow scenarios to evaluate the changes in the residence time of water in Lake Houston. The two scenarios tested were an increased inflow of approximately 300 cubic feet per second for 1 month (May 2010) from two watersheds: the West Fork San Jacinto River and Luce Bayou. These scenarios were chosen to mimic the effects of possible small releases or diversions of water from outside the San Jacinto River Basin into the basin (or directly into the lake) on the residence time of water in Lake Houston. During the time of increased inflow for the two scenarios tested, maximum residence time decreased slightly from approximately 106 to 97 days.

Texas↗
USGS1994

Experimental studies of alunite: I. 18 O- 16 O and D-H fractionation factors between alunite and water at 250-450°C

We have determined oxygen and hydrogen isotope fractionation factors between alunite and water over a temperature range of 250-450??C by reacting synthetic natroalunite with 0.7 m K2SO4 -0.1 to 0.65 m H2SO4 solutions to produce K-rich alunite. From 88 to 95% alkali and isotope exchange were observed in most of these experiments, and the partial equilibrium method was used to compute equilibrium fractionation factors. Least-squares fits of the data give 103 In ??alunite(so4)-H2O = 3.09 ( 106 T2 (K)) - 2.94 and 103 In ??alunite(OH)-H2O = 2.28 ( 106 T2 (K)) - 3.90. The intramineral 18O- 16O fractionation factor 103 In ??alunite(so4-OH site) is given by the expression 0.8 ( 106 T2 (K)) + 0.96. The alunite-water D-H fractionation factor ranges from -19 at 450??C to -6 at 250??C and does not appear to be strongly dependent on temperature. Runs with alkali exchange in the opposite direction were used to obtain 18O- 16O and D-H fractionation factors between natroalunite (mol% Na = 70-75) and water at 350-450??C. These indicate that mol% Na has negligible effect on the fractionation factors over this temperature range. Measured 18O-16O and D-H fractionation factors between alunite and 1.0 m KCl -0.5 m H2SO4 fluids also agree within 2?? with the values obtained from the K2SO4-H2SO4 fluids. However, experiments with alunite and distilled water at 400??C gave a value of 103 In ??alunite(SO4)-H2O of 0.0, compared with a value of 3.9 obtained at this temperature with K2SO4- and H2SO4-bearing fluids. This suggests that changes in fluid composition can affect alunite-water 18O-16O fractionation factors. Reconnaissance experiments with fine-grained natural natroalunite demonstrate that alunite-water D-H exchange can occur by hydrogen diffusion, although this process is generally not significant in the experiments with coarser grained synthetic alunites. ?? 1994.

Geochimica et Cosmochimica Acta↗
USGS2026

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

Tennessee↗
USGS1986

Color images of Kansas subsurface geology from well logs

Modern wireline log combinations give highly diagnostic information that goes beyond the basic shale content, pore volume, and fluid saturation of older logs. Pattern recognition of geology from logs is made conventionally through either the examination of log overlays or log crossplots. Both methods can be combined through the use of color as a medium of information by setting the three color primaries of blue, green, and red light as axes of three dimensional color space. Multiple log readings of zones are rendered as composite color mixtures which, when plotted sequentially with depth, show lithological successions in a striking manner. The method is extremely simple to program and display on a color monitor. Illustrative examples are described from the Kansas subsurface. ?? 1986.

Computers & Geosciences↗
USGS2003

Finite element model predictions of static deformation from dislocation sources in a subduction zone: Sensitivities to homogeneous, isotropic, Poisson-solid, and half-space assumptions

Dislocation models can simulate static deformation caused by slip along a fault. These models usually take the form of a dislocation embedded in a homogeneous, isotropic, Poisson-solid half-space (HIPSHS). However, the widely accepted HIPSHS assumptions poorly approximate subduction zone systems of converging oceanic and continental crust. This study uses three-dimensional finite element models (FEMs) that allow for any combination (including none) of the HIPSHS assumptions to compute synthetic Green's functions for displacement. Using the 1995 Mw = 8.0 Jalisco-Colima, Mexico, subduction zone earthquake and associated measurements from a nearby GPS array as an example, FEM-generated synthetic Green's functions are combined with standard linear inverse methods to estimate dislocation distributions along the subduction interface. Loading a forward HIPSHS model with dislocation distributions, estimated from FEMs that sequentially relax the HIPSHS assumptions, yields the sensitivity of predicted displacements to each of the HIPSHS assumptions. For the subduction zone models tested and the specific field situation considered, sensitivities to the individual Poisson-solid, isotropy, and homogeneity assumptions can be substantially greater than GPS. measurement uncertainties. Forward modeling quantifies stress coupling between the Mw = 8.0 earthquake and a nearby Mw = 6.3 earthquake that occurred 63 days later. Coulomb stress changes predicted from static HIPSHS models cannot account for the 63-day lag time between events. Alternatively, an FEM that includes a poroelastic oceanic crust, which allows for postseismic pore fluid pressure recovery, can account for the lag time. The pore fluid pressure recovery rate puts an upper limit of 10-17 m2 on the bulk permeability of the oceanic crust. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS1984

The evolution of the southern California uplift, 1955 through 1976

The southern California uplift culminated in 1974 as a 150- km-wide crustal swell that extended about 600 km eastward and east-southeastward from Point Arguello to the Colorado River and Salton Sea, respectively; it was characterized by remarkably uniform height changes between 1959 and 1974 of 0.30-0.35 m over at least half of its 60,000-70,000 km 2 area. At its zenith, the uplift included virtually the entire Transverse Ranges geologic province and parts of the Coast Ranges, San Joaquin Valley, Sierra Nevada, Basin and Range, Mojave Desert, Peninsular Ranges, and Salton Trough provinces. The alinement of the western part of the uplift closely paralleled the east-trending Transverse Ranges, whereas the southern flank of the eastern lobe roughly coincided with the west-northwest-trending San Andreas fault. The position and configuration of the uplift associate it with a singularly complex section of the boundary between the North American and Pacific plates that has certainly sustained major modification during the past 5 million years and probably during the past 1 million years. Surface deformation can be categorized as tectonic or nontectonic. Nontectonic vertical displacements associated with the activities of man have overwhelmed natural compaction and areally significant soil expansion in the southern California area. Because tectonic displacements are implicitly defined as those that cannot be otherwise explained, those vertical movements that can be reasonably attributed to artificial processes have been subtracted from our reconstructed configurations of the uplift. Hence this reconstruction has necessarily included the assembly and evaluation of an enormous volume of data on oil-field operations, changes in ground-water levels, and measured subsidence (or rebound) associated with changes in the underground fluid regimen. Measured changes in height at various stages in the evolution of the uplift have been based chiefly on first-order levelings carried out between 1953 and 1976. Exceptions to this generalization consist largely of the results of pre-1953 surveys through the western Transverse Ranges and the eastern Mojave Desert. Errors in measured height differences derive from blunders, systematic survey errors, random survey errors, improperly formulated orthometric corrections, and intrasurvey movement; the last of these has created the most serious problems encountered in our reconstruction of the basic data. A variety of independent tests indicate that survey error associated with the utilized levelings was generally small and fell largely within the predicted random-error range. Moreover, the redundancy and coherence displayed by the entire data set provide convincing evidence of survey accuracy and the virtual absence of height- and slope-dependent error in particular. Our reconstructions of the changing configuration of the uplift derive chiefly from comparisons among sequentially developed observed elevations along the same route. Most of the observed elevations from which the vertical displacements were computed have been reconstructed with respect to bench mark Tidal 8, San Pedro, as invariant in height. Because the San Pedro tide station has been characterized by a history of modest relative uplift, vertical displacements referred to this station are biased slightly toward the appearance of subsidence. Where the observed elevations cannot be conveniently tied to Tidal 8, they have been referred to secondary control points whose history with respect to Tidal 8 can be independently established. Each of the lines of observed elevation changes provides, accordingly, a section athwart or along the axis of the uplift from which the changes in the configuration of the uplift can be roughly generalized. Because relatively few surveys were run in 1955, which we choose as a representative temporal datum, we have commonly incorporated the results of earlier or of somewhat later levelings as the equivalents of 1955 surveys. Although this procedure introduces a certain subjectivity, the probable equivalence between the results of these earlier or later surveys with those that would have been obtained had this leveling been carried out in 1955, usually can be independently tested. Wherever the calculated vertical displacements are based on comparisons between the results of levelings over different routes, the observed elevations have been orthometrically corrected to agree with those that would have been produced had each of these surveys been along the same route. The growth of the southern California uplift consisted of two well-defined spasms of positive movement, the second of which was closely followed by partial collapse. Our reconstruction, although it clearly errs in detail, indicates that the uplift, together with marginal and apparently ephemeral tectonic subsidence, nucleated in the west-central Transverse Ranges near Ozena, sometime between the spring of 1959 and the spring of 1960. The uplift expanded rapidly eastward (and probably westward as well), and by the fall of 1961 much of the Transverse Ranges and the Mojave Desert at least as far east as Twentynine Palms had risen by as much as 0.25 m. Between 1962 and 1972 the area included by the initially developed (1959-61) uplift sustained additional but clearly decelerating uplift accompanied locally by oscillatory displacements. Between 1972/ 73 and 1974 a second crustal spasm extended the uplift eastward to the Colorado River and elevated much of the eastern Mojave Desert by values that equaled or exceeded those developed within the western lobe. Between 1974 and 1976, at least the central part of the uplift sustained partial collapse that nowhere amounted to less than 50 percent of the cumulative uplift since 1959. Whether this collapse affected the entire uplift is conjectural, but we now recognize well-defined evidence of major down-to-the-north tilting that must have occurred within the eastern part of the uplift at some time between 1974 and 1976. Accumulating evidence indicates that nearly all the area included with the southern California uplift underwent similar uplift and partial collapse during the early part of the 20th century. Thus we infer that the recent uplift represents but a single event in an ongoing, more or less cyclic deformational process characterized by a period of about 50 years. Even though less than two full cycles are expressed in the geodetic record, the cumulative rate of uplift near the center of the recent uplift probably has averaged about 5 mm/yr, a value that is roughly consistent with the uplift rates that have been deduced for the late Quaternary emergent marine terraces along the south flank of the Transverse Ranges. Although the evolution of the recent uplift is relatively well defined, its correlation with the regional seismicity is poorly defined. A comparison between the occurrence of southern California earthquakes of magnitude ≥4 during the period 1932 to 1976 with the 1974 configuration of the uplift demonstrates the existence of (1) relatively aseismic areas within the western lobe of the uplift (in the western Transverse Ranges), in the central part of the uplift (in the western Mojave Desert), and along an east-trending zone that extends into the eastern Mojave athwart the south flank of the uplift (north of the Salton Sea) and (2) localized concentrations of seismic activity along the flanks of the uplift. Moreover, 9 of the 10 largest earthquakes recorded within or around the area of the southern California uplift during the period 1932 to 1976 (the 1933 Long Beach, the 1941 Santa Barbara, the 1946 Walker Pass, the 1947 Manix, the 1948 Desert Hot Springs, and the four major 1952 Kern County shocks) occurred before the inception of the uplift in 1959 or 1960. The area embraced by the southern California uplift has been identified with geodetically defined horizontal strain, part of which may have accumulated as a major north-south contractional event that roughly coincided with the first spasm of uplift. Nonetheless, continuing contractional strain associated with regionally developed partial collapse argues that the uplift cannot be fully explained simply as the vertical expression of continuing north-south compression. Consideration of the two well-defined historical episodes of uplift and partial collapse indicate that the southern California uplift may be the product of decoupling and viscous flow beneath the seismogenic zone, presumably driven by continuing motion between the irregularly margined plates south of the great bend of the San Andreas fault. Because the magnitude of the maximum uplift associated with each episode was approximately the same, there may be some threshold value above which collapse (viscous flow) may ensue; the absence of total collapse may be a function of precollapse strain hardening within the postulated subseismogenic viscoelastic layer.

California↗
USGS1998

Geohistory and thermal maturation in the Cherokee Basin (Mid-Continent, U.S.A.): results from modeling

The Cherokee basin in southeastern Kansas contains a stratigraphic section consisting mostly of Permian-Pennsylvanian alternating clastics and thin carbonates overlying carbonates of Mississippian and Cambrian-Ordovician age on a Precambrian crytalline basement. Based on a conceptual model of events of deposition, nondeposition, and erosion, a burial history model for (1) noncompaction, and a series of models for (2) compaction are computed for a borehole location in the south-central part of the basin. The models are copled with the calculation of nonsteady-state geothermal conditions. Maximum temperatures during basin evolution of about 70??C at the base of the organic-rich Pennsylvanian are predicted by our models, assuming pure heat conduction and a heat flow from the basement of 60 m W/m2. The maturation of organic matter as indicated by three different vitrinite reflectance (Ro) models is on the order og 0.3-0.5% Ro for Pennsylvanian rocks and 0.6% Ro for the Devonian-Mississippian Cattanooga Shale. Vitrinite reflectance was measured on subsurface smaples from three wells. The measured values correlate in the upper part of the sequence with modeled data, but diverge slightly in the Lower Pennsylvanian and Cattanooga Shale. The differences in maturation may be a result of differing local geological conditions within the basin. The relatively high Ro-depth gradients observed in one borehole may be explained by conditions in the Teeter oil field, which is a typical plains-type anticline that has been affected by fluid flow through vertical faults. Higher Ro values correlate positively with the grade of sulfidfe mineralization in the sediment, which may be a hint of fluid impact. The high Ro values relative to the shallow depth of the Mississippian and the Chattanooga Shale in the Brown well are on the order of Ro values modeled for the same stratigraphic units at present-day greater depths and may reflect uplift of the Ozark dome, located further east, affecting the eastern side of the Cherokee Basin.Based on a concept model of deposition, nondeposition and erosion, a burial history model for noncompaction, and a series of models for compaction are developed for a borehole location in a south-central part of the Cherokee basin in southeastern Kansas. Coupled with the calculation of nonsteady state-state geothermal conditions, the models predict maximum temperatures during evolution of about 70 ??C at the base of the organic-rich Pennsylvanian. A difference in organic matter maturation in the Pennsylvanian and the Chattanooga shale exhibited by vitrinite reflectance models indicate probably differing local geological conditions within the basin.

American Association of Petroleum Geologists Bulle↗
USGS1983

Geotherm: the U.S. geological survey geothermal information system

GEOTHERM is a comprehensive system of public databases and software used to store, locate, and evaluate information on the geology, geochemistry, and hydrology of geothermal systems. Three main databases address the general characteristics of geothermal wells and fields, and the chemical properties of geothermal fluids; the last database is currently the most active. System tasks are divided into four areas: (1) data acquisition and entry, involving data entry via word processors and magnetic tape; (2) quality assurance, including the criteria and standards handbook and front-end data-screening programs; (3) operation, involving database backups and information extraction; and (4) user assistance, preparation of such items as application programs, and a quarterly newsletter. The principal task of GEOTHERM is to provide information and research support for the conduct of national geothermal-resource assessments. The principal users of GEOTHERM are those involved with the Geothermal Research Program of the U.S. Geological Survey. Information in the system is available to the public on request. ?? 1983.

Computers & Geosciences↗
USGS2013

Colloid transport in saturated porous media: Elimination of attachment efficiency in a new colloid transport model

A colloid transport model is introduced that is conceptually simple yet captures the essential features of colloid transport and retention in saturated porous media when colloid retention is dominated by the secondary minimum because an electrostatic barrier inhibits substantial deposition in the primary minimum. This model is based on conventional colloid filtration theory (CFT) but eliminates the empirical concept of attachment efficiency. The colloid deposition rate is computed directly from CFT by assuming all predicted interceptions of colloids by collectors result in at least temporary deposition in the secondary minimum. Also, a new paradigm for colloid re-entrainment based on colloid population heterogeneity is introduced. To accomplish this, the initial colloid population is divided into two fractions. One fraction, by virtue of physiochemical characteristics (e.g., size and charge), will always be re-entrained after capture in a secondary minimum. The remaining fraction of colloids, again as a result of physiochemical characteristics, will be retained “irreversibly” when captured by a secondary minimum. Assuming the dispersion coefficient can be estimated from tracer behavior, this model has only two fitting parameters: (1) the fraction of the initial colloid population that will be retained “irreversibly” upon interception by a secondary minimum, and (2) the rate at which reversibly retained colloids leave the secondary minimum. These two parameters were correlated to the depth of the Derjaguin-Landau-Verwey-Overbeek (DLVO) secondary energy minimum and pore-water velocity, two physical forces that influence colloid transport. Given this correlation, the model serves as a heuristic tool for exploring the influence of physical parameters such as surface potential and fluid velocity on colloid transport.

Water Resources Research↗
USGS2011

A distribution-based parameterization for improved tomographic imaging of solute plumes

Difference geophysical tomography (e.g. radar, resistivity and seismic) is used increasingly for imaging fluid flow and mass transport associated with natural and engineered hydrologic phenomena, including tracer experiments, in situ remediation and aquifer storage and recovery. Tomographic data are collected over time, inverted and differenced against a background image to produce ‘snapshots’ revealing changes to the system; these snapshots readily provide qualitative information on the location and morphology of plumes of injected tracer, remedial amendment or stored water. In principle, geometric moments (i.e. total mass, centres of mass, spread, etc.) calculated from difference tomograms can provide further quantitative insight into the rates of advection, dispersion and mass transfer; however, recent work has shown that moments calculated from tomograms are commonly biased, as they are strongly affected by the subjective choice of regularization criteria. Conventional approaches to regularization (Tikhonov) and parametrization (image pixels) result in tomograms which are subject to artefacts such as smearing or pixel estimates taking on the sign opposite to that expected for the plume under study. Here, we demonstrate a novel parametrization for imaging plumes associated with hydrologic phenomena. Capitalizing on the mathematical analogy between moment-based descriptors of plumes and the moment-based parameters of probability distributions, we design an inverse problem that (1) is overdetermined and computationally efficient because the image is described by only a few parameters, (2) produces tomograms consistent with expected plume behaviour (e.g. changes of one sign relative to the background image), (3) yields parameter estimates that are readily interpreted for plume morphology and offer direct insight into hydrologic processes and (4) requires comparatively few data to achieve reasonable model estimates. We demonstrate the approach in a series of numerical examples based on straight-ray difference-attenuation radar monitoring of the transport of an ionic tracer, and show that the methodology outlined here is particularly effective when limited data are available.

Geophysical Journal International↗
USGS2014

A ternary age-mixing model to explain contaminant occurrence in a deep supply well

The age distribution of water from a public-supply well in a deep alluvial aquifer was estimated and used to help explain arsenic variability in the water. The age distribution was computed using a ternary mixing model that combines three lumped parameter models of advection-dispersion transport of environmental tracers, which represent relatively recent recharge (post- 1950s) containing volatile organic compounds (VOCs), old intermediate depth groundwater (about 6500 years) that was free of drinking-water contaminants, and very old, deep groundwater (more than 21,000 years) containing arsenic above the USEPA maximum contaminant level of 10 µg/L. The ternary mixing model was calibrated to tritium, chloroflorocarbon-113, and carbon-14 ( 14 C) concentrations that were measured in water samples collected on multiple occasions. Variability in atmospheric 14 C over the past 50,000 years was accounted for in the interpretation of 14 C as a tracer. Calibrated ternary models indicate the fraction of deep, very old groundwater entering the well varies substantially throughout the year and was highest following long periods of nonoperation or infrequent operation, which occurred during the winter season when water demand was low. The fraction of young water entering the well was about 11% during the summer when pumping peaked to meet water demand and about 3% to 6% during the winter months. This paper demonstrates how collection of multiple tracers can be used in combination with simplified models of fluid flow to estimate the age distribution and thus fraction of contaminated groundwater reaching a supply well under different pumping conditions.

Groundwater↗
USGS2001

Differential rates of feldspar weathering in granitic regoliths

Differential rates of plagioclase and K-feldspar weathering commonly observed in bedrock and soil environments are examined in terms of chemical kinetic and solubility controls and hydrologic permeability. For the Panola regolith, in the Georgia Piedmont Province of southeastern United States, petrographic observations, coupled with elemental balances and 87Sr/86Sr ratios, indicate that plagioclase is being converted to kaolinite at depths > 6 m in the granitic bedrock. K-feldspar remains pristine in the bedrock but subsequently weathers to kaolinite at the overlying saprolite. In contrast, both plagioclase and K-feldspar remain stable in granitic bedrocks elsewhere in Piedmont Province, such as Davis Run, Virginia, where feldspars weather concurrently in an overlying thick saprolite sequence. Kinetic rate constants, mineral surface areas, and secondary hydraulic conductivities are fitted to feldspar losses with depth in the Panola and Davis Run regoliths using a time-depth computer spreadsheet model. The primary hydraulic conductivities, describing the rates of meteoric water penetration into the pristine granites, are assumed to be equal to the propagation rates of weathering fronts, which, based on cosmogenic isotope dating, are 7 m/106 yr for the Panola regolith and 4 m/106 yr for the Davis Run regolith. Best fits in the calculations indicate that the kinetic rate constants for plagioclase in both regoliths are factors of two to three times faster than K-feldspar, which is in agreement with experimental findings. However, the range for plagioclase and K-feldspar rates (kr = 1.5 x 10-17 to 2.8 x 10-16 mol m-2 s-1) is three to four orders of magnitude lower than for that for experimental feldspar dissolution rates and are among the slowest yet recorded for natural feldspar weathering. Such slow rates are attributed to the relatively old geomorphic ages of the Panola and Davis Run regoliths, implying that mineral surface reactivity decreases significantly with time. Differential feldspar weathering in the low-permeability Panola bedrock environment is more dependent on relative feldspar solubilities than on differences in kinetic reaction rates. Such weathering is very sensitive to primary and secondary hydraulic conductivities (qp and qs), which control both the fluid volumes passing through the regolith and the thermodynamic saturation of the feldspars. Bedrock permeability is primarily intragranular and is created by internal weathering of networks of interconnected plagioclase phenocrysts. Saprolite permeability is principally intergranular and is the result of dissolution of silicate phases during isovolumetric weathering. A secondary to primary hydraulic conductivity ratio of qs/qp = 150 in the Panola bedrock results in kinetically controlled plagioclase dissolution but thermodynamically inhibited K-feldspar reaction. This result is in accord with calculated chemical saturation states for groundwater sampled in the Panola Granite. In contrast, greater secondary conductivities in the Davis Run saprolite, qs/qp = 800, produces both kinetically controlled plagioclase and K-feldspar dissolution. Faster plagioclase reaction, leading to bedrock weathering in the Panola Granite but not at Davis Run, is attributed to a higher anorthite component of the plagioclase and a wetter and warmer climate. In addition, the Panola Granite has an abnormally high content of disseminated calcite, the dissolution of which precedes the plagioclase weathering front, thus creating additional secondary permeability. Copyright ?? 2001 Elsevier Science Ltd.

Geochimica et Cosmochimica Acta↗
USGS1980

Transition of basaltic lava from pahoehoe to aa, Kilauea Volcano, Hawaii: Field observations and key factors

Nearly all Hawaiian basaltic lava erupts as pahoehoe, and some changes to aa during flowage and cooling; factors governing the transition involve certain critical relations between viscosity and rate of shear strain. If the lava slows, cools, and stops in direct response to concomitant increase in viscosity before these critical relations are reached, it remains pahoehoe. But, if flow mechanics (flow rate, flow dimensions, slope, momentum, etc.) impel the lava to continue to move and deform even after it has become highly viscous, the critical relations may be reached and the lava changes to aa. Typical modes of transition from pahoehoe to aa include: (1) spontaneous formation of relatively stiff clots in parts of the flowing lava where shear rate is highest; these clots grow into discrete, rough, sticky masses to which the remaining fluid lava incrementally adheres; (2) fragmentation and immersion of solid or semi-solid surface crusts of pahoehoe by roiling movements of the flow, forming cores of discrete, tacky masses; (3) sudden renewed movement of lava stored and cooled within surface reservoirs to form clots. The masses, fragments, and clots in these transition modes are characterized by spinose, granulated surfaces; as flow movement continues, the masses and fragments aggregate, fracture, and grind together, completing the transition to aa. Observations show that the critical relation between viscosity and rate of shear strain is inverse: if viscosity is low, a high rate of shear is required to begin the transition to aa; conversely, if viscosity is high, a much lower rate of shear will induce the transition. These relations can be demonstrated qualitatively with simple graphs, which can be used to examine the flow history of any selected finite lava element by tracing the path represented by its changing viscosity and shear rate. A broad, diffuse “transition threshold zone” in these graphs portrays the inverse critical relation between viscosity and shear rate; the transition to aa is represented by the path of the lava element crossing this zone. Moving lava flows can be regarded as natural viscometers, by which shear stress and rate of shear strain at selected points can be determined and viscosity can be computed. By making such determinations under a wide range of conditions on pahoehoe, aa, and transitional flow types, the critical relations that control the pahoehoe-aa transition can be quantified.

Hawaii↗
USGS2012

The mechanics and three-dimensional internal structure of active magmatic systems: Kilauea volcano, Hawaii

Interpretation of abundant seismic data suggests that Kilauea's primary conduit within the upper mantle is concentrically zoned to about 34-km depth. This zoned structure is inferred to contain a central core region of relatively higher permeability, surrounded by numerous dikes that are in intermittent hydraulic communication with each other and with the central core. During periods of relatively high magma transport, the entire cross section of the conduit is utilized. During periods of relatively low to moderate transport, however, only the central core is active. As the conduit penetrates the oceanic crust and enters the volcanic shield, it simultaneously supplies the deeper sections of the rift zones (6-to 10-km depth) and the roots of the summit reservoir with picritic magma. The rift zones at depth are inferred to be almost wholly molten and to possess a high degree of fluid continuity from Heiheiahulu in the East Rift Zone, 45 km westward through the roots of the summit magma reservoir, and well into the Southwest Rift Zone. Higher in the shield, the subcaldera magma reservoir and the shallow rift zones occupy the 2-to 4-km depth interval. Summit-differentiated olivine tholeiite (ρ ≈ 2.62 g cm −3 ) is periodically injected laterally along a horizon of neutral buoyancy within the rift zones, where the density of the magma is just balanced by the in situ density of the shield (Ryan, 1987 a , b ). Deep rift zone intrusions push seaward the deep tectonic blocks of the volcano's south flank. Shallow rift intrusions build a sheeted dike complex, inferred to be in isostatic equilibrium with the higher-density deep rift cores below. General finite element analyses are presented for the deformation and stress fields surrounding such dikes in the horizontal and vertical planes. The dike tip in two and three dimensions is surrounded by a tubular core of tensile (σ 1 , σ 2 ) and shear stress (τ max ). The displacement field is characterized by counterrotating cells on either side of the dike tip which, in vertical orientation, produce the characteristic subsidence above the dike complex, with uplift on either side, forming a ridge-trough-ridge structure. A finite element model of Kilauea's shield computes the displacement fields and principal stress (σ 1 ) distributions resulting from intrusive activity on each or both of the rift zones. Within the summit region, tensile stress lobes produced by the three-dimensional upward extension of the intrusions superpose constructively to produce calderawide regimes of tensile stress, conducive to caldera development. Parametric studies of (1) intrusion in the East Rift Zone only, (2) intrusion in the Southwest Rift Zone only, and (3) intrusion in both rift zones demonstrate their unique kinematic contributions. For case 1, the caldera undergoes a counterclockwise rotation (torque up state) conducive to the development of rightstepping en echelon eruptive fissures, as exemplified by the August 14, 1971, eruption. For case 2, the caldera undergoes a clockwise rotation (torque down state) conducive to the development of left-stepping eruptive fissures, as occurred during the December 31, 1974, eruption. For case 3, the caldera substructure is driven due southward, producing the southward migration of the upper portions of the summit magma reservoir.

Journal of Geophysical Research Solid Earth↗
USGS2010

Geology, geochemistry, and genesis of the Greens Creek massive sulfide deposit, Admiralty Island, southeastern Alaska

In 1996, a memorandum of understanding was signed by representatives of the U.S. Geological Survey and Kennecott Greens Creek Mining Company to initiate a cooperative applied research project focused on the Greens Creek massive sulfide deposit in southeastern Alaska. The goals of the project were consistent with the mandate of the U.S. Geological Survey Mineral Resources Program to maintain a leading role in national mineral deposits research and with the need of Kennecott Greens Creek Mining Company to further development of the Greens Creek deposit and similar deposits in Alaska and elsewhere. The memorandum enumerated four main research priorities: (1) characterization of protoliths for the wall rocks, and elucidation of their alteration histories, (2) determination of the ore mineralogy and paragenesis, including metal residences and metal zonation within the deposit, (3) determination of the ages of events important to ore formation using both geochronology and paleontology, and (4) development of computer models that would allow the deposit and its host rocks to be examined in detail in three dimensions. The work was carried out by numerous scientists of diverse expertise over a period of several years. The written results, which are contained in this Professional Paper, are presented by 21 authors: 13 from the U.S. Geological Survey, 4 from Kennecott Greens Creek Mining Company, 2 from academia, and 2 from consultants. The Greens Creek deposit (global resource of 24.2 million tons at an average grade of 13.9 percent zinc, 5.1 percent lead, 0.15 troy ounce per ton gold, and 19.2 troy ounces per ton silver at zero cutoff) formed in latest Triassic time during a brief period of rifting of the Alexander terrane. The deposit exhibits a range of syngenetic, diagenetic, and epigenetic features that are typical of volcanogenic (VMS), sedimentary exhalative (SEDEX), and Mississippi Valley-type (MVT) genetic models. In the earliest stages of rifting, formation of precious-metal-rich silica-barite-carbonate white ores began at low temperature in a shallow, subaqueous setting, probably a thin carbonate shelf on the flanks of the Alexander landmass. Epigenetic carbonate replacement textures in the footwall dolostones are overlain by stratiform silica-carbonate-barite-rich ores and indicate that early mineralization formed at and just beneath the paleo sea floor by mixing of a reduced, precious-metal-rich, base-metal-poor hydrothermal fluid with oxygenated seawater. As rifting intensified, the shelf was downfaulted and isolated as a graben. Isolation of the basin and onset of starved-basin shale sedimentation was concurrent with emplacement of mafic-ultramafic intrusives at shallow levels in the rift, resulting in an increasingly higher temperature and progressively more anoxic ore-forming environment. The formation of the main stage of massive sulfide ores began as the supply of bacterially reduced sulfur increased in the accumulating shales. As the main-stage mineralization intensified, shale sedimentation inundated the hydrothermal system, eventually forming a cap. Biogenic sulfate reduction supplied reduced sulfur to the base of the shales where mixing occurred with hot, base-metal-rich hydrothermal fluids. Ore deposition continued by destruction and epigenetic replacement of the early white ores in proximal areas and by inflation and diagenetic replacement of unlithified shale at the interface between the white ores and the base of the shale cap. Ore deposition waned as the shales became lithified and as the supply of bacterially reduced sulfur to the site of ore deposition ceased. The final stages of rifting resulted in the emplacement of mafic-ultramafic intrusive rocks into the Greens Creek system and extrusion of voluminous basaltic flows at the top of the Triassic section. Greenschist facies metamorphism during the Jurassic-Cretaceous accretion of the Alexander terrane to the continental margin resulted in recrystalli

Professional Paper↗
USGS2018

Modeled inundation limits of potential lahars from Mount Adams in the White Salmon River Valley, Washington

Lahars large enough to reach populated areas are a hazard at Mount Adams, a massive volcano in the southern Cascade Range of Washington State (fig. 1). It is considered to be still active and has the potential to erupt again. By definition, lahars are gravity-driven flows of water-saturated mixtures of mud and rock (plus or minus ice, wood, and other debris), which originate from volcanoes and have a variety of potential triggering mechanisms (Vallance, 2000; Vallance and Iverson, 2015). Flowing mixtures can range in fluid consistency from something like a milkshake to something more like wet concrete, and they behave like flash floods, in that they can appear suddenly in river channels with little warning and commonly have boulder- or log-choked flow fronts. Lahars are hazardous because they can flow rapidly in confined valleys (commonly 20–35 miles per hour [mph] or 9–16 meters per second [m/s]), can travel more than 100 miles (mi) (161 kilometers [km]) from a source volcano, and can move with incredible destructive force, carrying multi-ton boulders and logs that can act as battering rams (Pierson, 1998). The biggest threats from lahars to downstream communities are present during eruptive activity, and impacts to communities can be dire. For example, a very large eruption-triggered lahar in Colombia in 1985 surprised and killed more than 20,000 people in a large town located about 45 mi (72 km) downstream and out of sight of the volcano that produced it (Pierson and others, 1990). Mount Adams, one of the largest volcanoes in the Cascade Range, is a composite stratocone composed primarily of andesite lava flows. It has been the most continuously active volcano within the 480-mi 2 Mount Adams volcanic field—a region covering parts of Klickitat, Skamania, Yakima, andLewis Counties and part of the Yakama Nation Reservation in Washington State (Hildreth and Fierstein,1995, 1997). About 500,000 years in age, Mount Adams reached its present size by about 15,000 years ago, primarily through the episodic effusion of lava flows; it has not had a history of major explosive eruptions like Mount St. Helens, its neighbor to the west. Timing of the most recent eruptive activity (recorded by four thin tephra layers) is on the order of 1,000 years ago; the tephras are bracketed by 2,500-year-old and 500-year-old ash layers from Mount St. Helens (Hildreth and Fierstein, 1995, 1997). Mount Adams currently shows no signs of renewed unrest. Eruptive history does not tell us everything we need to know about hazards at Mount Adams, however, which are fully addressed in the volcano hazard assessment for Mount Adams (W.E. Scott and others, 1995). This volcano has had a long-active hydrothermal system that circulated acidic hydrothermal fluids, formed by the solution of volcanic gases in heated groundwater, through fractures and permeable zones into upper parts of the volcanic cone. Acid sulfate leaching of rocks in the summit area may still be occurring, but chemical and thermal evidence suggests that the main hydrothermal system is no longer active at Mount Adams (Nathenson and Mariner, 2013). However, these rock-weakening chemical reactions have operated long enough to change about 0.4 cubic miles (mi 3 ) (1.7 cubic kilometers [km 3 ]) of the hard lava rock in the volcano’s upper cone to a much weaker clay-rich rock, thus significantly reducing rock strength and thereby slope stability in parts of the cone (Finn and others, 2007). The two largest previous lahars from Mount Adams were triggered by landslides of hydrothermally altered rock from the upper southwestern flank of the cone, and any future large lahars are likely to be triggered by the same mechanism. Mount Rainier also has had extensive hydrothermal alteration of rock in its upper edifice, and it also has a history of large landslides that transform into lahars (K.M. Scott and others, 1995; Vallance and Scott, 1997; Reid and others, 2001). The spatial depiction of modeled lahar inundation zones accompanying this report, shown in two different map perspectives, is intended to augment (not replace) the existing hazard maps for Mount Adams (W.E. Scott and others, 1995; Vallance, 1999). The maps in this report show potential areas of inundation by lahars of different initial volumes, which are determined by a computer model, LAHARZ (Iverson and others, 1998; Schilling, 1998). One map sheet presents LAHARZ-determined inundation areas on a normal plan-view shaded-relief map of the study area; the other gives an oblique perspective of the landscape with raised topography, as if one were viewing the landscape at an angle from an aircraft (Jenny and Patterson, 2007). LAHARZ was developed after the original hazard maps (based only on mapping of geologic deposits) were made. Predicted inundation zones on these maps provide an alternative approach to estimation of areas that could be inundated as lahars of different volumes pass through the valley. However, there is considerable uncertainty in the exact location of the hazard-zone boundaries shown on these maps, as well as on earlier maps.

Washington↗
USGS2026

Multiple-well monitoring site adjacent to the Midway- Sunset and Buena Vista Oil Fields, Kern County, California

Groundwater quality in and around oil fields in the Southern San Joaquin Valley is of interest to many California residents that rely heavily on groundwater for domestic, commercial, and agricultural use. To help assess the effects of historical oil-field activities and natural geologic sources on groundwater near the southwest margins of the Kern County Groundwater Subbasin, a multiple-well monitoring site was installed near the administrative boundary between the Midway-Sunset and Buena Vista Oil Fields in Kern County, California. The installation of the Midway-Sunset Buena Vista multiple-well monitoring site (MSBV) supports regional analysis of the relations of oil and gas sources to groundwater quality by providing information about the geology, hydrology, geophysical properties, and water quality of the alluvial and upper Tulare aquifers in areas where groundwater data were limited. Data collected from the site included drill cuttings, whole core samples, sidewall core samples, mud-gas analysis, borehole geophysical logs, depth to water measurements, and water quality samples. Whole cores were scanned using dual energy computed tomography. Subsamples of selected cores were analyzed for density, porosity, specific retention, and bulk minerology. Thin sections of the subsamples were prepared, photographed, and examined. Two samples were analyzed using scanning electron microscope technology to examine the microporosity of diatomite laden sediment. Instrumentation installed in the wells collect hourly depth to water measurements. Analysis of the data show there is 355 feet of alluvium overlying the Tulare Formation at the well site. The contact between the two formations is an aquitard resulting in a perched aquifer in the alluvium and unconfined aquifer in the Tulare Formation. The alluvium is more heterogenous and finer grained than the Tulare Formation resulting in markedly higher porosity in the alluvium compared to the Tulare Formation. Higher specific retention observed in the alluvium is attributed to the finer grained sediment and greater abundance of reworked diatomite (as represented by opal-CT [cristobalite-tridymite]) compared to the Tulare Formation. Total dissolved solids (TDS) approached or exceeded 10,000 milligrams per liter (mg/L) in the alluvium from approximately 176 to 242 feet below land surface and at the top of the Amnicola clay at approximately 670 feet below land surface within the Tulare Formation. Elevated TDS, chloride, and boron concentrations in the alluvium and on top of the Amnicola clay likely reflect groundwater that is mixed with oil-field water. Water chemistry and modern-aged groundwater in the alluvial monitoring well (MSBV #3) are consistent with the oil-field water in the alluvium being derived from documented historical surface disposal of oil-field water upslope (northwest) of the site. Water chemistry and pre-modern groundwater age in the deeper Tulare monitoring well (MSBV #1) on top of the Amnicola clay are consistent with oil-field fluids derived from upslope natural geologic sources or old oil wells that leak in the subsurface. Shallow groundwater in the Tulare (MSBV #2) is not affected by mixing with oil-field sources.

EarthArXiv↗
USGS2017

Using CO 2 Prophet to estimate recovery factors for carbon dioxide enhanced oil recovery

Introduction The Oil and Gas Journal’s enhanced oil recovery (EOR) survey for 2014 (Koottungal, 2014) showed that gas injection is the most frequently applied method of EOR in the United States and that carbon dioxide (CO 2 ) is the most commonly used injection fluid for miscible operations. The CO 2 -EOR process typically follows primary and secondary (waterflood) phases of oil reservoir development. The common objective of implementing a CO 2 -EOR program is to produce oil that remains after the economic limit of waterflood recovery is reached. Under conditions of miscibility or multicontact miscibility, the injected CO 2 partitions between the gas and liquid CO2 phases, swells the oil, and reduces the viscosity of the residual oil so that the lighter fractions of the oil vaporize and mix with the CO 2 gas phase (Teletzke and others, 2005). Miscibility occurs when the reservoir pressure is at least at the minimum miscibility pressure (MMP). The MMP depends, in turn, on oil composition, impurities of the CO 2 injection stream, and reservoir temperature. At pressures below the MMP, component partitioning, oil swelling, and viscosity reduction occur, but the efficiency is increasingly reduced as the pressure falls farther below the MMP. CO 2 -EOR processes are applied at the reservoir level, where a reservoir is defined as an underground formation containing an individual and separate pool of producible hydrocarbons that is confined by impermeable rock or water barriers and is characterized by a single natural pressure system. A field may consist of a single reservoir or multiple reservoirs that are not in communication but which may be associated with or related to a single structural or stratigraphic feature (U.S. Energy Information Administration [EIA], 2000). The purpose of modeling the CO 2 -EOR process is discussed along with the potential CO 2 -EOR predictive models. The data demands of models and the scope of the assessments require tradeoffs between reservoir-specific data that can be assembled and simplifying assumptions that allow assignment of default values for some reservoir parameters. These issues are discussed in the context of the CO 2 Prophet EOR model, and their resolution is demonstrated with the computation of recovery-factor estimates for CO 2 -EOR of 143 reservoirs in the Powder River Basin Province in southeastern Montana and northeastern Wyoming.

Scientific Investigations Report↗