USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geothermics”

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 991 records · Page 55Linked to original sources

Global volcanic hazards and risk

An estimated 800 million people live within 100 km of an active volcano in 86 countries and additional overseas territories worldwide [see Chapter 4 and Appendix B]1. Volcanoes are compelling evidence that the Earth is a dynamic planet characterised by endless change and renewal. Humans have always found volcanic activity fascinating and have often chosen to live close to volcanoes, which commonly provide favourable environments for life. Volcanoes bring many benefits to society: eruptions fertilise soils; elevated topography provides good sites for infrastructure (e.g. telecommunications on elevated ground); water resources are commonly plentiful; volcano tourism can be lucrative; and volcanoes can acquire spiritual, aesthetic or religious significance. Some volcanoes are also associated with geothermal resources, making them a target for exploration and a potential energy resource. Much of the time volcanoes are not a threat because they erupt very infrequently or because communities have become resilient to frequently erupting volcanoes. However, there is an everpresent danger of a long-dormant volcano re-awakening or of volcanoes producing anomalously large or unexpected eruptions. Volcanic eruptions can cause loss of life and livelihoods in exposed communities, damage or disrupt critical infrastructure and add stress to already fragile environments. Their impacts can be both short-term, e.g. physical damage, and long-term, e.g. sustained or permanent displacement of populations. The risk from volcanic eruptions and their attendant hazards is often underestimated beyond areas within the immediate proximity of a volcano. For example, volcanic ash hazards can have effects hundreds of kilometres away from the vent and have an adverse impact on human and animal health, infrastructure, transport, agriculture and horticulture, the environment and economies. The products of volcanism and their impacts can extend beyond country borders, to be regional and even global in scale. Although known historical loss of life from volcanic eruptions (since 1600 AD about 280,000 fatalities are recorded, Auker et al. (2013)) is modest compared to other major natural hazards, volcanic eruptions can be catastrophic for exposed communities. In 1985 the town of Armero in Colombia was buried by lahars (volcanic mudflows) with more than 21,000 fatalities due to relatively small explosive eruptions at the summit of Nevado del Ruiz volcano that partially melted a glacier (Voight, 1990).

Book chapter↗

Virtual special issue of recent advances on gas hydrates scientific drilling in Alaska

Gas hydrate refers to a non-stoichiometric clathrate that forms spontaneously in the natural environment whenever sufficient quantities of gases of appropriate size (most commonly methane) interact with abundant water under specific conditions of temperature and pressure. (1,2) Such conditions occur wherever the shallow geothermal gradient has been suppressed by either deepwater or thick permafrost, allowing for relatively low temperatures to coexist with elevated pressures. The volume of gas hydrate on Earth is difficult to constrain, (3) but it is sufficient that gas hydrate is a meaningful potential component of (1) the long-term natural cycling of carbon, (2) the nearer term environmental changes in response to warming climates, (4) and (3) future energy supply systems. Since the initial recognition of gas hydrate as an abundant component in nature in the late 1960s, a series of scientific drilling expeditions conducted by both the Integrated Ocean Discovery Program (and successors) and national research and development programs in Canada, Japan, China, the United States, South Korea, and others (5) have explored the occurrence and nature of gas hydrates. In particular, the desire for expanded energy supply options to support the economic development and energy security for nations around the globe is currently motivating a broad range of laboratory and numerical simulation studies in support of an ongoing series of field-based scientific tests of the potential commercial viability of gas extraction from natural gas hydrate deposits. (6)

Journal of Energy & Fuels↗

Salinity trends in a groundwater system supplemented by 50 years of imported Colorado River water

The Indio subbasin of the Coachella Valley is a desert area of southern California where a growing population depends primarily on groundwater for drinking and agricultural uses. The aquifer system has been supplemented with Colorado River water through managed recharge and widespread irrigation since the mid-20th century. We use a combination of geochemical modeling and trend analysis to identify changes in total dissolved solids through time, elucidate the sources of dissolved solids, and quantify the extent of contributions from those sources throughout the Indio subbasin. We conclude that recharged Colorado River water is the primary source and driver of increasing salinity, particularly in areas immediately downgradient from the recharge locations and in the eastern part of the subbasin away from the recharge ponds due to irrigation using imported water. Other contributions of dissolved solids to groundwater resources include geothermal waters, wastewater effluent, and agricultural return flow, although their effects are more localized. This study presents a broadly applicable framework for identifying sources of dissolved solids in groundwater wells and salinity trends at a regional scale in a large data set.

California↗

New Method for Electrical Conductivity Temperature Compensation

Electrical conductivity (κ) measurements of natural waters are typically referenced to 25 °C (κ 25 ) using standard temperature compensation factors (α). For acidic waters (pH < 4), this can result in a large κ 25 error (δκ 25 ). The more the sample temperature departs from 25 °C, the larger the potential δκ 25 . For pH < 4, the hydrogen ion transport number becomes substantial and its mode of transport is different from most other ions resulting in a different α. A new method for determining α as a function of pH and temperature is presented. Samples with varying amounts of H 2 SO 4 and NaCl were used to develop the new α, which was then applied to 65 natural water samples including acid mine waters, geothermal waters, seawater, and stream waters. For each sample, the κ and pH were measured at several temperatures from 5 to 90 °C and κ 25 was calculated. The δκ 25 ranged from −11 to 9% for the new method as compared to −42 to 25% and −53 to 27% for the constant α (0.019) and ISO-7888 methods, respectively. The new method for determining α is a substantial improvement for acidic waters and performs as well as or better than the standard methods for circumneutral waters.

Environmental Science & Technology↗

Hydrologic processes in deep vadose zones in interdrainage arid environments

A unifying theory for the hydrology of desert vadose zones is particularly timely considering the rising population and water stresses in arid and semiarid regions. Conventional models cannot reconcile the apparent discrepancy between upward flow indicated by hydraulic gradient data and downward flow suggested by environmental tracer data in deep vadose zone profiles. A conceptual model described here explains both hydraulic and tracer data remarkably well by incorporating the hydrologic role of desert plants that encroached former juniper woodland 10 to 15 thousand years ago in the southwestern United States. Vapor transport also plays an important role in redistributing moisture through deep soils, particularly in coarse-grained sediments. Application of the conceptual model to several interdrainage arid settings reproduces measured matric potentials and chloride accumulation by simulating the transition from downward flow to upward flow just below the root zone initiated by climate and vegetation change. Model results indicate a slow hydraulic drying response in deep vadose zones that enables matric potential profiles to be used to distinguish whether precipitation episodically percolated below the root zone or was completely removed via evapotranspiration during the majority of the Holocene. Recharge declined dramatically during the Holocene in interdrainage basin floor settings of arid and semiarid basins. Current flux estimates across the water table in these environmental settings, are on the order of 0.01 to 0.1 mm yr -1 and may be recharge (downward) or discharge (upward) depending on vadose zone characteristics, such as soil texture, geothermal gradient, and water table depth. In summary, diffuse recharge through the basin floor probably contributes only minimally to the total recharge in arid and semiarid basins.

Water Science and Application↗

Seismovolcanic signals at Deception Island volcano, Antarctica: Wave field analysis and source modeling

The seismovolcanic signals associated with the volcanic activity of Deception Island (Antarctica), recorded during three Antarctic summers (1994–1995, 1995–1996 and 1996–1997), are analyzed using a dense small-aperture (500 m) seismic array. The visual and spectral classification of the seismic events shows the existence of long-period and hybrid isolated seismic events, and of low-frequency, quasi-monochromatic and spasmodic continuous tremors. All spectra have the highest amplitudes in the frequency band between 1 and 4 Hz, while hybrids and spasmodic tremors have also significant amplitudes in the high-frequency band (4–10 Hz). The array analysis indicates that almost all the well-correlated low-frequency signals share similar array parameters (slowness and back azimuth) and have the same source area, close to the array site. The polarization analysis shows that phases at high-frequency are mostly composed of P waves, and those phases dominated by low frequencies can be interpreted as surface waves. No clear shear waves are evidenced. From the energy evaluation, we have found that the reduced displacement values for surface and body waves are confined in a narrow interval. Volcano-tectonic seismicity is located close to the array, at a depth shallower than 1 km. The wave-field properties of the seismovolcanic signals allow us to assume a unique source model, a shallow resonating fluid-filled crack system at a depth of some hundreds of meters. All of the seismic activity is interpreted as the response of a reasonably stable stationary geothermal process. The differences observed in the back azimuth between low and high frequencies are a near-field effect. A few episodes of the degassification process in an open conduit were observed and modeled with a simple organ pipe.

Journal of Geophysical Research: Solid Earth↗

New constraints on mechanisms of remotely triggered seismicity at Long Valley Caldera

Regional-scale triggering of local earthquakes in the crust by seismic waves from distant main shocks has now been robustly documented for over a decade. Some of the most thoroughly recorded examples of repeated triggering of a single site from multiple, large earthquakes are measured in geothermal fields of the western United States like Long Valley Caldera. As one of the few natural cases where the causality of an earthquake sequence is apparent, triggering provides fundamental constraints on the failure processes in earthquakes. We show here that the observed triggering by seismic waves is inconsistent with any mechanism that depends on cumulative shaking as measured by integrated energy density. We also present evidence for a frequency-dependent triggering threshold. On the basis of the seismic records of 12 regional and teleseismic events recorded at Long Valley Caldera, long-period waves (>30 s) are more effective at generating local seismicity than short-period waves of comparable amplitude. If the properties of the system are stationary over time, the failure threshold for long-period waves is ~0.05 cm/s vertical shaking. Assuming a phase velocity of 3.5 km/s and an elastic modulus of 3.5 x 10 10 Pa, the threshold in terms of stress is 5 kPa. The frequency dependence is due in part to the attenuation of the surface waves with depth. Fluid flow through a porous medium can produce the rest of the observed frequency dependence of the threshold. If the threshold is not stationary with time, pore pressures that are >99.5% of lithostatic and vary over time by a factor of 4 could explain the observations with no frequency dependence of the triggering threshold.

California↗

Debris-bed friction of hard-bedded glaciers

[1] Field measurements of debris-bed friction on a smooth rock tablet at the bed of Engabreen, a hard-bedded, temperate glacier in northern Norway, indicated that basal ice containing 10% debris by volume exerted local shear traction of up to 500 kPa. The corresponding bulk friction coefficient between the dirty basal ice and the tablet was between 0.05 and 0.08. A model of friction in which nonrotating spherical rock particles are held in frictional contact with the bed by bed-normal ice flow can account for these measurements if the power law exponent for ice flowing past large clasts is 1. A small exponent (n < 2) is likely because stresses in ice are small and flow is transient. Numerical calculations of the bed-normal drag force on a sphere in contact with a flat bed using n = 1 show that this force can reach values several hundred times that on a sphere isolated from the bed, thus drastically increasing frictional resistance. Various estimates of basal friction are obtained from this model. For example, the shear traction at the bed of a glacier sliding at 20 m a-1 with a geothermally induced melt rate of 0.006 m a-1 and an effective pressure of 300 kPa can exceed 100 kPa. Debris-bed friction can therefore be a major component of sliding resistance, contradicting the common assumption that debris-bed friction is negligible. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

Tectonic stressing in California modeled from GPS observations

What happens in the crust as a result of geodetically observed secular motions? In this paper we find out by distorting a finite element model of California using GPS-derived displacements. A complex model was constructed using spatially varying crustal thickness, geothermal gradient, topography, and creeping faults. GPS velocity observations were interpolated and extrapolated across the model and boundary condition areas, and the model was loaded according to 5-year displacements. Results map highest differential stressing rates in a 200-km-wide band along the Pacific-North American plate boundary, coinciding with regions of greatest seismic energy release. Away from the plate boundary, GPS-derived crustal strain reduces modeled differential stress in some places, suggesting that some crustal motions are related to topographic collapse. Calculated stressing rates can be resolved onto fault planes: useful for addressing fault interactions and necessary for calculating earthquake advances or delays. As an example, I examine seismic quiescence on the Garlock fault despite a calculated minimum 0.1-0.4 MPa static stress increase from the 1857 M???7.8 Fort Tejon earthquake. Results from finite element modeling show very low to negative secular Coulomb stress growth on the Garlock fault, suggesting that the stress state may have been too low for large earthquake triggering. Thus the Garlock fault may only be stressed by San Andreas fault slip, a loading pattern that could explain its erratic rupture history.

Journal of Geophysical Research B: Solid Earth↗

Anisotropy, repeating earthquakes, and seismicity associated with the 2008 eruption of Okmok Volcano, Alaska

We use shear wave splitting (SWS) analysis and double-difference relocation to examine temporal variations in seismic properties prior to and accompanying magmatic activity associated with the 2008 eruption of Okmok volcano, Alaska. Using bispectrum cross-correlation, a multiplet of 25 earthquakes is identified spanning five years leading up to the eruption, each event having first motions compatible with a normal fault striking NE–SW. Cross-correlation differential times are used to relocate earthquakes occurring between January 2003 and February 2009. The bulk of the seismicity prior to the onset of the eruption on 12 July 2008 occurred southwest of the caldera beneath a geothermal field. Earthquakes associated with the onset of the eruption occurred beneath the northern portion of the caldera and started as deep as 13 km. Subsequent earthquakes occurred predominantly at 3 km depth, coinciding with the depth at which the magma body has been modeled using geodetic data. Automated SWS analysis of the Okmok catalog reveals radial polarization outside the caldera and a northwest-southeast polarization within. We interpret these polarizations in terms of a magma reservoir near the center of the caldera, which we model with a Mogi point source. SWS analysis using the same input processing parameters for each event in the multiplet reveals no temporal changes in anisotropy over the duration of the multiplet, suggesting either a short-term or small increase in stress just before the eruption that was not detected by GPS, or eruption triggering by a mechanism other than a change of stress in the system.

Alaska↗

Remotely triggered microearthquakes and tremor in central California following the 2010 M w 8.8 Chile earthquake

We examine remotely triggered microearthquakes and tectonic tremor in central California following the 2010 M w 8.8 Chile earthquake. Several microearthquakes near the Coso Geothermal Field were apparently triggered, with the largest earthquake ( M l 3.5) occurring during the large-amplitude Love surface waves. The Chile mainshock also triggered numerous tremor bursts near the Parkfield-Cholame section of the San Andreas Fault (SAF). The locally triggered tremor bursts are partially masked at lower frequencies by the regionally triggered earthquake signals from Coso, but can be identified by applying high-pass or matched filters. Both triggered tremor along the SAF and the M l 3.5 earthquake in Coso are consistent with frictional failure at different depths on critically-stressed faults under the Coulomb failure criteria. The triggered tremor, however, appears to be more phase-correlated with the surface waves than the triggered earthquakes, likely reflecting differences in constitutive properties between the brittle, seismogenic crust and the underlying lower crust.

California↗

Generation and evolution of hydrothermal fluids at Yellowstone: Insights from the Heart Lake Geyser Basin

We sampled fumaroles and hot springs from the Heart Lake Geyser Basin (HLGB), measured water and gas discharge, and estimated heat and mass flux from this geothermal area in 2009. The combined data set reveals that diverse fluids share an origin by mixing of deep solute-rich parent water with dilute heated meteoric water, accompanied by subsequent boiling. A variety of chemical and isotopic geothermometers are consistent with a parent water that equilibrates with rocks at 205°C ± 10°C and then undergoes 21% ± 2% adiabatic boiling. Measured diffuse CO 2 flux and fumarole compositions are consistent with an initial dissolved CO 2 concentration of 21 ± 7 mmol upon arrival at the caldera boundary and prior to southeast flow, boiling, and discharge along the Witch Creek drainage. The calculated advective flow from the basin is 78 ± 16 L s −1 of parent thermal water, corresponding to 68 ± 14 MW, or &ndash;1% of the estimated thermal flux from Yellowstone. Helium and carbon isotopes reveal minor addition of locally derived crustal, biogenic, and meteoric gases as this fluid boils and degasses, reducing the He isotope ratio (Rc/Ra) from 2.91 to 1.09. The HLGB is one of the few thermal areas at Yellowstone that approaches a closed system, where a series of progressively boiled waters can be sampled along with related steam and noncondensable gas. At other Yellowstone locations, steam and gas are found without associated neutral Cl waters (e.g., Hot Spring Basin) or Cl-rich waters emerge without significant associated steam and gas (Upper Geyser Basin).

Wyoming↗

Changes in permeability caused by transient stresses: field observations, experiments, and mechanisms

Oscillations in stress, such as those created by earthquakes, can increase permeability and fluid mobility in geologic media. In natural systems, strain amplitudes as small as 10&ndash;6 can increase discharge in streams and springs, change the water level in wells, and enhance production from petroleum reservoirs. Enhanced permeability typically recovers to prestimulated values over a period of months to years. Mechanisms that can change permeability at such small stresses include unblocking pores, either by breaking up permeability-limiting colloidal deposits or by mobilizing droplets and bubbles trapped in pores by capillary forces. The recovery time over which permeability returns to the prestimulated value is governed by the time to reblock pores, or for geochemical processes to seal pores. Monitoring permeability in geothermal systems where there is abundant seismicity, and the response of flow to local and regional earthquakes, would help test some of the proposed mechanisms and identify controls on permeability and its evolution.

Reviews of Geophysics↗

Rethinking the use of seabed sediment temperature profiles to trace submarine groundwater flow

Submarine groundwater fluxes across the seafloor facilitate important hydrological and biogeochemical exchanges between oceans and seabed sediment, yet few studies have investigated spatially distributed groundwater fluxes in deep‐ocean environments such as continental slopes. Heat has been previously applied as a submarine groundwater tracer using an analytical solution to a heat flow equation assuming steady state conditions and homogeneous thermal conductivity. These assumptions are often violated in shallow seabeds due to ocean bottom temperature changes or sediment property variations. Here heat tracing analysis techniques recently developed for terrestrial settings are applied in concert to examine the influences of groundwater flow, ocean temperature changes, and seabed thermal conductivity variations on deep‐ocean sediment temperature profiles. Temperature observations from the sediment and bottom ocean water on the Scotian Slope off eastern Canada are used to demonstrate how simple thermal methods for tracing groundwater can be employed if more comprehensive techniques indicate that the simplifying assumptions are valid. The spatial distribution of the inferred groundwater fluxes on the slope suggests a downward groundwater flow system with recharge occurring over the upper‐middle slope and discharge on the lower slope. We speculate that the downward groundwater flow inferred on the Scotian Slope is due to density‐driven processes arising from underlying salt domes, in contrast with upward slope systems driven by geothermal convection. Improvements in the design of future submarine hydrogeological studies are proposed for thermal data collection and groundwater flow analysis, including new equations that quantify the minimum detectable flux magnitude for a given sensor accuracy and profile length.

Water Resources Research↗

Interseismic ground deformation and fault slip rates in the greater San Francisco Bay Area from two decades of space geodetic data

The detailed spatial variations of strain accumulation and creep on major faults in the northern San Francisco Bay Area (North Bay), which are important for seismic potential and evaluation of natural hazards, remain poorly understood. Here we combine interferometric synthetic aperture radar data from the ERS‐1/2 and Envisat satellites between 1992 and 2010 with continuous and campaign GPS data to obtain a high spatial and temporal coverage of ground deformation of the North Bay. The SAR data from both ascending and descending orbits are combined to separate horizontal and vertical components of the deformation. We jointly invert the horizontal component of the mean velocities derived from these data to infer the deep strike‐slip rates on major locked faults. We use the estimated deep rates to simulate the long‐wavelength deformation due to interseismic elastic strain accumulation along these locked faults. After removing the long‐wavelength signal from the InSAR horizontal mean velocity field, we estimate fault‐parallel surface creep rates of up to 2 mm/year along the central section of the Rodgers Creek fault and surface creep rates ranging between 2 and 4 mm/year along the Concord fault. No surface creep is geodetically resolved along the West Napa and Green Valley fault zones. We identified characteristically repeating earthquakes on the Rodgers Creek fault, the West Napa fault, the Green Valley fault, and the Concord fault. Nontectonic deformation in the Geysers geothermal field and in Late Cenozoic basins (Rohnert Park and Sonoma basins) are also observed, likely due to hydrological and sediment‐compaction processes, respectively.

Journal of Geophysical Research B: Solid Earth↗

Submarine permafrost map in the arctic modelled using 1D transient heat flux (SuPerMAP)

Offshore permafrost plays a role in the global climate system, but observations of permafrost thickness, state, and composition are limited to specific regions. The current global permafrost map shows potential offshore permafrost distribution based on bathymetry and global sea level rise. As a first‐order estimate, we employ a heat transfer model to calculate the subsurface temperature field. Our model uses dynamic upper boundary conditions that synthesize Earth System Model air temperature, ice mass distribution and thickness, and global sea level reconstruction and applies globally distributed geothermal heat flux as a lower boundary condition. Sea level reconstruction accounts for differences between marine and terrestrial sedimentation history. Sediment composition and pore water salinity are integrated in the model. Model runs for 450 ka for cross‐shelf transects were used to initialize the model for circumarctic modeling for the past 50 ka. Preindustrial submarine permafrost (i.e., cryotic sediment), modeled at 12.5‐km spatial resolution, lies beneath almost 2.5 ×10 6 km 2 of the Arctic shelf. Our simple modeling approach results in estimates of distribution of cryotic sediment that are similar to the current global map and recent seismically delineated permafrost distributions for the Beaufort and Kara seas, suggesting that sea level is a first‐order determinant for submarine permafrost distribution. Ice content and sediment thermal conductivity are also important for determining rates of permafrost thickness change. The model provides a consistent circumarctic approach to map submarine permafrost and to estimate the dynamics of permafrost in the past.

Journal of Geophysical Research C: Oceans↗

Radiocarbon dating of silica sinter and postglacial hydrothermal activity in the El Tatio geyser field

The El Tatio geothermal field in the Chilean Altiplano contains hydrothermal silica sinter deposits overlaying glacial and volcanic units, providing an opportunity to constrain the timing of deglaciation and volcanic activity in an area with sparse absolute chronologies. We obtained 51 new radiocarbon ages and δ 13 C values on the organic material trapped in these sinter deposits. Based on the δ 13 C values, we exclude 29 samples for possible contamination with bacterial mats that incorporate old carbon. We infer that hydrothermal activity initiated ~27 ka ago and has been nearly continuous ever since. The ages of the oldest sinter deposits coincide with ages of moraines that stabilized after the most recent deglaciation. Whereas late Pleistocene sinters are broadly distributed in the field, Holocene deposits are found around active hydrothermal features. Although recent volcanism is absent in the vicinity of El Tatio, persistent hydrothermal discharge implies a long‐lived magmatic heat source.

El Tatio Geothermal Field↗

Characteristics of frequent dynamic triggering of microearthquakes in Southern California

Dynamic triggering of earthquakes has been reported at various fault systems. The triggered earthquakes are thought to be caused either directly by dynamic stress changes due to the passing seismic waves, or indirectly by other nonlinear processes that are initiated by the passing waves. Distinguishing these physical mechanisms is difficult because of the general lack of high‐resolution earthquake catalogs and robust means to quantitatively evaluate triggering responses, particularly, delayed responses. Here we use the high‐resolution Quake Template Matching catalog in Southern California to systematically evaluate teleseismic dynamic triggering patterns in the San Jacinto Fault Zone and the Salton Sea Geothermal Field from 2008 to 2017. We develop a new statistical approach to identify triggered cases, finding that approximately 1 out of every 5 global M w ≥ 6 earthquakes dynamically trigger microearthquakes in Southern California. The triggering responses include both instantaneous and delayed triggering, showing a highly heterogeneous pattern and indicating possible evolving triggering thresholds. We do not observe a clear peak ground velocity triggering threshold that can differentiate triggering earthquakes from non‐triggering events, but there are subtle differences in the frequency content that may possibly differentiate the earthquakes. In contrast to the depth distribution of background seismicity, the identified triggered earthquakes tend to concentrate at the edges of the seismogenic zones. Although instantaneously triggered earthquakes are likely a result of dynamic Coulomb stress changes, the cases of delayed dynamic triggering are best explained by nonlinear triggering processes, including cyclic material fatigue, accelerated transient creep, and stochastic frictional heterogeneities.

California↗