Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysical Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7Linked to original sources
USGS2018

KG²B, a collaborative benchmarking exercise for estimating the permeability of the Grimsel granodiorite - Part 1: measurements, pressure dependence and pore-fluid effects

Measuring the permeability of tight rocks remains a challenging task. In addition to the traditional sources of errors that affect more permeable formations (e.g. sample selection, non-representative specimens, disturbance introduced during sample acquisition and preparation), tight rocks can be particularly prone to solid–fluid interactions and thus more sensitive to the methods, procedures and techniques used to measure permeability. To address this problem, it is desirable to collect, for a single material, measurements obtained by different methods and pore-fluids. For that purpose a collaborative benchmarking exercise involving 24 laboratories was organized for measuring the permeability of a single low permeability material, the Grimsel granodiorite, at a common effective confining pressure (5 MPa). The objectives of the benchmark were: (i) to compare the results for a given method, (ii) to compare the results between different methods, (iii) to analyze the accuracy of each method, (iv) to study the influence of experimental conditions (especially the nature of pore fluid), (v) to discuss the relevance of indirect methods and models and finally (vi) to suggest good practice for low permeability measurements. In total 39 measurements were collected that allowed us to discuss the influence of (i) pore-fluid, (ii) measurement method, (iii) sample size and (iv) pressure sensitivity. Discarding some outliers from the bulk data set (4 out of 39) an average permeability of 1.11 × 10 −18 m² with a standard deviation of 0.57 × 10 −18 m² was obtained. The most striking result was the large difference in permeability for gas measurements compared to liquid measurements. Regardless of the method used, gas permeability was higher than liquid permeability by a factor approximately 2 ( k gas = 1.28 × 10 −18 m² compared to k liquid = 0.65 × 10 −18 m²). Possible explanations are that (i) liquid permeability was underestimated due to fluid-rock interactions (ii) gas permeability was overestimated due to insufficient correction for gas slippage and/or (iii) gases and liquids do not probe exactly the same porous networks. The analysis of Knudsen numbers shows that the gas permeability measurements were performed in conditions for which the Klinkenberg correction is sufficient. Smaller samples had a larger scatter of permeability values, suggesting that their volume were below the Representative Elementary Volume. The pressure dependence of permeability was studied by some of the participating teams in the range 1–30 MPa and could be fitted to an exponential law k = k o .exp(– γP eff ) with γ = 0.093 MPa −1 . Good practice rules for measuring permeability in tight materials are also provided.

Geophysical Journal International↗
USGS2018

KG²B, a collaborative benchmarking exercise for estimating the permeability of the Grimsel granodiorite - Part 2: modeling, microstructures and complementary data

Measuring and modelling the permeability of tight rocks remains a challenging task. In addition to the traditional sources of errors that affect more permeable formations (e.g. sample selection, non-representative specimens, disturbance introduced during sample acquisition and preparation), tight rocks can be particularly prone to solid–fluid interactions and thus more sensitive to the methods, procedures and techniques used to measure permeability. To address this problem, it is desirable to collect, for a single material, measurements obtained by different methods and pore fluids. For that purpose, a benchmarking exercise involving 24 laboratories was organized for measuring and modelling the permeability of a single low-permeability material, the Grimsel granodiorite. The objectives of the benchmark were: (i) to compare the results for a given method, (ii) to compare the results between different methods, (iii) to analyse the accuracy of each method, (iv) to study the influence of experimental conditions (especially the nature of pore fluid), (v) to discuss the relevance of indirect methods and models and finally (vi) to suggest good practice for low-permeability measurements. To complement the data set of permeability measurements presented in a companion paper, we focus here on (i) quantitative analysis of microstructures and pore size distribution, (ii) permeability modelling and (iii) complementary measurements of permeability anisotropy and poroelastic parameters. Broad ion beam—scanning electron microscopy, micro-computerized tomography, mercury injection capillary pressure (MICP) and nuclear magnetic resonance (NMR) methods were used to characterize the microstructures and provided the input parameters for permeability modelling. Several models were used: (i) basic statistical models, (ii) 3-D pore network and effective medium models, (iii) percolation model using MICP data and (iv) free-fluid model using NMR data. The models were generally successful in predicting the actual range of measured permeability. Statistical models overestimate the permeability because they do not adequately account for the heterogeneity of the crack network. Pore network and effective medium models provide additional constraints on crack parameters such as aspect ratio, aperture, density and connectivity. MICP and advanced microscopy techniques are very useful tools providing important input data for permeability estimation. Permeability measured—orthogonal to foliation is lower that—parallel to foliation. Combining the experimental and modelling results provide a unique and rich data set.

Geophysical Journal International↗
USGS2018

Rayleigh wave ellipticity measurement uncertainty across the IRIS/USGS and New China Digital Seismograph Networks

Long-period Rayleigh wave horizontal to vertical amplitude (H/V) ratios at a station provide information about local earth structure that is complementary to phase velocity. However, a number of studies have observed that significant scatter appears in these measurements making it difficult to use H/V ratio measurements to resolve earth structure. Some of the scatter in these measurements has been attributed to local geological structure while some has remained unaccounted for. Most Global Seismographic Network (GSN) stations contain two nearby high-quality broad-band seismometers (e.g. in the same vault, but on different piers or in different boreholes). For each broad-band sensor in the IRIS/USGS component of the GSN, we estimate H/V ratios of fundamental mode Rayleigh waves using M > 6.5 earthquakes from 2001 to 2018 (around 19 000 measurements). We compute these ratios at a number of discrete periods (25, 50, 75, 100 and 150 s) and find that for well-isolated Rayleigh waves (windows where the correlation coefficients between radial and the phase-shifted vertical components are greater than 0.9) significant scatter in H/V ratios occurs between colocated sensors (greater than 25 per cent at 100 s period). This suggests the scatter in H/V ratio measurements can be at least partially attributed to extremely local phenomena such as sensor emplacement in the vault. We also find that H/V ratios can vary as a function of event backazimuth, indicating that care must be taken when computing average ratios for a station, as a large number of events from a given region could bias H/V ratio measurements at a station.

Geophysical Journal International↗
USGS2019

Inversion of airborne EM data with an explicit choice of prior model

Inversion of airborne electromagnetic (AEM) data is an under-determined inverse problem, in that infinitely many resistivity models exist that will be able to explain the observed data, within measurement errors. Therefore, additional information or constraints must be taken into account to solve the inverse problem. In deterministic approaches, the goal is to locate one optimal model that can be obtained by using some form of smoothness constraints implied through a number of regularization choices. This model, however, will not necessarily represent realistic geological features. Probabilistic methods offer an alternative in which the solution is not one model, but a collection of models, whose variability represents the uncertainty. The probabilistic approach can also rely on implicit model assumptions, representing prior information (a type of regularization information) that may or may not be consistent with the actual available information. Here, we present an approach for AEM inversion in which the prior model is explicitly chosen by a user, preferably selected based on actual prior information available and then integrated with AEM data using a general Monte Carlo based sampling approach. This approach leads to a new workflow to AEM inversion in which geological prior information is independently and explicitly chosen before inversion is carried out. The main benefit of this approach is that each model obtained will, by construction, be consistent with prior (geological) information as well as geophysical data. Through examples based on synthetic and real AEM data, we will demonstrate the methodology, not least that the choice of prior information cannot be avoided: Either it is done explicitly, or it will be chosen implicitly by the choice of method used to invert the AEM data.

Geophysical Journal International↗
USGS2019

Variability in synthetic earthquake ground motions caused by source variability and errors in wave propagation models

Numerical simulations of earthquake ground motions are used both to anticipate the effects of hypothetical earthquakes by forward simulation and to infer the behaviour of the real earthquake source ruptures by the inversion of recorded ground motions. In either application it is necessary to assume some Earth structure that is necessarily inaccurate and to use a computational method that is also inaccurate for simulating the wavefield Green's functions. We refer to these two sources of error as ‘propagation inaccuracies’, which might be considered to be epistemic. We show that the variance of the Fourier spectrum of the synthetic earthquake seismograms caused by propagation inaccuracies is related to the spatial covariance on the rupture surface of errors in the computed Green's functions, which we estimate for the case of the 2009 L'Aquila, Italy, earthquake by comparing erroneous computed Green's functions with observed L'Aquila aftershock seismograms (empirical Green's functions). We further show that the variance of the synthetic seismograms caused by the rupture variability (aleatory uncertainty) is related to the spatial covariance on the rupture surface of aleatory variations in the rupture model, and we investigate the effect of correlated variations in Green's function errors and variations in rupture models. Thus, we completely characterize the variability of synthetic earthquake seismograms induced by errors in propagation and variability in the rupture behaviour. We calculate the spectra of the variance of the ground motions of the L'Aquila main shock caused by propagation inaccuracies for two specific broad-band stations, the AQU and the FIAM stations. These variances are distressingly large, being comparable or in some cases exceeding the data amplitudes, suggesting that the best-fitting L'Aquila rupture model significantly overfits the data and might be seriously in error. If these computed variances are typical, the accuracy of many other rupture models for past earthquakes may need to be reconsidered. The results of this work might be useful in seismic hazard estimation because the variability of the computed ground motion, caused both by propagation inaccuracies and variations in the rupture model, can be computed directly, not requiring laborious consideration of multiple Earth structures.

Geophysical Journal International↗
USGS2019

A combinatorial approach to determine earthquake magnitude distributions on a variable slip-rate fault

Combinatorial methods are used to determine the spatial distribution of earthquake magnitudes on a fault whose slip rate varies along strike. Input to the problem is a finite sample of earthquake magnitudes that span 5 kyr drawn from a truncated Pareto distribution. The primary constraints to the problem are maximum and minimum values around the target slip-rate function indicating where feasible solutions can occur. Two methods are used to determine the spatial distribution of earthquakes: integer programming and the greedy-sequential algorithm. For the integer-programming method, the binary decision vector includes all possible locations along the fault where each earthquake can occur. Once a set of solutions that satisfy the constraints is found, the cumulative slip misfit on the fault is globally minimized relative to the target slip-rate function. The greedy algorithm sequentially places earthquakes to locally optimize slip accumulation. As a case study, we calculate how earthquakes are distributed along the megathrust of the Nankai subduction zone, in which the slip rate varies significantly along strike. For both methods, the spatial distribution of magnitudes depends on slip rate, except for the largest magnitude earthquakes that span multiple sections of the fault. The greedy-sequential algorithm, previously applied to this fault (Parsons et al., 2012), tends to produce smoother spatial distributions and fewer lower magnitude earthquakes in the low slip-rate section of the fault compared to the integer-programming method. Differences in results from the two methods relate to how much emphasis is placed on minimizing the misfit to the target slip rate (integer programming) compared to finding a solution within the slip-rate constraints (greedy sequential). Specifics of the spatial distribution of magnitudes also depend on the shape of the target slip-rate function: i.e. stepped at the section boundaries versus a smooth function. This study isolates the effects of slip-rate variation along a single fault in determining the spatial distribution of earthquake magnitudes, helping to better interpret results from more complex, interconnected fault systems.

Geophysical Journal International↗
USGS2019

The effect of stress changes on time-dependent earthquake probabilities for the central Wasatch Fault Zone, Utah, USA

Static and quasi-static Coulomb stress changes produced by large earthquakes can modify the probability of occurrence of subsequent events on neighboring faults. This approach is based on physical (Coulomb stress changes) and statistical (probability calculations) models, which are influenced by the quality and quantity of data available in the study region. Here, we focus on the Wasatch Fault Zone (WFZ), a well-studied active normal fault system having abundant geologic and paleoseismological data. Paleoseismological trench investigations of the WFZ indicate that at least 24 large, surface-faulting earthquakes have ruptured the fault’s five central, 35–59-km long segments since ~7 ka. Our goal is to determine if the stress changes due to the youngest paleoevents have significantly modified the present-day probability of occurrence of large earthquakes on each of the segments. For each segment, we modeled the cumulative (coseismic + postseismic) Coulomb stress changes (∆CFScum) due to earthquakes younger than the most recent event on the segment in question and applied the resulting values to the time-dependent probability calculations. Results from the Coulomb stress modeling suggest that the Brigham City, Salt Lake City, and Provo segments have accumulated ∆CFScum larger than 10 bars, whereas the Weber segment has experienced a stress decrease of 5 bars, in the scenario of recent rupture of the Great Salt Lake fault to the west. Probability calculations predict high probability of occurrence for the Brigham City and Salt Lake City segments, due to their long elapsed times (>1-2 ka) when compared to the Weber, Provo, and Nephi segments (< 1 ka). The range of calculated coefficients of variation (CV) has a large influence on the final probabilities, mostly in the case of the Brigham City segment. Finally, when the Coulomb stress and the probability models are combined, our results indicate that the ∆CFScum resulting from earthquakes postdating the youngest events on each of the five segments significantly affects the probability calculations for three of the segments: Brigham City, Salt Lake City, and Provo. The probability of occurrence of a large earthquake in the next 50 years on these three segments may therefore be underestimated if a time-independent approach, or a time-dependent approach that does not consider ∆CFS, is adopted.

Utah↗
USGS2019

Improvements in seismic resolution and current limitations in the Global Seismographic Network

Station noise levels play a fundamental limitation in our ability to detect seismic signals. These noise levels are frequency-dependent and arise from a number of physically different drivers. At periods greater than 100 s, station noise levels are often limited by the self-noise of the instrument as well as the sensitivity of the instrument to non-seismic noise sources. Recently, station operators in the Global Seismographic Network (GSN) have deployed several Streckeisen STS-6A very broadband borehole seismometers. These sensors provide a potential replacement for the no-longer-produced Streckeisen STS-1 seismometer and the GeoTech KS-54000 borehole seismometer. Along with showing some of the initial observational improvements from installing modern very broadband seismometers at depth, we look at current limitations in the seismic resolution from Earth tide periods 100,000 s (0.01 mHz) to Nyquist at most GSN sites (0.02 s or 50 Hz). Finally, we show the potential for improved observations of continuously excited horizontal Earth hum as well as the splitting of very long-period torsional modes as a result of installing instruments at depth. Both of these observations make use of the low horizontal noise levels which are obtained by installing very broadband borehole seismometers at depth with noise levels similar to the Streckeisen STS-1.

Geophysical Journal International↗
USGS2010

Viscosity and finite strength of the mantle as determined from water and ice loads

Some recent examples of transient Earth loads (Lake Bonneville, Utah; Glacier Bay, Alaska; northeast Greenland) indicate that both the viscosity and finite strength of the mantle are lower than is commonly presumed. A time constant (1/ e ) of 4000 years is estimated for Lake Bonneville, and of 1000 years for northeast Greenland. A strain rate of 10 −14 is typical. These figures imply viscosities in a homogeneous half space ranging from 10 20 to 10 21 poises. An upper limit of finite strength is set by Lake Bonneville at a few times 10 6 dyn/cm 2 . If mountain ranges like the Sierra Nevada or Himalaya are regarded as dynamically supported rather than static systems, this low value is not incompatible with other geologic observations.

Alaska, Utah↗
USGS1969

Confidence limits for the precision parameter κ

Confidence limits are calculated for the precision parameter κ used in the analysis of palaeomagnetic data and for the angular standard deviation σ. A set of tables for 95 per cent and 99 per cent confidence limits is presented.

Geophysical Journal International↗
USGS1970

Latitude dependence of the angular dispersion of the geomagnetic field

Changes in the direction of the Earth's magnetic field at a given site are produced in part by wobble of the main geomagnetic dipole, in part by fluctuations in the intensity and direction of the non-dipole field, and in part by changes in the intensity of the main dipole field. These three processes combine to produce an angular variance that is strongly latitude dependent. A method is presented for isolating the contribution due to variation with latitude of the average intensity of the non-dipole field.

Geophysical Journal International↗
USGS1971

Palaeomagnetism and potassium-argon ages of volcanic rocks of Ngorongoro caldera, Tanzania

Natural remanent magnetization (NRM) has been measured in 21 lava flows and 1 tuff in the south-west wall of Ngorongoro caldera, Tanzania. The lowest three lavas are normally magnetized, the next two have intermediate directions, and the remainder are reversed; potassium-argon dating places the reversal at 2.45 My, corresponding to the Gauss Matuyama epoch boundary. Alternating-field demagnetization was effective in removing the ubiquitous lightning-produced secondary magnetizations, but partial thermal demagnetization was not. Curie temperatures of the lavas fell into two groups: 250°–350°C and 570°–585°C, probably corresponding to unoxidized and oxidized titanomagnetites respectively. NRM directions in the lavas occurred in several distinct and stratigraphically continuous groups, suggesting intermittent eruption of groups of lavas. No correspondence was found between NRM groups and other parameters, such as Curie temperatures. The dispersion of mean NRM directions is large in comparison with late Quaternary lavas from low latitudes in the Pacific basin, but probably is not large relative to the proportion of the geomagnetic field in Africa that is non-dipolar. The mean palaeomagnetic pole is offset 10.5° from the geographic pole, but because of the large dispersion, the difference is not considered significant.

Ngorongoro caldera↗
USGS1973

Critically refracted waves in a spherically symmetric radially heterogeneous Earth model

A theoretical analysis of acoustic waves refracted by a spherical boundary across which velocity and density increase abruptly and below which velocity and density may either increase or decrease continuously with depth is formulated in terms of waves generated at a harmonic point source and scattered by a radially heterogeneous spherical body. Through the application of an Earth-flattening transformation on the radial solution and the Watson transform on the sum over eigenfunctions, the solution to the spherical problem for high frequencies is expressed as an integral for the corresponding half-space problem in which the effect of boundary curvature maps into an effective positive velocity gradient with depth. The results of both analytical and numerical evaluation of this integral can be summarized as follows for body waves in the crust and upper mantle: (1) In the special case of a critical velocity gradient (a gradient equal and opposite to the effective curvature gradient), waves interacting with the boundary at the critical angle of incidence have the same form as the classical head wave for flat, homogeneous layers. (2) For gradients more negative than critical, the amplitude of waves incident at the critical angle decay more rapidly with distance than the classical head wave. (3) For gradients that are positive, null, and less negative than critical, the amplitude of waves near the critical angle decays less rapidly with distance than the classical head wave, and at sufficiently large distances, the refracted wave field can be adequately described in terms of ray-theoretical diving waves. At intermediate distances from the critical point, the spectral amplitude of the refracted wave is scalloped due to multiple diving wave interference.

Geophysical Journal of the Royal Astronomical Soci↗
USGS1974

Palaeomagnetism and magnetic–polarity zonation in some Oligocene volcanic rocks of the San Juan Mountains, south–western Colorado

Palaeomagnetic results have been obtained from thirty sites in intrusive and extrusive rocks of Oligocene age from the San Juan Mountains, south-western Colorado. All specimens from each site were subjected to af demagnetization, and the reliability of each site determined. Twenty-three sites gave reliable results. Because five sites from the thick intracaldera part of the La Jara Canyon Member of the Treasure Mountain Tuff appear to have become magnetized during the same small interval of geological time, their results were combined and their mean pole and direction used in the final calculations. The eighteen remaining reliable sites yielded an average Oligocene palaeomagnetic pole at 85·6° N and 298·0° E (δ p = 8°, δ m = 11°, k = 24·1). A stratigraphic sequence is given for the major San Juan ash-flow sheets and selected additional units with their corresponding magnetic polarities and mean K-Ar ages.

Colorado↗
USGS1980

Driving- stress waveform and the determination of rock internal friction by the stress-strain curve method.

Harmonic distortion in the stress-time function applied to rock specimens affects the measurement of rock internal friction in the seismic wave periods by the stress-strain hysteresis loop method. If neglected, the harmonic distortion can cause measurements of rock internal friction to be in error by 3O% in the linear range. The stress-time function therefore must be recorded and Fourier analysed for correct interpretation of the experimental data. Such a procedure would also yield a value for internal friction at the higher harmonic frequencies.-Author

Geophysical Journal of the Royal Astronomical Soci↗