USGS Science⌕ Search

SEARCH · USGS Science

Results for “The Earth Observer”

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 865 records · Page 48Linked to original sources
USGS2006

The sand seas of titan: Cassini RADAR observations of longitudinal dunes

The most recent Cassini RADAR images of Titan show widespread regions (up to 1500 kilometers by 200 kilometers) of near-parallel radar-dark linear features that appear to be seas of longitudinal dunes similar to those seen in the Namib desert on Earth. The Ku-band (2.17-centimeter wavelength) images show ∼100-meter ridges consistent with duneforms and reveal flow interactions with underlying hills. The distribution and orientation of the dunes support a model of fluctuating surface winds of ∼0.5 meter per second resulting from the combination of an eastward flow with a variable tidal wind. The existence of dunes also requires geological processes that create sand-sized (100- to 300-micrometer) particulates and a lack of persistent equatorial surface liquids to act as sand traps.

Science↗
USGS2001

hypoDD-A Program to Compute Double-Difference Hypocenter Locations

HypoDD is a Fortran computer program package for relocating earthquakes with the double-difference algorithm of Waldhauser and Ellsworth (2000). This document provides a brief introduction into how to run and use the programs ph2dt and hypoDD to compute double-difference (DD) hypocenter locations. It gives a short overview of the DD technique, discusses the data preprocessing using ph2dt, and leads through the earthquake relocation process using hypoDD. The appendices include the reference manuals for the two programs and a short description of auxiliary programs and example data. Some minor subroutines are presently in the c language, and future releases will be in c. Earthquake location algorithms are usually based on some form of Geiger’s method, the linearization of the travel time equation in a first order Taylor series that relates the difference between the observed and predicted travel time to unknown adjustments in the hypocentral coordinates through the partial derivatives of travel time with respect to the unknowns. Earthquakes can be located individually with this algorithm, or jointly when other unknowns link together the solutions to indivdual earthquakes, such as station corrections in the joint hypocenter determination (JHD) method, or the earth model in seismic tomography. The DD technique (described in detail in Waldhauser and Ellsworth, 2000) takes advantage of the fact that if the hypocentral separation between two earthquakes is small compared to the event-station distance and the scale length of velocity heterogeneity, then the ray paths between the source region and a common station are similar along almost the entire ray path (Fréchet, 1985; Got et al., 1994). In this case, the difference in travel times for two events observed at one station can be attributed to the spatial offset between the events with high accuracy. DD equations are built by differencing Geiger’s equation for earthquake location. In this way, the residual between observed and calculated travel-time difference (or double-difference) between two events at a common station are a related to adjustments in the relative position of the hypocenters and origin times through the partial derivatives of the travel times for each event with respect to the unknown. HypoDD calculates travel times in a layered velocity model (where velocity depends only on depth) for the current hypocenters at the station where the phase was recorded. The double-difference residuals for pairs of earthquakes at each station are minimized by weighted least squares using the method of singular value decomposition (SVD) or the conjugate gradients method (LSQR, Paige and Saunders, 1982). Solutions are found by iteratively adjusting the vector difference between nearby hypocentral pairs, with the locations and partial derivatives being updated after each iteration. Details about the algorithm can be found in Waldhauser and Ellsworth (2000). When the earthquake location problem is linearized using the double-difference equations, the common mode errors cancel, principally those related to the receiver-side structure. Thus we avoid the need for station corrections or high-accuracy of predicted travel times for the portion of the raypath that lies outside the focal volume. This approach is especially useful in regions with a dense distribution of seismicity, i.e. where distances between neighboring events are only a few hundred meters. The improvement of double-difference locations over ordinary JHD locations is shown in Figure 1 for about 10,000 earthquakes that occurred during the 1997 seismic crisis in the Long Valley caldera, California. While the JHD locations (left panel) show a diffuse picture of the seismicity, double-difference locations (right panel) bring structural details such as the location of active fault planes into sharp focus.

Open-File Report↗
USGS2025

Petrogenesis and mineralization potential of spinifex komatiitic basalts in the Bradley Peak greenstone terrane, Wyoming Province

Komatiitic volcanic rocks are important hosts of Ni sulfide mineralization and record early Earth evolution; however, those in the well-studied Archean Wyoming Province have received little attention. Here, we elucidate the timing and petrogenesis of the Bradley Peak komatiitic volcanic rocks using field and textural observations, geochronology, and geochemistry. Detrital and igneous zircon U-Pb ages for two samples from previously undated units support published age determinations, placing the eruption age at 2.72 Ga. Stratigraphy of the volcanic flows was mapped and 36 samples including cumulates, greenschists, and spinifex-textured rocks were collected. Whole-rock geochemistry was used to classify the spinifex-textured samples as Al-undepleted komatiitic basalts (11–17 wt% MgO). Platinum-group element concentrations (n = 25) are like those in global Al-undepleted komatiitic basalts, and PGE/Ti ratios do not indicate the volcanic flows likely host sulfide mineralization. Initial ε Nd values of −0.5 to +4.7 (n = 16), indicate that these lavas were derived from a depleted mantle source and have negligible evolved crust contamination. The primary magma to the komatiitic basalt flows is estimated to have had 19 wt% MgO and be derived from ∼15 to 25 % mantle partial melting at 3–4 GPa. Trace element chemistry and thermodynamic modeling suggest the primary melt assimilated local banded iron formation. Although the Bradley Peak komatiitic basalts do not contain positive evidence of magmatic sulfide deposits, depleted Au in the flows suggests they could be source rocks for nearby orogenic gold deposits.

Wyoming↗
USGS2019

Dome formation on Ceres by sold-state flow analogous to terrestrial salt tectonics

The dwarf planet Ceres’s outer crust is a complex, heterogeneous mixture of ice, clathrates, salts and silicates. Numerous large domes on Ceres’s surface indicate a degree of geological activity. These domes have been attributed to cryovolcanism, but that is difficult to reconcile with Ceres’s small size and lack of long-lived heat sources. Here we alternatively propose that Ceres’s domes form by solid-state flow within the compositionally heterogeneous crust, a mechanism directly analogous to salt tectonics on Earth. We use numerical simulations to illustrate that differential loading of a crust with compositional heterogeneity on a scale of tens of kilometres can produce dome-like features of scale similar to those observed. The mechanism requires the presence of low-viscosity and low-density, possibly ice-rich, material in the upper 1–10 km of the subsurface. Such substantial regional heterogeneity in Ceres’s crustal composition is consistent with observations from the National Aeronautics and Space Administration’s Dawn mission. We conclude that deformation analogous to that in terrestrial salt tectonics is a viable alternative explanation for the observed surface morphologies, and is consistent with Ceres being both cold and geologically active.

Nature Geoscience↗
USGS2013

Projected future changes in vegetation in western North America in the 21st century

Rapid and broad-scale forest mortality associated with recent droughts, rising temperature, and insect outbreaks has been observed over western North America (NA). Climate models project additional future warming and increasing drought and water stress for this region. To assess future potential changes in vegetation distributions in western NA, the Community Earth System Model (CESM) coupled with its Dynamic Global Vegetation Model (DGVM) was used under the future A2 emissions scenario. To better span uncertainties in future climate, eight sea surface temperature (SST) projections provided by phase 3 of the Coupled Model Intercomparison Project (CMIP3) were employed as boundary conditions. There is a broad consensus among the simulations, despite differences in the simulated climate trajectories across the ensemble, that about half of the needleleaf evergreen tree coverage (from 24% to 11%) will disappear, coincident with a 14% (from 11% to 25%) increase in shrubs and grasses by the end of the twenty-first century in western NA, with most of the change occurring over the latter half of the twenty-first century. The net impact is a ~6 GtC or about 50% decrease in projected ecosystem carbon storage in this region. The findings suggest a potential for a widespread shift from tree-dominated landscapes to shrub and grass-dominated landscapes in western NA because of future warming and consequent increases in water deficits. These results highlight the need for improved process-based understanding of vegetation dynamics, particularly including mortality and the subsequent incorporation of these mechanisms into earth system models to better quantify the vulnerability of western NA forests under climate change.

Journal of Climate↗
USGS2012

Casual instrument corrections for short-period and broadband seismometers

Of all the filters applied to recordings of seismic waves, which include source, path, and site effects, the one we know most precisely is the instrument filter. Therefore, it behooves seismologists to accurately remove the effect of the instrument from raw seismograms. Applying instrument corrections allows analysis of the seismogram in terms of physical units (e.g., displacement or particle velocity of the Earth’s surface) instead of the output of the instrument (e.g., digital counts). The instrument correction can be considered the most fundamental processing step in seismology since it relates the raw data to an observable quantity of interest to seismologists. Complicating matters is the fact that, in practice, the term “instrument correction” refers to more than simply the seismometer. The instrument correction compensates for the complete recording system including the seismometer, telemetry, digitizer, and any anti‐alias filters. Knowledge of all these components is necessary to perform an accurate instrument correction. The subject of instrument corrections has been covered extensively in the literature (Seidl, 1980; Scherbaum, 1996). However, the prospect of applying instrument corrections still evokes angst among many seismologists—the authors of this paper included. There may be several reasons for this. For instance, the seminal paper by Seidl (1980) exists in a journal that is not currently available in electronic format and cannot be accessed online. Also, a standard method for applying instrument corrections involves the programs TRANSFER and EVALRESP in the Seismic Analysis Code (SAC) package (Goldstein et al. , 2003). The exact mathematical methods implemented in these codes are not thoroughly described in the documentation accompanying SAC.

Seismological Research Letters↗
USGS1999

Earthquake probabilities in the San Francisco Bay region: 2000 to 2030 - A summary of findings

The San Francisco Bay region sits astride a dangerous “earthquake machine,” the tectonic boundary between the Pacific and North American Plates. The region has experienced major and destructive earthquakes in 1838, 1868, 1906, and 1989, and future large earthquakes are a certainty. The ability to prepare for large earthquakes is critical to saving lives and reducing damage to property and infrastructure. An increased understanding of the timing, size, location, and effects of these likely earthquakes is a necessary component in any effective program of preparedness. This study reports on the probabilities of occurrence of major earthquakes in the San Francisco Bay region (SFBR) for the three decades 2000 to 2030. The SFBR extends from Healdsberg on the northwest to Salinas on the southeast and encloses the entire metropolitan area, including its most rapidly expanding urban and suburban areas. In this study a “major” earthquake is defined as one with M≥6.7 (where M is moment magnitude). As experience from the Northridge, California (M6.7, 1994) and Kobe, Japan (M6.9, 1995) earthquakes has shown us, earthquakes of this size can have a disastrous impact on the social and economic fabric of densely urbanized areas. To reevaluate the probability of large earthquakes striking the SFBR, the U.S. Geological Survey solicited data, interpretations, and analyses from dozens of scientists representing a wide crosssection of the Earth-science community (Appendix A). The primary approach of this new Working Group (WG99) was to develop a comprehensive, regional model for the long-term occurrence of earthquakes, founded on geologic and geophysical observations and constrained by plate tectonics. The model considers a broad range of observations and their possible interpretations. Using this model, we estimate the rates of occurrence of earthquakes and 30-year earthquake probabilities. Our study considers a range of magnitudes for earthquakes on the major faults in the region—an innovation over previous studies of the SFBR that considered only a small number of potential earthquakes of fixed magnitude.

California↗
USGS2017

Climatic history of the northeastern United States during the past 3000 years

Many ecosystem processes that influence Earth system feedbacks, including vegetation growth, water and nutrient cycling, and disturbance regimes, are strongly influenced by multi-decadal to millennial-scale variations in climate that cannot be captured by instrumental climate observations. Paleoclimate information is therefore essential for understanding contemporary ecosystems and their potential trajectories under a variety of future climate conditions. With the exception of fossil pollen records, there are a limited number of northeastern US (NE US) paleoclimate archives that can provide constraints on its temperature and hydroclimate history. Moreover, the records that do exist have not been considered together. Tree-ring data indicate that the 20th century was one of the wettest of the past 500 years in the eastern US (Pederson et al., 2014), and lake-level records suggest it was one of the wettest in the Holocene (Newby et al., 2014); how such results compare with other available data remains unclear, however. Here we conduct a systematic review, assessment, and comparison of paleotemperature and paleohydrological proxies from the NE US for the last 3000 years. Regional temperature reconstructions are consistent with the long-term cooling trend (1000 BCE–1700 CE) evident in hemispheric-scale reconstructions, but hydroclimate reconstructions reveal new information, including an abrupt transition from wet to dry conditions around 550–750 CE. NE US paleo data suggest that conditions during the Medieval Climate Anomaly were warmer and drier than during the Little Ice Age, and drier than today. There is some evidence for an acceleration over the past century of a longer-term wetting trend in the NE US, and coupled with the abrupt shift from a cooling trend to a warming trend from increased greenhouse gases, may have wide-ranging implications for species distributions, ecosystem dynamics, and extreme weather events. More work is needed to gather paleoclimate data in the NE US, make inter-proxy comparisons, and improve estimates of uncertainty in the reconstructions.

Climate of the Past↗
USGS1999

Climatic impact of glacial cycle polar motion: Coupled oscillations of ice sheet mass and rotation pole position

Precessional motion of Earth's rotation axis relative to its orbit is a well-known source of long-period climatic variation. It is less well appreciated that growth and decay of polar ice sheets perturb the symmetry of the global mass distribution enough that the geographic location of the rotation axis will change by at least 15 km and possibly as much as 100 km during a single glacial cycle. This motion of the pole will change the seasonal and latitudinal pattern of temperatures. We present calculations, based on a diurnal average energy balance, which compare the summer and winter temperature anomalies due to a 1° decrease in obliquity with those due to a 1° motion of the rotation pole toward Hudson Bay. Both effects result in peak temperature perturbations of about 1° Celsius. The obliquity change primarily influences the amplitude of the seasonal cycle, while the polar motion primarily changes the annual mean temperatures. The polar motion induced temperature anomaly is such that it will act as a powerful negative feedback on ice sheet growth. We also explore the evolution of the coupled system composed of ice sheet mass and pole position. Oscillatory solutions result from the conflicting constraints of rotational and thermal stability. A positive mass anomaly on an otherwise featureless Earth is in rotational equilibrium only at the poles or the equator. The two polar equilibria are rotationally unstable, and the equatorial equilibrium, though rotationally stable, is thermally unstable. We find that with a plausible choice for the strength of coupling between the thermal and rotational systems, relatively modest external forcing can produce significant response at periods of 10 4 –10 6 years, but it strongly attenuates polar motion at longer periods. We suggest that these coupled oscillations may contribute to the observed dominance of 100 kyr glacial cycles since the mid-Pleistocene and will tend to stabilize geographic patterns that are suitable to glaciations.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

In-situ observations of the physical properties of the Martian surface

The physical properties of rocks and soils on the surface of Mars have been investigated by several landed spacecraft. Studies of these physical properties constrain interpretation of Martian geologic processes and provide engineering data for future mission planning. As on Earth, these properties vary considerably from place to place, and provide constraints on the origin and evolution of the surface materials. Martian soils commonly have thin surface crusts that may be caused by salts cementing grains together. Estimates of soil physical properties at the various landing sites are generally comparable, but rather uncertain in many cases. Rock physical properties, based on abrasion by the Mars Exploration Rover (MER) Rock Abrasion Tool (RAT) and other experiments, vary widely.

Book chapter↗
USGS2018

The Polar WRF downscaled historical and projected twenty-first century climate for the coast and foothills of Arctic Alaska

Climate change is most pronounced in the northern high latitude region. Yet, climate observations are unable to fully capture regional-scale dynamics due to the sparse weather station coverage, which limits our ability to make reliable climate-based assessments. A set of simulated data products was therefore developed for the North Slope of Alaska through a dynamical downscaling approach. The polar-optimized Weather Research and Forecast (Polar WRF) model was forced by three sources: The ERA-interim reanalysis data (for 1979–2014), the Community Earth System Model 1.0 (CESM1.0) historical simulation (for 1950–2005), and the CESM1.0 projected (for 2006–2100) simulations in two Representative Concentration Pathways (RCP4.5 and RCP8.5) scenarios. Climatic variables were produced in a 10-km grid spacing and a 3-h interval. The ERA-interim forced WRF (ERA-WRF) proves the value of dynamical downscaling, which yields more realistic topographical-induced precipitation and air temperature, as well as corrects underestimations in observed precipitation. In summary, dry and cold biases to the north of the Brooks Range are presented in ERA-WRF, while CESM forced WRF (CESM-WRF) holds wet and warm biases in its historical period. A linear scaling method allowed for an adjustment of the biases, while keeping the majority of the variability and extreme values of modeled precipitation and air temperature. CESM-WRF under RCP 4.5 scenario projects smaller increase in precipitation and air temperature than observed in the historical CESM-WRF product, while the CESM-WRF under RCP 8.5 scenario shows larger changes. The fine spatial and temporal resolution, long temporal coverage, and multi-scenario projections jointly make the dataset appropriate to address a myriad of physical and biological changes occurring on the North Slope of Alaska.

Alaska↗
USGS2014

Residual shear strength variability as a primary control on movement of landslides reactivated by earthquake-induced ground motion: Implications for coastal Oregon, U.S.

Most large seismogenic landslides are reactivations of preexisting landslides with basal shear zones in the residual strength condition. Residual shear strength often varies during rapid displacement, but the response of residual shear zones to seismic loading is largely unknown. We used a ring shear apparatus to perform simulated seismic loading tests, constant displacement rate tests, and tests during which shear stress was gradually varied on specimens from two landslides to improve understanding of coseismic landslide reactivation and to identify shear strength models valid for slow gravitational failure through rapid coseismic failure. The landslides we studied represent many along the Oregon, U.S., coast. Seismic loading tests resulted in (1) catastrophic failure involving unbounded displacement when stresses represented those for the existing landslides and (2) limited to unbounded displacement when stresses represented those for hypothetical dormant landslides, suggesting that coseismic landslide reactivation may be significant during future great earthquakes occurring near the Oregon Coast. Constant displacement rate tests indicated that shear strength decreased exponentially during the first few decimeters of displacement but increased logarithmically with increasing displacement rate when sheared at 0.001 cm s −1 or greater. Dynamic shear resistance estimated from shear strength models correlated well with stresses observed during seismic loading tests, indicating that displacement rate and amount primarily controlled failure characteristics. We developed a stress-based approach to estimate coseismic landslide displacement that utilizes the variable shear strength model. The approach produced results that compared favorably to observations made during seismic loading tests, indicating its utility for application to landslides.

Oregon↗
USGS1990

Textural development of clayey and quartzofeldspathic fault gouges relative to their sliding behavior

Many of the secondary fault structures developed during triaxial friction experiments have been generally correlated with the structures of natural fault zones. Therefore, any physical differences that can be found between laboratory samples that slide stably and those that show stick-slip motion may help to identify the cause of earthquakes. We have examined petrographically the run products of many triaxial friction experiments using clayey and quartzofeldspathic gouges, which comprise the principal types of natural fault gouge material. The examined samples were tested under a wide range of temperature, confining and fluid pressure, and velocity conditions. The clayey and quartzofeldspathic gouges show some textural differences, owing to their different mineral contents and grain sizes and shapes. In the clayey gouges, for example, a clay mineral fabric and kink band sets are commonly developed, whereas in the quartzofeldspathic gouges fracturing and crushing of the predominately quartz and feldspar grains are important processes. For both types of gouge, however, and whatever the pressure-temperature-velocity conditions of the experiments, the transition from stable sliding to stick-slip motion is correlated with: (i) a change from pervasive deformation of the gouge layer to localized slip in subsidiary shears; and (ii) an increase in the angle betweem the shears that crosscut the gouge layer (Riedel shears) and ones that form along the gouge-rock cylinder boundaries (boundary shears). This suggests that the localization of shear within a fault zone combined with relatively high Riedel-shear angles are somehow connected with earthquakes. Secondary fracture sets similar to Riedel shears have been identified at various scales in major strike-slip faults such as the San Andreas of the western United States (Wallace, 1973) and the Luhuo and Fuyun earthquake faults of China (Deng and Zhang, 1984; Deng et al. , 1986). The San Andreas also contains locked and creeping sections that correspond to the stick-slip and stably sliding experimental samples, respectively. We plan to study the physical structure of the San Andreas fault, to see if the experimentally observed differences related to sliding behavior can also be distinguished in the field.

Physics and Chemistry of the Earth↗
USGS1998

Poroelastic rebound along the Landers 1992 earthquake surface rupture

Maps of surface displacement following the 1992 Landers, California, earthquake, generated by interferometric processing of ERS-1 synthetic aperture radar (SAR) images, reveal effects of various postseismic deformation processes along the 1992 surface rupture. The large-scale pattern of the postseismic displacement field includes large lobes, mostly visible on the west side of the fault, comparable in shape with the lobes observed in the coseismic displacement field. This pattern and the steep displacement gradient observed near the Emerson-Camp Rock fault cannot be simply explained by afterslip on deep sections of the 1992 rupture. Models show that horizontal slip occurring on a buried dislocation in a Poisson's material produces a characteristic quadripole pattern in the surface displacement field with several centimeters of vertical motion at distances of 10–20 km from the fault, yet this pattern is not observed in the postseismic interferograms. As previously proposed to explain local strain in the fault step overs [ Peltzer et al ., 1996b], we argue that poroelastic rebound caused by pore fluid flow may also occur over greater distances from the fault, compensating the vertical ground shift produced by fault afterslip. Such a rebound is explained by the gradual change of the crustal rocks' Poisson's ratio value from undrained (coseismic) to drained (postseismic) conditions as pore pressure gradients produced by the earthquake dissipate. Using the Poisson's ratio values of 0.27 and 0.31 for the drained and undrained crustal rocks, respectively, elastic dislocation models show that the combined contributions of afterslip on deep sections of the fault and poroelastic rebound can account for the range change observed in the SAR data and the horizontal displacement measured at Global Positioning System (GPS) sites along a 60-km-long transect across the Emerson fault [ Savage and Svarc , 1997]. Using a detailed surface slip distribution on the Homestead Valley, Kickapoo, and Johnson Valley faults, we modeled the poroelastic rebound in the Homestead Valley pull apart. A Poisson's ratio value of 0.35 for the undrained gouge rocks in the fault zone is required to account for the observed surface uplift in the 3.5 years following the earthquake. This large value implies a seismic velocity ratio V p /V s of 2.1, consistent with the observed low V s values of fault zone guided waves at shallow depth [ Li et al ., 1997]. The SAR data also reveal postseismic creep along shallow patches of the Eureka Peak and Burnt Mountain faults with a characteristic decay time of 0.8 years. Coseismic, dilatant hardening (locking process) followed by post-seismic, pore pressure controlled fault creep provide a plausible mechanism to account for the decay time of the observed slip rate along this section of the fault.

Journal of Geophysical Research B: Solid Earth↗
USGS2023

Steady-state forms of channel profiles shaped by debris flow and fluvial processes

Debris flows regularly traverse bedrock channels that dissect steep landscapes, but our understanding of bedrock erosion by debris flows and their impact on steepland morphology is still rudimentary. Quantitative models of steep bedrock channel networks are based on geomorphic transport laws designed to represent erosion by water-dominated flows. To quantify the impact of debris flow erosion on steep channel network form, it is first necessary to develop methods to estimate spatial variations in bulk debris flow properties (e.g., flow depth, velocity) throughout the channel network that can be integrated into landscape evolution models. Here, we propose and evaluate two methods to estimate spatial variations in bulk debris flow properties along the length of a channel profile. We incorporate both methods into a model designed to simulate the evolution of longitudinal channel profiles that evolve in response to debris flow and fluvial processes. To explore this model framework, we propose a general family of debris flow erosion laws where erosion rate is a function of debris flow depth and channel slope. Model results indicate that erosion by debris flows can explain the occurrence of a scaling break in the slope–area curve at low-drainage areas and that upper-network channel morphology may be useful for inferring catchment-averaged erosion rates in quasi-steady landscapes. Validating specific forms of a debris flow incision law, however, would require more detailed model–data comparisons in specific landscapes where input parameters and channel morphometry can be better constrained. Results improve our ability to interpret topographic signals within steep channel networks and identify observational targets critical for constraining a debris flow incision law.

Earth Surface Dynamics↗
USGS2015

Crustal stress and structure at Kīlauea Volcano inferred from seismic anisotropy

Seismic anisotropy, measured through shear wave splitting (SWS) analysis, can be indicative of the state of stress in Earth's crust. Changes in SWS at Kīlauea Volcano, Hawai‘i, associated with the onset of summit eruptive activity in 2008 hint at the potential of the technique for tracking volcanic activity. To use SWS observations as a monitoring tool, however, it is important to understand the cause of seismic anisotropy at the volcano throughout the eruptive cycle. To address this need, we analyzed SWS results from across Kīlauea in combination with macroscopic surface structures (mapped fractures, faults, and fissures) and stress orientations inferred from fault plane solutions. Seismic anisotropy seems to be due to pervasive aligned structures in most regions of the volcano. The upper East and Southwest Rift Zones, however, show a bimodality in stress and SWS, suggesting a stress discontinuity with depth, perhaps related to magma conduits that trend obliquely to the dominant structure. Other areas in and around Kīlauea Caldera display principal stresses of similar magnitudes, indicating that small stress perturbations can rotate the maximum horizontal compressive stress direction by up to 90 ° . In these locations, static structures generally control SWS, but dynamic conditions due to magmatic activity can override the structural control. Monitoring of SWS may therefore provide important signs of impending volcanism.

Hawaii↗
USGS2026

Modeling the seasonality of wind-driven hydrocarbon waves in Titan’s polar lakes

Titan, the only body in the solar system aside from Earth with standing liquids on its surface, has polar hydrocarbon lakes and seas. As Titan’s atmosphere generates light winds, there should be waves on the surface of these lakes and seas, yet, direct wave observations are scant. We introduce and use PlanetWaves , an open source 4D spectral wave model, to study Titan’s waves and create seasonal maps of wave shape and propagation on Ontario Lacus and Ligeia Mare. Titan’s modeled waves grow up to 30 times larger than terrestrial waves for the same wind speed, are seasonally present and are largest in the spring and summer when winds are strongest. Average daily winds almost never exceed the wave generation threshold of 0.5–0.7 m/s. Average storm winds (∼1.5 m/s) generate waves 15–48 cm in height with a period ranging 6–10.5 s while maximum storm winds (∼4 m/s) generate waves 2.7–3.2 m in height with a period up to 32 s. Titan’s waves become fetch-independent at ∼40 km for average storm winds occurring ∼1% of a Titan year and ∼100 kilometers for maximum storm winds occurring 2-3 times per Titan decade. On Ontario Lacus, storm winds blow nearly parallel to the eastern shore, potentially driving wave modification of the smooth eastern shoreline. On Ligeia Mare, waves rarely propagate toward a hypothesized wave modified shoreline suggesting that another process, such as tectonics, may contribute to a straight shoreline morphology.

JGR Planets↗
USGS1997

Correlated errors in geodetic time series: Implications for time-dependent deformation

Analysis of frequent trilateration observations from the two-color electronic distance measuring networks in California demonstrate that the noise power spectra are dominated by white noise at higher frequencies and power law behavior at lower frequencies. In contrast, Earth scientists typically have assumed that only white noise is present in a geodetic time series, since a combination of infrequent measurements and low precision usually preclude identifying the time-correlated signature in such data. After removing a linear trend from the two-color data, it becomes evident that there are primarily two recognizable types of time-correlated noise present in the residuals. The first type is a seasonal variation in displacement which is probably a result of measuring to shallow surface monuments installed in clayey soil which responds to seasonally occurring rainfall; this noise is significant only for a small fraction of the sites analyzed. The second type of correlated noise becomes evident only after spectral analysis of line length changes and shows a functional relation at long periods between power and frequency of 1/ƒ α , where ƒ is frequency and α≈2. With α=2, this type of correlated noise is termed random-walk noise, and its source is mainly thought to be small random motions of geodetic monuments with respect to the Earth's crust, though other sources are possible. Because the line length changes in the two-color networks are measured at irregular intervals, power spectral techniques cannot reliably estimate the level of 1/ƒ α noise. Rather, we also use here a maximum likelihood estimation technique which assumes that there are only two sources of noise in the residual time series (white noise and random-walk noise) and estimates the amount of each. From this analysis we find that the random-walk noise level averages about 1.3 mm/√yr and that our estimates of the white noise component confirm theoretical limitations of the measurement technique. In addition, the seasonal noise can be as large as 3 mm in amplitude but typically is less than 0.5 mm. Because of the presence of random-walk noise in these time series, modeling and interpretation of the geodetic data must account for this source of error. By way of example we show that estimating the time-varying strain tensor (a form of spatial averaging) from geodetic data having both random-walk and white noise error components results in seemingly significant variations in the rate of strain accumulation; spatial averaging does reduce the size of both noise components but not their relative influence on the resulting strain accumulation model.

Journal of Geophysical Research B: Solid Earth↗