Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Geophysical Research- Earth Surface”

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 235 records · Page 13Linked to original sources

Seismic and geodetic signatures of fault slip at the Slumgullion Landslide Natural Laboratory

We tested the hypothesis that the Slumgullion landslide is a useful natural laboratory for observing fault slip, specifically that slip along its basal surface and side-bounding strike-slip faults occurs with comparable richness of aseismic and seismic modes as along crustal- and plate-scale boundaries. Our study provides new constraints on models governing landslide motion. We monitored landslide deformation with temporary deployments of a 29-element prism array surveyed by a robotic theodolite and an 88-station seismic network that complemented permanent extensometers and environmental instrumentation. Aseismic deformation observations show that large blocks of the landslide move steadily at approximately centimeters per day, possibly punctuated by variations of a few millimeters, while localized transient slip episodes of blocks less than a few tens of meters across occur frequently. We recorded a rich variety of seismic signals, nearly all of which originated outside the monitoring network boundaries or from the side-bounding strike-slip faults. The landslide basal surface beneath our seismic network likely slipped almost completely aseismically. Our results provide independent corroboration of previous inferences that dilatant strengthening along sections of the side-bounding strike-slip faults controls the overall landslide motion, acting as seismically radiating brakes that limit acceleration of the aseismically slipping basal surface. Dilatant strengthening has also been invoked in recent models of transient slip and tremor sources along crustal- and plate-scale faults suggesting that the landslide may indeed be a useful natural laboratory for testing predictions of specific mechanisms that control fault slip at all scales.

Colorado

Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity

A new model of the subducted Juan de Fuca plate beneath western North America allows first-order correlations between the occurrence of Wadati-Benioff zone earthquakes and slab geometry, temperature, and hydration state. The geo-referenced 3D model, constructed from weighted control points, integrates depth information from earthquake locations and regional seismic velocity studies. We use the model to separate earthquakes that occur in the Cascadia forearc from those that occur within the underlying Juan de Fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes. Seismicity within the slab is most prevalent where the slab is warped beneath northwestern California and western Washington suggesting that slab flexure, in addition to expected metamorphic dehydration processes, promotes earthquake occurrence within the subducted oceanic plate. Earthquake patterns beneath western Vancouver Island are consistent with slab dehydration processes. Conversely, the lack of slab earthquakes beneath western Oregon is consistent with an anhydrous slab. Double-differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern California that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 M7.1 Mendocino earthquake. We infer that this earthquake ruptured a surface within the Cascadia accretionary margin above the Juan de Fuca plate. We further speculate that this earthquake is associated with a detached fragment of former Farallon plate. Other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes.

British Columbia;California;Oregon;Washington

Detection of aseismic slip and poroelastic reservoir deformation at the North Brawley Geothermal Field from 2009 to 2019

The North Brawley Geothermal Field, located within the Brawley Seismic Zone of Southern California, presents a case study for understanding seismic hazards linked to fluid injection and geothermal energy extraction. An earthquake swarm near the geothermal field in 2012 included two earthquakes with magnitudes greater than 5 and was potentially preceded by a years-long aseismic slip transient. To better understand ground deformation around the geothermal field, including its evolution with time and its physical mechanisms, we analyze deformation before, during, and after the swarm using ground- and satellite-based geodetic techniques between 2009 and 2019. We integrate observations from GNSS, Sentinel-1, TerraSAR-X, UAVSAR, and leveling surveys into a single deformation history. Modeling of this new collection of observations at the North Brawley Geothermal Field provides evidence for 80% more pre-swarm aseismic slip than previously recognized from 2009 to 2012. During the 2012 Brawley swarm, our geodetic slip inversions closely match the results of seismic waveform inversions from the swarm events. After the 2012 swarm, surface deformation is dominated by poroelastic deformation of a shallow fluid reservoir at <1 km depth rather than fault slip. The deformation history and seismicity catalogs at North Brawley suggest a cessation of fault-related slip during the ∼7 years after the 2012 earthquake swarm.

California

Uplift and subsidence associated with the great Aceh-Andaman earthquake of 2004

Rupture of the Sunda megathrust on 26 December 2004 produced broad regions of uplift and subsidence. We define the pivot line separating these regions as a first step in defining the lateral extent and the downdip limit of rupture during that great Mw ??? 9.2 earthquake. In the region of the Andaman and Nicobar islands we rely exclusively on the interpretation of satellite imagery and a tidal model. At the southern limit of the great rupture we rely principally on field measurements of emerged coral microatolls. Uplift extends from the middle of Simeulue Island, Sumatra, at ??? 2.5??N, to Preparis Island, Myanmar (Burma), at ??? 14.9??N. Thus the rupture is ??? 1600 km long. The distance from the pivot line to the trench varies appreciably. The northern and western Andaman Islands rose, whereas the southern and eastern portion of the islands subsided. The Nicobar Islands and the west coast of Aceh province, Sumatra, subsided. Tilt at the southern end of the rupture is steep; the distance from 1.5 m of uplift to the pivot line is just 60 km. Our method of using satellite imagery to recognize changes in elevation relative to sea surface height and of using a tidal model to place quantitative bounds on coseismic uplift or subsidence is a novel approach that can be adapted to other forms of remote sensing and can be applied to other subduction zones in tropical regions. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth

Implications of regional gravity for state of stress in the earth's crust and upper mantle

Topography is maintained by stress differences within the earth. Depending on the distribution of the stress we classify the support as either local or regional compensation. In general, the stresses implied in a regional compensation scheme are an order of magnitude larger than those corresponding to local isostasy. Gravity anomalies, a measure of the earth's departure from hydrostatic equilibrium, can be used to distinguish between the two compensation mechanisms and thus to estimate the magnitude of deviatoric stress in the crust and upper mantle. Topography created at an oceanic ridge crest or in a major continental orogenic zone appears to be locally compensated. Such features were formed on weak crust incapable of maintaining stress differences much greater than the stress from the applied load. Oceanic volcanoes formed on an already cooled, thickened lithosphere are regionally supported with elastic stresses. Similarly, the broad topographic rise seaward of subduction zones is elastically supported as the lithosphere is bent near the plate margin. Although the implied stress is to some degree dependent on the rheological model assumed, the gravity anomalies and surface deformation produced by these features demonstrate that the upper 30–40 km of the oceanic lithosphere is capable of regionally supporting stress differences in the 100-MPa range. Given certain conditions of load emplacement, continental crust can also support loads regionally over 100-m.y. time scales, but the effects of erosion only allow an estimate of a lower bound on stress. Data from space probes indicate that the upper layers of other terrestrial planets also support topographic-induced stress differences in excess of 100 MPa.

Journal of Geophysical Research Solid Earth

Near-surface versus fault zone damage following the 1999 Chi-Chi earthquake: Observation and simulation of repeating earthquakes

We observe crustal damage and its subsequent recovery caused by the 1999&thinsp; M 7.6 Chi-Chi earthquake in central Taiwan. Analysis of repeating earthquakes in Hualien region, ~70&thinsp;km east of the Chi-Chi earthquake, shows a remarkable change in wave propagation beginning in the year 2000, revealing damage within the fault zone and distributed across the near surface. We use moving window cross correlation to identify a dramatic decrease in the waveform similarity and delays in the S wave coda. The maximum delay is up to 59&thinsp;ms, corresponding to a 7.6% velocity decrease averaged over the wave propagation path. The waveform changes on either side of the fault are distinct. They occur in different parts of the waveforms, affect different frequencies, and the size of the velocity reductions is different. Using a finite difference method, we simulate the effect of postseismic changes in the wavefield by introducing S wave velocity anomaly in the fault zone and near the surface. The models that best fit the observations point to pervasive damage in the near surface and deep, along-fault damage at the time of the Chi-Chi earthquake. The footwall stations show the combined effect of near-surface and the fault zone damage, where the velocity reduction (2&ndash;7%) is twofold to threefold greater than the fault zone damage observed in the hanging wall stations. The physical models obtained here allow us to monitor the temporal evolution and recovering process of the Chi-Chi fault zone damage.

Journal of Geophysical Research B: Solid Earth

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth

The influence of formation material properties on the response of water levels in wells to Earth tides and atmospheric loading

The water level in an open well can change in response to deformation of the surrounding material, either because of applied strains (tidal or tectonic) or surface loading by atmospheric pressure changes. Under conditions of no vertical fluid flow and negligible well bore storage (static-confined conditions), the sensitivities to these effects depend on the elastic properties and porosity which characterize the surrounding medium. For a poroelastic medium, high sensitivity to applied areal strains occurs for low porosity, while high sensitivity to atmospheric loading occurs for high porosity; both increase with decreasing compressibility of the solid matrix. These material properties also influence vertical fluid flow induced by areally extensive deformation and can be used to define two types of hydraulic diffusivity which govern pressure diffusion, one for applied strain and one for surface loading. The hydraulic diffusivity which governs pressure diffusion in response to surface loading is slightly smaller than that which governs fluid flow in response to applied strain. Given the static-confined response of a water well to atmospheric loading and Earth tides, the in situ drained matrix compressibility and porosity (and hence the one-dimensional specific storage) can be estimated. Analysis of the static-confined response of five wells to atmospheric loading and Earth tides gives generally reasonable estimates for material properties.

Journal of Geophysical Research Solid Earth

Geodetic measurement of crustal deformation on the San Andreas, Hayward, and Calaveras faults near San Francisco, California

Analysis of a geodetic network of 115 lines crossing the San Andreas, Hayward, and Calaveras faults in the vicinity of San Francisco Bay and measured repeatedly between 1970 and 1980 has revealed details about the accommodation of relative plate motion in this area. The most striking result is that the deformation is not uniformly distributed across the area. In the east bay, along the Hayward and Calaveras faults, all motion appears to take place as slip directly on the fault, with no accumulation of strain in the adjacent crust. On both the Calaveras and the Hayward faults the rate obtained for the 1970 to 1980 period agrees with geologic rates spanning a few million years and with creep rates spanning a few decades. The Hayward fault slip rate is 7 ± 1 mm/yr. The Calaveras fault slip rate is 7 ± 1 mm/yr, with perhaps half of this slip distributed across a zone a few kilometers wide, probably as inelastic deformation of weak near‐surface material. The absence of strain accumulation in the east bay is surprising since the Hayward and Calaveras faults have been the site of large earthquakes in the past. A block located east of the Calaveras fault and south of the Las Positas fault has been rotating clockwise at a rate of 0.3 ± 0.1 μrad/yr with very little internal deformation. Along the San Francisco peninsula no detectable slip occurs (less than 1.5 mm/yr) at the surface, but appreciable strain is accumulating. Near fault shear strain rates are 0.6 ± 0.1 μstrain/yr (engineering) with direction N47°W ± 9. The slip rate near the San Andreas fault is 12.2 ± 3.9 mm/yr distributed across a broad zone. The relative motion across the whole region during the period 1970–1980 is 32.1 ± 7.4 mm/yr.

California

Subsidence rates at the southern Salton Sea consistent with reservoir depletion

Space geodetic measurements from the Envisat satellite between 2003 and 2010 show that subsidence rates near the southeastern shoreline of the Salton Sea in Southern California are up to 52mmyr −1 greater than the far-field background rate. By comparing these measurements with model predictions, we find that this subsidence appears to be dominated by poroelastic contraction associated with ongoing geothermal fluid production, rather than the purely fault-related subsidence proposed previously. Using a simple point source model, we suggest that the source of this proposed volumetric strain is at depths between 1.0 km and 2.4 km (95% confidence interval), comparable to generalized boundaries of the Salton Sea geothermal reservoir. We find that fault slip on two previously imaged tectonic structures, which are part of a larger system of faults in the Brawley Seismic Zone, is not an adequate predictor of surface velocity fields because the magnitudes of the best fitting slip rates are often greater than the full plate boundary rate and at least 2 times greater than characteristic sedimentation rates in this region. Large-scale residual velocity anomalies indicate that spatial patterns predicted by fault slip are incompatible with the observations.

California

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

Mineralogy and geochemistry of a sediment‐hosted hydrothermal sulfide deposit from the Southern Trough of Guaymas Basin, Gulf of California

Samples dredged from a 15‐m‐high hydrothermal mound atop the flat turbidite pond in the Southern Trough of Guaymas Basin consist of pyrrhotite‐rich massive sulfide, barite, barite + calcite, talc, and opaline silica as well as substrate material composed of fossiliferous, clay‐rich ooze. An 11‐m‐long sediment core taken near the dredge site shows increasing hydrothermal alteration with depth; anhydrite‐filled fractures near the base of the core appear to be channels for hydrothermal discharge. Oxidation of the sulfide‐rich samples to an assemblage of geothite, lepidocrocite, and amorphous Fe oxyhydroxide is ubiquitous. Compared to other massive sulfide deposits on sediment‐starved oceanic ridges, the hydrothermal deposit dredged in Guaymas Basin has a high pyrrhotite/pyrite ratio, a low Zn sulfide and combined ore metal (Cu + Zn + Pb + Ag + Cd) content, and a greater abundance of sulfate, carbonate, and silicate phases. Venting hydrothermal solutions are alkaline with moderately high p H; high Ca, Ba, and SiO 2 content; low ƒS 2 and ƒo 2 ; and very low transition metal content. Disequilibrium assemblages of pyrrhotite and sulfate minerals form during rapid mixing of this evolved vent fluid with ambient bottom waters at the discharge site. Talc is formed at a temperature near 270°C by mixing or entrainment of Mg‐rich bottom water or pore fluid with upwelling hydrothermal fluid that is saturated with silica. Calcite may precipitate from the alkaline, Ca‐rich fluid during degassing of CO 2 . The minimum temperature range for sulfide and nonsulfide deposition is approximately 190°–326°C. The composition of hydrothermal deposits, vent solutions, and altered sediment requires that circulating fluids evolve during deep penetration into the basaltic basement complex, further interaction with the organic‐and carbonate‐rich sediment pile, and near‐surface mixing with ambient seawater. Although the stable assemblage albite‐epidote‐clinochlore present at depth in the sediment pile requires very low dissolved Mg and Fe in the altering fluid, the addition of Mg to deeply buried sediment indicates significant recharge of the system by Guaymas Basin bottom water.

Southern Trough of Guaymas Basin

The role of viscous magma mush spreading in volcanic flank motion at Kīlauea Volcano, Hawai‘i

Multiple mechanisms have been suggested to explain seaward motion of the south flank of Kīlauea Volcano, Hawai‘i. The consistency of flank motion during both waxing and waning magmatic activity at Kīlauea suggests that a continuously acting force, like gravity body force, plays a substantial role. Using finite element models, we test whether gravity is the principal driver of long-term motion of Kīlauea's flank. We compare our model results to geodetic data from Global Positioning System and interferometric synthetic aperture radar during a time period with few magmatic and tectonic events (2000-2003), when deformation of Kīlauea was dominated by summit subsidence and seaward motion of the south flank. We find that gravity-only models can reproduce the horizontal surface velocities if we incorporate a regional décollement fault and a deep, low-viscosity magma mush zone. To obtain quasi steady state horizontal surface velocities that explain the long-term seaward motion of the flank, we find that an additional weak zone is needed, which is an extensional rift zone above the magma mush. The spreading rate in our model is mainly controlled by the magma mush viscosity, while its density plays a less significant role. We find that a viscosity of 2.5 × 1017–2.5 × 1019 Pa s for the magma mush provides an acceptable fit to the observed horizontal surface deformation. Using high magma mush viscosities, such as 2.5 × 1019 Pa s, the deformation rates remain more steady state over longer time scales. These models explain a significant amount of the observed subsidence at Kīlauea's summit. Some of the remaining subsidence is probably a result of magma withdrawal from subsurface reservoirs.

Hawai'i

New constraints on mechanisms of remotely triggered seismicity at Long Valley Caldera

Regional-scale triggering of local earthquakes in the crust by seismic waves from distant main shocks has now been robustly documented for over a decade. Some of the most thoroughly recorded examples of repeated triggering of a single site from multiple, large earthquakes are measured in geothermal fields of the western United States like Long Valley Caldera. As one of the few natural cases where the causality of an earthquake sequence is apparent, triggering provides fundamental constraints on the failure processes in earthquakes. We show here that the observed triggering by seismic waves is inconsistent with any mechanism that depends on cumulative shaking as measured by integrated energy density. We also present evidence for a frequency-dependent triggering threshold. On the basis of the seismic records of 12 regional and teleseismic events recorded at Long Valley Caldera, long-period waves (>30 s) are more effective at generating local seismicity than short-period waves of comparable amplitude. If the properties of the system are stationary over time, the failure threshold for long-period waves is ~0.05 cm/s vertical shaking. Assuming a phase velocity of 3.5 km/s and an elastic modulus of 3.5 x 10 10 Pa, the threshold in terms of stress is 5 kPa. The frequency dependence is due in part to the attenuation of the surface waves with depth. Fluid flow through a porous medium can produce the rest of the observed frequency dependence of the threshold. If the threshold is not stationary with time, pore pressures that are >99.5% of lithostatic and vary over time by a factor of 4 could explain the observations with no frequency dependence of the triggering threshold.

California

Location of seismic events and eruptive fissures on the Piton de la Fournaise volcano using seismic amplitudes

We present a method for locating the source of seismic events on Piton de la Fournaise. The method is based on seismic amplitudes corrected for station site effects using coda site amplification factors. Once corrected, the spatial distribution of amplitudes shows smooth and simple contours for many types of events, including rockfalls, long-period events and eruption tremor. On the basis of the simplicity of these distributions we develop inversion methods for locating their origins. To achieve this, the decrease of the amplitude as a function of the distance to the source is approximated by the decay either of surface or body waves in a homogeneous medium. The method is effective for locating rockfalls, long-period events, and eruption tremor sources. The sources of eruption tremor are usually found to be located at shallow depth and close to the eruptive fissures. Because of this, our method is a useful tool for locating fissures at the beginning of eruptions.

Journal of Geophysical Research B: Solid Earth

Externally triggered renewed bubble nucleation in basaltic magma: the 12 October 2008 eruption at Halema‘uma‘u Overlook vent, Kīlauea, Hawai‘i, USA

From October 2008 until present, dozens of small impulsive explosive eruptions occurred from the Overlook vent on the southeast side of Halema&lsquo;uma&lsquo;u Crater, at Kīlauea volcano, USA. These eruptions were triggered by rockfalls from the walls of the volcanic vent and conduit onto the top of the lava column. Here we use microtextural observations and data from clasts erupted during the well-characterized 12 October 2008 explosive eruption at Halema&lsquo;uma&lsquo;u to extend existing models of eruption triggering. We present a potential mechanism for this eruption by combining microtextural observations with existing geophysical and visual data sets. We measure the size and number density of bubbles preserved in juvenile ejecta using 2D images and X-ray microtomography. Our data suggest that accumulations of large bubbles with diameters of >50 &mu; m to at least millimeters existed at shallow levels within the conduit prior to the 12 October 2008 explosion. Furthermore, a high number density of small bubbles <50 &mu; m is measured in the clasts, implying very rapid nucleation of bubbles. Visual observations, combined with preexisting geophysical data, suggest that the impact of rockfalls onto the magma free surface induces pressure changes over short timescales that (1) nucleated new additional bubbles in the shallow conduit leading to high number densities of small bubbles and (2) expanded the preexisting bubbles driving upward acceleration. The trigger of eruption and bubble nucleation is thus external to the degassing system.

Hawaii

Postseismic deformation following the 1989 (M = 7.1): Loma Prieta, California, earthquake

Postseismic deformation along a 90‐km profile bisecting the projected surface trace of the coseismic rupture of the 1989 Loma Prieta earthquake has been monitored by frequent GPS surveys for 3.3 years following the earthquake. In addition to the expected deformation associated with secular strain accumulation on the San Andreas and Calaveras faults, deformation associated with postseismic readjustment has been detected. Most of that deformation can be attributed to 1.5 m right‐lateral and 0.9 m reverse postseismic slip on a 5‐km‐wide downdip extension of the Loma Prieta rupture. In addition, there seems to be a 0.1 m postseismic collapse of the Loma Prieta rupture zone in the direction perpendicular to the plane of the rupture. The fault‐normal (N48°E) surface displacements plotted as a function of time exhibit a curvature suggesting a relaxation time of about 1.4 years. Similar plots of the fault‐parallel (N42°W) displacement components do not exhibit significant curvature. Presumably, the deformation shown in those plots is dominated by secular strain accumulation along the San Andreas and Calaveras faults rather than postseismic relaxation.

California

Detection and monitoring of H2O and CO2 ice clouds on Mars

We have developed an observational scheme for the detection and discrimination of Mars atmospheric H 2 O and CO 2 clouds using ground-based instruments in the near infrared. We report the results of our cloud detection and characterization study using Mars near IR images obtained during the 1990 and 1993 oppositions. We focused on specific wavelengths that have the potential, based on previous laboratory studies of H 2 O and CO 2 ices, of yielding the greatest degree of cloud detectability and compositional discriminability. We have detected and mapped absorption features at some of these wavelengths in both the northern and southern polar regions of Mars. Compositional information on the nature of these absorption features was derived from comparisons with laboratory ice spectra and with a simplified radiative transfer model of a CO 2 ice cloud overlying a bright surface. Our results indicate that both H 2 O and CO 2 ices can be detected and distinguished in the polar hood clouds. The region near 3.00 μm is most useful for the detection of water ice clouds because there is a strong H 2 O ice absorption at this wavelength but only a weak CO 2 ice band. The region near 3.33 μm is most useful for the detection of CO 2 ice clouds because there is a strong, relatively narrow CO 2 ice band at this wavelength but only broad “continuum” H 2 O ice absorption. Weaker features near 2.30 μm could arise from CO 2 ice at coarse grain sizes, or surface/dust minerals. Narrow features near 2.00 μm, which could potentially be very diagnostic of CO 2 ice clouds, suffer from contamination by Mars atmospheric CO 2 absorptions and are difficult to interpret because of the rather poor knowledge of surface elevation at high latitudes. These results indicate that future ground-based, Earth-orbital, and spacecraft studies over a more extended span of the seasonal cycle should yield substantial information on the style and timing of volatile transport on Mars, as well as a more detailed understanding of the role of CO 2 condensation in the polar heat budget.

Journal of Geophysical Research E: Planets