Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Oceans”

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

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

At least 415 records · Page 23Linked to original sources
USGS2002

Undrained poroelastic response of sandstones to deviatoric stress change

Deformation of porous crustal rock through diagenesis, tectonic loading, or other processes can change pore volume and affect fluid pressure. The largest stress-induced pore pressure changes occur when fluid is trapped in pores in an “undrained” condition. We have measured the undrained poroelastic response of two sandstones to changes in mean and deviatoric stress. Pore pressure was found to respond to mean stress σ m in the usual manner: Δ p = B Δσ m ( B ranging from 0.4 to 0.7), nearly independent of the ambient deviatoric stress state. However, variations in deviatoric stress (σ d = (σ 1 − σ 3 )/2) at constant mean stress were also found to induce a reversible (elastic) pore pressure response to stress levels up to and exceeding 80% failure stress (i.e., Δ p = = const.). The coefficient η became more negative with increasing deviatoric stress level in sandstone and Ottawa sand samples. That is, η represents a dilatant response where increased deviatoric stress causes a decrease in pore pressure. The poroelastic response to deviatoric stress is explained in terms of anisotropic matrix stiffening due to closure of crack-like pore space or flattening of grain contacts at high ambient stress levels and can be important in calculations of earthquake stress transfer.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Concordant paleolatitudes for Neoproterozoic ophiolitic rocks of the Trinity Complex, Klamath Mountains, California

New paleomagnetic results from the eastern Klamath Mountains of northern California show that Neoproterozoic rocks of the Trinity ophiolitic complex and overlying Middle Devonian volcanic rocks are latitudinally concordant with cratonal North America. Combining paleomagnetic data with regional geologic and faunal evidence suggests that the Trinity Complex and related terranes of the eastern Klamath plate were linked in some fashion to the North American craton throughout that time, but that distance between them may have varied considerably. A possible model that is consistent with our paleomagnetic results and the geologic evidence is that the Trinity Complex formed and migrated parallel to paleolatitude in the basin between Laurasia and Australia–East Antarctica as the Rodinian supercontinent began to break up. It then continued to move parallel to paleolatitude at least through Middle Devonian time. Although the eastern Klamath plate served as a nucleus against which more western components of the Klamath Mountains province amalgamated, the Klamath superterrane was not accreted to North America until Early Cretaceous time.

California↗
USGS2003

Importance of a stochastic distribution of floods and erosion thresholds in the bedrock river incision problem

Fluvial erosion of bedrock occurs during occasional flood events when boundary shear stress exceeds a critical threshold to initiate incision. Therefore efforts to model the evolution of topography over long timescales should include an erosion threshold and should be driven by a stochastic distribution of erosive events. However, most bedrock incision models ignore the threshold as a second‐order detail. In addition, climate is poorly represented in most landscape evolution models, so the quantitative relationship between erosion rate and measurable climatic variables has been elusive. Here we show that the presence of an erosion threshold, when combined with a well‐constrained, probabilistic model of storm and flood occurrence, has first‐order implications for the dynamics of river incision in tectonically active areas. First, we make a direct calculation of the critical shear stress required to pluck bedrock blocks for a field site in New York. Second, we apply a recently proposed stochastic, threshold, bedrock incision model to a series of streams in California, with known tectonic and climatic forcing. Previous work in the area has identified a weak relationship between channel gradient or relief and rock uplift rate that is not easily explained by simpler detachment‐limited models. The results with the stochastic threshold model show that even low erosion thresholds, which are exceeded in steep channels during high‐frequency flood events, fundamentally affect the predicted relationship between gradient and uplift rate in steady state rivers, in a manner consistent with the observed topography. This correspondence between theory and data is, however, nonunique; models in which a thin alluvial cover may act to inhibit channel incision in the low uplift rate zone also provide plausible explanations for the observed topography. Third, we explore the broader implications of the stochastic threshold model to the development of fluvial topography in active tectonic settings. We suggest that continued field applications of geomorphic models, including physically meaningful thresholds and stochastic climate distributions, are required to advance our knowledge of interactions among surficial, climatic, and crustal processes.

California↗
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / n , n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Fire dynamics and implications for nitrogen cycling in boreal forests

We used a dynamic, long-term mass balance approach to track cumulative carbon (C) and nitrogen (N) losses to fire in boreal Manitoba over the 6500 years since deglaciation. Estimated C losses to decomposition and fire, combined with measurements of N pools in mature and burned forest floors, suggest that loss of N by combustion has likely resulted in a long-term loss that exceeds the amount of N stored in soil today by 2 to 3 times. These estimates imply that biological N fixation rates could be as high as 5 to 10 times atmospheric deposition rates in boreal regions. At the site scale, the amount of N lost is due to N content of fuels, which varies by stand type and fire severity, which in turn vary with climate and fire dynamics. The interplay of fire frequency, fire severity, and N partitioning during regrowth are important for understanding rates and sustainability of nutrient and carbon cycling over millenia and over broad regions.

Manitoba↗
USGS2001

Moss and lichen cover mapping at local and regional scales in the boreal forest ecosystem of central Canada

Mosses and lichens are important components of boreal landscapes [Vitt et al., 1994; Bubier et al., 1997]. They affect plant productivity and belowground carbon sequestration and alter the surface runoff and energy balance. We report the use of multiresolution satellite data to map moss and lichens over the BOREAS region at a 10 m, 30 m, and 1 km scales. Our moss and lichen classification at the 10 m scale is based on ground observations of associations among soil drainage classes, overstory composition, and cover type among four broad classes of ground cover (feather, sphagnum, and brown mosses and lichens). For our 30 m map, we used field observations of ground cover-overstory associations to map mosses and lichens in the BOREAS southern study area (SSA). To scale up to a 1 km (AVHRR) moss map of the BOREAS region, we used the TM SSA mosaics plus regional field data to identify AVHRR overstory-ground cover associations. We found that: 1) ground cover, overstory composition and density are highly correlated, permitting inference of moss and lichen cover from satellite-based land cover classifications; 2) our 1 km moss map reveals that mosses dominate the boreal landscape of central Canada, thereby a significant factor for water, energy, and carbon modeling; 3) TM and AVHRR moss cover maps are comparable; 4) satellite data resolution is important; particularly in detecting the smaller wetland features, lakes, and upland jack pine sites; and 5) distinct regional patterns of moss and lichen cover correspond to latitudinal and elevational gradients. Copyright 2001 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Stable isotope composition of waters in the Great Basin, United States 1. Air-mass trajectories

Isentropic trajectories, calculated using the NOAA/Climate Monitoring and Diagnostics Laboratory's isentropic transport model, were used to determine air-parcel origins and the influence of air mass trajectories on the isotopic composition of precipitation events that occurred between October 1991 and September 1993 at Cedar City, Utah, and Winnemucca, Nevada. Examination of trajectories that trace the position of air parcels backward in time for 10 days indicated five distinct regions of water vapor origin: (1) Gulf of Alaska and North Pacific, (2) central Pacific, (3) tropical Pacific, (4) Gulf of Mexico, and (5) continental land mass. Deuterium (??D) and oxygen-18 (??18O) analyses were made of precipitation representing 99% of all Cedar City events. Similar analyses were made on precipitation representing 66% of the precipitation falling at Winnemucca during the same period. The average isotopic composition of precipitation derived from each water vapor source was determined. More than half of the precipitation that fell at both sites during the study period originated in the tropical Pacific and traveled northeast to the Great Basin; only a small proportion traversed the Sierra Nevada. The isotopic composition of precipitation is determined by air-mass origin and its track to the collection station, mechanism of droplet formation, reequilibration within clouds, and evaporation during its passage from cloud to ground. The Rayleigh distillation model can explain the changes in isotopic composition of precipitation as an air mass is cooled pseudo-adiabatically during uplift. However, the complicated processes that take place in the rapidly convecting environment of cumulonimbus and other clouds that are common in the Great Basin, especially in summer, require modification of this model because raindrops that form in the lower portion of those clouds undergo isotopic change as they are elevated to upper levels of the clouds from where they eventually drop to the ground.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Stable isotope compositions of waters in the Great Basin, United States 3. Comparison of groundwaters with modern precipitation

Groundwater samples from wells and springs, scattered over most of the Great Basin province, were collected and analyzed for their isotopic makeup. They were augmented by previously published isotopic data on groundwaters from southeast California and by several hundred unpublished isotopic analyses. The ratio of 2H (deuterium, D) to 1H, in water samples from valleys in parts of California, Idaho, Nevada, Oregon, and Utah, are here compared with the winter, summer, and annual isotopic compositions of precipitation falling in or near the sampled areas. The main goal of this study was to identify basins where the groundwaters have isotopic compositions that are "lighter" (depleted in the heavier isotope, D) relative to modern winter precipitation. Where these basins do not adjoin substantially higher terrain, we consider those light groundwaters to be of Pleistocene age and thus more than 10,000 years old. Where the groundwater is 10 to 19??? lighter than local winter precipitation, we consider it to be possibly an indication of Pleistocene water; where the ??D makeup is >20??? lighter, we consider it to be probably Pleistocene water. More than 80 sites underlain by waters of possible or probable Pleistocene age were identified.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Modeling soil thermal and carbon dynamics of a fire chronosequence in interior Alaska

In this study, the dynamics of soil thermal, hydrologic, and ecosystem processes were coupled to project how the carbon budgets of boreal forests will respond to changes in atmospheric CO 2 , climate, and fire disturbance. The ability of the model to simulate gross primary production and ecosystem respiration was verified for a mature black spruce ecosystem in Canada, the age-dependent pattern of the simulated vegetation carbon was verified with inventory data on aboveground growth of Alaskan black spruce forests, and the model was applied to a postfire chronosequence in interior Alaska. The comparison between the simulated soil temperature and field-based estimates during the growing season (May to September) of 1997 revealed that the model was able to accurately simulate monthly temperatures at 10 cm ( R > 0.93) for control and burned stands of the fire chronosequence. Similarly, the simulated and field-based estimates of soil respiration for control and burned stands were correlated ( R = 0.84 and 0.74 for control and burned stands, respectively). The simulated and observed decadal to century-scale dynamics of soil temperature and carbon dynamics, which are represented by mean monthly values of these variables during the growing season, were correlated among stands ( R = 0.93 and 0.71 for soil temperature at 20- and 10-cm depths, R = 0.95 and 0.91 for soil respiration and soil carbon, respectively). Sensitivity analyses indicate that along with differences in fire and climate history a number of other factors influence the response of carbon dynamics to fire disturbance. These factors include nitrogen fixation, the growth of moss, changes in the depth of the organic layer, soil drainage, and fire severity.

Alaska↗
USGS2002

Thermal infrared spectroscopy of experimentally shocked anorthosite and pyroxenite: Implications for remote sensing of Mars

The feldspar and pyroxene mineralogies on Mars revealed by the Thermal Emission Spectrometer (TES) on Mars Global Surveyor likely record a variety of shock effects, as suggested by petrologic analyses of the Martian meteorites and the abundance of impact craters on the planet's surface. To study the effects of shock pressures on thermal infrared spectra of these minerals, we performed shock recovery experiments on orthopyroxenite and anorthosite samples from the Stillwater Complex (Montana) over peak pressures from 17 to 63 GPa. We acquired emissivity and hemispherical reflectance spectra (350–1400 cm −1 ; ∼7–29 μm) of both coherent chips and fine-grained powders of shocked and unshocked samples. These spectra are more directly comparable to remotely sensed data of Mars (e.g., TES) than previously acquired absorption or transmission spectra of shocked minerals. The spectra of experimentally shocked feldspar show systematic changes with increasing pressure due to depolymerization of the silica tetrahedra. For the spectra of chips, this includes the disappearance of small bands in the 500–650 cm −1 region and a strong band at 1115 cm −1 , and changes in positions of a strong band near 940 cm −1 and the Christiansen feature near 1250 cm −1 . Spectra of the shocked powders show the gradual disappearance of a transparency feature near 830 cm −1 . Fewer changes are observed in the pyroxene spectra at pressures as high as 63 GPa. Spectra of experimentally shocked minerals will help identify more precisely the mineralogy of rocks and soils not only from TES but also from Mars instruments such as miniTES and THEMIS.

Journal of Geophysical Research E: Planets↗
USGS2003

Slip history and dynamic implications of the 1999 Chi-Chi, Taiwan, earthquake

[1] We investigate the rupture process of the 1999 Chi‐Chi, Taiwan, earthquake using extensive near‐source observations, including three‐component velocity waveforms at 36 strong motion stations and 119 GPS measurements. A three‐plane fault geometry derived from our previous inversion using only static data [ Ji et al. , 2001 ] is applied. The slip amplitude, rake angle, rupture initiation time, and risetime function are inverted simultaneously with a recently developed finite fault inverse method that combines a wavelet transform approach with a simulated annealing algorithm [ Ji et al. , 2002b ]. The inversion results are validated by the forward prediction of an independent data set, the teleseismic P and SH ground velocities, with notable agreement. The results show that the total seismic moment release of this earthquake is 2.7 × 10 20 N m and that most of the slip occurred in a triangular‐shaped asperity involving two fault segments, which is consistent with our previous static inversion. The rupture front propagates with an average rupture velocity of ∼2.0 km s −1 , and the average slip duration (risetime) is 7.2 s. Several interesting observations related to the temporal evolution of the Chi‐Chi earthquake are also investigated, including (1) the strong effect of the sinuous fault plane of the Chelungpu fault on spatial and temporal variations in slip history, (2) the intersection of fault 1 and fault 2 not being a strong impediment to the rupture propagation, and (3) the observation that the peak slip velocity near the surface is, in general, higher than on the deeper portion of the fault plane, as predicted by dynamic modeling.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Stress-induced, time-dependent fracture closure at hydrothermal conditions

Time-dependent closure of fractures in quartz was measured in situ at 22–530°C temperature and 0.1-150 MPa water pressure. Unlike previous crack healing and rock permeability studies, in this study, fracture aperture is monitored directly and continuously using a windowed pressure vessel, a long-working-distance microscope, and reflected-light interferometry. Thus the fracture volume and geometry can be measured as a function of time, temperature, and water pressure. Relatively uniform closure occurs rapidly at temperatures and pressures where quartz becomes significantly soluble in water. During closure the aperture is reduced by as much as 80% in a few hours. We infer that this closure results from the dissolution of small particles or asperities that prop the fracture open. The driving force for closure via dissolution of the prop is the sum of three chemical potential terms: (1) the dissolution potential, proportional to the logarithm of the degree of undersaturation of the solution; (2) the coarsening potential, proportional to the radius of curvature of the prop; and (3) the pressure solution potential, proportional to the effective normal stress at the contact between propping particles and the fracture wall. Our observations suggest that closure is controlled by a pressure solution-like process. The aperture of dilatant fractures and microcracks in the Earth that are similar to those in our experiments, such as ones generated from thermal stressing or brittle failure during earthquake rupture and slip, will decrease rapidly with time, especially if the macroscopic stress is nonhydrostatic.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Three-dimensional velocity structure of crust and upper mantle in southwestern China and its tectonic implications

Using P and S arrival times from 4625 local and regional earthquakes recorded at 174 seismic stations and associated geophysical investigations, this paper presents a three‐dimensional crustal and upper mantle velocity structure of southwestern China (21°–34°N, 97°–105°E). Southwestern China lies in the transition zone between the uplifted Tibetan plateau to the west and the Yangtze continental platform to the east. In the upper crust a positive velocity anomaly exists in the Sichuan Basin, whereas a large‐scale negative velocity anomaly exists in the western Sichuan Plateau, consistent with the upper crustal structure under the southern Tibetan plateau. The boundary between these two anomaly zones is the Longmen Shan Fault. The negative velocity anomalies at 50‐km depth in the Tengchong volcanic area and the Panxi tectonic zone appear to be associated with temperature and composition variations in the upper mantle. The Red River Fault is the boundary between the positive and negative velocity anomalies at 50‐km depth. The overall features of the crustal and the upper mantle structures in southwestern China are a low average velocity, large crustal thickness variations, the existence of a high‐conductivity layer in the crust or/and upper mantle, and a high heat flow value. All these features are closely related to the collision between the Indian and the Asian plates.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Basin-centered asperities in great subduction zone earthquakes: A link between slip, subsidence, and subduction erosion?

Published areas of high coseismic slip, or asperities, for 29 of the largest Circum-Pacific megathrust earthquakes are compared to forearc structure revealed by satellite free-air gravity, bathymetry, and seismic profiling. On average, 71% of an earthquake's seismic moment and 79% of its asperity area occur beneath the prominent gravity low outlining the deep-sea terrace; 57% of an earthquake's asperity area, on average, occurs beneath the forearc basins that lie within the deep-sea terrace. In SW Japan, slip in the 1923, 1944, 1946, and 1968 earthquakes was largely centered beneath five forearc basins whose landward edge overlies the 350??C isotherm on the plate boundary, the inferred downdip limit of the locked zone. Basin-centered coseismic slip also occurred along the Aleutian, Mexico, Peru, and Chile subduction zones but was ambiguous for the great 1964 Alaska earthquake. Beneath intrabasin structural highs, seismic slip tends to be lower, possibly due to higher temperatures and fluid pressures. Kilometers of late Cenozoic subsidence and crustal thinning above some of the source zones are indicated by seismic profiling and drilling and are thought to be caused by basal subduction erosion. The deep-sea terraces and basins may evolve not just by growth of the outer arc high but also by interseismic subsidence not recovered during earthquakes. Basin-centered asperities could indicate a link between subsidence, subduction erosion, and seismogenesis. Whatever the cause, forearc basins may be useful indicators of long-term seismic moment release. The source zone for Cascadia's 1700 A.D. earthquake contains five large, basin-centered gravity lows that may indicate potential asperities at depth. The gravity gradient marking the inferred downdip limit to large coseismic slip lies offshore, except in northwestern Washington, where the low extends landward beneath the coast. Transverse gravity highs between the basins suggest that the margin is seismically segmented and could produce a variety of large earthquakes. Published in 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Magma supply dynamics at Westdahl volcano, Alaska, modeled from satellite radar interferometry

A group of satellite radar interferograms that span the time period from 1991 to 2000 shows that Westdahl volcano, Alaska, deflated during its 1991-1992 eruption and is reinflating at a rate that could produce another eruption within the next several years. The rates of inflation and deflation are approximated by exponential decay functions having time constants of about 6 years and a few days, respectively. This behavior is consistent with a deep, constant-pressure magma source connected to a shallow reservoir by a magma-filled conduit. An elastic deformation model indicates that the reservoir is located about 6 km below sea level and beneath Westdahl Peak. We propose that the magma flow rate through the conduit is governed by the pressure gradient between the deep source and the reservoir. The pressure gradient, and hence the flow rate, are greatest immediately after eruptions. Pressurization of the reservoir decreases both the pressure gradient and the flow rate, but eventually the reservoir ruptures and an eruption or intrusion ensues. The eruption rate is controlled partly by the pressure gradient between the reservoir and surface, and therefore it, too, decreases with time. When the supply of eruptible magma is exhausted, the eruption stops, the reservoir begins to repressurize at a high rate, and the cycle repeats. This model might also be appropriate for other frequently active volcanoes with stable magma sources and relatively simple magma storage systems.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Oxygen-18 concentrations in recent precipitation and ice cores on the Tibetan Plateau

A detailed study of the climatic significance of δ 18 O in precipitation was completed on a 1500 km southwest-northeast transect of the Tibetan Plateau in central Asia. Precipitation samples were collected at four meteorological stations for up to 9 years. This study shows that the gradual impact of monsoon precipitation affects the spatial variation of δ 18 O-T relationship along the transect. Strong monsoon activity in the southern Tibetan Plateau results in high precipitation rates and more depleted heavy isotopes. This depletion mechanism is described as a precipitation “amount effect” and results in a poor δ 18 O-T relationship at both seasonal and annual scales. In the middle of the Tibetan Plateau, the effects of the monsoon are diminished but continue to cause a reduced correlation of δ 18 O and temperature at the annual scale. At the monthly scale, however, a significant δ 18 O-T relationship does exist. To the north of the Tibetan Plateau beyond the extent of the effects of monsoon precipitation, δ 18 O in precipitation shows a strong temperature dependence. δ 18 O records from two shallow ice cores and historic air temperature data were compared to verify the modern δ 18 O-T relationship. δ 18 O in Dunde ice core was positively correlated with air temperature from a nearby meteorological station in the north of the plateau. The δ 18 O variation in an ice core from the southern Plateau, however, was inversely correlated with precipitation amount at a nearby meteorological station and also the accumulation record in the ice core. The long-term variation of δ 18 O in the ice core record in the monsoon regions of the southern Tibetan Plateau suggest past monsoon seasons were probably more expansive. It is still unclear, however, how changes in large-scale atmosphere circulation might influence summer monsoon precipitation on the Tibetan Plateau.

Journal of Geophysical Research D: Atmospheres↗
USGS2003

Multiple oxygen and sulfur isotopic analyses on water-soluble sulfate in bulk atmospheric deposition from the southwestern United States

Sulfate is a major component of bulk atmospheric deposition (including dust, aerosol, fog, and rain). We analyzed sulfur and oxygen isotopic compositions of water-soluble sulfate from 40 sites where year-round dust traps collect bulk atmospheric deposition in the southwestern United States. Average sulfur and oxygen isotopic compositions (δ 34 S and δ 18 O) are 5.8 ± 1.4 (CDT) and 11.2 ± 1.9 (SMOW) (n = 47), respectively. Samples have an oxygen 17 anomaly (Δ 17 O), with an average value of 1.0 ± 0.6‰. Except for a weak positive correlation between δ 18 O and Δ 17 O values (r 2 ≈ 0.4), no correlation exists for δ 18 O versus δ 34 S, Δ 17 O versus δ 34 S, or any of the three isotopic compositions versus elevation of the sample site. Exceptional positive Δ 17 O values (up to 4.23‰) are found in samples from sites in the vicinity of large cities or major highways, and near-zero Δ 17 O values are found in samples close to dry lakes. Comparison of isotopic values of dust trap sulfate and desert varnish sulfate from the region reveals that varnish sulfate has average isotopic values that are ∼4.8‰ lower for δ 18 O, ∼2.1‰ higher for δ 34 S, and ∼0.3‰ lower for Δ 17 O than those of the present-day bulk deposition sulfate. Although other factors could cause the disparity, this observation suggests a possibility that varnish sulfate may have recorded a long-term atmospheric sulfate deposition during the Holocene or Pleistocene, as well as the differences between sulfur and oxygen isotopic compositions of the preindustrial bulk deposition sulfate and those of the industrial era.

Journal of Geophysical Research D: Atmospheres↗
USGS2003

Orogenic tectonism on Io

We catalog 143 Ionian mountains (montes) and mountain‐like features (mensae, tholi, plana, and small peaks) in order to investigate orogenic tectonism on Io. From this comprehensive list, we select 96 mountains for which there are sufficient coverage and resolution to discern spatial relationships with surrounding geologic features. Three of the 96 mountains are probably volcanoes, 92 appear to be tectonic massifs, and 1 is ambiguous. Of the 92 tectonic mountains, 38 abut paterae (volcanic or volcano‐tectonic craters with irregular or scalloped margins). This juxtaposition is unlikely to be a coincidence as the probability of it occurring by chance is ∼0.1%. We propose instead that orogenic faults may act as conduits for magma ascent, thus fueling patera formation near mountains. As resurfacing buries a shell of material from Io's surface to the base of the lithosphere, its effective radius is reduced and it heats up. We calculate the lithospheric volume change due to subsidence and thermal expansion as a function of lithospheric thickness. Conservation of volume dictates that this material must be uplifted at Io's surface. By estimating the total volume of the mountains, we are able to place a lower limit of 12 km on Io's lithospheric thickness. We hypothesize that, in some cases, mountain formation may be facilitated by asthenospheric diapirs impinging on the base of the lithosphere. The resulting lithospheric swell could focus the compressive stresses that drive orogenic tectonism. This model is one of several possible mechanisms for uplifting isolated mountains such as are observed on Io.

Journal of Geophysical Research E: Planets↗