USGS Science⌕ Search

SEARCH · USGS Science

Results for “earthquake hazards”

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 1,711 records · Page 95Linked to original sources

Regression models for estimating coseismic landslide displacement

Newmark's sliding-block model is widely used to estimate coseismic slope performance. Early efforts to develop simple regression models to estimate Newmark displacement were based on analysis of the small number of strong-motion records then available. The current availability of a much larger set of strong-motion records dictates that these regression equations be updated. Regression equations were generated using data derived from a collection of 2270 strong-motion records from 30 worldwide earthquakes. The regression equations predict Newmark displacement in terms of (1) critical acceleration ratio, (2) critical acceleration ratio and earthquake magnitude, (3) Arias intensity and critical acceleration, and (4) Arias intensity and critical acceleration ratio. These equations are well constrained and fit the data well (71% < R2 < 88%), but they have standard deviations of about 0.5 log units, such that the range defined by the mean ?? one standard deviation spans about an order of magnitude. These regression models, therefore, are not recommended for use in site-specific design, but rather for regional-scale seismic landslide hazard mapping or for rapid preliminary screening of sites. ?? 2007 Elsevier B.V. All rights reserved.

Engineering Geology↗

Map and database of Quaternary faults and folds in Colombia and its offshore regions

As part of the International Lithosphere Program’s “World Map of Major Active Faults,” the U.S. Geological Survey (USGS) is assisting in the compilation of a series of digital maps of Quaternary faults and folds in Western Hemisphere countries. The maps show the locations, ages, and activity rates of major earthquake-related features such as faults and fault-related folds. They are accompanied by databases that describe these features and document current information on their activity in the Quaternary. Top date, the project has published fault and fold maps for Costa Rica (Montero and others, 1998), Panama (Cowan and others, 1998), Venezuela (Audemard and others, 2000), Bolovia/Chile (Lavenu, and others, 2000), and Argentina (Costa and others, 2000). The project is a key part of the Global Seismic Hazards Assessment Program (ILP Project II-0) for the International Decade for Natural Hazard Disaster Reduction.

Open-File Report↗

The thermal regime in the resurgent dome of Long Valley Caldera, California: Inferences from precision temperature logs in deep wells

Long Valley Caldera in eastern California formed 0.76Ma ago in a cataclysmic eruption that resulted in the deposition of 600km3 of Bishop Tuff. The total current heat flow from the caldera floor is estimated to be ~290MW, and a geothermal power plant in Casa Diablo on the flanks of the resurgent dome (RD) generates ~40MWe. The RD in the center of the caldera was uplifted by ~80cm between 1980 and 1999 and was explained by most models as a response to magma intrusion into the shallow crust. This unrest has led to extensive research on geothermal resources and volcanic hazards in the caldera. Here we present results from precise, high-resolution, temperature-depth profiles in five deep boreholes (327-1,158m) on the RD to assess its thermal state, and more specifically 1) to provide bounds on the advective heat transport as a guide for future geothermal exploration, 2) to provide constraints on the occurrence of magma at shallow crustal depths, and 3) to provide a baseline for future transient thermal phenomena in response to large earthquakes, volcanic activity, or geothermal production. The temperature profiles display substantial non-linearity within each profile and variability between the different profiles. All profiles display significant temperature reversals with depth and temperature gradients <50??C/km at their bottom. The maximum temperature in the individual boreholes ranges between 124.7??C and 129.5??C and bottom hole temperatures range between 99.4??C and 129.5??C. The high-temperature units in the three Fumarole Valley boreholes are at the approximate same elevation as the high-temperature unit in borehole M-1 in Casa Diablo indicating lateral or sub-lateral hydrothermal flow through the resurgent dome. Small differences in temperature between measurements in consecutive years in three of the wells suggest slow cooling of the shallow hydrothermal flow system. By matching theoretical curves to segments of the measured temperature profiles, we calculate horizontal groundwater velocities in the hydrothermal flow unit under the RD that range from 1.9 to 2.8m/yr, which corresponds to a maximum power flowing through the RD of 3-4MW. The relatively low temperatures and large isothermal segments at the bottom of the temperature profiles are inconsistent with the presence of magma at shallow crustal levels. ?? 2010.

Journal of Volcanology and Geothermal Research↗

Monitoring a supervolcano in repose: Heat and volatile flux at the yellostone caldera

Although giant calderas ("supervolcanoes") may slumber for tens of thousands of years between eruptions, their abundant earthquakes and crustal deformation reveal the potential for future upheaval. Any eventual supereruption could devastate global human populations, so these systems must be carefully scrutinized. Insight into dormant but restless calderas can be gained by monitoring their output of heat and gas. At Yellowstone, the large thermal and CO2 fluxes require massive input of basaltic magma, which continues to invade the lower to mid-crust, sustains the overlying high-silica magma reservoir, and may result in volcanic hazard for millennia to come. The high flux of CO2 may contribute to the measured deformation of the caldera floor and can also modify the pressure, thermal, and chemical signals emitted from the magma. In order to recognize precursors to eruption, we must scrutinize the varied signals emerging from restless calderas with the goal of discriminating magmatic, hydrothermal, and hybrid phenomena.

Elements↗

Improving near‐real‐time coseismic landslide models: Lessons learned from the 2016 Kaikōura, New Zealand, earthquake

The U.S. Geological Survey (USGS) is developing near‐real‐time global earthquake‐triggered‐landslide products to augment the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) system. The 14 November 2016 M w "> M w Mw 7.8 Kaikōura, New Zealand, earthquake provided a test case for evaluating the performance and near‐real‐time response applicability of three published global seismically induced landslide models. All three models obtain shaking estimates from the USGS ShakeMap, which is updated and sometimes changes significantly in the hours to days after an earthquake. The Kaikōura earthquake is a particularly valuable event that helps us better understand how changes to the ShakeMap affect the landslide models because the ShakeMap evolved significantly over several weeks as multifault rupture and seismic data were incorporated. We used the detailed landslide inventory available for this event for qualitative landslide model evaluation. We found that once a point source was replaced with an approximate rupture extent in ShakeMap, the landslide models were all successful at roughly identifying the area of highest hazard. This is notable, given that the models are relatively simple, coarse in resolution, and are based solely on input proxies that are globally available. However, all of the models dramatically overpredicted the hazard level, which indicates that improvements can be made. Subsequent updates to the ShakeMap resulted in improvements to model performance by some metrics and declining performance by others. In all cases, details of the ShakeMap strongly controlled the spatial pattern, even when those details were erroneous, such as the inclusion of a fault segment that did not rupture. If maps of landslide hazard are to be used effectively for rapid response, then we need to understand and clearly communicate the control that ShakeMap has over the models and how that typically evolves with time and is (or is not) reflected in reported uncertainties.

Bulletin of the Seismological Society of America↗

The need for the International Decade of Natural Hazard Reduction

Natural disasters affect countires large and small, rich or poor, whatever their political persuasion. The toll exacted by natural calamities each year drains the human and economic resources of every nation and stands as one of the formidable barriers to national, regional, and world development. Over the last 20 years, natural disasters have killed nearly 3 million people and disrupted the lives of over 800 million others. In 2 years there were more than 50 serious natural disasters, including landslides in Italy, France, and Colombia; a typhoon in Korea; wildfires in China and the United State; a windstorm in England; grasshopper plagues in Africa's horns and the Sahel; tornadoes in Canada; devastating earthquakes in Soviet Armenia and Tadzhikistan; infestations in Africa; landslides in Brazil; and tornadoes in the United States.

Earthquakes & Volcanoes (USGS)↗

Event sedimentation in low-latitude deep-water carbonate basins, Anegada passage, northeast Caribbean

The Virgin Islands and Whiting basins in the Northeast Caribbean are deep, structurally controlled depocentres partially bound by shallow-water carbonate platforms. Closed basins such as these are thought to document earthquake and hurricane events through the accumulation of event layers such as debris flow and turbidity current deposits and the internal deformation of deposited material. Event layers in the Virgin Islands and Whiting basins are predominantly thin and discontinuous, containing varying amounts of reef- and slope-derived material. Three turbidites/sandy intervals in the upper 2 m of sediment in the eastern Virgin Islands Basin were deposited between ca. 2000 and 13 600 years ago, but do not extend across the basin. In the central and western Virgin Islands Basin, a structureless clay-rich interval is interpreted to be a unifite. Within the Whiting Basin, several discontinuous turbidites and other sand-rich intervals are primarily deposited in base of slope fans. The youngest of these turbidites is ca. 2600 years old. Sediment accumulation in these basins is low (&minus;1) for basin adjacent to carbonate platform, possibly due to limited sediment input during highstand sea-level conditions, sediment trapping and/or cohesive basin walls. We find no evidence of recent sediment transport (turbidites or debris flows) or sediment deformation that can be attributed to the ca. M7.2 1867 Virgin Islands earthquake whose epicentre was located on the north wall of the Virgin Islands Basin or to recent hurricanes that have impacted the region. The lack of significant appreciable pebble or greater size carbonate material in any of the available cores suggests that submarine landslide and basin-wide blocky debris flows have not been a significant mechanism of basin margin modification in the last several thousand years. Thus, basins such as those described here may be poor recorders of past natural hazards, but may provide a long-term record of past oceanographic conditions in ocean passages.

Basin Research↗

Sustained long-period seismicity at Shishaldin Volcano, Alaska

From September 1999 through April 2004, Shishaldin Volcano, Aleutian Islands, Alaska, exhibited a continuous and extremely high level of background seismicity. This activity consisted of many hundreds to thousands of long-period (LP; 1–2 Hz) earthquakes per day, recorded by a 6-station monitoring network around Shishaldin. The LP events originate beneath the summit at shallow depths (0–3 km). Volcano tectonic events and tremor have rarely been observed in the summit region. Such a high rate of LP events with no eruption suggests that a steady state process has been occurring ever since Shishaldin last erupted in April–May 1999. Following the eruption, the only other signs of volcanic unrest have been occasional weak thermal anomalies and an omnipresent puffing volcanic plume. The LP waveforms are nearly identical for time spans of days to months, but vary over longer time scales. The observations imply that the spatially close source processes are repeating, stable and non-destructive. Event sizes vary, but the rate of occurrence remains roughly constant. The events range from magnitude ∼0.1 to 1.8, with most events having magnitudes <1.0. The observations suggest that the conduit system is open and capable of releasing a large amount of energy, approximately equivalent to at least one magnitude 1.8–2.6 earthquake per day. The rate of observed puffs (1 per minute) in the steam plume is similar to the typical seismic rates, suggesting that the LP events are directly related to degassing processes. However, the source mechanism, capable of producing one LP event about every 0.5–5 min, is still poorly understood. Shishaldin's seismicity is unusual in its sustained high rate of LP events without accompanying eruptive activity. Every indication is that the high rate of seismicity will continue without reflecting a hazardous state. Sealing of the conduit and/or change in gas flux, however, would be expected to change Shishaldin's behavior.

Alaska↗

Searching for seismic precursors - The Barry Landslide hazard clean up

The Barry Landslide, located in Barry Arm of Prince William Sound, Alaska, poses a major hazard due to its steep, unstable slopes and the potential for a massive landslide-generated tsunami. With an estimated volume of 500–700 million cubic meters, the Barry Landslide could trigger highly destructive waves. In this study, we focus on seismic signals from the Barry Landslide, which are critical for providing timely tsunami warnings. Since the summer of 2020, the region has been instrumented to monitor the landslide, but the seismic record is complicated by the presence of nearby glaciers and frequent regional earthquakes. Among these signals, we analyze a specific class of short-duration, high-frequency seismic events that exhibit strong seasonal variability, increasing in rate from late summer to mid-winter before ceasing abruptly in late winter or early spring. Our analysis suggests that the source of these signals is likely near or beneath Cascade Glacier, adjacent to the landslide, rather than within the landslide mass itself. We apply detection algorithms to construct a time history for this signal type, which we then compare with environmental factors like precipitation, temperature, and slope displacement data from ground-based radar and remote sensing. Correlations indicate that these seismic events may be driven by seasonal hydrological changes, particularly the freeze-up of subglacial water pathways. While these events are not directly linked to landslide motion, they serve as indirect markers of subsurface hydrological conditions that influence slope stability. Our findings highlight the complex interplay between glaciers, groundwater, and landslide dynamics, emphasizing the need for multi-parameter monitoring to assess evolving geohazards in the region.

Alaska↗

Source character of microseismicity in the San Francisco Bay block, California, and implications for seismic hazard

We examine relocated seismicity within a 30-km-wide crustal block containing San Francisco Bay and bounded by two major right-lateral strike-slip fault systems, the Hayward and San Andreas faults, to determine seismicity distribution, source character, and possible relationship to proposed faults. Well-located low-level seismicity ( M d ≦ 3.0) has occurred persistently within this block throughout the recording interval (1969 to 1995), with the highest levels of activity occurring along or directly adjacent to (within ∼5 km) the bounding faults and falling off toward the long axis of the bay. The total seismic moment release within the interior of the Bay block since 1969 is equivalent to one M L 3.8 earthquake, one to two orders of magnitude lower than activity along and within 5 km of the bounding faults. Focal depths of reliably located events within the Bay block are generally less than 13 km with most seismicity in the depth range of 7 to 12 km, similar to focal depths along both the adjacent portions of the San Andreas and Hayward faults. Focal mechanisms for M d 2 to 3 events within the Bay block mimic focal mechanisms along the adjacent San Andreas fault zone and in the East Bay, suggesting that Bay block is responding to a similar regional stress field. Two potential seismic source zones have been suggested within the Bay block. Our hypocentral depths and focal mechanisms suggest that a proposed subhorizontal detachment fault 15 to 18 km beneath the Bay is not seismically active. Several large-scale linear NW-trending aeromagnetic anomalies within the Bay block were previously suggested to represent large through-going subvertical fault zones. The two largest earthquakes (both M d 3.0) in the Bay block since 1969 occur near two of these large-scale linear aeromagnetic anomalies; both have subvertical nodal planes with right-lateral slip subparallel to the magnetic anomalies, suggesting that structures related to the anomalies may be capable of brittle failure. Geodetic, focal mechanism and seismicity data all suggest the Bay block is responding elastically to the same regional stresses affecting the bounding faults; however, continuous Holocene reflectors across the proposed fault zones suggest that if the magnetic anomalies represent basement fault zones, then these faults must have recurrence times one to several orders of magnitude longer than on the bounding faults.

California↗

Effect of structural heterogeneity and slip distribution on coseismic vertical displacement from rupture on the Seattle Fault

Workshops in 2001 and 2002 were convened to determine critical issues in the development of tsunami inundation maps for the Puget Sound region. The Tsunami Inundation Mapping Effort (TIME) is conducted under the multi-agency National Tsunami Hazard Mitigation Program (NTHMP). The Puget Sound Tsunami/Landslide Workshop in 2001 focused on integrated tsunami research involving a wide range of research studies and tsunami hazard mitigation issues. The 2002 Puget Sound Tsunami Sources workshop (González et al., 2003) made specific recommendations for tsunami source modeling and improving our state of knowledge for sources in the Puget Sound region. One of the recommendations stated in González et al. (2003) is "Develop methods to assess the sensitivity of coastal areas to tsunami inundation, based on multiple simulations that reflect the possible range of variations in the source parameters." Tsunami inundation models rely heavily on the imposed initial conditions which, for an earthquake source, is the coseismic vertical displacement field. For example, Koshimura et al. (2002) use the geologic uplift observations (Buknam et al., 1992) to constrain the slip distribution for the event that occurred 1100 years ago, resulting in an average slip of 3.7 m and a magnitude of 7.6. Walsh et al. (2003) develop a tsunami inundation map for Elliot Bay based on a M 7.3 earthquake and the geologic uplift observations from the 1100 y.b.p. event as in Koshimura et al. (2002), though they use a constant fault dip of 60° rather than different dips for deep and shallow segments. The objective of this report is to examine how coseismic vertical displacement from a smaller M 6.5 Seattle Fault earthquake (as in Hartzell et al., 2002) is affected by structural heterogeneity and different slip distribution patterns.

Washington↗

A ground‐motion prediction model for shallow crustal earthquakes in Greece

Using a recently completed database of uniformly processed strong‐motion data recorded in Greece, we derive a ground‐motion prediction model (GMPM) for horizontal‐component peak ground velocity, peak ground acceleration, and 5% damped pseudoacceleration response spectra, at 105 periods ranging from 0.01 to 10 s. The equations were developed by modifying a global GMPM, to account for more rapid attenuation and weaker magnitude scaling in the Greek ground motions than in the global GMPM. Our GMPM is calibrated using the Greek data for distances up to 300 km, magnitudes from 4.0 to 7.0, and time‐averaged 30 m shear‐wave velocities from 150 to 1200 &#x2009;&#x2009; m / s "> 1200 m / s 1200 m/s ⁠ . The GMPM has important attributes for hazard applications including magnitude scaling that extends the range of applicability to M "> M M 8.0 and nonlinear site response. These features are possible because they are well constrained by data in the global GMPM from which our model is derived. An interesting feature of the Greek data, also observed previously in studies of mid‐magnitude events (6.1–6.5) in Italy, is that they are substantially overpredicted by the global GMPM, which may be a repeatable regional feature, but may also be influenced by soil–structure interaction. This bias is an important source of epistemic uncertainty that should be considered in hazard analysis.

Bulletin of the Seismological Society of America↗

Book review: Three great tsunamis: Lisbon (1755), Sumatra-Andaman (2004), and Japan (2011)

&ldquo;Three Great Tsunamis: Lisbon (1755), Sumatra&ndash;Andaman (2004), and Japan (2011)&rdquo; is published in Springer&rsquo;s new series SpringerBriefs. According to Springer&rsquo;s website, the SpringBriefs volumes are intended to provide &ldquo;concise summaries of cutting-edge research and practical applications across a wide spectrum of fields&rdquo;. Among the several categories considered for SpringerBriefs are in-depth case studies, for which this volume is most closely aligned. This 89-page volume comprises six chapters. Chapter 1 describes the fundamentals of tsunamis at a basic level, including a selective list of major historical tsunamis. Chapter 2 describes the 1755 Lisbon earthquake and tsunami, although no tectonic background is given for this event, as it is for the other two events considered in this volume. Included in this chapter, however, is a brief but fascinating account of how this earthquake affected scientific and philosophical thought during the Age of Enlightenment. The authors&rsquo; description of the 2004 Sumatra&ndash;Andaman earthquake and tsunami in Chapter 3 includes instrumental and post-tsunami measurements of the tsunami, ground shaking and subsidence, and other effects such as mud volcanoes. Primarily, tsunami observations are described for the 2011 Japan event, including a brief summary of the nuclear disaster in the Fukushima region. The fifth chapter is on tsunami warning systems, including new and emerging systems around the world in the twenty-first century. The final chapter provides perspectives on why these events were so destructive and how future hazards may be mitigated. Review info: Three Great Tsunamis: Lisbon (1755), Sumatra&ndash;Andaman (2004), and Japan (2011). By Harsh K. Gupta and Vineet K. Gahalaut, 2013. ISBN: 978-9400765757

Pure and Applied Geophysics↗

The rupture process of the 2018 Mw 6.9 Hawaiʻi earthquake as imaged by a genetic algorithm-based back-projection technique

An episode of unrest began at Kīlauea in April 2018 that produced both significant volcanic output and high rates of seismicity, including a M w 6.9 earthquake on 4 May 2018. In this study, we image the rupture process of this earthquake using a genetic algorithm-based back-projection technique. The dominant feature of the earthquake is a slowly propagating western rupture, which shares similar characteristics with the region's largest recorded event in 1975 ( M w 7.7). The location of this western segment suggests that small asperities on this section of the décollement that frequently fail as slow slip events may achieve seismic slip rates when rupture is initiated on adjacent sections of the fault. Given the interaction between volcanic and seismic activity in this region, imaging the rupture properties of these events can improve our understanding of future geologic hazards in this region.

Hawaii↗

Earthquake magnitude distributions on northern Caribbean faults from combinatorial optimization models

On-fault earthquake magnitude distributions are calculated for northern Caribbean faults using estimates of fault slip and regional seismicity parameters. Integer programming, a combinatorial optimization method, is used to determine the optimal spatial arrangement of earthquakes sampled from a truncated Gutenberg-Richter distribution that minimizes the global misfit in slip rates on a complex fault system. Slip rates and their uncertainty on major faults are derived from a previously published GPS block model for the region, with fault traces determined from offshore geophysical mapping and previously published onshore studies. The optimal spatial arrangement of the sampled earthquakes is compared with the 500-year history of earthquake observations. Rupture segmentation of the subduction interface along the Hispaniola-Puerto Rico Trench (PRT) fault and seismic coupling on the PRT fault appear to exert the primary control over this spatial arrangement. Introducing a rupture barrier for the Hispaniola-PRT fault northwest of Mona Passage, based on geophysical and seismicity observations, and assigning a low slip rate of 2 mm/yr on the PRT fault are most consistent with historical earthquakes in the region. The addition of low slip-rate secondary faults as well as segmentation of the Hispaniola and Septentrional strike-slip fault improves the consistency with historical seismicity. An important observation from the modeling is that varying the slip rate on the PRT fault and different segmentation scenarios result in significant changes to the optimal magnitude distribution on faults farther away. In general, optimal on-fault magnitude distributions are more complex and inter-dependent than is typically assumed in probabilistic seismic hazard analysis and probabilistic tsunami hazard analysis.

JGR Solid Earth↗

Earthquake impact scale

With the advent of the USGS prompt assessment of global earthquakes for response (PAGER) system, which rapidly assesses earthquake impacts, U.S. and international earthquake responders are reconsidering their automatic alert and activation levels and response procedures. To help facilitate rapid and appropriate earthquake response, an Earthquake Impact Scale (EIS) is proposed on the basis of two complementary criteria. On the basis of the estimated cost of damage, one is most suitable for domestic events; the other, on the basis of estimated ranges of fatalities, is generally more appropriate for global events, particularly in developing countries. Simple thresholds, derived from the systematic analysis of past earthquake impact and associated response levels, are quite effective in communicating predicted impact and response needed after an event through alerts of green (little or no impact), yellow (regional impact and response), orange (national-scale impact and response), and red (international response). Corresponding fatality thresholds for yellow, orange, and red alert levels are 1, 100, and 1,000, respectively. For damage impact, yellow, orange, and red thresholds are triggered by estimated losses reaching \$1M, \$100M, and \$1B, respectively. The rationale for a dual approach to earthquake alerting stems from the recognition that relatively high fatalities, injuries, and homelessness predominate in countries in which local building practices typically lend themselves to high collapse and casualty rates, and these impacts lend to prioritization for international response. In contrast, financial and overall societal impacts often trigger the level of response in regions or countries in which prevalent earthquake resistant construction practices greatly reduce building collapse and resulting fatalities. Any newly devised alert, whether economic- or casualty-based, should be intuitive and consistent with established lexicons and procedures. Useful alerts should also be both specific (although allowably uncertain) and actionable. In this analysis, an attempt is made at both simple and intuitive color-coded alerting criteria; yet the necessary uncertainty measures by which one can gauge the likelihood for the alert to be over- or underestimated are preserved. The essence of the proposed impact scale and alerting is that actionable loss information is now available in the immediate aftermath of significant earthquakes worldwide on the basis of quantifiable loss estimates. Utilizing EIS, PAGER’s rapid loss estimates can adequately recommend alert levels and suggest appropriate response protocols, despite the uncertainties; demanding or awaiting observations or loss estimates with a high level of accuracy may increase the losses.

Natural Hazards Review↗

Holocene surface ruptures on the Salinas Fault and southeastern Great Southern Puerto Rico Fault Zone, South Coastal Plain of Puerto Rico

We analyzed light detection and ranging (lidar) data and aerial photography to locate active faults near the south coast of Puerto Rico and excavated paleoseismic trenches across the Salinas fault and the Great Southern Puerto Rico fault zone (GSPRFZ). We document evidence for two Holocene surface‐rupturing earthquakes along both faults. Two earthquakes on the Salinas fault occurred after the deposition of sediments that are 7400–10,400 yrs old. No quantitative ages constrain the timing of the two earthquakes on the southeast GSPRFZ, but we interpret both events to have occurred during the Holocene or latest Pleistocene, based on the similarities in the characteristics of the faulted sediment and the development of soils exposed in both trenches. Stratigraphic and geomorphic evidence suggests components of both vertical and lateral slip on both faults. These results show that onshore active faults in Puerto Rico are more common than previously recognized and highlight the need for additional study to search for Holocene active faults and for an updated seismic hazard analysis for the island.

Puerto Rico↗

The Wallula fault and tectonic framework of south-central Washington, as interpreted from magnetic and gravity anomalies

The Yakima fold and thrust belt (YFTB) in central Washington has accommodated regional, mostly north-directed, deformation of the Cascadia backarc since prior to emplacement of Miocene flood basalt of the Columbia River Basalt Group (CRBG). The YFTB consists of two structural domains. Northern folds of the YFTB strike eastward and terminate at the western margin of a 20-mGal negative gravity anomaly, the Pasco gravity low, straddling the North American continental margin. Southern folds of the YFTB strike southeastward, form part of the Olympic–Wallowa lineament (OWL), and pass south of the Pasco gravity low as the Wallula fault zone. An upper crustal model based on gravity and magnetic anomalies suggests that the Pasco gravity low is caused in part by an 8-km-deep Tertiary basin, the Pasco sub-basin, abutting the continental margin and concealed beneath CRBG. The Pasco sub-basin is crossed by north-northwest-striking magnetic anomalies caused by dikes of the 8.5 Ma Ice Harbor Member of the CRBG. At their northern end, dikes connect with the eastern terminus of the Saddle Mountains thrust of the YFTB. At their southern end, dikes are disrupted by the Wallula fault zone. The episode of NE–SW extension that promoted Ice Harbor dike injection apparently involved strike-slip displacement on the Saddle Mountains and Wallula faults. The amount of lateral shear on the OWL impacts the level of seismic hazard in the Cascadia region. Ice Harbor dikes, as mapped with aeromagnetic data, are dextrally offset by the Wallula fault zone a total of 6.9 km. Assuming that dike offsets are tectonic in origin, the Wallula fault zone has experienced an average dextral shear of 0.8 mm/y since dike emplacement 8.5 Ma, consistent with right-lateral stream offsets observed at other locations along the OWL. Southeastward, the Wallula fault transfers strain to the north-striking Hite fault, the possible location of the M 5.7 Milton-Freewater earthquake in 1936.

Washington↗