USGS Science⌕ Search

SEARCH · USGS Science

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.

637 records · Page 36Linked to original sources
USGS1977

Application of gold compositional analyses to mineral exploration in the United States

Native gold is a mineral composed of Au, Ag and Cu in solid solution and it usually contains one or more trace metals as lattice impurities, as mineral inclusions, in grain boundaries or in surface coatings. Alloy proportions of Au, Ag and Cu, together with certain other elements, can be thought of as constituting a gold “signature”. Gold is associated with a great variety of ore deposits and has characteristic signatures for each of several types of ore deposits. Signatures for gold derived from igneous-metamorphic, hypothermal, mesothermal and epithermal deposits reflect conditions of ore formation by their content of Ag, Cu and characteristic associated elements. At higher temperatures of ore formation, gold has low Ag and high Cu content, and Bi and Pb are the most abundant trace elements. But at lower temperatures of ore formation, Ag is high, Cu is low, and Pb is the most abundant trace element. The same trend in gold signatures is observable in gold mining districts, such as Central City, Colorado, where zoning as shown by mineral assemblages indicates ore deposition at progressively lower temperatures as the distance from a central high-temperature zone increases. The signatures of gold may be useful in searching for porphyry Cu deposits. Signatures from Butte (Montana), Mineral Park (Arizona) and Cala Abajo (Puerto Rico), on the basis of limited sampling, are similar and distinctive. They are characterized by a similar assemblage of trace elements and are relatively high in both Ag and Cu. Another application of gold compositional data is in tracing placer gold to its bedrock source. For example, the Ag content of placer gold in the Tarryall district of Colorado differed from that of nearly all of the bedrock sources of gold found by early prospectors. However, one lightly prospected area peripheral to the Tertiary quartz monzonite stock at Montgomery Gulch contains gold with a Ag content similar to that of the placer gold. This area is the most likely source of the gold in the productive placers and may be a potential exploration target. Gold signatures may be useful in prospecting for metals other than gold. Several metals of low crustal abundance — notably Sn, W, Mo and the Pt group metals — are detected in analyses of some gold samples and may indicate economic deposits of these metals.

Journal of Geochemical Exploration↗
USGS2014

Seismological analyses of the 2010 March 11, Pichilemu, Chile Mw 7.0 and Mw 6.9 coastal intraplate earthquakes

On 2010 March 11, a sequence of large, shallow continental crust earthquakes shook central Chile. Two normal faulting events with magnitudes around M w 7.0 and M w 6.9 occurred just 15 min apart, located near the town of Pichilemu. These kinds of large intraplate, inland crustal earthquakes are rare above the Chilean subduction zone, and it is important to better understand their relationship with the 2010 February 27, M w 8.8, Maule earthquake, which ruptured the adjacent megathrust plate boundary. We present a broad seismological analysis of these earthquakes by using both teleseismic and regional data. We compute seismic moment tensors for both events via a W-phase inversion, and test sensitivities to various inversion parameters in order to assess the stability of the solutions. The first event, at 14 hr 39 min GMT, is well constrained, displaying a fault plane with strike of N145°E, and a preferred dip angle of 55°SW, consistent with the trend of aftershock locations and other published results. Teleseismic finite-fault inversions for this event show a large slip zone along the southern part of the fault, correlating well with the reported spatial density of aftershocks. The second earthquake (14 hr 55 min GMT) appears to have ruptured a fault branching southward from the previous ruptured fault, within the hanging wall of the first event. Modelling seismograms at regional to teleseismic distances (Δ > 10°) is quite challenging because the observed seismic wave fields of both events overlap, increasing apparent complexity for the second earthquake. We perform both point- and extended-source inversions at regional and teleseismic distances, assessing model sensitivities resulting from variations in fault orientation, dimension, and hypocentre location. Results show that the focal mechanism for the second event features a steeper dip angle and a strike rotated slightly clockwise with respect to the previous event. This kind of geological fault configuration, with secondary rupture in the hanging wall of a large normal fault, is commonly observed in extensional geological regimes. We propose that both earthquakes form part of a typical normal fault diverging splay, where the secondary fault connects to the main fault at depth. To ascertain more information on the spatial and temporal details of slip for both events, we gathered near-fault seismological and geodetic data. Through forward modelling of near-fault synthetic seismograms we build a kinematic k −2 earthquake source model with spatially distributed slip on the fault that, to first-order, explains both coseismic static displacement GPS vectors and short-period seismometer observations at the closest sites. As expected, the results for the first event agree with the focal mechanism derived from teleseismic modelling, with a magnitude M w 6.97. Similarly, near-fault modelling for the second event suggests rupture along a normal fault, M w 6.90, characterized by a steeper dip angle (dip = 74°) and a strike clockwise rotated (strike = 155°) with respect to the previous event.

Pichilemu↗
USGS2010

Vadose zone attenuation of organic compounds at a crude oil spill site - Interactions between biogeochemical reactions and multicomponent gas transport

Contaminant attenuation processes in the vadose zone of a crude oil spill site near Bemidji, MN have been simulated with a reactive transport model that includes multicomponent gas transport, solute transport, and the most relevant biogeochemical reactions. Dissolution and volatilization of oil components, their aerobic and anaerobic degradation coupled with sequential electron acceptor consumption, ingress of atmospheric O 2 , and the release of CH 4 and CO 2 from the smear zone generated by the floating oil were considered. The focus of the simulations was to assess the dynamics between biodegradation and gas transport processes in the vadose zone, to evaluate the rates and contributions of different electron accepting processes towards vadose zone natural attenuation, and to provide an estimate of the historical mass loss. Concentration distributions of reactive (O 2 , CH 4 , and CO 2 ) and non-reactive (Ar and N 2 ) gases served as key constraints for the model calibration. Simulation results confirm that as of 2007, the main degradation pathway can be attributed to methanogenic degradation of organic compounds in the smear zone and the vadose zone resulting in a contaminant plume dominated by high CH 4 concentrations. In accordance with field observations, zones of volatilization and CH 4 generation are correlated to slightly elevated total gas pressures and low partial pressures of N 2 and Ar, while zones of aerobic CH 4 oxidation are characterized by slightly reduced gas pressures and elevated concentrations of N 2 and Ar. Diffusion is the most significant transport mechanism for gases in the vadose zone; however, the simulations also indicate that, despite very small pressure gradients, advection contributes up to 15% towards the net flux of CH 4 , and to a more limited extent to O 2 ingress. Model calibration strongly suggests that transfer of biogenically generated gases from the smear zone provides a major control on vadose zone gas distributions and vadose zone carbon balance. Overall, the model was successful in capturing the complex interactions between biogeochemical reactions and multicomponent gas transport processes. However, despite employing a process-based modeling approach, honoring observed parameter ranges, and generally obtaining good agreement between field observations and model simulations, accurate quantification of natural attenuation rates remains difficult. The modeling results are affected by uncertainties regarding gas phase saturations, tortuosities, and the magnitude of CH 4 and CO 2 flux from the smear zone. These findings highlight the need to better delineate gas fluxes at the model boundaries, which will help constrain contaminant degradation rates, and ultimately source zone longevity.

Minnesota↗
USGS1998

Earthquake triggering by transient and static deformations

Observational evidence for both static and transient near-field and far-field triggered seismicity are explained in terms of a frictional instability model, based on a single degree of freedom spring-slider system and rate- and state-dependent frictional constitutive equations. In this study a triggered earthquake is one whose failure time has been advanced by Δ t (clock advance) due to a stress perturbation. Triggering stress perturbations considered include square-wave transients and step functions, analogous to seismic waves and coseismic static stress changes, respectively. Perturbations are superimposed on a constant background stressing rate which represents the tectonic stressing rate. The normal stress is assumed to be constant. Approximate, closed-form solutions of the rate-and-state equations are derived for these triggering and background loads, building on the work of Dieterich [1992, 1994]. These solutions can be used to simulate the effects of static and transient stresses as a function of amplitude, onset time t 0 , and in the case of square waves, duration. The accuracies of the approximate closed-form solutions are also evaluated with respect to the full numerical solution and t 0 . The approximate solutions underpredict the full solutions, although the difference decreases as t 0 approaches the end of the earthquake cycle. The relationship between Δ t and t 0 differs for transient and static loads: a static stress step imposed late in the cycle causes less clock advance than an equal step imposed earlier, whereas a later applied transient causes greater clock advance than an equal one imposed earlier. For equal Δ t , transient amplitudes must be greater than static loads by factors of several tens to hundreds depending on t 0 . We show that the rate-and-state model requires that the total slip at failure is a constant, regardless of the loading history. Thus a static load applied early in the cycle, or a transient applied at any time, reduces the stress at the initiation of failure, whereas static loads that are applied sufficiently late raise it. Rate-and-state friction predictions differ markedly from those based on Coulomb failure stress changes (ΔCFS) in which Δ t equals the amplitude of the static stress change divided by the background stressing rate. The ΔCFS model assumes a stress failure threshold, while the rate-and-state equations require a slip failure threshold. The complete rate-and-state equations predict larger Δ t than the ΔCFS model does for static stress steps at small t 0 , and smaller Δ t than the ΔCFS model for stress steps at large t 0 . The ΔCFS model predicts nonzero Δ t only for transient loads that raise the stress to failure stress levels during the transient. In contrast, the rate-and-state model predicts nonzero Δ t for smaller loads, and triggered failure may occur well after the transient is finished. We consider heuristically the effects of triggering on a population of faults, as these effects might be evident in seismicity data. Triggering is manifest as an initial increase in seismicity rate that may be followed by a quiescence or by a return to the background rate. Available seismicity data are insufficient to discriminate whether triggered earthquakes are “new” or clock advanced. However, if triggering indeed results from advancing the failure time of inevitable earthquakes, then our modeling suggests that a quiescence always follows transient triggering and that the duration of increased seismicity also cannot exceed the duration of a triggering transient load. Quiescence follows static triggering only if the population of available faults is finite.

Journal of Geophysical Research B: Solid Earth↗
USGS1990

Dolomite dissolution rates and possible Holocene dedolomitization of water-bearing units in the Edwards aquifer, south-central Texas

Rates of dolomite dissolution can be used to test the concept, based on geomorphologic evidence, that a major part of the Edwards aquifer could have formed within the Holocene, a timeframe of approximately 10,000 years. During formation of the aquifer in the Edwards limestone (Cretaceous, Albian) of the Balcones fault zone, dolomite dissolution and porosity development were synchronous and the result of mixing-zone dedolomitization. Initiation of the mixing zone in the early Holocene (∼11,000 years before present) is suggested by the maximum age of formation of major discharge sites that allowed the influx of meteoric water into brine-filled, dolomitic preaquifer units. Dedolomitization, the dissolution of dolomite and net precipitation of calcite, has left aquifer units that are calcitic, and 40 vol.% interconnected pore space. The mass of dolomite missing is obtained by comparison of stratigraphically equivalent altered and unaltered units. One dissolution rate (1.76 × 10 −4 mmol dolomite kg H2O −1 yr −1 ) is determined from this mass, 10 4 yr reaction time, and a log-linear function describing the increase in mass discharge (three orders of magnitude) during aquifer formation. The second estimated dissolution rate is obtained from the mass transfer of dolomite to solution calculated from the increase in magnesium in pore fluids selected from the modern aquifer to represent a typical flowpath during aquifer formation. A reaction time of 10 4 yr for this mass transfer yields a rate of 0.56 × 10 −4 mmol dolomite kg H2O −1 yr −1 Both of these rates are comparable to modern rates of dolomite dissolution (0.3 to 4.5 × 10 −4 mmol dolomite kg H2O −1 yr −1 ) calculated from measured reaction times in the Tertiary Floridan aquifer system in Florida and the Madison aquifer in the Mississippian Madison Limestone of the Northern Great Plains. Similarity of these rates to the estimated paleo-rates of dolomite dissolution supports a 10 4 yr reaction timeframe. Both of these rates are comparable to modern rates of dolomite dissolution (0.3 to 4.5 × 10 −4 mmol dolomite kg H2O −1 yr −1 ) calculated from measured reaction times in the Tertiary Floridan aquifer system in Florida and the Madison aquifer in the Mississippian Madison Limestone of the Northern Great Plains. Similarity of these rates to the estimated paleo-rates of dolomite dissolution supports a 10 4 yr reaction timeframe. The Holocene reaction time also can be compared to a series of reaction times calculated by assuming that the mass of dolomite missing from the Edwards was removed at rates observed in the Floridan and Madison aquifers. These reaction times (for complete removal of dolomite) range from 2700 to 58,500 yr and span the Pleistocene-Holocene boundary. Finally, an estimated dolomite reaction rate during dedolomitization of the Edwards aquifer based on surface area of exposed dolomite [mmol cm −2 s −1 (millimoles per square centimeter per second)] may be approximated from reaction times. This rate is directly a function of the mass of dolomite removed and the surface area exposed per pore volume passing through the rock. The surface area is available from the observed dolomite rhomb size in unaltered rock. The rate of pore fluid movement is obtained from the averaged annual discharge. Rates during formation of the Edwards aquifer calculated from all reaction times range from 10 −13 to 10 −14 mmol dolomite cm −2 s −1 . These rates are faster than rates (10 −18 mmol cm −2 s −1 ), measured in the pure laboratory system, CaMg(CO 3 ) 2 CO 2 H 2 O, but slower than rates determined in an alpine stream study (10 −10 to 10 −11 mmol cm −2 s −1 ) where cold glacial melt water flows over dolostone. Dolomite dissolution rates from both the Edwards and other aquifers support the concept that a major part of the Edwards aquifer could have formed within the Holocene.

Texas↗
USGS1996

Long-term geochemical surveillance of fumaroles at Showa-Shinzan dome, Usu volcano, Japan

This study investigates 31 years of fumarole gas and condensate (trace elements) data from Showa-Shinzan, a dacitic dome-cryptodome complex that formed during the 1943-1945 eruption of Usu volcano. Forty-two gas samples were collected from the highest-temperature fumarole, named A-1, from 1954 (800??C) to 1985 (336??C), and from lower-temperature vents. Condensates were collected contemporaneously with the gas samples, and we reanalyzed ten of these samples, mostly from the A-1 vent, for 32 cations and three anions. Modeling using the thermochemical equilibrium program, SOLVGAS, shows that the gas samples are mild disequilibrium mixtures because they: (a) contain unequilibrated sedimentary CH4 and NH3; (b) have unequilibrated meteoric water; or (c) lost CO, either by air oxidation or by absorption by the sodium hydroxide sampling solution. SOLVGAS also enabled us to restore the samples by removing these disequilibrium effects, and to estimate their equilibrium oxygen fugacities and amounts of S2 and CH4. The restored compositions contain > 98% H2O with minor to trace amounts of CO2, H2, HCl, SO2, HF, H2S, CO, S2 and CH4. We used the restored gas and condensate data to test the hypotheses that these time-series compositional data from the dome's fumaroles provide: (1) sufficient major-gas data to analyze long-term degassing trends of the dome's magma-hydrothermal system without the influence of sampling or contamination effects; (2) independent oxygen fugacity-versus-temperature estimates of the Showa-Shinzan dacite; (3) the order of release of trace elements, especially metals, from magma; and (4) useful information for assessing volcanic hazards. The 1954-1985 restored A-1 gas compositions confirm the first hypothesis because they are sufficient to reveal three long-term degassing trends: (1) they became increasingly H2O-rich with time due to the progressive influx of meteoric water into the dome; (2) their C/S and S/Cl ratios decreased dramatically while their Cl/F ratios stayed roughly constant, indicating the progressive outgassing of less soluble components (F ??? Cl > S > C) from the magma reservoir; and (3) their H2O/H2, CO2/CO and H2S/SO2 ratios increased significantly in concert with equilibrium changes expected for the ??? 500??C temperature drop. When plotted against reciprocal temperature, the restored-gas log oxygen fugacities follow a tight linear trend from < NNO + 0.5 at > 800??C to NNO + 2.5 at ??? 400??C. This trend largely disproves the second hypothesis because the oxygen fugacities for the < 800??C restored gases can only be explained by mixing of hot magmatic gases with ??? 350??C steam from superheated meteoric water. But above 800??C this trend intersects the opposing linear trend for other Usu eruptive products, implying a log oxygen fugacity of -11.45 at 902??C for the Showa-Shinzan magma. The time-series trace-element data also disprove the third hypothesis because rock- and incrustation-particle contaminants in the condensates account for most of the trace-element variation. Nonetheless, highly volatile elements like B and As are relatively unaffected by this particle contamination, and they show similar time-series trends as Cl and F. Finally, except for infrequent sampling around the 1977 Usu eruption, the results generally confirm the fourth hypothesis, since the time-series trends for the major gases and selected trace elements indicate that, with time, the system cooled, degassed and was infiltrated by meteoric water, all of which are positive signs that volcanic activity declined over the 31-year history. This study also suggests that second boiling of shallow magma within and possibly beneath the cryptodome sustained magmatic degassing for at least 20 years after emplacement.

Journal of Volcanology and Geothermal Research↗
USGS2012

Modification of wave-cut and faulting-controlled landforms

From a casual observation that the form of degraded fault scarps resembles the error function, this investigation proceeds through an elementary diffusion equation representation of landform evolution to the application of the resulting equations to the modern topography of scarplike landforms. The morphologic observations can be analyzed either in the form of one or more cross-strike elevation profiles or in the form of the slope-offset plot, a point plot of maximum scarp slope versus scarp offset. Working with either or both of these data representations for nine geologic structures, which range in age from 3 to 400 ka B.P. and in offset from 1 to 50 m, we apply analytical solutions for the vertical initial value scarp, the vertical continuous offset scarp, and the finite slope, initial value scarp. The model calculations are intrinsically ambiguous, yielding as the final answer only the product κ t (in the case of the initial value problem) or the product κ A −1 (in the case of the repeated faulting problem); here t is the age of a single scarp-forming event, 2 A is the vertical slip rate, and κ is the “mass diffusivity.” A single profile across three sea cliffs along the Santa Cruz, California, coast is analyzed as three separate initial value problems. A reasonably constrained age for the sea cliff standing above the Highway 1 platform returns κ = 11 GKG (1 GKG = 1 m 2 /ka). With this κ, we can date the two older sea cliffs. In fact, we do the converse: age estimates for these two older sea cliffs based on a uniform rate of uplift both yield the same κ as for the lower sea cliff. We treat a single profile of the Raymond fault in Pasadena/San Marino in terms of the repeated faulting problem; for it the uplift rate of R. Crook and others yields κ = 16 GKG. The very substantial preexisting offset across the Raymond fault must have been buried/leveled some 230 ka B.P., when the modern topography began to form. Our analysis of the Lake Bonneville shoreline scarps reveals a dependence of κ t on 2a, suggestive of nonlinear modification processes. This appearance is treated with the finite slope initial value scarp model to determine κ=1.1 GKG for the Lake Bonneville shoreline scarps. The suggestion of M. N. Machette that approximately 100,000-year-old, meter-high scarps are “unobservable” in weakly consolidated alluvial terranes of the Basin and Range and Rio Grande Rift Valley provinces can be formulated as κ ≳ 1 GKG. The coincidence between this inequality and the Lake Bonneville shoreline κ is striking, and it suggests that the value of κ = 1 GKG may be generally applicable, as a good first approximation, to the modification of alluvial terranes within the semiarid regions of the western United States. The Lake Bonneville shoreline κ is the basis for dating four sets of fault scarps in west-central Utah. The Drum Mountains fault scarps can be modeled in several different circumstances, but the most likely interpretation is that these fault scarps formed as the result of a single episode of normal faulting 3.6 to 5.7 ka B.P. The younger age is associated with quite low initial slope angles (25°). The other three sets of fault scarps show no evidence for finite initial value slopes. Fault scarps along the eastern base of the Fish Springs Range are very young, 3 ka B.P. We estimate the age of fault scarps along the western flank of the Oquirrh Mountains to be 32 ka B.P., which meets the weak geologic constraint that they be older than the Lake Bonneville shoreline. Fault scarps along the northeastern margin of the Sheeprock Mountains are even older, 53 ka B.P. An intriguing consequence of our single-event analysis of these scarps is that an 11.5-m offset occurred in a single earthquake.

Journal of Geophysical Research Solid Earth↗