SEARCH · Geology Reports

Results for “Interfaces”

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

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27Linked to original sources

USGS2013

Quantifying potential earthquake and tsunami hazard in the Lesser Antilles subduction zone of the Caribbean region

In this study, we quantify the seismic and tsunami hazard in the Lesser Antilles subduction zone, focusing on the plate interface offshore of Guadeloupe. We compare potential strain accumulated via GPS-derived plate motions to strain release due to earthquakes that have occurred over the past 110 yr, and compute the resulting moment deficit. Our results suggest that enough strain is currently stored in the seismogenic zone of the Lesser Antilles subduction arc in the region of Guadeloupe to cause a large and damaging earthquake of magnitude M w ∼ 8.2 ± 0.4. We model several scenario earthquakes over this magnitude range, using a variety of earthquake magnitudes and rupture areas, and utilizing the USGS ShakeMap and PAGER software packages. Strong ground shaking during the earthquake will likely cause loss of life and damage estimated to be in the range of several tens to several hundreds of fatalities and hundreds of millions to potentially billions of U.S. dollars of damage. In addition, such an event could produce a significant tsunami. Modelled tsunamis resulting from these scenario earthquakes predict meter-scale wave amplitudes even for events at the lower end of our magnitude range (M 7.8), and heights of over 3 m in several locations with our favoured scenario (M 8.0, partially locked interface from 15–45 km depth). In all scenarios, only short lead-times (on the order of tens of minutes) would be possible in the Caribbean before the arrival of damaging waves.

Lesser Antilles↗

USGS1986

Shuttle imaging radar: Physical controls on signal penetration and subsurface scattenng in the Eastern Sahara

SIR-A signal penetration and subsurface backscatter within the upper meter or so of the sediment blanket in the Eastern Sahara of southern Egypt and northern Sudan are enhanced both by radar sensor parameters and by the physical and chemical characteristics of eolian and alluvial materials. Interpretation of SIR-A images by McCauley et al. [1], [2] dramatically changed previous concepts of the role that fluvial processes have played over the past 10000 to 30000000 years in shaping this now extremely flat, featureless, and hyperarid landscape. In the present paper we summarize the near-surface stratigraphy, the electrical properties of materials, and the types of radar interfaces found to be responsible for different classes of SIR-A tonal response. The dominant factors related to efficient microwave signal penetration into the sediment blanket include a) favorable distribution of particle sizes, b) extremely low moisture content and c) reduced geometric scattering at the SIR-A frequency (1.3 GHz). The depth of signal penetration that results in a recorded backscatter, here called "radar imaging depth," was documented in the field to be a maximum of 1.5 m, or 0.25 of the calculated "skin depth," for the sediment blanket. Radar imaging depth is estimated to be between 2 and 3 m for active sand dune materials. Diverse permittivity interfaces and volume scatterers within the shallow subsurface are responsible for most of the observed backscatter not directly attributable to grazing outcrops.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2011

A trans-dimensional Bayesian Markov chain Monte Carlo algorithm for model assessment using frequency-domain electromagnetic data

A meaningful interpretation of geophysical measurements requires an assessment of the space of models that are consistent with the data, rather than just a single, ‘best’ model which does not convey information about parameter uncertainty. For this purpose, a trans-dimensional Bayesian Markov chain Monte Carlo (MCMC) algorithm is developed for assessing frequencydomain electromagnetic (FDEM) data acquired from airborne or ground-based systems. By sampling the distribution of models that are consistent with measured data and any prior knowledge, valuable inferences can be made about parameter values such as the likely depth to an interface, the distribution of possible resistivity values as a function of depth and non-unique relationships between parameters. The trans-dimensional aspect of the algorithm allows the number of layers to be a free parameter that is controlled by the data, where models with fewer layers are inherently favoured, which provides a natural measure of parsimony and a signiﬁcant degree of ﬂexibility in parametrization. The MCMC algorithm is used with synthetic examples to illustrate how the distribution of acceptable models is affected by the choice of prior information, the system geometry and conﬁguration and the uncertainty in the measured system elevation. An airborne FDEM data set that was acquired for the purpose of hydrogeological characterization is also studied. The results compare favorably with traditional least-squares analysis, borehole resistivity and lithology logs from the site, and also provide new information about parameter uncertainty necessary for model assessment.

Geophysical Journal International↗

USGS2011

Evaluation of MODFLOW-LGR in connection with a synthetic regional-scale model

This work studies costs and benefits of utilizing local‐grid refinement (LGR) as implemented in MODFLOW‐LGR to simulate groundwater flow in a buried tunnel valley interacting with a regional aquifer. Two alternative LGR methods were used: the shared‐node (SN) method and the ghost‐node (GN) method. To conserve flows the SN method requires correction of sources and sinks in cells at the refined/coarse‐grid interface. We found that the optimal correction method is case dependent and difficult to identify in practice. However, the results showed little difference and suggest that identifying the optimal method was of minor importance in our case. The GN method does not require corrections at the models' interface, and it uses a simpler head interpolation scheme than the SN method. The simpler scheme is faster but less accurate so that more iterations may be necessary. However, the GN method solved our flow problem more efficiently than the SN method. The MODFLOW‐LGR results were compared with the results obtained using a globally coarse (GC) grid. The LGR simulations required one to two orders of magnitude longer run times than the GC model. However, the improvements of the numerical resolution around the buried valley substantially increased the accuracy of simulated heads and flows compared with the GC simulation. Accuracy further increased locally around the valley flanks when improving the geological resolution using the refined grid. Finally, comparing MODFLOW‐LGR simulation with a globally refined (GR) grid showed that the refinement proportion of the model should not exceed 10% to 15% in order to secure method efficiency.

Ground Water↗

USGS1991

Strain measurements and the potential for a great subduction earthquake off the coast of Washington

Geodetic measurements of deformation in northwestern Washington indicate that strain is accumulating at a rate close to that predicted by a model of the Cascadia subduction zone in which the plate interface underlying the continental slope and outer continental shelf is currently locked but the remainder of the interface slips continuously. Presumably this locked segment will eventually rupture in a great thrust earthquake with a down-dip extent greater than 100 kilometers.

Washington↗

USGS1966

Tropical lakes, copropel, and oil shale

During a long-continued study of the lacustrine beds of the Eocene Green River Formation, I have tried to interpret past events from observation of present-day processes. After a search of some 40 years, four lakes have been found that are producing a kind of organic ooze judged to be a modern analogue of the precursors of rich oil shale. Two of the lakes are in central Africa and two are in Florida. All four are shallow. The ooze in all four is predominantly algal, entirely in the form of minute fecal pellets, and does not decay in warm, wet, oxidizing environments. Several of the most unusual, mummified microorganisms found in the oil shale of the Green River are illustrated as testament to the inference that the Eocene organic oozes also were resistant to decay. Studies to determine why these algal oozes do not decay are in progress, but as yet no satisfactory explanation is available. The ooze from Mud Lake, Florida, contains very few living bacteria but a great many bacterial spores, suggesting some active inhibitor. Gentle and slow anaerobic decay takes place in the ooze 1 foot or more below the mud-water interface. The algal ooze accumulates slowly. That at a depth of 3 feet below the mud-water interface has a C 14 age of 2280 ± 200 years. If compacted, this 3-foot layer would amount to a layer only about 0.5 inch thick. The air-dried algal mud (from Mud Lake) looks much like oil shale and has a C-H ratio essentially like that of the organic matter in oil shale. The oxygen of the mud, however, is roughly 5 times as high as in oil shale. The calorific value of dried Mud Lake algal ooze is about 6600 cal/g, whereas the organic matter from Green River oil shale averages about 9500 cal/g. Analyses show that the dried algal ooze from Mud Lake contains small quantities of higher fatty acids (C 12 –C 34 ), with C 16 being dominant. It also contains some unsaturated fatty acids and about 0.3 per cent of n-alkanes, predominantly odd carbons, with C 29 being dominant. A few qualitative analyses show that carotenoid pigments and terpenes are also present. A possible source of hydrocarbons is the long branched side chain of the chlorophyll molecule. Small crustaceans liberate this in the form of phytol, which, by dehydration, can go over into a series of saturated and unsaturated hydrocarbons, including phytanc and pristane, both of which are common in Green River oil shale. Search for other precursors of hydrocarbons continues.

GSA Bulletin↗

USGS1970

Quartz gabbro and anorthositic gabbro: Markers of offset along the San Andreas fault in the California Coast Ranges

Large-scale lateral movement on the San Andreas fault zone is suggested by the distribution of gabbroic rocks that may be slivered remnants of oceanic crust. Distinctive and unusual hornblende quartz gabbro and anorthositic gabbro that are virtually identical both petrographically and chemically are exposed at Logan and Gold Hill in the Coast Ranges along the San Andreas fault. The hornblende quartz gabbro is made up of labradorite to bytownite plagioclase, hornblende, and quartz, with very minor biotite and pyroxene, and accessory metallic opaques and apatite. The coarse-grained anorthositic gabbro with anorthositic layers and associated gneiss is made up mainly of labradorite to bytownite plagioclase, hornblende, lesser clinopyroxene, and, locally, orthopyroxene. The present 100 mi of separation between these two gabbro bodies is probably due to displacement along the San Andreas fault of what was originally one gabbro mass. Somewhat similar quartz gabbro and anorthositic gabbro associated with ultramafic rocks near Eagle Rest Peak in the San Emigdio Mountains may represent a source for the Logan and Gold Hill slivers. This suggests a minimum right-lateral movement of about 200 mi on the San Andreas fault zone. Anorthositic gabbro clasts from Cretaceous conglomerate in the Gualala area have strong resemblance to gabbroic rocks at Logan, Gold Hill, and Eagle Rest Peak and speculatively suggest 350 mi of right-lateral movement on the fault zone. It is tentatively suggested that the gabbro clasts of Gualala may have been shed from a large area of exposed gabbroic oceanic crust, pieces of which occur as fault slivers at Gold Hill and Logan. It is further suggested that the Eagle Rest Peak locality may be a relatively in-place exposure of this gabbroic oceanic crust that lies near a continental-oceanic interface. This interface, thought to represent a fossil subduction zone between Franciscan and Sierran basements, appears to be overlain depositionally by Eocene sedimentary rocks. If these Eocene deposits do “seal over” the possible subduction zone, it suggests that movement on such a zone ceased before the Eocene deposition, and that the new and different pattern of right-lateral movement on the San Andreas fault zone probably began sometime later.

California↗

USGS2010

Slow slip phenomena in Cascadia from 2007 and beyond: a review

Recent technological advances combined with more detailed analyses of seismologic and geodetic observations have fundamentally changed our understanding of the ways in which tectonic stresses arising from plate motions are accommodated by slip on faults. The traditional view that relative plate motions are accommodated by a simple cycle of stress accumulation and release on “locked” plate-boundary faults has been revolutionized by the serendipitous discovery and recognition of the significance of slow-slip phenomena, mostly in the deeper reaches of subduction zones. The Cascadia subduction zone, located in the Pacific Northwest of the conterminous United States and adjacent Canada, is an archetype of exploration and learning about slow-slip phenomena. These phenomena are manifest as geodetically observed aseismic transient deformations accompanied by a previously unrecognized class of seismic signals. Although secondary failure processes may be involved in generating the seismic signals, the primary origins of both aseismic and seismic phenomena appear to be episodic fault slip, probably facilitated by fluids, on a plate interface that is critically stressed or weakened. In Cascadia, this transient slip evolves more slowly and over more prolonged durations relative to the slip in earthquakes, and it occurs between the 30- and 40-km-depth contours of the plate interface where information was previously elusive. Although there is some underlying organization that relaxes nearly all the accrued plate-motion stresses along the entirety of Cascadia, we now infer that slow slip evolves in complex patterns indicative of propagating stress fronts. Our new understanding provides key constraints not only on the region where the slow slip originates, but also on the probable characteristics of future megathrust earthquakes in Cascadia. Herein, we review the most significant scientific issues and progress related to understanding slow-slip phenomena in Cascadia and highlight some of their societal implications. We provide a comprehensive review, from the big picture as inferred from studies of regional-scale monitoring data to the details revealed by innovative, focused experiments and new instrumentation. We focus on what has been learned largely since 2007, when several major investments in monitoring and temporary deployments dramatically increased the quality and quantity of available data.

Oregon↗

USGS2022

Melt surges, flow differentiation, and remobilization of crystal-rich mushes in response to unloading: Observations from Kīlauea Iki lava lake, Hawaii

Field and laboratory studies of the 1959 Kīlauea Iki lava lake have provided insight into differentiation processes in mafic magma chambers. This paper explores how partially molten basaltic mushes responded to unloading as a consequence of drilling. Most holes drilled from 1967 to 1979 terminated in a melt-rich internal differentiate with a sharp crust-melt interface. These interfaces were not stable, so the boreholes were backfilled by melt-rich (<5% crystal) ooze. This process, with melt ascent rates of 1.3–4.2 m/s, occurred within minutes of intersecting the bodies, mimicking volcanic eruptions, albeit on a small scale. One borehole (KI79-1), which did not encounter such a discontinuity, was backfilled over a period of 16 days by upward flow of crystal-rich mushes rather than melt-rich ooze. The first interval of ooze recovered had undergone extensive internal differentiation. Its most conspicuous feature was production of melt-rich layers by lateral migration of interstitial melt from the wallrock into the rising crystal-rich mush. In addition, two smaller-scale processes occurred within the rising mush: segregation of melt into discrete blebs within the rising mush column and aggregation of groundmass crystals into crystal-rich clumps formed adjacent to coarser olivine crystals. The upper parts of the ooze are enriched in melt relative to deeper samples, which suggests that the melt blebs rose relative to their olivine-rich matrix. Similar melt blebs and crystal-rich clumps are observed in naturally occurring diapiric bodies within the lava lake. These processes appear to be intrinsic to the upwelling of narrow cylindrical mush bodies whether constrained within a borehole (like the oozes) or unconstrained (as were the diapirs in the lava lake). The most striking behavior observed during repeated reentry of KI79-1 was a sharp change in rheology during the second and third re-entries of the borehole. The shift in behavior observed was that the oozes rose up the borehole, with ascent rates of 1.0–1.7 m/s, which are comparable to the rates of the crystal-poor oozes from melt-rich internal differentiates. These oozes contain more melt than the original core at equivalent depths, presumably because melt moved relative to crystals down the pressure gradient created by the open borehole. Groundmass textures in these inflated mushes show erosion of crystal outlines, especially of grain-to-grain contacts between different phases, so that the tenuous crystalline network observed in the original core samples was replaced by rounded crystals in continuous melt at crystallinities of 55–65 vol%. The transition from stable coherent mush to inflatable mush occurred at 25–28 vol% melt. This behavior appears similar to certain types of reactive transport observed in other studies.

Hawaii↗

USGS2016

Extending Alaska's plate boundary: tectonic tremor generated by Yakutat subduction

The tectonics of the eastern end of the Alaska-Aleutian subduction zone are complicated by the inclusion of the Yakutat microplate, which is colliding into and subducting beneath continental North America at near-Pacific-plate rates. The interaction among these plates at depth is not well understood, and further east, even less is known about the plate boundary or the source of Wrangell volcanism. The drop-off in Wadati-Benioff zone (WBZ) seismicity could signal the end of the plate boundary, the start of aseismic subduction, or a tear in the downgoing plate. Further compounding the issue is the possible presence of the Wrangell slab, which is faintly outlined by an anemic, eastward-dipping WBZ beneath the Wrangell volcanoes. In this study, I performed a search for tectonic tremor to map slow, plate-boundary slip in south-central Alaska. I identified ∼11,000 tremor epicenters, which continue 85 km east of the inferred Pacific plate edge marked by WBZ seismicity. The tremor zone coincides with the edges of the downgoing Yakutat terrane, and tremors transition from periodic to continuous behavior as they near the aseismic Wrangell slab. I interpret tremor to mark slow, semicontinuous slip occurring at the interface between the Yakutat and North America plates. The slow slip region lengthens the megathrust interface beyond the WBZ and may provide evidence for a connection between the Yakutat slab and the aseismic Wrangell slab.

Geology↗

USGS2019

Melt equilibration depths as sensors of lithospheric thickness during Eurasia-Arabia collision and the uplift of the Anatolian Plateau

A co-investigation of mantle melting conditions and seismic structure revealed an evolutionary record of mantle dynamics accompanying the transition from subduction to collision along the Africa-Eurasia margin and the >1 km uplift of the Anatolian Plateau. New 40 Ar/ 39 Ar dates of volcanic rocks from the Eastern Taurides (southeast Turkey) considerably expand the known spatial extent of Miocene-aged mafic volcanism following a magmatic lull over much of Anatolia that ended at ca. 20 Ma. Mantle equilibration depths for these chemically diverse basalts are interpreted to indicate that early to middle Miocene lithospheric thickness in the region varied from ∼50 km or less near the Bitlis suture zone to ∼80 km near the Inner Tauride suture zone. This southward-tapering lithospheric base could be a vestige of the former interface between the subducted (and now detached) portion of the Arabian plate and the overriding Eurasian plate, and/or a reflection of mantle weakening associated with greater mantle hydration trenchward prior to collision. Asthenospheric upwelling driven by slab tearing and foundering along this former interface, possibly accompanied by convective removal of the lithosphere, could have led to renewed volcanic activity after 20 Ma. Melt equilibration depths for late Miocene and Pliocene basalts together with seismic imaging of the present lithosphere indicate that relatively invariant lithospheric thicknesses of 60–70 km have persisted since the middle Miocene. Thus, no evidence is found for large-scale (tens of kilometers) Miocene delamination of the lower lithosphere from the overriding plate, which has been proposed elsewhere to account for late Miocene and younger uplift of Anatolia.

Anatolian Plate↗

USGS2019

The Shumagin seismic gap structure and associated tsunami hazards, Alaska convergent margin

The potential for a major earthquake in the Shumagin seismic gap, and the tsunami it could generate, was reported in 1971. However, while potentially tsunamigenic splay faults in the adjacent Unimak and Semidi earthquake segments are known, such features along the Shumagin segment were undocumented until recently. To investigate margin structure and search for splay faults, we reprocessed six legacy seismic records and also processed seismic data acquired by RV Langseth during the ALEUT project (cf. Bécel et al., 2017). All records show splay faults separating the frontal prism from the margin framework. A ridge uplifted by the splay fault hanging wall extends along the entire segment. At the plate interface, the splay fault cuts across subducted sediment strata in some images, whereas in others, the plate interface sediment cuts across the fault. Splay fault zones are commonly associated with subducting lower-plate relief. Along the upper slope, beneath a sediment cover, major normal faults dipping landward and seaward border a ridge of basement rock. The faults displace a regional unconformity that elsewhere received Oligocene–Miocene sediment. Low seafloor scarps above some normal faults indicate recent tectonism. The buried ridge is a continuation of the Unimak Ridge structure that extends NE of the Unimak/Shumagin segment boundary. Some geological characteristics of the Shumagin segment differ from those of other Alaskan earthquake segments, but a causal link to the proposed Shumagin creeping seismic behavior is equivocal.

Alaska↗

USGS2017

Development of a coastal drought index using salinity data

A critical aspect of the uniqueness of coastal drought is the effects on the salinity dynamics of creeks, rivers, and estuaries. The location of the freshwater–saltwater interface along the coast is an important factor in the ecological and socioeconomic dynamics of coastal communities. Salinity is a critical response variable that integrates hydrologic and coastal dynamics including sea level, tides, winds, precipitation, streamflow, and tropical storms. The position of the interface determines the composition of freshwater and saltwater aquatic communities as well as the freshwater availability for water intakes. Many definitions of drought have been proposed, with most describing a decline in precipitation having negative impacts on the water supply. Indices have been developed incorporating data such as rainfall, streamflow, soil moisture, and groundwater levels. These water-availability drought indices were developed for upland areas and may not be ideal for characterizing coastal drought. The availability of real-time and historical salinity datasets provides an opportunity for the development of a salinity-based coastal drought index. An approach similar to the standardized precipitation index (SPI) was modified and applied to salinity data obtained from sites in South Carolina and Georgia. Using the SPI approach, the index becomes a coastal salinity index (CSI) that characterizes coastal salinity conditions with respect to drought periods of higher-saline conditions and wet periods of higher-freshwater conditions. Evaluation of the CSI indicates that it provides additional coastal response information as compared to the SPI and the Palmer hydrologic drought index, and the CSI can be used for different estuary types and for comparison of conditions along coastlines.

Bulletin of the American Meteorological Society↗

USGS2021

Spectral damping scaling factors for horizontal components of ground motions from subduction earthquakes using NGA-Subduction data

This article develops global models of damping scaling factors (DSFs) for subduction zone earthquakes that are functions of the damping ratio, spectral period, earthquake magnitude, and distance. The Next Generation Attenuation for subduction earthquakes (NGA-Sub) project has developed the largest uniformly processed database of recorded ground motions to date from seven subduction regions: Alaska, Cascadia, Central America and Mexico, South America, Japan, Taiwan, and New Zealand. NGA-Sub used this database to develop new ground motion models (GMMs) at a reference 5% damping ratio. We worked with the NGA-Sub project team to develop an extended database that includes pseudo-spectral accelerations (PSA) for 11 damping ratios between 0.5% and 30%. We use this database to develop parametric models of DSF for both interface and intraslab subduction earthquakes that can be used to adjust any subduction GMM from a reference 5% damping ratio to other damping ratios. The DSF is strongly influenced by the response spectral shape and the duration of motion; therefore, in addition to the damping ratio, the median DSF model uses spectral period, magnitude, and distance as surrogate predictor variables to capture the effects of the spectral shape and the duration of motion. We also develop parametric models for the standard deviation of DSF. The models presented in this article are for the RotD50 horizontal component of PSA and are compared with the models for shallow crustal earthquakes in active tectonic regions. Some noticeable differences arise from the considerably longer duration of interface records for very large magnitude events and the enriched high-frequency content of intraslab records, compared with shallow crustal earthquakes. Regional differences are discussed by comparing the proposed global models with the data from each subduction region along with recommendations on the applicability of the models.

Earthquake Spectra↗

USGS2021

NGA-subduction global ground motion models with regional adjustment factors

We develop semi-empirical ground motion models (GMMs) for peak ground acceleration, peak ground velocity, and 5%-damped pseudo-spectral accelerations for periods from 0.01 to 10 s, for the median orientation-independent horizontal component of subduction earthquake ground motion. The GMMs are applicable to interface and intraslab subduction earthquakes in Japan, Taiwan, Mexico, Central America, South America, Alaska, the Aleutian Islands, and Cascadia. The GMMs are developed using a combination of data inspection, data regression with respect to physics-informed functions, ground-motion simulations, and geometrical constraints for certain model components. The GMMs capture observed differences in source and path effects for interface and intraslab events, conditioned on moment magnitude, rupture distance, and hypocentral depth. Site effect and aleatory variability models are shared between event types. Regionalized GMM components include the model constant (that controls ground motion amplitude), anelastic attenuation, magnitude-scaling break point, linear site response, and sediment depth terms. We develop models for the aleatory between-event variability ( τ ) "> ( τ ) , within-event variability ( ϕ ) "> ( ϕ ) , single-station within-event variability ( ϕ SS ) "> ( ϕ SS ) , and site-to-site variability ( ϕ S 2 S ) "> ( ϕ S 2 S ) . Ergodic analyses should use the median GMM and aleatory variability computed using the between-event and within-event variability models. An analysis incorporating non-ergodic site response should use the median GMM at the reference shear-wave velocity condition, a site-specific site response model, and aleatory variability computed using the between-event and single-station within-event variability models. Epistemic uncertainty in the median model is represented by standard deviations on the regional model constants, which facilitates scaled-backbone representations of model uncertainty in hazard analyses.

Earthquake Spectra↗

USGS2016

A decision framework for identifying models to estimate forest ecosystem services gains from restoration

Restoring degraded forests and agricultural lands has become a global conservation priority. A growing number of tools can quantify ecosystem service tradeoffs associated with forest restoration. This evolving “tools landscape” presents a dilemma: more tools are available, but selecting appropriate tools has become more challenging. We present a Restoration Ecosystem Service Tool Selector (RESTS) framework that describes key characteristics of 13 ecosystem service assessment tools. Analysts enter information about their decision context, services to be analyzed, and desired outputs. Tools are filtered and presented based on five evaluative criteria: scalability, cost, time requirements, handling of uncertainty, and applicability to benefit-cost analysis. RESTS uses a spreadsheet interface but a web-based interface is planned. Given the rapid evolution of ecosystem services science, RESTS provides an adaptable framework to guide forest restoration decision makers toward tools that can help quantify ecosystem services in support of restoration.

Forest Ecosystems↗

USGS1994

Large-explosive source, wide-recording aperture, seismic profiling on the Columbia Plateau, Washington

Clear subsurface seismic images have been obtained at low cost on the Columbia Plateau, Washington. The Columbia Plateau is perhaps the most notorious of all 'bad-data' areas because large impedance contrasts in surface flood basalts severely degrade the seismic wavefield. This degradation was mitigated in this study via a large-explosive source, wide-recording aperture shooting method.The shooting method emphasizes the wide-angle portion of the wavefield, where Fermat's principle guarantees reverberation will not interfere with the seismic manifestations of crucial geologic interfaces. The basalt diving wave, normally discarded in standard common midpoint (CMP) seismic profiling, can be used to image basalt velocity structure via traveltime inversion. Maximum depth-penetration of the diving wave tightly constrains basalt-sediment interface depth. An arrival observed only at shot-receiver offsets greater than 15 km can be used to determine the velocity and geometry of basement via simultaneous inversion.The results from this study suggest that previous geologic hypotheses and hydrocarbon play concepts for the Columbia Plateau may have been in error.

Geophysics↗

USGS1994

An image of the Columbia Plateau from inversion of high‐resolution seismic data

We use a method of traveltime inversion of high‐resolution seismic data to provide the first reliable images of internal details of the Columbia River Basalt Group (CRBG), the subsurface basalt/sediment interface, and the deeper sediment/basement interface. Velocity structure within the basalts, delineated on the order of 1 km horizontally and 0.2 km vertically, is constrained to within ±0.1 km/s for most of the seismic profile. Over 5000 observed traveltimes fit our model with an rms error of 0.018 s. The maximum depth of penetration of the basalt diving waves (truncated by underlying low‐velocity sediments) provides a reliable estimate of the depth to the base of the basalt, which agrees with well‐log measurements to within 0.05 km (165 ft). We use image blurring, calculated from the resolution matrix, to estimate the aspect ratio of imaged velocity anomaly widths to true widths for velocity features within the basalt. From our calculations of image blurring, we interpret low velocity zones (LVZ) within the basalts at Boylston Mountain and the Whiskey Dick anticline to have widths of 4.5 and 3 km, respectively, within the upper 1.5 km of the model. At greater depth, the widths of these imaged LVZs thin to approximately 2 km or less. We interpret these linear, subparallel, low‐velocity zones imaged adjacent to anticlines of the Yakima Fold Belt to be brecciated fault zones. These fault zones dip to the south at angles between 15 to 45 degrees.

Geophysics↗