Search USGS⌕ Search

SEARCH · Search USGS

Results for “Inverse Problems”

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 109 records · Page 6Linked to original sources
USGS2024

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications↗
USGS1998

Singularity and Nonnormality in the Classification of Compositional Data

Geologists may want to classify compositional data and express the classification as a map. Regionalized classification is a tool that can be used for this purpose, but it incorporates discriminant analysis, which requires the computation and inversion of a covariance matrix. Covariance matrices of compositional data always will be singular (noninvertible) because of the unit-sum constraint. Fortunately, discriminant analyses can be calculated using a pseudo-inverse of the singular covariance matrix; this is done automatically by some statistical packages such as SAS. Granulometric data from the Darss Sill region of the Baltic Sea is used to explore how the pseudo-inversion procedure influences discriminant analysis results, comparing the algorithm used by SAS to the more conventional Moore-Penrose algorithm. Logratio transforms have been recommended to overcome problems associated with analysis of compositional data, including singularity. A regionalized classification of the Darss Sill data after logratio transformation is different only slightly from one based on raw granulometric data, suggesting that closure problems do not influence severely regionalized classification of compositional data.

Mathematical Geology↗
USGS2013

Effects of error covariance structure on estimation of model averaging weights and predictive performance

When conducting model averaging for assessing groundwater conceptual model uncertainty, the averaging weights are often evaluated using model selection criteria such as AIC, AICc, BIC, and KIC (Akaike Information Criterion, Corrected Akaike Information Criterion, Bayesian Information Criterion, and Kashyap Information Criterion, respectively). However, this method often leads to an unrealistic situation in which the best model receives overwhelmingly large averaging weight (close to 100%), which cannot be justified by available data and knowledge. It was found in this study that this problem was caused by using the covariance matrix, C E , of measurement errors for estimating the negative log likelihood function common to all the model selection criteria. This problem can be resolved by using the covariance matrix, C ek , of total errors (including model errors and measurement errors) to account for the correlation between the total errors. An iterative two-stage method was developed in the context of maximum likelihood inverse modeling to iteratively infer the unknown C ek from the residuals during model calibration. The inferred C ek was then used in the evaluation of model selection criteria and model averaging weights. While this method was limited to serial data using time series techniques in this study, it can be extended to spatial data using geostatistical techniques. The method was first evaluated in a synthetic study and then applied to an experimental study, in which alternative surface complexation models were developed to simulate column experiments of uranium reactive transport. It was found that the total errors of the alternative models were temporally correlated due to the model errors. The iterative two-stage method using C ek resolved the problem that the best model receives 100% model averaging weight, and the resulting model averaging weights were supported by the calibration results and physical understanding of the alternative models. Using C ek obtained from the iterative two-stage method also improved predictive performance of the individual models and model averaging in both synthetic and experimental studies.

Water Resources Research↗
USGS2010

Time-dependent seismic tomography

Of methods for measuring temporal changes in seismic-wave speeds in the Earth, seismic tomography is among those that offer the highest spatial resolution. 3-D tomographic methods are commonly applied in this context by inverting seismic wave arrival time data sets from different epochs independently and assuming that differences in the derived structures represent real temporal variations. This assumption is dangerous because the results of independent inversions would differ even if the structure in the Earth did not change, due to observational errors and differences in the seismic ray distributions. The latter effect may be especially severe when data sets include earthquake swarms or aftershock sequences, and may produce the appearance of correlation between structural changes and seismicity when the wave speeds are actually temporally invariant. A better approach, which makes it possible to assess what changes are truly required by the data, is to invert multiple data sets simultaneously, minimizing the difference between models for different epochs as well as the rms arrival-time residuals. This problem leads, in the case of two epochs, to a system of normal equations whose order is twice as great as for a single epoch. The direct solution of this system would require twice as much memory and four times as much computational effort as would independent inversions. We present an algorithm, tomo4d, that takes advantage of the structure and sparseness of the system to obtain the solution with essentially no more effort than independent inversions require. No claim to original US government works Journal compilation ?? 2010 RAS.

Geophysical Journal International↗
USGS1979

Monitoring massive fracture growth at 2-km depths using surface tiltmeter arrays

Tilt due to massive hydraulic fractures induced in sedimentary rocks at depths of up to 2.2 km have been recorded by surface tiltmeters. Injection of fluid volumes up to 4 · 10 5 liters and masses of propping agent up to 5 · 10 5 kg is designed to produce fractures approximately 1 km long, 50–100 m high and about 1 cm wide. The surface tilt data adequately fit a dislocation model of a tensional fault in a half-space. Theoretical and observational results indicate that maximum tilt occurs at a distance off the strike of the fracture equivalent to 0.4 of the depth to the fracture. Azimuth and extent of the fracture deduced from the geometry of the tilt field agree with other kinds of geophysical measurements. Detailed correlation of the tilt signatures with pumping parameters (pressure, rate, volume, mass) have provided details on asymmetry in geometry and growth rate. Whereas amplitude variations in tilt vary inversely with the square of the depth, changes in flow rate or pressure gradient can produce a cubic change in width. These studies offer a large-scale experimental approach to the study of problems involving fracturing, mass transport, and dilatancy processes.

Tectonophysics↗
USGS2005

Constraints on fault slip rates of the southern California plate boundary from GPS velocity and stress inversions

We use Global Positioning System (GPS) velocities and stress orientations inferred from seismicity to invert for the distribution of slip on faults in the southern California plate-boundary region. Of particular interest is how long-term slip rates are partitioned between the Indio segment of the San Andreas fault (SAF), the San Jacinto fault (SJF) and the San Bernardino segment of the SAE We use two new sets of constraints to address this problem. The first is geodetic velocities from the Southern California Earthquake Center's (SCEC) Crustal Motion Map (version 3 by Shen et al.), which includes significantly more data than previous models. The second is a regional model of stress-field orientations at seismogenic depths, as determined from earthquake focal mechanisms. While GPS data have been used in similar studies before, this is the first application of stress-field observations to this problem. We construct a simplified model of the southern California fault system, and estimate the interseismic surface velocities using a backslip approach with purely elastic strain accumulation, following Meade et al. In addition, we model the stress orientations at seismogenic depths, assuming that crustal stress results from the loading of active faults. The geodetically derived stressing rates are found to be aligned with the stress orientations from seismicity. We therefore proceed to invert simultaneously GPS and stress observations for slip rates of the faults in our network. We find that the regional patterns of crustal deformation as imaged by both data sets can be explained by our model, and that joint inversions lead to better constrained slip rates. In our preferred model, the SJF accommodates ???15 mm yr-1 and the Indio segment of the SAF ???23 mm yr-1 of right-lateral motion, accompanied by a low slip rate on the San Bernardino segment of the SAF 'Anomalous' fault segments such as around the 1992 Mw = 7.3 Landers surface rupture can be detected. There, observed stresses deviate strongly from the long-term loading as predicted by our simple model. Evaluation of model misfits together with information from palaeoseismology may provide further insights into the time dependence of strain accumulation along the San Andreas system. ?? 2004 RAS.

Geophysical Journal International↗
USGS1982

Earthquake location in island arcs

A comprehensive data set of selected teleseismic P-wave arrivals and local-network P- and S-wave arrivals from large earthquakes occurring at all depths within a small section of the central Aleutians is used to examine the general problem of earthquake location in island arcs. Reference hypocenters for this special data set are determined for shallow earthquakes from local-network data and for deep earthquakes from combined local and teleseismic data by joint inversion for structure and location. The high-velocity lithospheric slab beneath the central Aleutians may displace hypocenters that are located using spherically symmetric Earth models; the amount of displacement depends on the position of the earthquakes with respect to the slab and on whether local or teleseismic data are used to locate the earthquakes. Hypocenters for trench and intermediate-depth events appear to be minimally biased by the effects of slab structure on rays to teleseismic stations. However, locations of intermediate-depth events based on only local data are systematically displaced southwards, the magnitude of the displacement being proportional to depth. Shallow-focus events along the main thrust zone, although well located using only local-network data, are severely shifted northwards and deeper, with displacements as large as 50 km, by slab effects on teleseismic travel times. Hypocenters determined by a method that utilizes seismic ray tracing through a three-dimensional velocity model of the subduction zone, derived by thermal modeling, are compared to results obtained by the method of joint hypocenter determination (JHD) that formally assumes a laterally homogeneous velocity model over the source region and treats all raypath anomalies as constant station corrections to the travel-time curve. The ray-tracing method has the theoretical advantage that it accounts for variations in travel-time anomalies within a group of events distributed over a sizable region of a dipping, high-velocity lithospheric slab. In application, JHD has the practical advantage that it does not require the specification of a theoretical velocity model for the slab. Considering earthquakes within a 260 km long by 60 km wide section of the Aleutian main thrust zone, our results suggest that the theoretical velocity structure of the slab is presently not sufficiently well known that accurate locations can be obtained independently of locally recorded data. Using a locally recorded earthquake as a calibration event, JHD gave excellent results over the entire section of the main thrust zone here studied, without showing a strong effect that might be attributed to spatially varying source-station anomalies. We also calibrated the ray-tracing method using locally recorded data and obtained results generally similar to those obtained by JHD. ?? 1982.

Physics of the Earth and Planetary Interiors↗
USGS2011

Mechanical analysis of fault slip data: Implications for paleostress analysis

Stress inversions are a useful and popular tool for structural geologist and seismologist alike. These methods were first introduced by Wallace (1951) and Bott (1959) and subsequent studies continue to be based on their assumptions: the remote stress tensor is spatially uniform for the rock mass containing the faults and temporally constant over the history of faulting in that region, and the slip on each fault surface has the same direction and sense as the maximum shear stress resolved on that surface from the remote stress tensor. Furthermore, successful implementation requires that slip accumulates on faults of diverse orientation. Many studies employ these methods on isolated faults or on fault systems with limited ranges of orientations, which can lead to erroneous results. We propose a new method that incorporates the effects of mechanical interaction of the entire fault or fault system, and solves the complete mechanical problem rather than employing empirical relationships between slip and stress or strain (or strain rate). The method requires knowledge of the fault geometry and information on at least one slip vector component along portions of the known fault geometry. For example, if throw is known, the strike-slip component can be solved for. We test the method using a single synthetic fault with anisotropic roughness similar to that measured at fault outcrops. While the orientation of remote stress may be determined precisely, the lack of diverse fault orientations introduces a systematic error in the remote stress ratio. We further test the effect of diversity of fault orientations and find that Wallace–Bott type inversions do not perform as well for limited ranges of orientations when compared to the proposed method. Finally, we use published data from the 1999 Chi-Chi, Taiwan, earthquake, and find that the method using surface data only, and surface data with subsurface focal mechanisms, produce similar results. The resulting stress orientations are in good agreement with results from Wallace–Bott inversions. Furthermore, the slip distribution is in general agreement with kinematic slip inversions using coseismic surface deformation. Stress inversion methods using fault slip data can thus be improved upon, significantly in some cases, by solving a mechanical boundary value problem that takes into account the geometry of faults or fault systems. As a bonus, the solution provides the stress, strain, and displacement fields throughout the region and the slip distributions on the faults.

Journal of Structural Geology↗
USGS2009

Obtaining parsimonious hydraulic conductivity fields using head and transport observations: A Bayesian geostatistical parameter estimation approach

Flow path delineation is a valuable tool for interpreting the subsurface hydrogeochemical environment. Different types of data, such as groundwater flow and transport, inform different aspects of hydrogeologic parameter values (hydraulic conductivity in this case) which, in turn, determine flow paths. This work combines flow and transport information to estimate a unified set of hydrogeologic parameters using the Bayesian geostatistical inverse approach. Parameter flexibility is allowed by using a highly parameterized approach with the level of complexity informed by the data. Despite the effort to adhere to the ideal of minimal a priori structure imposed on the problem, extreme contrasts in parameters can result in the need to censor correlation across hydrostratigraphic bounding surfaces. These partitions segregate parameters into facies associations. With an iterative approach in which partitions are based on inspection of initial estimates, flow path interpretation is progressively refined through the inclusion of more types of data. Head observations, stable oxygen isotopes ( 18 O/ 16 O ratios), and tritium are all used to progressively refine flow path delineation on an isthmus between two lakes in the Trout Lake watershed, northern Wisconsin, United States. Despite allowing significant parameter freedom by estimating many distributed parameter values, a smooth field is obtained.

Water Resources Research↗
USGS2013

Assessing grain-size correspondence between flow and deposits of controlled floods in the Colorado River, USA

Flood-deposited sediment has been used to decipher environmental parameters such as variability in watershed sediment supply, paleoflood hydrology, and channel morphology. It is not well known, however, how accurately the deposits reflect sedimentary processes within the flow, and hence what sampling intensity is needed to decipher records of recent or long-past conditions. We examine these problems using deposits from dam-regulated floods in the Colorado River corridor through Marble Canyon–Grand Canyon, Arizona, U.S.A., in which steady-peaked floods represent a simple end-member case. For these simple floods, most deposits show inverse grading that reflects coarsening suspended sediment (a result of fine-sediment-supply limitation), but there is enough eddy-scale variability that some profiles show normal grading that did not reflect grain-size evolution in the flow as a whole. To infer systemwide grain-size evolution in modern or ancient depositional systems requires sampling enough deposit profiles that the standard error of the mean of grain-size-change measurements becomes small relative to the magnitude of observed changes. For simple, steady-peaked floods, 5–10 profiles or fewer may suffice to characterize grain-size trends robustly, but many more samples may be needed from deposits with greater variability in their grain-size evolution.

Journal of Sedimentary Research↗
USGS2007

Physically based estimation of soil water retention from textural data: General framework, new models, and streamlined existing models

Numerous models are in widespread use for the estimation of soil water retention from more easily measured textural data. Improved models are needed for better prediction and wider applicability. We developed a basic framework from which new and existing models can be derived to facilitate improvements. Starting from the assumption that every particle has a characteristic dimension R associated uniquely with a matric pressure ψ and that the form of the ψ– R relation is the defining characteristic of each model, this framework leads to particular models by specification of geometric relationships between pores and particles. Typical assumptions are that particles are spheres, pores are cylinders with volume equal to the associated particle volume times the void ratio, and that the capillary inverse proportionality between radius and matric pressure is valid. Examples include fixed-pore-shape and fixed-pore-length models. We also developed alternative versions of the model of Arya and Paris that eliminate its interval-size dependence and other problems. The alternative models are calculable by direct application of algebraic formulas rather than manipulation of data tables and intermediate results, and they easily combine with other models (e.g., incorporating structural effects) that are formulated on a continuous basis. Additionally, we developed a family of models based on the same pore geometry as the widely used unsaturated hydraulic conductivity model of Mualem. Predictions of measurements for different suitable media show that some of the models provide consistently good results and can be chosen based on ease of calculations and other factors.

Vadose Zone Journal↗
USGS2018

Gas hydrate quantification using full-waveform inversion of sparse ocean-bottom seismic data: A case study from Green Canyon Block 955, Gulf of Mexico

We present a case study of gas hydrate quantification using dense short-offset multichannel seismic (MCS) and sparse long-offset ocean-bottom-seismometer (OBS) data in lease block Green Canyon 955 (GC955), Gulf of Mexico (GOM), where the presence of gas hydrate was interpreted using logging while drilling (LWD) data acquired by the GOM Gas Hydrate Joint Industry Project Leg II expedition. We use frequency-domain full-waveform inversion (FWI) of seven OBS gathers to invert for a P-wave velocity model of an approximately 7 km long MCS profile connecting two LWD sites, GC955-H and GC955-Q. We build the starting model for FWI using traveltime inversion (TI) of the MCS and OBS data. In addition, we use the TI model for depth migrating the MCS stack. At the LWD sites, we constrain the hydrate saturation (Sgh) using sonic and resistivity logs. Unfortunately, as is typical of seismic quantification problems, the FWI model resolution is not sufficient to extrapolate the LWD-based Sgh. Therefore, we apply Backus averaging to the sonic log, at 60 m wavelength, bringing it within approximately 8% of the FWI model and make the assumption that averaging the sonic log is same as redistributing the gas hydrate within the Backus wavelength. In this manner, instead of Sgh, the FWI model is able to estimate the total gas hydrate volume. In the end, we use the FWI model and the migrated stack to constrain the locations and bulk volumes of free gas and gas hydrate. Our results demonstrate that with careful processing, reasonable estimates on locations and bulk volumes of submarine gas hydrate accumulations can be achieved even with sparse seismic data that are not adequate for amplitude-based assessments.

Geophysics↗
USGS2023

Using DC resistivity ring array surveys to resolve conductive structures around tunnels or mine-workings

In underground environments, conventional direct current (DC) resistivity surveys with a single linear array of electrodes produce fundamentally non-unique inversions. These non-uniqueness and model resolution issues stem from limitations placed on the location of transmitters (TXs) and receivers (RXs) by the geometry of existing tunnels and boreholes. Poor excitation and/or sampling of the region of interest (ROI) can create artifacts and reduce the resolution of the recovered model. To address these problems we propose the use of an ensemble of ring arrays, which are created by placing one or more electrodes in each face (sidewalls, floor, and ceiling) of the tunnel to form a ring of electrodes at each along-tunnel location. Using a series of increasingly complex synthetic models, we assess the benefits of ring arrays and show that they can be used to better constrain the location and shape of anomalous bodies around the tunnel. Although ring arrays significantly improve the resolution of the recovered model, the size of the comprehensive ring array survey increases rapidly with the number of electrodes used. To balance model resolution and survey size, we developed a physics-based survey design methodology. In this methodology, TXs are selected based upon secondary charge accumulations on a test block that is moved throughout the ROI. Although this survey design methodology does not produce a strictly optimal survey, it balances model resolution and survey size in a practical and computationally efficient manner. Since the ring array more accurately estimates the around-tunnel location of targets and ensures that targets on all sides of the tunnel are detected, it is ideally suited to tunnel-based environments. Our results show that only about 6 % "> 6% of the possible TXs and 0.5 % "> 0.5% of the RXs in the comprehensive ring array survey are needed to retain the improvements in resolution. Therefore, economical ring array surveys can be designed for both reconnaissance and target characterization. Following the inversion of the reconnaissance dataset, additional rings can be added to reduce the inter-ring spacing or off-tunnel boreholes can be added to the region around identified anomalies to increase resolution as required.

Journal of Applied Geophysics↗
USGS1995

Analysis of nutrient and ancillary water-quality data for surface and ground water of the Willamette Basin, Oregon, 1980-90

An analysis of historical water-quality data for surface and ground water collected in the Willamette and Sandy River Basins during the 1980-90 water years was performed. For surface water, most data were concentrated at sites on the main stem Willamette River or near the mouths of major tributaries. All seasons were represented. Data for nitrogen and phosphorus species were readily available, but simultaneously collected discharge measurements frequently were not. Seven primary sites were used for a quantitative analysis of nutrient data in surface water. At six of the seven primary sites, median concentrations of nitrite plus nitrate and total phosphorus were less than 0.5 milligrams per liter (mg/L) as nitrogen (N) and 0.1 mg/L as phosphorus (P), respectively. These concen- trations were lower than national median concen- trations for basins with similar land uses. At the Pudding River site, which received significant point-source and agricultural nonpoint-source inputs, median values for nitrite plus nitrate and total phosphorus were 1.4 mg/L as N and 0.2 mg/L as P, respectively. Nitrite-plus-nitrate concentrations at this site were seasonally dependent, with the highest concentrations occurring during winter. Forested sites had significantly lower water temperatures and nutrient concentrations than urban or agricultural sites. Evidence of diel variations in dissolved oxygen concentrations and pH values suggest that, at some sites, low dissolved oxygen concentrations may be a problem during early morning hours in the summer. Historical nutrient data for ground water were limited primarily to nitrate and nitrite-plus-nitrate determinations for wells completed in the basin-fill and alluvial aquifer. Nitrate and nitrite-plus-nitrate concentrations in this aquifer showed a weak inverse relationship with depth.

Water-Resources Investigations Report↗
USGS2002

Groundwater contamination downstream of a contaminant penetration site. II. Horizontal penetration of the contaminant plume

Part I of this study (Rubin, H.; Buddemeier, R.W. Groundwater Contamination Downstream of a Contaminant Penetration Site Part 1: Extension-Expansion of the Contaminant Plume. J. of Environmental Science and Health Part A (in press).) addressed cases, in which a comparatively thin contaminated region represented by boundary layers (BLs) developed within the freshwater aquifer close to contaminant penetration site. However, at some distance downstream from the penetration site, the top of the contaminant plume reaches the top or bottom of the aquifer. This is the location of the "attachment point," which comprises the entrance cross section of the domain evaluated by the present part of the study. It is shown that downstream from the entrance cross section, a set of two BLs develop in the aquifer, termed inner and outer BLs. It is assumed that the evaluated domain, in which the contaminant distribution gradually becomes uniform, can be divided into two sections, designated: (a) the restructuring section, and (b) the establishment section. In the restructuring section, the vertical concentration gradient leads to expansion of the inner BL at the expense of the outer BL, and there is almost no transfer of contaminant mass between the two layers. In the establishment section, each of the BLs occupies half of the aquifer thickness, and the vertical concentration gradient leads to transfer of contaminant mass from the inner to the outer BL. By use of BL approximations, changes of salinity distribution in the aquifer are calculated and evaluated. The establishment section ends at the uniformity point, downstream from which the contaminant concentration profile is practically uniform. The length of the restructuring section, as well as that of the establishment section, is approximately proportional to the aquifer thickness squared, and is inversely proportional to the transverse dispersivity. The study provides a convenient set of definitions and terminology that are helpful in visualizing the gradual development of uniform contaminant concentration distribution in an aquifer subject to contaminant plume penetration. The method developed in this study can be applied to a variety of problems associated with groundwater quality, such as initial evaluation of field data, design of field data collection, the identification of appropriate boundary conditions for numerical models, selection of appropriate numerical modeling approaches, interpretation and evaluation of field monitoring results, etc.

Journal of Environmental Science and Health - Part↗
USGS2015

Point spread functions for earthquake source imaging: An interpretation based on seismic interferometry

Recently, various methods have been proposed and applied for earthquake source imaging, and theoretical relationships among the methods have been studied. In this study, we make a follow-up theoretical study to better understand the meanings of earthquake source imaging. For imaging problems, the point spread function (PSF) is used to describe the degree of blurring and degradation in an obtained image of a target object as a response of an imaging system. In this study, we formulate PSFs for earthquake source imaging. By calculating the PSFs, we find that waveform source inversion methods remove the effect of the PSF and are free from artifacts. However, the other source imaging methods are affected by the PSF and suffer from the effect of blurring and degradation due to the restricted distribution of receivers. Consequently, careful treatment of the effect is necessary when using the source imaging methods other than waveform inversions. Moreover, the PSF for source imaging is found to have a link with seismic interferometry with the help of the source-receiver reciprocity of Green’s functions. In particular, the PSF can be related to Green’s function for cases in which receivers are distributed so as to completely surround the sources. Furthermore, the PSF acts as a low-pass filter. Given these considerations, the PSF is quite useful for understanding the physical meaning of earthquake source imaging.

Geophysical Journal International↗
USGS2010

Estimating 14C groundwater ages in a methanogenic aquifer

This paper addresses the problem of 14 C age dating of groundwaters in a confined regional aquifer affected by methanogenesis. Increasing CH 4 concentrations along the groundwater flow system and 13 C and 14 C isotopic data for dissolved inorganic carbon, dissolved organic carbon, and CH 4 clearly show the effect of methanogenesis on groundwater chemistry. Inverse reaction path modeling using NETPATH indicates the predominant geochemical reactions controlling the chemical evolution of groundwater in the aquifer are incongruent dissolution of dolomite, ion exchange, methanogenesis, and oxidation of sedimentary organic matter. Modeling of groundwater 14 C ages using NETPATH indicates that a significant part of groundwater in the Alliston aquifer is less than 13,000 years old; however, older groundwater in the range of 15,000–23,000 years is also present in the aquifer. This paper demonstrates that 14C ages calculated using NETPATH, incorporating the effects of methanogenesis on the carbon pools, provide reasonable groundwater ages that were not possible by other isotopic methods.

Water Resources Research↗
USGS2016

The Earthquake‐Source Inversion Validation (SIV) Project

Finite‐fault earthquake source inversions infer the (time‐dependent) displacement on the rupture surface from geophysical data. The resulting earthquake source models document the complexity of the rupture process. However, multiple source models for the same earthquake, obtained by different research teams, often exhibit remarkable dissimilarities. To address the uncertainties in earthquake‐source inversion methods and to understand strengths and weaknesses of the various approaches used, the Source Inversion Validation (SIV) project conducts a set of forward‐modeling exercises and inversion benchmarks. In this article, we describe the SIV strategy, the initial benchmarks, and current SIV results. Furthermore, we apply statistical tools for quantitative waveform comparison and for investigating source‐model (dis)similarities that enable us to rank the solutions, and to identify particularly promising source inversion approaches. All SIV exercises (with related data and descriptions) and statistical comparison tools are available via an online collaboration platform, and we encourage source modelers to use the SIV benchmarks for developing and testing new methods. We envision that the SIV efforts will lead to new developments for tackling the earthquake‐source imaging problem.

Seismological Research Letters↗