Search USGS⌕ Search

SEARCH · Search USGS

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 325 records · Page 18Linked to original sources
USGS2019

A combinatorial approach to determine earthquake magnitude distributions on a variable slip-rate fault

Combinatorial methods are used to determine the spatial distribution of earthquake magnitudes on a fault whose slip rate varies along strike. Input to the problem is a finite sample of earthquake magnitudes that span 5 kyr drawn from a truncated Pareto distribution. The primary constraints to the problem are maximum and minimum values around the target slip-rate function indicating where feasible solutions can occur. Two methods are used to determine the spatial distribution of earthquakes: integer programming and the greedy-sequential algorithm. For the integer-programming method, the binary decision vector includes all possible locations along the fault where each earthquake can occur. Once a set of solutions that satisfy the constraints is found, the cumulative slip misfit on the fault is globally minimized relative to the target slip-rate function. The greedy algorithm sequentially places earthquakes to locally optimize slip accumulation. As a case study, we calculate how earthquakes are distributed along the megathrust of the Nankai subduction zone, in which the slip rate varies significantly along strike. For both methods, the spatial distribution of magnitudes depends on slip rate, except for the largest magnitude earthquakes that span multiple sections of the fault. The greedy-sequential algorithm, previously applied to this fault (Parsons et al., 2012), tends to produce smoother spatial distributions and fewer lower magnitude earthquakes in the low slip-rate section of the fault compared to the integer-programming method. Differences in results from the two methods relate to how much emphasis is placed on minimizing the misfit to the target slip rate (integer programming) compared to finding a solution within the slip-rate constraints (greedy sequential). Specifics of the spatial distribution of magnitudes also depend on the shape of the target slip-rate function: i.e. stepped at the section boundaries versus a smooth function. This study isolates the effects of slip-rate variation along a single fault in determining the spatial distribution of earthquake magnitudes, helping to better interpret results from more complex, interconnected fault systems.

Geophysical Journal International↗
USGS1956

Ceroid in fish

Since the original description of ceroid in rats, many papers have appeared on the etiology and characteristics of this pigment. It was first seen as a yellow, granular pigment in hematoxylin and eosin sections of the cirrhotic livers of choline deficient rats. The pigment was more fully characterized by Endicott and Lillie, and additional stainging reactions were summarized recently by Lillie. The pigment is sudanophilic in paraffin sections, acid-fast, basophilic, isotropic, iron negative, and highly resistant to solution in water, alcohol, fat solvents, and dilute aqueous acids and alkalis. It is stained by Mallory's hemofuscin stain and Weigert's myelin stain. It reduces osmium tetraoxide and diamine silver carbonate but not ferric ferricyanide. The Gmelin reactions for bile pigments is negative. It has a greenish yellow fluorescence at 3650-3660 Å. It is Schiff positive with or without antecedent diastase digestion after performic or periodic acid oxidation.

The American Journal of Pathology↗
USGS2006

Post-seismic relaxation following the great 2004 Sumatra-Andaman earthquake on a compressible self-gravitating Earth

he Mw ??? 9.0 2004 December 26 Sumatra-Andaman and Mw =8.7 2005 March 28 Nias earthquakes, which collectively ruptured approximately 1800 km of the Andaman and Sunda subduction zones, are expected to be followed by vigorous viscoelastic relaxation involving both the upper and lower mantle. Because of these large spatial dimensions it is desirable to fully account for gravitational coupling effects in the relaxation process. We present a stable method of computing relaxation of a spherically-stratified, compressible and self-gravitating viscoelastic Earth following an impulsive moment release event. The solution is cast in terms of a spherical harmonic expansion of viscoelastic normal modes. For simple layered viscoelastic models, which include a low-viscosity oceanic asthenosphere, we predict substantial post-seismic effects over a region several 100s of km wide surrounding the eastern Indian Ocean. We compare observed GPS time-series from ten regional sites (mostly in Thailand and Indonesia), beginning in 2004 December, with synthetic time-series that include the coseismic and post-seismic effects of the 2004 December 26 and 2005 March 28 earthquakes. A viscosity structure involving a biviscous (Burgers body) rheology in the asthenosphere explains the pattern and amplitude of post-seismic offsets remarkably well. ?? 2006 The Authors Journal compilation ?? 2006 RAS.

Geophysical Journal International↗
USGS1991

A model for hydrostatic consolidation of Pierre shale

This paper presents closed-form solutions for consolidation of transversely isotropic porous media under hydrostatic stress. The solutions are applied to model the time variation of pore pressure, volume strain and strains parallel and normal to bedding, and to obtain coefficients of consolidation and permeability, as well as other properties, and the bulk modulus resulting from hydrostatic consolidation of Pierre shale. It is found that the coefficients consolidation and permeability decrease and the bulk moduli increase with increasing confining pressure, reflecting the closure of voids in the rock. ?? 1991.

International Journal of Rock Mechanics and Mining↗
USGS2001

Model coupling intraparticle diffusion/sorption, nonlinear sorption, and biodegradation processes

Diffusion, sorption and biodegradation are key processes impacting the efficiency of natural attenuation. While each process has been studied individually, limited information exists on the kinetic coupling of these processes. In this paper, a model is presented that couples nonlinear and nonequilibrium sorption (intraparticle diffusion) with biodegradation kinetics. Initially, these processes are studied independently (i.e., intraparticle diffusion, nonlinear sorption and biodegradation), with appropriate parameters determined from these independent studies. Then, the coupled processes are studied, with an initial data set used to determine biodegradation constants that were subsequently used to successfully predict the behavior of a second data set. The validated model is then used to conduct a sensitivity analysis, which reveals conditions where biodegradation becomes desorption rate-limited. If the chemical is not pre-equilibrated with the soil prior to the onset of biodegradation, then fast sorption will reduce aqueous concentrations and thus biodegradation rates. Another sensitivity analysis demonstrates the importance of including nonlinear sorption in a coupled diffusion/sorption and biodegradation model. While predictions based on linear sorption isotherms agree well with solution concentrations, for the conditions evaluated this approach overestimates the percentage of contaminant biodegraded by as much as 50%. This research demonstrates that nonlinear sorption should be coupled with diffusion/sorption and biodegradation models in order to accurately predict bioremediation and natural attenuation processes. To our knowledge this study is unique in studying nonlinear sorption coupled with intraparticle diffusion and biodegradation kinetics with natural media.

Journal of Contaminant Hydrology↗
USGS2009

Thermal conductivity of hydrate-bearing sediments

A thorough understanding of the thermal conductivity of hydrate-bearing sediments is necessary for evaluating phase transformation processes that would accompany energy production from gas hydrate deposits and for estimating regional heat flow based on the observed depth to the base of the gas hydrate stability zone. The coexistence of multiple phases (gas hydrate, liquid and gas pore fill, and solid sediment grains) and their complex spatial arrangement hinder the a priori prediction of the thermal conductivity of hydrate-bearing sediments. Previous studies have been unable to capture the full parameter space covered by variations in grain size, specific surface, degree of saturation, nature of pore filling material, and effective stress for hydrate-bearing samples. Here we report on systematic measurements of the thermal conductivity of air dry, water- and tetrohydrofuran (THF)-saturated, and THF hydrate–saturated sand and clay samples at vertical effective stress of 0.05 to 1 MPa (corresponding to depths as great as 100 m below seafloor). Results reveal that the bulk thermal conductivity of the samples in every case reflects a complex interplay among particle size, effective stress, porosity, and fluid-versus-hydrate filled pore spaces. The thermal conductivity of THF hydrate–bearing soils increases upon hydrate formation although the thermal conductivities of THF solution and THF hydrate are almost the same. Several mechanisms can contribute to this effect including cryogenic suction during hydrate crystal growth and the ensuing porosity reduction in the surrounding sediment, increased mean effective stress due to hydrate formation under zero lateral strain conditions, and decreased interface thermal impedance as grain-liquid interfaces are transformed into grain-hydrate interfaces.

Journal of Geophysical Research B: Solid Earth↗
USGS1990

Plane-strain shear dislocations moving steadily in linear elastic diffusive solids

This paper derives the stress and pore pressure fields induced by a plane-strain shear (gliding edge) dislocation moving steadily at a constant speed V in a linear elastic, fluid-infiltrated (Biot) solid. Solutions are obtained for the limiting cases in which the plane containing the moving dislocation (y = 0) is permeable and impermeable to the diffusing species. Although the solutions for the permeable and impermeable planes are required to agree with each other and with the ordinary elastic solution in the limits of V = 0 (corresponding to drained response) and V = ∞ (corresponding to undrained response), the stress and pore pressure fields differ considerably for finite nonzero velocities. For the dislocation on the impermeable plane, the pore pressure is discontinuous on y = 0 and attains values which are equal in magnitude and opposite in sign as y = 0 is approached from above and below. The solution reveals the surprising result that the pore pressure on the impermeable plane is zero everywhere behind the moving dislocation (x < 0). For the dislocation on the permeable plane, the pore pressure is zero on y = 0 and attains its maximum at about (2c/V, 2c/V) where c is the diffusivity, and the origin of the coordinate system coincides with the dislocation. For the impermeable plane, the largest pore pressure change occurs at the origin.

Journal of Applied Mechanics↗
USGS2008

Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System

Systematic improvements in algorithmic design of regional ocean circulation models have led to significant enhancement in simulation ability across a wide range of space/time scales and marine system types. As an example, we briefly review the Regional Ocean Modeling System, a member of a general class of three-dimensional, free-surface, terrain-following numerical models. Noteworthy characteristics of the ROMS computational kernel include: consistent temporal averaging of the barotropic mode to guarantee both exact conservation and constancy preservation properties for tracers; redefined barotropic pressure-gradient terms to account for local variations in the density field; vertical interpolation performed using conservative parabolic splines; and higher-order, quasi-monotone advection algorithms. Examples of quantitative skill assessment are shown for a tidally driven estuary, an ice-covered high-latitude sea, a wind- and buoyancy-forced continental shelf, and a mid-latitude ocean basin. The combination of moderate-order spatial approximations, enhanced conservation properties, and quasi-monotone advection produces both more robust and accurate, and less diffusive, solutions than those produced in earlier terrain-following ocean models. Together with advanced methods of data assimilation and novel observing system technologies, these capabilities constitute the necessary ingredients for multi-purpose regional ocean prediction systems.

Journal of Computational Physics↗
USGS2000

Fracture process zone in granite

In uniaxial compression tests performed on Aue granite cores (diameter 50 mm, length 100 mm), a steel loading plate was used to induce the formation of a discrete shear fracture. A zone of distributed microcracks surrounds the tip of the propagating fracture. This process zone is imaged by locating acoustic emission events using 12 piezoceramic sensors attached to the samples. Propagation velocity of the process zone is varied by using the rate of acoustic emissions to control the applied axial force. The resulting velocities range from 2 mm/s in displacement-controlled tests to 2 μm/s in tests controlled by acoustic emission rate. Wave velocities and amplitudes are monitored during fault formation. P waves transmitted through the approaching process zone show a drop in amplitude of 26 dB, and ultrasonic velocities are reduced by 10%. The width of the process zone is ∼9 times the grain diameter inferred from acoustic data but is only 2 times the grain size from optical crack inspection. The process zone of fast propagating fractures is wider than for slow ones. The density of microcracks and acoustic emissions increases approaching the main fracture. Shear displacement scales linearly with fracture length. Fault plane solutions from acoustic events show similar orientation of nodal planes on both sides of the shear fracture. The ratio of the process zone width to the fault length in Aue granite ranges from 0.01 to 0.1 inferred from crack data and acoustic emissions, respectively. The fracture surface energy is estimated from microstructure analysis to be ∼2 J. A lower bound estimate for the energy dissipated by acoustic events is 0.1 J.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

Source characterization for an explosion during the 2009 eruption of Redoubt Volcano from very-long-period seismic waves

The 2009 eruption of Redoubt produced several very-long-period (VLP) signals associated with explosions. We invert for the source location and mechanism of an explosion at Redoubt volcano using waveform methods applied to broadband recordings. Such characterization of the source carries information on the geometry of the conduit and the physics of the explosion process. Inversions are carried out assuming the volcanic source can be modeled as a point source, with mechanisms described by a) a set of 3 orthogonal forces, b) a moment tensor consisting of force couples, and c) both forces and moment tensor components. We find that the source of the VLP seismic waves during the explosion is well-described by either a combined moment/force source located northeast of the crater and at an elevation of 1.6 km ASL or a moment source at an elevation of 800 m to the southwest of the crater. The moment tensors for the solutions with moment and force and moment-only share similar characteristics. The source time functions for both moment tensors begin with inflation (pressurization) and execute two cycles of deflation-reinflation (depressurization–repressurization). Although the moment/force source provides a better fit to the data, we find that owing to the limited coverage of the broadband stations at Redoubt the moment-only source is the more robust and reliable solution. Based on the moment-only solution, we estimate a volume change of 19,000 m 3 and a pressure change of 7 MPa in a dominant sill and an out-of-phase volume change of 5000 m 3 and pressure change of 1.8 MPa in a subdominant dike at the source location. These results shed new light on the magmatic plumbing system beneath Redoubt and complement previous studies on Vulcanian explosions at other volcanoes.

Alaska↗
USGS1986

Predicting two-dimensional steady-state soil freezing fronts using the CVBEM

The complex variable boundary element method (CVBEM) is used instead of a real variable boundary element method due to the available modeling error evaluation techniques developed. The modeling accuracy is evaluated by the model-user in the determination of an approximative boundary upon which the CVBEM provides an exact solution. Although inhomogeneity (and anisotropy) can be included in the CVBEM model, the resulting fully populated matrix system quickly becomes large. Therefore in this paper, the domain is assumed homogeneous and isotropic except for differences in frozen and thawed conduction parameters on either side of the freezing front. The example problems presented were obtained by use of a popular 64K microcomputer (the current version of the program used in this study has the capacity to accommodate 30 nodal points).

Journal of Heat Transfer↗
USGS1987

Comparative efficacy of 16 anesthetic chemicals on rainbow trout

Presently there are no legally registered fish anesthetics that allow for the release of fish or use of the fish for food soon after they have been anesthetized. MS-222 (tricaine), the only anesthetic registered for use on fish in the United States, cannot be used within 21 d of harvesting the fish for food. As the start in a search for an anesthetic that can be used with little or no withdrawal period, we tested the efficacy of 16 chemicals as anesthetics on rainbow trout Salmo gairdneri. Efficacy was defined by the fish (1) becoming handleable (quiet enough to be manipulated and handled readily) in 3 min or less, (2) recovering in 10 min or less, and (3) showing no mortality after 15 min in the anesthetic solution. Four chemicals--MS-222, quinaldine sulfate, benzocaine, and 2-phenoxyethanol--met these criteria for efficacy. Chemicals that yielded excessive induction or recovery times or caused excessive mortality were methylpentynol, chlorobutanol, etomidate, metomidate, Piscaine, propanidid, carbon dioxide, nicotine, salt, Halothane, Metofane, and Biotal. Because carbon dioxide leaves no residues and requires no withdrawal period, it may be an acceptable alternative for fishery workers who can tolerate somewhat shallower anesthesia and longer induction and recovery times.

North American Journal of Fisheries Management↗
USGS1977

Measuring total antimony in geothermal waters by flame atomic absorption spectrometry

A flame atomic absorption procedure utilizing an electrodeless discharge lamp is described for determining total solute antimony in silica-rich geothermal waters. Following NaNO 2 oxidation of Sb +3 , SbCl 5 is extracted from 6 N HCl solution by using MIBK (methyl isobutyl ketone) ; silica in the organic layer is removed by centrifugation. The analytical detection limit for antimony is about 6 micrograms per liter in the original sample; the coefficient of variation is 4 percent at the 250 μ g/L Sb level. Comparable total antimony levels were found in sample splits which had been either filtered and acidified (to pH less than 1.5) with HNO 3 or left untreated at the time of sample collection.

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

Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M = 9.2 Sumatra earthquake

The 2004 M=9.2 Sumatra-Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ???1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress-strain relation; an aspherical perturbation in viscoelastic structure; a 'static'mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra-Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation. ?? Journal compilation ?? 2008 RAS.

Geophysical Journal International↗
USGS2013

Nitrate in watersheds: straight from soils to streams?

Human activities are rapidly increasing the global supply of reactive N and substantially altering the structure and hydrologic connectivity of managed ecosystems. There is long-standing recognition that N must be removed along hydrologic flowpaths from uplands to streams, yet it has proven difficult to assess the generality of this removal across ecosystem types, and whether these patterns are influenced by land-use change. To assess how well upland nitrate (NO 3 - ) loss is reflected in stream export, we gathered information from >50 watershed biogeochemical studies that reported nitrate concentrations ([NO 3 - ]) for stream water and for either upslope soil solution or groundwater NO 3 - to examine whether stream export of NO3- accurately reflects upland NO 3 - losses. In this dataset, soil solution and streamwater [NO 3 - ] were correlated across 40 undisturbed forest watersheds, with streamwater [NO 3 - ] typically half (median = 50%) soil solution [NO 3 - ]. A similar relationship was seen in 10 disturbed forest watersheds. However, for 12 watersheds with significant agricultural or urban development, the intercept and slope were both significantly higher than the relationship seen in forest watersheds. Differences in concentration between soil solution or groundwater and stream water may be attributed to biological uptake, microbial processes including denitrification, and/or preferential flow routing. The results of this synthesis are consistent with the hypotheses that undisturbed watersheds have a significant capacity to remove nitrate after it passes below the rooting zone and that land use changes tend to alter the efficiency or the length of watershed flowpaths, leading to reductions in nitrate removal and increased stream nitrate concentrations.

Journal of Geophysical Research G: Biogeosciences↗
USGS1987

Estimating threshold values for the land disposal of organic solvent-contaminated wastes

This paper describes conceptual and mathematical modeling approaches to evaluate environmentally suitable waste management practices that prevent ground water pollution from landfill disposal of solid wastes contaminated with organic solvents. Threshold values for disposal of solid wastes contaminated with acetone, nitrobenzene, chlorobenzene, and carbon tetrachloride were estimated using a disposal site model using a two-dimensional advection-,dispersion solute transport computer code. For example, the threshold value for disposal of carbon tetrachloride at the site model was approximately 600 g/acre, while the value for acetone was limited only by the characteristics of the waste to prevent free drainage of liquid solvent. The study indicated that the toxicity, mass loading rate, and the mass of solvent initially disposed are critical parameters determining the relative success of a given site to attenuate contaminants to environmentally acceptable levels.

Journal of Hazardous Materials↗
USGS2014

Bayesian historical earthquake relocation: an example from the 1909 Taipei earthquake

Locating earthquakes from the beginning of the modern instrumental period is complicated by the fact that there are few good-quality seismograms and what traveltimes do exist may be corrupted by both large phase-pick errors and clock errors. Here, we outline a Bayesian approach to simultaneous inference of not only the hypocentre location but also the clock errors at each station and the origin time of the earthquake. This methodology improves the solution for the source location and also provides an uncertainty analysis on all of the parameters included in the inversion. As an example, we applied this Bayesian approach to the well-studied 1909 M w 7 Taipei earthquake. While our epicentre location and origin time for the 1909 Taipei earthquake are consistent with earlier studies, our focal depth is significantly shallower suggesting a higher seismic hazard to the populous Taipei metropolitan area than previously supposed.

Taipei↗
USGS2011

Role of back diffusion and biodegradation reactions in sustaining an MTBE/TBA plume in alluvial media

A methyl tert-butyl ether (MTBE) / tert-butyl alcohol (TBA) plume originating from a gasoline spill in late 1994 at Vandenberg Air Force Base (VAFB) persisted for over 15 years within 200 feet of the original spill source. The plume persisted until 2010 despite excavation of the tanks and piping within months after the spill and excavations of additional contaminated sediments from the source area in 2007 and 2008. The probable history of MTBE concentrations along the plume centerline at its source was estimated using a wide variety of available information, including published details about the original spill, excavations and monitoring by VAFB consultants, and our own research data. Two-dimensional reactive transport simulations of MTBE along the plume centerline were conducted for a 20-year period following the spill. These analyses suggest that MTBE diffused from the thin anaerobic aquifer into the adjacent anaerobic silts and transformed to TBA in both aquifer and silt layers. The model reproduces the observation that after 2004 TBA was the dominant solute, diffusing back out of the silts into the aquifer and sustaining plume concentrations much longer than would have been the case in the absence of such diffusive exchange. Simulations also suggest that aerobic degradation of MTBE or TBA at the water table in the overlying silt layer significantly affected concentrations of MTBE and TBA by limiting the chemical mass available for back diffusion to the aquifer.

California↗