USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Conversation”

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,225 records · Page 68Linked to original sources

Effects of landscape characteristics on annual survival of Lesser Prairie-Chickens

Agriculture and development have caused landscape change throughout the southwestern Great Plains in the range of the lesser prairie-chicken ( Tympanuchus pallidicinctus ). Landscape alteration within the lesser prairie-chicken range may contribute to range contraction and population losses through decreases in survival rates. Our objectives were to determine if: (1) landscape configuration ( i.e ., the spatial arrangement of habitat) or composition ( i.e ., the amount of habitat), at the study-site scale, affected annual survival of females, (2) relationships exist between landscape context ( i.e ., landscape configuration and composition) and weekly survival to assess effects of landscape composition and configuration on lesser prairie-chicken populations, and (3) anthropogenic features influenced daily mortality risk. We captured 170 female lesser prairie-chickens and used very-high-frequency and GPS (Global Positioning System) transmitters to track their movement and survival for 2 y. We used known-fate survival models to test if landscape configuration or composition within three sites in Kansas were related to differences in female survival among sites. In addition we tested for relationships between weekly survival and landscape configuration or composition within home ranges. Finally, we used Andersen-Gill models to test the influence of distance to anthropogenic features on daily mortality risk. Differences in survival were evident between sites with differing landscape compositions as annual survival in Northwestern Kansas (Ŝ = 0.27) was half that of Clark County, Kansas (Ŝ = 0.56), which corresponded with 41.9% more grassland on the landscape in Clark County; landscape configuration did not measurably differ among sites. Survival was greater for prairie-chickens with home-ranges that had greater patch richness and in areas with 30% crop and 57% grassland. Female lesser prairie-chickens also experienced greater mortality risk closer to fences at patch edges. Further conversion of grassland landscapes occupied by lesser prairie-chickens should be avoided to reduce habitat loss and fragmentation thresholds that could affect survival. We suggest continued encouragement of Conservation Reserve Program enrollment in western areas of the lesser prairie-chicken range to maintain or increase the amount of grassland to increase annual survival.

Kansas↗

Fine scale movements and habitat use of black brant during the flightless Wing Molt in Arctic Alaska

Thousands of Black Brant ( Branta bernicla nigricans ) migrate annually to the Teshekpuk Lake Special Area (TLSA), Alaska, to undergo the flightless wing molt on tundra lakes and wetlands. GPS transmitters were attached to Brant over two summers (2007–2008) to examine patterns of movement and habitat use of molting Brant, including variation by habitat type, year and body mass. Molting Brant were located an average of 31 ± 1 m (SE) from shore and this distance did not vary across any of the explanatory variables. Brant moved an average of 123 ± 3 m hr -1 while flightless. Movement rates varied by year, averaging 22 ± 12 m hr -1 faster in 2008, and across habitat types, averaging 22 ± 13 m hr -1 faster in inland versus coastal and estuarine habitats. Two kernel home ranges were estimated: entire home range, which encompassed the complete 95% probability contour, and shoreline home range, which included only shoreline areas used by molting Brant. Entire home range (x bar = 15.1 ± 2.2 km 2 ) was negatively correlated with body mass, suggesting that heavier individuals have more body reserves to contribute to feather growth and thereby require less food and smaller home ranges. Conversely, shoreline home range (x bar = 4.3 ± 0.6 km 2 ) did not vary by body mass, but rather by habitat type, being larger in estuarine habitats. The complex shorelines and numerous deltaic islands of estuarine habitats offer more shoreline per area than either coastal or inland habitats. Brant appear to have limited ability to adjust their home range size or forage further from shore in response to variable food resources across years or habitats, instead altering their movement rate. Given this apparent lack of behavioral flexibility, Brant may be sensitive to development-related disturbances or habitat losses at molt sites in the TLSA.

Alaska↗

Fidelity and persistence of Ring-billed ( Larus delawarensis ) and Herring ( Larus argentatus ) gulls to wintering sites

While the breeding ecology of gulls (Laridae) has been well studied, their movements and spatial organization during the non-breeding season is poorly understood. The seasonal movements, winter-site fidelity, and site persistence of Ring-billed ( Larus delawarensis ) and Herring ( L. argentatus ) gulls to wintering areas were studied from 2008–2012. Satellite transmitters were deployed on Ring-billed Gulls ( n = 21) and Herring Gulls ( n = 14). Ten Ring-billed and six Herring gulls were tracked over multiple winters and > 300 wing-tagged Ring-billed Gulls were followed to determine winter-site fidelity and persistence. Home range overlap for individuals between years ranged between 0–1.0 (95% minimum convex polygon) and 0.31–0.79 (kernel utilization distributions). Ringbilled and Herring gulls remained at local wintering sites during the non-breeding season from 20–167 days and 74–161 days, respectively. The probability of a tagged Ring-billed Gull returning to the same site in subsequent winters was high; conversely, there was a low probability of a Ring-billed Gull returning to a different site. Ring-billed and Herring gulls exhibited high winter-site fidelity, but exhibited variable site persistence during the winter season, leading to a high probability of encountering the same individuals in subsequent winters.

Waterbirds↗

Preliminary evaluation of behavioral response of nesting waterbirds to small unmanned aircraft flight

Small unmanned aircraft systems present an emerging technology with the potential to survey colonial waterbird populations while reducing disturbance in comparison to traditional ground counts. Recent research with these systems has been performed on some colonially nesting avian species; however, none have focused on wading bird species. During 2015–2016, this study tested the behavioral response of a mixed-species rookery (Cattle Egret ( Bubulcus ibis ), Snowy Egret ( Egretta thula ), Glossy Ibis ( Plegadis falcinellus ) and a groundnesting colony of Common Terns ( Sterna hirundo )) in shrub habitat to small unmanned aircraft system flights at 12 m, 15 m, 30 m, and 50 m. Even at the lowest altitudes, the birds either showed no reaction or acclimated within 60 sec of the fly-over. Conversely, physically entering the colony to conduct ground surveys resulted in all Common Terns flushing from their nests beginning when the observer was 50 m away and required significantly more time in the colony overall: ~30–60 min vs. ~3–7 min with the small unmanned aircraft system. While this study focuses only on the behavioral response of nesting birds and not comparison of count estimates, these results provide preliminary evidence that small unmanned aircraft systems provide the potential to monitor colonial nesting bird populations while minimizing disturbance to the colony.

Waterbirds↗

Temporal and spatial shifts in habitat use by Black Brant immediately following flightless molt

Each year thousands of Pacific Black Brant ( Branta bernicla nigricans ) undergo flightless wing molt in the Teshekpuk Lake Special Area (TLSA), Alaska, in two distinct habitats: inland, freshwater lakes and coastal, brackish wetlands. Brant lose body mass during wing molt and likely must add reserves upon regaining flight to help fuel their 2,500 km migration to autumn staging areas. We characterized movements and habitat use by Brant during post-molt (the period immediately following the recovery of flight) by (1) marking individual Brant with GPS (global positioning system) transmitters, and (2) conducting a series of replicate aerial surveys. Individuals molting in inland habitats promptly abandoned their molt wetland during the post-molt and moved into coastal habitats. Consequently, inland habitats were nearly deserted by early August when Brant had regained flight, a decrease of >5,000 individuals from the flightless period of early July. Conversely, coastal molting Brant largely remained in coastal habitats during the post-molt and many coastal wetlands were occupied by large flocks (>1,000 birds). Our results indicate that inland, freshwater wetlands were less suitable post-molt habitats for Brant, while coastal wetlands were preferred as they transitioned from flightless molt. The immediacy with which Brant vacated inland habitats upon regaining flight suggests that food may be limiting during molt and they are not selecting inland molt sites strictly for food resources, but rather a balance of factors including predator avoidance and acquisition of protein for feather growth. Our data clearly demonstrate that patterns of habitat use by Brant in the TLSA change over the course of the molt season, an important consideration for management of future resource development activities in this area.

Alaska↗

Strong motion from surface waves in deep sedimentary basins

It is widely recognized that long-period surface waves generated by conversion of body waves at the boundaries of deep sedimentary basins make an important contribution to strong ground motion. The factors controlling the amplitude of such motion, however, are not widely understood. A study of pseudovelocity response spectra of strong-motion records from the Los Angeles Basin shows that late-arriving surface waves with group velocities of about 1 km/sec dominate the ground motion for periods of 3 sec and longer. The rate of amplitude decay for these waves is less than for the body waves and depends significantly on period, with smaller decay for longer periods. The amplitude can be modeled by the equation log y = f(M, RE) + c + bRB where y is the pseudovelocity response, f(M, RE) is an attenuation relation based on a general strong-motion data set, M is moment magnitude, RE is the distance from the source to the edge of the basin, RB is the distance from the edge of the basin to the recording site, and b and c are parameters fit to the data. The equation gives values larger by as much as a factor of 3 than given by the attenuation relationships based on general strong-motion data sets for the same source-site distance. It is clear that surface waves need to be taken into account in the design of long-period structures in deep sedimentary basins. The ground-motion levels specified by the earthquake provisions of current building codes, in California at least, accommodate the long-period ground motions from basin-edge-generated surface waves for periods of 5 sec and less and earthquakes with moment magnitudes of 7.5 or less located more than 20 km outside the basin. There may be problems at longer periods and for earthquakes located closer to the basin edge. The results of this study suggest that anelastic attenuation may need to be included in attempts to model long-period motion in deep sedimentary basins. To obtain better data on surface waves in the future, operators of strong-motion networks should take special care for the faithful recording of the long-period components of ground motion. It will also be necessary to insure that at least some selected recorders, once triggered, continue to operate for a time sufficient for the surface waves to traverse the basin. With velocities of about 1 km/sec, that time will be as long as 100 sec for a basin the size of the Los Angeles Basin.

Bulletin of the Seismological Society of America↗

Viscoelasticity, postseismic slip, fault interactions, and the recurrence of large earthquakes

The Brownian Passage Time (BPT) model for earthquake recurrence is modified to include transient deformation due to either viscoelasticity or deep post seismic slip. Both of these processes act to increase the rate of loading on the seismogenic fault for some time after a large event. To approximate these effects, a decaying exponential term is added to the BPT model's uniform loading term. The resulting interevent time distributions remain approximately lognormal, but the balance between the level of noise (e.g., unknown fault interactions) and the coefficient of variability of the interevent time distribution changes depending on the shape of the loading function. For a given level of noise in the loading process, transient deformation has the effect of increasing the coefficient of variability of earthquake interevent times. Conversely, the level of noise needed to achieve a given level of variability is reduced when transient deformation is included. Using less noise would then increase the effect of known fault interactions modeled as stress or strain steps because they would be larger with respect to the noise. If we only seek to estimate the shape of the interevent time distribution from observed earthquake occurrences, then the use of a transient deformation model will not dramatically change the results of a probability study because a similar shaped distribution can be achieved with either uniform or transient loading functions. However, if the goal is to estimate earthquake probabilities based on our increasing understanding of the seismogenic process, including earthquake interactions, then including transient deformation is important to obtain accurate results. For example, a loading curve based on the 1906 earthquake, paleoseismic observations of prior events, and observations of recent deformation in the San Francisco Bay region produces a 40% greater variability in earthquake recurrence than a uniform loading model with the same noise level.

Bulletin of the Seismological Society of America↗

Explosion source strong ground motions in the Mississippi embayment

Two strong-motion arrays were deployed for the October 2002 Embayment Seismic Excitation Experiment to study the spatial variation of strong ground motions in the deep, unconsolidated sediments of the Mississippi embayment because there are no comparable strong-motion data from natural earthquakes in the area. Each linear array consisted of eight three-component K2 accelerographs spaced 15 m apart situated 1.2 and 2.5 kin from 2268-kg and 1134-kg borehole explosion sources, respectively. The array data show distinct body-wave and surface-wave arrivals that propagate within the thick, unconsolidated sedimentary column, the high-velocity basement rocks, and small-scale structure near the surface. Time-domain coherence of body-wave and surface-wave arrivals is computed for acceleration, velocity, and displacement time windows. Coherence is high for relatively low-frequency verticalcomponent Rayleigh waves and high-frequency P waves propagating across the array. Prominent high-frequency PS conversions seen on radial components, a proxy for the direct S wave from earthquake sources, lose coherence quickly over the 105-m length of the array. Transverse component signals are least coherent for any ground motion and appear to be highly scattered. Horizontal phase velocity is computed by using the ratio of particle velocity to estimates of the strain based on a plane-wave-propagation model. The resulting time-dependent phase-velocity map is a useful way to infer the propagation mechanisms of individual seismic phases and time windows of three-component waveforms. Displacement gradient analysis is a complementary technique for processing general spatial-array data to obtain horizontal slowness information.

Arkansas, Missouri, Tennessee↗

Sedimentary basin effects in Seattle, Washington: Ground-motion observations and 3D simulations

Seismograms of local earthquakes recorded in Seattle exhibit surface waves in the Seattle basin and basin-edge focusing of S waves. Spectral ratios of Swaves and later arrivals at 1 Hz for stiff-soil sites in the Seattle basin show a dependence on the direction to the earthquake, with earthquakes to the south and southwest producing higher average amplification. Earthquakes to the southwest typically produce larger basin surface waves relative to S waves than earthquakes to the north and northwest, probably because of the velocity contrast across the Seattle fault along the southern margin of the Seattle basin. S to P conversions are observed for some events and are likely converted at the bottom of the Seattle basin. We model five earthquakes, including the M 6.8 Nisqually earthquake, using 3D finite-difference simulations accurate up to 1 Hz. The simulations reproduce the observed dependence of amplification on the direction to the earthquake. The simulations generally match the timing and character of basin surface waves observed for many events. The 3D simulation for the Nisqually earth-quake produces focusing of S waves along the southern margin of the Seattle basin near the area in west Seattle that experienced increased chimney damage from the earthquake, similar to the results of the higher-frequency 2D simulation reported by Stephenson et al. (2006). Waveforms from the 3D simulations show reasonable agreement with the data at low frequencies (0.2-0.4 Hz) for the Nisqually earthquake and an M 4.8 deep earthquake west of Seattle.

Bulletin of the Seismological Society of America↗

On baseline corrections and uncertainty in response spectrafor baseline variations commonly encountered in digital accelerograph records

Most digital accelerograph recordings are plagued by long-period drifts, best seen in the velocity and displacement time series obtained from integration of the acceleration time series. These drifts often result in velocity values that are nonzero near the end of the record. This is clearly unphysical and can lead to inaccurate estimates of peak ground displacement and long-period spectral response. The source of the long-period noise seems to be variations in the acceleration baseline in many cases. These variations could be due to true ground motion (tilting and rotation, as well as local permanent ground deformation), instrumental effects, or analog-to-digital conversion. Very often the trends in velocity are well approximated by a linear trend after the strong shaking subsides. The linearity of the trend in velocity implies that no variations in the baseline could have occurred after the onset of linearity in the velocity time series. This observation, combined with the lack of any trends in the pre-event motion, allows us to compute the time interval in which any baseline variations could occur. We then use several models of the variations in a Monte Carlo procedure to derive a suite of baseline-corrected accelerations for each noise model using records from the 1999 Chi-Chi earthquake and several earthquakes in Turkey. Comparisons of the mean values of the peak ground displacements, spectral displacements, and residual displacements computed from these corrected accelerations for the different noise models can be used as a guide to the accuracy of the baseline corrections. For many of the records considered here the mean values are similar for each noise model, giving confidence in the estimation of the mean values. The dispersion of the ground-motion measures increases with period and is noise-model dependent. The dispersion of inelastic spectra is greater than the elastic spectra at short periods but approaches that of the elastic spectra at longer periods. The elastic spectra from the most basic processing, in which only the pre-event mean is removed from the acceleration time series, do not diverge from the baseline-corrected spectra until periods of 10-20 sec or more for the records studied here, implying that for many engineering purposes elastic spectra can be used from records with no baseline correction or filtering.

Bulletin of the Seismological Society of America↗

Seismic imaging of the Waltham Canyon fault, California: comparison of ray‐theoretical and Fresnel volume prestack depth migration

Near‐vertical faults can be imaged using reflected refractions identified in controlled‐source seismic data. Often theses phases are observed on a few neighboring shot or receiver gathers, resulting in a low‐fold data set. Imaging can be carried out with Kirchhoff prestack depth migration in which migration noise is suppressed by constructive stacking of large amounts of multifold data. Fresnel volume migration can be used for low‐fold data without severe migration noise, as the smearing along isochrones is limited to the first Fresnel zone around the reflection point. We developed a modified Fresnel volume migration technique to enhance imaging of steep faults and to suppress noise and undesired coherent phases. The modifications include target‐oriented filters to separate reflected refractions from steep‐dipping faults and reflections with hyperbolic moveout. Undesired phases like multiple reflections, mode conversions, direct P and S waves, and surface waves are suppressed by these filters. As an alternative approach, we developed a new prestack line‐drawing migration method, which can be considered as a proxy to an infinite frequency approximation of the Fresnel volume migration. The line‐drawing migration is not considering waveform information but requires significantly shorter computational time. Target‐oriented filters were extended by dip filters in the line‐drawing migration method. The migration methods were tested with synthetic data and applied to real data from the Waltham Canyon fault, California. The two techniques are applied best in combination, to design filters and to generate complementary images of steep faults.

California↗

Optimizing earthquake early warning alert distance strategies using the July 2019 Mw6.4 and Mw7.1 Ridgecrest, California, earthquakes

The ShakeAlert earthquake early warning system aims to alert people who experience modified Mercalli intensity (MMI) IV+ shaking during an earthquake using source estimates (magnitude and location) to estimate median‐expected peak ground motions with distance, then using these ground motions to determine median‐expected MMI and thus the extent of MMI IV shaking. Because median ground motions are used, even if magnitude and location are correct, there will be people outside the alert region who experience MMI IV shaking but do not receive an alert (missed alerts). We use 91,000 “Did You Feel It?” survey responses to the July 2019 Mw 6.4 and Mw 7.1 Ridgecrest, California, earthquakes to determine which ground‐motion to intensity conversion equation (GMICE) best fits median MMI with distance. We then explore how incorporating uncertainty from the ground‐motion prediction equation and the GMICE in the alert distance calculation can produce more accurate MMI IV alert regions for a desired alerting strategy (e.g., aiming to alert 95% of people who experience MMI IV+ shaking), assuming accurate source characterization. Without incorporating ground‐motion uncertainties, we find MMI IV alert regions using median‐expected ground motions alert fewer than 20% of the population that experiences MMI IV+ shaking. In contrast, we find >94% of the people who experience MMI IV+ shaking can be included in the MMI IV alert region when two standard deviations of ground‐motion uncertainty are included in the alert distance computation. The optimal alerting strategy depends on the false alert tolerance of the community due to the trade‐off between minimizing missed and false alerts. This is especially the case for situations like the Mw 6.4 earthquake when alerting 95% of the 5 million people who experience MMI IV+ also results in alerting 14 million people who experience shaking below this level and do not need to take protective action.

California↗

Revisiting California’s past great earthquakes and long-term earthquake rate

In this study, we revisit the three largest historical earthquakes in California—the 1857 Fort Tejon, 1872 Owens Valley, and 1906 San Francisco earthquakes—to review their published moment magnitudes, and compare their estimated shaking distributions with predictions using modern ground‐motion models (GMMs) and ground‐motion intensity conversion equations. Currently accepted moment magnitude estimates for the three earthquakes are 7.9, 7.6, and 7.8, respectively. We first consider the extent to which the intensity distributions of all three earthquakes are consistent with a moment magnitude toward the upper end of the estimated range. We then apply a GMM‐based method to estimate the magnitudes of large historical earthquakes. The intensity distribution of the 1857 earthquake is too sparse to provide a strong constraint on magnitude. For the 1872 earthquake, consideration of all available constraints suggests that it was a high stress‐drop event, with a magnitude on the higher end of the range implied by scaling relationships, that is, higher than moment magnitude 7.6. For the 1906 earthquake, based on our analysis of regional intensities and the detailed intensity distribution in San Francisco, along with other available constraints, we estimate a preferred moment magnitude of 7.9, consistent with the published estimate based on geodetic and instrumental seismic data. These results suggest that, although there can be a tendency for historical earthquake magnitudes to be overestimated, the accepted catalog magnitudes of California’s largest historical earthquakes could be too low. Given the uncertainties of the magnitude estimates, the seismic moment release rate between 1850 and 2019 could have been either higher or lower than the average over millennial time scales. It is further not possible to reject the hypothesis that California seismicity is described by an untruncated Gutenberg–Richter distribution with a b "> b ‐value of 1.0 for moment magnitudes up to 8.0.

California↗

Structure and Qp-Qs relations in the Seattle and Tualatin basins from converted seismic phases

We use converted body‐wave phases from local earthquakes to constrain depth to basement and average attenuation relations for the Seattle basin in Washington and the Tualatin basin in Oregon. P ‐, P‐to‐S ‐( Ps ), S‐to‐P‐ ( Sp ), and S ‐wave arrivals are present in three‐component recordings of magnitude 2.5–4.0 earthquakes at seismic stations located in these basins. Based on their relative travel times, these phases are attributed to body‐wave conversions at the basement‐to‐basin contact or to high‐impedance interfaces within the basins. Depth to basement values are calculated using the differential travel times between direct and converted phases, as well as average P ‐ and S ‐wave velocity values. We also identify a high‐impedance layer in the Tualatin basin that likely represents a laterally extensive deposit of volcanic materials embedded between the basement contact and the Columbia River Basalt Group. In addition, the average Q P – Q S "> Q P - Q S attenuation relation is calculated for each station by taking the spectral ratio of converted phases to their parent body‐wave arrivals. For the Seattle basin, our analysis yields an average Q P "> Q P value of 73 and an average Q S "> Q S value of 60 for seismic waves with frequencies between 2 and 25 Hz. In the Tualatin basin, a much reduced Q P – Q S "> Q P – Q S relation suggests that average body‐wave attenuation is likely higher than in the Seattle basin. The converted phase techniques presented here provide a reliable way to develop estimates of basin depth and attenuation structure for undercharacterized regions using simple passive source seismic records.

Oregon, Washington↗

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S ‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

central Apennines↗

Capturing, preserving and digitizing legacy seismic data from the Montserrat Volcano Observatory analog seismic network, July 1995 – December 2004

An eruption of the Soufrière Hills Volcano (SHV) on the eastern Caribbean island of Montserrat began on 18 July 1995 and continued until February 2010. Within nine days of the eruption onset, an existing four‐station analog seismic network (ASN) was expanded to 10 sites. Telemetered data from this network were recorded, processed, and archived locally using a system developed by scientists from the U.S. Geological Survey (USGS) Volcano Disaster Assistance Program (VDAP). In October 1996, a digital seismic network (DSN) was deployed with the ability to capture larger amplitude signals across a broader frequency range. These two networks operated in parallel until December 2004, with separate telemetry and acquisition systems (analysis systems were merged in March 2001). Although the DSN provided better quality data for research, the ASN featured superior real‐time monitoring tools and captured valuable data including the only seismic data from the first 15 months of the eruption. These successes of the ASN have been rather overlooked. This article documents the evolution of the ASN, the VDAP system, the original data captured, and the recovery and conversion of more than 230,000 seismic events from legacy SUDS, Hypo71, and Seislog formats into Seisan database with waveform data in miniSEED format. No digital catalog existed for these events, but students at the University of South Florida have classified two-thirds of the 40,000 events that were captured between July 1995 and October 1996. Locations and magnitudes were recovered for ~10,000 of these events. Real-time seismic amplitude measurement, seismic spectral amplitude measurement, and tiltmeter data were also captured. The result is that the ASN seismic dataset is now more discoverable, accessible, and reusable, in accordance with FAIR data principles. These efforts could catalyze new research on the 1995–2010 SHV eruption. Furthermore, many observatories have data in these same legacy data formats and might benefit from procedures and codes documented here.

Soufrière Hills Volcano↗

Imaging the tectonic grain of the Northern Cordillera orogen using Transportable Array receiver functions

Azimuthal variations in receiver function conversions can image lithospheric structural contrasts and anisotropic fabrics that together compose tectonic grain. We apply this method to data from EarthScope Transportable Array in Alaska and additional stations across the northern Cordillera. The best‐resolved quantities are the strike and depth of dipping fabric contrasts or interfaces. We find a strong geographic gradient in such anomalies, with large amplitudes extending inboard from the present‐day subduction margin, the Aleutian arc, and an influence of flat‐slab subduction of the Yakutat microplate north of the Denali fault. An east–west band across interior Alaska shows low‐amplitude crustal anomalies. Anomaly amplitudes correlate with structural intensity (density of aligned geological elements), but are the highest in areas of strong Cenozoic deformation, raising the question of an influence of current stress state. Imaged subsurface strikes show alignment with surface structures. We see concentric strikes around arc volcanoes implying dipping magmatic structures and fabric into the middle crust. Regions with present‐day weaker deformation show lower anomaly amplitudes but structurally aligned strikes, suggesting pre‐Cenozoic fabrics may have been overprinted or otherwise modified. We observe general coherence of the signal across the brittle‐plastic transition. Imaged crustal fabrics are aligned with major faults and shear zones, whereas intrafault blocks show imaged strikes both parallel to and at high angles to major block‐bounding faults. High‐angle strikes are subparallel to neotectonic deformation, seismicity, fault lineaments, and prominent metallogenic belts, possibly due to overprinting and/or co‐evolution with fault‐parallel fabrics. We suggest that the underlying tectonic grain in the northern Cordillera is broadly distributed rather than strongly localized. Receiver functions thus reveal key information about the nature and continuity of tectonic fabrics at depth and can provide unique insights into the deformation history and distribution of regional strain in complex orogenic belts.

Alaska↗

Six decades of seismology at South Pole, Antarctica: Current limitations and future opportunities to facilitate new geophysical observations

Seismograms from the South Pole have been important for seismological observations for over six decades by providing (until 2007) the only continuous seismic records from the interior of the Antarctic continent. The South Pole, Antarctica station has undergone many updates over the years, including conversion to a digital recording station as part of the Global Seismographic Network (GSN) in 1991 and being relocated to multiple deep ( ⁠ &gt; 250 &#x2009;&#x2009; m "> > 250 m >250 m ⁠ ) boreholes 8 km away from the station in 2003 (and renamed to Quiet South Pole, Antarctica [QSPA]). Notably, QSPA is the second most used GSN station by the National Earthquake Information Center to pick phases used to rapidly detect and locate earthquakes globally, and has been used for a variety of glaciological and oceanography studies. In addition, it is the only seismic station on the Earth where low‐frequency ( ⁠ &lt; 5 &#x2009;&#x2009; mHz "> < 5 mHz <5 mHz ⁠ ), normal‐mode oscillations of the planet excited by large earthquakes can be recorded without influence from Earth’s rotation, and most of the direct effects of the solid Earth tide vanish. However, the current sensors are largely 1980s vintage, and, while able to make some lower‐frequency observations from earthquakes, the borehole sensors appear unable to resolve ambient ground motions at frequencies lower than 25 mHz due to instrument noise and contamination from magnetic field variations. Recently developed borehole sensors offer the potential to extend background noise observations to below 3 mHz, which would substantially improve the fidelity and scientific value of seismic observations at South Pole. Through collaboration with the IceCube Neutrino Observatory, the opportunity exists to emplace a modern very broadband seismometer near the base ( ⁠ &gt; 2 &#x2009;&#x2009; km "> > 2 km >2 km depth) of the Antarctic ice cap, which could lead to unprecedented seismic observations at long periods and facilitate a broad spectrum of Earth science studies.

Seismological Research Letters↗