USGS ScienceSearch

Geology topics

Eric A. Bergman

Publications and source records attributed to Eric A. Bergman.

11 recordsLinked to original sources

Warm places, warm years, and warm seasons increase parasitizing of moose by winter ticks

Observed links between parasites, such as ticks, and climate change have aroused concern for human health, wildlife population dynamics, and broader ecosystem effects. The one-host life history of the winter tick ( Dermacentor albipictus ) links each annual cohort to environmental conditions during three specific time periods when they are predictably vulnerable: spring detachment from hosts, summer larval stage, and fall questing for hosts. We used mixed-effects generalized linear models to investigate the drivers of tick loads carried by moose ( Alces alces ) relative to these time periods and across 750 moose, 10 years, and 16 study areas in the western United States. We tested for the effects of biotic factors (moose density, shared winter range, vegetation, migratory behavior) and weather conditions (temperature, snow, humidity) during each seasonal period when ticks are vulnerable and off-host. We found that warm climatic regions, warm seasonal periods across multiple partitions of the annual tick life cycle, and warm years relative to long-term averages each contributed to increased tick loads. We also found important effects of snow and other biotic factors such as host density and vegetation. Tick loads in the western United States were, on average, lower than those where tick-related die-offs in moose populations have occurred recently, but loads carried by some individuals may be sufficient to cause mortality. Lastly, we found interannual variation in tick loads to be most correlated with spring snowpack, suggesting this environmental component may have the highest potential to induce change in tick load dynamics in the immediate future of this region.

Ecosphere

A global catalog of calibrated earthquake locations

We produced a globally distributed catalog of earthquakes and nuclear explosions with calibrated hypocenters, referred to as the Global Catalog of Calibrated Earthquake Locations (GCCEL). This dataset currently contains 18,782 events in 289 clusters with >3.2 million arrival times observed at 19,258 stations. The term “calibrated” refers to the property that the hypocenters are minimally biased by unknown Earth structure. In addition, we calculate uncertainties using empirically determined variability of the arrival‐time data itself, specific to each calibrated cluster of hypocenters. Outliers in the arrival‐time dataset are removed based on measured variability of the data. In each cluster, we estimate the empirically determined uncertainty for each set of station‐phase arrival times. We use a version of the hypocentroidal decomposition multiple event relocation algorithm specifically adapted for calibrated relocations of clusters of seismic events. Most clusters are calibrated by fitting the subset of direct crustal first arrivals ( Pg and Sg ) with a locally appropriate travel‐time model to estimate the cluster hypocentroid. A few clusters are calibrated by aligning the pattern of relative locations in space and time with one or more events for which a ground‐truth hypocenter is available from an independent source with known uncertainty, such as a nuclear explosion. Epicentral uncertainties in GCCEL typically range from 1 to 5 km with a 90% confidence interval. Most events have depth constraint from one or more sources, usually with an uncertainty of ≤5 km. GCCEL is a significant resource for research at local, regional, and global scales because it provides minimally biased absolute hypocenters, meaningful associated error estimates, and curated arrival times as a reference dataset that can be used as prior constraints in the development of new regional, national, and global earthquake catalogs; validation of new location techniques; and the generation of advanced Earth models.

Seismological Research Letters

Structural control on megathrust rupture and slip behavior: Insights from the 2016 Mw 7.8 Pedernales Ecuador earthquake

The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 Mw 7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906 Mw 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a Mw 7.7 earthquake.

Journal of Geophysical Research: Earth Surface

The August 2018 Kaktovik earthquakes: Active tectonics in northeastern Alaska revealed With InSAR and seismology

The largest earthquakes recorded in northern Alaska (M w 6.4 and M w 6.0) occurred ~6 hours apart on August 12, 2018 in the northeastern Brooks Range. The earthquakes were captured by Sentinel-1 InSAR satellites and Earthscope Transportable Array seismic data, giving insight into the little-known active tectonic processes of Arctic Alaska, obscured until recently by sparse data availability. In this study, InSAR modelling, teleseismic back projections, calibrated hypocentral relocations and regional moment tensor solutions resolve two previously unknown, SSW-dipping right-lateral fault segments. These are the first active faults identified as conjugate to the NE-trending sinistral Canning Displacement Zone directly to the west, which is therefore a more complex zone of diffuse faulting than previously thought. The northeastern Brooks Range has been characterized as an area of low to moderate seismic hazard, but these earthquakes illustrate the potential for larger, possibly destructive events in a region earmarked for rapid resource development.

Alaska

The 12 November 2017 Mw 7.3 Ezgeleh–Sarpolzahab (Iran) earthquake and active tectonics of the Lurestan arc

The 12 November 2017 M w 7.3 Ezgeleh‐Sarpolzahab earthquake is the largest instrumentally recorded earthquake in the Zagros Simply Folded Belt by a factor of ∼10 in seismic moment. Exploiting local, regional, and teleseismic data and synthetic aperture radar interferometry imagery, we characterize the rupture, its aftershock sequence, background seismicity, and regional tectonics. The mainshock ruptured slowly (∼2 km/s), unilaterally southward, for ∼40 km along an oblique (dextral‐thrust) fault that dips ∼14°E beneath the northwestern Lurestan arc. Slip is confined to basement depths of ∼12–18 km, resolvably beneath the sedimentary cover which is ∼8 km thick in this area. The gentle dip angle and basement location allow for a broad slip area, explaining the large magnitude relative to earthquakes in the main Fars arc of the Zagros, where shallower, steeper faults are limited in rupture extent by weak sedimentary layers. Early aftershocks concentrate around the southern and western edges of the mainshock slip area and therefore cluster in the direction of rupture propagation, implying a contribution from dynamic triggering. A cluster of events ∼100 km to the south near Mandali (Iraq) reactivated the ∼50° dipping Zagros Foredeep Fault. The basement fault responsible for the Ezgeleh‐Sarpolzahab earthquake probably accounts for the ∼1 km elevation contrast between the Lurestan arc and the Kirkuk embayment but is distinct from sections of the Mountain Front Fault that define frontal escarpments elsewhere in the Zagros. It may be related to a seismic interface underlying the central and southern Lurestan arc, and a key concern is whether or not the more extensive regional structure is also seismogenic.

Journal of Geophysical Research: Solid Earth

2017 Valparaíso earthquake sequence and the megathrust patchwork of central Chile

In April 2017, a sequence of earthquakes offshore Valparaíso, Chile, raised concerns of a potential megathrust earthquake in the near future. The largest event in the 2017 sequence was a M 6.9 on 24 April, seemingly colocated with the last great-sized earthquake in the region—a M 8.0 in March 1985. The history of large earthquakes in this region shows significant variation in rupture size and extent, typically highlighted by a juxtaposition of large ruptures interspersed with smaller magnitude sequences. We show that the 2017 sequence ruptured an area between the two main slip patches during the 1985 earthquake, rerupturing a patch that had previously slipped during the October 1973 M 6.5 earthquake sequence. A significant gap in historic ruptures exists directly to the south of the 2017 sequence, with large enough moment deficit to host a great-sized earthquake in the near future, if it is locked.

Valparaiso

RMT focal plane sensitivity to seismic network geometry and faulting style

Modern tectonic studies often use regional moment tensors (RMTs) to interpret the seismotectonic framework of an earthquake or earthquake sequence; however, despite extensive use, little existing work addresses RMT parameter uncertainty. Here, we quantify how network geometry and faulting style affect RMT sensitivity. We examine how data-model fits change with fault plane geometry (strike and dip) for varying station configurations. We calculate the relative data fit for incrementally varying geometries about a best-fitting solution, applying our workflow to real and synthetic seismograms for both real and hypothetical station distributions and earthquakes. Initially, we conduct purely observational tests, computing RMTs from synthetic seismograms for hypothetical earthquakes and a series of well-behaved network geometries. We then incorporate real data and station distributions from the International Maule Aftershock Deployment (IMAD), which recorded aftershocks of the 2010 M W 8.8 Maule earthquake, and a set of regional stations capturing the ongoing earthquake sequence in Oklahoma and southern Kansas. We consider RMTs computed under three scenarios: (1) real seismic records selected for high data quality; (2) synthetic seismic records with noise computed for the observed source-station pairings and (3) synthetic seismic records with noise computed for all possible station-source pairings. To assess RMT sensitivity for each test, we observe the ‘fit falloff’, which portrays how relative fit changes when strike or dip varies incrementally; we then derive the ranges of acceptable strikes and dips by identifying the span of solutions with relative fits larger than 90 per cent of the best fit. For the azimuthally incomplete IMAD network, Scenario 3 best constrains fault geometry, with average ranges of 45° and 31° for strike and dip, respectively. In Oklahoma, Scenario 3 best constrains fault dip with an average range of 46°; however, strike is best constrained by Scenario 1, with a range of 26°. We draw two main conclusions from this study. (1) Station distribution impacts our ability to constrain RMTs using waveform time-series; however, in some tectonic settings, faulting style also plays a significant role and (2) increasing station density and data quantity (both the number of stations and the number of individual channels) does not necessarily improve RMT constraint. These results may be useful when organizing future seismic deployments (e.g. by concentrating stations in alignment with anticipated nodal planes), and in computing RMTs, either by guiding a more rigorous data selection process for input data or informing variable weighting among the selected data (e.g. by eliminating the transverse component when strike-slip mechanisms are expected).

Geophysical Journal International

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

The Mw 5.8 Mineral, Virginia, earthquake of August 2011 and aftershock sequence: constraints on earthquake source parameters and fault geometry

The M w 5.8 earthquake of 23 August 2011 (17:51:04 UTC) (moment, M0 5.7×10 17 N·m) occurred near Mineral, Virginia, within the central Virginia seismic zone and was felt by more people than any other earthquake in United States history. The U.S. Geological Survey (USGS) received 148,638 felt reports from 31 states and 4 Canadian provinces. The USGS PAGER system estimates as many as 120,000 people were exposed to shaking intensity levels of IV and greater, with approximately 10,000 exposed to shaking as high as intensity VIII. Both regional and teleseismic moment tensor solutions characterize the earthquake as a northeast‐striking reverse fault that nucleated at a depth of approximately 7±2 km. The distribution of reported macroseismic intensities is roughly ten times the area of a similarly sized earthquake in the western United States (Horton and Williams, 2012). Near‐source and far‐field damage reports, which extend as far away as Washington, D.C., (135 km away) and Baltimore, Maryland, (200 km away) are consistent with an earthquake of this size and depth in the eastern United States (EUS). Within the first few days following the earthquake, several government and academic institutions installed 36 portable seismograph stations in the epicentral region, making this among the best‐recorded aftershock sequences in the EUS. Based on modeling of these data, we provide a detailed description of the source parameters of the mainshock and analysis of the subsequent aftershock sequence for defining the fault geometry, area of rupture, and observations of the aftershock sequence magnitude–frequency and temporal distribution. The observed slope of the magnitude–frequency curve or b‐value for the aftershock sequence is consistent with previous EUS studies (b=0.75), suggesting that most of the accumulated strain was released by the mainshock. The aftershocks define a rupture that extends between approximately 2–8 km in depth and 8–10 km along the strike of the fault plane. Best‐fit modeling of the geometry of the aftershock sequence defines a rupture plane that strikes N36°E and dips to the east‐southeast at 49.5°. Moment tensor solutions of the mainshock and larger aftershocks are consistent with the distribution of aftershock locations, both indicating reverse slip along a northeast–southwest striking southeast‐dipping fault plane.

Virginia

Toward a comprehensive catalog of global historical seismicity

The United States Geological Survey (USGS) and the Colorado School of Mines (CSM) have initiated a project to locate more accurately all earthquakes recorded by instruments during the period 1900 to 1963. Seismicity for this period (hereafter referred to as historical seismicity, following Lee et al . [1988]) is still poorly understood, even for basic parameters such as earthquake locations (see Figure 1). In some cases this is the result of inherent limitations in the distribution, response characteristics, and timing of the instruments. However, locations for most of the pre-1964 earthquakes are poorly determined simply because modern data analysis techniques have yet to be applied to the available arrival-time observations, which are mainly preserved as printed bulletins and not in a computer-ready digital format. The arduous task of hand-entering these data has prevented the systematic analysis and relocation of historical seismicity.

Eos, Transactions, American Geophysical Union

The international seismological observing period in Africa

The International Seismological Observing Period (ISOP) is a specific time interval designated for enhanced international cooperation in the collection and dissemination of observatory measurements from the global seismographic network. The primary purpose of the ISOP is to strengthen the international infrastructure that supports current seismological practice and increase the cooperation among nations that operate seismological observatories. Measurements, reported by the existing global network and compiled by agencies such as the International Seismological Centre (ISC), are providing new information about earthquakes and the structure of the Earth of fundamental importance to the Earth sciences. However, these data represent but a small fraction of the information contained in the seismograms. One of the goals of the ISOP is to collect improved sets of data. In particular, the measurement and reporting of later-arriving phases, during a fixed ISOP period, from earthquakes selected for detailed observation by the cooperating stations will be encouraged. The use of advanced, digital instrumentation provides an unprecedented opportunity for enhancing the methods of seismogram interpretation and seismic parameter extraction, by the implementation of digital processing methods at seismic observatories worldwide. It must be ensured that this new information will be available to the entire seismological community. It is believed that this purpose is best served with an ISOP that promotes increased on-site processing at digital stations in Africa and elsewhere. Improvements in seismology require truly international cooperation and the educational aspects of seismological practice form one of the goals of the ISOP. Thus, workshops will be needed in Africa to train analysts in ISOP procedures and to introduce them to modern techniques and applications of the data. Participants will, thus, benefit from theoretical results and practical experience that are of direct relevance to their own work.

Tectonophysics