USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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,477 records · Page 82Linked to original sources

A seismic refraction survey of the Imperial Valley Region, California

The U.S. Geological Survey conducted an extensive seismic refraction survey in the Imperial Valley region of California in 1979. The Imperial Valley is located in the Salton Trough, an active rift between the Pacific and North American plates. Forty shots fired at seven shot points were recorded by 100 portable seismic instruments at typical spacing of 0.5–1 km. More than 1300 recording locations were occupied, and more than 3000 usable seismograms were obtained. We analyzed five profiles using a standard ray‐tracing program, constructed a contour map of reduced travel times from our most widely recorded shot point, and modeled an existing gravity profile across the Salton Trough. Results are itemized: (1) All models have in common a sedimentary layer ( V p = 1.8–5.0 km/s), a “transition zone” ( V p = 5.0–5.65 km/s), a basement ( V p = 5.65 km/s in the Imperial Valley, 5.9 km/s on the bordering mesas), and subbasement ( V p = 7.2 km/s). (2) The sedimentary layer ranges in thickness along the axis of the Salton Trough from 3.7 km (Salton Sea) to 4.8 km (U.S.‐Mexican border). On the bordering mesas it is quite variable in thickness. (3) The “transition” zone is about 1 km thick in most places. In the Imperial Valley there are no marked velocity discontinuities in this zone between the sedimentary layer and basement. On the bordering mesas, however, there is a discontinuity at the top of this zone. (4) There are apparently two types of basement. On the bordering mesas, basement is crystalline igneous and metamorphic rocks. In the Imperial Valley, basement is mostly lower‐greenshist‐facies sedimentary rocks, based primarily on the smooth transition in character from sediment to basement arrivals, the low value of basement velocity, and the fact that deep (4 km) wells in the valley penetrate only the upper part of the known Cenozoic stratigraphic column for the Salton Trough. (5) The subbasement, or intermediate crustal layer, ranges in depth along the axis of the Salton Trough from 16 km (Salton Sea) to 10 km (U.S.‐Mexican border). Gravity modeling requires that this layer deepen and/or pinch out beneath the bordering mesas and mountain ranges. Based on its high velocity and the presence of intrusive basaltic rocks in the sedimentary section in the Imperial Valley, the subbasement is thought to be a mafic intrusive complex similar to oceanic middle crust. (6) Several structures are seen that affect basement, transition zone, and deeper parts of the sedimentary layer. They include a scarp along the Imperial fault, as much as 1 km down to the northeast, and a scarp passing roughly along the topographic boundary between the Imperial Valley and the bordering mesa to the west, as much as 3½ km down to the east. We interpret the latter scarp to be the suture, or rift boundary, between the older crystalline basement on the mesa and the younger metasedimentary basement in the Imperial Valley. (7) On a contour map of reduced travel time from our most widely recorded shot point, subtle patches of early arrivals among otherwise late arrivals in the central Imperial Valley correlate well with known geothermal resource areas having reservoir temperatures of more than 150°C. Apparently the Salton Trough is a location where new crust is being generated. As the rift opens, mafic intrusive rocks fill it from below as sedimentary rocks fill it from above. Rifting and intrusion produce high heat flow that metamorphoses the sedimentary rocks to shallow depth (metasedimentary basement in the Imperial Valley) and thus consolidates the new crust.

California↗

Local gravity anomalies produced by dislocation sources

Rundle (1978) and Walsh and Rice (1979) have shown that the change in the vertical component of gravity is proportional to uplift for a spherical source of dilatation and for slip on an infinitely long dip-slip fault. In the first case, no free air gravity anomaly is produced and in the second case no Bouguer gravity anomaly. Gravity anomalies due to other dislocation sources in three dimensions are as follows: For strike-slip faulting the ratio of the gravity change to uplift depends upon position; however, the gravity change contours are roughly similar to those corresponding to a zero free air gravity anomaly. Nor is the ratio constant for dip-slip faulting except for the two special cases of dip slip on a vertical fault and horizontal slip on a horizontal fault, neither of which produce a Bouguer anomaly. The Bouguer anomaly produced by an open horizontal crack is the same as would be produced had the material within the crack been mined out without deforming the solid. If the horizontal crack were filled with material of density equal to the host rock (a good approximation to sill formation), no Bouguer anomaly is produced. For cracks of other inclinations the ratio of gravity change to uplift is not constant. Thus, dilatancy, in general, does not correspond to the absence of a free air anomaly, as might be suggested by the special case of a spherical source of dilatation. For two-dimensional models a cylindrical source of dilatation produces no free air gravity anomaly, dip-slip faulting produces no Bouguer anomaly, and open cracks produce a Bouguer anomaly equal to that which would be produced had the material within the crack been mined out without deforming the solid. A two-dimensional crack filled with material of density equal to that of the host rock would produce no Bouguer anomaly. Jachens et al. (1983) have reported temporal changes in gravity, elevation, and area strain along the San Andreas fault in southern California such that the Bouguer anomaly apparently remains unchanged and the uplift-to-strain ratio is about −100 km. Several dislocation mechanisms are proposed that fulfill these constraints, but these mechanisms appear to be rather contrived and are not regarded as satisfactory explanations.

Journal of Geophysical Research Solid Earth↗

Regional deformation near Palmdale, California, 1973-1983 (USA)

The Tehachapi trilateration network spans the intersection of the San Andreas and Garlock faults in southern California in the “Big Bend” region of the San Andreas fault. Analysis of data from 1973–1983 shows strain differences between the northwest and southeast regions of the network and slip at depth on both faults. The Palmdale network, spanning the San Andreas fault entirely within the Tehachapi network, showed increases of about 1 μstrain in both east‐west and north‐south extension in late 1979. The Tehachapi strains also jumped at this time, but the magnitude of the increase was only about one third that of Palmdale. The principal strain rates for Tehachapi over the time interval 1973–1983 are μstrain/yr and μstrain/yr, with the 1 axis directed N76°E. Strains were also computed for two subregions. The principal strain rates for the southeast Tehachapi region are μstrain/yr and μStrain/yr, with the 1 axis directed N73°E. This result differs significantly from the principal strain rates at Palmdale ( μstrain/yr, μstrain/yr, with the 1 axis directed N71°E). For the northwest Tehachapi region, μstrain/yr and μstrain/yr, with the 1 axis directed N93°E.

California↗

A seismic refraction study of the Oregon Cascades

A 275‐km‐long reversed refraction profile in the Oregon Cascades, two shallow earthquakes of magnitude 5 in southern Washington, a shallow earthquake of magnitude 4.6 in northern California, and a previously published analysis of the Bouguer gravity field are used to develop a crustal P wave velocity model for the Oregon Cascades. Travel time analysis of the refraction profile indicates a crustal structure characterized by surface layers with P wave velocities that vary from 2.9 to 5.2 km/s and thicknesses that vary from 2.5 to 5.0 km, upper crustal velocities of 6.1 to 6.5 km/s between the depths of 3 and 29 km, lower crustal velocities near 7.0 km/s between the depths of 29 and 44 km, and a mantle reflector at a depth of 44 km. Comparison with synthetic seismograms supports this structure and shows that the lower crustal and Mono transitions can be better modeled by continuous velocity gradients than by first‐order discontinuities. Arrival times from two shallow earthquakes in southern Washington in 1981 across 14 telemetry stations in the Oregon Cascades show an apparent velocity of 7.62 km/s. On the basis of a previously published gravity analysis, the mean north‐south component of the dip to the crust‐mantle boundary is estimated at 1° down dip to the south. The apparent velocity from the earthquakes and the dip estimate from gravity indicate that the true Pn velocity is 7.70 km/s. Arrival times from a shallow earthquake in northern California in 1978 across telemetry stations in the Washington Cascades are consistent with this model and indicate that the upper mantle velocity reaches 8.20 km/s at a depth of approximately 100 km.

Oregon↗

Deformation of clinopyroxenite: Evidence for a transition in flow mechanisms and semibrittle behavior

A systematic suite of constant strain rate experiments was performed on a vacuum-dried, high-purity, fine-grained clinopyroxenite using NaCl and NaF as confining media in a Griggs-type piston-cylinder apparatus. The experiments were carried out over a range of temperatures from 400° to 1100°C, strain rates from 10 −3 to 10 −7 s −1 , and confining pressures from 170 to 1990 MPa. At T = 600°C and = 1.1 × 10 −5 s −1 , three modes of deformation occur with increasing confining pressure: (1) Macroscopic faulting associated with low strength and stress drops, (2) stable microfracturing and plastic deformation associated with pressure-dependent strength, and (3) plastic deformation (mechanical twinning and 〈001〉 slip) with high strengths which are insensitive to pressure variations. In experiments at P = 1500 MPa, within this high-pressure plastic mode, two regimes of flow are clearly defined. At low to intermediate temperatures and high strain rates, flow strengths are insensitive to changes in strain rate and temperature. Optical and transmission electron microscope observations indicate that plastic strain is accomplished by mechanical twinning on (100) and (001) and by {100}〈001〉 slip. In contrast, at high temperatures and low strain rates the flow stress is strongly dependent on temperature and strain rate. Specimens deformed in this regime show evidence of recovery, multiple slip, and recrystallization; and plastic strain is much more homogeneous. The flow data within each regime can be satisfactorily fit to thermally activated power laws. In the low-temperature regime n (the stress exponent) = 83 ± 16 and E * (the activation energy for flow) = 220 ± 40 kJ/mol. We believe that these parameters reflect flow dominated by the kinetics of dislocation glide associated with mechanical twinning and (100)〈001〉 slip. In the high-temperature regime, n = 5.3 ± 1.1 and E * = 380 ± 30 kJ/mol. These parameters describe creep by multiple slip accompanied by increased rates of diffusion and recovery.

Journal of Geophysical Research Solid Earth↗

Permeability of fault gouge under confining pressure and shear stress

The permeability of both clay-rich and non-clay gouges, as well as several pure clays, was studied as a function of confining pressures from 5 to 200 MPa and shear strain to 10. Permeability ranged over four orders of magnitude, from around 10 −22 to 10 −18 m 2 (1 darcy = 0.987 × 10 −12 m 2 ). The lowest values were characteristic of the montmorillonite-rich and finer grained non-clay gouges. Illite, kaolinite, and chlorite had intermediate permeabilities, while the highest values were typical of the serpentine and coarser grained non-clay gouges. Grain size was an important factor in determining permeability, particularly for the clay-rich samples. The coarse grained gouges were the most permeable and decreased in permeability after shearing. Conversely, the fine grained gouges had characteristically lower permeabilities that did not vary significantly after various amounts of shearing. The permeabilities of the non-clay samples were not significantly different than those of the clays. Therefore, comminuted rock flours can be equally as effective in reducing the flow of water as the characteristically low permeability clay gouges. The strengths of the samples were quite variable. The non-clay gouges were consistently the strongest, with yield points (beginning of nonelastic behavior) around 850 MPa, while montmorillonite had an anomalously low strength in relation to all the other gouges at 250 MPa. Strength of the saturated samples under drained (low pore pressure) conditions did not correlate with high or low permeability. However, the low permeabilities of these gouges could be a factor in the measured low shear stresses along fault regions if excess pore pressures were created as a result of shearing or compaction, and this pressure was unable to dissipate through a thick section of the material.

Journal of Geophysical Research Solid Earth↗

Origin of Hawaiian tholeiite: A metasomatic model

Two voluminous magma types generated in the mantle underlying the Pacific plate are mid‐ocean ridge tholeiite (MORB) erupted at the East Pacific Rise spreading center and Hawaiian tholeiite (HT) erupted above the Hawaiian hot spot or melting anomaly. MORB has low initial 87 Sr/ 86 Sr ratios and low amounts of all incompatible trace elements including rare earths; chondrite‐normalized patterns are depleted in light rare earths. HT, by contrast, has higher initial 87 Sr/ 86 Sr and higher amounts of incompatible trace elements; chondrite‐nor‐malized patterns are enriched in the middle and light rare earths. HT is generally poorer in CaO and Al 2 O 3 and much richer in total iron and TiO 2 compared with MORB having the same MgO content. Primary magma compositions for the two volcanic systems are calculated in Fe‐Mg equilibrium with residual olivine (Fo 92 ). MORB is generated by partial melting of a trace element depleted Iherzolite source leaving a residual assemblage dominated by olivine and orthopyroxene. The percentage of partial melting for a primary magma containing 15% MgO is calculated to be 35–42% in a source mantle having a heavy rare earth content of 3×chondrite and 33–35% MgO. HT, represented by Kilauea tholeiite, is generated by partial melting of a mixture of unmelted and residual mantle for MORB which has been modified by metasomatic addition of a nephelinitic fluid, amphibole, and minor amounts of apatite and Fe‐bearing phases such as sulfide and magnetite/ilmenite. This model yields a picritic magma in equilbrium with magnesian dunite at high (>40%) degrees of partial melting. The source also has 35% MgO before partial melting. Melting in both systems in polyvariant and not controlled by lower‐temperature invariant equilibria. The low‐velocity zone is considered to be the source of metasomatic fluids that are driven upward into the lowermost lithosphere in response to a thermal plume. Picritic primary magmas are produced by shear melting, localized in the zone of thinned and metasomatized lithosphere beneath Hawaii. Melt extraction is rapid and episodic at intervals of months to decades; magma is not stored in the mantle but passes upward to a plexus of storage reservoirs located 2–6 km beneath the surface of Kilauea. Kilauea primary magmas fractionate olivine during upward transport to reach bulk compositions of 13–14% MgO in storage. Different magma batches erupted to the surface, distinguished by different major and minor element compositons compared at similar MgO content, represent combinations of differing degrees of metasomatic enrichment, differing degrees of partial melting, and some effects of premelting mantle heterogeneity.

Hawaii↗

Introduction and digest to the Special Issue on Chemical Effects of Water on the Deformation and Strengths of Rocks

The important role of pore pressure in promoting such brittle processes as cataclasis, hydraulic fracturing, large‐scale faulting, and earthquakes within the crust is widely accepted in geology and geophysics [, 1957; , 1959; , 1958; ., 1963; , 1968; ., 1968; ., 1976; , 1973, 1980; , 1981]. Provided that fluid pressure is fully communicated with rock pore space, the effective normal stresses that control crack growth, macroscopic fracture, and friction are reduced by the magnitude of the fluid pressui'e. Beyond this physical effect of pore fluids, there are chemical effects of water on the strength of rocks that are also important in governing differential stresseg and flow in the continental crust. Some of these chemical effects of water on rock deformation have long been recognized.

Journal of Geophysical Research B: Solid Earth↗

Experimental deformation of topaz crystals: Possible embrittlement by intracrystalline water

Crystallographically oriented single‐crystal prisms of gem quality topaz (composition AlSiO (OH− F) where x = 0.04 ± 0.01) were deformed at a confining pressure of 1.50 GPa, a temperature of 800°C, and a strain rate of 2×10 s. Under nearly identical conditions, all crystals of anhydrous rock‐forming minerals that have been tested to date, such as olivine, quartz, feldspars, pyroxenes, and refractory oxides, deform plastically; in contrast, our topaz crystals failed by brittle fracture regardless of the orientation of the compression direction. No optical evidence for plastic deformation was detected. Another suite of experiments with compression perpendicular to the (001) cleavage at = 100°–950°C and a strain rate of 2×10 s displayed two regimes of behavior: (1) at >400°C, fracture strength was independent of temperature, and fracture occurred on one or two surfaces parallel to {103}; (2) at <400°C, the fracture strength increased rapidly with decreasing temperature, no macroscopic stress drop was observed, and many closely spaced conjugate fractures formed on (103) and (103). The anomalous brittleness of topaz compared to anhydrous silicate and oxide crystals indicates that intracrystalline “water” plays a role in the embrittlement. We suggest that water within the topaz crystals promotes fracture in ways similar to the mechanisms of slow crack growth aided by environmental moisture.

Journal of Geophysical Research B: Solid Earth↗

Modification of wave-cut and faulting-controlled landforms

From a casual observation that the form of degraded fault scarps resembles the error function, this investigation proceeds through an elementary diffusion equation representation of landform evolution to the application of the resulting equations to the modern topography of scarplike landforms. The morphologic observations can be analyzed either in the form of one or more cross-strike elevation profiles or in the form of the slope-offset plot, a point plot of maximum scarp slope versus scarp offset. Working with either or both of these data representations for nine geologic structures, which range in age from 3 to 400 ka B.P. and in offset from 1 to 50 m, we apply analytical solutions for the vertical initial value scarp, the vertical continuous offset scarp, and the finite slope, initial value scarp. The model calculations are intrinsically ambiguous, yielding as the final answer only the product κ t (in the case of the initial value problem) or the product κ A −1 (in the case of the repeated faulting problem); here t is the age of a single scarp-forming event, 2 A is the vertical slip rate, and κ is the “mass diffusivity.” A single profile across three sea cliffs along the Santa Cruz, California, coast is analyzed as three separate initial value problems. A reasonably constrained age for the sea cliff standing above the Highway 1 platform returns κ = 11 GKG (1 GKG = 1 m 2 /ka). With this κ, we can date the two older sea cliffs. In fact, we do the converse: age estimates for these two older sea cliffs based on a uniform rate of uplift both yield the same κ as for the lower sea cliff. We treat a single profile of the Raymond fault in Pasadena/San Marino in terms of the repeated faulting problem; for it the uplift rate of R. Crook and others yields κ = 16 GKG. The very substantial preexisting offset across the Raymond fault must have been buried/leveled some 230 ka B.P., when the modern topography began to form. Our analysis of the Lake Bonneville shoreline scarps reveals a dependence of κ t on 2a, suggestive of nonlinear modification processes. This appearance is treated with the finite slope initial value scarp model to determine κ=1.1 GKG for the Lake Bonneville shoreline scarps. The suggestion of M. N. Machette that approximately 100,000-year-old, meter-high scarps are “unobservable” in weakly consolidated alluvial terranes of the Basin and Range and Rio Grande Rift Valley provinces can be formulated as κ ≳ 1 GKG. The coincidence between this inequality and the Lake Bonneville shoreline κ is striking, and it suggests that the value of κ = 1 GKG may be generally applicable, as a good first approximation, to the modification of alluvial terranes within the semiarid regions of the western United States. The Lake Bonneville shoreline κ is the basis for dating four sets of fault scarps in west-central Utah. The Drum Mountains fault scarps can be modeled in several different circumstances, but the most likely interpretation is that these fault scarps formed as the result of a single episode of normal faulting 3.6 to 5.7 ka B.P. The younger age is associated with quite low initial slope angles (25°). The other three sets of fault scarps show no evidence for finite initial value slopes. Fault scarps along the eastern base of the Fish Springs Range are very young, 3 ka B.P. We estimate the age of fault scarps along the western flank of the Oquirrh Mountains to be 32 ka B.P., which meets the weak geologic constraint that they be older than the Lake Bonneville shoreline. Fault scarps along the northeastern margin of the Sheeprock Mountains are even older, 53 ka B.P. An intriguing consequence of our single-event analysis of these scarps is that an 11.5-m offset occurred in a single earthquake.

Journal of Geophysical Research Solid Earth↗

Deformation in the White Mountain seismic gap, California-Nevada, 1972-1982

A 100×40 km trilateration network extending from Bishop, California, to near Hawthorne, Nevada, crosses the east end of the Long Valley caldera, site of renewed magma inflation in the 1979–1980 interval, and spans most of the White Mountain seismic gap. The network was surveyed in 1972, 1973, 1976, 1979, 1980, and 1982. The 1980 survey may be contaminated by a scale error. In addition, leveling surveys across the caldera have been run in 1932, 1957, 1975, 1980, 1982, and 1983. Interpretation of the deformation is complicated by the occurrence of the May 1980 Mammoth Lakes earthquake sequence (four earthquakes M L ≥6) at the south edge of the caldera as well as other moderate earthquakes within the White Mountain seismic gap. The vertical deformation is largely accounted for by 0.10‐ to 0.15‐km 3 expansion of a spherical magma chamber 8–10 km beneath the resurgent dome within the Long Valley caldera sometime between July 1979 and September 1980 with an additional expansion of perhaps 0.05 km 3 between September 1980 and July 1982. Some additional sources of deformation within the aftershock zone of the Mammoth Lakes earthquakes seem to be required to explain the horizontal deformation. We show that right‐lateral slip on vertical faults extending WNW from each of the three largest earthquakes in the Mammoth Lakes sequence provides the required additional deformation, but this solution is by no means unique. There are simply too few data to define the rather complex deformation that apparently occurred within the aftershock zone. There is little doubt, however, that inflation of a magma chamber beneath the resurgent dome within the Long Valley caldera was involved in the deformation.

California↗

Earthquake swarm in Long Valley caldera, California, January 1983: Evidence for dike inflation

The 1982–1983 deformation observed by trilateration and leveling surveys across the Long Valley caldera is apparently related to the 8.5‐km‐long by 8‐km‐deep vertical rupture surface defined by the January 1983 earthquake swarm that occurred in the south moat of the caldera. The observed deformation can be explained as follows. In late 1982, 0.03 km 3 of magma was injected into a dike that dips 30° northward from the bottom of the rupture surface. The downdip dimension of this dike is 8 km. The dike inflation accounts for the uplift observed across the caldera as well as some of the horizontal deformation. Inflation of the dike generated a tension of about 3 bars across the vertical plane that was to become the rupture surface of the January swarm. This reduced the frictional stress on the rupture plane and perhaps triggered the slip that caused the January swarm. Right‐lateral slip (0.22 m) on the uppermost 2 km of the rupture plane during and after the January swarm accounts for the additional horizontal deformation observed. The model can be improved marginally if strike slip is admitted over the entire rupture surface and 0.006 km 3 of magma is injected along that surface in the depth interval 3–8 km. The improvement in the model fit, however, is not sufficient to require shallow injection of magma. Thus we conclude that inflation of a dike at depth (8–12 km) dipping northward beneath the resurgent dome plus shallow right‐lateral slip on the rupture surface is a simple, but not unique, explanation of the observed deformation and seismicity.

California↗

Catastrophic isotopic modification of rhyolitic magma at times of caldera subsidence, Yellowstone Plateau Volcanic Field

The Yellowstone Plateau volcanic field has undergone repeated eruption of rhyolitic magma strongly depleted in 18 O. Large calderas subsided 2.0, 1.3, and 0.6 Ma ago, on eruption of ash flow sheets that represent at least 2500, 280, and 1000 km 3 of zoned magma. More than 60 other rhyolite lavas and tuffs permit reconstruction of the long-term chemical and isotopic evolution of the silicic system. Narrow δ 18 O ranges in the ash flow sheets contrast with wide δ 18 O variations in postcaldera lavas of the first and third caldera cycles. Earliest postcollapse lavas are 3 to 6‰ lighter than the preceding ash flow sheets. The O 18 depletions were short-lived events that immediately followed caldera subsidence; hundreds of cubic kilometers of magma were drastically 18 O depleted and thousands were depleted by 1–2‰. Sequences of postcaldera lavas record partial recovery toward precaldera δ 18 O values; secular trends between collapse events thus reflect gradual reenrichment of the roofmost magma in δ 18 O. Much of the subcaldera reservoir was affected, because lavas that erupted as far apart as 115 km reflect the same pattern of depletion and partial recovery. Contemporaneous extracaldera rhyolites have the highest δ 18 O values in the volcanic field and show no effects of the repeated depletions. Sr and Pb isotope ratios of intracaldera rhyolites jump to more radiogenic values at times of caldera formation and show a longterm zigzag pattern like that of δ 18 O. Although some contamination by foundering roof rocks seenis to be required, water was probably the predominant contaminant. Even if roof rocks had been strongly depleted in O 18 before engulfment, their assimilation would have been far from sufficient to account for the large O 18 shift. The low- O 18 lavas contain no xenocrysts and show no trace element or phenocryst evidence of massive contamination. Their Fe-Ti-oxide temperatures indicate no cooling relative to the caldera-forming ash flow magma, and their whole-rock, glass, and phenoeryst chemistry suggests compositional continuity with the ash flow sequence. Oxygen exchange between the magma and a mass of low-O 18 water greatly exceeding solubility limits may require (1) recurrent explosive activity to sustain access and mixing of water with the magma and (2) convection of the magma reservoir to prevent local saturation.

Journal of Geophysical Research Solid Earth↗

Modeling crater topography and albedo from monoscopic Viking orbiter images: 1. Methodology

A new photoclinometric technique for extraction of topographic data from single planetary images is presented that overcomes many previous limitations of photoclinometry. The procedure fully compensates for oblique viewing geometry prevalent in spacecraft images. Albedo variations have been one of the most serious obstacles in the application of photoclinometry to planetary surfaces. This problem is overcome in the topographic solution by simultaneously utilizing brightness data from a pair of profiles; both segments are assumed to have the same topographic and albedo variations along their lengths. Profile directions are chosen where the orientation of downslope or upslope is obvious, thus resolving a major ambiguity in photoclinometry. This requirement is particularly easy to satisfy for craters and not very difficult for many irregular features. An additional procedure is presented that eliminates even the requirement of topographic symmetry along the pair of profiles. If two profiles have the same relief but their shapes are very different, another method can be used in an iterative process to derive topographic profiles; however, this procedure does, require that the albedo not vary along the profiles. Test results indicate that both procedures have an accuracy and precision of approximately 2° for slopes of typical bowl-shaped craters, which translates to approximately 5% for depths.

Journal of Geophysical Research Solid Earth↗

Propagation and linkage of oceanic ridge segments

We have investigated the propagation of spreading ridges and the development of structures that link ridge segments using an analogy between ridges and cracks in elastic plates. The ridge-propagation force and a path factor that controls propagation direction were calculated for echelon ridge segments propagating toward each other. The ridge-propagation force increases as ridge ends approach but then declines sharply as the ends pass, so ridge segments may overlap somewhat. The sign of the path factor changes as ridge ends approach and pass, so the overlapping ridge ends may diverge and then converge following a hook-shaped path. The magnitudes of shear stresses in the plane of the plate and orientations of maximum shear planes between adjacent ridge segments were calculated to study transform faulting. For different loading conditions simulating ridge push, plate pull, and ridge suction, we identify a zone of intense mechanical interaction between adjacent ridge ends in which stresses are concentrated. For all conditions, the shear stress in the interaction zone increases as ends approach and remains large as the ends overlap; thus crust in this zone may fracture and weaken in preparation for the formation of a through-going transform fault. The calculated shear planes rotate toward an orientation about 90° from the strike of ridges as the ends pass, thus favoring the orthogonal arrangement of ridges and transforms. The magnitudes of mean stresses in the plane of the plate and orientations of principal stress planes were also calculated. The mean stress is tensile in the interaction zone, so basins may form there, except in the case of ridge push loading. The planes across which the maximum tension acts are oblique to ridges, thus favoring obliquely oriented normal faults bounding the transform valley.

Journal of Geophysical Research Solid Earth↗

Interpretation of gravity data in a complex volcano-tectonic setting, southwestern Nevada

This regional gravity study, based on an irregular 2-km data grid, was conducted during the past few years at Yucca Mountain, southern Nye County, Nevada, as part of a program to locate a suitable repository for high-level nuclear waste. About 100 surface rock samples, three borehole gamma-gamma logs, and one borehole gravity study provide excellent density control. A nearly linear increase in density of 0.26 g/cm 3 per kilometer of depth is indicated in the thick tuff sequences that underlie the mountain. Isostatic and 2.0-g/cm 3 Bouguer corrections were applied to the observed gravity values to remove regional gradients and topographic effects, respectively. The Bare Mountain gravity high, with an isostatic anomaly maximum of 48 mGal, is connected with a greater gravity high over the Funeral Mountains, to the southwest; together, these highs result from a continuous block of dense, metamorphosed Precambrian and Paleozoic rocks that stretches across much of the Walker Lane from the east edge of Death Valley to Bare Mountain. The Calico Hills gravity high appears more likely to originate from a northeast trending buried ridge of Paleozoic rocks that extends southwestward beneath Busted Butte, 5 km southeast of the proposed repository, where two- and three-dimensional modeling indicates that the pre-Cenozoic rocks lie less than 1000 m beneath the surface. Tuff, at least 4000 m thick, fills a large steep-sided depression in the pretuff rocks beneath Yucca Mountain and Crater Flat. The gravity low and the thick tuff section lie within a large collapse area that includes the Crater Flat-Timber Mountain-Silent Canyon caldera complexes. Gravity lows in Crater Flat itself are interpreted to coincide with the source areas of the Prow Pass Member, the Bullfrog Member, and the Tram Member of the Crater Flat Tuff; these source areas add nearly 350 km 2 to the previously recognized extent of the local caldera complexes. Southward extension of the broad gravity low associated with Crater Flat into the Amargosa Desert is evidence for sector graben-type collapse segments related to the formation of the Timber Mountain caldera and superimposed on the other volcanic and extensional structures within Crater Flat.

Journal of Geophysical Research Solid Earth↗

Ice and debris in the fretted terrain, Mars

Viking moderate- and high-resolution images along the northern highland margin were studied monoscopically and stereoscopically to contribute to an understanding of the development of fretted terrain. Results support the hypothesis that the fretting process involved flow facilitated by interstitial ice. The process apparently continued for a long period of time, and debris-apron formation shaped the fretted terrain in the past as well as the present. Interstitial ice in debris aprons is most likely derived from ground ice obtained by sapping or scarp collapse. Debris aprons could have been removed by sublimation if they consisted mostly of ice, or by deflation if they consisted mostly of debris. To remove the debris, wind erosion was either very intense early in martian history, or was intermittent, perhaps owing to climatic cycles.

Journal of Geophysical Research↗

Lu-Hf constraints on the evolution of lunar basalts

Very low Ti basalts and green glass samples from the moon show high Lu/Hf ratios and low Hf concentrations. Low-Ti lunar basalts show high and variable Lu/Hf ratios and higher Hf concentrations, whereas high-Ti lunar basalts show low Lu/Hf ratios and high Hf concentrations. KREEP basalts have constant Lu/Hf ratios and high but variable Hf concentrations. Using the Lu-Hf behavior as a constraint, we propose a model for the mare basalts evolution. This constraint requires extensive crystallization of the primary lunar magma ocean prior to formation of the lunar mare basalt sources and the KREEP basalts. Mare basalts are produced by the melting of the cumulate rocks, and KREEP basalts represent the residual liquid of the magma ocean. Lu and Hf concentrations and the Hf isotopic data of lunar rocks suggest that assimilation cannot be accepted as a major process to explain the diversity of the lunar mare basalts. The urKREEP hypothesis is also unnecessary. Both high- and low-Ti basalts show enough iron enrichment to be regarded as melting products of the last stage cumulate rocks from the lunar magma ocean. The KREEP basalts are also rich in iron and may be regarded as the final, residual liquid left after the crystallization of the major portion of the primary lunar magma ocean.

Journal of Geophysical Research Solid Earth↗