Search USGS⌕ Search

SEARCH · Search USGS

Results for “Applied Acoustics”

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

Rapid earthquake magnitude classification via P-wave strains from borehole strainmeters and Distributed Acoustic Sensing

Distributed Acoustic Sensing (DAS) offers a promising approach for earthquake early warning (EEW) in settings where seismic networks are costly to maintain. By repurposing fiber-optic cables as dense strainmeter arrays, DAS enables real-time earthquake detection wherever those fibers are accessible. However, poor azimuthal coverage and challenges in estimating magnitude from strain measurements remain key hurdles in applying for earthquake monitoring. Here, we develop a machine learning method to distinguish large (M≥5.4) earthquakes from smaller ones within the first 4 seconds of a strain waveform after a P-wave arrival without determining location. Using ensemble decision tree models trained on borehole strainmeter data (3.5≤M≤7.1) and tested on onshore DAS waveforms (including the 2024 M7 Offshore Cape Mendocino earthquake), we find that low-frequency (0.2–0.5 Hz) continuous wavelet transform coefficients are the strongest predictors of magnitude, in addition to strain amplitude. Both DAS and borehole strainmeters effectively capture long-period strain signals, making these findings valuable for EEW systems. Our method shows high precision compared to the real-time EEW system, ShakeAlert®, supporting the position that DAS is a viable technology for earthquake monitoring and magnitude classification.

California↗

Interpretation of well logs in a carbonate aquifer

This report describes the log analysis of the Randolph and Sabial core holes in the Edwards aquifer in Texas, with particular attention to the principles that can be applied generally to any carbonate system. The geologic and hydrologic data were obtained during the drilling of the two holes, from extensive laboratory analysis of the cores, and from numerous geophysical logs run in the two holes. Some logging methods are inherently superiors to others for the analysis of limestone and dolomite aquifers. Three such systems are the dentistry, neutron, and acoustic-velocity (sonic) logs. Most of the log analysis described here is based on the interpretation of suites of logs from these three systems. In certain instances, deeply focused resistivity logs can be used to good advantage in carbonate rock studies; this technique is used to computer the water resistivity in the Randolph core hole. The rocks penetrated by the Randolph core hole are typical of those carbonates that have undergone very little solution by recent ground-water circulation. There are few large solutional openings; the water is saline; and the rocks are dark, dolomitic, have pore space that is interparticle or intercrystalline, and contain unoxidized organic material. The total porosity of rocks in the saline zone is higher than that of rocks in the fresh-water aquifer; however, the intrinsic permeability is much less in the saline zone because there are fewer large solutional openings. The Sabinal core hole penetrates a carbonate environment that has experienced much solution by ground water during recent geologic time. The rocks have high secondary porosities controlled by sedimentary structures within the rock; the water is fresh; and the dominant rock composition is limestone. The relative percentages of limestone and dolomite, the average matrix (grain) densities of the rock mixtures , and the porosity of the rock mass can be calculated from density, neutron, and acoustic logs. With supporting data from resistivity logs, the formation water quality can be estimated, as well as the relative cementation or tortuosity of the rock. Many of these properties calculated from logs can be verified by analysis of the core available from test holes drilled in the saline and fresh water zones.

Texas↗

Defining winter trophic habitat of juvenile Gulf Sturgeon in the Suwannee and Apalachicola rivermouth estuaries, acoustic telemetry investigations

Three automated listening post-telemetry studies were undertaken in the Suwannee and Apalachicola estuaries to gain knowledge of habitats use by juvenile Gulf Sturgeons ( Acipenser oxyrinchus desotoi ) on winter feeding grounds. A simple and reliable method for external attachment of small acoustic tags to the dorsal fin base was developed using shrink-tubing. Suspending receivers on masts below anchored buoys improved reception and facilitated downloading; a detection range of 500–2500 m was realized. In the Apalachicola estuary, juvenile GS stayed in shallow water (< 2 m) within the estuarine transition zone all winter in the vicinity of the Apalachicola River mouth. Juvenile GS high-use areas did not coincide with high density benthic macrofauna areas from the most recent (1999) benthos survey. In the Suwannee estuary, juveniles ranged widely and individually throughout oligohaline to mesohaline subareas of the estuary, preferentially using mesohaline subareas seaward of Suwannee Reef (52% of acoustic detections). The river mouth subarea was important only in early and late winter, during the times of adult Gulf Sturgeon migrations (41% of detections). Preferred winter feeding subareas coincided spatially with known areas of dense macrofaunal benthos concentrations. Following a dramatic drop in air and water temperatures, juvenile GS left the river mouth and estuary, subsequently being detected 8 km offshore in polyhaline open Gulf of Mexico waters, before returning to the estuary. Cold-event offshore excursions demonstrate that they can tolerate full-salinity polyhaline waters in the open Gulf of Mexico, for at least several days at a time. For juvenile sturgeons, the stress and metabolic cost of enduring high salinity (Jarvis et al., 2001; McKenzie et al., 2001; Singer and Ballantyne, 2002) for short periods in deep offshore waters seems adaptively advantageous relative to the risk of cold-event mortality in shallow inshore waters of lower salinity. Thus, while juveniles can tolerate high salinities for days to weeks to escape cold events, they appear to make only infrequent use of open polyhaline waters. Throughout the winter foraging period, juvenile GS stayed primarily within the core area of Suwannee River mouth influence, extending about 12 km north and south of the river mouth, and somewhat seaward of Suwannee Reef (< 5 km offshore). None were detected departing the core area past either of the northern or southern acoustic gates, located 66 and 52 km distant from the river mouth, respectively.

Florida↗

Evaluation of a mass-balance approach to determine consumptive water use in northeastern Illinois

A principal component of evaluating and managing water use is consumptive use. This is the portion of water withdrawn for a particular use, such as residential, which is evaporated, transpired, incorporated into products or crops, consumed by humans or livestock, or otherwise removed from the immediate water environment. The amount of consumptive use may be estimated by a water (mass)-balance approach; however, because of the difficulty of obtaining necessary data, its application typically is restricted to the facility scale. The general governing mass-balance equation is: Consumptive use = Water supplied - Return flows . This study explored a mass-balance field-based computation of consumptive use in a residential setting at the scale of a sanitary sewer service area (sewershed). In addition, the feasibility (cost and difficulty) and relative uncertainties (accuracies) associated with applying the approach at this scale were evaluated. The study was conducted during 2011&ndash;13 within a 3.5-square mile (mi 2 ) sewershed confined to a predominantly residential area of Elk Grove Village, Illinois. Following background evaluation of the geohydrologic setting, sewershed infrastructure, and possible components of supplied and returned water, the identified primary components were 1. public water deliveries by the Elk Grove Village Department of Public Works, 2. self-served groundwater withdrawals in an included unincorporated neighborhood with public sanitary sewer service, 3. return flows to the sanitary sewer system, and 4. direct return of water discharged from swimming pools to Salt Creek. Water volumes principally were reported for deliveries, measured for sanitary sewer returns by using an acoustic Doppler current-velocity meter, and estimated for domestic withdrawals and swimming pool discharges to storm sewers. All water volumes required some degree of estimation. Observation wells were installed adjacent to sewer pipelines (lines) to determine the depth of the water table relative to that of the sewer lines and to collect water samples for detection of optical brighteners, as they are routinely discharged as clotheswashing waste to sanitary sewers. These data provided qualitative information on gains (inflow and infiltration) and losses (exfiltration) of sewer flow by pipe leakage, which might otherwise not be considered in the sewer flow return measurements. Hydrographs of sewer flow also were evaluated to identify and estimate storm-associated inputs to sewer flow. The volume of sanitary sewer return flow (778 million gallons per year [Mgal/yr]) was determined to substantially exceed the volume of supplied water (566 Mgal/yr), thus, for this study setting, voiding the utility of the applied mass-balance approach for estimating consumptive water use. Mass-balance components, including sanitary sewer flow and supplied-water use, were estimated within reasonable limits of uncertainty. Evidence of a water table that is typically shallower than the area&rsquo;s sewer lines, yet is sometimes depressed near more deeply buried sewer lines, suggests groundwater infiltration into the sewers contributes to the excess volume of return flow. Technical obstacles and project resources precluded accurate quantification of infiltration volumes and other gains and losses to sanitary sewer flow. As estimated from various simplified methods, a minimum of 26 percent of return flow measured in the sanitary sewer represented groundwater infiltration and stormwater inflow; separately, about 2 percent of return flow was estimated as inflow. On the basis of the alternative winter base-rate method, consumptive use in the sewershed was estimated as 13 percent, which compares favorably with that used by the State of Illinois for Lake Michigan allocation accounting (10 percent) and other States and Canadian Provinces in the Great Lakes region (generally 10-15 percent). The study also provided other findings considered useful to studies of water use and to performance evaluation of sanitary sewer infrastructure. In urban residential settings, the comparatively small volumes of nonpublic sources of water (self-supplied) and direct (nonsanitary) return flow potentially can be ignored in the estimation of consumptive use. An acoustic Doppler current-velocity meter can be used in sanitary sewers to accurately measure discharge and reasonably estimate storm-associated inflows. Hourly to daily patterns of water use can be readily identified and quantified in the return flow record for the sanitary sewers. Relative volumes of infiltration gains (and exfiltration losses) can be substantial, even in sewer systems of communities making significant investments in system upgrades to limit sewer line leakage. Monitoring of optical brighteners in groundwater (and potentially in sanitary sewer flow) can provide a useful means of identifying probable leakage from (and to) sewer lines. Accurate quantification of gains and losses to sanitary sewer flow at the sewershed scale will require additional research effort and technical advances. Under ideal conditions, accurate quantification of consumptive use at the sewershed scale by the described mass-balance approach might be possible. Under most prevailing conditions, quantification likely would be more costly and time consuming than that of the present study, given the freely contributed technical support of the host community and relatively appropriate conditions of the study area. Essentials to quantification of consumptive use are a fully cooperative community, storm and sanitary sewers that are separate, and newer sewer infrastructure and (or) a robust program for limiting infiltration, exfiltration, and inflow.

Illinois↗

Measurement of bedload transport in sand-bed rivers: A look at two indirect sampling methods

Sand-bed rivers present unique challenges to accurate measurement of the bedload transport rate using the traditional direct sampling methods of direct traps (for example the Helley-Smith bedload sampler). The two major issues are: 1) over sampling of sand transport caused by &ldquo;mining&rdquo; of sand due to the flow disturbance induced by the presence of the sampler and 2) clogging of the mesh bag with sand particles reducing the hydraulic efficiency of the sampler. Indirect measurement methods hold promise in that unlike direct methods, no transport-altering flow disturbance near the bed occurs. The bedform velocimetry method utilizes a measure of the bedform geometry and the speed of bedform translation to estimate the bedload transport through mass balance. The bedform velocimetry method is readily applied for the estimation of bedload transport in large sand-bed rivers so long as prominent bedforms are present and the streamflow discharge is steady for long enough to provide sufficient bedform translation between the successive bathymetric data sets. Bedform velocimetry in small sandbed rivers is often problematic due to rapid variation within the hydrograph. The bottom-track bias feature of the acoustic Doppler current profiler (ADCP) has been utilized to accurately estimate the virtual velocities of sand-bed rivers. Coupling measurement of the virtual velocity with an accurate determination of the active depth of the streambed sediment movement is another method to measure bedload transport, which will be termed the &ldquo;virtual velocity&rdquo; method. Much research remains to develop methods and determine accuracy of the virtual velocity method in small sand-bed rivers.

Missouri↗

Overview of selected surrogate technologies for continuous suspended-sediment monitoring

Surrogate technologies for inferring selected characteristics of suspended sediments in surface waters are being tested by the U.S. Geological Survey and several partners with the ultimate goal of augmenting or replacing traditional monitoring methods. Optical properties of water such as turbidity and optical backscatter are the most commonly used surrogates for suspended-sediment concentration, but use of other techniques such as those based on acoustic backscatter, laser diffraction, digital photo-optic, and pressure-difference principles is increasing for concentration and, in some cases, particle-size distribution and flux determinations. The potential benefits of these technologies include acquisition of automated, continuous, quantifiably accurate data obtained with increased safety and at less expense. When suspended-sediment surrogate data meet consensus accuracy criteria and appropriate sediment-record computation techniques are applied, these technologies have the potential to revolutionize the way fluvial-sediment data are collected, analyzed, and disseminated.

Conference Paper↗

Acoustic and optical borehole-wall imaging for fractured-rock aquifer studies

Imaging with acoustic and optical televiewers results in continuous and oriented 360° views of the borehole wall from which the character, relation, and orientation of lithologic and structural planar features can be defined for studies of fractured-rock aquifers. Fractures are more clearly defined under a wider range of conditions on acoustic images than on optical images including dark-colored rocks, cloudy borehole water, and coated borehole walls. However, optical images allow for the direct viewing of the character of and relation between lithology, fractures, foliation, and bedding. The most powerful approach is the combined application of acoustic and optical imaging with integrated interpretation. Imaging of the borehole wall provides information useful for the collection and interpretation of flowmeterand other geophysical logs, core samples, and hydraulic and water-quality data from packer testing and monitoring.

Journal of Applied Geophysics↗

How accurately does eDNA reflect the spatial distribution of cold-water fish? Field validation from a temperate lake

Applications of environmental DNA (eDNA) based detection technology to evaluate the distribution of aquatic organisms are increasing; yet field validations of eDNA are important to measure accuracy in study systems. To successfully apply this technology to species conservation, it is critical to understand how both species biology and environmental conditions affect the accuracy of inference from eDNA detection data. We implemented a field assessment of the accuracy and spatial resolution of eDNA-based species distributions for a native cold-water, schooling fish, cisco Coregonus artedi , that has been reintroduced to a deep temperate lake. We leveraged a combination of acoustic telemetry, providing known spatial locations of tagged fish, and lake-wide eDNA sampling to infer their distribution in Keuka Lake, New York, USA. Sub-surface (12 m and 18 m depths) eDNA samples were collected to accommodate the diel vertical migration behaviour of this fish species. The results of this study validated the accuracy of positive eDNA detections with the distribution of tagged fish to coarse spatial scales. Yet, several fine-scale locations revealed a mismatch between eDNA and acoustic telemetry detections; consistent with rapid transport of genetic material via lake currents. Empirical measurements of lake currents using drifters found cisco eDNA detections could deviate from specimens' source locations by as much as 3.3 km at 12 m depth or 1.5 km at 18 m depth over a 24 h transport period. Our study indicates that accurate species distributions estimated from eDNA sampling in lakes may require further understanding of transport mechanisms and persistence of environmental genetic material to relate point detections to source animal locations. Integrating eDNA sampling with additional data collection of species biology and environmental conditions will increase the spatial resolution of fish distribution assessments.

New York↗

First documented occurrences of the shortnose sturgeon (Acipenser brevirostrum) in the Saco River, Maine, USA

During sampling efforts to study the more abundant Atlantic sturgeon, Acipenser oxyrinchus oxyrinchus, between May of 2009 and November of 2011, four shortnose sturgeon were captured in gill nets near the mouth of the Saco River, Maine. Two of these individuals were tagged with acoustic transmitters to monitor their movement within the Saco River. Additionally, six shortnose sturgeon that had been tagged with acoustic transmitters in the Merrimack River, Massachusetts were detected on the acoustic array deployed within the Saco River and its estuary over this time period. These incidences represent the first verified documentation of shortnose sturgeon within this estuary.

Maine↗

Bibliography of borehole geophysics as applied to ground-water hydrology

Most of the references on borehole geophysics that are relevant to ground-water hydrology are contained in this bibliography, but it does not include every reference that is available under each subject heading; the literature is much too extensive to compile a complete listing. Some of the references may appear under more than one subject heading because the references commonly relate to more than one main topic. Many articles have been cross-referenced in order to assist the reader in locating an article. For example, the article entitled, 'Application of the acoustic televiewer to the characterization of hydraulic fractures in geothermal wells' is listed under both 'Acoustic televiewer,' and 'Geothermal'. The bibliography is intended to lead the reader to other articles on borehole-geophysical logging and related subjects, because each article cited also will have a list of references, which may be more specialized, covering many subjects with related applications, such as physics, mathematics, chemistry, geology, electronics, acoustics, hydrology, and surface geophysics. However, not all of these related subject headings could be included in this bibliography.

Circular↗

Sharp structural variability of the Gorda slab imaged by a fiber array

Structural variability at the subduction interface imposes stress and strength heterogeneity that affects the behavior of megathrust earthquakes. Characterizing this variability can be challenging because the slab interface is often deep and seismically quiet. Seismic wavefields from intraslab earthquakes can contain high-frequency phase conversions from the megathrust fault that are delayed and distorted by interface topography and velocity heterogeneity. Distributed acoustic sensing (DAS) allows for long-term deployments of dense seismic arrays that can yield many observations of these converted phases at high spatial resolution over large areas. Here, we develop a technique to identify and enhance slab-converted phases in DAS data, and we pick the differential arrival times and relative polarities between the direct and converted phases at short wavelength along the fiber optic cable. We apply this framework to earthquake wavefields measured by a DAS array at the southern end of the Cascadia Subduction Zone, above the locked zone, which suggests kilometer-scale interface roughness and variations in crustal thickness and absolute depth of the Gorda slab. Converted phase wavefields also suggest sharp variability in the presence and position of fluids in the slab and forearc crust.

California↗

Accuracy of flowmeters measuring horizontal groundwater flow in an unconsolidated aquifer simulator.

Borehole flowmeters that measure horizontal flow velocity and direction of groundwater flow are being increasingly applied to a wide variety of environmental problems. This study was carried out to evaluate the measurement accuracy of several types of flowmeters in an unconsolidated aquifer simulator. Flowmeter response to hydraulic gradient, aquifer properties, and well-screen construction was measured during 2003 and 2005 at the U.S. Geological Survey Hydrologic Instrumentation Facility in Bay St. Louis, Mississippi. The flowmeters tested included a commercially available heat-pulse flowmeter, an acoustic Doppler flowmeter, a scanning colloidal borescope flowmeter, and a fluid-conductivity logging system. Results of the study indicated that at least one flowmeter was capable of measuring borehole flow velocity and direction in most simulated conditions. The mean error in direction measurements ranged from 15.1 degrees to 23.5 degrees and the directional accuracy of all tested flowmeters improved with increasing hydraulic gradient. The range of Darcy velocities examined in this study ranged 4.3 to 155 ft/d. For many plots comparing the simulated and measured Darcy velocity, the squared correlation coefficient (r 2 ) exceeded 0.92. The accuracy of velocity measurements varied with well construction and velocity magnitude. The use of horizontal flowmeters in environmental studies appears promising but applications may require more than one type of flowmeter to span the range of conditions encountered in the field. Interpreting flowmeter data from field settings may be complicated by geologic heterogeneity, preferential flow, vertical flow, constricted screen openings, and nonoptimal screen orientation.

Ground Water Monitoring and Remediation↗

A numerical simulation of magma motion, crustal deformation, and seismic radiation associated with volcanic eruptions

The finite difference method is used to calculate the magma dynamics, seismic radiation, and crustal deformation associated with a volcanic eruption. The model geometry consists of a cylindrical reservoir and narrow cylindrical conduit embedded in a homogeneous crust. We consider two models of eruption. In the first model, a lid caps the vent and the magma is overpressurized prior to the eruption. The eruption is triggered by the instantaneous removal of the lid, at which point the exit pressure becomes equal to the atmospheric pressure. In the second model, a plug at the reservoir outlet allows pressurization of only the magmatic fluid in the reservoir before the eruption. Magma transfer between the reservoir and conduit is triggered by the instantaneous removal of the plug, and the eruption occurs when the pressure at the conduit orifice exceeds the material strength of the lid capping the vent. In both models, magma dynamics are expressed by the equations of mass and momentum conservation in a compressible fluid, in which fluid expansion associated with depressurization is accounted for by a constitutive law relating pressure and density. Crustal motions are calculated from the equations of elastodynamics. The fluid and solid are dynamically coupled by applying the continuity of wall velocities and normal stresses across the conduit and reservoir boundaries. Free slip is allowed at the fluid-solid boundary. Both models predict the gradual depletion of the magma reservoir, which causes crustal deformation observed as a long-duration dilatational signal. Superimposed on this very-long-period (VLP) signal generated by mass transport are long-period (LP) oscillations of the magma reservoir and conduit excited by the acoustic resonance of the reservoir-conduit system during the eruption. The volume of the reservoir, vent size, and magma properties control the duration of VLP waves and dominant periods of LP oscillations. The second model predicts that when the magmatic fluid reaches the vent, a high-pressure pulse occurs at this location in accordance with the basic theory of compressible fluid dynamics. This abrupt pressure increase just beneath the vent is consistent with observed seismograms in which pulse-like Rayleigh waves excited by a shallow source are dominant. The strength of the lid plays an important role in the character of the seismograms and in defining the type of eruption observed.

Geophysical Journal International↗

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth↗

Submarine slope failures near Seward, Alaska, during the M9.2 1964 earthquake

Following the 1964 M9.2 megathrust earthquake in southern Alaska, Seward was the only town hit by tsunamis generated from both submarine landslides and tectonic sources. Within 45 seconds of the start of the earthquake, a 1.2-km-long section of waterfront began sliding seaward, and soon after, ~6-8-m high waves inundated the town. Studies soon after the earthquake concluded that submarine landslides along the Seward waterfront generated the tsunamis that occurred immediately after the earthquake. We analyze pre- and post-earthquake bathymetry data to assess the location and extent of submarine mass failures and sediment transport. New NOAA multibeam bathymetry shows the morphology of the entire fjord at 15 m resolution. We also assembled all older soundings from smooth sheets for comparison to the multibeam dataset. We gridded the sounding data, applied corrections for coseismic subsidence, post-seismic rebound, unrecovered co-seismic subsidence, sea-level rise (vertical datum shift), and measurement errors. The difference grids show changes resulting from the 1964 earthquake. We estimate the total volume of slide material to be about 211 million m3. Most of this material was transported to a deep, flat area, which we refer to as &ldquo;the bathtub&rdquo;, about 6 to 13 km south of Seward. Sub-bottom profiling of the bathtub shows an acoustically transparent unit, which we interpret as a sediment flow deposit resulting from the submarine landslides. The scale of the submarine landslides and the distance over which sediment was transported is much larger than previously appreciated.

Alaska↗

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗

Local explosion detection and infrasound localization by reverse time migration using 3-D finite-difference wave propagation

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (<15 km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ∼20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.

Alaska↗

An initial investigation of multidimensional flow and transverse mixing characteristics of the Ohio River near Cincinnati, Ohio

Two-dimensional hydrodynamic and transport models were applied to a 34-mile reach of the Ohio River from Cincinnati, Ohio, upstream to Meldahl Dam near Neville, Ohio. The hydrodynamic model was based on the generalized finite-element hydrodynamic code RMA2 to simulate depth-averaged velocities and flow depths. The generalized water-quality transport code RMA4 was applied to simulate the transport of vertically mixed, water-soluble constituents that have a density similar to that of water. Boundary conditions for hydrodynamic simulations included water levels at the U.S. Geological Survey water-level gaging station near Cincinnati, Ohio, and flow estimates based on a gate rating at Meldahl Dam. Flows estimated on the basis of the gate rating were adjusted with limited flow-measurement data to more nearly reflect current conditions. An initial calibration of the hydrodynamic model was based on data from acoustic Doppler current profiler surveys and water-level information. These data provided flows, horizontal water velocities, water levels, and flow depths needed to estimate hydrodynamic parameters related to channel resistance to flow and eddy viscosity. Similarly, dye concentration measurements from two dye-injection sites on each side of the river were used to develop initial estimates of transport parameters describing mixing and dye-decay characteristics needed for the transport model. A nonlinear regression-based approach was used to estimate parameters in the hydrodynamic and transport models. Parameters describing channel resistance to flow (Manning’s “n”) were estimated in areas of deep and shallow flows as 0.0234, and 0.0275, respectively. The estimated RMA2 Peclet number, which is used to dynamically compute eddy-viscosity coefficients, was 38.3, which is in the range of 15 to 40 that is typically considered appropriate. Resulting hydrodynamic simulations explained 98.8 percent of the variability in depth-averaged flows, 90.0 percent of the variability in water levels, 93.5 percent of the variability in flow depths, and 92.5 percent of the variability in velocities. Estimates of the water-quality-transport-model parameters describing turbulent mixing characteristics converged to different values for the two dye-injection reaches. For the Big Indian Creek dye-injection study, an RMA4 Peclet number of 37.2 was estimated, which was within the recommended range of 15 to 40, and similar to the RMA2 Peclet number. The estimated dye-decay coefficient was 0.323. Simulated dye concentrations explained 90.2 percent of the variations in measured dye concentrations for the Big Indian Creek injection study. For the dye-injection reach starting downstream from Twelvemile Creek, however, an RMA4 Peclet number of 173 was estimated, which is far outside the recommended range. Simulated dye concentrations were similar to measured concentration distributions at the first four transects downstream from the dye-injection site that were considered vertically mixed. Farther downstream, however, simulated concentrations did not match the attenuation of maximum concentrations or cross-channel transport of dye that were measured. The difficulty of determining a consistent RMA4 Peclet was related to the two-dimension model assumption that velocity distributions are closely approximated by their depth-averaged values. Analysis of velocity data showed significant variations in velocity direction with depth in channel reaches with curvature. Channel irregularities (including curvatures, depth irregularities, and shoreline variations) apparently produce transverse currents that affect the distribution of constituents, but are not fully accounted for in a two-dimensional model. The two-dimensional flow model, using channel resistance to flow parameters of 0.0234 and 0.0275 for deep and shallow areas, respectively, and an RMA2 Peclet number of 38.3, and the RMA4 transport model with a Peclet number of 37.2, may have utility for emergency-planning purposes. Emergency-response efforts would be enhanced by continuous streamgaging records downstream from Meldahl Dam, real-time water-quality monitoring, and three-dimensional modeling. Decay coefficients are constituent specific.

Kentucky, Ohio↗