Search USGS⌕ Search

SEARCH · Search USGS

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 361 records · Page 20Linked to original sources
USGS2000

Pore fluid pressure, apparent friction, and Coulomb failure

Many recent studies of stress-triggered seismicity rely on a fault failure model with a single free parameter, the apparent coefficient of friction, presumed to be a material constant with possible values 0 ≤ μ′ ≤ 1. These studies may present a misleading view of fault strength and the role of pore fluid pressure in earthquake failure. The parameter μ′ is intended to incorporate the effects of both friction and pore pressure, but is a material constant only if changes in pore fluid pressure induced by changes in stress are proportional to the normal stress change across the potential failure plane. Although specific models of fault zones permit such a relation, neither is it known that fault zones within the Earth behave this way, nor is this behavior expected in all cases. In contrast, for an isotropic homogeneous poroelastic model the pore pressure changes are proportional to changes in mean stress, μ′ is not a material constant, and −∞ ≤ μ′ ≤ +∞. Analysis of the change in Coulomb failure stress for tectonically loaded reverse and strike-slip faults shows considerable differences between these two pore pressure models, suggesting that such models might be distinguished from one another using observations of triggered seismicity (e.g., aftershocks). We conclude that using the constant apparent friction model exclusively in studies of Coulomb failure stress is unwise and could lead to significant errors in estimated stress change and seismic hazard.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Seismic mapping of shallow fault zones in the San Gabriel Mountains from the Los Angeles Region Seismic Experiment, southern California

During the Los Angeles Region Seismic Experiment (LARSE), a reflection/refraction survey was conducted along a profile (line 1) extending from Seal Beach, California, northeastward to the Mojave Desert and crossing the Los Angeles and San Gabriel Valley basins and San Gabriel Mountains. In most shot gathers from the southern and central San Gabriel Mountains, clear secondary arrivals are seen that merge, or appear to merge, with first arrivals at three locations, including the location of the Vincent thrust fault, an exposed late Mesozoic/early Cenozoic megathrust. These secondary arrivals are interpretable as reflections in the shallow crust (<5 km depth) from a concave-upward interface that projects to the surface in the north near the Vincent thrust fault, is offset in its central part at the San Gabriel fault (an old branch of the San Andreas fault), and terminates in the south at 1 to 2 km depth at the southern mountain front. The velocity structure above and below this interface strongly suggests it is the Vincent thrust fault: intermediate velocities (6.2 km/s), consistent with mylonites overlying the Vincent thrust fault, are observed above it; lower velocities (5.8 km/s), consistent with the Pelona Schist underlying the Vincent thrust fault, are observed below it. Problems arise, however, in attempting to match this reflector to the exposed Vincent thrust fault, which is seen in outcrops east of line 1. The Vincent thrust fault is shallower than the reflector in most places. An unmapped structure (steep fault, monocline, or thrust fault) is required between line 1 and the outcrops that either drops the Vincent thrust fault down to the depths of the reflector or repeats the Vincent thrust fault beneath line 1 in the footwall of another thrust fault. An alternative interpretation of the reflector is a deep greenstone horizon within the Pelona Schist, although this alternative is not favored by the velocity structure. Copyright 2001 by the American Geophysical Union.

California↗
USGS2000

Fracture process zone in granite

In uniaxial compression tests performed on Aue granite cores (diameter 50 mm, length 100 mm), a steel loading plate was used to induce the formation of a discrete shear fracture. A zone of distributed microcracks surrounds the tip of the propagating fracture. This process zone is imaged by locating acoustic emission events using 12 piezoceramic sensors attached to the samples. Propagation velocity of the process zone is varied by using the rate of acoustic emissions to control the applied axial force. The resulting velocities range from 2 mm/s in displacement-controlled tests to 2 μm/s in tests controlled by acoustic emission rate. Wave velocities and amplitudes are monitored during fault formation. P waves transmitted through the approaching process zone show a drop in amplitude of 26 dB, and ultrasonic velocities are reduced by 10%. The width of the process zone is ∼9 times the grain diameter inferred from acoustic data but is only 2 times the grain size from optical crack inspection. The process zone of fast propagating fractures is wider than for slow ones. The density of microcracks and acoustic emissions increases approaching the main fracture. Shear displacement scales linearly with fracture length. Fault plane solutions from acoustic events show similar orientation of nodal planes on both sides of the shear fracture. The ratio of the process zone width to the fault length in Aue granite ranges from 0.01 to 0.1 inferred from crack data and acoustic emissions, respectively. The fracture surface energy is estimated from microstructure analysis to be ∼2 J. A lower bound estimate for the energy dissipated by acoustic events is 0.1 J.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Earthquake stress drop and laboratory-inferred interseismic strength recovery

We determine the scaling relationships between earthquake stress drop and recurrence interval t r that are implied by laboratory-measured fault strength. We assume that repeating earthquakes can be simulated by stick-slip sliding using a spring and slider block model. Simulations with static/kinetic strength, time-dependent strength, and rate- and state-variable-dependent strength indicate that the relationship between loading velocity and recurrence interval can be adequately described by the power law V L ∝ t r n where n≈−1. Deviations from n=−1 arise from second order effects on strength, with n>−1 corresponding to apparent time-dependent strengthening and n<−1 corresponding to weakening. Simulations with rate and state-variable equations show that dynamic shear stress drop Δτ d scales with recurrence as d Δτ d / d ln t r ≤σ e ( b-a ), where σ e is the effective normal stress, μ=τ/σ e , and ( a-b )= d μ ss / d ln V is the steady-state slip rate dependence of strength. In addition, accounting for seismic energy radiation, we suggest that the static shear stress drop Δτ s scales as d Δτ s / d ln t r ≤σ e (1 +ζ)( b-a ), where ζ is the fractional overshoot. The variation of Δτ s with ln t r for earthquake stress drops is somewhat larger than implied by room temperature laboratory values of ζ and b-a . However, the uncertainty associated with the seismic data is large and the discrepancy between the seismic observations and the rate of strengthening predicted by room temperature experiments is less than an order of magnitude.

Journal of Geophysical Research B: Solid Earth↗
USGS2000

Acoustic properties of a crack containing magmatic or hydrothermal fluids

We estimate the acoustic properties of a crack containing magmatic or hydrothermal fluids to quantify the source properties of long-period (LP) events observed in volcanic areas assuming that a crack-like structure is the source of LP events. The tails of synthetic waveforms obtained from a model of a fluid-driven crack are analyzed by the Sompi method to determine the complex frequencies of one of the modes of crack resonance over a wide range of the model parameters α/ a and ρ f /ρ s , where αis the P wave velocity of the rock matrix, a is the sound speed of the fluid, and ρ f and p s are the densities of the fluid and rock matrix, respectively. The quality factor due to radiation loss ( Q r ) for the selected mode almost monotonically increases with increasing α/ a , while the dimensionless frequency ( v ) of the mode decreases with increasing α/ a and ρ f /ρ s . These results are used to estimate Q and v for a crack containing various types of fluids (gas-gas mixtures, liquid-gas mixtures, and dusty and misty gases) for values of a , ρ f , and quality factor due to intrinsic losses ( Q i ) appropriate for these types of fluids, in which Q is given by Q −l = Q r −l + Q i −1 . For a crack containing such fluids, we obtain Q ranging from almost unity to several hundred, which consistently explains the wide variety of quality factors measured in LP events observed at various volcanoes. We underscore the importance of dusty and misty gases containing small-size particles with radii around 1 μm to explain long-lasting oscillations with Q significantly larger than 100. Our results may provide a basis for the interpretation of spatial and temporal variations in the observed complex frequencies of LP events in terms of fluid compositions beneath volcanoes.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Modeling consolidation and dewatering near the toe of the northern Barbados accretionary complex

At the toe of the northern Barbados accretionary complex, temperature and pore water chemistry data indicate that fluid flow is channeled along the décollement and other shallow thrust faults. We examine mechanisms that may prevent consolidation and maintain high permeability over large sections of the décollement. High-resolution bulk density data from five boreholes show that the décollement is well consolidated at some sites while other sites remain underconsolidated. Underconsolidated décollement behavior is associated with kilometer-scale negative-polarity seismic reflections from the décollement plane that have been interpreted to be fluid conduits. We use a coupled fluid flow/consolidation model to simulate the loading response of a 10-km-long by 680-m-thick slice of sediment as it enters the accretionary complex. The simulations capture 185 ka (5 km) of subduction, with a load function representing the estimated effective stress of the overriding accretionary prism (3.8° taper angle). Simulation results of bulk density in the décollement 3.2 km arcward of the deformation front are compared with observations. The results show that persistent high pore pressures at the arcward edge of the simulation domain can explain underconsolidated behavior. The scenario is consistent with previous modeling results showing that high pore pressures can propagate intermittently along the décollement from deeper in the complex. Simulated seaward fluxes in the décollement (1–14 cm yr −1 ) lie between previous estimates from modeling studies of steady state (<1 cm yr −1 ) and transient (>1 m yr −1 ) flow. Maximum simulated instantaneous fluid sources (2.5×10 −13 s −1 ) are comparable to previous estimates. The simulations show minor swelling of incoming sediments (fluid sources ∼−3×10 −15 s −1 ) up to 3 km before subduction that may help to explain small-scale shearing and normal faulting proximal to the protodécollement.

Atlantic Ocean, Northern Barbados Accretionary Com↗
USGS2001

Travel times of P and S from the global digital seismic networks: Implications for the relative variation of P and S velocity in the mantle

We present new data sets of P and S arrival times which have been handpicked from long-period vertical and transverse component recordings of the various global seismic networks. Using events which occurred from 1976 to 1994 results in ∼38,000 globally well-distributed measurements of teleseismic P and ∼41,000 measurements of S . These data are particularly useful for looking at the relative variation of S and P velocities in the lower mantle. We describe both the measurement techniques and the gross characteristics of the data sets. The size of our data sets allows us to exploit the internal consistency of the data to identify outliers using a summary ray analysis. Since the polarity of each arrival is also known, we can construct fault plane solutions and/or compare with polarities predicted by the Harvard centroid moment tensor solutions to further diagnose phase misidentification. This analysis results in ∼5% of the data being identified as outliers. An analysis of variance indicates that the S residual travel times are dominated by the effects of three-dimensional structure but the P data have comparable contributions from noise and source mislocation effects. The summary ray analysis reveals the basic character of lower mantle structure, and there are large-scale patterns in both the S and P data sets that correlate quite well with each other. This analysis suggests that on average, d ln v s / d ln v p is an increasing function of depth in the mantle going from a value of ∼1.7 at the top of the lower mantle to an apparent value of 4 near the base of the mantle. This latter extreme value of R seems to result mainly from data which sample one region in the lowermost mantle under the central Pacific, where large positive S residuals are associated with very small P residuals. Such an anomaly cannot be thermal in origin.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Deformation across the Pacific-North America plate boundary near San Francisco, California

We have detected a narrow zone of compression between the Coast Ranges and the Great Valley, and we have estimated slip rates for the San Andreas, Rodgers Creek, and Green Valley faults just north of San Francisco. These results are based on an analysis of campaign and continuous Global Positioning System (GPS) data collected between 1992 and 2000 in central California. The zone of compression between the Coast Ranges and the Great Valley is 25 km wide. The observations clearly show 3.8±1.5 mm yr −1 of shortening over this narrow zone. The strike slip components are best fit by a model with 20.8±1.9 mm yr −1 slip on the San Andreas fault, 10.3±2.6 mm yr −1 on the Rodgers Creek fault, and 8.1±2.1 mm yr −1 on the Green Valley fault. The Pacific-Sierra Nevada-Great Valley motion totals 39.2±3.8 mm yr −1 across a zone that is 120 km wide (at the latitude of San Francisco). Standard deviations are one σ. The geodetic results suggest a higher than geologic rate for the Green Valley fault. The geodetic results also suggest an inconsistency between geologic estimates of the San Andreas rate and seismologic estimates of the depth of locking on the San Andreas fault. The only convergence observed is in the narrow zone along the border between the Great Valley and the Coast Ranges.

California↗
USGS2001

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

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

Thermal thickness and evolution of Precambrian lithosphere: A global study

The thermal thickness of Precambrian lithosphere is modeled and compared with estimates from seismic tomography and xenolith data. We use the steady state thermal conductivity equation with the same geothermal constraints for all of the Precambrian cratons (except Antarctica) to calculate the temperature distribution in the stable continental lithosphere. The modeling is based on the global compilation of heat flow data by Pollack et al. [1993] and more recent data. The depth distribution of heat‐producing elements is estimated using regional models for ∼300 blocks with sizes varying from 1°×1° to about 5°×5° in latitude and longitude and is constrained by laboratory, seismic and petrologic data and, where applicable, empirical heat flow/heat production relationships. Maps of the lateral temperature distribution at depths 50, 100, and 150 km are presented for all continents except Antarctica. The thermal thickness of the lithosphere is calculated assuming a conductive layer overlying the mantle with an adiabat of 1300°C. The Archean and early Proterozoic lithosphere is found to have two typical thicknesses, 200–220 km and 300–350 km. In general, thin (∼220 km) roots are found for Archean and early Proterozoic cratons in the Southern Hemisphere (South Africa, Western Australia, South America, and India) and thicker (>300 km) roots are found in the Northern Hemisphere (Baltic Shield, Siberian Platform, West Africa, and possibly the Canadian Shield). We find that the thickness of continental lithosphere generally decreases with age from >200 km beneath Archean cratons to intermediate values of 200±50 km in early Proterozoic lithosphere, to about 140±50 km in middle and late Proterozoic cratons. Using known crustal thickness, our calculated geotherms, and assuming that isostatic balance is achieved at the base of the lithosphere, we find that Archean and early Proterozoic mantle lithosphere is 1.5% less dense (chemically depleted) than the underlying asthenosphere, while middle and late Proterozoic subcrustal lithosphere should be depleted by ∼0.6–0.7%. Our results suggest three contrasting stages of lithosphere formation at the following ages: >2.5 Ga, 2.5–1.8 Ga, and <1.8 Ga. Ages of komatiites, greenstone belts, and giant dike swarms broadly define similar stages and apparently reflect secular changes in mantle temperature and, possibly, convection patterns.

Journal of Geophysical Research B: Solid Earth↗
USGS2001

Creation of residual flows in a partially stratified estuary

The creation of residual flows in estuaries is examined using acoustic Doppler current profiler data sets from northern San Francisco Bay. The data sets are analyzed using principal component analysis to examine the temporal variability of the flows which create the residual circulation. It is seen that in this periodically and partially stratified estuary the residual flows are created through a series of pulses with strong variability at the 24‐hour timescale, through the interaction of shear, stratification and mixing. This interaction is captured through the use of a dimensionless number, the horizontal Richardson number ( Ri x ), which is developed to examine the local balance between the stratifying and destratifying forces at the tidal timescale. It is seen that Ri x is a valuable parameter in predicting the onset of the residual‐creating events, with a threshold value of ≈3 on ebb tides. This critical value is argued to be a threshold, above which the stratification and shear flow create a feedback effect, each further intensifying the other. This feedback results in a highly variable exchange flow which creates the estuarine residual in intermittent pulses rather than as a steady flow. Although typically attributed to baroclinic forcing, an argument is made that these pulses of residual‐creating exchange flow could be created by barotropic forcing in the presence of variable stratification which is asymmetric between flood and ebb tides. This result poses a great challenge for turbulence modeling, as the timing and magnitude of stratification and shear must be correctly simulated on the tidal timescale in order to reproduce the effects seen in the data sets presented.

Journal of Geophysical Research C: Oceans↗
USGS2002

Across-shelf sediment transport: Interactions between suspended and bed sediment

We use a two-dimensional, time-dependent sediment-transport model to quantify across-shelf transport, deposition, and sorting during wave-driven resuspension events characteristic of those that dominate sediment transport on many continental shelves. Decreases in wave-orbital velocities as water depth increases, and the resulting cross-shelf gradient in bed shear stress favor a net offshore transport of sediment. On wide, flat shelves (slopes ∼0.1%percnt;), these gradients are low, and the depth to which the seabed is reworked depends mainly on bottom shear stress and local sediment availability. On narrow, steep shelves (slopes ∼0.5%percnt;), however, the gradient in bottom stress generates significant cross-shelf suspended sediment flux gradients that create regions of net erosion and deposition. While the magnitude of waves generally determines the water depth to which sediment can be resuspended, erosional and depositional patterns on narrow shelves are sensitive to cross-shelf gradients in wave energy, nonlocal sediment availability, and the direction and magnitude of the cross-shelf current. During energetic waves, cross-shelf divergence of suspended sediment flux can create a coarsened, erosional area on the inner shelf that abuts a region of fine-grained sediment deposition on the mid-to-outer shelf. If currents are strongly shoreward, however, flux divergence leads to erosion over the entire shelf.

Journal of Geophysical Research C: Oceans↗
USGS2000

Numerical model of frazil ice and suspended sediment concentrations and formation of sediment laden ice in the Kara Sea

A one-dimensional (vertical) numerical model of currents, mixing, frazil ice concentration, and suspended sediment concentration has been developed and applied in the shallow southeastern Kara Sea. The objective of the calculations is to determine whether conditions suitable for turbid ice formation can occur during times of rapid cooling and wind- and wave-induced sediment resuspension. Although the model uses a simplistic approach to ice particles and neglects ice-sediment interactions, the results for low-stratification, shallow (∼20-m) freeze-up conditions indicate that the coconcentrations of frazil ice and suspended sediment in the water column are similar to observed concentrations of sediment in turbid ice. This suggests that wave-induced sediment resuspension is a viable mechanism for turbid ice formation, and enrichment mechanisms proposed to explain the high concentrations of sediment in turbid ice relative to sediment concentrations in underlying water may not be necessary in energetic conditions. However, salinity stratification found near the Ob' and Yenisey Rivers damps mixing between ice-laden surface water and sediment-laden bottom water and probably limits incorporation of resuspended sediment into turbid ice until prolonged or repeated wind events mix away the stratification. Sensitivity analyses indicate that shallow (≤20 m), unstratified waters with fine bottom sediment (settling speeds of ∼1 mm s −1 or less) and long open water fetches (>25 km) are ideal conditions for resuspension.

Kara Sea↗
USGS2001

Primary centers and secondary concentrations of tectonic activity through time in the western hemisphere of Mars

Five main stages of radial and concentric structures formed around Tharsis from the Noachian through the Amazonian as determined by geologic mapping of 24,452 structures within the stratigraphic framework of Mars and by testing their radial and concentric orientations. Tectonic activity peaked in the Noachian (stage 1) around the largest center, Claritas, an elongate center extending more than 20° in latitude and defined by about half of the total grabens which are concentrated in the Syria Planum, Thaumasia, and Tempe Terra regions. During the Late Noachian and Early Hesperian (stage 2), extensional structures formed along the length of present-day Valles Marineris and in Thaumasia (with a secondary concentration near Warrego Vallis) radial to a region just to the south of the central margin of Valles Marineris. Early Hesperian (stage 3) radial grabens in Pavonis, Syria, Ulysses, and Tempe Terra and somewhat concentric wrinkle ridges in Lunae and Solis Plana and in Thaumasia, Sirenum, Memnonia, and Amazonis are centered northwest of Syria with secondary centers at Thaumasia, Tempe Terra, Ulysses Fossae, and western Valles Marineris. Late Hesperian/Early Amazonian (stage 4) structures around Alba Patera, the northeast trending alignment of Tharsis Montes, and Olympus Mons appears centered on Alba Patera. Stage 5 structures (Middle-Late Amazonian) represent the last pulse of Tharsis-related activity and are found around the large shield volcanoes and are centered near Pavonis Mons. Tectonic activity around Tharsis began in the Noachian and generally decreased through geologic time to the Amazonian. Statistically significant radial distributions of structures formed during each stage, centered at different locations within the higher elevations of Tharsis. Secondary centers of radial structures during many of the stages appear related to previously identified local magmatic centers that formed at different times and locations throughout Tharsis.

Journal of Geophysical Research E: Planets↗
USGS2001

Mars' "White Rock" feature lacks evidence of an aqueous origin: Results from Mars Global Surveyor

The “White Rock” feature on Mars has long been viewed as a type example for a Martian playa largely because of its apparent high albedo along with its location in a topographic basin (a crater). Data from the Mars Global Surveyor Thermal Emission Spectrometer (TES) demonstrate that White Rock is not anomalously bright relative to other Martian bright regions, reducing the significance of its albedo and weakening the analogy to terrestrial playas. Its thermal inertia value indicates that it is not mantled by a layer of loose dust, nor is it bedrock. The thermal infrared spectrum of White Rock shows no obvious features of carbonates or sulfates and is, in fact, spectrally flat. Images from the Mars Orbiter Camera show that the White Rock massifs are consolidated enough to retain slopes and allow the passage of saltating grains over their surfaces. Material appears to be shed from the massifs and is concentrated at the crests of nearby bedforms. One explanation for these observations is that White Rock is an eroded accumulation of compacted or weakly cemented aeolian sediment.

Journal of Geophysical Research E: Planets↗
USGS2001

Mountains on Io: High-resolution Galileo observations, initial interpretations, and formation models

During three close flybys in late 1999 and early 2000 the Galileo spacecraft acquired new observations of the mountains that tower above Io's surface. These images have revealed surprising variety in the mountains' morphologies. They range from jagged peaks several kilometers high to lower, rounded structures. Some are very smooth, others are covered by numerous parallel ridges. Many mountains have margins that are collapsing outward in large landslides or series of slump blocks, but a few have steep, scalloped scarps. From these observations we can gain insight into the structure and material properties of Io's crust as well as into the erosional processes acting on Io. We have also investigated formation mechanisms proposed for these structures using finite-element analysis. Mountain formation might be initiated by global compression due to the high rate of global subsidence associated with Io's high resurfacing rate; however, our models demonstrate that this hypothesis lacks a mechanism for isolating the mountains. The large fraction (∼40%) of mountains that are associated with paterae suggests that in some cases these features are tectonically related. Therefore we have also simulated the stresses induced in Io's crust by a combination of a thermal upwelling in the mantle with global lithospheric compression and have shown that this can focus compressional stresses. If this mechanism is responsible for some of Io's mountains, it could also explain the common association of mountains with paterae.

Journal of Geophysical Research E: Planets↗