Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Computational Physics”

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 19 recordsLinked to original sources
USGS1978

P-Wave to Rayleigh-wave conversion coefficients for wedge corners; model experiments

An analytic solution is not available for the diffraction of elastic waves by wedges; however, numerical solutions of finite-difference type are available for selected wedge angles. The P- to Rayleigh-wave conversion coefficients at wedge tips have been measured on two-dimensional seismic models for stress-free wedges with wedge angles, ??0, of 10, 30, 60, 90 and 120??. The conversion coefficients show two broad peaks and a minimum as a function of the angle between the wedge face and the direction of the incident P-wave. The minimum occurs for the P wave incident parallel to the wedge face and one maximum is near an incidence angle of 90?? to the wedge face. The amplitude of this maximum, relative to the other, decreases as the wedge angle increases. The asymmetry of the conversion coefficients, CPR(??; ??0), relative to parallel incidence (?? = 0) increases as the wedge angle increases. The locations of the maxima and the minimum as well as the asymmetry can be explained qualitatively. The conversion coefficients are measured with an accuracy of ??5% in those regions where there are no interfering waves. A comparison of the data for the 10?? wedge with the theoretical results for a half plane (0?? wedge) shows good correlation. ?? 1978.

Journal of Computational Physics↗
USGS1984

Computation of free oscillations of the earth

Although free oscillations of the Earth may be computed by many different methods, numerous practical considerations have led us to use a Rayleigh-Ritz formulation with piecewise cubic Hermite spline basis functions. By treating the resulting banded matrix equation as a generalized algebraic eigenvalue problem, we are able to achieve great accuracy and generality and a high degree of automation at a reasonable cost.

Journal of Computational Physics↗
USGS2008

Ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the Regional Ocean Modeling System

Systematic improvements in algorithmic design of regional ocean circulation models have led to significant enhancement in simulation ability across a wide range of space/time scales and marine system types. As an example, we briefly review the Regional Ocean Modeling System, a member of a general class of three-dimensional, free-surface, terrain-following numerical models. Noteworthy characteristics of the ROMS computational kernel include: consistent temporal averaging of the barotropic mode to guarantee both exact conservation and constancy preservation properties for tracers; redefined barotropic pressure-gradient terms to account for local variations in the density field; vertical interpolation performed using conservative parabolic splines; and higher-order, quasi-monotone advection algorithms. Examples of quantitative skill assessment are shown for a tidally driven estuary, an ice-covered high-latitude sea, a wind- and buoyancy-forced continental shelf, and a mid-latitude ocean basin. The combination of moderate-order spatial approximations, enhanced conservation properties, and quasi-monotone advection produces both more robust and accurate, and less diffusive, solutions than those produced in earlier terrain-following ocean models. Together with advanced methods of data assimilation and novel observing system technologies, these capabilities constitute the necessary ingredients for multi-purpose regional ocean prediction systems.

Journal of Computational Physics↗
USGS2020

Dimensional effects of inter-phase mass transfer on attenuation of structurally trapped gaseous carbon dioxide in shallow aquifers

Based on experimental evidence and using mathematical modeling, inter-phase mass transfer processes of CO 2 exsolving from and dissolving into water in heterogeneous porous media are investigated under two fundamentally different flow conditions: in a quasi one dimensional vertical column and in a two-dimensional tank with a lateral background water flow, both at laboratory scale. In both cases, the CO 2 dissolved in water under a given overpressure is injected for a certain period at the bottom of the tank, exsolves, and migrates upwards. A layer of fine sand is present in the tanks designed to mimic geological scenarios of accumulation and trapping of exsolved CO 2 in shallow aquifers. Then, clean water is injected and the accumulated CO 2 is dissolved back into the flowing water. The study aims to point out the differences in the mass transfer processes between the quasi-1D and 2D cases using a mathematical model of two-phase compositional flow in heterogeneous porous media calibrated to the experimental datasets, and expose strategies that should be explored in future research. Additionally, temperature variations observed during the 2D experiments allow for analysis of isothermal versus non-isothermal effects on the processes of multiphase CO 2 evolution. The mathematical model is discretized and solved using the mixed hybrid finite element method in 2D that allows for the simulation of both advection- and diffusion-dominated processes accurately.

Journal of Computational Physics↗
USGS2025

Total uncertainty quantification in inverse solutions with deep learning surrogate models

We propose an approximate Bayesian method for quantifying the total uncertainty in inverse partial differential equation (PDE) solutions obtained with machine learning surrogate models, including operator learning models. The proposed method accounts for uncertainty in the observations, PDE, and surrogate models. First, we use the surrogate model to formulate a minimization problem in the reduced space for the maximum a posteriori (MAP) inverse solution. Then, we randomize the MAP objective function and obtain samples of the posterior distribution by minimizing different realizations of the objective function. We test the proposed framework by comparing it with the iterative ensemble smoother and deep ensembling methods for a nonlinear diffusion equation with an unknown space-dependent diffusion coefficient. Among other applications, this equation describes the flow of groundwater in an unconfined aquifer. Depending on the training dataset and ensemble sizes, the proposed method provides similar or more descriptive posteriors of the parameters and states than the iterative ensemble smoother method. Deep ensembling underestimates uncertainty and provides less-informative posteriors than the other two methods. Our results show that, despite inherent uncertainty, surrogate models can be used for parameter and state estimation as an alternative to the inverse methods relying on (more accurate) numerical PDE solvers.

Journal of Computational Physics↗
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain 2. Numerical predictions and experimental tests

Numerical solutions of the equations describing flow of variably fluidized Coulomb mixtures predict key features of dry granular avalanches and water-saturated debris flows measured in physical experiments. These features include time-dependent speeds, depths, and widths of flows as well as the geometry of resulting deposits. Threedimensional (3-D) boundary surfaces strongly influence flow dynamics because transverse shearing and cross-stream momentum transport occur where topography obstructs or redirects motion. Consequent energy dissipation can cause local deceleration and deposition, even on steep slopes. Velocities of surge fronts and other discontinuities that develop as flows cross 3-D terrain are predicted accurately by using a Riemann solution algorithm. The algorithm employs a gravity wave speed that accounts for different intensities of lateral stress transfer in regions of extending and compressing flow and in regions with different degrees of fluidization. Field observations and experiments indicate that flows in which fluid plays a significant role typically have high-friction margins with weaker interiors partly fluidized by pore pressure. Interaction of the strong perimeter and weak interior produces relatively steep-sided, flat-topped deposits. To simulate these effects, we compute pore pressure distributions using an advection-diffusion model with enhanced diffusivity near flow margins. Although challenges remain in evaluating pore pressure distributions in diverse geophysical flows, Riemann solutions of the depthaveraged 3-D Coulomb mixture equations provide a powerful tool for interpreting and predicting flow behavior. They provide a means of modeling debris flows, rock avalanches, pyroclastic flows, and related phenomena without invoking and calibrating Theological parameters that have questionable physical significance.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Biomechanical factors contributing to self-organization in seagrass landscapes

Field observations have revealed that when water flow is consistently from one direction, seagrass shoots align in rows perpendicular to the primary axis of flow direction. In this study, live Zostera marina shoots were arranged either randomly or in rows perpendicular to the flow direction and tested in a seawater flume under unidirectional flow and waves to determine if shoot arrangement: a) influenced flow-induced force on individual shoots, b) differentially altered water flow through the canopy, and c) influenced light interception by the canopy. In addition, blade breaking strength was compared with flow-induced force to determine if changes in shoot arrangement might reduce the potential for damage to shoots. Under unidirectional flow, both current velocity in the canopy and force on shoots were significantly decreased when shoots were arranged in rows as compared to randomly. However, force on shoots was nearly constant with downstream distance, arising from the trade-off of shoot bending and in-canopy flow reduction. The coefficient of drag was higher for randomly-arranged shoots at low velocities (< 30 cm s- 1) but converged rapidly among the two shoot arrangements at higher velocities. Shoots arranged in rows tended to intercept slightly more light than those arranged randomly. Effects of shoot arrangement under waves were less clear, potentially because we did not achieve the proper plant size?row spacing ratio. At this point, we may only suggest that water motion, as opposed to light capture, is the dominant physical mechanism responsible for these shoot arrangements. Following a computation of the Environmental Stress Factor, we concluded that even photosynthetically active blades may be damaged or broken under frequently encountered storm conditions, irrespective of shoot arrangement. We hypothesize that when flow is generally from one direction, seagrass bed patterns over multiple scales of consideration may arise as a cumulative effect of individual shoot self-organization driven by reduced force and drag on the shoots and somewhat improved light capture.

Journal of Experimental Marine Biology and Ecology↗
USGS2012

Measurements of storm-generated bottom stresses on the continental shelf

Large values of bottom friction velocity, u *, and roughness length, z 0 , determined from burst-averaged speed data taken on the continental shelf in outer Norton Sound, Alaska, with the GEOPROBE tripod during a storm in September 1977 are correlated with extremely large values of near-bottom concentration of total suspended particulate matter (TSM). Combined wind-driven and tidal currents exceeding 30 cm/s at 1 m above the bottom and intense oscillatory bottom currents with maxima above 45 cm/s were associated with the largest measured values of TSM at 2 m above the sea floor. The values of u * and z 0 obtained from the ‘law of the wall’ velocity-depth relationship are diminished substantially throughout the storm period (average reduction of 44%) when the turbulence reducing effects of the vertical concentration gradient of TSM are considered. The form of the latter correction was adapted from Smith and McLean (1977 a ). Values of the mean u * computed from the theory of Grant and Madsen (1979), which predicts an enhanced shear stress due to nonlinear wave-current interactions, compare favorably with the u * values determined from the measured velocity profiles. The measured values of z 0 , however, are considerably larger than any of the estimates based on (1) the actual scales of the physical roughness elements; (2) the apparent roughness of Grant and Madsen (1979); or (3) the thickness of the bed-load layer as formulated by Smith and McLean (1977 a ).

Journal of Geophysical Research - Oceans↗
USGS2016

Comment on “The reduction of friction in long-runout landslides as an emergent phenomenon” by Brandon C. Johnson et al.

Results from a highly idealized, 2-D computational model indicate that dynamic normal-stress rarefactions might cause friction reduction in long-runout landslides, but the physical relevance of the idealized dynamics has not been confirmed by experimental tests. More importantly, the model results provide no evidence that refutes alternative hypotheses about friction reduction mechanisms. One alternative hypothesis, which is strongly supported by field evidence, experimental data, and the predictions of a well-constrained computational model, involves development of high pore fluid pressures in deforming landslide material or overridden bed material. However, no scientific basis exists for concluding that a universal mechanism is responsible for friction reduction in all long-runout landslides.

Journal of Geophysical Research F: Earth Surface↗
USGS2010

Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration

Shallow rainfall-induced landslides commonly occur under conditions of transient infiltration into initially unsaturated soils. In an effort to predict the timing and location of such landslides, we developed a model of the infiltration process using a two-layer system that consists of an unsaturated zone above a saturated zone and implemented this model in a geographic information system (GIS) framework. The model links analytical solutions for transient, unsaturated, vertical infiltration above the water table to pressure-diffusion solutions for pressure changes below the water table. The solutions are coupled through a transient water table that rises as water accumulates at the base of the unsaturated zone. This scheme, though limited to simplified soil-water characteristics and moist initial conditions, greatly improves computational efficiency over numerical models in spatially distributed modeling applications. Pore pressures computed by these coupled models are subsequently used in one-dimensional slope-stability computations to estimate the timing and locations of slope failures. Applied over a digital landscape near Seattle, Washington, for an hourly rainfall history known to trigger shallow landslides, the model computes a factor of safety for each grid cell at any time during a rainstorm. The unsaturated layer attenuates and delays the rainfall-induced pore-pressure response of the model at depth, consistent with observations at an instrumented hillside near Edmonds, Washington. This attenuation results in realistic estimates of timing for the onset of slope instability (7 h earlier than observed landslides, on average). By considering the spatial distribution of physical properties, the model predicts the primary source areas of landslides.

Journal of Geophysical Research F: Earth Surface↗
USGS1988

Range indices of geomagnetic activity

The simplest index of geomagnetic activity is the range in nT from maximum to minimum value of the field in a given time interval. The hourly range R was recommended by IAGA for use at observatories at latitudes greater than 65??, but was superceded by AE. The most used geomagnetic index K is based on the range of activity in a 3 h interval corrected for the regular daily variation. In order to take advantage of real time data processing, now available at many observatories, it is proposed to introduce a 1 h range index and also a 3 h range index. Both will be computed hourly, i.e. each will have a series of 24 per day, the 3 h values overlapping. The new data will be available as the range (R) of activity in nT and also as a logarithmic index (I) of the range. The exponent relating index to range in nT is based closely on the scale used for computing K values. The new ranges and range indices are available, from June 1987, to users in real time and can be accessed by telephone connection or computer network. Their first year of production is regarded as a trial period during which their value to the scientific and commercial communities will be assessed, together with their potential as indicators of regional and global disturbances' and in which trials will be conducted into ways of eliminating excessive bias at quiet times due to the rate of change of the daily variation field. ?? 1988.

Journal of Atmospheric and Terrestrial Physics↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 2. Experimental tests

Scaling considerations indicate that miniature experiments can be used to test models of granular avalanches in which the effects of intergranular fluid and cohesion are negligible. To test predictions of a granular avalanche model described in a companion paper, we performed bench top experiments involving avalanches of dry sand across irregular basal topography that mimicked the complexity of natural terrain. The experiments employed a novel method of laser-assisted cartography to map the three-dimensional morphology of rapidly moving avalanches, thereby providing high-resolution data for comparison with model output. Model input consisted of two material properties (angles of internal and basal Coulomb friction of the sand), which were measured in independent tests, and of initial and boundary conditions that characterized the geometry of the experimental apparatus. Experimental results demonstrate that the model accurately predicts not only the gross behavior but also many details of avalanche motion, from initiation to deposition. We attribute this accuracy to a mathematical and computational formulation that conserves mass and momentum in three-dimensional physical space and satisfies the Coulomb equation in three-dimensional stress space. Our results support the hypothesis that a Coulomb proportionality between shear and normal stresses applies in moderately rapid granular flows and that complicated constitutive postulates are unnecessary if momentum conservation is strictly enforced in continuum avalanche models. Furthermore, predictions of our Coulomb continuum model contrast with those of a Coulomb point mass model, illustrating the importance of multidimensional modeling and model testing.

Journal of Geophysical Research Earth Surface↗
USGS2008

Physical property changes in hydrate-bearing sediment due to depressurization and subsequent repressurization

Physical property measurements of sediment cores containing natural gas hydrate are typically performed on material exposed, at least briefly, to non-in situ conditions during recovery. To examine the effects of a brief excursion from the gas-hydrate stability field, as can occur when pressure cores are transferred to pressurized storage vessels, we measured physical properties on laboratory-formed sand packs containing methane hydrate and methane pore gas. After depressurizing samples to atmospheric pressure, we repressurized them into the methane-hydrate stability field and remeasured their physical properties. Thermal conductivity, shear strength, acoustic compressional and shear wave amplitudes, and speeds of the original and depressurized/repressurized samples are compared. X– ray computed tomography images track how the gas-hydrate distribution changes in the hydrate-cemented sands owing to the depressurizaton/repressurization process. Because depressurization-induced property changes can be substantial and are not easily predicted, particularly in water-saturated, hydrate-bearing sediment, maintaining pressure and temperature conditions throughout the core recovery and measurement process is critical for using laboratory measurements to estimate in situ properties.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Partly cloudy with a chance of lava flows: Forecasting volcanic eruptions in the 21st century

A primary goal of volcanology is forecasting hazardous eruptive activity. Despite much progress over the last century, however, volcanoes still erupt with no detected precursors, lives and livelihoods are lost to eruptive activity, and forecasting the onsets of eruptions remains fraught with uncertainty. Long‐term forecasts are generally derived from the geological and historical records, from which recurrence intervals and styles of activity can be inferred, while shorter‐term forecasts are derived from patterns in monitoring data. Information from geology and monitoring data can be evaluated and combined using statistical analysis, expert elicitation, and conceptual and or mathematical models. Integrative frameworks, such as event trees, combine this diversity of information to produce probabilistic forecasts that can inform the style and scale of the societal response to a potential future eruption. Several developments show promise to revolutionize the utility and accuracy of these forecasts. These include growth in the quantity and quality of multidisciplinary monitoring data, coupled with increases in computing power; machine learning algorithms, which will allow far better utilization of this growing volume of data; and new physiochemical volcano models and data assimilation algorithms, which take advantage of a wide range of monitoring data and realistic physics to better predict the evolution of a given physical state. Although eruption forecasts may never be as generally reliable as weather forecasts, and great caution must be exercised when attempting to predict highly complex volcanic behavior, these and other innovations—particularly when combined in integrative, fully probabilistic forecasting frameworks—should help volcanologists to better issue warnings of volcanic activity on societally relevant time frames.

Hawaii↗
USGS1981

A five-collector system for the simultaneous measurement of argon isotope ratios in a static mass spectrometer

A system is described that utilizes five separate Faraday-cup collector assemblies, aligned along the focal plane of a mass spectrometer, to collect simultaneous argon ion beams at masses 36–40. Each collector has its own electrometer amplifier and analog-to-digital measuring channel, the outputs of which are processed by a minicomputer that also controls the mass spectrometer. The mass spectrometer utilizes a 90° sector magnetic analyzer with a radius of 23 cm, in which some degree of z -direction focussing is provided for all the ion beams by the fringe field of the magnet. Simultaneous measurement of the ion beams helps to eliminate mass-spectrometer memory as a significant source of measurement error during an analysis. Isotope ratios stabilize between 7 and 9 s after sample admission into the spectrometer, and thereafter changes in the measured ratios are linear, typically to within ±0.02%. Thus the multi-collector arrangement permits very short extrapolation times for computation of initial ratios, and also provides the advantages of simultaneous measurement of the ion currents in that errors due to variations in ion beam intensity are minimized. A complete analysis takes less than 10 min, so that sample throughput can be greatly enhanced. In this instrument, the factor limiting analytical precision now lies in short-term apparent variations in the interchannel calibration factors.

International Journal of Mass Spectrometry and Ion↗
USGS2011

Assessing first-order emulator inference for physical parameters in nonlinear mechanistic models

We present an approach for estimating physical parameters in nonlinear models that relies on an approximation to the mechanistic model itself for computational efficiency. The proposed methodology is validated and applied in two different modeling scenarios: (a) Simulation and (b) lower trophic level ocean ecosystem model. The approach we develop relies on the ability to predict right singular vectors (resulting from a decomposition of computer model experimental output) based on the computer model input and an experimental set of parameters. Critically, we model the right singular vectors in terms of the model parameters via a nonlinear statistical model. Specifically, we focus our attention on first-order models of these right singular vectors rather than the second-order (covariance) structure.

Journal of Agricultural, Biological, and Environme↗
USGS2021

Retention and dimensional changes of evergreen brush piles within a flood control reservoir

Brush piles (i.e., trees and large woody debris) are often installed in reservoirs to supplement fish habitat. The retention and dimensional change of brush piles after installation is important information that can be used to maximize the effectiveness of this management action. We evaluated the retention and dimensional change of 70 eastern red cedar Juniperus virginiana and bald cypress Taxodium distichum brush piles in an embayment of a drawdown reservoir up to four annual cycles of submergence and exposure. We used satellite imagery to supplement our onsite measurements of retention. We also examined spatial patterns of brush pile retention and dimensional change. Brush piles were lost at 10% per year, and their volume was lost at 14% per year. We compared our rates of brush pile retention and dimensional change with those from a holdout data set of 50 brush piles. Estimates between data sets did not differ statistically. Spatial patterns of retention and dimensional change coincided with morphological features in our study area, suggesting that retention and dimensional change is influenced by variable physical forces (e.g., wave action and flow) at installation locations. Our estimates of brush pile retention and dimensional change can be used to generally sustain desirable brush densities. For example, to maintain a fixed total volume of brush in our study embayment, roughly 23% of the total brush volume installed would need to be replaced annually. Similar research in reservoirs managed for other purposes is needed, as length and cycle of inundation could lead to variable rates of retention and dimensional change. Additionally, advancements into computer-assisted detection and volume estimation could reduce the time and effort needed to monitor brush piles.

Mississippi↗
USGS2017

A multimodal wave spectrum-based approach for statistical downscaling of local wave climate

Characterization of wave climate by bulk wave parameters is insufficient for many coastal studies, including those focused on assessing coastal hazards and long-term wave climate influences on coastal evolution. This issue is particularly relevant for studies using statistical downscaling of atmospheric fields to local wave conditions, which are often multimodal in large ocean basins (e.g. the Pacific). Swell may be generated in vastly different wave generation regions, yielding complex wave spectra that are inadequately represented by a single set of bulk wave parameters. Furthermore, the relationship between atmospheric systems and local wave conditions is complicated by variations in arrival time of wave groups from different parts of the basin. Here, we address these two challenges by improving upon the spatiotemporal definition of the atmospheric predictor used in statistical downscaling of local wave climate. The improved methodology separates the local wave spectrum into “wave families,” defined by spectral peaks and discrete generation regions, and relates atmospheric conditions in distant regions of the ocean basin to local wave conditions by incorporating travel times computed from effective energy flux across the ocean basin. When applied to locations with multimodal wave spectra, including Southern California and Trujillo, Peru, the new methodology improves the ability of the statistical model to project significant wave height, peak period, and direction for each wave family, retaining more information from the full wave spectrum. This work is the base of statistical downscaling by weather types, which has recently been applied to coastal flooding and morphodynamic applications.

Journal of Physical Oceanography↗