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,351 records · Page 75Linked to original sources
USGS2022

Western U.S. deformation models for the 2023 update to the U.S. National Seismic Hazard Model

This report describes geodetic and geologic information used to constrain deformation models of the 2023 update to the National Seismic Hazard Model (NSHM), a set of deformation models to interpret these data, and their implications for earthquake rates in the western United States. Recent updates provide a much larger data set of Global Positioning System crustal velocities than used in the 2014 NSHM, as well as hundreds of new faults considered as active sources for the 2023 NSHM. These data are interpreted by four geodetic models of deformation that estimate fault slip rates and their uncertainties together with off‐fault moment release rates. Key innovations in the 2023 NSHM relative to past practice include (1) the addition of two new (in addition to two existing) deformation models, (2) the revision and expansion of the geologic slip rate database, (3) accounting for fault creep through development of a creep‐rate model that is employed by the four deformation models, and (4) accounting for time‐dependent earthquake‐cycle effects through development of viscoelastic models of the earthquake cycle along the San Andreas fault and the Cascadia subduction zone. The effort includes development of a geologic deformation model that complements the four geodetic models. The current deformation models provide a new assessment of outstanding discrepancies between geologic and geodetic slip rates, at the same time highlighting the need for both geologic and geodetic slip rates to robustly inform the earthquake rate model.

Seismological Research Letters↗
USGS1996

Faults, lineaments, and earthquake epicenters digital map of the Pahute Mesa 30' x 60' quadrangle, Nevada

This map database, identified as Faults, lineaments, and earthquake epicenters digital map of the Pahute Mesa 30' X 60' quadrangle, Nevada, has been approved for release and publication by the Director of the USGS. Although this database has been subjected to rigorous review and is substantially complete, the USGS reserves the right to revise the data pursuant to further analysis and review. Furthermore, it is released on condition that neither the USGS nor the United States Government may be held liable for any damages resulting from its authorized or unauthorized use. This digital map compilation incorporates fault, air photo lineament, and earthquake epicenter data from within the Pahute Mesa 30' by 60' quadrangle, southern Nye County, Nevada (fig. 1). The compilation contributes to the U.S. Department of Energy's Yucca Mountain Project, established to determine whether or not the Yucca Mountain site is suitable for the disposal of high-level nuclear waste. Studies of local and regional faulting and earthquake activity, including the features depicted in this compilation, are carried out to help characterize seismic hazards and tectonic processes that may be relevant to the future stability of Yucca Mountain. The Yucca Mountain site is located in the central part of the Beatty 30' by 60' quadrangle approximately 15 km south of the south edge of the Pahute Mesa quadrangle (fig. 1). The U.S. Geological Survey participates in studies of the Yucca Mountain site under Interagency Agreement DE-AI08-78ET44802. The map compilation is only available on line as a digital database in ARC/INFO ASCII (Generate) and export formats. The database can be downloaded via 'anonymous ftp' from a USGS system named greenwood.cr.usgs.gov (136.177.48.5). The files are located in a directory named /pub/open-file-reports/ofr-96-0262. This directory contains a text document named 'README.1 ST' that contains database technical and explanatory documentation, including instructions for uncompressing the bundled (tar) file. In displaying the compilation it is important to note that the map data set is considered accurate when depicted at a scale of about 1:100,000; displaying the compilation at scales significantly larger than this may result in distortions and (or) mislocations of the data.

Nevada↗
USGS2019

The 2018 update of the US National Seismic Hazard Model: Overview of model and implications

During 2017–2018, the National Seismic Hazard Model for the conterminous United States was updated as follows: (1) an updated seismicity catalog was incorporated, which includes new earthquakes that occurred from 2013 to 2017; (2) in the central and eastern United States (CEUS), new ground motion models were updated that incorporate updated median estimates, modified assessments of the associated epistemic uncertainties and aleatory variabilities, and new soil amplification factors; (3) in the western United States (WUS), amplified shaking estimates of long-period ground motions at sites overlying deep sedimentary basins in the Los Angeles, San Francisco, Seattle, and Salt Lake City areas were incorporated; and (4) in the conterminous United States, seismic hazard is calculated for 22 periods (from 0.01 to 10 s) and 8 uniform V S30 maps (ranging from 1500 to 150 m/s). We also include a description of updated computer codes and modeling details. Results show increased ground shaking in many (but not all) locations across the CEUS (up to ~30%), as well as near the four urban areas overlying deep sedimentary basins in the WUS (up to ~50%). Due to population growth and these increased hazard estimates, more people live or work in areas of high or moderate seismic hazard than ever before, leading to higher risk of undesirable consequences from forecasted future ground shaking.

Earthquake Spectra↗
USGS2015

Earthquake hypocenters and focal mechanisms in central Oklahoma reveal a complex system of reactivated subsurface strike-slip faulting

The sharp increase in seismicity over a broad region of central Oklahoma has raised concern regarding the source of the activity and its potential hazard to local communities and energy industry infrastructure. Since early 2010, numerous organizations have deployed temporary portable seismic stations in central Oklahoma in order to record the evolving seismicity. In this study, we apply a multiple-event relocation method to produce a catalog of 3,639 central Oklahoma earthquakes from late 2009 through 2014. RMT source parameters were determined for 195 of the largest and best-recorded earthquakes. Combining RMT results with relocated seismicity enabled us to determine the length, depth and style-of-faulting occurring on reactivated subsurface fault systems. Results show that the majority of earthquakes occur on near vertical, optimally oriented (NE-SW and NW-SE), strike-slip faults in the shallow crystalline basement. These are necessary first order observations required to assess the potential hazards of individual faults in Oklahoma.

Oklahoma↗
USGS1999

Seismic slope-performance analysis: from hazard map to decision support system

In response to the growing recognition of engineers and decision-makers of the regional effects of earthquake-induced landslides, this paper presents a general approach to conducting seismic landslide zonation, based on the popular Newmark's sliding block analogy for modeling coherent landslides. Four existing models based on the sliding block analogy are compared. The comparison shows that the models forecast notably different levels of slope performance. Considering this discrepancy along with the limitations of static maps as a decision tool, a spatial decision support system (SDSS) for seismic landslide analysis is proposed, which will support investigations over multiple scales for any number of earthquake scenarios and input conditions. Most importantly, the SDSS will allow use of any seismic landslide analysis model and zonation approach. Developments associated with the SDSS will produce an object-oriented model for encapsulating spatial data, an object-oriented specification to allow construction of models using modular objects, and a direct-manipulation, dynamic user-interface that adapts to the particular seismic landslide model configuration.

Technical Council on Lifeline Earthquake Engineeri↗
USGS2003

Slip rate and earthquake recurrence along the central Septentrional fault, North American-Caribbean plate boundary, Dominican Republic

The Septentrional fault zone (SFZ) is the major North American-Caribbean, strike-slip, plate boundary fault at the longitude of eastern Hispaniola. The SFZ traverses the densely populated Cibao Valley of the Dominican Republic, forming a prominent scarp in alluvium. Our studies at four sites along the central SFZ are aimed at quantifying the late Quaternary behavior of this structure to better understand the seismic hazard it represents for the northeastern Caribbean. Our investigations of excavations at sites near Rio Cenovi show that the most recent ground-rupturing earthquake along this fault in the north central Dominican Republic occurred between A.D. 1040 and A.D. 1230, and involved a minimum of ???4 m of left-lateral slip and 2.3 m of normal dip slip at that site. Our studies of offset stream terraces at two locations, Rio Juan Lopez and Rio Licey, provide late Holocene slip rate estimates of 6-9 mm/yr and a maximum of 11-12 mm/yr, respectively, across the Septentrional fault. Combining these results gives a best estimate of 6-12 mm/yr for the slip rate across the SFZ. Three excavations, two near Tenares and one at the Rio Licey site, yielded evidence for the occurrence of earlier prehistoric earthquakes. Dates of strata associated with the penultimate event suggest that it occurred post-A.D. 30, giving a recurrence interval of 800-1200 years. These studies indicate that the SFZ has likely accumulated elastic strain sufficient to generate a major earthquake during the more than 800 years since it last slipped and should be considered likely to produce a destructive future earthquake.

Journal of Geophysical Research B: Solid Earth↗
USGS2021

Changes in liquefaction severity in the San Francisco Bay Area with sea-level rise

This paper studies the impacts of sea-level rise on liquefaction triggering and severity around the San Francisco Bay Area, California, for the M 7.0 “HayWired” earthquake scenario along the Hayward fault. This work emerged from stakeholder engagement for the US Geological Survey releases of the HayWired earthquake scenario and the Coastal Storm Modeling System projects, in which local planners and engineers asked where, why, and by how much liquefaction hazards may change due to sea-level rise in the future. We assess the impacts of sea-level rise on liquefaction by computing changes in liquefaction potential index (LPI) for over 400 cone penetration test (CPT) soundings around the San Francisco Bay for groundwater table models developed for current and increased sea levels of up to 5 m. For the M 7.0 HayWired earthquake scenario, we find that while the majority of sites are insensitive to sea-level changes of less than 1 m, some sites are highly sensitive to small changes in water levels. We then repeat these analyses for a uniform shaking scenario to isolate the effects of sea-level rise and we find similar patterns of change. For both earthquake scenarios, modest changes in overall LPI are expected for increases in sea level, but individual sites may see significant increases in liquefaction likelihood and severity.

Conference Paper↗
USGS2024

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington↗
USGS2018

Developing a global earthquake risk model

The understanding of earthquake risk is the first step towards the development and implementation of disaster risk reduction measures. However, in many countries, especially the countries of the developing world, earthquake risk models either do not exist or are publicly inaccessible. The Global Earthquake Model (GEM) Foundation and its partners have been supporting regional programmes and bilateral collaborations to develop a global earthquake risk model, due by the end of 2018. This paper describes how the main components (seismic hazard, exposure models and vulnerability functions) of this global effort are being collected, developed or improved. The calculations are being performed using the OpenQuake-engine, the open-source software for seismic hazard and risk calculations supported by GEM. This model will be able to provide estimates of critical risk metrics such as annualized average economic and human losses or aggregated losses for particular return periods, which are fundamental to the development of efficient and effective mitigation planning.

Conference Paper↗
USGS2022

Updates to and applications of the USGS National Crustal Model for seismic hazard studies

The U.S. Geological Survey (USGS) National Crustal Model (NCM) is being developed to assist in the modeling of seismic hazards across the conterminous United States. The NCM is composed of a grid of geophysical profiles, extending from the Earth’s surface into the upper mantle. It is constructed from a 3D geologic framework and geophysical rules defined by: (1) a petrologic and mineral physics database; (2) a 3D temperature model; and (3) a calibrated rock type- and age-dependent porosity model. Parameters needed to estimate site response for existing ground motion models (GMMs), including the time-averaged velocity in the upper 30 meters ( V S 30 ) and the depths to 1.0 and 2.5 km/s shear-wave velocity ( Z 1.0 and Z 2.5 ), can be extracted from the NCM. As GMMs develop, other metrics could also be extracted or derived from the NCM such as sediment thickness and travel times, fundamental frequency, a fully frequency-dependent site response function, or 3D geophysical volumes for wavefield simulations. Application of the NCM may also benefit other aspects of seismic hazard analysis including better accounting for path-dependent attenuation and geometric spreading, more accurate estimation of earthquake source properties such as hypocentral location and stress drop, and calculation of crustal strength profiles that inform estimates of the base of seismicity.

Conference Paper↗
USGS2024

Recent applications of the USGS National Crustal Model for Seismic Hazard Studies

The U.S. Geological Survey is developing the National Crustal Model (NCM) for seismic hazard studies to facilitate modeling site, path, and source components of seismic hazard across the conterminous United States. The NCM is composed of a 1km grid of geophysical profiles, extending from the Earth’s surface into the upper mantle. It is constructed from a threedimensional (3D) geologic framework and geophysical rules that use (1) a petrologic and mineral physics database; (2) a 3D temperature model; and (3) a calibrated rock type- and age-dependent porosity model. Parameters needed to estimate site response for existing ground motion models (GMMs), including the time-averaged velocity in the upper 30 meters (VS30), the depths to 1.0 and 2.5 km/s shear-wave velocity (Z1.0 and Z2.5), and sediment thickness, can be computed from the NCM. As GMMs continue to improve in the future, other metrics could also be extracted or derived from the NCM, such as fundamental period, site attenuation (ko), a fully frequency-dependent site response function, or 3D geophysical volumes for wavefield simulations. Application of the NCM may also benefit other aspects of seismic hazard analysis, including better accounting for path-dependent attenuation and geometric spreading, more accurate estimation of earthquake source properties such as hypocentral location and stress drop, and calculation of crustal strength profiles that inform estimates of the base of seismicity.

conterminous United States↗
USGS2017

Seismic displacement of gently-sloping coastal and marine sediment under multidirectional earthquake loading

Gentle sediment-laden slopes are typical of the onshore coastal zone and offshore continental shelf and slope. Coastal sediment are commonly young weakly consolidated materials that are well stratified, have low strength, and can mobilize shear displacements at low levels of stress. Seismically-driven plastic displacements of these sediment pose a hazard to coastal cities, buried onshore utilities, and offshore infrastructure like harbor protection and outfalls. One-dimensional rigid downslope-directed Newmark sliding block analyses have been used to predict earthquake deformations generally on steeper slopes that are modeled as frictional materials. This study probes the effect of multidirectional earthquake motions on inertial displacements of gently sloping ground of the coastal and offshore condition where soft-compliant soil is expected. Toward that objective, this investigation seeks to understand the effect on Newmark-type displacements of [1] multidirectional earthquake shaking and [2] soil compliance. In order to model multidirectional effects, the earthquake motions are rotated into the local slope strike- and dip-components. On gently sloping ground, including the strike component of motion always results in a larger and more accurate shear stress vector. Strike motions are found to contribute to downslope deformations on any declivity. Compliant response of the soil mass also influences the plastic displacements. The magnitude of seismic displacements can be estimated with a simplified model using only the estimated soil yield-acceleration ( k y ) and the peak ground velocity ( V max ) of the earthquake motions. Compliance effects can be effectively mapped using the concept of Plastic Displacement Response Spectra (PDRS).

Engineering Geology↗
USGS2009

Map of the Rinconada and Reliz Fault Zones, Salinas River Valley, California

The Rinconada Fault and its related faults constitute a major structural element of the Salinas River valley, which is known regionally, and referred to herein, as the 'Salinas Valley'. The Rinconada Fault extends 230 km from King City in the north to the Big Pine Fault in the south. At the south end of the map area near Santa Margarita, the Rinconada Fault separates granitic and metamorphic crystalline rocks of the Salinian Block to the northeast from the subduction-zone assemblage of the Franciscan Complex to the southwest. Northwestward, the Rinconada Fault lies entirely within the Salinian Block and generally divides this region into two physiographically and structurally distinct areas, the Santa Lucia Range to the west and the Salinas Valley to the east. The Reliz Fault, which continues as a right stepover from the Rinconada Fault, trends northwestward along the northeastern base of the Sierra de Salinas of the Santa Lucia Range and beyond for 60 km to the vicinity of Spreckels, where it is largely concealed. Aeromagnetic data suggest that the Reliz Fault continues northwestward another 25 km into Monterey Bay, where it aligns with a high-definition magnetic boundary. Geomorphic evidence of late Quaternary movement along the Rinconada and Reliz Fault Zones has been documented by Tinsley (1975), Dibblee (1976, 1979), Hart (1976, 1985), and Klaus (1999). Although definitive geologic evidence of Holocene surface rupture has not been found on these faults, they were regarded as an earthquake source for the California Geological Survey [formerly, California Division of Mines and Geology]/U.S. Geological Survey (CGS/USGS) Probabilistic Seismic Hazards Assessment because of their postulated slip rate of 1+-1 mm/yr and their calculated maximum magnitude of 7.3. Except for published reports by Durham (1965, 1974), Dibblee (1976), and Hart (1976), most information on these faults is unpublished or is contained in theses, field trip guides, and other types of reports. Therefore, the main purpose of this project is to compile and synthesize this body of knowledge into a comprehensive report for the geologic community. This report follows the format of Dibblee (1976) and includes discussions of the sections of the Rinconada Fault and of the Reliz Fault, as well as their Neogene history and key localities. Accompanying this report is a geologic map database of the faults, key localities, and earthquake epicenters, in ESRI shapefile format.

Scientific Investigations Map↗
USGS2017

Long-term afterslip of the M6.0, 2004 Parkfield, California, earthquake—Implications for forecasting amount and duration of afterslip on other major creeping faults

We present the longest record of surface afterslip on a continental strike‐slip fault for the 2004 M 6.0 Parkfield, California, earthquake, from which we can derive critical information about the duration and predictability of afterslip relevant to urban displacement hazard applications. Surface slip associated with this event occurred entirely postseismically along the interseismically creeping (0.6–1.5 cm/yr) main trace of the San Andreas fault. Using the first year of afterslip data, the program AFTER correctly predicted the cumulative surface afterslip (maximum ∼35 cm) eventually attained. By 1 yr postearthquake, observed afterslip had accumulated to only ∼74% of its modeled final value u f in units of length. The 6‐yr data suggested final slip would be reached everywhere by ∼6–12 yrs. Parkfield’s afterslip lasted much longer (∼6–12 yrs) than afterslip following a 2014 M 6.0 event in Napa, California, where no interseismic creep was known, and its afterslip neared completion (∼97% of u f ) by 1 yr. The uncertainty in u f for the Napa event fell to ≤2 cm in only three months, versus in 2 yrs for the Parkfield event, mostly because duration of the power‐law stage of afterslip at Parkfield is much longer, ∼1000 (493–1666) days versus ∼100 (35–421) days for Napa. Because the urban Hayward fault near San Francisco, California, like the Parkfield section, exhibits interseismic creep in a similar geological regime, significant afterslip might last for up to a decade following an anticipated M ≥6.7 earthquake, potentially delaying postearthquake recovery.

California↗
USGS2024

Precariously balanced rocks in northern New York and Vermont, U.S.A.: Ground-motion constraints and implications for fault sources

Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake‐source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground‐based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time‐independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site‐specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground‐motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M ∼6.5 earthquake‐source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake‐source faults in regions of cryptic faults.

Bulletin of the Seismological Society of America↗
USGS2018

Opportunities to enhance seismic demand parameters for future editions of the AS1170.4

Geoscience Australia has recently released its 2018 National Seismic Hazard Assessment (NSHA18). Results from the NSHA18 indicate significantly lower seismic hazard across almost all Australian localities at the 1/500 annual exceedance probability level relative to the factors adopted for the current Australian Standard AS1170.4–2007 (R2018). These new hazard estimates, coupled with larger probability factors (kp) for long return periods, have challenged notions of seismic hazard in Australia in terms of the recurrence of damaging ground motions. As a consequence, the new hazard estimates have raised questions over the appropriateness of the prescribed National Construction Code probability level as used in the AS1170.4 to determine appropriate seismic demands for the design of ordinary-use structures. Therefore, it is suggested that the ground-motion exceedance probability used in the current AS1170.4 be reviewed in light of the recent hazard assessment and the expected performance of modern buildings for rarer ground motions. Whilst adjusting the AS1170.4 exceedance probability level would be a major departure from previous earthquake loading standards, it would bring it into line with other international building codes in similar tectonic environments. Additionally, it would offer opportunities to further modernise how seismic demands are considered in Australian building design. In particular, the authors highlight the following additional opportunities: 1) the use of uniform hazard spectra to replace and simplify the spectral shape factors, which do not deliver uniform hazard across all natural periods; 2) updated site amplification factors to ensure continuity with modern ground-motion models, and; 3) the potential to define design ground motions in terms of uniform collapse risk rather than uniform hazard. Estimation of seismic hazard at any location is an uncertain science. However, as our knowledge improves, our estimates of the hazard will converge more closely to the actual – but unknowable – (time independent) hazard. It is therefore prudent to regularly update the estimates of the seismic demands in our building codes using the best available evidence-based methods and models.

Conference Paper↗
USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗