USGS Science⌕ Search

SEARCH · USGS Science

Results for “Interfaces”

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 307 records · Page 17Linked to original sources
USGS2012

Borehole geophysical, fluid, and hydraulic properties within and surrounding the freshwater/saline-water transition zone, San Antonio segment of the Edwards aquifer, south-central Texas, 2010-11

The freshwater zone of the San Antonio segment of the Edwards aquifer is used by residents of San Antonio and numerous other rapidly growing communities in south-central Texas as their primary water supply source. This freshwater zone is bounded to the south and southeast by a saline-water zone with an intermediate zone transitioning from freshwater to saline water, the transition zone. As demands on this water supply increase, there is concern that the transition zone could potentially move, resulting in more saline water in current supply wells. Since 1985, the U.S. Geological Survey (USGS), San Antonio Water System (SAWS), and other Federal and State agencies have conducted studies to better understand the transition zone. During 2010 and 2011, the USGS, in cooperation with SAWS, conducted a study to further assess the potential for movement of the transition zone in part of the San Antonio segment of the Edwards aquifer. Equivalent freshwater heads were computed to investigate the transition from saline to freshwater zones in the San Antonio segment and evaluate the potential for lateral flow at the freshwater/saline-water interface. Data were collected within and surrounding the transition zone from 13 wells in four transects (East Uvalde, Tri-County, Fish Hatchery, and Kyle). Hydraulic head and geophysical log data were used to calculate equivalent freshwater heads and then analyzed to identify possible horizontal gradients across the transition zone and thus flow. Unlike previous studies that used indirect methods to calculate fluid conductivity from fluid resistivity, in this study geophysical tools that directly measured fluid conductivity were used. Electromagnetic (EM) flowmeter logs were collected under both ambient and stressed (pumping) conditions and were processed to identify vertical flow zones within the borehole. The San Antonio segment of the Edwards aquifer (the study area) is about 175 miles long and extends from the western groundwater divide near Brackettville in Kinney County to the eastern groundwater divide near Kyle in Hays County. The four transects consist of two to five wells per transect and were configured approximately perpendicular to and across the expected trace of the freshwater/saline-water interface. The deep flow zone indicated by the EM flowmeter data for East Uvalde transect well EU2 corresponds directly with a large, negative deflection of the fluid logs, indicating an inflow of fresher water from the Devils River Limestone. To the southwest, towards the freshwater/saline-water interface, this same flow zone was observed in well EU1, but with a reduction of flow, and displayed no apparent fluid curve deflections. The highest observed transmissivity of the study area was observed in the saline zone of the Tri-County transect, at well TC3, which had a total transmissivity of 24,900 square feet per day. Zones of high transmissivity throughout the study site were observed to not be continuous and are likely caused by localized secondary porosity such as intersecting faults or karst features. Although analyses of daily mean equivalent freshwater heads for the East Uvalde transect indicated that the gradient across the freshwater/saline-water interface varied between into and out of the freshwater zone, the data indicate that there was a slightly longer period during which the gradient was out of the freshwater zone. Analyses of all daily mean equivalent freshwater heads for the Tri-County transect indicated that the lateral-head gradients across the freshwater/saline-water interface were typically mixed (not indicative of flow into or out of freshwater zone). Assessment of the daily mean equivalent freshwater heads indicated that, although the lateral-head gradient at the Kyle transect varied between into and out of the freshwater zone, the lateral-head gradient was typically from the transition zone into the freshwater zone.

Texas↗
USGS2015

Updated numerical model with uncertainty assessment of 1950-56 drought conditions on brackish-water movement within the Edwards aquifer, San Antonio, Texas

In 2010, the U.S. Geological Survey, in cooperation with the San Antonio Water System, began a study to assess the brackish-water movement within the Edwards aquifer (more specifically the potential for brackish-water encroachment into wells near the interface between the freshwater and brackish-water transition zones, referred to in this report as the transition-zone interface) and effects on spring discharge at Comal and San Marcos Springs under drought conditions using a numerical model. The quantitative targets of this study are to predict the effects of higher-than-average groundwater withdrawals from wells and drought-of-record rainfall conditions of 1950–56 on (1) dissolved-solids concentration changes at production wells near the transition-zone interface, (2) total spring discharge at Comal and San Marcos Springs, and (3) the groundwater head (head) at Bexar County index well J-17. The predictions of interest, and the parameters implemented into the model, were evaluated to quantify their uncertainty so the results of the predictions could be presented in terms of a 95-percent credible interval. The model area covers the San Antonio and Barton Springs segments of the Edwards aquifer; the history-matching effort was focused on the San Antonio segment. A previously developed diffuse-flow model of the Edwards aquifer, which forms the basis for the model in this assessment, is primarily based on a conceptualization in which flow in the aquifer is predominately through a network of numerous small fractures and openings. Primary updates to this model include an extension of the active area downdip, a conversion to an 8-layer SEAWAT variable-density flow and transport model to simulate dissolved-solids concentration effects on water density, history matching to 1999–2009 conditions, and parameter estimation in a highly parameterized context using automated methods in PEST (a model-independent Parameter ESTimation code). In addition to the best-fit parameter values derived from history matching, the uncertainty of model parameters was also estimated by using linear uncertainty analysis. Comparison of “prior” (before history matching) and “posterior” (after history matching) variances of parameters indicate that the information within the observation dataset used for history matching informs many parameters. The concentration threshold parameters were well-informed by the observation dataset as their posterior distributions were much narrower than their prior distributions. The transition-zone scaling parameters of hydraulic conductivity, effective porosity, and specific storage were all informed by the observation dataset, as evidenced by the difference between the prior and posterior variances. Saline-zone scaling parameters, alternatively, were not informed by the observation dataset for effective porosity and specific storage. Resulting posterior drier-month, wetter-month, and annual recharge multiplier parameter variances are important to understanding how well recharge is estimated and implemented within the model. The shifts of the posterior distributions left and right indicate that there were zones where less or more water was needed in the model. The widths of the distributions were not decreased substantially, indicating that many of the best-fit recharge parameters are not statistically different from the initial values specified in the history-matching effort. Recharge from rainfall is the driving force behind groundwater flow and heads in the aquifer; therefore, an increase in understanding of this process would benefit model development by potentially decreasing the uncertainty of this parameter. The history-matching effort was most helpful in informing the parameters in the model that control discharge at springs, namely, the spring orifice (drain) altitude and drain conductance parameters for each spring. The uncertainty assessment of the predictive model (a hypothetical recurrence of 1950–56 drought conditions and higher-than-average groundwater withdrawals from wells) provided insights into the potential effects of these conditions on dissolved-solids concentration changes at production wells near the transition-zone interface, discharges at Comal and San Marcos Springs, and heads at Bexar County index well J-17. Results at the 25 production wells near the transition-zone interface indicate that the uncertainty of model input parameters based on expert knowledge yielded an upper bound of the 95-percent credible interval of dissolved-solids concentrations that exceeds the secondary drinking water standards of 1,000 milligrams per liter (mg/L) of the Texas Commission on Environmental Quality (TCEQ) for many wells. However, the history-matching process provided key information to inform prediction-sensitive model parameters and therefore, contributed to a substantial decrease of the upper bound of the 95-percent credible interval to below the secondary drinking water standards. Reductions in dissolved-solids concentration changes were on the order of 400 mg/L to 1,300 mg/L. The reduction in uncertainty in regards to this prediction implies that this prediction of dissolved-solids concentration change can be made with some certainty using this current model and that those parameters that control this prediction are informed by the observation dataset. Even though predictive uncertainty was reduced for this prediction, dissolved-solids concentration changes were still greater than zero, indicating a minimal increase in concentration at these 25 production wells during the 7-year simulation period is likely. However, this minimal concentration increase indicates a small potential for movement of the brackish-water transition zone near these wells during the 7-year simulation period of drought-ofrecord (1950–56) rainfall conditions with higher-than-average groundwater withdrawals by wells. Predictive results of total spring discharge during the 7-year period, as well as head predictions at Bexar County index well J-17, were much different than the dissolved-solids concentration change results at the production wells. These upper bounds are an order of magnitude larger than the actual prediction which implies that (1) the predictions of total spring discharge at Comal and San Marcos Springs and head at Bexar County index well J-17 made with this model are not reliable, and (2) parameters that control these predictions are not informed well by the observation dataset during historymatching, even though the history-matching process yielded parameters to reproduce spring discharges and heads at these locations during the history-matching period. Furthermore, because spring discharges at these two springs and heads at Bexar County index well J-17 represent more of a cumulative effect of upstream conditions over a larger distance (and longer time), many more parameters (with their own uncertainties) are potentially controlling these predictions than the prediction of dissolved-solids concentration change at the prediction wells, and therefore contributing to a large posterior uncertainty.

Texas↗
USGS2019

Use of a Numerical Model to Simulate the Hydrologic System and Transport of Contaminants Near Joint Base Cape Cod, Western Cape Cod, Massachusetts

Historical training and operational activities at Joint Base Cape Cod (JBCC) on western Cape Cod, Massachusetts, have resulted in the release of contaminants into an underlying glacial aquifer that is the sole source of water to the surrounding communities. Remedial systems have been installed to contain and remove contamination from the aquifer. Groundwater withdrawals for public supply are expected to increase as the region continues to urbanize. Increases in water-supply withdrawals and wastewater return flow likely will affect the hydrologic system around JBCC and could affect the transport of any contamination that may remain in the aquifer following remediation of contamination from the JBCC. The U.S. Geological Survey, in cooperation with the Air Force Civil Engineer Center, developed a numerical, steady-state regional model of the Sagamore flow lens on western Cape Cod and evaluated the potential effects of future (2030) groundwater withdrawals on water levels, streamflows, hydraulic gradients, and advective transport near the JBCC. The aquifer consists generally of sandy sediments underlain by impermeable bedrock and is bounded laterally by a freshwater/saltwater interface. Data on the altitude of the bedrock surface, position of the freshwater/saltwater interface, lithology of the aquifer, spatial distribution of recharge, and hydrologic boundaries were incorporated into the three-dimensional, finite-difference groundwater flow model. Some inputs into the numerical model—aquifer properties, leakances, and recharge—are represented as parameters to facilitate estimation of optimal parameter values in an inverse calibration. A hybrid parameterization scheme, with both zones of piecewise constancy and pilot points, is used to represent hydraulic conductivity; other adjustable parameters include recharge, boundary leakance, and porosity. Data on water levels, the distribution of subsurface contamination, and groundwater ages were compiled, evaluated, and used to develop observations of long-term average hydraulic gradients and advective-transport patterns. These observations of steady-state hydrologic conditions were combined with the parameterized groundwater model in an inverse calibration to estimate model parameters that best fit the observations. Current (2010) and future (2030) conditions were simulated in the calibrated model to characterize the groundwater flow system and to determine potential effects of increased groundwater withdrawals on advective-transport patterns at the JBCC. Groundwater flow and advective transport are radially outward from a water-table divide in the northern part of the JBCC; flow diverges from the divide toward all points of the compass. Most groundwater flow and contaminant transport occur in shallow parts of the aquifer. On average, about one-half of the groundwater flux occurs in the shallowest 20 percent of the saturated thickness; shallow flow is even more predominant near streams and lakes. Projected (2030) increases in groundwater withdrawals decrease water levels by a maximum of about 1.2 feet in the northern part of the JBCC; drawdowns exceeding 1 foot generally are limited to areas near the largest increases in withdrawals, such as in the northern part of the JBCC, near Long Pond in Falmouth, and in eastern Barnstable. Streamflow decreases average about 6 percent; the largest decreases are in areas with the largest drawdowns. Changes in hydraulic-gradient directions at the water table exceed 1 degree in about 13 percent of the aquifer, generally near groundwater divides where gradient magnitudes are small and near large groundwater withdrawals. Predictions of advective transport from randomly selected locations at the water table are similar for current (2010) and future (2030) groundwater withdrawals. The results indicate that projected increases in groundwater withdrawals affect water levels and streamflows, but effects on hydraulic gradients and advective transport at the JBCC likely are small. Several underlying assumptions inherent in the model, including observations and weights used in the calibration, representation of local-scale heterogeneity, and simulation of the freshwater/saltwater interface, could affect model calibration and predictions; these assumptions were evaluated with alternative models and alternative inverse calibrations. Eight alternative calibrations were performed in which different, but reasonable, observations and weights were used. The preferred calibrated model had the best overall fit to the observations. Fine-grained silty sediments occur in many parts of the aquifer, and silt lenses can locally affect hydraulic gradients. A set of alternative models in which silts were represented with different correlation distances and hydraulic conductivities indicated that explicitly representing silt lenses could affect model calibration but that the implicit representation of local-scale heterogeneity may be sufficient at the regional scale to represent regional-scale hydraulic gradients. For the coastal boundary, two alternative models representing silty and sandy seabeds and their associated interface positions were developed to test the importance of the assumed coastal-boundary condition. The two alternative models resulted in different predictions of streamflow—streamflows increase with smaller (silty) seabed leakances. However, predictions of advective transport, particularly near the JBCC, generally were similar between the alternative and preferred calibrated models, indicating that the seabed leakance and associated interface position at the coastal boundary does not affect simulations of advective transport in inland parts of the aquifer.

Massachusetts↗
USGS2001

Simulation of ground-water discharge to Biscayne Bay, southeastern Florida

As part of the Place-Based Studies Program, the U.S. Geological Survey initiated a project in 1996, in cooperation with the U.S. Army Corps of Engineers, to quantify the rates and patterns of submarine ground-water discharge to Biscayne Bay. Project objectives were achieved through field investigations at three sites (Coconut Grove, Deering Estate, and Mowry Canal) along the coastline of Biscayne Bay and through the development and calibration of variable-density, ground-water flow models. Two-dimensional, vertical cross-sectional models were developed for steady-state conditions for the Coconut Grove and Deering Estate transects to quantify local-scale ground-water discharge patterns to Biscayne Bay. A larger regional-scale model was developed in three dimensions to simulate submarine ground-water discharge to the entire bay. The SEAWAT code, which is a combined version of MODFLOW and MT3D, was used to simulate the complex variable-density flow patterns. Field data suggest that ground-water discharge to Biscayne Bay relative to the shoreline is restricted to within 300 meters at Coconut Grove, 600 to 1,000 meters at Deering Estate, and 100 meters at Mowry Canal. The vertical cross-sectional models, which were calibrated to the field data using the assumption of steady state, tend to focus ground-water discharge to within 50 to 200 meters of the shoreline. With homogeneous distributions for aquifer parameters and a constant-concentration boundary for Biscayne Bay, the numerical models could not reproduce the lower ground-water salinities observed beneath the bay, which suggests that further research may be necessary to improve the accuracy of the numerical simulations. Results from the cross-sectional models, which were able to simulate the approximate position of the saltwater interface, suggest that longitudinal dispersivity ranges between 1 and 10 meters, and transverse dispersivity ranges from 0.1 to 1 meter for the Biscayne aquifer. The three-dimensional, regional-scale model was calibrated to ground-water heads, canal baseflow, and the general position of the saltwater interface for nearly a 10-year period from 1989 to 1998. The mean absolute error between observed and simulated head values is 0.15 meter. The mean absolute error between observed and simulated baseflow is 3 x 105 cubic meters per day. The position of the simulated saltwater interface generally matches the position observed in the field, except for areas north of the Miami Canal where the simulated saltwater interface is located about 5 kilometers inland of the observed saltwater interface. Results from the regional-scale model suggest that the average rate of fresh ground-water discharge to Biscayne Bay for the 10-year period (1989-98) is about 2 x 105 cubic meters per day for 100 kilometers of coastline. This simulated discharge rate is about 6 percent of the measured surface-water discharge to Biscayne Bay for the same period. The model also suggests that nearly 100 percent of the fresh ground-water discharge is to the northern half of Biscayne Bay, north of the Cutler Drain Canal. South of the Cutler Drain Canal, coastal lowlands prevent the water table from rising high enough to drive measurable quantities of ground water to Biscayne Bay. Annual variations in sea-level elevation, which can be as large as 0.3 meter, have a substantial effect on rates of ground-water discharge. During 1989-98, simulated rates of ground-water discharge to Biscayne Bay generally are highest when sea level is relatively low.

Florida↗
USGS1988

Freshwater supply potential of the Atlantic Intracoastal Waterway near Myrtle Beach, South Carolina

A study was conducted to determine the low-flow frequency of freshwater flow in the Atlantic Intracoastal Waterway (AICW) near Myrtle Beach, South Carolina and to determine the effects of proposed freshwater withdrawals of 45 cu ft/sec at the location of the saltwater-freshwater interface. Discharges simulated in the AICW for 1982-86 using BRANCH one-dimensional flow model were used to establish a relation of 7-day average flows in the AICW to summed 7-day average flows of four tributary streams. This relation was used with the tributary records for 1954-86 climatic years to generate 7-day minimum flows of the AICW, which were then used to develop a low-flow frequency relation. The relation indicated that the 7-day, 10-year flow of the Atlantic Intracoastal Waterway is 192 ct ft/s. A relation of the mile position of the saltwater-freshwater interface to recorded specific conductances at Vereen 's Marina was established. The 1982-85 period of record of specific conductance was used to simulate interface positions which were then used to establish a relation of 7-day average interface position to 7-day average discharge of the AICW. This relation indicated that the 7-day average interface position would be at mile 355.5 for the 7Q10 and at mile 356.2 if 45 cu ft were withdrawn during the 7Q10. The analysis indicates that the AICW can provide a reliable supply of freshwater at the proposed withdrawal location at mi 363.3 in the vicinity of Myrtle Beach, even during the 7Q10 low-flow conditions. (USGS)

South Carolina↗
USGS2000

Flow and salinity characteristics of the upper Suwannee River Estuary, Florida

Continuous stage and salinity data were recorded from August 1995 to December 1997 at four gages located in the upper Suwannee River Estuary. Continuous velocity data were recorded at two of the four gages and continuous discharge data were computed for these two gages. Additional salinity data were collected at 15 monitoring sites from November 1992 to October 1997. Wind-speed data collected at Cedar Key, Florida, during the study period were utilized in the regression analysis. Correlations were developed to describe the longitudinal extent of the saltwater/freshwater interface (defined as 0.5 parts per thousand (ppt) salinity) and salinity distribution in the upper Suwannee River Estuary. On East Pass, the median of difference between daily maximum and daily minimum stage ranged from 2.92 feet for a gage at river mile 3.8 to 3.33 feet for a gage at river mile 1.2. Velocities tended to be unidirectional with some instances of bilateral flow. Reversal in flow direction was common and coincided with rising tides. Monthly mean discharges for the Suwannee River near Wilcox, Florida, during the study period typically were lower than the average for the period of record (1931.97). Discharge near Wilcox averaged 4,000 cubic feet per second (ft3/s) lower than the long-term average from June to September 1996. An El Ni�o event induced precipitation that was responsible for higher than average monthly mean discharge measured near Wilcox during November and December 1997. The maximum observed salinity concentrations for the study period ranged from 28.20 ppt at river mile 3.8 to 31.00 ppt at river mile 1.9. Median daily fluctuations of salinity at river miles 3.8 and 1.2 were 0.12 and 11.31 ppt, respectively. The maximum daily upstream extent of the saltwater/freshwater interface was at or upstream from river mile 4.0 for about 50 percent of the study period. The interface was at or upstream from river mile 3.8 and river mile 2.8 40 and 57 percent of the time. The interface was downstream from river mile 1.2 and river mile 1.9 11 and 21 percent of the time, respectively. The median daily maximum salinity for the four gages ranged from 0.22 ppt at river mile 3.8 to 11.50 ppt at river mile 1.2. Multiple linear-regression models were developed to determine the isohaline location for 0.5, 2, 5, 10, 15, and 20 ppt salinity, and to predict the maximum daily salinity concentrations at gages as a function of stage, river discharge, and wind. The salinity at a location was inversely proportional to the daily mean discharge at the Suwannee River near Wilcox. Under extreme low-flow conditions (3,500 ft3/s), the regression models predicted that the interface would occur at river mile 7.2, upstream from the Gopher River confluence with the Suwannee River. Wind speed did not have a substantial influence on model predictions. The period of record for the Suwannee River at Wilcox was applied to appropriate regression models to produce a synthetic record of historical salinity distributions. Two withdrawal scenarios, a 10-percent diversion and a 1,000 ft3/s diversion, were evaluated relative to high-, medium-, and low-flow conditions and compared to actual salinity distributions. The 10-percent and 1,000 ft3/s withdrawals scenarios resulted in the isohaline of 0.5 ppt migrating 0.6 and 1.58 miles upstream from the actual isohaline location for a low-flow condition of 4,500 ft3/s, and migrating 0.14 and 0.65 miles upstream from the actual isohaline location for a high-flow conditions of 20,300 ft3/s for Wadley Pass.

Florida↗
USGS1991

New hydrologic instrumentation in the U.S. Geological Survey

New water-level sensing and recording instrumentation is being used by the U.S. Geological Survey for monitoring water levels, stream velocities, and water-quality characteristics. Several of these instruments are briefly described. The Basic Data Recorder (BDR) is an electronic data logger, that interfaces to sensor systems through a serial-digital interface standard (SDI-12), which was proposed by the data-logger industry; the Incremental Shaft Encoder is an intelligent water-level sensor, which interfaces to the BDR through the SDI-12; the Pressure Sensor is an intelligent, nonsubmersible pressure sensor, which interfaces to the BDR through the SDI-12 and monitors water levels from 0 to 50 feet; the Ultrasonic Velocity Meter is an intelligent, water-velocity sensor, which interfaces to the BDR through the SDI-12 and measures the velocity across a stream up to 500 feet in width; the Collapsible Hand Sampler can be collapsed for insertion through holes in the ice and opened under the ice to collect a water sample; the Lighweight Ice Auger, weighing only 32 pounds, can auger 6- and 8-inch holes through approximately 3.5 feet of ice; and the Ice Chisel has a specially hardened steel blade and 6-foot long, hickory D-handle.

Conference Paper↗
USGS2002

Simulations of seismic hazard for the Pacific Northwest of the United States from earthquakes associated with the Cascadia subduction zone

We investigate the impact of different rupture and attenuation models for the Cascadia subduction zone by simulating seismic hazard models for the Pacific Northwest of the U.S. at 2% probability of exceedance in 50 years. We calculate the sensitivity of hazard (probabilistic ground motions) to the source parameters and the attenuation relations for both intraslab and interface earthquakes and present these in the framework of the standard USGS hazard model that includes crustal earthquakes. Our results indicate that allowing the deep intraslab earthquakes to occur anywhere along the subduction zone increases the peak ground acceleration hazard near Portland, Oregon by about 20%. Alternative attenuation relations for deep earthquakes can result in ground motions that differ by a factor of two. The hazard uncertainty for the plate interface and intraslab earthquakes is analyzed through a Monte-Carlo logic tree approach and indicates a seismic hazard exceeding 1 g (0.2 s spectral acceleration) consistent with the U.S. National Seismic Hazard Maps in western Washington, Oregon, and California and an overall coefficient of variation that ranges from 0.1 to 0.4. Sensitivity studies indicate that the paleoseismic chronology and the magnitude of great plate interface earthquakes contribute significantly to the hazard uncertainty estimates for this region. Paleoseismic data indicate that the mean earthquake recurrence interval for great earthquakes is about 500 years and that it has been 300 years since the last great earthquake. We calculate the probability of such a great earthquake along the Cascadia plate interface to be about 14% when considering a time-dependent model and about 10% when considering a time-independent Poisson model during the next 50-year interval.

Pure and Applied Geophysics↗
USGS1987

Analysis of saltwater upconing beneath a pumping well

Aquifer systems that contain freshwater and saltwater are usually stratified, with the more dense saltwater underlying the freshwater. A groundwater well discharging from the freshwater zone causes the saltwater to move upwards towards the well. This phenomenon is known as saltwater upconing. Two methods of analysis, the sharp-interface method and the fluid-density-dependent solute-transport method, are used to simulate saltwater upconing. Numerical experiments including comparisons of the two methods indicate: (1) for low to moderate pumpages the 50% isochlor and sharp interface correlate well; (2) the well can discharge significant concentrations of saltwater, even though a stable cone (according to the sharp-interface method) exists below the well screen; (3) an almost linear relationship exists between the well discharge rate and the concentration of the discharge at low pumping rates that maintain a stable cone; and (4) upconing is sensitive to transverse dispersivity, whereas it is insensitive to longitudinal dispersivity. A simulation of upconing at Test Site No. 4, Truro, Cape Cod, Massachusetts, indicates that the appropriate field value of transverse dispersivity is very small. This supports the validity of the sharp-interface assumption for analyzing the behavior of systems with thin saltwater-freshwater transition zones.

Journal of Hydrology↗
USGS2011

Modules based on the geochemical model PHREEQC for use in scripting and programming languages

The geochemical model PHREEQC is capable of simulating a wide range of equilibrium reactions between water and minerals, ion exchangers, surface complexes, solid solutions, and gases. It also has a general kinetic formulation that allows modeling of nonequilibrium mineral dissolution and precipitation, microbial reactions, decomposition of organic compounds, and other kinetic reactions. To facilitate use of these reaction capabilities in scripting languages and other models, PHREEQC has been implemented in modules that easily interface with other software. A Microsoft COM (component object model) has been implemented, which allows PHREEQC to be used by any software that can interface with a COM server—for example, Excel®, Visual Basic®, Python, or MATLAB". PHREEQC has been converted to a C++ class, which can be included in programs written in C++. The class also has been compiled in libraries for Linux and Windows that allow PHREEQC to be called from C++, C, and Fortran. A limited set of methods implements the full reaction capabilities of PHREEQC for each module. Input methods use strings or files to define reaction calculations in exactly the same formats used by PHREEQC. Output methods provide a table of user-selected model results, such as concentrations, activities, saturation indices, and densities. The PHREEQC module can add geochemical reaction capabilities to surface-water, groundwater, and watershed transport models. It is possible to store and manipulate solution compositions and reaction information for many cells within the module. In addition, the object-oriented nature of the PHREEQC modules simplifies implementation of parallel processing for reactive-transport models. The PHREEQC COM module may be used in scripting languages to fit parameters; to plot PHREEQC results for field, laboratory, or theoretical investigations; or to develop new models that include simple or complex geochemical calculations.

Computers & Geosciences↗
USGS2010

Controls on rind thickness on basaltic andesite clasts weathering in Guadeloupe

A clast of low porosity basaltic andesite collected from the B horizon of a soil developed on a late Quaternary volcaniclastic debris flow in the Bras David watershed on Basse-Terre Island, Guadeloupe, exhibits weathering like that observed in many weathered clasts of similar composition in other tropical locations. Specifically, elemental profiles measured across the core-rind interface document that primary minerals and glass weather to Fe oxyhydroxides, gibbsite and minor kaolinite in the rind. The earliest reaction identified in the core is oxidation of Fe in pyroxene but the earliest reaction that creates significant porosity is plagioclase dissolution. Elemental loss varies in the order Ca???Na>K???Mg>Si>Al>Fe???P??Ti, consistent with the relative reactivity of phases in the clast from plagioclase???pyroxene???glass>apatite>ilmenite. The rind surrounds a core of unaltered material that is more spherical than the original clast. The distance from the core-rind boundary to a visually prominent rind layer, L, was measured as a proxy for the rind thickness at 36 locations on a slab cut vertically through the nominal center of the clast. This distance averaged 24.4??3.1mm. Maximum and minimum values for L, 35.8 and 20.6mm, were observed where curvature of the core-rind boundary is greatest (0.12mm-1) and smallest (0.018mm-1) respectively. Extrapolating from other rinds in other locations, the rate of rind formation is estimated to vary by a factor of about 2 (from ~4 to 7??10-14ms-1) from low to high curvature. The observation of a higher rate of rind formation for a higher curvature interface is consistent with a diffusion-limited model for weathering rind formation. The diffusion-limited model predicts that, like rind thickness, values of the thickness of the reaction front (h) for a given reaction, defined as the zone over which a parent mineral such as plagioclase completely weathers to rind material, should also increase with curvature. Values of h were quantified as a function of interface curvature using bulk chemical analysis (500<h<2000??m). Values of h were also quantified by measuring loss of matrix glass and increase in porosity as a function of curvature. In contrast to rind thickness, h shows no consistent increase with curvature. This contradiction is attributed to the mm-scale roughness of the interface which is related to phenocryst grain size. Therefore, the overall rind formation rate is strongly affected by curvature measured at the scale of the clast, while mineral reaction rates documented by reaction front thickness are strongly affected by curvature at the scale of phenocrysts. Similarly, the weathering advance rate (ms-1) for the entire Bras David watershed can be extrapolated from the clast weathering rate if roughness at the watershed scale equals values of approximately 400-800. ?? 2010 Elsevier B.V.

Chemical Geology↗
USGS2023

The presence of silicate melt may enhance rates of cation diffusion in olivine

Olivine is commonly used as a ‘crystal clock’ to extract timescales relevant to pre-eruptive perturbations within mafic magmatic systems. Diffusion chronometry applications require accurate calibrations for the rates at which Fe-Mg or other commonly measured elements like Ni, Mn, and Ca diffuse through the crystal lattice. In the past, these rates have been mainly characterized using solid-solid diffusion couple experiments involving olivine single crystals, thin films, or powder sources. Despite the presence of melt surrounding olivine in natural magmatic systems, very few experiments involving magma have been performed, largely because controlling interface reactions is difficult. For this study, we carried out olivine-melt diffusion experiments as a test of the diffusion chronometry method, and to determine whether the presence of melt influences the calculated timescales. To approximate a natural system, we incorporated small natural Kīlauea and San Carlos olivine seeds within a natural Kīlauea basalt and tracked diffusive re-equilibration through time. To better control interface reactions, after some equilibration period at an initial superliquidus temperature of 1290 °C, the runs were rapidly cooled to form a rim and left to dwell at various final temperatures (1200, 1220, 1240, 1255 °C) for 6–84 h. Concentration gradients for Fe-Mg, Mn, Ni, Ca were measured, and the step-wise nature of the core-rim transition was ascertained using slow diffusing elements like P or Al. When these gradients are modeled using published diffusivities, the timescales retrieved are typically 10 times longer than the actual experiment durations. Thus, measured diffusivities are an order of magnitude faster than those previously obtained in olivine-solid source experiments, but they are in excellent agreement with the only two other melt-olivine datasets. We explore reasons for why melt-bearing olivine diffusion experiments tend to yield faster rates. The possible effects of (1) growth during diffusion, (2) diffusion during any initial dissolution step, and (3) extended tube or planar defects at the interface on calculated diffusivities are all considered but found to be inconsequential. Instead, we argue that additional point defects (vacancies) are likely created at the interface by higher concentrations in elements like Al or H in the basalt melt compared to other solid couple diffusant sources. Future applications of diffusion chronometry in olivine may require a complete re-evaluation of published diffusivities using melt-bearing experimental configurations.

Earth and Planetary Science Letters↗
USGS2020

Characterizing the diverse hydrogeology underlying rivers and estuaries using new floating transient electromagnetic methodology

The hydrogeology below large surface water features such as rivers and estuaries is universally under-informed at the long reach to basin scales (tens of km+). This challenge inhibits the accurate modeling of fresh/saline groundwater interfaces and groundwater/surface water exchange patterns at management-relevant spatial extents. Here we introduce a towed, floating transient electromagnetic (TEM) system (i.e. FloaTEM) for rapid (up to 15 km/h) high resolution electrical mapping of the subsurface below large water bodies to depths often a factor of 10 greater than other towed instruments. The novel FloaTEM system is demonstrated at a range of diverse 4th through 6th-order riverine settings across the United States including 1) the Farmington River, near Hartford, Connecticut; 2) the Upper Delaware River near Barryville, New York; 3) the Tallahatchie River near Shellmound, Mississippi; and, 4) the Eel River estuary, on Cape Cod, near Falmouth, Massachusetts. Airborne frequency-domain electromagnetic and land-based towed TEM data are also compared at the Tallahatchie River site, and streambed geologic scenarios are explored with forward modeling. A range of geologic structures and pore water salinity interfaces were identified. Process-based interpretation of the case study data indicated FloaTEM can resolve varied sediment-water interface materials, such as the accumulation of fines at the bottom of a reservoir and permeable sand/gravel riverbed sediments that focus groundwater discharge. Bedrock layers were mapped at several sites, and aquifer confining units were defined at comparable resolution to airborne methods. Terrestrial fresh groundwater discharge with flowpaths extending hundreds of meters from shore was also imaged below the Eel River estuary, improving on previous hydrogeological characterizations of that nutrient-rich coastal exchange zone. In summary, the novel FloaTEM system fills a critical gap in our ability to characterize the hydrogeology below surface water features and will support more accurate prediction of groundwater/surface water exchange dynamics and fresh-saline groundwater interfaces.

Science of the Total Environment↗
USGS2002

Hydrogeologic processes in saline systems: Playas, sabkhas, and saline lakes

Pans, playas, sabkhas, salinas, saline lakes, and salt flats are hydrologically similar, varying only in their boundary conditions. Thus, in evaluating geochemical processes in these systems, a generic water and solute mass-balance approach can be utilized. A conceptual model of a coastal sabkha near the Arabian Gulf is used as an example to illustrate the various water and solute fluxes. Analysis of this model suggests that upward flux of ground water from underlying formations could be a major source of solutes in the sabkha, but contribute only a small volume of the water. Local rainfall is the main source of water in the modeled sabkha system with a surprisingly large recharge-to-rainfall ratio of more than 50%. The contribution of seawater to the solute budget depends on the ratio of the width of the supratidal zone to the total width and is generally confined to a narrow zone near the shoreline of a typical coastal sabkha. Because of a short residence time of water, steady-state flow is expected within a short time (<100 years), while steady state for solutes may take much longer (>50,000 years). The solute composition of the brine in a closed saline system depends largely on the original composition of the input water. The high total ion content in the brine limits the efficiency of water-rock interaction and absorption. Because most natural systems are hydrologically open, the chemistry of the brines and the associated evaporite deposits may be significantly different than that predicted for hydrologically closed systems. Seasonal changes in temperature of the unsaturated zone cause precipitation of minerals in saline systems undergoing evaporation. Thus, during the hot dry season months, minerals exhibit retrograde solubility so that gypsum, anhydrite and calcite precipitate. Evaporation near the surface is also a major process that causes mineral precipitation in the upper portion of the unsaturated zone (e.g. halite and carnallite), provided that the relative humidity of the atmosphere is less than the activity of water. The slope of the fresh/brine-water interface in saline lake systems is shallower than in fresh/seawater interface because of the greater density difference between the fresh/brine-water bodies. The interface between sabkha brines and seawater slopes seaward, unlike normal marine-fresh water systems that slope landward. Moreover, the brine/seawater interface does not achieve steady state because it is pushed toward the sea by the sabkha's brine. ?? 2002 Elsevier Science B.V. All rights reserved.

Earth-Science Reviews↗
USGS2018

Time-dependent pore filling

Capillarity traps fluids in porous media during immiscible fluid displacement. Most field situations involve relatively long time scales, such as hydrocarbon migration into reservoirs, resource recovery, nonaqueous phase liquid remediation, geological CO 2 storage, and sediment‐atmosphere interactions. Yet laboratory studies and numerical simulations of capillary phenomena rarely consider the impact of time on these processes. We use time‐lapse microphotography to record the evolution of saturation in air‐ or hydrocarbon‐filled capillary tubes submerged in water to investigate long‐term pore filling phenomena beyond imbibition. Microphotographic sequences capture a lively pore filling history where various concurrent physical phenomena coexist. Dissolution and diffusion play a central role. Observations indicate preferential transport of the wetting liquid along corners, vapor condensation, capillary flow induced by asymmetrical interfaces, and interface pinning that defines the diffusion length. Other processes include internal snap‐offs, fluid redistribution, and changes in wettability as fluids dissolve into each other. Overall, the rate of pore filling is diffusion‐controlled for a given interfacial configuration; diffusive transport takes place at a constant rate for pinned interfaces and is proportional to the square root of time for free interfaces where the diffusion length increases with time.

Water Resources Research↗
USGS1997

Wrinkle-like slip pulse on a fault between different materials

Pulses of slip velocity can propagate on a planar interface governed by a constant coefficient of friction, where the interface separates different elastic materials. Such pulses have been found in two-dimensional plane strain finite difference calculations of slip on a fault between elastic media with wave speeds differing by 20%. The self-sustaining propagation of the slip pulse arises from interaction between normal and tangential deformation that exists only with a material contrast. These calculations confirm the prediction of Weertman [1980] that a dislocation propagating steadily along a material interface has a tensile change of normal traction with the same pulse shape as slip velocity. The self-sustaining pulse is associated with a rapid transition from a head wave traveling along the interface with the S wave speed of the faster material, to an opposite polarity body wave traveling with the slower S speed. Slip occurs during the reversal of normal particle velocity. The pulse can propagate in a region with constant coefficient of friction and an initial stress state below the frictional criterion. Propagation occurs in only one direction, the direction of slip in the more compliant medium, with rupture velocity near the slower S wave speed. Displacement is larger in the softer medium, which is displaced away from the fault during the passage of the slip pulse. Motion is analogous to a propagating wrinkle in a carpet. The amplitude of slip remains approximately constant during propagation, but the pulse width decreases and the amplitudes of slip velocity and stress change increase. The tensile change of normal traction increases until absolute normal traction reaches zero. The pulse can be generated as a secondary effect of a drop of shear stress in an asperity. The pulse shape is unstable, and the initial slip pulse can change during propagation into a collection of sharper pulses. Such a pulse enables slip to occur with little loss of energy to friction, while at the same time increasing irregularity of stress and slip at the source.

Journal of Geophysical Research B: Solid Earth↗
USGS2026

Sharp structural variability of the Gorda slab imaged by a fiber array

Structural variability at the subduction interface imposes stress and strength heterogeneity that affects the behavior of megathrust earthquakes. Characterizing this variability can be challenging because the slab interface is often deep and seismically quiet. Seismic wavefields from intraslab earthquakes can contain high-frequency phase conversions from the megathrust fault that are delayed and distorted by interface topography and velocity heterogeneity. Distributed acoustic sensing (DAS) allows for long-term deployments of dense seismic arrays that can yield many observations of these converted phases at high spatial resolution over large areas. Here, we develop a technique to identify and enhance slab-converted phases in DAS data, and we pick the differential arrival times and relative polarities between the direct and converted phases at short wavelength along the fiber optic cable. We apply this framework to earthquake wavefields measured by a DAS array at the southern end of the Cascadia Subduction Zone, above the locked zone, which suggests kilometer-scale interface roughness and variations in crustal thickness and absolute depth of the Gorda slab. Converted phase wavefields also suggest sharp variability in the presence and position of fluids in the slab and forearc crust.

California↗
USGS2007

Spectral element modelling of fault-plane reflections arising from fluid pressure distributions

The presence of fault-plane reflections in seismic images, besides indicating the locations of faults, offers a possible source of information on the properties of these poorly understood zones. To better understand the physical mechanism giving rise to fault-plane reflections in compacting sedimentary basins, we numerically model the full elastic wavefield via the spectral element method (SEM) for several different fault models. Using well log data from the South Eugene Island field, offshore Louisiana, we derive empirical relationships between the elastic parameters (e.g. P-wave velocity and density) and the effective-stress along both normal compaction and unloading paths. These empirical relationships guide the numerical modelling and allow the investigation of how differences in fluid pressure modify the elastic wavefield. We choose to simulate the elastic wave equation via SEM since irregular model geometries can be accommodated and slip boundary conditions at an interface, such as a fault or fracture, are implemented naturally. The method we employ for including a slip interface retains the desirable qualities of SEM in that it is explicit in time and, therefore, does not require the inversion of a large matrix. We performa complete numerical study by forward modelling seismic shot gathers over a faulted earth model using SEM followed by seismic processing of the simulated data. With this procedure, we construct post-stack time-migrated images of the kind that are routinely interpreted in the seismic exploration industry. We dip filter the seismic images to highlight the fault-plane reflections prior to making amplitude maps along the fault plane. With these amplitude maps, we compare the reflectivity from the different fault models to diagnose which physical mechanism contributes most to observed fault reflectivity. To lend physical meaning to the properties of a locally weak fault zone characterized as a slip interface, we propose an equivalent-layer model under the assumption of weak scattering. This allows us to use the empirical relationships between density, velocity and effective stress from the South Eugene Island field to relate a slip interface to an amount of excess pore-pressure in a fault zone. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗