USGS Science⌕ Search

SEARCH · USGS Science

Results for “Energy”

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 1,513 records · Page 84Linked to original sources
USGS2000

A compartmentalized solute transport model for redox zones in contaminated aquifers: 1. Theory and development

This paper, the first of two parts [see Abrams and Loague , this issue], takes the compartmentalized approach for the geochemical evolution of redox zones presented by Abrams et al . [1998] and embeds it within a solute transport framework. In this paper the compartmentalized approach is generalized to facilitate the description of its incorporation into a solute transport simulator. An equivalent formulation is developed which removes any discontinuities that may occur when switching compartments. Rate‐limited redox reactions are modeled with a modified Monod relationship that allows either the organic substrate or the electron acceptor to be the rate‐limiting reactant. Thermodynamic constraints are used to inhibit lower‐energy redox reactions from occurring under infeasible geochemical conditions without imposing equilibrium on the lower‐energy reactions. The procedure used allows any redox reaction to be simulated as being kinetically limited or thermodynamically limited, depending on local geochemical conditions. Empirical reaction inhibition methods are not needed. The sequential iteration approach (SIA), a technique which allows the number of solute transport equations to be reduced, is adopted to solve the coupled geochemical/solute transport problem. When the compartmentalized approach is embedded within the SIA, with the total analytical concentration of each component as the dependent variable in the transport equation, it is possible to reduce the number of transport equations even further than with the unmodified SIA. A one‐dimensional, coupled geochemical/solute transport simulation is presented in which redox zones evolve dynamically in time and space. The compartmentalized solute transport (COMPTRAN) model described in this paper enables the development of redox zones to be simulated under both kinetic and thermodynamic constraints. The modular design of COMPTRAN facilitates the use of many different, preexisting solute transport and geochemical codes. The companion paper [ Abrams and Loague , this issue] presents examples of the application of COMPTRAN to field‐scale problems.

Water Resources Research↗
USGS2001

Stress drop with constant, scale independent seismic efficiency and overshoot

To model dissipated and radiated energy during earthquake stress drop, I calculate dynamic fault slip using a single degree of freedom spring-slider block and a laboratory-based static/kinetic fault strength relation with a dynamic stress drop proportional to effective normal stress. The model is scaled to earthquake size assuming a circular rupture; stiffness varies inversely with rupture radius, and rupture duration is proportional to radius. Calculated seismic efficiency, the ratio of radiated to total energy expended during stress drop, is in good agreement with laboratory and field observations. Predicted overshoot, a measure of how much the static stress drop exceeds the dynamic stress drop, is higher than previously published laboratory and seismic observations and fully elasto-dynamic calculations. Seismic efficiency and overshoot are constant, independent of normal stress and scale. Calculated variation of apparent stress with seismic moment resembles the observational constraints of McGarr [1999].

Geophysical Research Letters↗
USGS2001

Moss and lichen cover mapping at local and regional scales in the boreal forest ecosystem of central Canada

Mosses and lichens are important components of boreal landscapes [Vitt et al., 1994; Bubier et al., 1997]. They affect plant productivity and belowground carbon sequestration and alter the surface runoff and energy balance. We report the use of multiresolution satellite data to map moss and lichens over the BOREAS region at a 10 m, 30 m, and 1 km scales. Our moss and lichen classification at the 10 m scale is based on ground observations of associations among soil drainage classes, overstory composition, and cover type among four broad classes of ground cover (feather, sphagnum, and brown mosses and lichens). For our 30 m map, we used field observations of ground cover-overstory associations to map mosses and lichens in the BOREAS southern study area (SSA). To scale up to a 1 km (AVHRR) moss map of the BOREAS region, we used the TM SSA mosaics plus regional field data to identify AVHRR overstory-ground cover associations. We found that: 1) ground cover, overstory composition and density are highly correlated, permitting inference of moss and lichen cover from satellite-based land cover classifications; 2) our 1 km moss map reveals that mosses dominate the boreal landscape of central Canada, thereby a significant factor for water, energy, and carbon modeling; 3) TM and AVHRR moss cover maps are comparable; 4) satellite data resolution is important; particularly in detecting the smaller wetland features, lakes, and upland jack pine sites; and 5) distinct regional patterns of moss and lichen cover correspond to latitudinal and elevational gradients. Copyright 2001 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2004

Crystallographic controls on the frictional behavior of dry and water-saturated sheet structure minerals

We compare the frictional strengths of 17 sheet structure mineral powders, measured under dry and water-saturated conditions, to identify the factors that cause many of them to be relatively weak. The dry coefficient of friction μ ranges upward from 0.2 for graphite, leveling off at 0.8 for margarite, clintonite, gibbsite, kaolinite, and lizardite. The values of μ (dry) correlate directly with calculated (001) interlayer bond strengths of the minerals. This correlation occurs because shear becomes localized along boundary and Riedel shears and the platy minerals in them rotate into alignment with the shear planes. For those gouges with μ (dry) < 0.8, shear occurs by breaking the interlayer bonds to form new cleavage surfaces. Where μ (dry) = 0.8, consistent with Byerlee's law, the interlayer bonds are sufficiently strong that other frictional processes dominate. The transition in dry friction mechanisms corresponds to calculated surface energies of 2–3 J/m 2 . Adding water causes μ to decrease for every mineral tested except graphite. If the minerals are separated into groups with similar crystal structures, μ (wet) increases with increasing interlayer bond strength within each group. This relationship also holds for the swelling clay montmorillonite, whose water-saturated strength is consistent with the strengths of nonswelling clays of similar crystal structure. Water in the saturated gouges forms thin, structured films between the plate surfaces. The polar water molecules are bonded to the plate surfaces in proportion to the mineral's surface energy, and μ (wet) reflects the stresses required to shear through the water films.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

International heat flow commission celebrates 40 years

The outflow of heat from the Earth's interior is, in terms of energy, the most impressive terrestrial phenomenon. Its present rate of about 1021 joules per year is order‐of‐magnitudes greater than the energy dissipation of earthquakes or heat loss from volcanic eruptions. The study of the Earth's internal heat plays an important role in understanding the Earth's origin, internal constitution, and plate tectonics. The thermal structure of the Earth, generally referred to as the geothermics, or also as the (terrestrial) heat flow, has been studied for a long time. However, modern geothermics, one of the fundamental geophysical disciplines, is relatively young. The International Heat Flow Commission (IHFC), under which the academic geothermal research on the international scale is organized, was created only in 1963. The IHFC, operating under the International Association of Seismology and Physics of the Earth's Interior (IASPEI), covers a broad scope of geophysical studies, and links the activities of other associations of the International Union of Geodesy and Geophysics (IUGG).

Eos, Earth and Space Science News↗
USGS2004

Application of near real-time radial semblance to locate the shallow magmatic conduit at Kilauea Volcano, Hawaii

Radial Semblance is applied to broadband seismic network data to provide source locations of Very-Long-Period (VLP) seismic energy in near real time. With an efficient algorithm and adequate network coverage, accurate source locations of VLP energy are derived to quickly locate the shallow magmatic conduit system at Kilauea Volcano, Hawaii. During a restart in magma flow following a brief pause in the current eruption, the shallow magmatic conduit is pressurized, resulting in elastic radiation from various parts of the conduit system. A steeply dipping distribution of VLP hypocenters outlines a region extending from sea level to about 550 m elevation below and just east of the Halemaumau Pit Crater. The distinct hypocenters suggest the shallow plumbing system beneath Halemaumau consists of a complex plexus of sills and dikes. An unconstrained location for a section of the conduit is also observed beneath the region between Kilauea Caldera and Kilauea Iki Crater.

Hawaii↗
USGS2005

Nearshore wave-induced cyclical flexing of sea cliffs

[1] Evolution of a tectonically active coast is driven by geomorphically destructive energy supplied by ocean waves. Wave energy is episodic and concentrated; sea cliffs are battered by the geomorphic wrecking ball every 4-25 s. We measure the response of sea cliffs to wave assault by sensing the ground motion using near-coastal seismometers. Sea cliffs respond to waves in two distinct styles. High-frequency motion (20 Hz) reflects the natural frequency of the sea cliff as it rings in response to direct wave impact. Low-frequency motion in the 0.1-0.05 Hz (10-20 s) band consistently agrees with the dominant nearshore wave period. Integrating microseismic velocities suggests 50 ??m and 10 ??m displacements in horizontal and vertical directions, respectively. Displacement ellipsoids exhibit simultaneous downward and seaward sea cliff motion with each wave. Video footage corroborates the downward sea cliff flex in response to the imposed water load on the wave cut platform. Gradients in displacement amplitudes documented using multiple seismometers suggest longitudinal and shear strain of the flexing sea cliff on the order of 0.5-4 ?? strains during each wave loading cycle. As this sea cliff flexure occurs approximately 3 million times annually, it has the potential to fatigue the rock through cyclical loading. Local sea cliff retreat rates of 10 cm/yr imply that a given parcel of rock is flexed through roughly 109 cycles of increasing amplitude before exposure to direct wave attack at the cliff face. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

Surface drifter derived circulation in the northern and middle Adriatic Sea: Response to wind regime and season

More than 120 satellite-tracked drifters were deployed in the northern and middle Adriatic (NMA) Sea between September 2002 and November 2003, with the purpose of studying the surface circulation at mesoscale to seasonal scale in relation to wind forcing, river runoff, and bottom topography. Pseudo-Eulerian and Lagrangian statistics were calculated from the low-pass-filtered drifter velocity data between September 2002 and December 2003. The structure of the mean circulation is determined with unprecedented high horizontal resolution by the new data. In particular, mean currents, velocity variance, and kinetic energy levels are shown to be maximal in the Western Adriatic Current (WAC). Separating data into seasons, we found that the mean kinetic energy is maximal in fall, with high values also in winter, while it is significantly weaker in summer. High-resolution Local Area Model Italy winds were used to relate the drifter velocities to the wind fields. The surface currents appear to be significantly influenced by the winds. The mean flow during the northeasterly bora regime shows an intensification of the across-basin recirculating currents. In addition, the WAC is strongly intensified both in intensity and in its offshore lateral extension. In the southeasterly sirocco regime, northward flow without recirculation dominates in the eastern half of the basin, while during northwesterly maestro the WAC is enhanced. Separating the data into low and high Po River discharge rates for low-wind conditions shows that the WAC and the velocity fluctuations in front of the Po delta are stronger for high Po River runoff. Lagrangian covariance, diffusivity, and integral time and space scales are larger in the along-basin direction and are maximal in the southern portion of the WAC. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2006

Mountain hydrology of the western United States

Climate change and climate variability, population growth, and land use change drive the need for new hydrologic knowledge and understanding. In the mountainous West and other similar areas worldwide, three pressing hydrologic needs stand out: first, to better understand the processes controlling the partitioning of energy and water fluxes within and out from these systems; second, to better understand feedbacks between hydrological fluxes and biogeochemical and ecological processes; and, third, to enhance our physical and empirical understanding with integrated measurement strategies and information systems. We envision an integrative approach to monitoring, modeling, and sensing the mountain environment that will improve understanding and prediction of hydrologic fluxes and processes. Here extensive monitoring of energy fluxes and hydrologic states are needed to supplement existing measurements, which are largely limited to streamflow and snow water equivalent. Ground‐based observing systems must be explicitly designed for integration with remotely sensed data and for scaling up to basins and whole ranges.

Water Resources Research↗
USGS2007

Integrating laboratory creep compaction data with numerical fault models: A Bayesian framework

[1] We developed a robust Bayesian inversion scheme to plan and analyze laboratory creep compaction experiments. We chose a simple creep law that features the main parameters of interest when trying to identify rate-controlling mechanisms from experimental data. By integrating the chosen creep law or an approximation thereof, one can use all the data, either simultaneously or in overlapping subsets, thus making more complete use of the experiment data and propagating statistical variations in the data through to the final rate constants. Despite the nonlinearity of the problem, with this technique one can retrieve accurate estimates of both the stress exponent and the activation energy, even when the porosity time series data are noisy. Whereas adding observation points and/or experiments reduces the uncertainty on all parameters, enlarging the range of temperature or effective stress significantly reduces the covariance between stress exponent and activation energy. We apply this methodology to hydrothermal creep compaction data on quartz to obtain a quantitative, semiempirical law for fault zone compaction in the interseismic period. Incorporating this law into a simple direct rupture model, we find marginal distributions of the time to failure that are robust with respect to errors in the initial fault zone porosity.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

The source of infrasound associated with long-period events at mount St. Helens

During the early stages of the 2004-2008 Mount St. Helens eruption, the source process that produced a sustained sequence of repetitive long-period (LP) seismic events also produced impulsive broadband infrasonic signals in the atmosphere. To assess whether the signals could be generated simply by seismic-acoustic coupling from the shallow LP events, we perform finite difference simulation of the seismo-acoustic wavefield using a single numerical scheme for the elastic ground and atmosphere. The effects of topography, velocity structure, wind, and source configuration are considered. The simulations show that a shallow source buried in a homogeneous elastic solid produces a complex wave train in the atmosphere consisting of P/SV and Rayleigh wave energy converted locally along the propagation path, and acoustic energy originating from , the source epicenter. Although the horizontal acoustic velocity of the latter is consistent with our data, the modeled amplitude ratios of pressure to vertical seismic velocity are too low in comparison with observations, and the characteristic differences in seismic and acoustic waveforms and spectra cannot be reproduced from a common point source. The observations therefore require a more complex source process in which the infrasonic signals are a record of only the broadband pressure excitation mechanism of the seismic LP events. The observations and numerical results can be explained by a model involving the repeated rapid pressure loss from a hydrothermal crack by venting into a shallow layer of loosely consolidated, highly permeable material. Heating by magmatic activity causes pressure to rise, periodically reaching the pressure threshold for rupture of the "valve" sealing the crack. Sudden opening of the valve generates the broadband infrasonic signal and simultaneously triggers the collapse of the crack, initiating resonance of the remaining fluid. Subtle waveform and amplitude variability of the infrasonic signals as recorded at an array 13.4 km to the NW of the volcano are attributed primarily to atmospheric boundary layer propagation effects, superimposed upon amplitude changes at the source.

Washington↗
USGS2010

Inelastic off-fault response and three-dimensional dynamics of earthquake rupture on a strike-slip fault

Large dynamic stress off the fault incurs an inelastic response and energy loss, which contributes to the fracture energy, limiting the rupture and slip velocity. Using an explicit finite element method, we model three-dimensional dynamic ruptures on a vertical strike-slip fault in a homogeneous half-space. The material is subjected to a pressure-dependent Drucker-Prager yield criterion. Initial stresses in the medium increase linearly with depth. Our simulations show that the inelastic response is confined narrowly to the fault at depth. There the inelastic strain is induced by large dynamic stresses associated with the rupture front that overcome the effect of the high confining pressure. The inelastic zone increases in size as it nears the surface. For material with low cohesion (~5 MPa) the inelastic zone broadens dramatically near the surface, forming a "flowerlike" structure. The near-surface inelastic strain occurs in both the extensional and the compressional regimes of the fault, induced by seismic waves ahead of the rupture front under a low confining pressure. When cohesion is large (~10 MPa), the inelastic strain is significantly reduced near the surface and confined mostly to depth. Cohesion, however, affects the inelastic zone at depth less significantly. The induced shear microcracks show diverse orientations near the surface, owing to the low confining pressure, but exhibit mostly horizontal slip at depth. The inferred rupture-induced anisotropy at depth has the fast wave direction along the direction of the maximum compressive stress.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Introduction to special section on impacts of land use change on water resources

Changes in land use have potentially large impacts on water resources, yet quantifying these impacts remains among the more challenging problems in hydrology. Water, food, energy, and climate are linked through complex webs of direct and indirect effects and feedbacks. Land use is undergoing major changes due not only to pressures for more efficient food, feed, and fiber production to support growing populations but also due to policy shifts that are creating markets for biofuel and agricultural carbon sequestration. Hydrologic systems embody flows of water, solutes, sediments, and energy that vary even in the absence of human activity. Understanding land use impacts thus necessitates integrated scientific approaches. Field measurements, remote sensing, and modeling studies are shedding new light on the modes and mechanisms by which land use changes impact water resources. Such studies can help deconflate the interconnected influences of human actions and natural variations on the quantity and quality of soil water, surface water, and groundwater, past, present, and future.

Water Resources Research↗
USGS2010

Analysis of nonvolcanic tremor on the San Andreas Fault near Parkfield, CA using U.S. Geological Survey Parkfield Seismic Array

Reports by Nadeau and Dolenc (2005) that tremor had been detected near Cholame Valley spawned an effort to use UPSAR (U. S. Geological Survey Parkfield Seismic Array) to study characteristics of tremor. UPSAR was modified to record three channels of velocity at 40–50 sps continuously in January 2005 and ran for about 1 month, during which time we recorded numerous episodes of tremor. One tremor, on 21 January at 0728, was recorded with particularly high signal levels as well as another episode 3 days later. Both events were very emergent, had a frequency content between 2 and 8 Hz, and had numerous high-amplitude, short-duration arrivals within the tremor signal. Here using the first episode as an example, we discuss an analysis procedure, which yields azimuth and apparent velocity of the tremor at UPSAR. We then provide locations for both tremor episodes. The emphasis here is how the tremor episode evolves. Twelve stations were operating at the time of recording. Slowness of arrivals was determined using cross correlation of pairs of stations; the same method used in analyzing the main shock data from 28 September 2004. A feature of this analysis is that 20 s of the time series were used at a time to calculate correlation; the longer windows resulted in more consistent estimates of slowness, but lower peak correlations. These values of correlation (peaks of about 0.25), however, are similar to that obtained for the S wave of a microearthquake. Observed peaks in slowness were traced back to source locations assumed to lie on the San Andreas fault. Our inferred locations for the two tremor events cluster near the locations of previously observed tremor, south of the Cholame Valley. Tremor source depths are in the 14–24 km range, which is below the seismogenic brittle zone, but above the Moho. Estimates of error do not preclude locations below the Moho, however. The tremor signal is very emergent but contains packets that are several times larger than the background tremor signal and lasts about 5 s. These impulsive wavelets are similar to low-frequency earthquakes signals seen in Japan but appear to be broader band rather than just higher in low-frequency energy. They may be more appropriately called high-energy tremor (HET). HET signals at UPSAR correlate well with the record of this event from station GHIB of the HRSN borehole array at Parkfield and HETs typically have a higher cross-correlation coefficient than the rest of the tremor event. The amplitudes of a large HET are consistent with a magnitude of 0.1 when compared with a M2.3 event that had about the same epicenter. Polarizations of the tremor episode at UPSAR are mostly just north of east. Both linearity and azimuth evolve over time suggesting a change in tremor source location over time and linearity is typically higher at the HETs.

California↗
USGS2011

Volcanic plume height measured by seismic waves based on a mechanical model

In August 2008 an unmonitored, largely unstudied Aleutian volcano, Kasatochi, erupted catastrophically. Here we use seismic data to infer the height of large eruptive columns such as those of Kasatochi based on a combination of existing fluid and solid mechanical models. In so doing, we propose a connection between a common, observable, short-period seismic wave amplitude to the physics of an eruptive column. To construct a combined model, we estimate the mass ejection rate of material from the vent on the basis of the plume height, assuming that the height is controlled by thermal buoyancy for a continuous plume. Using the estimated mass ejection rate, we then derive the equivalent vertical force on the Earth through a momentum balance. Finally, we calculate the far-field surface waves resulting from the vertical force. The model performs well for recent eruptions of Kasatochi and Augustine volcanoes if v , the velocity of material exiting the vent, is 120-230 m s -1 . The consistency between the seismically inferred and measured plume heights indicates that in these cases the far-field ~1 s seismic energy radiated by fluctuating flow in the volcanic jet during the eruption is a useful indicator of overall mass ejection rates. Thus, use of the model holds promise for characterizing eruptions and evaluating ash hazards to aircraft in real time on the basis of far-field short-period seismic data. This study emphasizes the need for better measurements of eruptive plume heights and a more detailed understanding of the full spectrum of seismic energy radiated coeruptively.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Hydrodynamics of a bathymetrically complex fringing coral reef embayment: Wave climate, in situ observations, and wave prediction

This paper examines the relationship between offshore wave climate and nearshore waves and currents at Hanalei Bay, Hawaii, an exposed bay fringed with coral reefs. Analysis of both offshore in situ data and numerical hindcasts identify the predominance of two wave conditions: a mode associated with local trade winds and an episodic pattern associated with distant source long-period swells. Analysis of 10 months of in situ data within the bay show that current velocities are up to an order of magnitude greater during long-period swell episodes than during trade wind conditions; overall circulation patterns are also fundamentally different. The current velocities are highly correlated with incident wave heights during the swell episodes, while they are not during the modal trade wind conditions. A phase-averaged wave model was implemented with the dual purpose of evaluating application to bathymetrically complex fringing reefs and to examine the propagation of waves into the nearshore in an effort to better explain the large difference in observed circulation during the two offshore wave conditions. The prediction quality of this model was poorer for the episodic condition than for the lower-energy mode, however, it illustrated how longer-period swells are preferentially refracted into the bay and make available far more nearshore wave energy to drive currents compared to waves during modal conditions. The highly episodic circulation, the nature of which is dependent on complex refraction patterns of episodic, long-period swell has implications for flushing and sediment dynamics for incised fringing reef-lined bays that characterize many high islands at low latitudes around the world.

Hawai'i↗
USGS2011

The influence of the Atlantic Warm Pool on the Florida panhandle sea breeze

In this paper we examine the variations of the boreal summer season sea breeze circulation along the Florida panhandle coast from relatively high resolution (10 km) regional climate model integrations. The 23 year climatology (1979–2001) of the multidecadal dynamically downscaled simulations forced by the National Centers for Environmental Prediction–Department of Energy (NCEP‐DOE) Reanalysis II at the lateral boundaries verify quite well with the observed climatology. The variations at diurnal and interannual time scales are also well simulated with respect to the observations. We show from composite analyses made from these downscaled simulations that sea breezes in northwestern Florida are associated with changes in the size of the Atlantic Warm Pool (AWP) on interannual time scales. In large AWP years when the North Atlantic Subtropical High becomes weaker and moves further eastward relative to the small AWP years, a large part of the southeast U.S. including Florida comes under the influence of relatively strong anomalous low‐level northerly flow and large‐scale subsidence consistent with the theory of the Sverdrup balance. This tends to suppress the diurnal convection over the Florida panhandle coast in large AWP years. This study is also an illustration of the benefit of dynamic downscaling in understanding the low‐frequency variations of the sea breeze.

Florida↗
USGS2012

Relating stick-slip friction experiments to earthquake source parameters

Analytical results for parameters, such as static stress drop, for stick-slip friction experiments, with arbitrary input parameters, can be determined by solving an energy-balance equation. These results can then be related to a given earthquake based on its seismic moment and the maximum slip within its rupture zone, assuming that the rupture process entails the same physics as stick-slip friction. This analysis yields overshoots and ratios of apparent stress to static stress drop of about 0.25. The inferred earthquake source parameters static stress drop, apparent stress, slip rate, and radiated energy are robust inasmuch as they are largely independent of the experimental parameters used in their estimation. Instead, these earthquake parameters depend on C, the ratio of maximum slip to the cube root of the seismic moment. C is controlled by the normal stress applied to the rupture plane and the difference between the static and dynamic coefficients of friction. Estimating yield stress and seismic efficiency using the same procedure is only possible when the actual static and dynamic coefficients of friction are known within the earthquake rupture zone.

Geophysical Research Letters↗