SEARCH · Geology Reports

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 721 records · Page 40Linked to original sources

USGS2016

Landsat 8: The plans, the reality, and the legacy

Landsat 8, originally known as the Landsat Data Continuity Mission (LDCM), is a National Aeronautics and Space Administration (NASA)-U.S. Geological Survey (USGS) partnership that continues the legacy of continuous moderate resolution observations started in 1972. The conception of LDCM to the reality of Landsat 8 followed an arduous path extending over nearly 13 years, but the successful launch on February 11, 2013 ensures the continuity of the unparalleled Landsat record. The USGS took over mission operations on May 30, 2013 and renamed LCDM to Landsat 8. Access to Landsat 8 data was opened to users worldwide. Three years following launch we evaluate the science and applications impact of Landsat 8. With a mission objective to enable the detection and characterization of global land changes at a scale where differentiation between natural and human-induced causes of change is possible, LDCM promised incremental technical improvements in capabilities needed for Landsat scientific and applications investigations. Results show that with Landsat 8, we are acquiring more data than ever before, the radiometric and geometric quality of data are generally technically superior to data acquired by past Landsat missions, and the new measurements, e.g., the coastal aerosol and cirrus bands, are opening new opportunities. Collectively, these improvements are sparking the growth of science and applications opportunities. Equally important, with Landsat 7 still operational, we have returned to global imaging on an 8-day cycle, a capability that ended when Landsat 5 ceased operational Earth imaging in November 2011. As a result, the Landsat program is on secure footings and planning is underway to extend the record for another 20 or more years.

Remote Sensing of Environment↗

USGS2009

Gas hydrate saturations estimated from fractured reservoir at Site NGHP-01-10, Krishna-Godavari Basin, India

During the Indian National Gas Hydrate Program Expedition 01 (NGHP-Ol), one of the richest marine gas hydrate accumulations was discovered at Site NGHP-01-10 in the Krishna-Godavari Basin. The occurrence of concentrated gas hydrate at this site is primarily controlled by the presence of fractures. Assuming the resistivity of gas hydratebearing sediments is isotropic, th?? conventional Archie analysis using the logging while drilling resistivity log yields gas hydrate saturations greater than 50% (as high as ???80%) of the pore space for the depth interval between ???25 and ???160 m below seafloor. On the other hand, gas hydrate saturations estimated from pressure cores from nearby wells were less than ???26% of the pore space. Although intrasite variability may contribute to the difference, the primary cause of the saturation difference is attributed to the anisotropic nature of the reservoir due to gas hydrate in high-angle fractures. Archie's law can be used to estimate gas hydrate saturations in anisotropic reservoir, with additional information such as elastic velocities to constrain Archie cementation parameters m and the saturation exponent n. Theory indicates that m and n depend on the direction of the measurement relative to fracture orientation, as well as depending on gas hydrate saturation. By using higher values of m and n in the resistivity analysis for fractured reservoirs, the difference between saturation estimates is significantly reduced, although a sizable difference remains. To better understand the nature of fractured reservoirs, wireline P and S wave velocities were also incorporated into the analysis.

Journal of Geophysical Research B: Solid Earth↗

USGS1978

Megalineament in southeastern Alaska marks southwest edge of Coast Range batholithic complex

The Coast Range megalineament is a prominent, nearly continuous topographic and structural feature that extends southeastward about 550 km (330 mi) from its junction with the Chatham Strait – Lynn Canal fault at Point Sherman to Tongass Passage near the mouth of Pearse Canal where it leaves southeastern Alaska. It probably extends still further southeastward into British Columbia along Work Channel and Chatham Sound – Grenville Channel.The megalineament is a zone a few hundred metres to 10 km (6 mi) wide in which closely spaced joints, foliation, compositional layering, and small faults define the megalineament trend. The zone usually coincides with topographic depressions apparently caused by selective fluvial and glacial erosion of the less resistant rocks of the zone.Studies in the Juneau, Endicott Arm, and Behm Canal areas indicate that the megalineament (1) is locally the site of lateral and (or) vertical separations of no greater than several kilometres; (2) does not mark a major structural or metamorphic discontinuity in the near-surface rocks; (3) may be located near a pre-metamorphic and pre-intrusive discontinuity; (4) is consistently associated with and parallel to steep gradients in both the gravity and aeromagnetic fields; and (5) probably is the surface expression of the western contact, at depth, of the intrusive rocks and gneisses of the Coast Range batholithic complex with the schists to the southwest.

Alaska↗

USGS2008

Probabilistic seismic hazard in the San Francisco Bay area based on a simplified viscoelastic cycle model of fault interactions

We construct a viscoelastic cycle model of plate boundary deformation that includes the effect of time-dependent interseismic strain accumulation, coseismic strain release, and viscoelastic relaxation of the substrate beneath the seismogenic crust. For a given fault system, time-averaged stress changes at any point (not on a fault) are constrained to zero; that is, kinematic consistency is enforced for the fault system. The dates of last rupture, mean recurrence times, and the slip distributions of the (assumed) repeating ruptures are key inputs into the viscoelastic cycle model. This simple formulation allows construction of stress evolution at all points in the plate boundary zone for purposes of probabilistic seismic hazard analysis (PSHA). Stress evolution is combined with a Coulomb failure stress threshold at representative points on the fault segments to estimate the times of their respective future ruptures. In our PSHA we consider uncertainties in a four-dimensional parameter space: the rupture peridocities, slip distributions, time of last earthquake (for prehistoric ruptures) and Coulomb failure stress thresholds. We apply this methodology to the San Francisco Bay region using a recently determined fault chronology of area faults. Assuming single-segment rupture scenarios, we find that fature rupture probabilities of area faults in the coming decades are the highest for the southern Hayward, Rodgers Creek, and northern Calaveras faults. This conclusion is qualitatively similar to that of Working Group on California Earthquake Probabilities, but the probabilities derived here are significantly higher. Given that fault rupture probabilities are highly model-dependent, no single model should be used to assess to time-dependent rupture probabilities. We suggest that several models, including the present one, be used in a comprehensive PSHA methodology, as was done by Working Group on California Earthquake Probabilities.

Journal of Geophysical Research B: Solid Earth↗

USGS2020

Plastic faulting in ice

Plastic faulting is a brittle‐like failure phenomenon exhibited by water ice and several other rock types under confinement. It is suspected to be the mechanism of deep earthquakes and extreme cases of shear localization in shallow rocks. Unlike ordinary Coulombic failure, plastic faulting is characterized by a pressure‐independent failure strength and fault plane oriented 45° to maximum principal stress. To research the question of how the instability initiates, we conducted over 50 constant‐displacement‐rate experiments on polycrystalline ice (phases Ih and II) near the brittle‐to‐ductile (B‐D) transition, at confining pressures P = 0–300 MPa, applied strain rates = 5 × 10 −5 – 7 × 10 −3 s −1 , temperatures T = 105–233 K, and mean grain sizes d = 0.25–1.18 mm. We find that (1) the width of the B‐D transition in variable space is vanishingly narrow, to the point of appearing as a crossover, (2) a plastic fault plane, once formed, is not a zone of subsequent weakness, (3) distributed ice I→II phase transformation in small amounts (<1 vol%) shows no causal relationship to subsequent failure, and (4) plastic faulting also occurs in ice II. We hypothesize that the elusive nucleating “trigger” parallels that of metals and ceramics undergoing severe plastic deformation, wherein transient local structural rearrangement occurs, in turn causing material strength to drop to a level sufficiently low, in a volume sufficiently large, that adiabatic instability is nucleated. Our results do not require and often are inconsistent with phase transformation. Plastic faulting may therefore be available to all solids undergoing severe deformation, and its appearance in so few is simply the result of insufficiently extreme conditions.

Journal of Geophysical Research- Solid Earth↗

USGS1976

Tectonic implications of space-time patterns of Cenozoic magmatism in the western United States

Locations of 2,100 radiometrically dated igneous rocks were plotted on a series of 20 maps, each representing an interval within the period 80 m.y. B.P. to present. Derivative maps showing the distributions in space and time of dated granitic intrusive rocks, silicic lavas and domes, ash-flow tuffs, andesitic-dacitic rocks, and basalts depict well the two main petrogenetic assemblages noted previously by others: (1) mainly intermediate andesitic-dacitic suites, including associated granitic intrusive rocks, silicic extrusive rocks, and minor basaltic lavas, are interpreted as reflecting plate interactions related to subduction along the continental margin; and (2) bimodal suites, dominantly basaltic but with minor silicic extrusive rocks, are interpreted as reflecting extensional tectonics. Space-time distribution of the two assemblages suggests that magmatic arcs extended continously parallel to the continental margin from Canada to Mexico in latest Mesozoic and in Oligocene times. An early Cenozoic null in magmatism in the Great Basin may delineate the region where subduction was arrested temporarily by development of the proto-San Andreas fault as a transform in coastal California or, alternatively, may reflect complex subsurface configurations of subducted plates. The late Cenozoic transition from subduction-related magmatism to extention-related basaltic volcanism in the southern Cordillera occurred at different times in different areas in harmony with current concepts about the migration of the Mendocino triple junction as the modern San Andreas transform fault was formed. The plots also reveal the existence of several discrete magmatic loci where igneous activity of various kinds was characteristically more intense and long-lived than elsewhere.

Earth and Planetary Science Letters↗

USGS2022

Complex magmatic-tectonic interactions during the 2020 Makushin Volcano, Alaska, earthquake swarm

On June 15, 2020, at 21:16 UTC, a locally-felt earthquake of magnitude 4.2 struck Unalaska Island, Alaska, ∼15 km west of the town of Unalaska and the large fishing port of Dutch Harbor. The event was followed by a M4.1 earthquake at 00:34 UTC and several M3+ aftershocks, initiating a prolific sequence with hundreds of earthquakes recorded into late December. The earthquakes all locate about 12 km southeast of the summit of Makushin Volcano at 7 to 10 km depth. To date, no eruptive activity or other surface changes have been observed at the volcano in webcam images, GPS or InSAR. Seismic bursts close to volcanoes are often associated with the onset of unrest that can lead to eruption. However, determining whether seismicity reflects magmatic rather than tectonic stresses is often challenging, although critical for hazard assessments and risk management strategies. To investigate the triggering mechanisms of the recent Makushin seismicity, we integrate information from space-time patterns of the earthquake hypocenters with their fault-plane solutions. We relocate the swarm events using double-difference relocation techniques and a 3D velocity model and find that the earthquakes, although they seem to follow two predominant orientations (NW-SE and SW-NE), do not show clear clustering into preferred alignments. Similarly, we do not observe pronounced migration in time and space. Fault-plane solutions (FPS) for all but one M2.5+ earthquakes have P-axis orientations consistent with subhorizontal NW-SE oriented regional maximum compression, whereas many of the lower-magnitude earthquakes have P-axes perpendicular to regional maximum compression. This provides evidence for the presence of a local stress field likely induced by magma intrusion. Results from Coulomb stress modeling are also consistent with dike inflation modulated by stresses induced by the M4+ earthquakes. The seismic swarm is thus likely linked to a superposition of driving stresses from both magmatic and tectonic processes on pre-existing faults. The case of the 2020 Makushin swarm, with its unusual characteristics, challenges traditional swarm classification schemes and suggests that a reconsideration of the definition of seismic swarms as having the maximum magnitude event in the middle of the swarm is warranted.

Alaska↗

USGS2011

InSAR observations of aseismic slip associated with an earthquake swarm in the Columbia River flood basalts

In 2009 a swarm of small shallow earthquakes occurred within the basalt flows of the Columbia River Basalt Group (CRBG). The swarm occurred within a dense seismic network in the U.S. Department of Energys Hanford Site. Data from the seismic network along with interferometric synthetic aperture radar (InSAR) data from the European Space Agencys (ESA) ENVISAT satellite provide insight into the nature of the swarm. By modeling the InSAR deformation data we constructed a model that consists of a shallow thrust fault and a near horizontal fault. We suggest that the near horizontal lying fault is a bedding-plane fault located between basalt flows. The geodetic moment of the modeled fault system is about eight times the cumulative seismic moment of the swarm. Precise location estimates of the swarm earthquakes indicate that the area of highest slip on the thrust fault, ???70mm of slip less than ???0.5km depth, was not located within the swarm cluster. Most of the slip on the faults appears to have progressed aseismically and we suggest that interbed sediments play a central role in the slip process. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Deformation from 1973 to 1987 in the epicentral area of the 1959 Hebgen Lake, Montana, earthquake (Ms = 7.5)

A 40‐km aperture trilateration network centered on the 1959 Hebgen Lake earthquake epicenter has been surveyed in 1973, 1974, 1976, 1978, 1981, 1984, and 1987. The deformation inferred from those surveys is described roughly by a uniaxial, 0.266 ± 0.014 μstrain/yr, N15°E ± 1°extension that is uniform in both time and space. That extension is orthogonal to the strike (N78°W ± 5°) of the 1959 rupture plane inferred from focal mechanism solutions. The absence of strain accumulation in the N75°W direction suggests that the source of deformation must be greatly elongated in that direction. Looked at in greater detail, the deformation is found to be concentrated in the northern part of the network. The zone of concentrated deformation coincides with a part of a well‐defined trend in seismicity that extends at least 100 km N75°W from the Sour Creek resurgent dome in the Yellowstone caldera. Extension perpendicular to that trend is indicated by focal mechanism solutions. Thus, the zone of rifting identified within the Hebgen Lake network apparently extends east southeast to the Yellowstone caldera. About 8 mm/yr extension is observed across a 20‐km width within the zone.

Montana, Wyoming↗

USGS2012

Petrology and geochronology of metamorphosed volcanic rocks and a middle Cretaceous volcanic neck in the east-central Sierra Nevada, California

Metamorphosed Mesozoic volcanic rocks from the east-central Sierra Nevada range in composition from basalt to rhyolite and have ages, based on whole rock Rb-Sr and U-Pb zircon dating, of about 237–224, 185, 163, 134, and 100 m.y. The major plutons of the batholith in this area are of Triassic (215–200 m.y.) and Cretaceous (94–80 m.y.) ages. Initial 87 Sr/ 86 Sr values for the metamorphosed volcanic rocks of the area are in the range from 0.7042 to 0.7058 and are generally different from the values for the surrounding batholithic rocks (0.7056–0.7066). A circular, zoned granitic pluton, with an outcrop area of 2.5 km 2 , similar in appearance to a ring dike complex, was apparently a conduit for some or possibly all of the middle-Cretaceous metamorphosed volcanic rocks exposed about 5 km to the south in the western part of the Ritter Range. Samples from the metamorphosed volcanic rocks and the pluton yield a Rb-Sr whole rock isochron age of 99.9 ± 2.2 m.y. with an initial 87 Sr/ 86 Sr of 0.7048 ± 0.0001. Major element variation diagrams of the pluton and volcanic rocks define coincident compositional trends. The pluton, intruded into Early Jurassic metamorphosed volcanic rocks (Rb-Sr whole rock isochron age of 185 ± 6 m.y., initial 87 / 86 Sr = 0.70474 ± 0.00009), is composed of a central body of granite that is intruded into and almost completely surrounded by a crescent-shaped outer rim of quartz monzodiorite. Contact relations between the granitic rocks, intrusive in some places and completely gradational at other places, indicate that the quartz monzodiorite was not completely crystallized prior to the intrusion of the granite. Modal minerology within the pluton is also consistent with this inferred crystallization history. The granitic rock textures are cataclastic, characterized by a groundmass of feldspar and quartz that surrounds larger grains of feldspar, hornblende, and biotite. Widely spaced aplite dikes from the granite intrude the quartz monzodiorite along the contact between the granitic rocks. The aplite dikes and interstitial alkali feldspar associated with myrmekite in the pluton, being undeformed, indicate crystallization of a small amount of silicate liquid after cataclasis. Aplite compositions suggest a minimum crystallization pressure of 0.5 kb (1–2 km); this depth is consistent with crystallization near the base of a large stratovolcano. The pluton is characterized by steep, linear features, defined by mineral clots and inclusions, that probably formed during the upward movement of the nearly crystallized magma in the volcanic conduit. The ages of volcanic events relative to the ages of the major intrusive epochs and the major element and isotopic compositions of the volcanic rocks relative to the major plutons indicate that the volcanic rocks are not simply or directly related to the major plutons in the Sierra Nevada.

Journal of Geophysical Research Solid Earth↗

USGS2009

How the continents deform: The evidence from tectonic geodesy

Space geodesy now provides quantitative maps of the surface velocity field within tectonically active regions, supplying constraints on the spatial distribution of deformation, the forces that drive it, and the brittle and ductile properties of continental lithosphere. Deformation is usefully described as relative motions among elastic blocks and is block-like because major faults are weaker than adjacent intact crust. Despite similarities, continental block kinematics differs from global plate tectonics: blocks are much smaller, typically ∼100–1000 km in size; departures from block rigidity are sometimes measurable; and blocks evolve over ∼1–10 Ma timescales, particularly near their often geometrically irregular boundaries. Quantitatively relating deformation to the forces that drive it requires simplifying assumptions about the strength distribution in the lithosphere. If brittle/elastic crust is strongest, interactions among blocks control the deformation. If ductile lithosphere is the stronger, its flow properties determine the surface deformation, and a continuum approach is preferable.

Annual Review of Earth and Planetary Sciences↗

USGS2009

Kilauea slow slip events: Identification, source inversions, and relation to seismicity

Several slow slip events beneath the south flank of Kilauea Volcano, Hawaii, have been inferred from transient displacements in daily GPS positions. To search for smaller events that may be close to the noise level in the GPS time series, we compare displacement fields on Kilauea's south flank with displacement patterns in previously identified slow slip events. Matching displacement patterns are found for several new candidate events, although displacements are much smaller than previously identified events. One of the candidates, 29 May 2000, is coincident with a microearthquake swarm, as are all of the previously identified slow slip events. The microearthquakes follow the onset of slow slip, implying that they are triggered by stress changes during slip. The new slow slip event brings the total number of events on Kilauea, between 1997 and 2007, to eight, the smallest having M W = 5.3, and the largest having M W = 6.0. While the recurrence time between the four largest events is 2.11 ± 0.01 years, the repeat time for all eight events is 0.9 ± 0.6 years. We invert for the fault geometry and distribution of slip during the slow slip events. The optimal source depths of 5 km, assuming uniform slip dislocations in an elastic half‐space, are considerably shallower than the accompanying swarm earthquakes (6.5–8.5 km), which would place the earthquakes in a zone of decreased Coulomb stress. Inversions including the effects of topography and layered elastic structure in the forward models favor depths comparable to microearthquake depths, such that the earthquakes are located in a region of increased Coulomb stress. We also invert for time‐dependent fault slip directly from the 30 s GPS phase observations, constraining the source to the optimal uniform slip geometry. On the basis of these inversions, the larger events last between 1.5–2.2 days. The data are unable to resolve migration of slip along the fault. The temporal pattern of accompanying microearthquakes is consistent with the fault slip history assuming a seismicity rate theory based on rate and state‐friction, making the swarm earthquakes coshocks and aftershocks of the slow slip events.

Hawaii↗

USGS2023

Scaling microseismic cloud shape during hydraulic stimulation using in-situ stress and permeability

Forecasting microseismic cloud shape as a proxy of stimulated rock volume may improve the design of an energy extraction system. The microseismic cloud created during hydraulic stimulation of geothermal reservoirs is known empirically to extend in the general direction of the maximum principal stress. However, this empirical relationship is often inconsistent with reported results, and the cloud growth process remains poorly understood. This study investigates microseismic cloud growth using data obtained from a hydraulic stimulation project in Basel, Switzerland, and explores its correlation with measured in situ stress. We applied principal component analysis to a time series of microseismicity for macroscopic characterization of microseismic cloud growth in two- and three-dimensional space. The microseismic cloud, in addition to extending in the general direction of maximum principal stress, expanded in the direction of intermediate principal stress. The orientation of the least microseismic cloud growth was stable and almost identical to the minimum principal stress direction. Further, microseismic cloud shape ratios showed good agreement when compared with in situ stress magnitude ratios. The permeability tensor estimated from microseismicity also provided a good correlation in terms of direction and magnitude with the microseismic cloud growth. We show that in situ stress plays a dominant role by controlling the permeability of each existing fracture in the reservoir fracture system. Consequently, microseismic cloud growth can be scaled by in situ stress as a first-order approximation if there is sufficient variation in the orientation of existing faults.

Basel↗

USGS2011

ASTER Global Digital Elevation Model Version 2 - summary of validation results

On June 29, 2009, NASA and the Ministry of Economy, Trade and Industry (METI) of Japan released a Global Digital Elevation Model (GDEM) to users worldwide at no charge as a contribution to the Global Earth Observing System of Systems (GEOSS). This “version 1” ASTER GDEM (GDEM1) was compiled from over 1.2 million scenebased DEMs covering land surfaces between 83°N and 83°S latitudes. A joint U.S.-Japan validation team assessed the accuracy of the GDEM1, augmented by a team of 20 cooperators. The GDEM1 was found to have an overall accuracy of around 20 meters at the 95% confidence level. The team also noted several artifacts associated with poor stereo coverage at high latitudes, cloud contamination, water masking issues and the stacking process used to produce the GDEM1 from individual scene-based DEMs (ASTER GDEM Validation Team, 2009). Two independent horizontal resolution studies estimated the effective spatial resolution of the GDEM1 to be on the order of 120 meters. A second version of the ASTER GDEM (GDEM2) is scheduled for release by NASA and METI in mid-October, 2011. Improvements in the GDEM2 result from acquiring 260,000 additional scenes to improve coverage, a smaller correlation kernel to yield higher spatial resolution, and improved water masking. As with the GDEM1, the GDEM2 validation was performed by the U.S. and Japanese partners. Vertical accuracy assessments included a comparison of the GDEM2 against absolute geodetic references over the Conterminous US (CONUS), against national elevation grids over the US and Japan, against the Shuttle Radar Topography Mission (SRTM) 1 arc-second elevation grids over the US and 20 sites around the globe, and against space borne laser altimeter data globally. Horizontal accuracy assessments were conducted as part of the Japan and the global SRTM studies, and horizontal resolution studies were conducted in both Japan and the US. Each group documented changes in artifacts in GDEM2 due to processing improvements. The absolute vertical accuracy study found the GDEM2 to be within -0.20 meters on average when compared against 18,000 geodetic control points over the CONUS, with an accuracy of 17 meters at the 95% confidence level. The Japan study noted the GDEM2 differed from the 10-meter national elevation grid by -0.7 meters over bare areas, and by 7.4 meters over forested areas. Similarly, the CONUS study noted the GDEM2 to be about 8 meters above the 1 arc-second NED over most forested areas, and more than a meter below NED over bare areas. The global altimeter study found the GDEM2 to be on average within 3 meters of altimeter-derived control, and also documented sensitivity to tree canopy height. The Japan study noted that the horizontal displacement in GDEM1 of 0.95 pixels was reduced to 0.23 pixels in GDEM2. Both teams noted improvements in horizontal resolution, between 71 and 82 meters, comparable to the SRTM 1 arc second elevation model, but at the cost of some increased noise. The number of voids and artifacts noted in GDEM1 were substantially reduced in GDEM2, and in some areas virtually eliminated. Based on these findings, the GDEM validation team recommends the release of the GDEM2 to the public, acknowledging that, while vastly improved, some artifacts still exist which could affect its utility in certain applications.

Report↗

USGS2021

The birth of a Hawaiian fissure eruption

Most basaltic explosive eruptions intensify abruptly, allowing little time to document processes at the start of eruption. One opportunity came with the initiation of activity from fissure 8 (F8) during the 2018 eruption on the lower East Rift Zone of Kīlauea, Hawaii. F8 erupted in four episodes. We recorded 28 min of high‐definition video during a 51‐min period, capturing the onset of the second episode on 5 May. From the videos, we were able to analyze the following in‐flight parameters: frequency and duration of explosions; ejecta heights; pyroclast exit velocities; in‐flight total mass and estimated mass eruption rates; and the in‐flight total grain size distributions. The videos record a transition from initial pulsating outgassing, via spaced, but increasingly rapid, discrete explosions, to quasisustained, unsteady fountaining. This transition accompanied waxing intensity (mass flux) of the F8 eruption. We infer that all activity was driven by a combination of the ascent of a coupled mixture of small bubbles and melt, and the buoyant rise of decoupled gas slugs and/or pockets. The balance between these two types of concurrent flow determined the exact form of the eruptive activity at any point in time, and changes to their relative contributions drove the transition we observed at early F8. Qualitative observations of other Hawaiian fountains at Kīlauea suggest that this physical model may apply more generally. This study demonstrates the value of in‐flight parameters derived from high‐resolution videos, which offer a rapid and highly time‐sensitive alternative to measurements based on sampling of deposits posteruption.

Island of Hawai'i↗

USGS2004

User's Guide for the MapImage Reprojection Software Package, Version 1.01

Scientists routinely accomplish small-scale geospatial modeling in the raster domain, using high-resolution datasets (such as 30-m data) for large parts of continents and low-resolution to high-resolution datasets for the entire globe. Recently, Usery and others (2003a) expanded on the previously limited empirical work with real geographic data by compiling and tabulating the accuracy of categorical areas in projected raster datasets of global extent. Geographers and applications programmers at the U.S. Geological Survey's (USGS) Mid-Continent Mapping Center (MCMC) undertook an effort to expand and evolve an internal USGS software package, MapImage, or mapimg, for raster map projection transformation (Usery and others, 2003a). Daniel R. Steinwand of Science Applications International Corporation, Earth Resources Observation Systems Data Center in Sioux Falls, S. Dak., originally developed mapimg for the USGS, basing it on the USGS's General Cartographic Transformation Package (GCTP). It operated as a command line program on the Unix operating system. Through efforts at MCMC, and in coordination with Mr. Steinwand, this program has been transformed from an application based on a command line into a software package based on a graphic user interface for Windows, Linux, and Unix machines. Usery and others (2003b) pointed out that many commercial software packages do not use exact projection equations and that even when exact projection equations are used, the software often results in error and sometimes does not complete the transformation for specific projections, at specific resampling resolutions, and for specific singularities. Direct implementation of point-to-point transformation with appropriate functions yields the variety of projections available in these software packages, but implementation with data other than points requires specific adaptation of the equations or prior preparation of the data to allow the transformation to succeed. Additional constraints apply to global raster data. It appears that some packages use the USGS's GCTP or similar point transformations without adaptation to the specific characteristics of raster data (Usery and others, 2003b). It is most common for programs to compute transformations of raster data in an inverse fashion. Such mapping can result in an erroneous position and replicate data or create pixels not in the original space. As Usery and others (2003a) indicated, mapimg performs a corresponding forward transformation to ensure the same location results from both methods. The primary benefit of this function is to mask cells outside the domain. MapImage 1.01 is now on the Web. You can download the User's Guide, source, and binaries from the following site: http://mcmcweb.er.usgs.gov/carto_research/projection/acc_proj_data.html

Open-File Report↗

USGS2026

Cajon Pass and the southern San Andreas Fault System: Earthquake cycle stress accumulation and present-day loading

With over a century since the last major rupture affecting the wider Los Angeles region, tectonic stress has steadily built along the southern San Andreas and San Jacinto fault systems, raising concerns of an imminent large earthquake. Cajon Pass, located at the junction of these faults, represents a critical site for potential through-going ruptures in Southern California. We constructed new 4D earthquake cycle simulations using a 1000-year paleoseismic rupture history of the San Andreas Fault System (SAFS) to assess spatial and temporal variations in stress. A semi-analytic Fourier transform model was used to compute stress from 3D dislocations in an elastic plate overlying a Maxwell viscoelastic half-space, assuming a complete coseismic reset of resolved shear stress on ruptured elements. Results show highest stress accumulation north of Cajon Pass (∼1.8 MPa/100 years) due to greater slip rates, and lower rates south of Cajon Pass (∼1.0–1.5 MPa/100 years). By 2025, Coulomb stress is estimated at 2.8 MPa on the Mojave South (MOS) segment, 1.8 MPa on the North San Bernardino (NSB1) segment and 3.6 MPa on the San Jacinto Bernardino (SJB) segment. Segments accumulate stresses with characteristic ranges of pre-event stress interpreted as failure thresholds: 1.2–2.7 MPa for MOS, 0.4–1.6 MPa for NSB1, and 1.2–2.9 MPa for SJB. When the stress disparity between segments SJB and MOS narrows, the faults appear to rupture jointly, suggesting that stress levels may control how Cajon Pass acts as an earthquake gate. These results may inform seismic hazard assessments by linking stress evolution to fault interactions.

California↗

USGS1997

Three-dimensional upper crustal velocity structure beneath San Francisco Peninsula, California

This paper presents new seismic data from, and crustal models of the San Francisco Peninsula. In much of central California the San Andreas fault juxtaposes the Cretaceous granitic Salinian terrane on its west and the Late Mesozoic/Early Tertiary Franciscan Complex on its east. On San Francisco Peninsula, however, the present-day San Andreas fault is completely within a Franciscan terrane, and the Pilarcitos fault, located southwest of the San Andreas, marks the Salinian-Franciscan boundary. This circumstance has evoked two different explanations: either the Pilarcitos is a thrust fault that has pushed Franciscan rocks over Salinian rocks or the Pilarcitos is a transform fault that has accommodated significant right-lateral slip. In an effort to better resolve the subsurface structure of the peninsula faults, we established a temporary network of 31 seismographs arrayed across the San Andreas fault and the subparallel Pilarcitos fault at ∼1–2 km spacings. These instruments were deployed during the first 6 months of 1995 and recorded local earthquakes, air gun sources set off in San Francisco Bay, and explosive sources. Travel times from these sources were used to augment earthquake arrival times recorded by the Northern California Seismic Network and were inverted for three-dimensional velocity structure. Results show lateral velocity changes at depth (∼0.5–7 km) that correlate with downward vertical projections of the surface traces of the San Andreas and Pilarcitos faults. We thus interpret the faults as high-angle to vertical features (constrained to a 70°–110° dip range). From this we conclude that the Pilarcitos fault is probably an important strike-slip fault that accommodated much of the right-lateral plate boundary strain on the peninsula prior to the initiation of the modern-day San Andreas fault in this region sometime after about 3.0 m.y. ago.

Journal of Geophysical Research B: Solid Earth↗