SEARCH · Search USGS

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 415 records · Page 23Linked to original sources

USGS2012

Comparison of soil thickness in a zero-order basin in the Oregon Coast Range using a soil probe and electrical resistivity tomography

Accurate estimation of the soil thickness distribution in steepland drainage basins is essential for understanding ecosystem and subsurface response to infiltration. One important aspect of this characterization is assessing the heavy and antecedent rainfall conditions that lead to shallow landsliding. In this paper, we investigate the direct current (DC) resistivity method as a tool for quickly estimating soil thickness over a steep (33–40°) zero-order basin in the Oregon Coast Range, a landslide prone region. Point measurements throughout the basin showed bedrock depths between 0.55 and 3.2 m. Resistivity of soil and bedrock samples collected from the site was measured for degrees of saturation between 40 and 92%. Resistivity of the soil was typically higher than that of the bedrock for degrees of saturation lower than 70%. Results from the laboratory measurements and point-depth measurements were used in a numerical model to evaluate the resistivity contrast at the soil-bedrock interface. A decreasing-with-depth resistivity contrast was apparent at the interface in the modeling results. At the field site, three transects were surveyed where coincident ground truth measurements of bedrock depth were available, to test the accuracy of the method. The same decreasing-with-depth resistivity trend that was apparent in the model was also present in the survey data. The resistivity contour of between 1,000 and 2,000 Ωm that marked the top of the contrast was our interpreted bedrock depth in the survey data. Kriged depth-to-bedrock maps were created from both the field-measured ground truth obtained with a soil probe and interpreted depths from the resistivity tomography, and these were compared for accuracy graphically. Depths were interpolated as far as 16.5 m laterally from the resistivity survey lines with root mean squared error (RMSE) = 27 cm between the measured and interpreted depth at those locations. Using several transects and analysis of the subsurface material properties, the direct current (DC) resistivity method is shown to be able to delineate bedrock depth trends within the drainage basin.

Oregon↗

USGS2012

Simulating the effects of ground-water withdrawals on streamflow in a precipitation-runoff model

Precipitation-runoff models are used to assess the effects of water use and management alternatives on streamflow. Often, ground-water withdrawals are a major water-use component that affect streamflow, but the ability of surface-water models to simulate ground-water withdrawals is limited. As part of a Hydrologic Simulation Program-FORTRAN (HSPF) precipitation-runoff model developed to analyze the effect of ground-water and surface-water withdrawals on streamflow in the Ipswich River in northeastern Massachusetts, an analytical technique (STRMDEPL) was developed for calculating the effects of pumped wells on streamflow. STRMDEPL is a FORTRAN program based on two analytical solutions that solve equations for ground-water flow to a well completed in a semi-infinite, homogeneous, and isotropic aquifer in direct hydraulic connection to a fully penetrating stream. One analytical method calculates unimpeded flow at the stream-aquifer boundary and the other method calculates the resistance to flow caused by semipervious streambed and streambank material. The principle of superposition is used with these analytical equations to calculate time-varying streamflow depletions due to daily pumping. The HSPF model can readily incorporate streamflow depletions caused by a well or surface-water withdrawal, or by multiple wells or surface-water withdrawals, or both, as a combined time-varying outflow demand from affected channel reaches. These demands are stored as a time series in the Watershed Data Management (WDM) file. This time-series data is read into the model as an external source used to specify flow from the first outflow gate in the reach where these withdrawals are located. Although the STRMDEPL program can be run independently of the HSPF model, an extension was developed to run this program within GenScn, a scenario generator and graphical user interface developed for use with the HSPF model. This extension requires that actual pumping rates for each well be stored in a unique WDM dataset identified by an attribute that associates each well with the model reach from which water is withdrawn. Other attributes identify the type and characteristics of the data. The interface allows users to easily add new pumping wells, delete exiting pumping wells, or change properties of the simulated aquifer or well. Development of this application enhanced the ability of the HSPF model to simulate complex water-use conditions in the Ipswich River Basin. The STRMDEPL program and the GenScn extension provide a valuable tool for water managers to evaluate the effects of pumped wells on streamflow and to test alternative water-use scenarios. Copyright ASCE 2004.

Conference Paper↗

USGS2022

The evolution of rock friction is more sensitive to slip than elapsed time, even at near-zero slip rates

For many decades, frictional strength increase at low slip rates has been ascribed to time-dependent contact-area growth across the sliding interface. As a result, phenomenological models that correctly predict contact-area growth, as observed in laboratory experiments, have also been widely assumed to be appropriate descriptors of frictional strength evolution. We present experiments that impose more than 5-orders-of-magnitude slip-rate reductions on granite to show that frictional strength evolution in these rocks unequivocally refutes such models. Instead, the data suggest that, even at subnanometric slip rates, frictional strength dominantly evolves with accrued slip. This remarkable slip-sensitivity of friction requires changes of intrinsic strength of the interface with slip that are absent from popular conceptual models of friction at the microscopic contact scale.

Proceedings of the National Academy of Sciences↗

USGS2018

The map as knowledge base

This paper examines the concept and implementation of a map as a knowledge base. A map as a knowledge base means that the visual map is not only the descriptive compilation of data and design principles, but also involves a compilation of semantic propositions and logical predicates that create a body of knowledge organized as a map. The digital product of a map as knowledge base can be interpreted by machines, as well as humans, and can provide access to the knowledge base through interfaces to select features and other information from the map. The design of maps as a knowledge base involves technical approaches and a system architecture to support a knowledge base. This paper clarifies how a map as a knowledge base differs from earlier map theory models by investigating the knowledge-based concepts of implementation through logical modelling, a knowledge repository, user interfaces for information access, and cartographic visualization. The paper ends with proof of concepts for two types of cartographic data query.

International Journal of Cartography↗

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↗

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↗