Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Solutions Journal”

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 487 records · Page 27Linked to original sources
USGS1973

Determination of the association and dissociation of humic acid fractions by small angle X-ray scattering

A procedure has been devised for the fractionation of humic acid samples from different environments. This procedure involves fractionation of the sample by adsorption chromatography on a Sephadex G-50 column followed by chromatography on either a G-25 or a G-100 column. The fractions of the solutions are then examined by small angle X-ray scattering. Three different types of behavior have been detected among the humic acid fractions: (1) Some fractions show very little change in aggregation at pH values above 3.5. (2) One fraction forms aggregates at pH values above and below pH 7, but at pH 7 it is completely dissociated. (3) In some fractions the degree of aggregation decreases with increasing pH. However, even at pH values as high as 11.5 some large particles are still present. These differences in association behavior are due to the interaction of different attractive and repulsive forces. In many aggregating systems only one type of attractive force is dominant; however in humic acid systems hydrogen bonding, π bonding between planar aromatic moieties, and other coulombic interactions apparently all play a role in the formation of molecular aggregates.

Journal of Research of the U.S. Geological Survey↗
USGS2009

Compression-cuticle relationship of seed ferns: Insights from liquid-solid states FTIR (Late Palaeozoic-Early Mesozoic, Canada-Spain-Argentina)

Cuticles have been macerated from suitably preserved compressed fossil foliage by Schulze's process for the past 150 years, whereas the physical-biochemical relationship between the "coalified layer" with preserved cuticle as a unit has hardly been investigated, although they provide complementary information. This relationship is conceptualized by an analogue model of the anatomy of an extant leaf: "vitrinite (mesophyll) + cuticle (biomacropolymer) = compression". Alkaline solutions from Schulze's process as a proxy for the vitrinite, are studied by means of liquid-solid states Fourier transform infrared spectroscopy (FTIR). In addition, cuticle-free coalified layers and fossilized cuticles of seed ferns mainly from Canada, Spain and Argentina of Late Pennsylvanian-Late Triassic age are included in the study sample. Infrared data of cuticle and alkaline solutions differ which is primarily contingent on the mesophyll +biomacropolymer characteristics. The compression records two pathways of organic matter transformation. One is the vitrinized component that reflects the diagenetic-post-diagenetic coalification history parallel with the evolution of the associated coal seam. The other is the cuticle that reflects the sum-total of evolutionary pathway of the biomacropolymer, its monomeric, or polymeric fragmentation, though factors promoting preservation include entombing clay minerals and lower pH conditions. Caution is advised when interpreting liquid-state-based FTIR data, as some IR signals may have resulted from the interaction of Schulze's process with the cuticular biochemistry. A biochemical-study course for taphonomy is suggested, as fossilized cuticles, cuticle-free coalified layers, and compressions are responses to shared physicogeochemical factors. ?? 2009 Elsevier B.V. All rights reserved.

International Journal of Coal Geology↗
USGS1998

Nonlinear refraction and reflection travel time tomography

We develop a rapid nonlinear travel time tomography method that simultaneously inverts refraction and reflection travel times on a regular velocity grid. For travel time and ray path calculations, we apply a wave front method employing graph theory. The first-arrival refraction travel times are calculated on the basis of cell velocities, and the later refraction and reflection travel times are computed using both cell velocities and given interfaces. We solve a regularized nonlinear inverse problem. A Laplacian operator is applied to regularize the model parameters (cell slownesses and reflector geometry) so that the inverse problem is valid for a continuum. The travel times are also regularized such that we invert travel time curves rather than travel time points. A conjugate gradient method is applied to minimize the nonlinear objective function. After obtaining a solution, we perform nonlinear Monte Carlo inversions for uncertainty analysis and compute the posterior model covariance. In numerical experiments, we demonstrate that combining the first arrival refraction travel times with later reflection travel times can better reconstruct the velocity field as well as the reflector geometry. This combination is particularly important for modeling crustal structures where large velocity variations occur in the upper crust. We apply this approach to model the crustal structure of the California Borderland using ocean bottom seismometer and land data collected during the Los Angeles Region Seismic Experiment along two marine survey lines. Details of our image include a high-velocity zone under the Catalina Ridge, but a smooth gradient zone between. Catalina Ridge and San Clemente Ridge. The Moho depth is about 22 km with lateral variations. Copyright 1998 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Imaging with cross-hole seismoelectric tomography

We propose a cross-hole imaging approach based on seismoelectric conversions (SC) associated with the transmission of seismic waves from seismic sources located in a borehole to receivers (electrodes) located in a second borehole. The seismoelectric (seismic-to-electric) problem is solved using Biot theory coupled with a generalized Ohm's law with an electrokinetic streaming current contribution. The components of the displacement of the solid phase, the fluid pressure, and the electrical potential are solved using a finite element approach with Perfect Match Layer (PML) boundary conditions for the seismic waves and boundary conditions mimicking an infinite material for the electrostatic problem. We develop an inversion algorithm using the electrical disturbances recorded in the second borehole to localize the position of the heterogeneities responsible for the SC. Because of the ill-posed nature of the inverse problem (inherent to all potential-field problems), regularization is used to constrain the solution at each time in the SC-time window comprised between the time of the seismic shot and the time of the first arrival of the seismic waves in the second borehole. All the inverted volumetric current source densities are aggregated together to produce an image of the position of the heterogeneities between the two boreholes. Two simple synthetic case studies are presented to test this concept. The first case study corresponds to a vertical discontinuity between two homogeneous sub-domains. The second case study corresponds to a poroelastic inclusion (partially saturated by oil) embedded into an homogenous poroelastic formation. In both cases, the position of the heterogeneity is recovered using only the electrical disturbances associated with the SC. That said, a joint inversion of the seismic and seismoelectric data could improve these results.

Geophysical Journal International↗
USGS1977

The age of groundwater in the Lincolnshire Limestone, England and its relevance to the flow mechanism

Groundwater samples from the Lincolnshire Limestone have been analysed for tritium, radiocarbon, and the stable-isotope ratios 13 C 12 C "> 13 C 12 C , 18 O 16 O "> 18 O 16 O and D/H. The age of the water increases in a downgradient direction below overlying confining deposits and reaches a maximum age greater than 25,000 years within 15 km of the outcrop. The δ 13 C ratio ultimately attains the unusually low negative values of −2‰ in a downgradient direction; this is approaching that of the aquifer matrix which is +2.35‰. The reason is believed to be exchange of carbon between the groundwater and the matrix by a continuous process of precipitation and further solution of calcium carbonate. The δ 18 O and δ D ratios imply that recharge of the aquifer during the late Pleistocene took place in the spring and autumn rather than the winter as at present. The data are interpreted by assuming that movement of water through the saturated zone is a combination of flow in fissures and “piston flow” through the micro-fissures and pores of the rock. The mechanism of water movement in the saturated zone is dominated by relatively rapid flow in fissures, but the fissure flow includes a contribution of much older water from “intergranular” storage which enters the fissures from the rock matrix by pressure differentials in the fissure distribution system and by diffusion. The distribution of water of different ages in the aquifer is closely related to recent groundwater abstraction patterns.

Journal of Hydrology↗
USGS2012

Development and application of methods used to source prehistoric Southwestern maize: a review

Archaeological cobs free of mineral contaminants should be used to source the soils in which they were grown. Mineral contaminants often contain much higher concentrations of metals than vegetal materials and can alter a cob’s apparent metal and heavy-isotope content. Cleaning a cob via immersion in an acid solution for more than a few minutes will result in the incongruent and sometimes complete leaching of metals, including strontium (Sr), from the cob. When using 87 Sr/ Sr to determine the location of potential agriculture fields, it is best to either integrate several depth-integrated soil samples or to integrate several vegetation samples from individual fields. Biologically labile Sr in semi-arid Southwestern soils largely originates from eolian source or sources and usually is not derived from underlying bedrock. Existing Sr-isotope data indicate that archaeological cobs from Aztec Ruins came from either the Mesa Verde-McElmo Dome or Totah areas, that Pueblo Bonito and Chetro Ketl cobs, from Chaco Canyon that predate A.D. 1130, probably came from the Rio Chaco corridor, and that cobs from Chaco Canyon, that postdate A.D. 1130, probably came from either the Totah or Zuni areas.

Journal of Archaeological Science↗
USGS2024

Unravelling spatial heterogeneity of inundation pattern domains for 2D analysis of fluvial landscapes and drainage networks

Fluvial landscape analysis is an essential part of geomorphology, hydrology, ecology, and cartography. It is traditionally focused on the transition between hillslopes and channel domain, in which the network drainage is represented by static flow lines. However, the natural fluctuations of the processes occurring in the watershed induce lateral and longitudinal expansions and contractions in the drainage patterns and variations of stream surface area. These dynamics can be better understood by introducing a two-dimensional (2D) view of catchment hydrography, in which river width and floodplain are included in the analysis. The novelty introduced in this work is the development of a hydrodynamic hierarchical framework (HHF) to analyse the transitions among geomorphic and hydrographic features of the fluvial landscape, distinguishing hillslope, unchanneled valleys, floodplains, and single/multithreads channels. HHF is based on the estimation of nested inundation pattern domains (IPDs) from digital elevation models and 2D hydrodynamic modeling. IPDs are defined by scaling laws that characterize log–log relations between watershed drainage density and unit discharge thresholds extracted from a 2D direct rainfall method (DRM) under steady state solutions. The physical significance of the IPDs is analysed within the context of both the physiographic features of the fluvial landscape and the rainfall rates employed as input for the modeling approach. Initially, the spatial heterogeneity of the IPDs is used to derive stream width metrics as a function of the rainfall rate. Then, a spatial index, representative of the IPDs' heterogeneity, is introduced as a measure of the susceptibility of the drainage network surface area to expansion and contraction. Finally, the consistency of the results is assessed in comparison to another hydrodynamic-based method for fluvial landscape analysis recently proposed in the literature. The proposed approach is analysed using challenging mountain and low-relief environments, characterized by multithread channels, meander cut-offs, oxbow lakes, and extreme landscapes that feature glacial outwash, permafrost, and peatlands.

Alaska↗
USGS2012

Chemistry of hydrothermal solutions from the southern Juan de Fuca Ridge

Fluids from three vent fields on the southern Juan de Fuca Ridge were sampled in September 1984 using the DSRV Alvin. The fluids are uniquely high in both chloride, which ranges up to twice the seawater value, and in metal content. Simple vapor-liquid phase separation could not have produced both the high chlorinity and gas concentrations observed in these fluids. The cause of the elevated chlorinity can not be uniquely identified but may be the result of either or a combination of two processes: (1) subsurface formation of a degassed brine and subsequent mixing of a small amount of this brine with a hydrothermal seawater which has not undergone a phase separation and (2) dissolution of a chloride-rich phase combined with a possible small loss of gas during sampling procedures. Although measured temperatures were all less than 300°C, quartz geothermometry suggests that the fluids have equilibrated at greater than 340°C. Quartz geobarometry is also in agreement with geophysical estimates of depth to the local magma chamber.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Detection of a locked zone at depth on the Parkfield, California, segment of the San Andreas fault

The Parkfield, California, segment of the San Andreas fault is transitional in character between the creeping segment of the fault to the northwest and the locked Carrizo Plain segment to the southeast. The rate of shallow fault slip decreases from 25–30 mm/yr northwest of the epicenter of the 1966 Parkfield earthquake to zero at the southeastern end of the 1966 rupture zone. Data from a network of trilateration lines spanning the San Andreas fault near Parkfield and extending to the Pacific coast near San Luis Obispo shed light on the rate of fault slip at depth since the 1966 earthquake. In this study, average rates of line length change and shallow fault slip were inverted to determine the slip rate at depth on the Parkfield fault segment. The fault is taken to be a vertical surface with unknown distribution of strike-slip displacement in an elastic half-space. A striking result of the inversions is that all solutions providing acceptable fits to the data exhibit a locked zone essentially coincident with the rupture surface of the 1966 Parkfield earthquake. The data require that the locked zone extend nearly as far north as the 1966 epicenter; however, the vertical extent of the locked zone is not well resolved. Over much of the Parkfield segment the fault is slipping faster at the earth's surface than it is at seismogenic depths. In order to fit the trilateration measurements it is necessary to include a component of contraction normal to the trend of the San Andreas. The inversion results suggest a spatially uniform normal strain of −0.06 μstrain/yr. The orientation of the contraction is compatible with geologic and seismic evidence of active folding and reverse faulting in the region. The magnitude of the contraction is consistent with convergence rates inferred from global plate motion models.

Journal of Geophysical Research Solid Earth↗
USGS1981

Geochemical evidence for modern sediment accumulation on the continental shelf off southern New England

An area of fine-grained sediment approximately 170 km x 74 km in size, located in water depths between 60 m and 150 m, south of Martha's Vineyard, Mass., is a site of modern sediment deposition. The 14 C ages systematically increase with sediment depth from about 1,300 years B.P. at the surface to 8,000-10,000 years B.P. at the depth of maximum core penetration. The old age for the surface sediments probably results from a combination of deposition of old carbon and faunal mixing. In the finest sediments, the sedimentation rates were approximately 130 cm/1,000 yrs when deposition began and have decreased to about 25 cm/1,000 yrs. The decreasing sedimentation rate reflects a diminishing source of fine sediments, which presumably came from the Georges Bank and Nantucket Shoals area. Inventories of excess 210 Pb in undisturbed cores average 70 dpm/cm 2 (disintegrations per minute per square centimeter), more than two times higher than the flux of 210 Pb from the atmosphere and from 226 Ra decay in the overlying water. This additional influx of 210 Pb either must be with new fine-grained sediment material or from solutions that are stripped of their 210 Pb by particulates in the bottom nepheloid layer. Stable Pb concentrations in surface sediments are about 28 ppm, as much as two times higher than concentrations at depth. The high accumulation rates, 210 Pb inventories, and trace-metal profiles imply that this area is a modern sink for fine-grained sediments and for pollutants associated with particulate matter in the water column. To our knowledge, this is the only site of present-day natural deposition on the Continental Shelf off the eastern United States, exclusive of the Gulf of Maine. Because the net currents on the outer half of this Continental Shelf flow from northeast to southwest, this fine-grained deposit may receive its sediments and possible contaminants from the Nantucket Shoals and Georges Bank regions.

Massachusetts↗
USGS2008

Temporal response of hydraulic head, temperature, and chloride concentrations to sea-level changes, Floridan aquifer system, USA

Three-dimensional density-dependent flow and transport modeling of the Floridan aquifer system, USA shows that current chloride concentrations are not in equilibrium with current sea level and, second, that the geometric configuration of the aquifer has a significant effect on system responses. The modeling shows that hydraulic head equilibrates first, followed by temperatures, and then by chloride concentrations. The model was constructed using a modified version of SUTRA capable of simulating multi-species heat and solute transport, and was compared to pre-development conditions using hydraulic heads, chloride concentrations, and temperatures from 315 observation wells. Three hypothetical, sinusoidal sea-level changes occurring over 100,000 years were used to evaluate how the simulated aquifer responds to sea-level changes. Model results show that hydraulic head responses lag behind sea-level changes only where the Miocene Hawthorn confining unit is thick and represents a significant restriction to flow. Temperatures equilibrate quickly except where the Hawthorn confining unit is thick and the duration of the sea-level event is long (exceeding 30,000 years). Response times for chloride concentrations to equilibrate are shortest near the coastline and where the aquifer is unconfined; in contrast, chloride concentrations do not change significantly over the 100,000-year simulation period where the Hawthorn confining unit is thick.

Hydrogeology Journal↗
USGS2025

Total uncertainty quantification in inverse solutions with deep learning surrogate models

We propose an approximate Bayesian method for quantifying the total uncertainty in inverse partial differential equation (PDE) solutions obtained with machine learning surrogate models, including operator learning models. The proposed method accounts for uncertainty in the observations, PDE, and surrogate models. First, we use the surrogate model to formulate a minimization problem in the reduced space for the maximum a posteriori (MAP) inverse solution. Then, we randomize the MAP objective function and obtain samples of the posterior distribution by minimizing different realizations of the objective function. We test the proposed framework by comparing it with the iterative ensemble smoother and deep ensembling methods for a nonlinear diffusion equation with an unknown space-dependent diffusion coefficient. Among other applications, this equation describes the flow of groundwater in an unconfined aquifer. Depending on the training dataset and ensemble sizes, the proposed method provides similar or more descriptive posteriors of the parameters and states than the iterative ensemble smoother method. Deep ensembling underestimates uncertainty and provides less-informative posteriors than the other two methods. Our results show that, despite inherent uncertainty, surrogate models can be used for parameter and state estimation as an alternative to the inverse methods relying on (more accurate) numerical PDE solvers.

Journal of Computational Physics↗
USGS2012

Rock-forming metals and Pb in modern Alaskan snow

Metal concentrations in annual and subannual increments of snowpack from the accumulation zone of a south central Alaska glacier indicate that the deposition of Pb with and upon snow is decoupled from that of rock dusts. Rock dusts accumulate, apparently as dry deposition, on the topmost, exposed surfaces of snowpacks in spring and summer, whereas Pb does not. Pb concentration is elevated throughout the latest one third of an annual snowpack, whereas that of rock dusts is not. For whole-year snowpacks, there is a generally sympathetic relationship among concentration of Pb, concentration of rock dust, degree of dominance of rock dusts over ocean solutes, and ferromagnesian character of the rock dusts; however, the fractional abundance of Pb in whole year samples may decrease when rock dust masses become large and/or when rock dusts dominate most strongly over salts. The metal suite chosen to characterize rock dusts and to distinguish them from ocean solutes gives detailed information about rock type of dust source areas and about the nature of the degraded rock products that are taken up, transported, and deposited by the atmosphere. Rock dusts are present at concentrations of only about 300 nanograms (ng) of dust per gram of snow in the Alaskan snowpacks. Concentrations of Pb in the Alaska snow samples are moderate, ranging from 0.1 to 0.3 ng Pb/g snow. This contrasts with larger Pb concentrations of 0.4 to 0.9 ng Pb/g snow in whole-year snowpack samples from the Sierra Nevada, California; with similar to smaller concentrations from north and south Greenland of about 0.04 ng Pb/g snow or less, and about 0.2 ng Pb/g snow or less, respectively, and with much smaller concentrations from Antarctica, now believed to range from a minimum of about 0.001 to a maximum of 0.005 (or 0.01) ng Pb/g snow.

Journal of Geophysical Research D: Atmospheres↗
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↗
USGS1995

Fluxes of water and solute in a coastal wetland sediment. 2. Effect of macropores on solute exchange with surface water

Chloride was highly concentrated relative to seawater in matrix porewater but was comparatively dilute in macropores. Concentration differences in pore-size classes declined with depth until indistinguishable below 10 cm. The segregated chloride distribution can be explained if recharge to the sediment occurred by downward infiltration in macropores and discharge occurred by an upward flux in matrix pores to satisfy evapotranspiration. Without disturbance by the downward infiltration flux in macropores, upward advection of chloride in matrix pores and evapoconcentration increased chloride concentrations in matrix pores to a level well above the concentration in seawater. The resulting high concentrations of chloride in matrix pores induced a large diffusive efflux of chloride into surface water that was sufficient to balance new input of chloride by infiltration of seawater in macropores (0.085 mmol Cl cm -2 day-1). Transport models that were constrained by water balance measurements at the field site explained both the exponential form of the vertical distribution of chloride in matrix pores and the rate of change in storage of chloride in sediment porewater over a one month period.

Journal of Hydrology↗
USGS2017

Using surface velocities to calculate ice thickness and bed topography: A case study at Columbia Glacier, Alaska, USA

Information about glacier volume and ice thickness distribution is essential for many glaciological applications, but direct measurements of ice thickness can be difficult and costly. We present a new method that calculates ice thickness via an estimate of ice flux. We solve the familiar continuity equation between adjacent flowlines, which decreases the computational time required compared to a solution on the whole grid. We test the method on Columbia Glacier, a large tidewater glacier in Alaska, USA, and compare calculated and measured ice thicknesses, with favorable results. This shows the potential of this method for estimating ice thickness distribution of glaciers for which only surface data are available. We find that both the mean thickness and volume of Columbia Glacier were approximately halved over the period 1957–2007, from 281m to 143 m, and from 294 km 3 to 134 km 3 , respectively. Using bedrock slope and considering how waves of thickness change propagate through the glacier, we conduct a brief analysis of the instability of Columbia Glacier, which leads us to conclude that the rapid portion of the retreat may be nearing an end.

Alaska↗
USGS2012

Modeling crater topography and albedo from monoscopic Viking orbiter images: 1. Methodology

A new photoclinometric technique for extraction of topographic data from single planetary images is presented that overcomes many previous limitations of photoclinometry. The procedure fully compensates for oblique viewing geometry prevalent in spacecraft images. Albedo variations have been one of the most serious obstacles in the application of photoclinometry to planetary surfaces. This problem is overcome in the topographic solution by simultaneously utilizing brightness data from a pair of profiles; both segments are assumed to have the same topographic and albedo variations along their lengths. Profile directions are chosen where the orientation of downslope or upslope is obvious, thus resolving a major ambiguity in photoclinometry. This requirement is particularly easy to satisfy for craters and not very difficult for many irregular features. An additional procedure is presented that eliminates even the requirement of topographic symmetry along the pair of profiles. If two profiles have the same relief but their shapes are very different, another method can be used in an iterative process to derive topographic profiles; however, this procedure does, require that the albedo not vary along the profiles. Test results indicate that both procedures have an accuracy and precision of approximately 2° for slopes of typical bowl-shaped craters, which translates to approximately 5% for depths.

Journal of Geophysical Research Solid Earth↗
USGS2018

Hydrothermal discharge from the El Tatio basin, Atacama, Chile

El Tatio in northern Chile is one of the best-studied geothermal fields in South America. However, there remain open questions about the mass and energy budgets, water recharge rates and residence time in the subsurface, origin of dissolved solutes, and processes affecting the phase and chemical composition of groundwater and surface water. We measured and sampled surface manifestations of the geothermal system (geysers perpetual spouters, mud pools/volcanoes, and non-eruptive hot springs) and meteoric water. From the isotopic composition we infer that the thermal water has a meteoric origin that is different from the composition of local meteoric water. The absence of detectable tritium in thermal waters indicates that most of the recharge occurred pre-1950. Boiling and steam separation from the deep reservoir appear to be the main subsurface processes affecting the thermal fluids. A large amount of heat is lost to the atmosphere by evaporation from surface water and by steam emitted from erupting geysers. Using the chloride inventory method, we estimate thermal water discharge to be 218 to 234 L/s, and the advective heat flow to be 120 to 170 MW.

El Tatio Basin↗