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,531 records · Page 85Linked to original sources
USGS2012

Calving seismicity from iceberg-sea surface interactions

Iceberg calving is known to release substantial seismic energy, but little is known about the specific mechanisms that produce calving icequakes. At Yahtse Glacier, a tidewater glacier on the Gulf of Alaska, we draw upon a local network of seismometers and focus on 80 hours of concurrent, direct observation of the terminus to show that calving is the dominant source of seismicity. To elucidate seismogenic mechanisms, we synchronized video and seismograms to reveal that the majority of seismic energy is produced during iceberg interactions with the sea surface. Icequake peak amplitudes coincide with the emergence of high velocity jets of water and ice from the fjord after the complete submergence of falling icebergs below sea level. These icequakes have dominant frequencies between 1 and 3 Hz. Detachment of an iceberg from the terminus produces comparatively weak seismic waves at frequencies between 5 and 20 Hz. Our observations allow us to suggest that the most powerful sources of calving icequakes at Yahtse Glacier include iceberg-sea surface impact, deceleration under the influence of drag and buoyancy, and cavitation. Numerical simulations of seismogenesis during iceberg-sea surface interactions support our observational evidence. Our new understanding of iceberg-sea surface interactions allows us to reattribute the sources of calving seismicity identified in earlier studies and offer guidance for the future use of seismology in monitoring iceberg calving.

Journal of Geophysical Research F: Earth Surface↗
USGS2018

Crustal structure and quaternary acceleration of deformation rates in central Washington revealed by stream profile inversion, potential field geophysics, and structural geology of the Yakima folds

Post‐Miocene tectonic uplift along fault‐cored anticlines within central Washington produced the Yakima Fold Province, a region of active NNE‐SSW shortening in the Cascadian backarc. The relative timing and rate of deformation along individual structures is coarsely defined yet imperative for seismic hazard assessment. In this work, we use geomorphic and geophysical mapping, stream profile inversion, and balanced cross‐section methods to constrain fault geometries and slip rates in the Yakima Canyon region. We extract stream profiles from LiDAR data and analytically solve for the rate of relative rock uplift along several active fault‐cored anticlines. To constrain the fault geometries at depth and the long‐term magnitude of deformation, we constructed two line‐balanced cross sections across the folds with forward‐modeled magnetic and gravity anomaly data. Our stream profile results indicate an increase of incision rates in the Pleistocene, and we infer the increase is tectonically controlled. We estimate modern slip rates between 0.4 and 0.6 mm/year accommodated on reverse faults that core the Manastash Ridge, Umtanum Ridge, and Selah Butte anticlines and establish that these faults reactivate and invert older normal faults in basement rocks. Finally, we calculate the time required to accumulate sufficient strain energy for a large magnitude earthquake (M ≥ 7) along individual structures in the Yakima Fold Province. Results show that the Yakima folds likely accommodate large magnitude earthquakes and that it takes several hundred to several thousand years to accumulate sufficient strain energy for an M ≥ 7 earthquake.

Washington↗
USGS2018

Climate change and future wildfire in the western USA: An ecological approach to nonstationarity

We developed ecologically based climate‐fire projections for the western United States. Using a finer ecological classification and fire‐relevant climate predictors, we created statistical models linking climate and wildfire area burned for ecosections, which are geographic delineations based on biophysical variables. The results indicate a gradient from purely fuel‐limited (antecedent positive water balance anomalies or negative energy balance anomalies) to purely flammability‐limited (negative water balance anomalies or positive energy balance anomalies) fire regimes across ecosections. Although there are other influences (such as human ignitions and management) on fire occurrence and area burned, seasonal climate significantly explains interannual fire area burned. Differences in the role of climate across ecosections are not random, and the relative dominance of climate predictors allows objective classification of ecosection climate‐fire relationships. Expected future trends in area burned range from massive increases, primarily in flammability limited systems near the middle of the water balance deficit distribution, to substantial decreases, in fuel‐limited nonforested systems. We predict increasing area burned in most flammability‐limited systems but predict decreasing area burned in primarily fuel‐limited systems with a flammability‐limited (“hybrid”) component. Compared to 2030–2059 (2040s), projected area burned for 2070–2099 (2080s) increases much more in the flammability and flammability‐dominated hybrid systems than those with equal control and continues to decrease in fuel‐limited hybrid systems. Exceedance probabilities for historical 95th percentile fire years are larger in exclusively flammability‐limited ecosections than in those with fuel controls. Filtering the projected results using a fire‐rotation constraint minimizes overprojection due to static vegetation assumptions, making projections more conservative.

Earth's Future↗
USGS2018

Hydrodynamics and sediment mobility processes over a degraded senile coral reef

Coral reefs can influence hydrodynamics and morphodynamics by dissipating and refracting incident wave energy, modifying circulation patterns, and altering sediment transport pathways. In this study, the sediment and hydrodynamic response of a senile (dead) barrier reef (Crocker Reef, located in the upper portion of the Florida Reef Tract) to storms and quiescent conditions was evaluated using field observations and the Coupled Ocean-Atmosphere-Wave-Sediment Transport model. Waves, circulation, and resultant sediment mobility were modeled across different reef zones. Sediment mobility during quiescent periods and the passage of far-field storms are driven by nonbreaking waves and, to a lesser degree, regional circulation. Spatial variability in these processes produces the present-day distribution of sediment grain size at Crocker Reef, wherein finer-grain material along a shallow central ridge is frequently mobilized (43% to 62% of the time), winnowed away, and deposited along the lower-energy flanks and in the fore reef where sand mobility occurs less frequently (32% to 43% and 1% to 22% of the time, respectively). Analysis of wave conditions for the period of 2006–2014 supports that wave heights rarely exceed the threshold for breaking (0.1% and 0.3% at the reef crest and at the reef flat, respectively), predominantly during the passage of tropical storms. There is a shift to a wave-breaking regime during near-field storms, creating the potential for mobilization of larger material and enhanced reef degradation. Sediment mobility can be enhanced due to wave skewness or the generation of free infragravity waves during periods of depth-induced wave breaking.

Florida↗
USGS2019

Induced seismicity reduces seismic hazard?

Earthquakes caused by human activities have been observed for decades. Often these are related to industrial activities pumping fluids into deep geologic formations, like with wastewater disposal. The simplest theory connecting these processes to earthquakes is straightforward: injection leads to fluid pressure changes that either reduce the strength of preexisting faults or generate new faults. In practice, the conditions that lead to induced earthquakes are not always clear in ways that can be generalized. Kao et al. (2018, https://doi.org/10.1029/2018GL079288) show how the distribution of induced earthquakes in Western Canada relate to natural rates of deformation in the crust. Using these new results, they discuss an intriguing paradox: induced seismicity can cause short-term increases in the seismic hazard that are followed by a period of reduced seismic hazard. Such hazard-reducing scenarios are plausible but hinge upon simplifying assumptions about how the crust stores and releases strain energy in the form of earthquakes.

Geophysical Research Letters↗
USGS2021

Permafrost promotes shallow groundwater flow and warmer headwater streams

The presence of permafrost influences the flow paths of water through Arctic landscapes and thereby has the potential to impact stream discharge and thermal regimes. Observations from eleven headwater streams in Alaska showed that July water temperatures were higher in catchments with more near‐surface permafrost. We apply a fully coupled cryohydrology model to investigate if the impact of permafrost on flow path depth could cause the same pattern in temperatures of groundwater discharging from hillslopes to streams. The model simulates surface energy and water balances, snow, and subsurface water and energy balances for two‐dimensional hillslope model cases with varying permafrost extent. We find that hillslopes with continuous permafrost have more shallow flow paths and twice as high rates of evapotranspiration, compared to hillslopes with no permafrost. For our simulated cases, 6.7 % of the horizontal water flux moves through the top organic soil layers when there is continuous permafrost, while only 0.5 % moves through organic layers without permafrost. The deeper flow paths in permafrost‐free simulations buffer seasonal temperature extremes, so that summer groundwater discharge temperatures are highest with continuous permafrost. Our results suggest that permafrost thawing alters groundwater flow paths and can lead to decreases in summer stream temperatures and reductions in evapotranspiration in headwater catchments. These changes are of potential importance for stream biotic components of ecosystems, however, the full impact remains unknown.

Alaska↗
USGS2021

Urbanization impacts on evapotranspiration across various spatio-temporal scales

Urbanization has been shown to locally increase the nighttime temperatures creating urban heat islands, which partly arise due to evapotranspiration (ET) reduction. It is unclear how the direction and magnitude of the change in local ET due to urbanization varies globally across different climatic regimes. This knowledge gap is critical, both for the key role of ET in the energy and water balance accounting for the majority of local precipitation, and for reducing the urban heat island effect. We explore and assess the impacts of urbanization on monthly and mean annual ET across a range of landscapes from local to global spatial scales. Remotely sensed land cover and ET available at 1 km resolution are used to quantify the differences in ET between urban and surrounding non-urban areas across the globe. The observed patterns show that the statistically significant difference between urban and non-urban ET can be estimated to first order as a function of local hydroclimate, with arid regions seeing increased ET, and humid regions showing decreased ET. Cities under cold climates also evaporate more than their non-urban surroundings during the winter, as the urban micro-climate has increased energy availability resulting from human activities. Increased ET in arid cities arises from municipal water withdrawals and increased irrigation during drought conditions. These results can help inform planners to improve the integration of environmental conditions into the design and management of urban landscapes.

Earth's Future↗
USGS2022

Multidisciplinary constraints on the thermal-chemical boundary between Earth's core and mantle

Heat flux from the core to the mantle provides driving energy for mantle convection thus powering plate tectonics, and contributes a significant fraction of the geothermal heat budget. Indirect estimates of core-mantle boundary heat flow are typically based on petrological evidence of mantle temperature, interpretations of temperatures indicated by seismic travel times, experimental measurements of mineral melting points, physical mantle convection models, or physical core convection models. However, previous estimates have not consistently integrated these lines of evidence. In this work, an interdisciplinary analysis is applied to co-constrain core-mantle boundary heat flow and test the thermal boundary layer (TBL) theory. The concurrence of TBL models, energy balance to support geomagnetism, seismology, and review of petrologic evidence for historic mantle temperatures supports Q CMB ∼15 TW, with all except geomagnetism supporting as high as ∼20 TW. These values provide a tighter constraint on core heat flux relative to previous work. Our work describes the seismic properties consistent with a TBL, and supports a long-lived basal mantle molten layer through much of Earth's history.

Geochemistry, Geophysics, Geosystems↗
USGS2021

Electrical properties of carbon dioxide hydrate: Implications for monitoring CO2 in the gas hydrate stability zone

CO 2 and CH 4 clathrate hydrates are of keen interest for energy and carbon cycle considerations. While both typically form on Earth as cubic structure I (sI), we find that pure CO 2 hydrate exhibits over an order of magnitude higher electrical conductivity ( σ ) than pure CH 4 hydrate at geologically relevant temperatures. The conductivity was obtained from frequency-dependent impedance ( Z ) measurements made on polycrystalline CO 2 hydrate (CO 2 ·6.0 ± 0.2H 2 O by methods here) with 25% gas-filled porosity, compared with CH 4 hydrate (CH 4 ·5.9H 2 O) formed and measured in the same apparatus and exhibiting closely matching grain characteristics. The conductivity of CO 2 hydrate is 6.5 × 10 −4 S/m at 273K with an activation energy ( E a ) of 46.5 kJ/mol at 260–281 K, compared with ∼5 × 10 −5 S/m and 34.8 kJ/m for CH 4 hydrate. Equivalent circuit modeling indicates that different pathways govern conduction in CO 2 versus CH 4 hydrate. Results show promise for use of electromagnetic methods in monitoring CO 2 hydrate formation in certain natural settings or in CO 2 /CH 4 exchange efforts.

Geophysical Research Letters↗
USGS2024

Laboratory hydrofractures as analogs to tectonic tremors

The fracture of Earth materials occurs over a wide range of time and length scales. Physical conditions, particularly the stress field and Earth material properties, may condition rupture in a specific fracture regime. In nature, fast and slow fractures occur concurrently: tectonic tremor events are fast enough to emit seismic waves and frequently accompany slow earthquakes, which are too slow to emit seismic waves and are referred to as aseismic slip events. In this study, we generate simultaneous seismic and aseismic processes in a laboratory setting by driving a penny-shaped crack in a transparent sample with pressurized fluid. We leverage synchronized high-speed imaging and high-frequency acoustic emission (AE) sensing to visualize and listen to the various sequences of propagation (breaks) and arrest (sticks) of a fracture undergoing stick-break instabilities. Slow radial crack propagation is facilitated by fast tangential fractures. Fluid viscosity and pressure regulate the fracture dynamics of slow and fast events, and control the inter-event time and the energy released during individual fast events. These AE signals share behaviors with observations of episodic tremors in Cascadia, United States; these include: (a) bursty or intermittent slow propagation, and (b) nearly linear scaling of radiated energy with area. Our laboratory experiments provide a plausible model of tectonic tremor as an indicative of hydraulic fracturing facilitating shear slip during slow earthquakes.

AGU Advances↗
USGS2026

Global primary and secondary microseism multi-decade geographic variation, secular intensification, and period lengthening

Earth's long-period background seismic wavefield is dominated by two distinct processes that couple ocean wave energy to a global microseism wavefield. We assess global microseism intensity in the secondary (4–10 s) and primary (14–20 s) bands, and across eight 2 s-wide period bands between 4 and 20 s. Robustly estimated primary and secondary secular amplitude trends are estimated at 73 globally distributed seismic station sites with continuous recording spanning at least 20 years, from as early as the late 1980s through October 2025. These trends are positive at 3⁢𝜎 significance for 61 (84%) and 46 (63%) stations with global average rates for vertical-component acceleration of 0.17 ± 0.04 and 0.11 ± 0.05%/yr, for the primary and secondary bands, respectively, with corresponding rates of energy increase of 0.27 ± 0.08 and 0.15 ± 0.09%/yr. Secular intensification is also observed within all 2 s period bands between 4 and 20 s. Amplitude histories for the longest primary microseism periods (18–20 s) correlate to near-antipodal distances, reflecting long-range teleconnections attributed to large-fetch storm systems, long-range swell and Rayleigh wave propagation, and geographically correlated El Niño Southern Oscillation and other geographically extensive atmospheric influences on storms and waves. The lower average rates of intensification for the secondary microseism suggest that crossing wave systems in remote regions are either under-observed or are intensifying more slowly than the primary microseism, possibly due to increasing swell unidirectionality. Secular intensification is greatest at the longest primary microseism periods. This is consistent with a broadening of the global ocean wave spectrum by approximately 0.01%/yr which may reflect an increasing occurrence of large storm systems.

JGR Solid Earth↗
USGS2026

Erosional and hydrological controls on the age and thermochemical stability of particulate organic carbon in an Arctic river

Understanding the mechanisms that drive the mobilization and fate of organic carbon (OC) in Arctic landscapes is important for modeling the feedbacks among permafrost thaw, carbon cycling, and climate change. While significant progress has been made toward measuring in situ OC decomposition in permafrost soils and bulk particulate organic carbon (POC) export from Arctic rivers, few studies have distinguished the source and lability of POC across Arctic landscapes, limiting our ability to predict whether mobilized POC will be oxidized to CO 2 and CH 4 or buried in downstream depositional environments. This study uses ramped pyrolysis/oxidation radiocarbon (RPO- 14 C) analyses to investigate spatial and temporal variations in the thermochemical stability and radiocarbon content of fluvial POC during downstream transport from mountains to the coast in the Canning River (North Slope, Alaska). Fluvial POC in the headwaters is predominantly comprised of high activation energy, thermally recalcitrant petrogenic OC (OC petro ) derived from shale bedrock. Moving into the foothills and low-relief coastal plains, river bank erosion primarily drives mobilization of labile, low activation energy, soil-derived OC (OC soil ). Fluvial POC in mountainous upstream reaches consisted of ∼70% OC petro and just ∼30% OC soil , while POC in the downstream coastal plain reaches comprised ∼85% OC soil and ∼15% OC petro . The high relative lability of POC exported to the coast indicates high susceptibility to oxidation and microbial decomposition, which could enhance CO 2 release as the Arctic hydrologic cycle intensifies. However, the persistence of refractory OC petro in the suspended load indicates the potential for long-term burial of rock organic carbon in marine sediments.

Alaska↗
USGS2006

National workshop on gas hydrates

The range of present knowledge on the subject of gas hydrates and related federal research programs was the topic of discussion at the National Workshop on Gas Hydrates, April 23–24. The intention of the meeting was to provide the impetus for an expanded and broader‐based national research program in both academia and government. Held at the U.S. Geological Survey National Center, Reston, Va., the workshop was organized by Michael D. Max, Naval Research Laboratory, Washington, D.C.; William P. Dillon, USGS, Woods Hole, Mass.; and Rodney D. Malone, U.S. Department of Energy, Morgantown Energy Technology Center, Morgantown, W.Va. The 33 attendees represented academia (33%), federal agencies (58%), and industry (9%).

Eos, Transactions, American Geophysical Union↗
USGS2012

An episode of reinflation of the Long Valley Caldera, eastern California: 1989-1991

Following the episodes of inflation of the resurgent dome associated with the May 1980 earthquake sequence (four M 6 earthquakes) and the January 1983 earthquake swarm (two M 5.2 events), 7 years of frequently repeated two-color geodimeter measurements spanning the Long Valley caldera document gradually decreasing extensional strain rates from 5 ppm/yr in mid-1983, when the measurements began, to near zero in mid-1989. The corresponding seismic activity within the caldera persisted at a low rate of fewer than 10 M ≥ 1.2 earthquakes per week from 1985 through November 1989 with no events exceeding M 3.0. Early October 1989 marked a change in activity when measurements of the two-color geodimeter network showed a significant increase in extensional strain rate (9 ppm/yr) across the caldera. The seismic activity began exceeding 10 M ≥ 1.2 per week in early December 1989 and rapidly increased to a sustained level of tens of M ≥ 1.2 per week with bursts having hundreds of events per day. Many events exceeded M 3.0 and the largest event was M ≈ 4. The 1989–1991 inflation episode is the first time that we have sufficient geodetic measurements in Long Valley to define the temporal relation between onset of an inflation episode and onset of brittle failure (earthquake swarm within the caldera). Here, the onset of deformation preceded the onset of increased earthquake activity by more than 2 months. The seismicity rate began to decrease in mid-July 1990, consistent with a gradually slowing of extension across the caldera as measured by the two-color geodimeter. The recent episode of inflation can be modeled by a single Mogi point source located about 7 km beneath the center of the resurgent dome. In contrast, the deformation pattern observed between mid-1983 and mid-1989 is best reproduced by fault slip in the south moat, inflation at 6.5 km depth near Casa Diablo Hot Springs and inflation beneath the resurgent dome. It appears that the 7-km source beneath the resurgent dome that was active for the earlier episodes is the primary source for the more recent episode. The model used to satisfy the line length observations predicts 7.5 cm of uplift along leveling route along highway 395 from mid-1983 to mid-1989 and an additional 11 cm through the end of 1991. To comp are with the energy release from seismicity, the modeled inflation from late 1989 through the end of 1991 has a moment that is a factor of 40 more than the cumulative seismic moment from earthquakes located within the caldera from the same period. Thus the recent inflation episode represents a significant portion of the observed geodetic deformation with only little seismic energy release.

Journal of Geophysical Research↗
USGS2010

Interpretation of surface flux measurements in heterogeneous terrain during the Monsoon '90 experiment

A network of 9-m-tall surface flux measurement stations were deployed at eight sparsely vegetated sites during the Monsoon '90 experiment to measure net radiation, Q , soil heat flux, G , sensible heat flux, H (using eddy correlation), and latent heat flux, λ E (using the energy balance equation). At four of these sites, 2-m-tall eddy correlation systems were used to measure all four fluxes directly. Also a 2-m-tall Bowen ratio system was deployed at one site. Magnitudes of the energy balance closure ( Q + G + H + λ E ) increased as the complexity of terrain increased. The daytime Bowen ratio decreased from about 10 before the monsoon season to about 0.3 during the monsoons. Source areas of the measurements are developed and compared to scales of heterogeneity arising from the sparse vegetation and the topography. There was very good agreement among simultaneous measurements of Q with the same model sensor at different heights (representing different source areas), but poor agreement among different brands of sensors. Comparisons of simultaneous measurements of G suggest that because of the extremely small source area, extreme care in sensor deployment is necessary for accurate measurement in sparse canopies. A recently published model to estimate fetch is used to interpret measurements of H at the 2 m and 9 m heights. Three sites were characterized by undulating topography, with ridgetops separated by about 200–600 m. At these sites, sensors were located on ridgetops, and the 9-m fetch included the adjacent valley, whereas the 2-m fetch was limited to the immediate ridgetop and hillside. Before the monsoons began, vegetation was mostly dormant, the watershed was uniformly hot and dry, and the two measurements of H were in close agreement. After the monsoons began and vegetation fully matured, the 2-m measurements of H were significantly greater than the 9-m measurements, presumably because the vegetation in the valleys was denser and cooler than on the ridgetops and hillsides. At one lowland site with little topographic relief, the vegetation was more uniform, and the two measurements of H were in close agreement during peak vegetation. Values of λ E could only be compared at two sites, but the 9-m values were greater than the 2-m values, suggesting λ E from the dense vegetation in the valleys was greater than elsewhere.

Water Environment Research↗
USGS1996

Modeling of fluidized ejecta emplacement over digital topography on Venus

The FLOW computer model of McEwen and Malin (1989) modified for application to the study of Venus fluidized ejecta blankets (FEBs) demonstrates that relatively low viscosities, yield strengths, and initial velocities are required to duplicate the observed flow paths of the outflow materials. The model calculates the velocities and simulated flow paths of gravity flows over Magellan topography. The model is formulated to determine flow movements from initial conditions, gravitational acceleration, and resistance to motion as described by Coulomb, viscous, and turbulent resistance forces. Successful duplication of observed FEB flow paths has been achieved for the FEB craters Addams, Isabella, and Cochran. When used as a simple energy‐line model, the model requires low coefficients of friction to extend FEBs to near their observed termini in the synthetic aperture radar (SAR) imagery, although the resulting straight flow lines do not follow the observed flow paths well. For Bingham flow, the model requires low values of viscosity and yield strength which are more similar to pyroclastic or debris flows than basaltic lavas. Flows of 100‐m depth require 1 to 2 orders of magnitude higher values of both viscosity and yield strength than 10‐m‐deep flows. The complicated nature of the flow lines for the low velocity model suggests that FEBs were probably emplaced under variably laminar and turbulent flow conditions, where underlying topography influenced both the direction and energy of flow materials.

Journal of Geophysical Research E: Planets↗
USGS2011

Seismograms offer insight into Oklahoma City bombing

The terrorist bombing of the Alfred P. Murrah Federal Building in Oklahoma City on April 19, 1995, generated seismic waves that were recorded on two permanent seismographs about 7 and 26 km away from the bombing. The seismogram recorded at 26 km shows two low-frequency wave trains, discrete sets of oscillatory signals, that begin about 10 s apart. Public release of this record prompted speculation that each wave train was caused by a different energy source. On May 23, 1995, the U.S. Geological Survey monitored the demolition of the bomb-ravaged Federal Building with portable seismographs (Figure 1). Two wave trains were picked up again. The recordings indicate that the wave trains during both the bombing and demolition represent seismic waves traveling at different velocities. We conclude that the two wave trains recorded during the bombing are consistent with a single impulsive energy source.

Oklahoma↗
USGS1997

Industrially induced changes in Earth structure at the geysers geothermal area, California

Industrial exploitation is causing clearly-measurable changes in Earth structure at The Geysers geothermal area, California. Production at The Geysers peaked in the late 1980s at ∼3.5 × 10³ kg s −1 of steam and 1800 MW of electricity. It subsequently decreased by about 10% per year [ Barker et al. , 1992] because of declining reservoir pressure. The steam reservoir coincides with a strong negative anomaly (∼0.16, ∼9%) in the compressional-to-shear seismic wave speed ratio V P / V S , consistent with the expected effects of low-pressure vapor-phase pore fluid [ Julian et al. , 1996]. Between 1991 and 1994 this anomaly increased in amplitude by up to about 0.07 (∼4%). This is consistent with the expected effects of continued pressure reduction and conversion of pore water to steam as a result of exploitation. These unique results show that V P / V S tomography can easily detect saturation changes caused by exploitation of reservoirs, and is a potentially valuable technique for monitoring environmental change. They also provide geophysical observational evidence that geothermal energy is not a renewable energy source.

Geophysical Research Letters↗