SEARCH · USGS Science

Results for “Earth in Space”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46Linked to original sources

USGS1998

Strain accumulation rates in the San Francisco Bay area, 1972–1989

Maps of the strain accumulation rate in the San Francisco Bay area have been constructed from trilateration observations extending from about 1972 until the Loma Prieta earthquake in late 1989. The observations were corrected to remove offsets imposed by shallow fault creep and by four M ∼6 earthquakes that occurred in the Bay area during that time interval. The Bay area was divided into 32 contiguous polygons, and the uniform (in both space and time) strain rates that best explain the changes in the corrected (earthquake and shallow fault creep offsets removed) distances within each polygon were calculated. In a coordinate system with the 1 axis directed N58°E and the 2 axis N32°W (perpendicular and parallel to the local tangent to the small circle drawn about the Pacific-Sierra Nevada pole of rotation) the averages of these 32 strain rates (each weighted by the area of the polygon) are ε11=9.2±7.4, ε12=−160.7±4.6, and ε22=8.2±6.2 nanostrain/yr, where extension is reckoned positive and quoted uncertainties are standard deviations. As expected from the Pacific-Sierra Nevada relative plate motion, the overall strain rate is predominantly right-lateral shear across a vertical plane striking N32°W. The net increase in the 12,225 km2 area of the trilateration network is only 212±110 m2/yr, which arises from almost equal extensions in the N32°W and N58°E directions. Within the network the strain rates vary from polygon to polygon. Significant areal dilatation rates are observed in almost ⅓ of the individual polygons and the N32°W extension rates tend to be negative to the west of the Hayward-Rodgers Creek fault trend and positive east of it. The N58E extension rate is positive in 22 out of 32 polygons, a proportion that is significantly larger than would be expected by chance if the N58°E extension rate were zero or negative. The pre-1989 strain accumulation across the eventual site of the Loma Prieta rupture involves fault normal contraction as well as right-lateral shear, consistent with the rupture mechanism.

California↗

USGS2000

Preliminary report on the 16 October 1999 M 7.1 Hector mine, California, earthquake

The M w 7.1 Hector Mine, California, earthquake occurred at 9:46 GMT on 16 October 1999. The event caused minimal damage because it was located in a remote, sparsely populated part of the Mojave Desert, approximately 47 miles east-southeast of Barstow, with epicentral coordinates 34.59°N 116.27°W and a hypocentral depth of 5 ± 3 km. Twelve foreshocks, M 1.9-3.8, preceded the mainshock during the previous twelve hours. All of these events were located close to the hypocenter of the mainshock. The Hector Mine earthquake occurred within the Eastern California Shear Zone (ECSZ). By virtue of its remote location, the societal impact of the Hector Mine earthquake was, fortunately, minimal in spite of the event's appreciable size. The ECSZ is characterized by high seismicity, a high tectonic strain rate, and a broad, distributed zone of north-northwest-trending faults (ECSZ; Figure 1 ; Dokka and Travis, 1990 ; Sauber et al., 1986 ; Sauber et al., 1994 ; Sieh et al., 1993 ). Data regarding the slip rates of faults within the ECSZ suggest that on the order of 15% of the Pacific-North American plate motion occurs along this zone ( Sauber et al., 1986 ; Wesnousky, 1986 ). Most of the faults in the ECSZ have low slip rates and long repeat times for major earthquakes, on the order of several thousands to tens of thousands of years. The occurrence of the Hector Mine earthquake within seven years and only about 30 km east of the 1992 M w 7.3 Landers earthquake suggests that the closely spaced surface faults in the ECSZ are mechanically related. The Hector Mine event involved rupture on two previously mapped fault zones—the Bullion Fault and an unnamed, more northerly-trending fault that is informally referred to in this paper as the Lavic Lake Fault (Dibblee, 1966 , 1967a , b ). Traces of the Bullion Fault exhibit evidence of Holocene displacement and were zoned as active in 1988 under California's Mquist-Priolo Earthquake Fault Zoning Act ( Hart and Bryant, 1997 ). The pattern of rupture along more than one named fault was also observed from the 1992 Landers earthquake ( Hauksson et al., 1993 ; Sieh et al., 1994). Much of the fault zone that produced the Hector Mine earthquake had been buried by relatively young stream deposits, and the fault scarps in bedrock have a subdued morphology. It appears that these faults have not experienced significant offset for perhaps 10,000 years or more ( Hart, 1987 ). Planned future investigations will refine the age of the last event on these faults. The portion of the Lavic Lake Fault that ruptured between the northern end of the Bullion Mountains and Lavic Lake had not previously been mapped. However, our field investigations have identified ancient, subdued fault scarps along portions of the 1999 rupture zone in this area. It thus appears that the entire segment of the Lavic Lake Fault that was involved in the 1999 event had ruptured in the past. As is typical for most faults within the Eastern California Shear Zone, the rate of movement along the Lavic Lake Fault may be quite slow (<1 mm/yr) and should produce earthquakes only infrequently. This event is a reminder that faults that have ruptured in late Quaternary time, but that lack evidence of Holocene displacement, can still produce earthquakes in this low-slip-rate tectonic setting. Additionally, the Hector Mine earthquake is noteworthy for a couple of other reasons. First, it clearly produced triggered seismicity over much of southern California, from the rupture zone toward the south-southwest in particular. Second, as we will discuss, the event may provide new data and insight into recently developed paradigms concerning earthquake interactions and the role of static stress changes. Questions such as these will, of course, be the subject of extensive detailed analyses in years to come. Fortunately, the Hector Mine sequence will provide one of the best data sets obtained to date for a significant earthquake in the United States. Because it occurred when major upgrades to both the regional seismic network (TriNet) and the regional geodetic network (SCIGN) were well underway, the Earth science community will have abundant high-quality data with which to explore the important and interesting questions that have been raised. In this paper, we present and discuss the basic data and preliminary results from the Hector Mine earthquake.

California↗

USGS1996

Cyclic stressing and seismicity at strongly coupled subduction zones

We use the finite element method to analyze stress variations in and near a strongly coupled subduction zone during an earthquake cycle. Deformation is assumed to be uniform along strike (plane strain on a cross section normal to the trench axis), and periodic earthquake slip is imposed consistent with the long-term rate of plate convergence and degree of coupling. Simulations of stress and displacement rate fields represent periodic fluctuations in time superimposed on an average field. The oceanic plate, descending slab, and continental lithosphere are assumed here to respond elastically to these fluctuations, and the remaining mantle under and between plates is assumed to respond as Maxwell viscoelastic. In the first part of the analysis we find that computed stress fluctuations in space and time are generally consistent with observed earthquake mechanism variations with time since a great thrust event. In particular, trench-normal extensional earthquakes tend to occur early in the earthquake cycle toward the outer rise but occur more abundantly late in the cycle in the subducting slab downdip of the main thrust zone. Compressional earthquakes, when they occur at all, have the opposite pattern. Our results suggest also that the actual timing of extensional outer rise events is controlled by the rheology of the shallow aseismic portion of the thrust interface. The second part of the analysis shows the effects of mantle relaxation on the rate of ground surface deformation during the earthquake cycle. Models without relaxation predict a strong overall compressional strain rate in the continental plate above the main thrust zone, with the strain rate constant between mainshocks. However with significant relaxation present, a localized region of unusually low compressional, or even slightly extensional, strain rate develops along the surface of the continental plate above and somewhat inland from the downdip edge of the locked main thrust zone. The low strain rate starts in the middle or late part of the cycle, depending on position. This result suggests that the negligible or small contraction measured on the Shumagin Islands, Alaska, during 1980 to 1991, may not invalidate an interpretation of that region as being a moderately coupled subduction zone. In contrast, mantle relaxation causes only modest temporal nonuniformity of uplift rates in the overriding plate and of extensional stress rates in the subducting plate, even when the Maxwell time is an order of magnitude less than the recurrence interval.

Journal of Geophysical Research B: Solid Earth↗

USGS1998

Waveform inversion of very long period impulsive signals associated with magmatic injection beneath Kilauea Volcano, Hawaii

We use data from broadband seismometers deployed around the summit of Kilauea Volcano to quantify the mechanism associated with a transient in the flow of magma feeding the east rift eruption of the volcano. The transient is marked by rapid inflation of the Kilauea summit peaking at 22 μrad 4.5 hours after the event onset, followed by slow deflation over a period of 3 days. Superimposed on the summit inflation is a series of sawtooth displacement pulses, each characterized by a sudden drop in amplitude lasting 5–10 s followed by an exponential recovery lasting 1–3 min. The sawtooth waveforms display almost identical shapes, suggesting a process involving the repeated activation of a fixed source. The particle motion associated with each sawtooth is almost linear, and its major swing shows compressional motion at all stations. Analyses of semblance and particle motion are consistent with a point source located 1 km beneath the northeast edge of the Halemaumau pit crater. To estimate the source mechanism, we apply a moment tensor inversion to the waveform data, assuming a point source embedded in a homogeneous half-space with compressional and shear wave velocities representative of the average medium properties at shallow depth under Kilauea. Synthetic waveforms are constructed by a superposition of impulse responses for six moment tensor components and three single force components. The origin times of individual impulses are distributed along the time axis at appropriately small, equal intervals, and their amplitudes are determined by least squares. In this inversion, the source time functions of the six tensor and three force components are determined simultaneously. We confirm the accuracy of the inversion method through a series of numerical tests. The results from the inversion show that the waveform data are well explained by a pulsating transport mechanism operating on a subhorizontal crack linking the summit reservoir to the east rift of Kilauea. The crack acts like a buffer in which a batch of fluid (magma and/or gas) accumulates over a period of 1–3 min before being rapidly injected into a larger reservoir (possibly the east rift) over a timescale of 5–10 s. The seismic moment and volume change associated with a typical batch of fluid are approximately 10 14 N m and 3000 m 3 , respectively. Our results also point to the existence of a single force component with amplitude of 10 9 N, which may be explained as the drag force generated by the flow of viscous magma through a narrow constriction in the flow path. The total volume of magma associated with the 4.5-hour-long activation of the pulsating source is roughly 500,000 m 3 in good agreement with the integrated volume flow rate of magma estimated near the eruptive site.

Journal of Geophysical Research B: Solid Earth↗

USGS2023

Summary of the history and research of the U.S. Geological Survey gas hydrate properties laboratory in Menlo Park, California, active from 1993 to 2022

The U.S. Geological Survey (USGS) Clathrate Hydrate Properties Project, active from 1993 to 2022 in Menlo Park, California, stemmed from an earlier project on the properties of planetary ices supported by the National Aeronautics and Space Administration’s (NASA’s) Planetary Geology and Geophysics Program. We took a material science approach in both projects, emphasizing chemical purity of samples, having controlled grain size and grain texture, and having verified crystal structures and phase relations. A foundational contribution from our USGS Gas Hydrate Properties Laboratory (GHPL) was in demonstrating the ability to reproducibly create such pure clathrate hydrate samples for study. Clathrate sample synthesis was achieved by heating sieved and weighed pure granular water ice in the presence of cold clathrate-forming gas or liquid. During heating, the ice melts at the grain scale and reacts with the gas to form clathrate. The resulting material has the desired uniformity and purity, with known intergranular porosity; our subsequent measurements showed that these clathrates exhibited the established clathrate structures and phase relations. This novel synthesis method was successful in creating clathrates of pure methane, ethane, propane, carbon dioxide, and multi-component gases. By mixing sand or silt with granular ice, we were also able to make clathrate-sediment aggregates with controlled grain textures. This simple method, adopted by many others in the community, permitted us to measure the physical and chemical properties of well-characterized and well-crystallized clathrates and clathrate/sediment aggregates. At about the same time, we adapted conventional scanning electron microscopy to cryogenic conditions for analysis of grain-scale characteristics of clathrates made in the GHPL as well as those collected from nature by drill core. The uniformity and reproducibility of our samples also allowed us to investigate how clathrates respond to environmental changes in chemistry, temperature, and pressure: we measured chemical exchange rates with dissolved gas species—such as noble gases and chlorofluorocarbons—as well as rates of clathrate dissolution and decomposition. These advances include the first accurate mapping of the conditions that promote the remarkable process of “anomalous preservation” at room pressure, a metastability that offers potential application for low-cost and safe transportation of natural gas from gas fields far from pipelines. Another advancement stemming from the GHPL was the compaction of as-synthesized porous clathrates to nearly full density by applying external pressure using three different techniques. Compaction allows for high-accuracy measurements of many fundamental physical and chemical properties of these materials, such as elastic wavespeeds and moduli, complete thermal properties, decomposition rates, thermal expansion, and clathrate equations of state. These properties and others, in turn, have helped USGS scientists to interpret geophysical well logs and active geophysical surveys, as well as model the rates of gas production from hydrate deposits in nature. Studying this class of icy minerals that occur in abundance on Earth and in the outer solar system has been a fascinating laboratory journey. Here, we summarize the history and major findings of the USGS GHPL in Menlo Park, including both in-house research as well as findings from the synergistic collaborations with other agencies and institutes that were key to the success of our laboratory. The Menlo Park GHPL was more formally incorporated within the USGS Gas Hydrates Project, a collaboration among multiple USGS Science Centers, in the early 2000s under the leadership of Deborah Hutchinson, and now under the leadership of Carolyn Ruppel and Timothy Collett.

Open-File Report↗

USGS2023

Correction to “The March 1940 superstorm: Geoelectromagnetic hazards and impacts on American communication and power systems”

In Love et al. (2023), 5 estimated geopotentials on long lines for the March 1989 storm are incorrect. This was the result of incorrect calculation of line integrals. This error affects the two paragraphs of Section 11 and Figure 10 of that paper. The most significant impact of this error is that we overestimated the average ratio of the geopotentials measured the long lines during the March 1940 storm to the estimated geopotentials for the March 1989 storm. Instead of being a factor of 9 higher, the average ratio is only 9% higher. This means the inductive geoeffectiveness of the March 1940 storm was comparable to, or just slightly greater than, the geoeffectiveness of the March 1989 storm. Corrected paragraphs and a corrected figure are provided here.

Space Weather↗

USGS2021

The making of the NEAM Tsunami Hazard Model 2018 (NEAMTHM18)

The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-eastern Atlantic, the Mediterranean, and connected seas (NEAM). NEAMTHM18 was designed as a three-phase project. The first two phases were dedicated to the model development and hazard calculations, following a formalized decision-making process based on a multiple-expert protocol. The third phase was dedicated to documentation and dissemination. The hazard assessment workflow was structured in Steps and Levels. There are four Steps: Step-1) probabilistic earthquake model; Step-2) tsunami generation and modeling in deep water; Step-3) shoaling and inundation; Step-4) hazard aggregation and uncertainty quantification. Each Step includes a different number of Levels. Level-0 always describes the input data; the other Levels describe the intermediate results needed to proceed from one Step to another. Alternative datasets and models were considered in the implementation. The epistemic hazard uncertainty was quantified through an ensemble modeling technique accounting for alternative models’ weights and yielding a distribution of hazard curves represented by the mean and various percentiles. Hazard curves were calculated at 2,343 Points of Interest (POI) distributed at an average spacing of ∼20 km. Precalculated probability maps for five maximum inundation heights (MIH) and hazard intensity maps for five average return periods (ARP) were produced from hazard curves. In the entire NEAM Region, MIHs of several meters are rare but not impossible. Considering a 2% probability of exceedance in 50 years (ARP≈2,475 years), the POIs with MIH >5 m are fewer than 1% and are all in the Mediterranean on Libya, Egypt, Cyprus, and Greece coasts. In the North-East Atlantic, POIs with MIH >3 m are on the coasts of Mauritania and Gulf of Cadiz. Overall, 30% of the POIs have MIH >1 m. NEAMTHM18 results and documentation are available through the TSUMAPS-NEAM project website (http://www.tsumaps-neam.eu/), featuring an interactive web mapper. Although the NEAMTHM18 cannot substitute in-depth analyses at local scales, it represents the first action to start local and more detailed hazard and risk assessments and contributes to designing evacuation maps for tsunami early warning.

Frontiers in Earth Science↗

USGS2016

Large along-strike variations in the onset of Subandean exhumation: Implications for Central Andean orogenic growth

Plate tectonics drives mountain building in general, but the space-time pattern and style of deformation is influenced by how climate, geodynamics, and basement structure modify the orogenic wedge. Growth of the Subandean thrust belt, which lies at the boundary between the arid, high-elevation Central Andean Plateau and its humid, low-elevation eastern foreland, figures prominently into debates of orogenic wedge evolution. We integrate new apatite and zircon (U-Th)/He thermochronometer data with previously published apatite fission-track data from samples collected along four Subandean structural cross-sections in Bolivia between 15° and 20°S. We interpret cooling ages vs. structural depth to indicate the onset of Subandean exhumation and signify the forward propagation of deformation. We find that Subandean growth is diachronous south (11 ± 3 Ma) vs. north (6 ± 2 Ma) of the Bolivian orocline and that Subandean exhumation magnitudes vary by more than a factor of two. Similar north-south contrasts are present in foreland deposition, hinterland erosion, and paleoclimate; these observations both corroborate diachronous orogenic growth and illuminate potential propagation mechanisms. Of particular interest is an abrupt shift to cooler, more arid conditions in the Altiplano hinterland that is diachronous in southern Bolivia (16-13 Ma) vs. northern Bolivia (10-7 Ma) and precedes the timing of Subandean propagation in each region. Others have interpreted the paleoclimate shift to reflect either rapid surface uplift due to lithosphere removal or an abrupt change in climate dynamics once orographic threshold elevations were exceeded. These mechanisms are not mutually exclusive and both would drive forward propagation of the orogenic wedge by augmenting the hinterland backstop, either through surface uplift or spatially variable erosion. In summary, we suggest that diachronous Subandean exhumation was driven by piecemeal hinterland uplift, orography, and the outward propagation of deformation.

Andes↗

USGS2022

Bathymetry retrieval from CubeSat image sequences with short time lags

The rapid expansion of CubeSat constellations could revolutionize the way inland and nearshore coastal waters are monitored from space. This potential stems from the ability of CubeSats to provide daily imagery with global coverage at meter-scale spatial resolution. In this study, we explore the unique opportunity to improve the retrieval of bathymetry offered by CubeSats, specifically those of the PlanetScope constellation. The orbital design of the PlanetScope constellation enables the acquisition of image sequences with short time lags (from seconds to hours). This characteristic allows multiple images to be captured during a short period of steady bathymetric conditions, especially in dynamic environments like rivers. We hypothesize that taking the ensemble mean of a CubeSat image sequence can enhance bathymetry retrieval compared to standard single-image analysis. Along with the existing optimal band ratio analysis (OBRA) algorithm, we also use a new neural network-based depth retrieval (NNDR) technique to infer bathymetry from both individual and time-averaged images. The two methodologies are evaluated using field data from five different river reaches with depths up to 15 m and both top-of-atmosphere (TOA) radiance and bottom-of-atmosphere (BOA) surface reflectance PlanetScope data products. Despite low spectral resolution and concerns about the radiometric quality of CubeSat imagery, accuracy assessment based on in-situ comparisons indicates the potential (0.52 < R 2 < 0.7 for the NNDR method) of PlanetScope imagery to retrieve depths up to ∼ 10 m in clear water conditions. The proposed image averaging consistently improves bathymetry retrieval over single image analysis. The NNDR technique was found to outperform OBRA, illustrating the importance of leveraging all spectral bands through machine learning approaches. TOA data provided more robust bathymetry results than BOA data for the OBRA technique, but the NNDR technique was minimally impacted by the type of data product.

International Journal of Applied Earth Observation↗

USGS2001

Viscoelastic shear zone model of a strike-slip earthquake cycle

I examine the behavior of a two-dimensional (2-D) strike-slip fault system embedded in a 1-D elastic layer (schizosphere) overlying a uniform viscoelastic half-space (plastosphere) and within the boundaries of a finite width shear zone. The viscoelastic coupling model of Savage and Prescott [1978] considers the viscoelastic response of this system, in the absence of the shear zone boundaries, to an earthquake occurring within the upper elastic layer, steady slip beneath a prescribed depth, and the superposition of the responses of multiple earthquakes with characteristic slip occurring at regular intervals. So formulated, the viscoelastic coupling model predicts that sufficiently long after initiation of the system, (1) average fault-parallel velocity at any point is the average slip rate of that side of the fault and (2) far-field velocities equal the same constant rate. Because of the sensitivity to the mechanical properties of the schizosphere-plastosphere system (i.e., elastic layer thickness, plastosphere viscosity), this model has been used to infer such properties from measurements of interseismic velocity. Such inferences exploit the predicted behavior at a known time within the earthquake cycle. By modifying the viscoelastic coupling model to satisfy the additional constraint that the absolute velocity at prescribed shear zone boundaries is constant, I find that even though the time-averaged behavior remains the same, the spatiotemporal pattern of surface deformation (particularly its temporal variation within an earthquake cycle) is markedly different from that predicted by the conventional viscoelastic coupling model. These differences are magnified as plastosphere viscosity is reduced or as the recurrence interval of periodic earthquakes is lengthened. Application to the interseismic velocity field along the Mojave section of the San Andreas fault suggests that the region behaves mechanically like a ???600-km-wide shear zone accommodating 50 mm/yr fault-parallel motion distributed between the San Andreas fault system and Eastern California Shear Zone. Copyright 2001 by the American Geophysical Union.

California↗

USGS2012

Regional deformation of the Sierra Nevada, California, on conjugate microfault sets

Strike slip microfaults are pervasive throughout the granitic rocks of the eastern Sierra Nevada. Offsets typically range from less than a millimeter to several tens of centimeters but exceed 100 m in some places. The spacing between microfaults varies from a few tens of centimeters to a few tens of meters throughout much of the high Sierra Nevada. Many of these microfaults are loci of slickensided, compact fault gouge, and they are commonly mineralized by quartz veinlets with minor amounts of epidote, chlorite, and rare sulphide minerals. The microfaults are oriented in two nearly vertical conjugate sets; a north to northeast striking set showing right lateral offset and an east to northeast striking set showing left lateral offset. Microfaults with left lateral offset are more common than microfaults with right lateral offset. Most lineaments visible on aerial photographs are microfaults. The age of this microfaulting is not precisely known. It developed after consolidation of the youngest granitic plutons in the Sierra (79 m.y. B.P.) and is known to cut a late Miocene volcanic dike in one area. Slickensides along the microfaults are subhorizontal but show a slight (about 3°) westward plunging average inclination, suggesting that much of the deformation occurred prior to the westward tilting of the Sierran block in late Tertiary time. The direction of maximum horizontal extensional strain (determined as the bisector of average microfault trends) changes systematically from north to south (WNW at 38.5°N; NW at 36.5°N). A pure shear constant volume solution based on a detailed study of microfaults at 37°20′N indicates a maximum extension of 2.3% in a N61°W direction. These extension directions are remarkably parallel to late Mesozoic to present‐day tectonic extension directions in the Basin and Range province. The pattern of microfaulting demonstrates that the supposedly monolithic Sierran terrane was also affected by the late Cenozoic and possibly earlier regional extension of western North America and provides an independent criterion for determining extensional strain directions.

California↗

USGS2012

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗

USGS2018

Wetlands inform how climate extremes influence surface water expansion and contraction

Effective monitoring and prediction of flood and drought events requires an improved understanding of how and why surface water expansion and contraction in response to climate varies across space. This paper sought to (1) quantify how interannual patterns of surface water expansion and contraction vary spatially across the Prairie Pothole Region (PPR) and adjacent Northern Prairie (NP) in the United States, and (2) explore how landscape characteristics influence the relationship between climate inputs and surface water dynamics. Due to differences in glacial history, the PPR and NP show distinct patterns in regards to drainage development and wetland density, together providing a diversity of conditions to examine surface water dynamics. We used Landsat imagery to characterize variability in surface water extent across 11 Landsat path/rows representing the PPR and NP (images spanned 1985–2015). The PPR not only experienced a 2.6-fold greater surface water extent under median conditions relative to the NP, but also showed a 3.4-fold greater change in surface water extent between drought and deluge conditions. The relationship between surface water extent and accumulated water availability (precipitation minus potential evapotranspiration) was quantified per watershed and statistically related to variables representing hydrology-related landscape characteristics (e.g., infiltration capacity, surface storage capacity, stream density). To investigate the influence stream connectivity has on the rate at which surface water leaves a given location, we modeled stream-connected and stream-disconnected surface water separately. Stream-connected surface water showed a greater expansion with wetter climatic conditions in landscapes with greater total wetland area, but lower total wetland density. Disconnected surface water showed a greater expansion with wetter climatic conditions in landscapes with higher wetland density, lower infiltration and less anthropogenic drainage. From these findings, we can expect that shifts in precipitation and evaporative demand will have uneven effects on surface water quantity. Accurate predictions regarding the effect of climate change on surface water quantity will require consideration of hydrology-related landscape characteristics including wetland storage and arrangement.

Prairie Pothole Regino↗

USGS2015

Coral 13C/12C records of vertical seafloor displacement during megathrust earthquakes west of Sumatra

The recent surge of megathrust earthquakes and tsunami disasters has highlighted the need for a comprehensive understanding of earthquake cycles along convergent plate boundaries. Space geodesy has been used to document recent crustal deformation patterns with unprecedented precision, however the production of long paleogeodetic records of vertical seafloor motion is still a major challenge. Here we show that carbon isotope ratios ( ) in the skeletons of massive Porites corals from west Sumatra record abrupt changes in light exposure resulting from coseismic seafloor displacements. Validation of the method is based on the coral response to uplift (and subsidence) produced by the March 2005 M w 8.6 Nias–Simeulue earthquake, and uplift further south around Sipora Island during a M&sim;8.4 megathrust earthquake in February 1797. At Nias, the average step-change in coral was 0.6±0.1&permil;/m for coseismic displacements of +1.8 m and −0.4 m in 2005. At Sipora, a distinct change in Porites microatoll growth morphology marks coseismic uplift of 0.7 m in 1797. In this shallow water setting, with a steep light attenuation gradient, the step-change in microatoll is 2.3&permil;/m , nearly four times greater than for the Nias Porites . Considering the natural variability in coral skeletal , we show that the lower detection limit of the method is around 0.2 m of vertical seafloor motion. Analysis of vertical displacement for well-documented earthquakes suggests this sensitivity equates to shallow events exceeding M w &sim;7.2 in central megathrust and back-arc thrust fault settings. Our findings indicate that the coral paleogeodesy technique could be applied to convergent tectonic margins throughout the tropical western Pacific and eastern Indian oceans, which host prolific coral reefs, and some of the world's greatest earthquake catastrophes. While our focus here is the link between coral , light exposure and coseismic crustal deformation, the same principles could be used to characterize interseismic strain during earthquake cycles over the last several millennia.

Earth and Planetary Science Letters↗

USGS2008

EAARL coastal topography - Northern Gulf of Mexico

These remotely sensed, geographically referenced elevation measurements of Lidar-derived coastal topography were produced as a collaborative effort between the U.S. Geological Survey (USGS), Florida Integrated Science Center (FISC), St. Petersburg, FL and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, VA. One objective of this research is to create techniques to survey areas for the purposes of geomorphic change studies following major storm events. The USGS Coastal and Marine Geology Program's National Assessment of Coastal Change Hazards project is a multi-year undertaking to identify and quantify the vulnerability of U.S. shorelines to coastal change hazards such as effects of severe storms, sea-level rise, and shoreline erosion and retreat. Airborne Lidar surveys conducted during periods of calm weather are compared to surveys collected following extreme storms in order to quantify the resulting coastal change. Other applications of high-resolution topography include habitat mapping, ecological monitoring, volumetric change detection, and event assessment. The purpose of this project is to provide highly detailed and accurate datasets of the northern Gulf of Mexico coastal areas, acquired on September 19, 2004, immediately following Hurricane Ivan. The datasets are made available for use as a management tool to research scientists and natural resource managers. An innovative airborne Lidar instrument originally developed at the NASA Wallops Flight Facility, and known as the Experimental Airborne Advanced Research Lidar (EAARL), was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532 nanometer) Lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive Lidar, a down-looking RGB (red-green-blue) digital camera, a high-resolution multi-spectral color infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers and an integrated miniature digital inertial measurement unit which provide for sub-meter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine Cessna 310 aircraft, but the instrument may be deployed on a range of light aircraft. A single pilot, a Lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system on September 19, 2004. The survey resulted in the acquisition of 3.2 gigabytes of data. The data were processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of Lidar data in an interactive or batch mode. Modules for pre-survey flight line definition, flight path plotting, Lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is routinely used to create maps that represent submerged or sub-aerial topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation from a point cloud of 'last return' elevations.

Alabama, Florida, Mississippi↗

USGS2022

Science mission requirements for a globally ranging, riserless drilling vessel for U.S. Scientific Ocean Drilling

Through the collection and analysis of shallow and deep subseafloor sediments, rocks, fluids, and life, scientific ocean drilling has enriched our understanding of the complex Earth system. Among other achievements, scientific ocean drilling has documented the history of Earth’s climate, the waxing and waning of polar ice sheets, the past changes in ocean and atmospheric circulation, the existence and function of microbial life in the subseafloor, the compositional variations in Earth’s crust and underlying mantle, and the physical and chemical processes acting at subduction zones, including those associated with tsunamigenic earthquakes. Over the decades, more than 12,000 articles that depend on analyses of scientific ocean drilling samples and geophysical data have been published, many detailing breakthrough contributions to global knowledge about the Earth system. Approximately 45% of these publications were led by U.S.-affiliated authors (International Ocean Discovery Program Publication Services, 2021). Since the mid-1980s, the workhorse of this multidisciplinary, international research effort has been the riserless D/V JOIDES Resolution, operated by Texas A&M University with funding from the U.S. National Science Foundation (NSF). D/V JOIDES Resolution has conducted the vast majority of scientific ocean drilling expeditions and collected most of the scientific cores over that period, including 82% of the expeditions and 93% of the cores in the last decade alone, despite being one of three platforms that is operated within the International Ocean Discovery Program. However, D/V JOIDES Resolution is approaching the end of its useful life. With a strong commitment to continue scientific ocean drilling beyond the end of the current phase, the community developed a document outlining the research frontiers that should be pursued. Exploring Earth by Scientific Ocean Drilling: 2050 Science Framework (Koppers and Coggon, 2020) describes seven scientific strategic objectives that focus on understanding interconnections within the Earth system and five flagship initiatives that integrate these objectives into long-term research efforts that address issues facing society. Additional elements in the 2050 Science Framework, including STEM education, workforce development, technology development, and innovative applications of data analytics, will advance the goals of scientific ocean drilling. Addressing the 2050 Science Framework also requires building partnerships with allied U.S. and international science programs and strengthening existing ones. To implement a significant portion of the 2050 Science Framework, the U.S. scientific community seeks to lease or acquire a newly built, globally ranging, state-of-the art, riserless drilling vessel. The many and varied technical and human resources requirements for successful accomplishment of scientific and educational goals summarized in this document and described in detail in the 2050 Science Framework require broad community input and careful consideration. Following receipt of NSF’s formal Request for Assistance to the United States Science Support Program (USSSP), the U.S. scientific ocean drilling community conducted a one-year exercise to identify its national scientific needs and priorities in order to determine the Science Mission Requirements (SMRs) presented here. This community effort included: (1) a U.S. community-wide survey to identify the specific operational and technical capabilities critical to addressing science in the 2050 Science Framework; (2) a series of online workshops focusing on critical capabilities identified by the survey; and (3) a large in-person workshop to synthesize the results of the survey and the virtual workshops (Appendix 1). The approach was designed to reach as many participants as possible. Overall, 278 survey responses were received from U.S. community members, representing 104 unique institutions from 39 states and the District of Columbia, and 137 unique individuals participated in the workshops (Appendix 2). The results of this effort comprise two classes of SMRs: Foundational Science Mission Requirements and Primary Science Mission Requirements. Foundational SMRs define minimum criteria for a new riserless drilling vessel that can address significant portions of the 2050 Science Framework. Primary SMRs build upon the Foundational SMRs and will create more robust science opportunities and data collection capabilities, will increase progress in addressing the 2050 Science Framework objectives, and will provide more real-time ship-to-shore interaction to improve science productivity, engagement, and outreach. Modern safety and environmental standards, including meeting standards to access protected waters such as exclusive economic zones, extended continental shelves, or high latitudes, while being cognizant of the vessel’s environmental footprint. Safe and efficient operations in global locations and in water depths from 70 m to 6000+ m, with total drill string length of at least 7000 m. High-quality core and data collection from a range of key subseafloor environments. Advanced heave compensation, dynamic positioning, and drill pipe stability. Modern mud and cement/casing systems. Critical onboard measurements for safety, operational decision-making, documentation of ephemeral properties, mission-specific science, and long-term science goals that extend beyond a single expedition. Designated and appropriate space for sample and data preservation. Highly skilled onboard personnel, including technical staff for curation and core handling; support for safety, time-sensitive, and critical shipboard measurements; computer support; equipment and instrument repair; application support; and data assurance. Primary Science Mission Requirements include: Flexible shipboard space for laboratories and on deck to ensure safe, successful implementation of diverse science objectives and operations. Minimizing contamination of recovered samples. Over-the-side capabilities for science-supporting technology (e.g., remotely operated vehicles, water- column sampling, sediment-water interface sampling). Downhole logging and measurements. Expanded borehole observatory capabilities. Reliable and consistent ship-to-shore communications. NSF’s investment in a new globally ranging, riserless drilling vessel will have a powerful economic multiplier effect, including the infusion of additional science support funds in the United States for training and research, the development of new technologies and tools, and the associated scientific and technical workforce development. The skills and knowledge gained through scientific ocean drilling are translatable to careers in fields such as sustainable energy development (e.g., geothermal and offshore wind), carbon sequestration, data management and cyberinfrastructure, biotechnology, communications, science education, policy, hazard mitigation, and environmental management. The United States is a leader in a well-established and internationally collaborative scientific ocean drilling community. A modern, globally ranging, riserless drilling vessel will allow the United States to expand its leadership position, address broad scientific questions that current capabilities preclude, and cultivate equitable international, multidisciplinary collaborations that will ensure scientific ocean drilling’s future success.

Report↗

USGS1996

Type of faulting and orientation of stress and strain as a function of space and time in Kilauea's south flank, Hawaii

Earthquake focal mechanisms of events occurring between 1972 and 1992 in the south flank of Kilauea volcano, Hawaii, are used to infer the state of stress and strain as a function of time and space. We have determined 870 fault plane solutions from P wave first motion polarities for events with magnitudes M L ≥ 2.5 and depth ranging between 6 and 12 km. Faulting is characterized by a mixture of decollement, reverse, and normal faults. Most large earthquakes with magnitude M < 7 slip on reverse faults striking NE at 40° and dipping SE between 60° and 70°. In Hawaii, the earthquakes with M > 7 rupture the decollement plane, since it is the only surface large enough to generate magnitude 7 or larger earthquakes. The percentage of reverse faulting events is high compared to the decollement and normal faulting mechanisms for the period 1972–1983. The percentage of decollement type focal mechanisms becomes dominant after 1983. This pattern of faulting activity suggests that pressure was building up within Kilauea's rift zone prior to the 1983 Puu'Oo eruption. Overall, a single stress orientation with the maximum compressive stress oriented SE perpendicular to the rift and dipping at 45° is compatible with the coeval existence of decollement, reverse, and normal faults. However, in a crustal volume east of longitude 155°10′W, we find a change of the orientation of σ 1 from nearly horizontal to plunging 45° SE occurring in 1979. This stress rotation suggests magma movements within the aseismic part of Kilauea's east rift zone. The strain and stress orientations are coaxial in the south flank except within the volume where the stress rotation is observed. We observe a change in the relationship between stress and strain directions caused either by the shifting of seismic activity from reverse faults to decollements, while stress stays constant, or by a rotation of stress, while strain remains constant. Assuming that the model of a noncohesive Coulomb wedge is appropriate for Kilauea's south flank, we find that high pore pressures are prevalent along the decollement and within the wedge for a coefficient of friction equal to 0.85.

Hawaii↗

USGS2024

Automatic identification and quantification of volcanic hotspots in Alaska using HotLINK: The hotspot learning and identification network

An increase in volcanic thermal emissions can indicate subsurface and surface processes that precede, or coincide with, volcanic eruptions. Space-borne infrared sensors can detect hotspots—defined here as localized volcanic thermal emissions—in near-real-time. However, automatic hotspot detection systems are needed to efficiently analyze the large quantities of data produced. While hotspots have been automatically detected for over 20 years with simple thresholding algorithms, new computer vision technologies, such as convolutional neural networks (CNNs), can enable improved detection capabilities. Here we introduce HotLINK: the Hotspot Learning and Identification Network, a CNN trained to detect hotspots with a dataset of −3,800 satellite-based, Visible Infrared Imaging Radiometer Suite (VIIRS) images from Mount Veniaminof and Mount Cleveland volcanoes, Alaska. We find that our model achieves an accuracy of 96% (F1-score 0.92) when evaluated on −1,700 unseen images from the same volcanoes, and 95% (F1-score 0.67) when evaluated on −3,000 images from six additional Alaska volcanoes (Augustine Volcano, Bogoslof Island, Okmok Caldera, Pavlof Volcano, Redoubt Volcano, Shishaldin Volcano). In comparison with an existing threshold-based hotspot detection algorithm, MIROVA (Coppola et al., Geological Society, London, Special Publications, 2016, 426, 181–205), our model detects 22% more hotspots and produces 12% fewer false positives. Additional testing on −700 labeled Moderate Resolution Imaging Spectroradiometer (MODIS) images from Mount Veniaminof demonstrates that our model is applicable to this sensor’s data as well, achieving an accuracy of 98% (F1-score 0.95). We apply HotLINK to 10 years of VIIRS data and 22 years of MODIS data for the eight aforementioned Alaska volcanoes and calculate the radiative power of detected hotspots. From these time series we find that HotLINK accurately characterizes background and eruptive periods, similar to MIROVA, but also detects more subtle warming signals, potentially related to volcanic unrest. We identify three advantages to our model over its predecessors: 1) the ability to detect more subtle volcanic hotspots and produce fewer false positives, especially in daytime images; 2) probabilistic predictions provide a measure of detection confidence; and 3) its transferability, i.e., the successful application to multiple sensors and multiple volcanoes without the need for threshold tuning, suggesting the potential for global application.

Alaska↗