Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Oceans”

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

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

Upper crustal structure of the Mount Hood, Oregon, region as revealed by time term analysis

Seismic refraction data with a dense areal distribution were collected to study the seismic structure of Mount Hood and the surrounding region. This area is typical of Cascade volcanoes and is geologically quite complex. The prime goals of this project were to search for velocity variations in the upper crustal rocks and to determine if the velocity of these rocks is anisotropic. A new system, including 100 remote recording units, was developed to facilitate the collection of data in this type of survey. The data collected in this study reveal a large variation in velocity and thickness of the uppermost crustal rocks that is probably typical of the High Cascade province. A regional structural pattern surrounding Mount Hood, where there is a marked thinning of low-velocity near-surface rocks, suggests that the present edifice of Mount Hood lies on top of a much larger structure, possibly the roof of a large batholith that was emplaced prior to the eruption of the volcanic rocks that form the modern mountain. A straightforward time term analysis of the data reveals this structure. When the time terms are used to correct the travel time observations for the variation in structure under each station, the remaining set of residuals indicates a variation of travel time with azimuth. While a systematic error in the time term reduction cannot be ruled out, this variation is probably the result of anisotropic velocity structure in the underlying refractor. The alignment of the direction of maximum velocity with the direction of maximum principal stress and the decrease in anisotropy with depth suggest that the anisotropic velocity structure is related to the opening and closing of fractures in the present stress field. The direction of maximum velocity is about N25°W, and the degree of anisotropy appears to vary with depth from about 3% near the surface to zero at depths below 8 km.

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

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↗
USGS2012

Source parameters of the 1980 Mammoth Lakes, California, earthquake sequence

From the more than 1500 Mammoth Lakes earthquakes recorded on three-component digital seismographs (Spudich et al., 1981), 150 were used in an analysis of the locations, mechanism, and source parameters. A composite fault plane solution of nine earthquakes 3.9 ≤ M ≤ 5.1 defines a right-lateral strike slip mechanism on a steeply dipping nearly east-west plane striking S75°E or left-lateral strike slip on a nearly north-south plane striking N10°E. Vertical cross sections of well-located aftershocks indicate possibly three east-west planes that coincide with the locations of the four largest earthquakes with M L ≥ 6.0. Using the spectral analysis of S waves (Brune, 1970), source parameters for 67 earthquakes were determined. Forty-eight had magnitudes greater than or equal to 3.0. Seismic moments ranged from 9.20×10 18 dyn cm to 2.33×10 24 dyn cm. Earthquakes with seismic moment greater than about 1.0×10 21 dyn cm had nearly constant stress drops (≃ 50 bars); earthquakes with seismic moment less than about 1.0×10 21 dyn cm had stress drops that apparently decrease as seismic moment decreases.

California↗
USGS2012

Uranium-lead isotopic ages from the Sierra Nevada Batholith, California

This study provides new information on the timing and distribution of Mesozoic magmatic events in the Sierra Nevada batholithic complex chiefly between 36° and 37°N. latitude. U-Pb ages have been determined for 133 zircon and 7 sphene separates from 82 samples of granitoid rocks. Granitoid rocks in this area range in age from 217 to 80 m.y. Triassic intrusions are restricted to the east side of the batholith; Jurassic plutons occur south of the Triassic plutons east of the Sierra Nevada, as isolated masses within the Cretaceous batholith, and in the western foothills of the range; Cretaceous plutons form a continuous belt along the axis of the batholith and occur as isolated masses east of the Sierra Nevada. No granitic intrusions were emplaced for 37 m.y. east of the Sierra Nevada following the end of Jurassic plutonism. However, following emplacement of the eastern Jurassic granitoids, regional extension produced a fracture system at least 350 km long into which the dominantly mafic, calc-alkalic Independence dike swarm was intruded 148 m.y. ago. The dike fractures probably represents a period of regional crustal extension caused by a redistribution of the regional stress pattern accompanying the Nevadan orogeny. Intrusion of Cretaceous granitic plutons began in large volume about 120 m.y. ago in the western Sierra Nevada and migrated steadily eastward for 40 m.y. at a rate of 2.7 mm/y. This slow and constant migration indicates remarkably uniform conditions of subduction with perhaps downward migration of parent magma generation or a slight flattening of the subduction zone. Such steady conditions could be necessary for the production of large batholithic complexes such as the Sierra Nevada. The abrupt termination of plutonism 80 m.y. ago may have resulted from an increased rate of convergence of the American and eastern Pacific plates and dramatic flattening of the subduction zone. U-Pb ages of the Giant Forest-alaskite sequence in Sequoia National Park are all in the range 99±3 m.y., indicating a relatively short period of emplacement and cooling for this nested group of plutons. U-Pb ages of a mafic inclusion and its host granodiorite indicate that both were derived from a common source or that the mafic inclusion was totally equilibrated with the granodioritic magma. Comparison of isotopic ages determined by different methods such as zircon U-Pb, sphene U-Pb, hornblende K-Ar, and biotite K-Ar suggests that zircon U-Pb ages generally approximate the emplacement age of a pluton. However, some plutons probably contain inherited or entrained old zircons, and the zircons of some samples are disturbed by younger thermal and metamorphic events. The ages reported here are consistent with U-Pb age determinations previously made on granitic rocks to the north [Stern et al., 1981], The age distribution of granitic belts determined here is in general agreement with those established by K-Ar dating [Evernden and Kistler, 1970] but does not differentiate the five epochs of plutonism determined in their study.

California↗
USGS2012

Contemporary block tectonics: California and Nevada

Well-determined fault plane solution and the gross pattern of late-Cenozoic faulting in California and Nevada show a systematic relation between the orientation of fault planes and slip directions. In general, normal faults have northerly strikes, reverse faults have easterly strikes, and dextral and sinstral strike slip faults have northwesterly and northeasterly strikes, respectively. Kinematically, this relation is consistent with the response of clusters of fault-bounded crustal blocks to a regional stress field generated by the relative motion between the Pacific and North American plates. In this stress field, the greatest and least principal (compressive) stresses are restricted to northerly and easterly striking vertical planes, respectively. Clusters of crustal blocks bound by faults having the above attitudes form three basic regimes: (1) a spreading regime with a combination of normal and conjugate, strike slip faults, (2) a converging regime with a combination of reverse and conjugate, strike slip faults, and (3) a wrench regime with sets of subparallel, throughgoing, northwest striking (dextral) faults or northeast striking (sinstral) faults. These three regimes are typified by faulting patterns in the Basin and Range province (spreading), the Transverse Ranges (converging), and the San Andreas-Garlock fault systems (wrench), respectively. The gross deformation of each regime resulting from relative displacements between individual blocks is characterized by north-south shortening and east-west extension with the ratio of extensional to shortening strains (and the areal dilitation) decreasing systematically from spreading to wrench to compressional regimes. The wrench regime involves a component of net rotational deformation (clockwise for dextral slip and counter-clockwise for sinstral slip), while deformation of the spreading and converging regimes is irrotational. Local deviations from regional kinematic directions are concentrated along the boundaries between regimes reflecting the mismatch in gross deformation fields between regimes. Maximum principal and shear stress magnitudes will increase systematically from spreading to wrench to converging regimes provided that fault slip is controlled by frictional strength (Byerlee's law) along preexisting fractures and that pore pressure in the brittle crust is laterally uniform. A minimum strength difference between active, block-bounding faults and block interiors is 15–30%. Simple arrangements of such block clusters mimic the gross kinematic pattern of Quaternary faulting in California and Nevada. Some implications for contemporary tectonics emphasized by this model involve the westward displacement of the Sierra Nevada block.with respect to the stable interior of the North American plates, oblique thrusting of the Salinian block over the Pacific plate, and a progressive increase in the offset of the San Andreas fault represented by the ‘big bend’ through the Transverse Ranges.

Journal of Geophysical Research Solid Earth↗
USGS2012

The sea slope problem revisited

Discrepancies in sea surface topography based on comparisons between the results of steric leveling and repeated geodetic levelings have identified what is known as the ‘sea slope problem.’ This problem is actually twofold: (1) the sea surface relief based on steric leveling differs significantly from that based on geodetic leveling along several generally north-south coasts, and (2) successively propagated levelings between several widely separated tide stations indicate that the stationary sea slope seemingly has been changing with time, whereas differenced sea level means between these stations indicate that the sea surface relief has remained virtually invariant during the same intervals. Reexamination of the three reported discrepancies between geodetically and sterically determined sea slopes indicates that the Australian example is based on leveling of a quality inappropriate to the comparison. The discrepancy developed along the Atlantic coast of the United States is limited to the reach between Portsmouth–Hampton Roads, Virginia, and Key West, Florida, where the accuracy of steric leveling may be especially vulnerable owing to the dynamic effects of the Gulf Stream. Reconsideration of the example along the Pacific coast of the United States indicates that the various discrepancies are due chiefly to intrasurvey movement and resultant distortion of geodetically defined height differences between tide stations. Agreement between the results of steric and geodetic leveling along tectonically inactive north-south coasts devoid of strong boundary currents is generally good. This observation supports the conclusions (1) that any directionally dependent systematic error in geodetic leveling is measurably insignificant and (2) that where allowance is made for the possible effects of major boundary currents or intrasurvey movement during levelings between tide stations, the sea slope problem tends to vanish.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Earthquakes of Loihi submarine volcano and the Hawaiian hot spot

Loihi is an active submarine volcano located 35 km south of the island of Hawaii and may eventually grow to be the next and southernmost island in the Hawaiian chain. The Hawaiian Volcano Observatory recorded two major earthquake swarms located there in 1971–1972 and 1975 which were probably associated with submarine eruptions or intrusions. The swarms were located very close to Loihi's bathymetric summit, except for earthquakes during the second stage of the 1971–1972 swarm, which occurred well onto Loihi's southwest flank. The flank earthquakes appear to have been triggered by the preceding activity and possible rifting along Loihi's long axis, similar to the rift-flank relationship at Kilauea volcano. Other changes accompanied the shift in locations from Loihi's summit to its flank, including a shift from burst to continuous seismicity, a rise in maximum magnitude, a change from small earthquake clusters to a larger elongated zone, a drop in b value, and a presumed shift from concentrated volcanic stresses to a more diffuse tectonic stress on Loihi's flank. The 1971–1972 swarm began at depths of 20–50 km about 1 month before the shallow swarm started and suggests an upward migration of seismicity as the first stage of the 13-month swarm. The seismic ‘root’ of Kilauea volcano is well defined by earthquakes which plunge to the south and southwest to a depth of 50–60 km, terminating in a region in which earthquakes are associated with deep harmonic tremor and magma. The zone of these deeper earthquakes and tremor lies between Kilauea, Loihi, and Mauna Loa volcanoes, may feed magma to all three volcanoes, and probably locates the Hawaiian hot spot.

Journal of Geophysical Research Solid Earth↗
USGS2012

Rapid intrusion of magma into wet rock: Groundwater flow due to pore pressure increases

Analytical and numerical solutions are developed to simulate the pressurization, expansion, and flow of groundwater contained within saturated, intact host rocks subject to sudden heating from the planar surface of an igneous intrusion. For most rocks, water diffuses more rapidly than heat, assuring that groundwater is not heated along a constant-volume pressure path and that thermal expansion and pressurization adjacent to the intrusion drives a flow that extends well beyond the heated region. The forcing parameter for pressurization and flow is α Δ T , where α is a thermal expansion coefficient reflecting the overall expansion of water heated through the temperature difference Δ T between the initial ambient and intrusive values. Pore pressure increases due to heating are greatest when the intrusion is emplaced rapidly and where the intrusive contact is impervious to groundwater contained in stiff, impermeable rocks with high thermal diffusivities and porosities. The maximum velocity of water flowing in pores decays with the inverse square root of time and is insensitive to hydraulic properties of the host rocks. Pressures are lessened and flow directions are reversed with the onset of hydrothermal convection. This occurs at times ranging from hours to weeks after onset of intrusion. As magma rises into near-surface rocks, steam can be generated. Solutions indicate that pressure increases and velocities are sensitive to the overall amount of expansion rather than the behavior of the water-steam transition. Both the overall thermal expansion coefficient α and the temperature difference Δ T are greater in shallow (<1 km) environments than in deep (∼5 km) ones. Thus, for rocks with similar transport properties, pressure increases due to heating are greatest in shallow environments. Although solutions can be applied to rocks with a wide variety of properties, pressure increases are calculated for compliant quartz-rich sedimentary rocks with a porosities between 1 and 20% and permeabilities between 1 darcy and 1 μdarcy, subject to temperature increases of 500 and 1000 K at depths ranging from 0.1 to 5 km in a region of hydrostatic pressures and normal geothermal gradient. Such rocks, with porosities greater than 5%, permeabilities less than a 0.1 mdarcy, and drained hydrostatic compressibilities of 10 −4 /MPa, undergo pressure increases greater than 10 MPa (100 bars)for conditions typical of water table depths of 2.5 km and heating to 500 K above ambient. Similar rocks, but with permeabilities less than 1 mdarcy, undergo pressure increases of 10 MPa for conditions typical of 1 km water table depth. Rocks commonly considered to be good aquifers undergo pressure increases of less than 1 MPa, primarily because of their high permeability. Although these estimates neglect the effects of fracturing and brecciation that may accompany such pressure increases, calculations indicate that pressure increases due to heating of cool groundwater can lead to failure of host rocks by a phreatic mechanism.

Journal of Geophysical Research Solid Earth↗
USGS2012

Extension in the Rio Grande rift

A positive gravity anomaly along the axis of the Rio Grande rift reflects a volume of anomalous mass added at the base of the crust and intruded into the crust. Part of this volume can be associated with vertical uplift of the crust. The remainder of this anomalous volume, plus the volume of surficial graben fill, can be associated with horizontal crustal extension. The volume of crustal uplift in the Rio Grande rift is unknown, but this term can be eliminated by means of an independent equation provided by assumption of generalized isostatic equilibrium. The volume and mass equations combined provide a solution for extension of the crust in terms of the following parameters: total anomalous mass deficiency in the mantle lithosphere, total anomalous mass excess in the crust and its density contrast, total anomalous mass deficiency of surficial graben fill and its density contrast, and the volume of material eroded from the uplift. Using standard density estimates and masses determined by equivalent-source modeling of gravity profiles, I obtained 1-km extension at 37°N (Colorado-New Mexico border), 13-km extension at 35°N (Albuquerque, New Mexico), and 24-km extension at 33°N in southern New Mexico. These estimates are of necessity very poorly constrained and may be as much as ±60% in error. Taking results at face value, the best fitting Euler pole occurs at about 41°N, in north central Colorado. The calculated opening angle is 1.32°. For an approximately 30-m.y. duration of the extensional system (late Oligocene to the present), calculated average angular velocity is 7.8×10 −4 rad/m.y.; calculated average spreading half-rate at 33°N (for example) is a stately 0.04 cm/yr. The ratio of extension to uplift increases southward.

Journal of Geophysical Research Solid Earth↗