USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Solid Earth”

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 307 records · Page 17Linked to original sources

Late Cenozoic volcanism, subduction, and extension in the Lassen region of California, southern Cascade Range

Some 537 volcanic vents younger than 7 Ma are identified and these are classified into five age intervals and five compositional categories based on SiO2 content. Maps of vents by age and composition illustrate regionally representative volcanic trends. By 2 Ma, the eastern limit of volcanism had contracted westward toward the late Quaternary arc. Late Quaternary volcanism is concentrated around and north of the silicic Lassen volcanic center. The belt of most recent volcanism (25-0 ka) has been active since at least 2 Ma. Most mafic volcanism is calcalkaline basalt and basaltic andesite. Normal faults and linear groups of vents are evidence of widespread crustal extension throughout most of the Lassen region. NNW orientation of maximum horizontal stress is indicated. The Lassen volcanic region is thought to occur above the subducting Gorda North plate but also within a broad zone of distributed extension that occurs in the North American lithosphere east and southeast of the present Cascadia subduction zone. The scarcity of volcanic rocks older than 7 Ma suggests that a more compressive lithospheric stress regime prior to the late Miocene extensional episode may have suppressed volcanism. -from Authors

Journal of Geophysical Research Solid Earth↗

Quaternary extrusion rates of the Cascade Range, northwestern United States and southern British Columbia

Quaternary (2–0 Ma) extrusion rates change significantly along the Cascade Range volcanic arc. The extrusion rate north of Mount Rainier is about 0.21 km 3 km −1 m.y. −1 ; the rate in southern Washington and northern Oregon south to Mount Hood is about 1.6 km 3 km −1 m.y. −1 ; in central Oregon the rate is 3–6 km 3 km −1 m.y. −1 ; and in northern California, the rate is 3.2 km 3 km −1 m.y. −1 . Eruption style also changes along the arc but at latitudes different from rate changes. At the ends of the arc, volcanism is focused at isolated intermediate to silicic composite volcanoes. The composite volcanoes represent ∼30% of the total volume of the arc. Mafic volcanic fields partly ring some composite volcanoes, especially in the south. In contrast, volcanism is diffused in the middle of the arc, where numerous overlapping mafic shields and a few composite volcanoes have built a broad ridge. Contrasting eruption style may signify diffuse versus focused heat sources or may reflect changes in permeability to ascending magma along the arc.

Journal of Geophysical Research Solid Earth↗

Geophysical constraints on Washington convergent margin structure

Gravity and magnetic maps of western Washington reveal the lateral structure and fabric of the Washington Coast Range, Puget Basin, and southern Washington Cascade Range. The magnetic and gravity maps show large amplitude positive anomalies associated with the shallow but largely buried section of Washington Coast Range mafic rocks which are separated by negative anomalies over deep sedimentary basins. The positive anomalies indicate that the Coast Range mafic basement extends farther east than previously thought, at least as far east as the longitude of Seattle. Linear and steep gravity and magnetic gradients indicate many unmapped, often buried faults in the Washington Coast Range Province. Magnetic highs are also associated with mapped batholiths in the Cascade arc. Several magnetic highs observed east of the Coast Range rocks and west of these batholiths may be associated with buried Tertiary plutons or ophiolites. Two-dimensional gravity and magnetic modeling constrained with geological and other geophysical data indicate that the Coast Range Province rocks are about 1 km thick at the coast, thickening to as much as 30 km near their postulated eastern edge. A maximum boundary on the average density of the upper 15–20 km of the rocks that compose the Coast Range Province of 2920 kg/m 3 was established by the modeling, suggesting a composition largely of basalt and gabbro with little interbedded sediments. Under these rocks may be mantle or a subduction complex composed of dense mafic, ultramafic, and sedimentary rocks like that proposed to underlie Vancouver Island. Previous gravity models of the Washington margin include lower densities for the proposed subduction complex than for Vancouver Island, suggesting a lower component of mafic and ultramafic rocks than the rocks underlying Vancouver Island. However, my Washington model requires that the proposed subduction complex be more dense than the trench sediments and, therefore, that material denser than sediments be incorporated within it. The absence of continental mantle and the modeled wedge shape of the Coast Range Province upper crust suggest that erosion of the bottom of the overriding plate by subduction processes may have occurred.

Washington↗

Stratigraphic, lithologic, and major element geochemical constraints on magmatic evolution at Lassen volcanic center, California

The Lassen volcanic center is the most recent of several long-lived volcanic centers in the southernmost Cascade Range. These centers have erupted products ranging from basaltic andesite to rhyolite and are superimposed on a background of regional basaltic to andesitic volcanism. The evolution of the Lassen volcanic center is described in three stages. Stages I and II comprise the Brokeoff volcano, and 80 km 3 andesitic stratocone, active from 600 to 400 ka. Brokeoff volcano is compositionally equivalent to the regional basaltic andesite to andesite volcanism in the Lassen region and is the result of structurally controlled focusing of the diffuse regional mafic magmatism. Stage III comprises a silicic dome field and adjacent area of hybrid andesites and has a total volume of about 100 km 3 . Volcanism during stage III was episodic and is subdivided into four sequences of lithologically and temporally distinct lavas. Stage III began at 400 ka with a rhyolitic, caldera-forming pyroclastic eruption and chemically related lavas. Additional sequences of dacite erupted between 250–200 ka and 100–0 ka. Hybrid andesites erupted adjacent to the silicic dome field between 300 and 0 ka. Porphyritic andesite and dacite with high Al 2 O 3 , low TiO 2 , medium K 2 O and FeO/MgO ratios of 1.5–2.0 are the most abundant rock types in the Lassen volcanic center. However, the single most voluminous unit is sparsely phyric rhyolite pumice. In general, the lavas of Lassen volcanic center form a single coherent trend on major element variation diagrams and in pseudo-quaternary phase space, consistant with an origin either by fractional crystallization or magma mixing. In detail, however, the lack of systematic temporal change in silica and subtly crossing trends indicate a complex origin. A variety of statistically successful fractional crystallization models can be constructed that derive Brokeoff andesites from regional magmas. An important conclusion of the modeling is that if fractional crystallization is the process responsible for generation of Brokeoff andesite, then the parent magma must be low to medium K in geochemical affinity in order to explain the variation in K 2 O. However, although major element variation can be modeled by fractional crystallization, petrographic and stratigraphic evidence indicates that magma mixing is an important but subtle process in Brokeoff lavas and suggests that lavas evolved in small independent batches. Lavas erupted during stage III, while predominantly silicic, range from 53 to 75% SiO 2 . Disequilibrium mineral assemblages in the stage III lavas indicate that they are not directly derived from Brokeoff andesite by fractional crystallization. Mixing of silicic magma with regional mafic magma and disaggregation of andesitic quenched magmatic inclusions play dominant roles in the compositional diversity of stage III lavas.

Journal of Geophysical Research Solid Earth↗

Trace element and isotopic constraints on magmatic evolution at Lassen volcanic center

Magmatic evolution at the Lassen volcanic center (LVC) is characterized by a transition from predominantly andesitic to predominantly silicic volcanism with time. Magmas of the andesitic, or “Brokeoff phase” of volcanism range in composition from basaltic andesite io dacite, whereas those of silicic, or “Lassen phase” range in composition from basaltic andesite to rhyolite. The compositions of magmas from each phase define well organized but distinct variation trends. Compared with Brokeoff-phase magmas of similar SiO 2 content, most Lassen-phase magmas contain lower concentrations of most incompatible minor and trace elements. Based on the behavior of both incompatible and compatible trace elements, the geochemical trends defined by the Brokeoff-phase magmas cannot be ascribed merely to fractional crystallization from a single or multiple mafic parental magmas, Moreover, the Lassen-phase magmas cannot be derived from the Brokeoff-phase magmas by fractional crystallization. Rather, the geochemical trends that characterize each volcanic phase define arrays that primarily indicate mixing between well-homogenized silicic and heterogeneous mafic magmas. The distinctive mixing-dominated arrays for each volcanic phase manifest the generation and evolution of two physically distinct, but genetically related magma systems. The LVC magmas have Sr, Nd, and Pb isotopic characteristics that approximate two-component mixing arrays. One isotopic component is similar in composition to that of NE Pacific Ocean ridge and seamount basalts (“MORB” component), the other to mafic Mesozoic granitoids sampled from the neighboring KSamath and Sierra Nevada provinces (“KSN” component). The isotopic compositions of the most silicic LVC magmas lie within the ranges defined by the more mafic LVC magmas, which in turn lie within broad ranges defined by primitive mafic lavas sampled from the Lassen region. The lack of a correlation between the major element and isotopic compositions of LVC magmas seriously limits any model for magmatic evolution that relies on assimilation of old middle to upper crust by isotopically homogeneous mafic magmas during their ascent through the crust. Alternatively, the isotopic and geochemical uniformity of the most silicic magmas of the Brokeoff and Lassen phases suggests that they are well-homogenized partial melts. The likely source region for these silicic melts is the lower crust, which we envision to consist primarily of mafic igneous rocks that are similar in geochemical and isotopic diversity to the regional mafic lavas. Magmatic evolution at LVC can be viewed in terms of a series of mantle melting events that subsequently stimulated meiting in a progressively increasing volume of the lower crust. In general, the LVC magmas represent slightly fractionated mixtures of the mantle-derived mafic magmas and silicic partial melts of the lower crust, the latter melts increasing in relative proportion over the history of the volcanic center. The voluminous rhyolitic lavas and pyroclastic materials erupted during the early Lassen phase represent lower crustal melts thai pooled into sufficient volumes to avoid significant blending and dilution with intruding mafic magmas. The geochemical and isotopic heterogeneity of both Brokeoff-phase andésites and mafic magmatic inclusions in the silicic Lassen-phase magmas must in part record the corresponding diversity of the mantle-derived contributions to arc magmagenesis in this region.

Journal of Geophysical Research Solid Earth↗

Post-11,000-year volcanism at Medicine Lake Volcano, Cascade Range, northern California

Eruptive activity during the past 11,000 years at Medicine Lake volcano has been episodic. Eight eruptions produced about 5.3 km 3 of basaltic lava during an interval of a few hundred years about 10,500 years B.P. After a hiatus of about 6000 years, eruptive activity resumed with a small andesite eruption at about 4300 years B.P. Approximately 2.5 km 3 of lava with compositions ranging from basalt to rhyolite vented in nine eruptions during an interval of about 3400 years in late Holocene time. The most recent eruption occurred about 900 years B.P. A compositional gap in SiO 2 values of erupted lavas occurs between 58 and 63%. The gap is spanned by chilled magmatic inclusions in late Holocene silicic lavas. Late Holocene andesitic to rhyolitic lavas were probably derived by fractionation, assimilation, and mixing from high-alumina basalt parental magma, possibly from basalt intruded into the volcano during the early mafic episode. Many basaltic to andesitic lavas contain iron-rich crystals and have high FeO*/MgO (characteristics caused by mixing of high-alumina basalt with ferrobasalt liquid produced by fractionation of parental high-alumina basalt). When ferrobasalt and high-alumina basalt are contaminated with a granitic crustal component, a calc-alkaline trend is produced. Some eruptions have produced both tholeiitic and calc-alkaline compositions. The eruptive activity is probably driven by intrusions of basalt that occur during east-west stretching of the crust in an extensional tectonic environment. Vents are typically aligned parallel or subparallel to major structural features, most commonly within 30° of north. Intruded magma should provide adequate heat for commercial geothermal development if sufficient fluids can be found. The nature and timing of future volcanic activity cannot be predicted from the observed pattern, but eruptions high on the edifice could produce high-silica products that might be accompanied by explosive activity, whereas eruptions lower on the flanks are likely to vent more fluid mafic lavas.

Journal of Geophysical Research Solid Earth↗

Constraints on the tectonics of the Mule Mountains thrust system, southeast California and southwest Arizona

The Mule Mountains thrust system crops out discontinuously over a 100-km-strike length in the Blythe-Quartzsite region of southeast California and southwest Arizona. Along the thrust system, middle and upper crustal metamorphic and plutonic rocks of Proterozoic and Mesozoic age are thrust north-northeastward (015° to 035°) over a lower plate metamorphic terrane that formed part of the Proterozoic North American craton, its Paleozoic sedimentary rock cover, overlying Mesozoic volcanic and sedimentary rocks, and the intruding Jurassic and Cretaceous granitic rocks. Stratigraphic, petrologic, and Pb isotopic ties for Jurassic granitoids and for Jurassic(?) and Cretaceous sedimentary rocks across the various parts of the thrust system indicate that related crustal blocks are superposed and preclude it from having large displacements. The thick-skinned thrust system is structurally symmetrical along its length with a central domain of synmetamorphic thrust faults that are flanked by western and eastern domains where lower plate synclines underlie the thrusts. Deformation occurred under low greenschist facies metamorphic conditions in the upper crust. Movement along the thrust system was probably limited to no more than a few tens of kilometers and occurred between 79±2 Ma and 70±4 Ma. The superposition of related rocks and the geometry of the thrust system preclude it from being a major tectonic boundary of post-Middle Jurassic age, as has been previously proposed. Rather, the thrust system forms the southern boundary of the narrow zone of Cretaceous intracratonic deformation, and it is one of the last tectonic events in the zone prior to regional cooling.

Journal of Geophysical Research Solid Earth↗

Comparison of ground motion from tremors and explosions in deep gold mines

Seismic body waves, from tamped chemical explosions, two with yields of 50 and one of 150 kg, were compared with corresponding data from three mining-induced tremors with a view to testing methods of discriminating between the two types of events. Detonated at depths of about 2 km, all three explosions generated P waves for which the low-frequency spectral asymptotes agree well with corresponding results scaled down from nuclear shots at the Nevada Test Site. For the two smaller explosions, recorded underground at hypocentral distances ranging from 234 to 871 m, the P wave corner frequencies, defined by the intersections of the low- and high-frequency asymptotes, are also in fine agreement with scaled-down results from the Nevada Test Site; for the 150-kg explosion, which was recorded only at the surface, it appears that 2.39 km of upward propagation caused the corner frequency to be at least a factor of 4 lower than anticipated due to attenuation. All three explosions generated S waves that appear to be a consequence of deviatoric stress release in the immediate environs of the explosions. The three tremors analyzed here generated P and S waves whose spectra and source parameters agree well with standard earthquake source models. With regard to discrimination based on P wave spectra, we conclude that for events of fixed low-frequency spectral asymptotes, the explosions typically have higher corner frequencies than tremors or earthquakes, although counterexamples certainly exist. Interestingly, the 150-kg explosion was identified as such on the basis of P and S wave polarities that are incompatible with the normally expected double-couple source model; instead, these initial motions are consistent with an explosion in conjunction with normal faulting. The body wave spectra of this explosion and those of a nearby tremor, however, were indistinguishable.

Journal of Geophysical Research Solid Earth↗

Infragravity waves over a natural barred profile

Measurements of cross-shore flow were made across the surf zone during a storm as a nearshore bar became better developed and migrated offshore. Measured infragravity band spectra were compared to synthetic spectra calculated numerically over the natural barred profile assuming a white run-up spectrum of leaky mode or high-mode edge waves. As in earlier studies, the spectra compared closely; however, for some frequencies the energy of the measured spectrum exceeded the energy of the synthetic spectrum, suggesting that the run-up spectrum was not white but had dominating frequencies. Utilizing cross-shore flow data and synthetic spectra from a number of cross-shore locations, an equivalent run-up spectrum was calculated for each day. On the first day of the storm, the equivalent run-up spectrum indicated a dominant wave that had a node in velocity reasonably close to the bar crest. Later during the storm, when the bar had migrated farther offshore, there was no evidence for a dominant motion having a velocity node at the bar crest. The structure of the equivalent run-up spectrum compared well with spectra of direct measurements of run-up obtained several hundred meters away. We have no clear evidence in support of the theory that infragravity waves might form or force the offshore migration of a bar. To confirm this finding, longer records obtained synoptically over a developing bar are required. The dominant wave observed early in the storm was consistent with Symond and Bowen's (1984) theoretical prediction of resonant amplification of discrete frequencies over a barred profile.

Journal of Geophysical Research Solid Earth↗

Effect of crustal layering upon dislocation modeling

Slip distribution at depth on a fault may be inferred from the deformation observed on the surface. In inverting the surface deformation data to obtain the slip distribution, the Earth is generally approximated by an elastic half‐space. Slip distributions inferred from a half‐space model may contain artifacts, including zones of reversed slip, due solely to effects of layering in the real Earth. This effect is demonstrated for a vertical strike‐slip fault in an Earth consisting of an elastic layer overlying an elastic half‐space. Slip on the fault is taken to be independent of the along‐strike coordinate (i.e., antiplane strain is assumed). For a given slip distribution in this model the slip distribution on a similar fault in an elastic half‐space is found that produces the identical surface deformation. Comparison of the two slip distributions reveals structure introduced into the half‐space equivalent slip profile by crustal layering. The comparisons suggest that low‐resolution inversion schemes (e.g., single screw dislocation models) are not drastically affected by Earth structure, but attempts at detailed inversion are likely to produce profiles contaminated by artifacts of Earth structure.

Journal of Geophysical Research B: Solid Earth↗

A migratory mantle plume on Venus: Implications for Earth?

A spatially fixed or at least internally rigid hotspot reference frame has been assumed for determining relative plate motions on Earth. Recent 1:5,000,000 scale mapping of Venus, a planet without terrestrial-style plate tectonics and ocean cover, reveals a systematic age and dimensional progression of corona-like arachnoids occurring in an uncinate chain. The nonrandom associations between arachnoids indicate they likely formed from a deep-seated mantle plume in a manner similar to terrestrial hotspot features. However, absence of expected convergent "plate" margin deformation suggests that the arachnoids are the surface expression of a migratory mantle plume beneath a stationary surface. If mantle plumes are not stationary on Venus, what if any are the implications for Earth?

Journal of Geophysical Research B: Solid Earth↗

Premonitory slip and tidal triggering of earthquakes

We have conducted a series of laboratory simulations of earthquakes using granite cylinders containing precut bare fault surfaces at 50 MPa confining pressure. Axial shortening rates between 10 −4 and 10 −6 mm/s were imposed to simulate tectonic loading. Average loading rate was then modulated by the addition of a small-amplitude sine wave to simulate periodic loading due to Earth tides or other sources. The period of the modulating signal ranged from 10 to 10,000 s. For each combination of amplitude and period of the modulating signal, multiple stick-slip events were recorded to determine the degree of correlation between the timing of simulated earthquakes and the imposed periodic loading function. Over the range of parameters studied, the degree of correlation of earthquakes was most sensitive to the amplitude of the periodic loading, with weaker dependence on the period of oscillations and the average loading rate. Accelerating premonitory slip was observed in these experiments and is a controlling factor in determining the conditions under which correlated events occur. In fact, some form of delayed failure is necessary to produce the observed correlations between simulated earthquake timing and characteristics of the periodic loading function. The transition from strongly correlated to weakly correlated model earthquake populations occurred when the amplitude of the periodic loading was approximately 0.05 to 0.1 MPa shear stress (0.03 to 0.06 MPa Coulomb failure function). Lower-amplitude oscillations produced progressively lower correlation levels. Correlations between static stress increases and earthquake aftershocks are found to degrade at similar stress levels. Typical stress variations due to Earth tides are only 0.001 to 0.004 MPa, so that the lack of correlation between Earth tides and earthquakes is also consistent with our findings. A simple extrapolation of our results suggests that approximately 1% of midcrustal earthquakes should be correlated with Earth tides. Triggered seismicity has been reported resulting from the passage of surface waves excited by the Landers earthquake. These transient waves had measured amplitudes in excess of 0.1 MPa at frequencies of 0.05 to 0.2 Hz in regions of notable seismicity increase. Similar stress oscillations in our laboratory experiments produced strongly correlated stick-slip events. We suggest that seemingly inconsistent natural observations of triggered seismicity and absence of tidal triggering indicate that failure is amplitude and frequency dependent. This is the expected result if, as in our laboratory experiments, the rheology of the Earth's crust permits delayed failure.

Journal of Geophysical Research B: Solid Earth↗

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 1. Strong motions

We develop a methodology to perform finite fault source inversions from strong motion data using Green's functions (GFs) calculated for a three-dimensional (3-D) velocity structure. The 3-D GFs are calculated numerically by inserting body forces at each of the strong motion sites and then recording the resulting strains along the target fault surface. Using reciprocity, these GFs can be recombined to represent the ground motion at each site for any (heterogeneous) slip distribution on the fault. The reciprocal formulation significantly reduces the required number of 3-D finite difference computations to at most 3NS, where NS is the number of strong motion sites used in the inversion. Using controlled numerical resolution tests, we have examined the relative importance of accurate GFs for finite fault source inversions which rely on near-source ground motions. These experiments use both 1-D and 3-D GFs in inversions for hypothetical rupture models in order (1) to analyze the ability of the 3-D methodology to resolve trade-offs between complex source phenomena and 3-D path effects, (2) to address the sensitivity of the inversion results to uncertainties in the 3-D velocity structure, and (3) to test the adequacy of the 1-D GF method when propagation effects are known to be three-dimensional. We find that given "data" from a prescribed 3-D Earth structure, the use of well-calibrated 3-D GFs in the inversion provides very good resolution of the assumed slip distribution, thus adequately separating source and 3-D propagation effects. In contrast, using a set of inexact 3-D GFs or a set of hybrid 1-D GFs allows only partial recovery of the slip distribution. These findings suggest that in regions of complex geology the use of well-calibrated 3-D GFs has the potential for increased resolution of the rupture process relative to 1-D GFs. However, realizing this full potential requires that the 3-D velocity model and associated GFs should be carefully validated against the true 3-D Earth structure before performing the inverse problem with actual data.

Journal of Geophysical Research B: Solid Earth↗

Geyser periodicity and the response of geysers to deformation

Numerical simulations of multiphase fluid and heat transport through a porous medium define combinations of rock properties and boundary conditions which lead to geyser‐like periodic discharge. Within the rather narrow range of conditions that allow geyser‐like behavior, eruption frequency and discharge are highly sensitive to the intrinsic permeabilities of the geyser conduit and the surrounding rock matrix, to the relative permeability functions assumed, and to pressure gradients in the matrix. In theory, heat pipes (concomitant upward flow of steam and downward flow of liquid) can exist under similar conditions, but our simulations suggest that the periodic solution is more stable. Simulated time series of geyser discharge are chaotic, but integrated quantities such as eruption frequency and mass discharge per eruption are free of chaos. These results may explain the observed sensitivity of natural geysers to small strains such as those caused by remote earthquakes, if ground motion is sufficient to induce permeability changes. Changes in geyser behavior caused by minor preseismic deformation, periodic surface loading, and Earth tides are more difficult to explain in the context of our current model.

Journal of Geophysical Research B: Solid Earth↗

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 2. Combining seismic and geodetic data

Using numerical tests for a prescribed heterogeneous earthquake slip distribution, we examine the importance of accurate Green's functions (GF) for finite fault source inversions which rely on coseismic GPS displacements and leveling line uplift alone and in combination with near-source strong ground motions. The static displacements, while sensitive to the three-dimensional (3-D) structure, are less so than seismic waveforms and thus are an important contribution, particularly when used in conjunction with waveform inversions. For numerical tests of an earthquake source and data distribution modeled after the 1994 Northridge earthquake, a joint geodetic and seismic inversion allows for reasonable recovery of the heterogeneous slip distribution on the fault. In contrast, inaccurate 3-D GFs or multiple 1-D GFs allow only partial recovery of the slip distribution given strong motion data alone. Likewise, using just the GPS and leveling line data requires significant smoothing for inversion stability, and hence, only a blurred vision of the prescribed slip is recovered. Although the half-space approximation for computing the surface static deformation field is no longer justifiable based on the high level of accuracy for current GPS data acquisition and the computed differences between 3-D and half-space surface displacements, a layered 1-D approximation to 3-D Earth structure provides adequate representation of the surface displacement field. However, even with the half-space approximation, geodetic data can provide additional slip resolution in the joint seismic and geodetic inversion provided a priori fault location and geometry are correct. Nevertheless, the sensitivity of the static displacements to the Earth structure begs caution for interpretation of surface displacements, particularly those recorded at monuments located in or near basin environments.

Journal of Geophysical Research B: Solid Earth↗

Mapping playa evaporite minerals and associated sediments in Death Valley, California, with multispectral thermal infrared images

Efflorescent salt crusts and associated sediments in Death Valley, California, were studied with remote-sensing data acquired by the NASA thermal infrared multispectral scanner (TIMS). Nine spectral classes that represent a variety of surface materials were distinguished, including several classes that reflect important aspects of the playa groundwater chemistry and hydrology. Evaporite crusts containing abundant thenardite (sodium sulfate) were mapped along the northern and eastern margins of the Cottonball Basin, areas where the inflow waters are rich in sodium. Gypsum (calcium sulfate) crusts were more common in the Badwater Basin, particularly near springs associated with calcic groundwaters along the western basin margin. Evaporite-rich crusts generally marked areas where groundwater is periodically near the surface and thus able to replenish the crusts though capillary evaporation. Detrital silicate minerals were prevalent in other parts of the salt pan where shallow groundwater does not affect the surface composition. The surface features in Death Valley change in response to climatic variations on several different timescales. For example, salt crusts on low-lying mudflats form and redissolve during seasonal-to-interannual cycles of wetting and desiccation. In contrast, recent flooding and erosion of rough-salt surfaces in Death Valley probably reflect increased regional precipitation spanning several decades. Remote-sensing observations of playas can provide a means for monitoring changes in evaporite facies and for better understanding the associated climatic processes. At present, such studies are limited by the availability of suitable airborne scanner data. However, with the launch of the Earth Observing System (EOS) AM-1 Platform in 1998, multispectral visible/near-infrared and thermal infrared remote-sensing data will become globally available.

Journal of Geophysical Research B: Solid Earth↗

CRUST 5.1: A global crustal model at 5° x 5°

We present a new global model for the Earth's crust based on seismic refraction data published in the period 1948-1995 and a detailed compilation of ice and sediment thickness. An extensive compilation of seismic refraction measurements has been used to determine the crustal structure on continents and their margins. Oceanic crust is modeled with both a standard model for normal oceanic crust, and variants for nonstandard regions, such as oceanic plateaus. Our model (CRUST 5.1) consists of 2592 5° × 5° tiles in which the crust and uppermost mantle are described by eight layers: (1) ice, (2) water, (3) soft sediments, (4) hard sediments, (5) crystalline upper, (6) middle, (7) lower crust, and (8) uppermost mantle. Topography and bathymetry are adopted from a Mandard database (ETOPO-5). Compressional wave velocity in each layer is based on field measurements, and shear wave velocity and density are estimated using recently published empirical V p -V s and V p -density relationships. The crustal model differs from previous models in that (1) the thickness and seismic/density structure of sedimentary basins is accounted for more completely, (2) the velocity structure of unmeasured regions is estimated using statistical averages that are based on a significantly larger database of crustal structure, (3) the compressional wave, shear wave, and density structure have been explicitly specified using newly available constraints from field and laboratory studies. Thus - this global crustal model is based on substantially more data than previous models and differs from them in many important respects. A new map of the thickness of the Earth's crust is presented, and we illustrate the application of this model by using it to provide the crustal correction for surface wave phase velocity maps. Love waves at 40 s are dominantly sensitive to crustal structure, and there is a very close correspondence between observed phase velocities at this period and those predicted by CRUST 5.1. We find that the application of crustal corrections to long-period (167 s) Rayleigh waves significantly increases the variance in the phase velocity maps and strengthens the upper mantle , velocity anomalies beneath stable continental regions. A simple calculation of crustal isostacy indicates significant lateral variations in upper mantle density. The model CRUST 5.1 provides a complete description of the physical properties of the Earth's crust at a scale of 5° × 5° and can be used for a wide range of seismological and nonseismological problems.

Journal of Geophysical Research B: Solid Earth↗

Apparent break in earthquake scaling due to path and site effects on deep borehole recordings

We reexamine the scaling of stress drop and apparent stress, rigidity times the ratio between seismically radiated energy to seismic moment, with earthquake size for a set of microearthquakes recorded in a deep borehole in Long Valley, California. In the first set of calculations, we assume a constant Q and solve for the corner frequency and seismic moment. In the second set of calculations, we model the spectral ratio of nearby events to determine the same quantities. We find that the spectral ratio technique, which can account for path and site effects or nonconstant Q, yields higher stress drops, particularly for the smaller events in the data set. The measurements determined from spectral ratios indicate no departure from constant stress drop scaling down to the smallest events in our data set (Mw 0.8). Our results indicate that propagation effects can contaminate measurements of source parameters even in the relatively clean recording environment of a deep borehole, just as they do at the Earth's surface. The scaling of source properties of microearthquakes made from deep borehole recordings may need to be reevaluated.

Journal of Geophysical Research B: Solid Earth↗