USGS Science⌕ Search

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.

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

At least 811 records · Page 45Linked to original sources
USGS2021

Projecting climate dependent coastal flood risk with a hybrid statistical dynamical model

Numerical models for tides, storm surge, and wave runup have demonstrated ability to accurately define spatially varying flood surfaces. However these models are typically too computationally expensive to dynamically simulate the full parameter space of future oceanographic, atmospheric, and hydrologic conditions that will constructively compound in the nearshore to cause both extreme event and nuisance flooding during the 21st century. A surrogate modeling framework of waves, winds, and tides is developed in this study to efficiently predict spatially varying nearshore and estuarine water levels contingent on any combination of offshore forcing conditions. The surrogate models are coupled with a time-dependent stochastic climate emulator that provides efficient downscaling for hypothetical iterations of offshore conditions. Together, the hybrid statistical-dynamical framework can assess present day and future coastal flood risk, including the chronological characteristics of individual flood and wave-induced dune overtopping events and their changes into the future. The framework is demonstrated at Naval Base Coronado in San Diego, CA, utilizing the regional Coastal Storm Modeling System (CoSMoS; composed of Delft3D and XBeach) as the dynamic simulator and Gaussian process regression as the surrogate modeling tool. Validation of the framework uses both in-situ tide gauge observations within San Diego Bay, and a nearshore cross-shore array deployment of pressure sensors in the open beach surf zone. The framework reveals the relative influence of large-scale climate variability on future coastal flood resilience metrics relevant to the management of an open coast artificial berm, as well as the stochastic nature of future total water levels.

California↗
USGS2023

Accuracy of finite fault slip estimates in subduction zone regions with topographic Green's functions and seafloor geodesy

Until recently, the lack of seafloor geodetic instrumentation and the use of unrealistically simple, half-space based forward models have resulted in poor resolution of near-trench slip in subduction zone settings. Here, we use a synthetic framework to investigate the impact of topography and geodetic data distribution on coseismic slip estimates in various subduction zone settings. We calculate surface displacements in two synthetic topographic domains that have topography similar to that of Chile and Japan, respectively. We then attempt to image target slip distributions by using a Bayesian approach to solve for slip with two sets of Green's functions—one that accounts for topography and one that does not—and five sets of 50 or more observation points selected from the synthetic surface displacements. Three of these sets of observation points are entirely onland, and two include 5–10 seafloor geodetic sites. We find that the use of topographic Green's functions always improves inferred slip models, and with seafloor geodetic data, it enables an almost perfect recovery of a target slip model, even in the near-trench region. Critically, our results demonstrate that it would be impossible for non-topographic Green's functions to properly recover the true slip distribution, particularly in the near-trench region. We also perform a parameter study with approximately 4,000 slip models estimated using a least-square approach, and find that topographic Green's functions yield significantly more accurate slip models in cases where good data (well distributed and reasonably dense) are available, even in the absence of seafloor geodetic sites.

Journal of Geophysical Research: Solid Earth↗
USGS1997

Post seismic deformation associated with the 1992 Mω = 7.3 Landers earthquake, southern California

Following the 1992 M ω =7.3 Landers earthquake, a linear array of 10 geodetic monuments at roughly 5‐km spacing was established across the Emerson fault segment of the Landers rupture. The array trends perpendicular to the local strike of the fault segment and extends about 30 km on either side of it. The array was surveyed by Global Positioning System 0.034, 0.048, 0.381, 1.27, 1.88, 2.60, and 3.42 years after the Landers earthquake to measure both the spatial and temporal character of the postearthquake relaxation. The temporal behavior is described roughly by a short‐term (decay time 84±23 days) exponential relaxation superimposed upon an apparently linear trend. Because the linear trend represents motions much more rapid than the observed preseismic motions, we attribute that trend to a slower (decay time greater than 5 years) postseismic relaxation, the curvature of which cannot be resolved in the short run (3.4 years) of postseismic data. About 100 mm of right‐lateral displacement and 50 mm of fault‐normal displacement accumulated across the geodetic array in the 3.4‐year interval covered by the postseismic surveys. Those displacements are attributed to postseismic, right‐lateral slip in the depth interval 10 to 30 km on the downward extension of the rupture trace. The right‐lateral slip amounted to about 1 m directly beneath the geodetic array, and the fault‐normal displacement is apparently primarily a consequence of the curvature of the rupture. These conclusions are based upon dislocation models fit to the observed deformation. However, no dislocation model was found with rms residuals as small as the expected observational error.

California↗
USGS2013

A new model for the growth of basaltic shields based on deformation of Fernandina volcano, Galápagos Islands

Space-geodetic measurements of surface deformation produced by the most recent eruptions at Fernandina – the most frequently erupting volcano in the Galápagos Archipelago – reveal that all have initiated with the intrusion of subhorizontal sills from a shallow magma reservoir. This includes eruptions from fissures that are oriented both radially and circumferentially with respect to the summit caldera. A Synthetic Aperture Radar (SAR) image acquired 1–2 h before the start of a radial fissure eruption in 2009 captures one of these sills in the midst of its propagation toward the surface. Galápagos eruptive fissures of all orientations have previously been presumed to be fed by vertical dikes, and this assumption has guided models of the origin of the eruptive fissure geometry and overall development of the volcanoes. Our findings allow us to reinterpret the internal structure and evolution of Galápagos volcanoes and of similar basaltic shields. Furthermore, we note that stress changes generated by the emplacement of subhorizontal sills feeding one type of eruption may control the geometry of subsequent eruptive fissures. Specifically, circumferential fissures tend to open within areas uplifted by sill intrusions that initiated previous radial fissure eruptions. This mechanism provides a possible explanation for the pattern of eruptive fissures that characterizes all the western Galápagos volcanoes, as well as the alternation between radial and circumferential fissure eruptions at Fernandina. The same model suggests that the next eruption of Fernandina will be from a circumferential fissure in the area uplifted by the 2009 sill intrusion, just southwest of the caldera rim.

Galápagos Islands↗
USGS2023

Geochemical evidence for evolving Proterozoic crustal thickness and orogenic styles in southwestern Laurentia

It has long been challenging for researchers to track the crustal thickness and mode(s) of crustal modification in ancient convergent margins, limiting evaluation of the tectonic styles and processes that modify continental crust during orogenesis. We present trace element igneous geochemical crustal thickness proxies that quantitatively track the crustal thickness evolution of the long-lived Proterozoic active margin in the southwestern U.S.A. We integrate these results with geobarometric data to constrain the mode of crustal modification. The data indicate a complex record of crustal thickness change in space and time and evolving orogenic styles. Geochemical proxies at 1.84–1.72 Ga are consistent with 20–40 km thick magmatic arcs that were locally thickened to ∼50 km during ∼1.75 Ga tectonism. During the Yavapai orogeny, 1.72–1.69 Ga, a ∼200-km-wide belt of 50-60 km thick crust extended from southern California to northern Colorado and was rapidly thinned and exhumed by ∼1.68 Ga. Crustal thickening and thinning during the Yavapai orogeny largely occurred by shortening and exhumation, respectively, in the upper 25 km of the crust. Subsequent 1.68–1.60 Ga tectonism involved crustal growth, local crustal thickening, and low-P, high-T metamorphism, consistent with extensional accretionary orogenesis. The 1.47–1.37 Ga Picuris orogeny was associated with 50–60 km thick crust across much of the Southwest and involved crustal shortening with ∼10 km of magmatic underplating. Advective heat from the emplacement of ferroan granites in the mid-crust likely contributed to elevated geothermal gradients and rheologically weakened the crust. Our results suggest evolving orogenic styles in the Southwest from 1.75–1.69 Ga short-lived crustal thickening associated with terrane accretion to 1.69–1.60 Ga largely extensional accretionary orogenesis, and regional, long-lived crustal thickening at 1.47–1.37 Ga involving extensive basaltic underplating. Contrasting with some recent hypotheses, our data document a complex middle Proterozoic record for the Southwest that was not orogenically quiescent or tectonically stagnant but involved complex mountain building styles.

Earth and Planetary Science Letters↗
USGS2012

Methods of practice and guidelines for using survey-grade global navigation satellite systems (GNSS) to establish vertical datum in the United States Geological Survey

Geodetic surveys have evolved through the years to the use of survey-grade (centimeter level) global positioning to perpetuate and post-process vertical datum. The U.S. Geological Survey (USGS) uses Global Navigation Satellite Systems (GNSS) technology to monitor natural hazards, ensure geospatial control for climate and land use change, and gather data necessary for investigative studies related to water, the environment, energy, and ecosystems. Vertical datum is fundamental to a variety of these integrated earth sciences. Essentially GNSS surveys provide a three-dimensional position x, y, and z as a function of the North American Datum of 1983 ellipsoid and the most current hybrid geoid model. A GNSS survey may be approached with post-processed positioning for static observations related to a single point or network, or involve real-time corrections to provide positioning "on-the-fly." Field equipment required to facilitate GNSS surveys range from a single receiver, with a power source for static positioning, to an additional receiver or network communicated by radio or cellular for real-time positioning. A real-time approach in its most common form may be described as a roving receiver augmented by a single-base station receiver, known as a single-base real-time (RT) survey. More efficient real-time methods involving a Real-Time Network (RTN) permit the use of only one roving receiver that is augmented to a network of fixed receivers commonly known as Continually Operating Reference Stations (CORS). A post-processed approach in its most common form involves static data collection at a single point. Data are most commonly post-processed through a universally accepted utility maintained by the National Geodetic Survey (NGS), known as the Online Position User Service (OPUS). More complex post-processed methods involve static observations among a network of additional receivers collecting static data at known benchmarks. Both classifications provide users flexibility regarding efficiency and quality of data collection. Quality assurance of survey-grade global positioning is often overlooked or not understood and perceived uncertainties can be misleading. GNSS users can benefit from a blueprint of data collection standards used to ensure consistency among USGS mission areas. A classification of GNSS survey qualities provide the user with the ability to choose from the highest quality survey used to establish objective points with low uncertainties, identified as a Level I, to a GNSS survey for general topographic control without quality assurance, identified as a Level IV. A Level I survey is strictly limited to post-processed methods, whereas Level II, Level III, and Level IV surveys integrate variations of a RT approach. Among these classifications, techniques involving blunder checks and redundancy are important, and planning that involves the assessment of the overall satellite configuration, as well as terrestrial and space weather, are necessary to ensure an efficient and quality campaign. Although quality indicators and uncertainties are identified in post-processed methods using CORS, the accuracy of a GNSS survey is most effectively expressed as a comparison to a local benchmark that has a high degree of confidence. Real-time and post-processed methods should incorporate these "trusted" benchmarks as a check during any campaign. Global positioning surveys are expected to change rapidly in the future. The expansion of continuously operating reference stations, combined with newly available satellite signals, and enhancements to the conterminous geoid, are all sufficient indicators for substantial growth in real-time positioning and quality thereof.

Techniques and Methods↗
USGS2012

Lithology and evolution of the crust-mantle boundary region in the southwestern Basin and Range Province

Mantle and crustal xenoliths from volcanic rocks in the southwestern Basin and Range province and Colorado Plateau Transition Zone reveal histories of episodic magmatism and deformation that have profoundly influenced the crustal structure of this region. Seismic transects in this area show a strongly reflective Moho of generally low relief, which, in the area of modern transects, consists of a thin zone (<2 km thick) of short reflectors. The upper mantle is transparent and has a P n of 7.8–8.0 km/s similar to much of the western United States. A lower crustal zone, 2–13 km thick, has variable internal reflectivity and a relatively low velocity of 6.6–6.8 km/s. Upper mantle peridotite xenoliths show both ductile and brittle deformational features and have structures and compositions affected by magmatic intrusion; intrusions form complex dike systems and extensive zones of grain boundary infiltration in peridotite xenoliths. Whereas melt infiltration preceded and followed ductile deformation, brittle deformation, represented by closely spaced joint systems and faults, followed ductile deformation and is related to the youngest magmatic episodes. These structural characteristics and high uppermost mantle temperature (∼1000°C) may combine to explain the relatively low P n . Alternating layers of ductily deformed and undeformed peridotites, with or without igneous intrusions, may contribute to the reflectivity of the Moho. Lower crustal xenoliths are dominantly igneous-textured pyroxenites and mafic to intermediate gabbros identical to the dikes in peridotite xenoliths. The crustal xenoliths also commonly are jointed, and in addition many show partial melting and have abundant cavities that probably were filled with CO 2 -rich fluids. These rocks are interpreted as products of underplated magmas that were fed through the mantle dike systems and may represent the lowest crustal unit identified in the seismic records. The mafic compositions and high densities of the crustal xenoliths indicate that the low velocity of the lower crust may be caused in part by fracture systems, partial melts, and high temperatures. Garnet granulite xenoliths from a locality with no mantle peridotite xenoliths probably represent crust of the region before late Miocene extension. Felsic granulite xenoliths from two localities have velocities like those of the two lower crustal units identified seismically and could be present in the modern crust as unequilibrated remnants of old crust. The preferred model for the evolution of the lower lithosphere is one in which extension affects the upper mantle as well as the crust and is overlapped in time by multiple magmatic episodes. The earliest magmatic events preceded extension, and later events accompanied and followed extension.

Arizona, California, Nevada↗
USGS2013

Current research issues related to post-wildfire runoff and erosion processes

Research into post-wildfire effects began in the United States more than 70 years ago and only later extended to other parts of the world. Post-wildfire responses are typically transient, episodic, variable in space and time, dependent on thresholds, and involve multiple processes measured by different methods. These characteristics tend to hinder research progress, but the large empirical knowledge base amassed in different regions of the world suggests that it should now be possible to synthesize the data and make a substantial improvement in the understanding of post-wildfire runoff and erosion response. Thus, it is important to identify and prioritize the research issues related to post-wildfire runoff and erosion. Priority research issues are the need to: (1) organize and synthesize similarities and differences in post-wildfire responses between different fire-prone regions of the world in order to determine common patterns and generalities that can explain cause and effect relations; (2) identify and quantify functional relations between metrics of fire effects and soil hydraulic properties that will better represent the dynamic and transient conditions after a wildfire; (3) determine the interaction between burned landscapes and temporally and spatially variable meso-scale precipitation, which is often the primary driver of post-wildfire runoff and erosion responses; (4) determine functional relations between precipitation, basin morphology, runoff connectivity, contributing area, surface roughness, depression storage, and soil characteristics required to predict the timing, magnitudes, and duration of floods and debris flows from ungaged burned basins; and (5) develop standard measurement methods that will ensure the collection of uniform and comparable runoff and erosion data. Resolution of these issues will help to improve conceptual and computer models of post-wildfire runoff and erosion processes.

Earth-Science Reviews↗
USGS2019

Size distributions of Arctic waterbodies reveal consistent relations in their statistical moments in space and time

Arctic lowlands are characterized by large numbers of small waterbodies, which are known to affect surface energy budgets and the global carbon cycle. Statistical analysis of their size distributions has been hindered by the shortage of observations at sufficiently high spatial resolutions. This situation has now changed with the high-resolution (<5 m) circum-Arctic Permafrost Region Pond and Lake (PeRL) database recently becoming available. We have used this database to make the first consistent, high-resolution estimation of Arctic waterbody size distributions, with surface areas ranging from 0.0001 km 2 (100 m 2 ) to 1 km 2 . We found that the size distributions varied greatly across the thirty study regions investigated and that there was no single universal size distribution function (including power-law distribution functions) appropriate across all of the study regions. We did, however, find close relationships between the statistical moments (mean, variance, and skewness) of the waterbody size distributions from different study regions. Specifically, we found that the spatial variance increased linearly with mean waterbody size ( R 2 = 0.97, p < 2.2e-16) and that the skewness decreased approximately hyperbolically. We have demonstrated that these relationships (1) hold across the 30 Arctic study regions covering a variety of (bio)climatic and permafrost zones, (2) hold over time in two of these study regions for which multi-decadal satellite imagery is available, and (3) can be reproduced by simulating rising water levels in a high-resolution digital elevation model. The consistent spatial and temporal relationships between the statistical moments of the waterbody size distributions underscore the dominance of topographic controls in lowland permafrost areas. These results provide motivation for further analyses of the factors involved in waterbody development and spatial distribution and for investigations into the possibility of using statistical moments to predict future hydrologic dynamics in the Arctic.

Alaska↗
USGS2012

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth↗
USGS1977

Voyager imaging experiment

The overall objective of this experiment is exploratory reconnaissance of Jupiter, Saturn, their satellites, and Saturn's rings. Such reconnaissance, at resolutions and phase angles unobtainable from Earth, can be expected to provide much new data relevant to the atmospheric and/or surface properties of these bodies. The experiment also has the following specific objectives: Observe and characterize the global circulation of the atmospheres of Jupiter and Saturn; Determine the horizontal and vertical structure of the visible clouds and establish their relationship to the belted appearance and dynamical properties of the planetary atmospheres; Determine the vertical structure of high, optically-thin, scattering layers on Jupiter and Saturn; Determine the nature of anomalous features such as the Great Red Spot, South Equatorial Belt disturbances, etc.; Characterize the nature of the colored material in the clouds of Jupiter and Saturn, and identify the nature and sources of chromophores on Io and Titan; Perform comparative geologic studies of many satellites at less than 15-km resolution; Map and characterize the geologic structure of several satellites at high resolution (???1 km); Investigate the existence and nature of atmospheres on the satellites; Determine the mass, size, and shape of many of the satellites by direct measurement; Determine the direction of the spin axes and periods of rotation of several satellites, and establish coordinate systems for the larger satellites; Map the radial distribution of material in Saturn's rings at high resolution; Determine the optical scattering properties of the primaries, rings, and satellites at several wavelengths and phase angles; Search for novel physical phenomena, e.g., phenomena associated with the Io flux tube, meteors, aurorae, lightning, or satellite shadows. ?? 1977 D. Reidel Publishing Company.

Space Science Reviews↗
USGS2017

Lithological influences on contemporary and long-term regolith weathering at the Luquillo Critical Zone Observatory

Lithologic differences give rise to the differential weatherability of the Earth’s surface and globally variable silicate weathering fluxes, which provide an important negative feedback on climate over geologic timescales. To isolate the influence of lithology on weathering rates and mechanisms, we compare two nearby catchments in the Luquillo Critical Zone Observatory in Puerto Rico, which have similar climate history, relief and vegetation, but differ in bedrock lithology. Regolith and pore water samples with depth were collected from two ridgetops and at three sites along a slope transect in the volcaniclastic Bisley catchment and compared to existing data from the granitic Río Icacos catchment. The depth variations of solid-state and pore water chemistry and quantitative mineralogy were used to calculate mass transfer (tau) and weathering solute profiles, which in turn were used to determine weathering mechanisms and to estimate weathering rates. Regolith formed on both lithologies is highly leached of most labile elements, although Mg and K are less depleted in the granitic than in the volcaniclastic profiles, reflecting residual biotite in the granitic regolith not present in the volcaniclastics. Profiles of both lithologies that terminate at bedrock corestones are less weathered at depth, near the rock-regolith interfaces. Mg fluxes in the volcaniclastics derive primarily from dissolution of chlorite near the rock-regolith interface and from dissolution of illite and secondary phases in the upper regolith, whereas in the granitic profile, Mg and K fluxes derive from biotite dissolution. Long-term mineral dissolution rates and weathering fluxes were determined by integrating mass losses over the thickness of solid-state weathering fronts, and are therefore averages over the timescale of regolith development. Resulting long-term dissolution rates for minerals in the volcaniclastic regolith include chlorite: 8.9 × 10 −14 mol m −2 s −1 , illite: 2.1 × 10 −14 mol m −2 s −1 and kaolinite: 4.0 × 10 −14 mol m −2 s −1 . Long-term weathering fluxes are several orders of magnitude lower in the granitic regolith than in the volcaniclastic, despite higher abundances of several elements in the granitic regolith. Contemporary weathering fluxes were determined from net (rain-corrected) solute profiles and thus represent rates over the residence time of water in the regolith. Contemporary weathering fluxes within the granitic regolith are similar to the long-term fluxes. In contrast, the long-term fluxes are faster than the contemporary fluxes in the volcaniclastic regolith. Contemporary fluxes in the granitic regolith are generally also slightly faster than in the volcaniclastic. The differences in weathering fluxes over space and time between these two watersheds indicate significant lithologic control of chemical weathering mechanisms and rates.

Bisley watersheds↗
USGS2022

Mapping a magnetic superstorm: March 1989 geoelectric hazards and impacts on United States power systems

A study is made of the relationships between geomagnetic and geoelectric field variation, Earth-surface impedance, and operational interference (anomalies) experienced on electric-power systems across the contiguous United States during the March 13-14, 1989 magnetic storm. For this, a 1-minute-resolution sequence of geomagnetic field maps is constructed from magnetometer time series acquired at ground-based observatories. Induced geoelectric field maps are calculated by convolving the geomagnetic maps with magnetotelluric impedance tensors. During the storm, anomalies were concentrated where the lithosphere is electrically resistive, and when and where geoelectric field amplitudes were high. This was particularly true in the Mid-Atlantic, Northeast, and the upper Midwest. Few anomalies were experienced in other parts of the Midwest and across much of the West, where the lithosphere is more conductive, and when and where geoelectric field amplitudes were low. Peak 1-minute-resolution geoelectric field amplitude ranged from 21.66 V/km in Maine and 19.02 V/km in Virginia to < 0.02 V/km in Idaho. Latitude-dependent organization of geoelectric hazards by auroral-zone electrojet currents is detectable, but it is much weaker than geographic organization due to surface impedance. Hazardous geoelectric fields were induced during different storm phases, at different local times, and, by inference, by a variety of ionospheric currents. Compared to geoelectric field amplitudes realized across the United States during March 1989, hazard maps used by utility companies to estimate systems exposure have much less geographic detail and a much smaller maximum-to-minimum range in geoelectric field amplitude. Future research will benet from denser geomagnetic monitoring, additional magnetotelluric surveying, and access to power-system impact data.

Space Weather↗
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↗