USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysics and Engineering”

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 73 records · Page 4Linked to original sources
USGS2011

Wave-current interaction in Willapa Bay

This paper describes the importance of wave-current interaction in an inlet-estuary system. The three-dimensional, fully coupled, Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system was applied in Willapa Bay (Washington State) from 22 to 29 October 1998 that included a large storm event. To represent the interaction between waves and currents, the vortex-force method was used. Model results were compared with water elevations, currents, and wave measurements obtained by the U.S. Army Corp of Engineers. In general, a good agreement between field data and computed results was achieved, although some discrepancies were also observed in regard to wave peak directions in the most upstream station. Several numerical experiments that considered different forcing terms were run in order to identify the effects of each wind, tide, and wave-current interaction process. Comparison of the horizontal momentum balances results identified that wave-breaking-induced acceleration is one of the leading terms in the inlet area. The enhancement of the apparent bed roughness caused by waves also affected the values and distribution of the bottom shear stress. The pressure gradient showed significant changes with respect to the pure tidal case. During storm conditions the momentum balance in the inlet shares the characteristics of tidal-dominated and wave-dominated surf zone environments. The changes in the momentum balance caused by waves were manifested both in water level and current variations. The most relevant effect on hydrodynamics was a wave-induced setup in the inner part of the estuary.

Washington↗
USGS2005

Accurate elevation and normal moveout corrections of seismic reflection data on rugged topography

The application of the seismic reflection method is often limited in areas of complex terrain. The problem is the incorrect correction of time shifts caused by topography. To apply normal moveout (NMO) correction to reflection data correctly, static corrections are necessary to be applied in advance for the compensation of the time distortions of topography and the time delays from near-surface weathered layers. For environment and engineering investigation, weathered layers are our targets, so that the static correction mainly serves the adjustment of time shifts due to an undulating surface. In practice, seismic reflected raypaths are assumed to be almost vertical through the near-surface layers because they have much lower velocities than layers below. This assumption is acceptable in most cases since it results in little residual error for small elevation changes and small offsets in reflection events. Although static algorithms based on choosing a floating datum related to common midpoint gathers or residual surface-consistent functions are available and effective, errors caused by the assumption of vertical raypaths often generate pseudo-indications of structures. This paper presents the comparison of applying corrections based on the vertical raypaths and bias (non-vertical) raypaths. It also provides an approach of combining elevation and NMO corrections. The advantages of the approach are demonstrated by synthetic and real-world examples of multi-coverage seismic reflection surveys on rough topography. ?? The Royal Society of New Zealand 2005.

New Zealand Journal of Geology and Geophysics↗
USGS2014

Methane hydrates in nature - Current knowledge and challenges

Recognizing the importance of methane hydrate research and the need for a coordinated effort, the United States Congress enacted the Methane Hydrate Research and Development Act of 2000. At the same time, the Ministry of International Trade and Industry in Japan launched a research program to develop plans for a methane hydrate exploratory drilling project in the Nankai Trough. India, China, the Republic of Korea, and other nations also have established large methane hydrate research and development programs. Government-funded scientific research drilling expeditions and production test studies have provided a wealth of information on the occurrence of methane hydrates in nature. Numerous studies have shown that the amount of gas stored as methane hydrates in the world may exceed the volume of known organic carbon sources. However, methane hydrates represent both a scientific and technical challenge, and much remains to be learned about their characteristics and occurrence in nature. Methane hydrate research in recent years has mostly focused on: (1) documenting the geologic parameters that control the occurrence and stability of methane hydrates in nature, (2) assessing the volume of natural gas stored within various methane hydrate accumulations, (3) analyzing the production response and characteristics of methane hydrates, (4) identifying and predicting natural and induced environmental and climate impacts of natural methane hydrates, (5) analyzing the methane hydrate role as a geohazard, (6) establishing the means to detect and characterize methane hydrate accumulations using geologic and geophysical data, and (7) establishing the thermodynamic phase equilibrium properties of methane hydrates as a function of temperature, pressure, and gas composition. The U.S. Department of Energy (DOE) and the Consortium for Ocean Leadership (COL) combined their efforts in 2012 to assess the contributions that scientific drilling has made and could continue to make to advance our understanding of methane hydrates in nature. COL assembled a Methane Hydrate Project Science Team with members from academia, industry, and government. This Science Team worked with COL and DOE to develop and host the Methane Hydrate Community Workshop, which surveyed a substantial cross section of the methane hydrate research community for input on the most important research developments in our understanding of methane hydrates in nature and their potential role as an energy resource, a geohazard, and/or as an agent of global climate change. Our understanding of how methane hydrates occur in nature is still growing and evolving, and it is known with certainty that field, laboratory, and modeling studies have contributed greatly to our understanding of hydrates in nature and will continue to be a critical source of the information needed to advance our understanding of methane hydrates.

Journal of Chemical and Engineering Data↗
USGS2021

Identifying geomorphic process domains in the synthetic landscapes of West Virginia, USA

Human activities such as mining, agriculture, and urbanization have resulted in severe, large‐scale alteration to landform organization and associated geomorphic processes. The mountaintop mining (MTM) region of West Virginia, USA has experienced dramatic topographic alteration, by removing steep slopes and introducing plateau‐like areas at ridgelines and benches on valley fills. The resulting engineered landforms create synthetic landscapes, disconnected from previous geomorphic processes. Invoking the process domain concept, we compare differences in slope‐area relations, cumulative area distributions (CADs), elevation, slope, upslope accumulated area, and a slope*area product before and after mining to adjacent unmined sub‐catchments in five study basins. Differences in the slope‐area relation include a 42% slope reduction in low drainage areas, corresponding to hillslopes, unchanneled valleys, and debris flow dominated channels, which may fall below thresholds required for debris flow processes. The curved slope‐area relation that represents valley incision by debris flows is replaced by slope‐area relations that resemble basins where gullying and the stream power law dominate. Extremely high chemical weathering of unconsolidated valley fills materials may facilitate process domain shifts from debris flows to gullying and fluvial erosion. The characteristic power law scaling break in CADs that represents the headward limit of the channelized network is subdued in post‐mined sites and may reflect headward channelized network extension in mined basins. Slope‐area relations and CADs present a unique topographic signature of MTM activity, potentially providing an analytical approach to assess impacts on underlying geomorphic processes for other synthetic landscapes such as cities or large‐scale agricultural production.

West Virginia↗
USGS2001

Geological characterization of remote field sites using visible and infrared spectroscopy: Results from the 1999 Marsokhod field test

Upcoming Mars Surveyor lander missions will include extensive spectroscopic capabilities designed to improve interpretations of the mineralogy and geology of landing sites on Mars. The 1999 Marsokhod Field Experiment (MFE) was a Mars rover simulation designed in part to investigate the utility of visible/near-infrared and thermal infrared field spectrometers to contribute to the remote geological exploration of a Mars analog field site in the California Mojave Desert. The experiment simultaneously investigated the abilities of an off-site science team to effectively analyze and acquire useful imaging and spectroscopic data and to communicate efficiently with rover engineers and an on-site field team to provide meaningful input to rover operations and traverse planning. Experiences gained during the MFE regarding effective communication between different mission operation teams will be useful to upcoming Mars mission teams. Field spectra acquired during the MFE mission exhibited features interpreted at the time as indicative of carbonates (both dolomitic and calcitic), mafic rocks and associated weathering products, and silicic rocks with desert varnish-like coatings. The visible/near-infrared spectra also suggested the presence of organic compounds, including chlorophyll in one rock. Postmission laboratory petrologic and spectral analyses of returned samples confirmed that all rocks identified as carbonates using field measurements alone were calc-silicates and that chlorophyll associated with endolithic organisms was present in the one rock for which it was predicted. Rocks classified from field spectra as silicics and weathered mafics were recognized in the laboratory as metamorphosed monzonites and diorite schists. This discrepancy was likely due to rock coatings sampled by the field spectrometers compared to fresh rock interiors analyzed petrographically, in addition to somewhat different surfaces analyzed by laboratory thermal spectroscopy compared to field spectra.

California↗
USGS2008

Coupled semivariogram uncertainty of hydrogeological and geophysical data on capture zone uncertainty analysis

This study investigates capture zone uncertainty that relates to the coupled semivariogram uncertainty of hydrogeological and geophysical data. Semivariogram uncertainty is represented by the uncertainty in structural parameters (range, sill, and nugget). We used the beta distribution function to derive the prior distributions of structural parameters. The probability distributions of structural parameters were further updated through the Bayesian approach with the Gaussian likelihood functions. Cokriging of noncollocated pumping test data and electrical resistivity data was conducted to better estimate hydraulic conductivity through autosemivariograms and pseudo-cross-semivariogram. Sensitivities of capture zone variability with respect to the spatial variability of hydraulic conductivity, porosity and aquifer thickness were analyzed using ANOVA. The proposed methodology was applied to the analysis of capture zone uncertainty at the Chicot aquifer in Southwestern Louisiana, where a regional groundwater flow model was developed. MODFLOW-MODPATH was adopted to delineate the capture zone. The ANOVA results showed that both capture zone area and compactness were sensitive to hydraulic conductivity variation. We concluded that the capture zone uncertainty due to the semivariogram uncertainty is much higher than that due to the kriging uncertainty for given semivariograms. In other words, the sole use of conditional variances of kriging may greatly underestimate the flow response uncertainty. Semivariogram uncertainty should also be taken into account in the uncertainty analysis. ?? 2008 ASCE.

Journal of Hydrologic Engineering↗
USGS2022

Measurements of Geologic Characteristics and Geophysical Properties of Sediments From the New England Mud Patch

The characterization of physical, geological, and geophysical properties of sediments within the New England Mud Patch (NEMP) was undertaken to provide a physical basis for acoustic inversions associated with the SeaBed Characterization EXperiment 2017 (SBCEX17). Using a suite of 89 sediment cores (piston/trigger, gravity [acoustic], and vibracore), a comprehensive database of laboratory-based sediment analyses, geophysical core logs, and the results of seismic reflection profiling, we formulate a three-layer lithostratigraphic model of the area within and immediately adjacent to the SBCEX17 focus area, referred to as the seabed experiment area (SEA). The uppermost lithostratigraphic unit, Unit 1, is relatively homogenous clayey- to sandy silt, with consistent downcore textural, mineralogical, and physical property attributes. Unit 2 is a variable-thickness transitional layer between Unit 1 and Unit 3, whose properties reflect a decrease in proximal erosion and transition to a lower energy depositional environment. Unit 3 is clean quartz sand containing abundant shells and shell fragments that was regionally deposited during Holocene sea-level rise. 210 Pb and 14 C radiocarbon geochronologies spanning the past 13 000 years are used to facilitate intercore comparison across the SEA. Analytical results and laboratory methods used in the derivation of those results are described in detail, serving as a reference for ongoing and future investigation of the SEA and entire NEMP. Although the derived lithostratigraphic model of the SEA is in good agreement with past evaluations of the regional sedimentology, comparisons of the lithostratigraphic and seismostratigraphic models highlight several significant incompatibilities that remain to be satisfactorily explained.

IEEE Journal of Ocean Engineering↗
USGS2017

Quantifying the heterogeneity of the tectonic stress field using borehole data

The heterogeneity of the tectonic stress field is a fundamental property which influences earthquake dynamics and subsurface engineering. Self-similar scaling of stress heterogeneities is frequently assumed to explain characteristics of earthquakes such as the magnitude-frequency relation. However, observational evidence for such scaling of the stress field heterogeneity is scarce. We analyze the local stress orientations using image logs of two closely spaced boreholes in the Coso Geothermal Field with sub-vertical and deviated trajectories, respectively, each spanning about 2 km in depth. Both the mean and the standard deviation of stress orientation indicators (borehole breakouts, drilling-induced fractures and petal-centerline fractures) determined from each borehole agree to the limit of the resolution of our method although measurements at specific depths may not. We find that the standard deviation in these boreholes strongly depends on the interval length analyzed, generally increasing up to a wellbore log length of about 600 m and constant for longer intervals. We find the same behavior in global data from the World Stress Map. This suggests that the standard deviation of stress indicators characterizes the heterogeneity of the tectonic stress field rather than the quality of the stress measurement. A large standard deviation of a stress measurement might be an expression of strong crustal heterogeneity rather than of an unreliable stress determination. Robust characterization of stress heterogeneity requires logs that sample stress indicators along a representative sample volume of at least 1 km.

Journal of Geophysical Research B: Solid Earth↗
USGS1997

Selection of the Mars Pathfinder landing site

The Mars Pathfinder spacecraft will land on a depositional fan near the mouth of the catastrophic outflow channel, Ares Vallis (19.5°N, 32.8°W). This site offers the prospect of analyzing a variety of rock types from the ancient cratered highlands, intermediate-age ridged plains, and reworked channel deposits. Analyses of these rocks by Pathfinder instruments will enable first-order scientific questions to be addressed, such as differentiation of the crust, the development of weathering products, and the nature of the early environment, as well as their subsequent evolution on Mars. Constraints imposed by (1) spacecraft and rover designs (which are robust), (2) entry, descent, and landing, (3) scientific potential at various sites, and (4) safety were important considerations in site selection. Engineering constraints require a 70 km by 200 km smooth, flat (low slopes) area located between 10° and 20°N that is below 0 km elevation, with average radar reflectivity, little dust, and moderate rock abundance. Three regions on Mars are between 10° and 20°N and below 0 km elevation: Chryse, Amazonis, and Isidis-Elysium. Science considerations favor sites at the mouths of outflow channels (grab bag sites offer an assay of rock types on Mars), highland sites (early crustal differentiation and climate), and sites covered with dark (unoxidized) material. Sites are considered safe if they are clearly below 0 km elevation, appear acceptably free of hazards in high-resolution (<50 m/pixel) Viking orbiter images and have acceptable reflectivity and roughness at radar wavelengths, thermal inertia, rock abundance, red to violet ratio, and albedo. Recent 3.5-cm wavelength radar observations were used to verify elevation, reflectivity, and roughness within the landing ellipses. Three sites meet all of these criteria: Ares Vallis, Tritonis Lacus, and Isidis. Although Isidis appears to be safer than Tritonis and Ares, the greater scientific potential at Ares Vallis resulted in its selection. Comparisons of the Grand Coulee (channel) and the depositional Ephrata Fan of the Channeled Scabland in eastern Washington, with Ares Vallis and its depositional fan also suggest the Ares Vallis landing site is safe and scientifically interesting.

Journal of Geophysical Research Planets↗
USGS2012

Strain accumulation across the Denali fault in the Delta River canyon, Alaska

Deformation along the Denali fault in the Delta River canyon was determined from geodetic surveys in 1941/1942, 1970, 1975, and 1979. The data were best for the 1975–1979 interval; in that period the average strain accumulation was essentially pure right lateral shear at a rate of 0.6 ± 0.1 μrad/a (a is years) (engineering shear) across a vertical plane striking N87°E. The plane of maximum shear is rotated about 30° counterclockwise from the local strike of the Denali fault but closely coincides with the strike of a major linear segment of the fault that begins 50 km farther west. The deformation between 1941/1942 and 1970 is consistent with a similar rate of strain accumulation if one removes the coseismic strain step contributed by the 1964 Alaska earthquake. The 1970–1975 deformation is poorly defined owing to uncertainties in the 1970 survey, but the strain accumulation during that period is certainly much less than during the 1975–1979 interval. The 1975–1979 strain accumulation is interpreted by means of a dislocation model which suggests that the Denali fault in the vicinity of the Delta River Canyon behaves as a leaky transform fault (i.e., a source of spreading as well as lateral slip). The block south of the Denali fault appears to be moving westward (parallel to the strike of the Denali fault west of the 147th meridian) at the rate of about 20 mm/a relative to the North American plate. Because the linear segment of the Denali fault east of the 147th meridian (along which the geodetic network is located) strikes N63°W, accommodation of the westward motion of the southern block requires some spreading within the fault zone as well as right lateral slip on the fault.

Journal of Geophysical Research Solid Earth↗
USGS2021

Movement of sediment through a burned landscape: Sediment volume observations and model comparisons in the San Gabriel Mountains, California, USA

Post-wildfire changes to hydrologic and geomorphic systems can lead to widespread sediment redistribution. Understanding how sediment moves through a watershed is crucial for assessing hazards, developing debris flow inundation models, engineering sediment retention solutions, and quantifying the role that disturbances play in landscape evolution. In this study, we used terrestrial and airborne lidar to measure sediment redistribution in the 2016 Fish Fire, in the San Gabriel Mountains in southern California, USA. The lidar areas are in two adjacent watersheds, at spatial scales of 900 m 2 to 4 km 2 , respectively. Terrestrial lidar data were acquired prior to rainfall, and two subsequent surveys show erosional change after rainstorms. Two airborne lidar flights occurred (1) 7 months before, and (2) 14 months after the fire ignition, capturing the erosional effects after rainfall. We found hillslope erosion dominated the overall sediment budget in the first rainy season after wildfire. Only 7% of the total erosion came from the active channel bed and channel banks, and the remaining 93% of eroded sediment was derived from hillslopes. Within the channelized portion of the watershed erosion/deposition could be generally described with topographic metrics used in a stream power equation. Observed sediment volumes were compared with four empirical models and one process-based model. We found that the best predictions of sediment volume were obtained from an empirical model developed in the same physiographic region. Moreover, this study showed that post-wildfire erosion rates in the San Gabriel Mountains attain the same magnitude as millennial time scale bedrock erosion rates.

California↗
USGS2019

Earthquake stress drop and Arias Intensity

We determine earthquake stress drops directly from the Arias intensity database of NGA-West2. Arias intensity (Arias, 1970) is an engineering measure proportional to the integral of the absolute value of acceleration squared, over the significant duration of the signal. As such, it is closely related to root-mean-square acceleration, and can readily be connected to earthquake stress drop (Hanks and McGuire, 1981). Arias intensity records out to 100 km yield stable stress drops for moderate-to-large magnitude earthquakes, M6.5+; for smaller events ~M4.5 – 6.5, only closer-in records yield stable results. For the 116 events considered, stress drops are about 35% larger for Class 1 mainshocks than for traditional on-fault Class 2 aftershocks, and smaller for those aftershocks close to the main fault plane. Aftershock stress drops show large variability, however, implying that on average they re-rupture weakened patches, but can also rupture intact rock or high-stress asperities. We observe an increase of stress drop with earthquake depth similar to that of other studies but do not find any significant faulting mechanism dependence. The variability of the Arias intensity-based stress drop is lower than that of eGf-based stress drops from Baltay et al. (2010, 2011), and nearly on par with variability seen in ground-motion prediction equations. The Arias intensity stress drop is a novel and promising method to estimate stress drop without the need for path and site corrections, and yields further insight into the connection between source physics and ground-motion.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

On extracting sediment transport information from measurements of luminescence in river sediment

Accurately quantifying sediment transport rates in rivers remains an important goal for geomorphologists, hydraulic engineers, and environmental scientists. However, current techniques for measuring long-time scale (102–106 years) transport rates are laborious, and formulae to predict transport are notoriously inaccurate. Here we attempt to estimate sediment transport rates by using luminescence, a property of common sedimentary minerals that is used by the geoscience community for geochronology. This method is advantageous because of the ease of measurement on ubiquitous quartz and feldspar sand. We develop a model from first principles by using conservation of energy and sediment mass to explain the downstream pattern of luminescence in river channel sediment. We show that the model can accurately reproduce the luminescence observed in previously published field measurements from two rivers with very different sediment transport styles. The model demonstrates that the downstream pattern of river sand luminescence should show exponential-like decay in the headwaters which asymptotes to a constant value with further downstream distance. The parameters from the model can then be used to estimate the time-averaged virtual velocity, characteristic transport lengthscale, storage time scale, and floodplain exchange rate of fine sand-sized sediment in a fluvial system. The sediment transport values predicted from the luminescence method show a broader range than those reported in the literature, but the results are nonetheless encouraging and suggest that luminescence demonstrates potential as a sediment transport indicator. However, caution is warranted when applying the model as the complex nature of sediment transport can sometimes invalidate underlying simplifications.

Journal of Geophysical Research F: Earth Surface↗
USGS1996

Inelastic models of lithospheric stress - I. Theory and application to outer-rise plate deformation

Outer-rise stress distributions determined in the manner that mechanical engineers evaluate inelastic stress distributions within conventional materials are contrasted with those predicted using simple elastic-plate models that are frequently encountered in studies of outer-rise seismicity. This comparison indicates that the latter are inherently inappropriate for studies of intraplate earthquakes, which are a direct manifestation of lithospheric inelasticity. We demonstrate that the common practice of truncating elastically superimposed stress profiles so that they are not permitted to exceed laboratory-based estimates of lithospheric yield strength will result in an accurate characterization of lithospheric stress only under relatively restrictive circumstances. In contrast to elastic-plate models, which predict that lithospheric stress distributions depend exclusively upon the current load, inelastic plate models predict that stress distributions are also significantly influenced by the plate-loading history, and, in many cases, this influence is the dominant factor in determining the style of potential seismicity (e.g. thrust versus normal faulting). Numerous ‘intuitive’ interpretations of outer-rise earthquakes have been founded upon the implicit assumption that a unique relationship exists between a specified combination of plate curvature and in-plane force, and the resulting lithospheric stress distribution. We demonstrate that the profound influence of deformation history often invalidates such interpretations. Finally, we examine the reliability of ‘yield envelope’ representations of lithospheric strength that are constructed on the basis of empirically determined frictional sliding relationships and silicate plastic-flow laws. Although representations of this nature underestimate the strength of some major interplate faults, such as the San Andreas, they appear to represent a reliable characterization of the strength of intraplate oceanic lithosphere.

Geophysical Journal International↗
USGS2003

Selection of the Mars Exploration Rover landing sites

The selection of Meridiani Planum and Gusev crater as the Mars Exploration Rover landing sites took over 2 years, involved broad participation of the science community via four open workshops, and narrowed an initial ∼155 potential sites (80–300 × 30 km) to four finalists based on science and safety. Engineering constraints important to the selection included (1) latitude (10°N–15°S) for maximum solar power, (2) elevation (less than −1.3 km) for sufficient atmosphere to slow the lander, (3) low horizontal winds, shear, and turbulence in the last few kilometers to minimize horizontal velocity, (4) low 10‐m‐scale slopes to reduce airbag spin‐up and bounce, (5) moderate rock abundance to reduce abrasion or strokeout of the airbags, and (6) a radar‐reflective, load‐bearing, and trafficable surface safe for landing and roving that is not dominated by fine‐grained dust. The evaluation of sites utilized existing as well as targeted orbital information acquired from the Mars Global Surveyor and Mars Odyssey. Three of the final four landing sites show strong evidence for surface processes involving water and appear capable of addressing the science objectives of the missions, which are to determine the aqueous, climatic, and geologic history of sites on Mars where conditions may have been favorable to the preservation of evidence of possible prebiotic or biotic processes. The evaluation of science criteria placed Meridiani and Gusev as the highest‐priority sites. The evaluation of the three most critical safety criteria (10‐m‐scale slopes, rocks, and winds) and landing simulation results indicated that Meridiani and Elysium Planitia are the safest sites, followed by Gusev and Isidis Planitia.

Journal of Geophysical Research E: Planets↗
USGS2012

Evidence for a new geomagnetic reversal from lava flows in Idaho: Discussion of short polarity reversals in the Brunhes and late Matuyama polarity chrons

K-Ar ages and paleomagnetic data for basalt samples from a new core hole (site E) at the Idaho National Engineering Laboratory (INEL) indicate that the age of the reversed polarity event recorded in Snake River Plain lavas is older than 465±50 ka (1000 years before present) reported previously by Champion et al. (1981). Nine basalt flows, eight with normal polarity and one with reversed polarity, were recognized in the site E core hole. The flows above and below the reversed flow have ages of 491±80 ka and 580±93 ka, respectively. The inclination of the paleomagnetic field direction of the reversed flow at site E agrees with the inclination of reversed flows elsewhere at INEL which have an age of 565±14 ka. These reversed flows were previously thought to be correlative with the Emperor event. We suggest that this polarity event is an older event which we name the Big Lost Reversed Polarity Subchronozone and Subchron. A review of data documenting short reversal records from volcanic and sedimentary rocks shows that there is evidence for eight polarity subchrons in the Brunhes and two besides the Jaramillo in the late Matuyama. These 10 short subchrons begin to indicate the many short events that Cox (1968) hypothesized must exist if polarity interval lengths have a Poisson distribution. These events are true subchrons, not excursions, and may or may not have low associated paleointensities, although low field strengths might explain why the reversal process aborts. The mean sustained polarity interval length since late Matuyama Chron time is 90,000 years. The similarity of this number with the 10 5 -year period of the Earth's orbital eccentricity suggests anew that linkage between geomagnetic, paleoclimatic, and possible underlying Earth orbital parameters should be evaluated.

Journal of Geophysical Research Solid Earth↗
USGS2003

An empirical approach to inversion of an unconventional helicopter electromagnetic dataset

A helicopter electromagnetic (HEM) survey acquired at the U.S. Idaho National Engineering and Environmental Laboratory (INEEL) used a modification of a traditional mining airborne method flown at low levels for detailed characterization of shallow waste sites. The low sensor height, used to increase resolution, invalidates standard assumptions used in processing HEM data. Although the survey design strategy was sound, traditional interpretation techniques, routinely used in industry, proved ineffective. Processed data and apparent resistivity maps were severely distorted, and hence unusable, due to low flight height effects, high magnetic permeability of the basalt host, and the conductive, three-dimensional nature of the waste site targets.To accommodate these interpretation challenges, we modified a one-dimensional inversion routine to include a linear term in the objective function that allows for the magnetic and three-dimensional electromagnetic responses in the in-phase data. Although somewhat ad hoc, the use of this term in the inverse routine, referred to as the shift factor, was successful in defining the waste sites and reducing noise due to the low flight height and magnetic characteristics of the host rock. Many inversion scenarios were applied to the data and careful analysis was necessary to determine the parameters appropriate for interpretation, hence the approach was empirical. Data from three areas were processed with this scheme to highlight different interpretational aspects of the method. Wastes sites were delineated with the shift terms in two of the areas, allowing for separation of the anthropomorphic targets from the natural one-dimensional host. In the third area, the estimated resistivity and the shift factor were used for geological mapping. The high magnetic content of the native soil enabled the mapping of disturbed soil with the shift term. Published by Elsevier Science B.V.

Journal of Applied Geophysics↗
USGS2011

A distribution-based parameterization for improved tomographic imaging of solute plumes

Difference geophysical tomography (e.g. radar, resistivity and seismic) is used increasingly for imaging fluid flow and mass transport associated with natural and engineered hydrologic phenomena, including tracer experiments, in situ remediation and aquifer storage and recovery. Tomographic data are collected over time, inverted and differenced against a background image to produce ‘snapshots’ revealing changes to the system; these snapshots readily provide qualitative information on the location and morphology of plumes of injected tracer, remedial amendment or stored water. In principle, geometric moments (i.e. total mass, centres of mass, spread, etc.) calculated from difference tomograms can provide further quantitative insight into the rates of advection, dispersion and mass transfer; however, recent work has shown that moments calculated from tomograms are commonly biased, as they are strongly affected by the subjective choice of regularization criteria. Conventional approaches to regularization (Tikhonov) and parametrization (image pixels) result in tomograms which are subject to artefacts such as smearing or pixel estimates taking on the sign opposite to that expected for the plume under study. Here, we demonstrate a novel parametrization for imaging plumes associated with hydrologic phenomena. Capitalizing on the mathematical analogy between moment-based descriptors of plumes and the moment-based parameters of probability distributions, we design an inverse problem that (1) is overdetermined and computationally efficient because the image is described by only a few parameters, (2) produces tomograms consistent with expected plume behaviour (e.g. changes of one sign relative to the background image), (3) yields parameter estimates that are readily interpreted for plume morphology and offer direct insight into hydrologic processes and (4) requires comparatively few data to achieve reasonable model estimates. We demonstrate the approach in a series of numerical examples based on straight-ray difference-attenuation radar monitoring of the transport of an ionic tracer, and show that the methodology outlined here is particularly effective when limited data are available.

Geophysical Journal International↗