USGS Science⌕ Search

SEARCH · USGS Science

Results for “Pure & Applied Geophysics”

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 163 records · Page 9Linked to original sources
USGS1999

Local tsunamis and distributed slip at the source

Variations in the local tsunami wave field are examined in relation to heterogeneous slip distributions that are characteristic of many shallow subduction zone earthquakes. Assumptions inherent in calculating the coseismic vertical displacement field that defines the initial condition for tsunami propagation are examined. By comparing the seafloor displacement from uniform slip to that from an ideal static crack, we demonstrate that dip-directed slip variations significantly affect the initial cross-sectional wave profile. Because of the hydrodynamic stability of tsunami wave forms, these effects directly impact estimates of maximum runup from the local tsunami. In most cases, an assumption of uniform slip in the dip direction significantly underestimates the maximum amplitude and leading wave steepness of the local tsunami. Whereas dip-directed slip variations affect the initial wave profile, strike-directed slip variations result in wavefront-parallel changes in amplitude that are largely preserved during propagation from the source region toward shore, owing to the effects of refraction. Tests of discretizing slip distributions indicate that small fault surface elements of dimensions similar to the source depth can acceptably approximate the vertical displacement field in comparison to continuous slip distributions. Crack models for tsunamis generated by shallow subduction zone earthquakes indicate that a rupture intersecting the free surface results in approximately twice the average slip. Therefore, the observation of higher slip associated with tsunami earthquakes relative to typical subduction zone earthquakes of the same magnitude suggests that tsunami earthquakes involve rupture of the seafloor, whereas rupture of deeper subduction zone earthquakes may be imbedded and not reach the seafloor.

Pure and Applied Geophysics↗
USGS1999

Seismic subduction of the Nazca Ridge as shown by the 1996-97 Peru earthquakes

By rupturing more than half of the shallow subduction interface of the Nazca Ridge, the great November 12, 1996 Peruvian earthquake contradicts the hypothesis that oceanic ridges subduct aseismically. The mainshock's rupture has a length of about 200 km and has an average slip of about 1.4 m. Its moment is 1.5 x 1028 dyne-cm and the corresponding M(w) is 8.0. The mainshock registered three major episodes of moment release as shown by a finite fault inversion of teleseismically recorded broadband body waves. About 55% of the mainshock's total moment release occurred south of the Nazca Ridge, and the remaining moment release occurred at the southern half of the subduction interface of the Nazca Ridge. The rupture south of the Nazca Ridge was elongated parallel to the ridge axis and extended from a shallow depth to about 65 km depth. Because the axis of the Nazca Ridge is at a high angle to the plate convergence direction, the subducting Nazca Ridge has a large southwards component of motion, 5 cm/yr parallel to the coast. The 900-1200 m relief of the southwards sweeping Nazca Ridge is interpreted to act as a 'rigid indenter,' causing the greatest coupling south of the ridge's leading edge and leading to the large observed slip. The mainshock and aftershock hypocenters were relocated using a new procedure that simultaneously inverts local and teleseismic data. Most aftershocks were within the outline of the Nazca Ridge. A three-month delayed aftershock cluster' occurred at the northern part of the subducting Nazca Ridge. Aftershocks were notably lacking at the zone of greatest moment release, to the south of the Nazca Ridge. However, a lone foreshock at the southern end of this zone, some 140 km downstrike of the mainshock's epicenter, implies that conditions existed for rupture into that zone. The 1996 earthquake ruptured much of the inferred source zone of the M(w) 7.9-8.2 earthquake of 1942, although the latter was a slightly larger earthquake. The rupture zone of the 1996 earthquake is immediately north of the seismic gap left by the great earthquakes (M(w) ~8.8-9.1) of 1868 and 1877. The M(w) 8.0 Antofagasta earthquake of 1995 occurred at the southern end of this great seismic gap. The M(w) 8.2 deep-focus Bolivian earthquake of 1994 occurred directly downdip of the 1868 portion of that gap. The recent occurrence of three significant earthquakes on the periphery of the great seismic gap of the 1868 and 1877 events, among other factors, may signal an increased seismic potential for that zone.

Pure and Applied Geophysics↗
USGS2000

Reevaluation of tsunami formation by debris avalanche at Augustine Volcano, Alaska

Debris avalanches entering the sea at Augustine Volcano, Alaska have been proposed as a mechanism for generating tsunamis. Historical accounts of the 1883 eruption of the volcano describe 6- to 9-meter-high waves that struck the coastline at English Bay (Nanwalek), Alaska about 80 kilometers east of Augustine Island. These accounts are often cited as proof that volcanigenic tsunamis from Augustine Volcano are significant hazards to the coastal zone of lower Cook Inlet. This claim is disputed because deposits of unequivocal tsunami origin are not evident at more than 50 sites along the lower Cook Inlet coastline where they might be preserved. Shallow water (<25 m) around Augustine Island, in the run-out zone for debris avalanches, limits the size of an avalanche-caused wave. If the two most recent debris avalanches, Burr Point (A.D. 1883) and West Island (<500 yr. B.P.) were traveling at velocities in the range of 50 to 100 meters per second, the kinetic energy of the avalanches at the point of impact with the ocean would have been between 1014 and 1015 joules. Although some of this energy would be dissipated through boundary interactions and momentum transfer between the avalanche and the sea, the initial wave should have possessed sufficient kinetic energy to do geomorphic work (erosion, sediment transport, formation of wave-cut features) on the coastline of lowwer Cook Inlet. Because widespread evidence of the effects of large waves cannot be found, it appears that the debris avalanches could not have been traveling very fast when they entered the sea, or they happened during low tide and displaced only small volumes of water. In light of these results, the hazard from volcanigenic tsunamis from Augustine Volcano appears minor, unless a very large debris avalanche occurs at high tide.

Pure and Applied Geophysics↗
USGS2001

Mapping crustal heterogeneity using Lg propagation efficiency throughout the Middle East, Mediterranean, Southern Europe and Northern Africa

In this paper we describe a technique for mapping the lateral variation of Lg characteristics such as Lg blockage, efficient Lg propagation, and regions of very high attenuation in the Middle East, North Africa, Europe and the Mediterranean regions. Lg is used in a variety of seismological applications from magnitude estimation to identification of nuclear explosions for monitoring compliance with the Comprehensive Nuclear-Test-Ban Treaty (CTBT). These applications can give significantly biased results if the Lg phase is reduced or blocked by discontinuous structure or thin crust. Mapping these structures using quantitative techniques for determining Lg amplitude attenuation can break down when the phase is below background noise. In such cases Lg blockage and inefficient propagation zones are often mapped out by hand. With our approach, we attempt to visually simplify this information by imaging crustal structure anomalies that significantly diminish the amplitude of Lg. The visualization of such anomalies is achieved by defining a grid of cells that covers the entire region of interest. We trace Lg rays for each event/ station pair, which is simply the great circle path, and attribute to each cell a value equal to the maximum value of the Lg/P-coda amplitude ratio for all paths traversing that particular cell. The resulting map, from this empirical approach, is easily interpreted in terms of crustal structure and can successfully image small blockage features often missed by analysis of raypaths alone. This map can then be used to screen out events with blocked Lg prior to performing Q tomography, and to avoid using Lg-based methods of event identification for the CTBT in regions where they cannot work. For this study we applied our technique to one of the most tectonically complex regions on the earth. Nearly 9000 earthquake/station raypaths, traversing the vast region comprised of the Middle East, Mediterranean, Southern Europe and Northern Africa, have been analyzed. We measured the amplitude of Lg relative to the P-coda and mapped the lateral variation of Lg propagation efficiency. With the relatively dense coverage provided by the numerous crossing paths we are able to map out the pattern of crustal heterogeneity that gives rise to the observed character of Lg propagation. We observe that the propagation characteristics of Lg within the region of interest are very complicated but are readily correlated with the different tectonic environments within the region. For example, clear strong Lg arrivals are observed for paths crossing the stable continental interiors of Northern Africa and the Arabian Shield. In contrast, weakened to absent Lg is observed for paths crossing much of the Middle East, and Lg is absent for paths traversing the Mediterranean. Regions that block Lg transmission within the Middle East are very localized and include the Caspian Sea, the Iranian Plateau and the Red Sea. Resolution is variable throughout the region and strongly depends on the distribution of seismicity and recording stations. Lg propagation is best resolved within the Middle East where regions of crustal heterogeneity on the order of 100 km are imaged (e.g., South Caspian Sea and Red Sea). Crustal heterogeneity is resolvable but is poorest in seismically quiescent Northern Africa.

Pure and Applied Geophysics↗
USGS2003

Simulation of ground motion using the stochastic method

A simple and powerful method for simulating ground motions is to combine parametric or functional descriptions of the ground motion's amplitude spectrum with a random phase spectrum modified such that the motion is distributed over a duration related to the earthquake magnitude and to the distance from the source. This method of simulating ground motions often goes by the name "the stochastic method." It is particularly useful for simulating the higher-frequency ground motions of most interest to engineers (generally, f>0.1 Hz), and it is widely used to predict ground motions for regions of the world in which recordings of motion from potentially damaging earthquakes are not available. This simple method has been successful in matching a variety of ground-motion measures for earthquakes with seismic moments spanning more than 12 orders of magnitude and in diverse tectonic environments. One of the essential characteristics of the method is that it distills what is known about the various factors affecting ground motions (source, path, and site) into simple functional forms. This provides a means by which the results of the rigorous studies reported in other papers in this volume can be incorporated into practical predictions of ground motion.

Pure and Applied Geophysics↗
USGS2006

Special issue: Terrestrial fluids, earthquakes and volcanoes: The Hiroshi Wakita volume I

Terrestrial Fluids, Earthquakes and Volcanoes: The Hiroshi Wakita Volume I is a special publication to honor Professor Hiroshi Wakita for his scientific contributions. This volume consists of 17 original papers dealing with various aspects of the role of terrestrial fluids in earthquake and volcanic processes, which reflect Prof. Wakita’s wide scope of research interests. Professor Wakita co-founded the Laboratory for Earthquake Chemistry in 1978 and served as its director from 1988 until his retirement from the university in 1997. He has made the laboratory a leading world center for studying earthquakes and volcanic activities by means of geochemical and hydrological methods. Together with his research team and a number of foreign guest researchers that he attracted, he has made many significant contributions in the above-mentioned scientific fields of interest. This achievement is a testimony to not only his scientific talent, but also his enthusiasm, his open mindedness, and his drive in obtaining both human and financial support.

Pure and Applied Geophysics↗
USGS2015

The Boulder magnetic observatory

The Boulder magnetic observatory has, since 1963, been operated by the Geomagnetism Program of the U.S. Geological Survey in accordance with Bureau and national priorities. Data from the observatory are used for a wide variety of scientific purposes, both pure and applied. The observatory also supports developmental projects within the Geomagnetism Program and collaborative projects with allied geophysical agencies.

Colorado↗
USGS2007

Integration of seafloor point data in usSEABED

Sediments of the beach, nearshore, and continental shelves record a complex interplay of processes including wave energy and direction , currents, beach erosion or accretion, bluff or cliff retreat, fluvial input, sediment longshore and cross-shelf transport processes, contaminant content and transport, sediment sources and sinks, and others. In turn, sediments and rocks modify wave patterns, affect recreation and tourism, and provide habitat for fish, epifauna, and infauna. Character of the surficial seafloor also influences navigation, commercial and recreational fishing and gathering of other food sources, communication, piplines, national defense, and provides geologic resources including sand and gravel aggregates, minerals, and real or potential energy sources. The beaches, nearshore, and continental margins fall under overlapping levels of managerial responsibility between Federal, State, regional, and local government agencies and consortia. In addition, universities and other academic institutions investigate these places for pure or applied scientific reasons. Mapping is usually the first step in understanding any issue and is often comprised of remotely gathered geophysical data such as bathymetry and backscatter imagery, and groundtruthing; that is, the collection of physical and virtual samples to tie the remotely gathered data to reality. The physical samples are described and (or) carefully analyzed for grain-size information -- which records both the site's physical conditions and geologic past -- and commonly, for constituent components such as mineral and rock types (to determine onland sources and in situ chemical processes), carbonate and organic content and microfossils (for biological and oceanographic influences), and structure such as layering and bioturbation (for physical influences). The samples may also be subjected to physical tests such as comp[action analyses, liquefaction or plasticity limits, ans other parameters important when considering construction of offshore structures. In recent years, virtual sampling of the seafloor has become popular, through the use of towed video or photographic equipment and the addition of camera to oceanographic equipment such as corers and tripods. Before about ten years ago, most maps were made by hand. Recently, with the advent of desktop GIS packages, map making and resource analysis can be done nearly "on-the-fly" if geographically located data exist. While the problems of projection, scale, and resolution of digitized paper maps are commonly known amongst GIS-users, access to the original underlying point data allows for maps to be regenerated for digital use using statistically proven methods, provides increasing data density by including multiple studies, as well as allows the point data to be used in other ways than just mapping. These point data may be available in raw or refined or in worded descriptions. Raw data such as granulometric analyses can be manipulated through the use of known equations or empirical relationships to provide information about other parameters of the sediment, such as mean grainsize, sorting, erodability, or rugosity. If refined data are presented such as gravel, sand, and mud percentages, the parameter noted earlier may be estimated. In the case of worded descriptions, values for geologic terms can be assigned, for example, "fine sand" equate to 0.2 mm sized particles, to provide numeric terms for GIS or modeling purposes.

Conference Paper↗
USGS1990

Selected translations of the Russian literature on the electrogeochemical sampling technique called CHIM (chastichnoe izvlechennye metallov)

As part of our research into new methods for the assessment of mineral deposits, the U.S. Geological Survey has recently begun investigation of the CHIM method. As part of our studies, translation of a Russian manual on the CHIM methodology and eight articles from the Russian literature were transit ted to provide background for our own research. The translations were done by Earth Science Translation Services of Albuquerque, New Mexico, and are presented as received, without editing on our part. Below is a bibliography of the translated articles. For approximately the past 20 years Russian geoscientists have been applying an electrogeochemical sampling technique given the Russian acronym CHIM, derived from Chastichnoe Izvlechennye Metallov which translates as "partial extraction of metals". In this technique a direct current is introduced into the earth through collector electrodes similar to "porous pots" used in electrical geophysical applications. The solution in the cathode is dilute nitric acid, and current is passed through the cathode for times ranging from 6 hours to 48 hours or more. Electrical connections to the nitric acid are made through an inner conductor that is typically spectroscopically pure graphite. At the cathode, mobile cations collect on the graphite or in the nitric acid solution, both of which serve as the geochemical sampling media. These media are then analyzed by appropriate methods for the ions of interest. In most applications of the CHIM method only mobile cations are sampled, although Russian literature does refer to collection of anions as well. More recently the CHIM method has been applied by the Peoples Republic of China and the Indian Geological Survey. The literature indicates that the method has advantages over other geochemical sampling techniques by providing increased sensitivity to the metals being searched for, especially where deposits are covered by substantial overburden. In some cases success has been claimed with overburden in excess of 500 meters. The technique appears to have been applied principally to exploration for base- and precious-metal deposits, but does not appear to be limited to these. References are made in the literature to its application in the search for nickel, cobalt, molybdenum, uranium, tin, REE, tungsten, berylium, and oil and gas.

Open-File Report↗
USGS1999

The color of the Martian sky and its influence on the illumination of the Martian surface

The dust in the atmosphere above the Mars Pathfinder landing site produced a bright, red sky that increases in redness toward the horizon at midday. There is also evidence for an absorption band in the scattered light from the sky at 860 nm. A model of the sky brightness has been developed [ Markiewicz et al. , this issue] and tested against Imager for Mars Pathfinder (IMP) observations of calibration targets on the lander. The resulting model has been used to quantify the total diffuse flux onto a surface parallel to the local level for several solar elevation angles and optical depths. The model shows that the diffuse illumination in shadowed areas is strongly reddened while areas illuminated directly by the Sun (and the blue forward scattering peak) see a more solar-type spectrum, in agreement with Viking and IMP observations. Quantitative corrections for the reddening in shadowed areas are demonstrated. It is shown quantitatively that the unusual appearance of the rock Yogi (the east face of which appeared relatively blue in images taken during the morning but relatively red during the afternoon) can be explained purely by the changing illumination geometry. We conclude that any spectrophotometric analysis of surfaces on Mars must take into account the diffuse flux. Specifically, the reflectances of surfaces viewed under different illumination geometries cannot be investigated for spectral diversity unless a correction has been applied which removes the influence of the reddened diffuse flux.

Journal of Geophysical Research E: Planets↗
USGS2011

High-frequency Born synthetic seismograms based on coupled normal modes

High-frequency and full waveform synthetic seismograms on a 3-D laterally heterogeneous earth model are simulated using the theory of coupled normal modes. The set of coupled integral equations that describe the 3-D response are simplified into a set of uncoupled integral equations by using the Born approximation to calculate scattered wavefields and the pure-path approximation to modulate the phase of incident and scattered wavefields. This depends upon a decomposition of the aspherical structure into smooth and rough components. The uncoupled integral equations are discretized and solved in the frequency domain, and time domain results are obtained by inverse Fourier transform. Examples show the utility of the normal mode approach to synthesize the seismic wavefields resulting from interaction with a combination of rough and smooth structural heterogeneities. This approach is applied to an ∼4 Hz shallow crustal wave propagation around the site of the San Andreas Fault Observatory at Depth (SAFOD).

Geophysical Journal International↗
USGS2011

Constraints on the long‐period moment‐dip tradeoff for the Tohoku earthquake

Since the work of Kanamori and Given (1981), it has been recognized that shallow, pure dip‐slip earthquakes excite long‐period surface waves such that it is difficult to independently constrain the moment ( M 0 ) and the dip ( δ ) of the source mechanism, with only the product M 0 sin(2 δ ) being well constrained. Because of this, it is often assumed that the primary discrepancies between the moments of shallow, thrust earthquakes are due to this moment‐dip tradeoff. In this work, we quantify how severe this moment‐dip tradeoff is depending on the depth of the earthquake, the station distribution, the closeness of the mechanism to pure dip‐slip, and the quality of the data. We find that both long‐period Rayleigh and Love wave modes have moment‐dip resolving power even for shallow events, especially when stations are close to certain azimuths with respect to mechanism strike and when source depth is well determined. We apply these results to USGS W phase inversions of the recent M9.0 Tohoku, Japan earthquake and estimate the likely uncertainties in dip and moment associated with the moment‐ dip tradeoff. After discussing some of the important sources of moment and dip error, we suggest two methods for potentially improving this uncertainty.

Tohoku↗
USGS2023

Summary of the history and research of the U.S. Geological Survey gas hydrate properties laboratory in Menlo Park, California, active from 1993 to 2022

The U.S. Geological Survey (USGS) Clathrate Hydrate Properties Project, active from 1993 to 2022 in Menlo Park, California, stemmed from an earlier project on the properties of planetary ices supported by the National Aeronautics and Space Administration’s (NASA’s) Planetary Geology and Geophysics Program. We took a material science approach in both projects, emphasizing chemical purity of samples, having controlled grain size and grain texture, and having verified crystal structures and phase relations. A foundational contribution from our USGS Gas Hydrate Properties Laboratory (GHPL) was in demonstrating the ability to reproducibly create such pure clathrate hydrate samples for study. Clathrate sample synthesis was achieved by heating sieved and weighed pure granular water ice in the presence of cold clathrate-forming gas or liquid. During heating, the ice melts at the grain scale and reacts with the gas to form clathrate. The resulting material has the desired uniformity and purity, with known intergranular porosity; our subsequent measurements showed that these clathrates exhibited the established clathrate structures and phase relations. This novel synthesis method was successful in creating clathrates of pure methane, ethane, propane, carbon dioxide, and multi-component gases. By mixing sand or silt with granular ice, we were also able to make clathrate-sediment aggregates with controlled grain textures. This simple method, adopted by many others in the community, permitted us to measure the physical and chemical properties of well-characterized and well-crystallized clathrates and clathrate/sediment aggregates. At about the same time, we adapted conventional scanning electron microscopy to cryogenic conditions for analysis of grain-scale characteristics of clathrates made in the GHPL as well as those collected from nature by drill core. The uniformity and reproducibility of our samples also allowed us to investigate how clathrates respond to environmental changes in chemistry, temperature, and pressure: we measured chemical exchange rates with dissolved gas species—such as noble gases and chlorofluorocarbons—as well as rates of clathrate dissolution and decomposition. These advances include the first accurate mapping of the conditions that promote the remarkable process of “anomalous preservation” at room pressure, a metastability that offers potential application for low-cost and safe transportation of natural gas from gas fields far from pipelines. Another advancement stemming from the GHPL was the compaction of as-synthesized porous clathrates to nearly full density by applying external pressure using three different techniques. Compaction allows for high-accuracy measurements of many fundamental physical and chemical properties of these materials, such as elastic wavespeeds and moduli, complete thermal properties, decomposition rates, thermal expansion, and clathrate equations of state. These properties and others, in turn, have helped USGS scientists to interpret geophysical well logs and active geophysical surveys, as well as model the rates of gas production from hydrate deposits in nature. Studying this class of icy minerals that occur in abundance on Earth and in the outer solar system has been a fascinating laboratory journey. Here, we summarize the history and major findings of the USGS GHPL in Menlo Park, including both in-house research as well as findings from the synergistic collaborations with other agencies and institutes that were key to the success of our laboratory. The Menlo Park GHPL was more formally incorporated within the USGS Gas Hydrates Project, a collaboration among multiple USGS Science Centers, in the early 2000s under the leadership of Deborah Hutchinson, and now under the leadership of Carolyn Ruppel and Timothy Collett.

Open-File Report↗
USGS1999

Steady-state flow of solid CO2: Preliminary results

To help answer the question of how much solid CO 2 exists in the Martian south polar cap, we performed a series of laboratory triaxial deformation experiments at constant displacement rate in compression on jacketed cylinders of pure, polycrystalline CO 2 . Test conditions were temperatures 150 < T < 190 K, hydrostatic confining pressures 5≤ P ≤40 MPa, and strain rates 4.5×10 −8 ≤ ε ≤4.3×10 −4 s −1 . Most of the measurements follow a constitutive law of the form ε = Aσ n exp(−Q/RT), where σ is the applied differential stress, R is the gas constant, and the other constants have values as follows: A = 10 3 86 MPa −n s −1 , n = 5.6, and Q = 33 kJ/mol. Solid CO 2 is markedly weaker than water ice. Our results suggest that the south polar cap on Mars is unlikely to be predominately solid CO 2 , because the elevation and estimated age of the cap is difficult to reconcile with the very weak rheology of the material.

Geophysical Research Letters↗
USGS2012

Friction of ice

The frictional strength of ice seems to be well below that for all other rocks. Triaxial testing of cylinders of pure water ice containing a 45° inclined sawcut, at temperatures of 77 ≤ T ≤ 115 K and confining pressures 0.1 ≤ P ≤ 250 MPa, reveals the frictional laws τ = 0.20 σ n + 8.3 MPa for P ≥ 10 MPa and τ = 0.55 σ n + 1.0 MPa for P ≤ 5 MPa, where τ and σ n are the shear and normal stresses on the sawcut surface. Friction is independent of T and, over the one order of magnitude tested, of average sliding velocity. The sliding behavior is invariably stick slip. Examination of sliding surfaces shows only minor gouge development. A curious arrangement of densely packed short vertical fractures occurs primarily in samples with anomalously low strength. These results apply to models of near‐surface tectonic activity on Jupiter's three large icy moons, Europa, Callisto, and Ganymede. They indicate, for example, that a global expansion on Ganymede of 3 linear percent will cause extensional movement on preexisting faults at depths to 7±3 km.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

3-D joint geodetic and strong-motion finite fault inversion of the 2008 May 12, Wenchuan, China Earthquake

We present a source inversion of the 2008 Wenchuan, China earthquake, using strong-motion waveforms and geodetic offsets together with three-dimensional synthetic ground motions. We applied the linear multiple time window technique considering geodetic and dynamic Green's functions computed with the finite element method and the reciprocity and Strain Green’s Tensor formalism. All ground motion estimates, valid up to 1 Hz, accounted for three-dimensional effects, including the topography and the geometry of the Beichuan and Pengguan faults. Our joint inversion has a higher moment (M0) than a purely geodetic inversion and the slip distribution presents differences when compared to one-dimensional model source inversions. The moment is estimated to be M0=1.2x1021 Nm, slightly larger than other works. Our results show that considering a complex 3D structure reduces the size of large areas of 10 m slip or greater by distributing it in wider zones, with reduced slips, in the central portion of the Beichuan and the Pengguan faults. Finally, we compare our source with a relocated aftershock catalog and conclude that the 4-5 m slip contours approximately bound the absence or presence of aftershocks.

Wenchuan↗
USGS2018

Impact of pore fluid chemistry on fine-grained sediment fabric and compressibility

Fines, defined here as grains or particles, less than 75 μm in diameter, exist nearly ubiquitously in natural sediment, even those classified as coarse. Macroscopic sediment properties, such as compressibility, which relates applied effective stress to the resulting sediment deformation, depend on the fabric of fines. Unlike coarse grains, fines have sizes and masses small enough to be more strongly influenced by electrical interparticle forces than by gravity. These electrical forces acting through pore fluids are influenced by pore fluid chemistry changes. Macroscopic property dependence on pore fluid chemistry must be accounted for in sediment studies involving subsurface flow and sediment stability analyses, as well as in engineered flow situations such as groundwater pollutant remediation, hydrocarbon migration, or other energy resource extraction applications. This study demonstrates how the liquid limit‐based electrical sensitivity index can be used to predict sediment compressibility changes due to pore fluid chemistry changes. Laboratory tests of electrical sensitivity, sedimentation, and compressibility illustrate mechanisms linking microscale and macroscale processes for selected pure, end‐member fines. A specific application considered here is methane extraction via depressurization of gas hydrate‐bearing sediment, which causes a dramatic pore water salinity drop concurrent with sediment being compressed by the imposed effective stress increase.

Journal of Geophysical Research↗
USGS2012

Surface heat flow and CO2 emissions within the Ohaaki hydrothermal field, Taupo Volcanic Zone, New Zealand

Carbon dioxide emissions and heat flow have been determined from the Ohaaki hydrothermal field, Taupo Volcanic Zone (TVZ), New Zealand following 20 a of production (116 MW e ). Soil CO 2 degassing was quantified with 2663 CO 2 flux measurements using the accumulation chamber method, and 2563 soil temperatures were measured and converted to equivalent heat flow (W m −2 ) using published soil temperature heat flow functions. Both CO 2 flux and heat flow were analysed statistically and then modelled using 500 sequential Gaussian simulations. Forty subsoil CO 2 gas samples were also analysed for stable C isotopes. Following 20 a of production, current CO 2 emissions equated to 111 ± 6.7 T/d. Observed heat flow was 70 ± 6.4 MW, compared with a pre-production value of 122 MW. This 52 MW reduction in surface heat flow is due to production-induced drying up of all alkali–Cl outflows (61.5 MW) and steam-heated pools (8.6 MW) within the Ohaaki West thermal area (OHW). The drying up of all alkali–Cl outflows at Ohaaki means that the soil zone is now the major natural pathway of heat release from the high-temperature reservoir. On the other hand, a net gain in thermal ground heat flow of 18 MW (from 25 MW to 43.3 ± 5 MW) at OHW is associated with permeability increases resulting from surface unit fracturing by production-induced ground subsidence. The Ohaaki East (OHE) thermal area showed no change in distribution of shallow and deep soil temperature contours despite 20 a of production, with an observed heat flow of 26.7 ± 3 MW and a CO 2 emission rate of 39 ± 3 T/d. The negligible change in the thermal status of the OHE thermal area is attributed to the low permeability of the reservoir beneath this area, which has limited production (mass extraction) and sheltered the area from the pressure decline within the main reservoir. Chemistry suggests that although alkali–Cl outflows once contributed significantly to the natural surface heat flow (∼50%) they contributed little (<1%) to pre-production CO 2 emissions due to the loss of >99% of the original CO 2 content due to depressurisation and boiling as the fluids ascended to the surface. Consequently, the soil has persisted as the major (99%) pathway of CO 2 release to the atmosphere from the high temperature reservoir at Ohaaki. The CO 2 flux and heat flow surveys indicate that despite 20 a of production the variability in location, spatial extent and magnitude of CO 2 flux remains consistent with established geochemical and geophysical models of the Ohaaki Field. At both OHW and OHE carbon isotopic analyses of soil gas indicate a two-stage fractionation process for moderate-flux (>60 g m −2 d −1 ) sites; boiling during fluid ascent within the underlying reservoir and isotopic enrichment as CO 2 diffuses through porous media of the soil zone. For high-flux sites (>300 g m −2 d −1 ), the δ 13 CO 2 signature (−7.4 ± 0.3‰ OHW and −6.5 ± 0.6‰ OHE) is unaffected by near-surface (soil zone) fractionation processes and reflects the composition of the boiled magmatic CO 2 source for each respective upflow. Flux thresholds of <30 g m −2 d −1 for purely diffusive gas transport, between 30 and 300 g m −2 d −1 for combined diffusive–advective transport, and ⩾300 g m −2 d −1 for purely advective gas transport at Ohaaki were assigned. δ 13 CO 2 values and cumulative probability plots of CO 2 flux data both identified a threshold of ∼15 g m −2 d −1 by which background (atmospheric and soil respired) CO 2 may be differentiated from hydrothermal CO 2 .

Taupo Volcanic Zone↗