Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research”

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

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

At least 1,459 records · Page 81Linked to original sources
USGS2012

Deep volcanic tremor and magma ascent mechanism under Kilauea, Hawaii

Deep harmonic tremor originating at depths around 40 km under Kilauea was studied using records accumulated since 1962 at the Hawaii Volcano Observatory of the U.S. Geological Survey. The deep source of the tremor was determined by onset times and confirmed by the relative amplitude across the island-wide network of seismometers. The period of tremor was conclusively shown to be determined by the source effect and not by the path or station site effect because the period would change considerably in time but maintained uniformity across the seismic net during the tremor episode. The tremor appeared to be primarily composed of P waves. We interpret the observed period and amplitude in terms of the stationary crack model of Aki et al. (1977) and find that the seismic moment rates for deep tremors are considerably larger than those for shallow-tremors suggesting more vigorous transport for the former. We propose a kinematic source model which may be more appropriate for deep tremor. According to this model, a measurable quantity called ‘reduced displacement’ is directly proportional to the rate of magma flow. A systematic search for deep tremor episodes was made for the period from 1962 through 1979, and the amplitude, period, and duration of the tremor were tabulated. We then constructed a cumulative reduced-displacement plot over the 18-year period. The result shows a generally steady process which does not seem to be significantly affected by major eruptions and large earthquakes near the surface. The total magma flow estimated from the reduced displacement is however, one order of magnitude smaller than that estimated by Swanson (1972). It may be that most channels transport magma aseismically, and only those with strong barriers generate tremor.

Hawaii↗
USGS2012

Magnetic signals from the core of the earth and secular variation

An oscillating, radial magnetic dipole source was assumed to exist in the core of the earth, 100 km beneath the core-mantle boundary. As an approximation, electromagnetic propagation was assumed in the core in lieu of hydromagnetic propagation, which could not be used because of unknown internal fields.Using Debye potentials, the radial and horizontal components of the surface fields were calculated using various assumed conductivity parameters in the core and in the mantle. It is concluded that most spherical harmonic models of the earth's magnetic field do not include enough terms to describe properly the field of core sources with periods of tens of years. These short-period variations are especially important in describing the secular variation. Because of this, a proper description of the secular variation requires more spherical harmonic terms than are required for the field itself. Inadequate representation of short-period variations in spherical harmonic models may contribute to the rapid deterioration of predictive models. Alternatives to spherical harmonic analysis for secular variation should be investigated; regardless of the method used, a much greater spatial distribution of high quality secular variation data is needed.

Journal of Geophysical Research Solid Earth↗
USGS2012

Quaternary and geodetically measured crustal movements in the Tokai District, Central Honshu, Japan

Geodetic measurements made in this century demonstrate that modern movements correlate well with deformation that has occurred during the past 2 m.y. Both show effects of subduction and of the Quaternary collision of Izu Peninsula with Central Honshu. Westward from Suruga Bay, deformation is greatest near the coast, decreases markedly towards the uplifting Akaishi Mountains, and then increases again in the conjugate block-faulted terranes of Chubu district. Quaternary faults, historic earthquakes, and the distribution of horizontal shear straining since ∼1880 all reflect this deformation pattern. Since 1900 the west coast of Suruga Bay has undergone secular subsidence, horizontal contraction and tilting towards Suruga Trough, all typical of the strain buildup process that occurs landward of subduction zones. However, futher westward, active mountain building is occuring, and the Akaishi Ranges are rising 3 mm/a. Unlike the episodic aseismic uplift of the Transverse Ranges of southern California, movements here have been relatively uniform. It seems likely that an increase in the tilt rate of the Suruga coast recognized since ∼1973 is part of a long-term readjustment to the great 1944 Tonankai earthquake. Data show that a monotonic increase has occurred since at least 1949, and similar effects are seen adjacent to the 1944 rupture zone over 200 km southwest along the Nankai Trough.

Journal of Geophysical Research Solid Earth↗
USGS2012

Gravity model and structural implications of the Goddard Pendant, Sierra Nevada, California

A subsurface model for the Goddard pendant is constructed from a residual gravity high of about 7 mGal over the pendant. The model, which is the simplest and most geologically reasonable possibility, shows a metamorphic block that tapers with depth and extends about 3.5 km below the surface. The structures in the Goddard pendant are similar in style and orientation to those in other Sierra Nevada pendants, indicating that the country rock was neither deformed nor rotated during pluton emplacement. Consequently, emplacement must have been a passive rather than a forceful process. The pendant itself represents a piece of country rock trapped between plutons which are dome shaped in cross section.

Journal of Geophysical Research Solid Earth↗
USGS2012

Minor and trace element geochemistry of volcanic rocks dredged from the Galapagos spreading center: Role of crystal fractionation and mantle heterogeneity

A wide range of rock types (abyssal tholeiite, Fe-Ti-rich basalt, andesite, and rhyodacite) were dredged from near 95°W and 85°W on the Galapagos spreading center. Computer modeling of major element compositions has shown that these rocks could be derived from common parental magmas by successive degrees of fractional crystallization. However, the P 2 O 5 /K 2 O ratio averages 0.83 at 95°W and 1.66 at 85°W and implies distinct mantle source compositions for the two areas. These source regions also have different rare earth element (REE) abundance patterns, with [La/Sm] EF = 0.67 at 95°W and 0.46 at 85°W. The sequence of fractionated lavas differs for the two areas and indicates earlier fractionation of apatite and titanomagnetite in the lavas from 95°W. The mantle source regions for these two areas are interpreted to be depleted in incompatible (and volatile?) elements, although the source region beneath 95°W is less severely depleted in La and K. Incompatible trace element abundances in 26 samples are used to infer that the range of Fe-Ti-rich basalt from 85°W represents 19 to 35% residual liquid following crystal fractionation of a mineral assemblage of plagioclase, clinopyroxene, and lesser olivine. The most highly differentiated samples have also had less than 1% titanomagnetite removed. Most samples from 85°W can be related to a common parental magma that contained approximately 9 wt % FeO*, 1 wt % TiO 2 , and had an Mg number (Mg# = 100 Mg/(Mg + Fe 2+ )) of about 65. Although the samples from 95°W cannot all be derived from a common parental magma, the inferred parental magmas may have been derived by varying degrees of partial melting of a common source. The fractionation sequence consists of two parts: an initial iron enrichment trend followed by a silica enrichment trend. We interpret the trace element data to indicate that the most iron rich lavas represent about 32% residual liquid derived by crystal fractionation of plagioclase, clinopyroxene, and lesser olivine from a parental magma with an Mg number of about 66. The silica enrichment trend results from crystallization of titanomagnetite and some apatite. Fractionation of pigeonite, which is a minor phase in the major element models, cannot be distinguished from clinopyroxene fractionation by using trace elements.

Journal of Geophysical Research Solid Earth↗
USGS2012

The Redskin granite: Evidence for thermogravitational diffusion in a Precambrian granite batholith

The 1000-Ma-old composite Pikes Peak batholith consists largely of medium- to coarse-grained biotite granite but is intruded by several smaller granite plutons. Two of these, the Tarryall Mountains batholith and the Redskin Granite (Hawley, 1969; Hawley and Wobus, 1977), are here interpreted as cupolas atop the main mass of the Pikes Peak batholith. Chemical analysis for 38 elements in rocks from the Redskin Granite, the Tarryall Mountains batholith, and the biotite granite of the Pikes Peak batholith support the contention that the compositions of the smaller plutons are unlikely to have evolved from that of the main mass of the batholith by crystal fractionation. It is believed that these compositions arose through liquid-state fractionation processes, possibly the convection-aided thermogravitational model of Shaw et al. (1976). It is important to recognize that the chemical signature of this fractionation process can be discerned in holocrystalline rocks as well as in their glassy equivalents because the process may be instrumental in the formation of some kinds of lithophile ore deposits.

Journal of Geophysical Research Solid Earth↗
USGS2012

Chemistry of rock-forming minerals of the Cretaceous-Paleocene batholith in southwestern Japan and implications for magma genesis

Petrographic descriptions and electron microprobe analyses of minerals are presented for 35 specimens from seven suites chosen to examine the transition from magnetite series to ilmenite series granitoids along two transects across the Cretaceous-Paleocene Inner Zone batholith of southwestern Japan. Regularities in chemical compositions of amphiboles, biotites, and feldspars suggest that fundamentally similar processes produced the magmas that formed the two series. Constant or decreasing Fe/(Fe+Mg) for biotites and amphiboles with increasing host rock silica content, coupled with the absence of early formed magnetite and sphene, suggest that magnetite series rocks may have become oxidized during crystallization near the level of intrusion, through the processes of second boiling and differential loss of hydrogen. For the Daito-Yokota, magnetite series suite, Fe/(Fe+Mg) for biotites decreased from 0.48 to 0.37 as SiO 2 content of the host rock increased from 55.3 to 75.5 wt %; for an ilmenite series suite from the Takanawa Peninsula, Fe/(Fe+Mg) for biotites increased from 0.51 to 0.77 with an increase in host rock SiO 2 from 53.4 to 75.5. Detailed consideration of amphibole chemistry shows predominance of edenitic and tschermakitic substitution schemes as well as coupling between substitutions of Ti in octahedral sites and Al IV . Interrelations between amphibole and biotite chemistry show that Fe/(Fe+Mg) and Mn contents can be interpreted in terms of equilibration, whereas Ti content cannot. The chemistry of chlorites correlates well with that of biotites; primary and secondary muscovites are distinct in composition. Plagioclase in all studied suites shows igneous zoning appropriate to host rock composition; perthitic alkali feldspars in all samples have lost albite component, and temperatures based on the two-feldspar geothermometer are low. The biotite-apatite geothermometer is also inoperative for this group of samples because fully fluorinated apatites typically occur in biotites of modest F content. Whereas magnetites have reequilibrated, analyses of ilmenites for the representative Daito-Yokota and Takanawa suites corroborate biotite compositional data and suggest that f O2 probably differed by 2 to 3 orders of magnitude during crystallization of silica-rich magnetite and ilmenite series granites. Whole-rock chemistry supports mineral chemistry in suggesting that the studied granitoids have crystallized from magmas generated in a lower crustal environment in which mantle-derived magmas partially melted source materials with igneous characteristics.

Journal of Geophysical Research Solid Earth↗
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↗
USGS2012

Petrogenesis of garnet two-mica granites in the Ruby Mountains, Nevada

Garnet-two-mica granites of Jurassic (160 ± 3 m.y.) and Cretaceous (83 ± 1.3 m.y.) ages intrude amphibolite facies metasedimentary rocks of Precambrian age and lower Paleozoic sedimentary rocks in the northern part of the Ruby Mountains, Nevada. High initial 87 Sr/ 86 Sr and high values of δ 18 O for minerals suggest a dominantly sedimentary source for the magmas. Total rare-earth contents are low and less fractionated (Ce N /Yb N = 10) than granitic rocks of similar composition in the Sierra Nevada and show a negative Eu anomaly (Eu/Eu* = 0.57). These data are consistent with melting in the stability field of plagioclase with little contribution to the melt from refractory minerals such as zircon, sphene, and apatite. Minor euhedral garnets are manganese-rich and magnesium-poor (al 60 sp 34 py 3.6 gr 2.8 ) compared to garnet (al 79 sp 6 py 13 gr 2 ) in the intruded metamorphic rocks and show an increase in manganese and a decrease in calcium in the outer few microns. Fractionation of magnesium-iron between garnet and biotite, and of sodium-potassium between K-feldspar and plagioclase in the granites suggest submagmatic equilibration temperatures, in the range of 365° to 505°C, whereas oxygen isotope equilibration temperatures for two granites are higher at 480° and 570°C. The compositions of the granites plot near minimum melting compositions in the water-saturated quartz-albite-orthoclase system. Comparison of muscovite plus quartz stability, the water-saturated granite solidus, and experimental garnet-melt equilibria suggests pressures of crystallization no lower than about 3.5 kbar. This pressure is compatible with that estimated from garnet-plagioclase-sillimanite-quartz equilibria for the intruded metasedimentary rocks.

Nevada↗
USGS2012

Two-mica granites of northeastern Nevada

The field settings are described and analytical data are presented for six two-mica granites from north-eastern Nevada. High δ 18 O and 87 Sr/ 86 Sr values indicate that all are S-type granites, derived from continental crust. The major element chemistry and accessory mineral contents of these rocks also are characteristic of S-type granites. Chemical, X ray, and other data are presented for the micas recovered from these granites. The muscovites are notably high in Fe 2 O 3 , FeO, and MgO. Except for one hydrobiotite, each of the biotites has an MgO content near 6.0 weight percent. Two different types of two-mica granites are recognized in the area of this study. One type is distinguished by the presence of many biotite euhedra within muscovite phenocrysts and by an unusual suite of accessory minerals completely devoid of opaque oxides. This type probably resulted from anatexis of late Precambrian argillites under conditions of relatively low oxygen fugacity, along a line that roughly coincides with the westward disappearance of continental basement. In the other textural type of two-mica granite the micas are equigranular and there is a greater variety of accessory minerals. The magmatic evolution of this type also appears to reflect the influence of late Precambrian argillites; there may be age differences between the two types of two-mica granites.

Journal of Geophysical Research Solid Earth↗
USGS2012

A stochastic fault model. 2. Time-dependent case

A random model of fault motion in an earthquake is formulated by assuming that the slip velocity is a random function of position and time truncated at zero, so that it does not have negative values. This random function is chosen to be self-affine; that is, on change of length scale, the function is multiplied by a scale factor but is otherwise unchanged statistically. A snapshot of slip velocity at a given time resembles a cluster of islands with rough topography; the final slip function is a smoother island or cluster of islands. In the Fourier transform domain, shear traction on the fault equals the slip velocity times an impedance function. The fact that this impedance function has a pole at zero frequency implies that traction and slip velocity cannot have the same spectral dependence in space and time. To describe stress fluctuations of the order of 100 bars when smoothed over a length of kilometers and of the order of kilobars at the grain size, shear traction must have a one-dimensional power spectrum is space proportional to the reciprocal wave number. Then the one-dimensional power spectrum for the slip velocity is proportional to the reciprocal wave number squared and for slip to its cube. If slip velocity has the same power law spectrum in time as in space, then the spectrum of ground acceleration will be flat (white noise) both on the fault and in the far field.

Journal of Geophysical Research Solid Earth↗
USGS2012

Geodetic measurement of crustal deformation on the San Andreas, Hayward, and Calaveras faults near San Francisco, California

Analysis of a geodetic network of 115 lines crossing the San Andreas, Hayward, and Calaveras faults in the vicinity of San Francisco Bay and measured repeatedly between 1970 and 1980 has revealed details about the accommodation of relative plate motion in this area. The most striking result is that the deformation is not uniformly distributed across the area. In the east bay, along the Hayward and Calaveras faults, all motion appears to take place as slip directly on the fault, with no accumulation of strain in the adjacent crust. On both the Calaveras and the Hayward faults the rate obtained for the 1970 to 1980 period agrees with geologic rates spanning a few million years and with creep rates spanning a few decades. The Hayward fault slip rate is 7 ± 1 mm/yr. The Calaveras fault slip rate is 7 ± 1 mm/yr, with perhaps half of this slip distributed across a zone a few kilometers wide, probably as inelastic deformation of weak near‐surface material. The absence of strain accumulation in the east bay is surprising since the Hayward and Calaveras faults have been the site of large earthquakes in the past. A block located east of the Calaveras fault and south of the Las Positas fault has been rotating clockwise at a rate of 0.3 ± 0.1 μrad/yr with very little internal deformation. Along the San Francisco peninsula no detectable slip occurs (less than 1.5 mm/yr) at the surface, but appreciable strain is accumulating. Near fault shear strain rates are 0.6 ± 0.1 μstrain/yr (engineering) with direction N47°W ± 9. The slip rate near the San Andreas fault is 12.2 ± 3.9 mm/yr distributed across a broad zone. The relative motion across the whole region during the period 1970–1980 is 32.1 ± 7.4 mm/yr.

California↗
USGS2012

Application of wave field continuation to the inversion of refraction data

Three examples of the inversion of refraction data by downward continuation illustrate the applicability of the method to field data. The first example is a refraction profile from the Mojave Desert, California. These data are spatially aliased and contain clear evidence of lateral inhomogeneity. The inversion in this case produces a broken image in the slowness‐depth domain due to the lateral inhomogeneity, but a useful average velocity model is still obtained. The second example is a shallow marine reflection profile. Here, the truncation effects due to the finite horizontal aperture of the recording cable produce artifacts in the slowness‐depth domain. The velocity model is, however, distinct from these artifacts, and the presence of strong precritical reflections aids in the inversion. The third example is another shallow marine reflection profile. The inversion of these data illustrate the utilization of constraints provided by multiples as well as primary arrivals.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Geomagnetic local and regional harmonic analyses

Procedures are developed for using rectangular and cylindrical harmonic analyses in local and regional areas. Both the linear least squares analysis, applicable when component data are available, and the nonlinear least squares analysis, applicable when only total field data are available, are treated. When component data are available, it is advantageous to work with residual fields obtained by subtracting components derived from a harmonic potential from the observed components. When only total field intensity data are available, they must be used directly. Residual values cannot be used. Cylindrical harmonic analyses are indicated when fields tend toward cylindrical symmetry; otherwise, rectangular harmonic analyses will be more advantageous. Examples illustrating each type of analysis are given.

Journal of Geophysical Research Solid Earth↗
USGS2012

The Yellowstone-Snake River Plain seismic profiling experiment: Crustal structure of the Eastern Snake River Plain

Seismic refraction profiles recorded along the eastern Snake River Plain (ESRP) in southeastern Idaho during the 1978 Yellowstone-Snake River Plain cooperative seismic profiling experiment are interpreted to infer the crustal velocity and attenuation (Q-1) structure of the ESRP. Travel-time and synthetic seismogram modeling of a 250 km reversed refraction profile as well as a 100 km detailed profile indicate that the crust of the ESRP is highly anomalous. Approximately 3 to 6 km of volcanic rocks (with some interbedded sediments) overlie an upper-crustal layer (compressional velocity ≅6.1 km/s) which thins southwestward along the ESRP from a thickness of 10 km near Island Park Caldera to 2 to 3 km beneath the central and southwestern portions of the ESRP. An intermediate-velocity (≅6.5 km/s) layer extends from ≅10 to ≅20 km depth. A thick (≅22 km) lower crust of compressional velocity 6.8 km/s, a total crustal thickness of ≅42 km, and a P n velocity of ≅7.9 km/s is observed in the ESRP, similar to the western Snake River Plain and the Rocky Mountains Provinces. High attenuation is evident on the amplitude corrected seismic data due to low-Q values in the volcanic rocks (Q p = 20 to 200) and throughout the crust (Q p = 160 to 300). Based on these characteristics of the crustal structure and volcanic-age progression data, it is suggested that the ESRP has resulted from an intensive period of intrusion of mantle-derived basaltic magma into the upper crust generating explosive silicic volcanism and associated regional uplift and caldera collapse. This activity began about 15 m.y. ago in southwestern Idaho and has migrated northeast to its present position at Yellowstone. Subsequent cooling of the intruded upper crust results in the 6.5 km/s velocity intermediate layer. Crustal subsidence and periodic basaltic volcanism as represented by the ESRP complete the sequence of crustal evolution.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Terranes and suture zones in east central Alaska

East central Alaska, with its 17 terranes, forms a part of the mosaic of allochthonous terranes that join the North American and Siberian plates. These terranes range from continental and continental margin, such as the Tatonduk with its thick well-bedded succession of marine shelf rocks, to seamount, arc, and ocean floor terranes. The Yukon crystalline terrane, the largest described here, is a composite of at least four subterranes juxtaposed across the Tintina fault with the Tatonduk terrane, a northwestern extension of the North American plate in Alaska. Inboard of the Yukon crystalline terrane are packets of closely appressed microterranes separated from the Tatonduk and other terranes belonging to North America by major suture zones. These microterranes lie between North America and the mosaic of accretionary terranes that form the more southerly part of Alaska. The most obviously allochthonous microterranes within the suture zones are the Woodchopper Canyon, an Early Devonian basaltic seamount, and the White Mountains, an Ordovician volcanic arc terrane capped by Silurian and Devonian carbonate bank deposits. The nearest counterpart of these terranes is the Alexander terrane in southeastern Alaska. The Tintina fault of Mesozoic and Cenozoic age, like the Denali fault, primarily follows old suture zones that separate terranes. Strike slip faulting developed after collision in places where further convergence was oblique to the terrane margins. Where terranes met head-on, their leading edges lie along a multiple set of high-angle faults that outline microterranes in accretion zones.

Alaska↗
USGS2012

Post 12 m.y. rotation of southwest Washington

Paleomagnetic field directions from the basalt of Pack Sack Lookout are compared to those from the Pomona Member of the Saddle Mountains Basalt of the Columbia River Basalt Group. The Pomona crops out over a wide region on the Columbia Plateau east of the Cascade Range, and the basalt of Pack Sack Lookout crops out well to the west of the Cascades about 30 to 60 km east of the Washington coast. Our paleomagnetic results support existing geologic and geochemical evidence that both these units are part of a single 12‐m.y.‐old flow that erupted in western Idaho and flowed to the Washington coast. The magnetic data further indicate that southwest Washington has undergone clockwise rotation with respect to the Columbia Plateau and stable North America. The data suggest that both a large‐scale regional rotation of ∼16° and locally complex small‐scale rotations exist, the two being present to different degrees in the eastern and western regions of southwest Washington. The Cascade Range appears to coincide with the tectonic boundary, separating rotated and unrotated regions of Washington state. Rotation of southwest Washington appears to have been associated with the rotation of large pieces of western Oregon and southern Washington, possibly as a result of either a ridge‐push force from the Basin and Range province or shear along the Pacific‐North America plate boundary.

Washington↗