SEARCH · Geology Reports

Results for “Earth Imaging Journal”

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

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11Linked to original sources

USGS2024

Vortex trapping of suspended sand grains over ripples

Coastal hydrodynamics and morphodynamics integrate the effects of small-scale fluid-sediment interactions; yet, these small-scale processes are not well understood. To investigate sediment trapping by turbulent coherent structures or vortices, the transport of coarse sand over ripples was analyzed in a small-oscillatory flow tunnel with phase-separated Particle Image and Tracking Velocimetry. Results from one of the first direct measurements of vortex-trapped sand grains under oscillatory flows are presented. The vortices mobilized sand grains along the ripple slopes just prior to flow reversal and transported the suspended sediment grains. During several flow cycles, some sand grains were temporarily trapped in the vortex, prescribing semi-circular trajectories off-center from the vortex core in quadrants of the vortex that were closest to the ripple slope, as illustrated by Nielsen (1992, https://doi.org/10.1142/1269 ). Comparisons of the horizontal sediment grain velocity with the horizontal fluid velocity yielded a linear relationship with a slope of 0.87. The vertical grain velocities also varied linearly with the vertical fluid velocity with a slope of approximately 1 and an offset of −0.08 m s −1 . The offset is close to the still water settling velocity for coarse sand grains, as hypothesized during vortex trapping. Additionally, estimates of the off-center distance, between the centers of the semi-circular sediment paths and vortex cores, compared well with the ratio of the settling velocity to the radian frequency of the vortex yielding a linear regression slope of 0.99. Improved understanding of vortex trapping effects on sediment dynamics may decrease uncertainty in model predictions of large-scale coastal hydrodynamics and sediment transport.

Journal of Geophysical Research: Earth Surface↗

USGS2005

Earthquake fracture energy inferred from kinematic rupture models on extended faults

We estimate fracture energy on extended faults for several recent earthquakes by retrieving dynamic traction evolution at each point on the fault plane from slip history imaged by inverting ground motion waveforms. We define the breakdown work (Wb) as the excess of work over some minimum traction level achieved during slip. Wb is equivalent to "seismological" fracture energy (G) in previous investigations. Our numerical approach uses slip velocity as a boundary condition on the fault. We employ a three-dimensional finite difference algorithm to compute the dynamic traction evolution in the time domain during the earthquake rupture. We estimate Wb by calculating the scalar product between dynamic traction and slip velocity vectors. This approach does not require specifying a constitutive law and assuming dynamic traction to be collinear with slip velocity. If these vectors are not collinear, the inferred breakdown work depends on the initial traction level. We show that breakdown work depends on the square of slip. The spatial distribution of breakdown work in a single earthquake is strongly correlated with the slip distribution. Breakdown work density and its integral over the fault, breakdown energy, scale with seismic moment according to a power law (with exponent 0.59 and 1.18, respectively). Our estimates of breakdown work range between 4 ?? 105 and 2 ?? 107 J/m2 for earthquakes having moment magnitudes between 5.6 and 7.2. We also compare our inferred values with geologic surface energies. This comparison might suggest that breakdown work for large earthquakes goes primarily into heat production. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Imaging the transition from Aleutian subduction to Yakutat collision in central Alaska, with local earthquakes and active source data

In southern and central Alaska the subduction and active volcanism of the Aleutian subduction zone give way to a broad plate boundary zone with mountain building and strike-slip faulting, where the Yakutat terrane joins the subducting Pacific plate. The interplay of these tectonic elements can be best understood by considering the entire region in three dimensions. We image three-dimensional seismic velocity using abundant local earthquakes, supplemented by active source data. Crustal low-velocity correlates with basins. The Denali fault zone is a dominant feature with a change in crustal thickness across the fault. A relatively high-velocity subducted slab and a low-velocity mantle wedge are observed, and high Vp/Vs beneath the active volcanic systems, which indicates focusing of partial melt. North of Cook Inlet, the subducted Yakutat slab is characterized by a thick low-velocity, high-Vp/Vs, crust. High-velocity material above the Yakutat slab may represent a residual older slab, which inhibits vertical flow of Yakutat subduction fluids. Alternate lateral flow allows Yakutat subduction fluids to contribute to Cook Inlet volcanism and the Wrangell volcanic field. The apparent northeast edge of the subducted Yakutat slab is southwest of the Wrangell volcanics, which have adakitic composition consistent with melting of this Yakutat slab edge. In the mantle, the Yakutat slab is subducting with the Pacific plate, while at shallower depths the Yakutat slab overthrusts the shallow Pacific plate along the Transition fault. This region of crustal doubling within the shallow slab is associated with extremely strong plate coupling and the primary asperity of the Mw 9.2 great 1964 earthquake. Copyright 2006 by the American Geophysical Union.

Alaska↗

USGS2004

Recalculated probability of M ≥ 7 earthquakes beneath the Sea of Marmara, Turkey

New earthquake probability calculations are made for the Sea of Marmara region and the city of Istanbul, providing a revised forecast and an evaluation of time-dependent interaction techniques. Calculations incorporate newly obtained bathymetric images of the North Anatolian fault beneath the Sea of Marmara [Le Pichon et al., 2001; Armijo et al., 2002]. Newly interpreted fault segmentation enables an improved regional A.D. 1500-2000 earthquake catalog and interevent model, which form the basis for time-dependent probability estimates. Calculations presented here also employ detailed models of coseismic and postseismic slip associated with the 17 August 1999 M = 7.4 Izmit earthquake to investigate effects of stress transfer on seismic hazard. Probability changes caused by the 1999 shock depend on Marmara Sea fault-stressing rates, which are calculated with a new finite element model. The combined 2004-2034 regional Poisson probability of M≥7 earthquakes is ~38%, the regional time-dependent probability is 44 ± 18%, and incorporation of stress transfer raises it to 53 ± 18%. The most important effect of adding time dependence and stress transfer to the calculations is an increase in the 30 year probability of a M ??? 7 earthquake affecting Istanbul. The 30 year Poisson probability at Istanbul is 21%, and the addition of time dependence and stress transfer raises it to 41 ± 14%. The ranges given on probability values are sensitivities of the calculations to input parameters determined by Monte Carlo analysis; 1000 calculations are made using parameters drawn at random from distributions. Sensitivities are large relative to mean probability values and enhancements caused by stress transfer, reflecting a poor understanding of large-earthquake aperiodicity.

Sea Of Marmara↗

USGS2012

Seismic reflection images of the crust of the northern part of the Chugach terrane, Alaska: Results of a survey for the Trans-Alaska Crustal Transect (TACT)

Deep crustal seismic reflection data show strong reflections from the middle and lower crust of the convergent continental margin near the eastern end of the Aleutian trench. These data were collected across the Border Ranges fault system, a major suture zone that separates the Peninsular and Chugach tectonostratigraphic terranes. The Chugach terrane consists of metamorphosed, strongly deformed flysch, subduction melange as well as blueschist- and greenschist-facies rocks. Even so, highly reflective rocks lie at great depth within or beneath this terrane. The shallow part of the seismic section (0–5 s) reveals only a few reflections, which undulate and interweave, indicating that some of them originate outside of the plane of the seismic section. Other undulating events could reveal early Cenozoic granitic plutons or antiformal stacks of rock imbricated along thrust faults. In contrast to this shallow reflection pattern, three reflection bands contain most of the events on the middle part of the seismic section (5–12 s or 12–34 km). We interpret divergent reflections within the uppermost band (5–5.5 s) as gently north dipping thrust faults, and subparallel reflections truncated at the boundaries of the upper and middle bands may image a metamorphic foliation or thrust faults. If the truncated reflections represent thrust faults, then the reflection geometry suggests that subparallel roof and floor thrust faults bound imbricated rocks, forming duplex structures. The upper and middle bands are separated by an area on the seismic section that reveals few reflections. This area correlates approximately in depth with a high velocity (7.5 km/s) layer evident in refraction models. The poorly reflective, high velocity unit might be volcanic rocks like those exposed in the structurally lowest part of the Chugach terrane. If so, then the bottom of the poorly reflective rocks indicates the base of the Chugach terrane. Rocks that cause the upper and middle reflection bands and the intervening, poorly reflective volcanic rocks may form a midcrustal shear zone that is about 10 km thick. The deepest band of subparallel reflections (30–35 km) correlates closely in depth with the top of the Wadati-Benioff zone associated with the underthrusting Pacific plate. The reflections probably stem from the decollement that separates lithospheric plates within the subduction zone. This reflection band increases abruptly in thickness, possibly indicating local underplating of subducted sediment. No events from the oceanic Moho have been recognized. Neither the Border Ranges fault system, the suture between the Chugach and Peninsular terranes, nor deep crustal layers of the Peninsular terrane are evident in these seismic reflection data.

Journal of Geophysical Research Solid Earth↗

USGS2003

Imaging the complexity of an active normal fault system: The 1997 Colfiorito (central Italy) case study

Six moderate magnitude earthquakes (5 < Mw < 6) ruptured normal fault segments of the southern sector of the North Apennine belt (central Italy) in the 1997 Colfiorito earthquake sequence. We study the progressive activation of adjacent and nearby parallel faults of this complex normal fault system using ???1650 earthquake locations obtained by applying a double-difference location method, using travel time picks and waveform cross-correlation measurements. The lateral extent of the fault segments range from 5 to 10 km and make up a broad, ???45 km long, NW trending fault system. The geometry of each segment is quite simple and consists of planar faults gently dipping toward SW with an average dip of 40??-45??. The fault planes are not listric but maintain a constant dip through the entire seismogenic volume, down to 8 km depth. We observe the activation of faults on the hanging wall and the absence of seismicity in the footwall of the structure. The observed fault segmentation appears to be due to the lateral heterogeneity of the upper crust: preexisting thrusts inherited from Neogene's compressional tectonic intersect the active normal faults and control their maximum length. The stress tensor obtained by inverting the six main shock focal mechanisms of the sequence is in agreement with the tectonic stress active in the inner chain of the Apennine, revealing a clear NE trending extension direction. Aftershock focal mechanisms show a consistent extensional kinematics, 70% of which are mechanically consistent with the main shock stress field.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

Volcanology: Lessons learned from Synthetic Aperture Radar imagery

Twenty years of continuous Earth observation by satellite SAR have resulted in numerous new insights into active volcanism, including a better understanding of subsurface magma storage and transport, deposition of volcanic materials on the surface, and the structure and development of volcanic edifices. This massive archive of data has resulted in fundamental leaps in our understanding of how volcanoes work – for example, identifying magma accumulation at supposedly quiescent volcanoes, even in remote areas or in the absence of ground-based data. In addition, global compilations of volcanic activity facilitate comparison of deformation behavior between different volcanic arcs and statistical evaluation of the strong link between deformation and eruption. SAR data are also increasingly used in timely hazard evaluation thanks to decreases in data latency and growth in processing and analysis techniques. The existing archive of SAR imagery is on the cusp of being enhanced by a new generation of satellite SAR missions, in addition to ground-based and airborne SAR systems, which will provide enhanced temporal and spatial resolution, broader geographic coverage, and improved availability of data to the scientific community. Now is therefore an opportune time to review the contributions of SAR imagery to volcano science, monitoring, and hazard mitigation, and to explore the future potential for SAR in volcanology. Provided that the ever-growing volume of SAR data can be managed effectively, we expect the future application of SAR data to expand from being a research tool for analyzing volcanic activity after the fact, to being a monitoring and research tool capable of imaging a wide variety of processes on different temporal and spatial scales as those processes are occurring. These data can then be used to develop new models of how volcanoes work and to improve quantitative forecasts of volcanic activity as a means of mitigating risk from future eruptions.

Journal of Volcanology and Geothermal Research↗

USGS2003

Fault trends on the seaward slope of the Aleutian Trench: Implications for a laterally changing stress field tied to a westward increase in oblique convergence

Normal faults along the seaward trench slope (STS) commonly strike parallel to the trench in response to bending of the oceanic plate into the subduction zone. This is not the circumstance for the Aleutian Trench, where the direction of convergence gradually changes westward, from normal to transform motion. GLORIA side-scan sonar images document that the Aleutian STS is dominated by faults striking oblique to the trench, west of 179??E and east of 172??W. These images also show a pattern of east-west trending seafloor faults that are aligned parallel to the spreading fabric defined by magnetic anomalies. The stress-strain field along the STS is divided into two domains west and east, respectively, of 179??E. Over the western domain, STS faults and nodal planes of earthquakes are oriented oblique (9??-46??) to the trench axis and (69??-90??) to the magnetic fabric. West of 179??E, STS fault strikes change by 36?? from the E-W trend of STS where the trench-parallel slip gets larger than its orthogonal component of convergence. This rotation indicates that horizontal stresses along the western domain of the STS are deflected by the increasing obliquity in convergence. An analytical model supports the idea that strikes of STS faults result from a superposition of stresses associated with the dextral shear couple of the oblique convergence and stresses caused by plate bending. For the eastern domain, most nodal planes of earthquakes strike parallel to the outer rise, indicating bending as the prevailing mechanism causing normal faulting. East of 172??W, STS faults strike parallel to the magnetic fabric but oblique (10??-26??) to the axis of the trench. On the basis of a Coulomb failure criterion the trench-oblique strikes probably result from reactivation of crustal faults generated by spreading. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2017

Magnetotelluric imaging of lower crustal melt and lithospheric hydration in the Rocky Mountain Front transition zone, Colorado, USA

We present an electrical resistivity model of the crust and upper mantle from two‐dimensional (2‐D) anisotropic inversion of magnetotelluric data collected along a 450 km transect of the Rio Grande rift, southern Rocky Mountains, and High Plains in Colorado, USA. Our model provides a window into the modern‐day lithosphere beneath the Rocky Mountain Front to depths in excess of 150 km. Two key features of the 2‐D resistivity model are (1) a broad zone (~200 km wide) of enhanced electrical conductivity (<20 Ωm) in the midcrust to lower crust that is centered beneath the highest elevations of the southern Rocky Mountains and (2) hydrated lithospheric mantle beneath the Great Plains with water content in excess of 100 ppm. We interpret the high conductivity region of the lower crust as a zone of partially molten basalt and associated deep‐crustal fluids that is the result of recent (less than 10 Ma) tectonic activity in the region. The recent supply of volatiles and/or heat to the base of the crust in the late Cenozoic implies that modern‐day tectonic activity in the western United States extends to at least the western margin of the Great Plains. The transition from conductive to resistive upper mantle is caused by a gradient in lithospheric modification, likely including hydration of nominally anhydrous minerals, with maximum hydration occurring beneath the Rocky Mountain Front. This lithospheric “hydration front” has implications for the tectonic evolution of the continental interior and the mechanisms by which water infiltrates the lithosphere.

Colorado↗

USGS2004

Transient volcano deformation sources imaged with interferometric synthetic aperture radar: Application to Seguam Island, Alaska

Thirty interferometric synthetic aperture radar (InSAR) images, spanning various intervals during 1992–2000, document coeruptive and posteruptive deformation of the 1992–1993 eruption on Seguam Island, Alaska. A procedure that combines standard damped least squares inverse methods and collective surfaces, identifies three dominant amorphous clusters of deformation point sources. Predictions generated from these three point source clusters account for both the spatial and temporal complexity of the deformation patterns of the InSAR data. Regularized time series of source strength attribute a distinctive transient behavior to each of the three source clusters. A model that combines magma influx, thermoelastic relaxation, poroelastic effects, and petrologic data accounts for the transient, interrelated behavior of the source clusters and the observed deformation. Basaltic magma pulses, which flow into a storage chamber residing in the lower crust, drive this deformational system. A portion of a magma pulse is injected into the upper crust and remains in storage during both coeruption and posteruption intervals. This injected magma degasses and the volatile products accumulate in a shallow poroelastic storage chamber. During the eruption, another portion of the magma pulse is transported directly to the surface via a conduit roughly centered beneath Pyre Peak on the west side of the island. A small amount of this magma remains in storage during the eruption, and posteruption thermoelastic contraction ensues. This model, made possible by the excellent spatial and temporal coverage of the InSAR data, reveals a relatively simple system of interrelated predictable processes driven by magma dynamics.

Journal of Geophysical Research B: Solid Earth↗

USGS1999

Ventifacts at the Pathfinder landing site

About half of the rocks at the Mars Pathfinder Ares Vallis landing site appear to be ventifacts, rocks abraded by windborne particles. Comparable resolution images taken by the Imager for Mars Pathfinder (IMP) camera and the Viking landers show that ventifacts are more abundant at the Pathfinder site. The ventifacts occur in several forms, including rocks with faceted edges, finger‐like projections, elongated pits, flutes, grooves, and possible rills. The trends of elongated pits, flutes, grooves, and rills cluster at ∼280–330° clockwise from north and generally dip 10–30° away from their trend direction. These orientations are indicative of southeast to northwest winds and differ from the trend of wind tails at the landing site, the direction of local wind streaks, and predictions of the Global Circulation Model, all of which indicate northeast to southwest winds. The disparity between these data sets strongly suggests that local circulation patterns have changed since the abrasion of the ventifacted rocks. The greater number of ventifacts at the Pathfinder site compared to either of the Viking sites is most easily explained as being due to a larger supply of abrading particles, composed of either sand‐sized grains or indurated dust aggregates, and higher surface roughness, which should increase the momentum of saltating grains. The Pathfinder ventifacts may have formed shortly after the deposition of outflow channel sediments nearly 2 Gry ago, when a large local supply of abrading particles should have been abundant and atmospheric conditions may have been more conducive to rock abrasion from saltating grains. Based on how ventifacts form on Earth, the several ventifact forms seen at the Pathfinder site and their presence on some rocks but not on others are probably due to local airflow conditions, original rock shape, exposure duration, rock movement, and to a lesser extent, rock lithology. The abundance of ventifacts at the Pathfinder site, together with other evidence of weathering, indicates that unaltered rock surfaces are rare on Mars.

Journal of Geophysical Research E: Planets↗

USGS2017

Terrestrial ecosystem model performance in simulating productivity and its vulnerability to climate change in the northern permafrost region

Realistic projection of future climate-carbon (C) cycle feedbacks requires better understanding and an improved representation of the C cycle in permafrost regions in the current generation of Earth system models. Here we evaluated 10 terrestrial ecosystem models for their estimates of net primary productivity (NPP) and responses to historical climate change in permafrost regions in the Northern Hemisphere. In comparison with the satellite estimate from the Moderate Resolution Imaging Spectroradiometer (MODIS; 246 ± 6 g C m −2 yr −1 ), most models produced higher NPP (309 ± 12 g C m −2 yr −1 ) over the permafrost region during 2000–2009. By comparing the simulated gross primary productivity (GPP) with a flux tower-based database, we found that although mean GPP among the models was only overestimated by 10% over 1982–2009, there was a twofold discrepancy among models (380 to 800 g C m −2 yr −1 ), which mainly resulted from differences in simulated maximum monthly GPP (GPP max ). Most models overestimated C use efficiency (CUE) as compared to observations at both regional and site levels. Further analysis shows that model variability of GPP and CUE are nonlinearly correlated to variability in specific leaf area and the maximum rate of carboxylation by the enzyme Rubisco at 25°C ( V c max_25 ), respectively. The models also varied in their sensitivities of NPP, GPP, and CUE to historical changes in climate and atmospheric CO 2 concentration. These results indicate that model predictive ability of the C cycle in permafrost regions can be improved by better representation of the processes controlling CUE and GPP max as well as their sensitivity to climate change.

Journal of Geophysical Research: Biogeosciences↗

USGS2005

Interferometric synthetic aperture radar study of Okmok volcano, Alaska, 1992-2003: Magma supply dynamics and postemplacement lava flow deformation

Okmok volcano, located in the central Aleutian arc, Alaska, is a dominantly basaltic complex topped with a 10-km-wide caldera that formed circa 2.05 ka. Okmok erupted several times during the 20th century, most recently in 1997; eruptions in 1945, 1958, and 1997 produced lava flows within the caldera. We used 80 interferometric synthetic aperture radar (InSAR) images (interferograms) to study transient deformation of the volcano before, during, and after the 1997 eruption. Point source models suggest that a magma reservoir at a depth of 3.2 km below sea level, located beneath the center of the caldera and about 5 km northeast of the 1997 vent, is responsible for observed volcano-wide deformation. The preeruption uplift rate decreased from about 10 cm yr −1 during 1992–1993 to 2 ∼ 3 cm yr −1 during 1993–1995 and then to about −1 ∼ −2 cm yr −1 during 1995–1996. The posteruption inflation rate generally decreased with time during 1997–2001, but increased significantly during 2001–2003. By the summer of 2003, 30 ∼ 60% of the magma volume lost from the reservoir in the 1997 eruption had been replenished. Interferograms for periods before the 1997 eruption indicate consistent subsidence of the surface of the 1958 lava flows, most likely due to thermal contraction. Interferograms for periods after the eruption suggest at least four distinct deformation processes: (1) volcano-wide inflation due to replenishment of the shallow magma reservoir, (2) subsidence of the 1997 lava flows, most likely due to thermal contraction, (3) deformation of the 1958 lava flows due to loading by the 1997 flows, and (4) continuing subsidence of 1958 lava flows buried beneath 1997 flows. Our results provide insights into the postemplacement behavior of lava flows and have cautionary implications for the interpretation of inflation patterns at active volcanoes.

Alaska↗

USGS2009

Constraining earthquake source inversions with GPS data: 2. A two-step approach to combine seismic and geodetic data sets

We present a new method to combine static and wavefield data to image earthquake ruptures. Our combined inversion is a two-step procedure, following the work of Hernandez et al. (1999), and takes into account the differences between the resolutions of the two data sets. The first step consists of an inversion of the static field, which yields a map of slip amplitude. This inversion exploits a special irregular grid that takes into account the resolution of the static data. The second step is an inversion of the radiated wavefield; it results in the determination of the time evolution of slip on the fault. In the second step, the slip amplitude is constrained to resemble the static slip amplitude map inferred from the GPS inversion. Using this combined inversion, we study the source process of the 2004 M6 Parkfield, California, earthquake. We conclude that slip occurred in two main regions of the fault, each of which displayed distinct rupture behaviors. Slip initiated at the hypocenter with a very strong bilateral burst of energy. Here, slip was localized in a narrow area approximately 10 km long, the rupture velocity was very fast (???3.5 km/s), and slip only lasted a short period of time (<1 s). Then the rupture proceeded to a wider region 12-20 km northwest of the hypocenter. Here, the earthquake developed in a more moderated way: the rupture velocity slowed to ???3.0 km/s and slip lasted longer (1-2 s). The maximum slip amplitude was 0.45 m. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2010

High‐resolution locations of triggered earthquakes and tomographic imaging of Kilauea Volcano's south flank

The spatiotemporal patterns of seismicity beneath Kilauea's south flank give insight to the structure and geometry of the decollement on which large, tsunamigenic earthquakes have occurred, and its relation to slow slip events (SSEs), which have been observed every 1 to 2 years since 1997. In order to record earthquakes triggered by a SSE that was predicted to occur in March 2007, a temporary network of 20 seismometers was deployed on Kilauea's south flank, termed the SEQ network. While the SSE did not occur until 17 June 2007, theSEQ network recorded over 3000 earthquakes, including those triggered by the SSE. We relocate hypocenters of volcano‐tectonic earthquakes and invert for P and S wave velocity structure using waveform cross‐correlation and double‐difference tomography using data from the SEQ network and the permanent Hawaii Volcano Observatory network (HVO) data, with additional data from other previous temporary arrays. The best‐constrained hypocenters, recorded by both the SEQ and HVO networks, indicate the decollement as a subhorizontal layer of seismicity at 8 km depth less than 1 km thick in most areas, with the western portion of the decollement dipping to the southeast. The seismicity triggered by the June 2007 SSE includes over 400 earthquakes overlapping with the southern edge of the decollement seismicity. A shallower swarm of earthquakes also occurred between 2 and 7 km depth in April 2007 near Apua Point, and may have been indirectly triggered by the M w 8.1 Solomon Islands earthquake at ∼6000 km distance, which occurred 48 h prior to the beginning of the swarm.

Hawaii↗

USGS1996

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↗

USGS1999

P wave crustal velocity structure in the greater Mount Rainier area from local earthquake tomography

We present results from a local earthquake tomographic imaging experiment in the greater Mount Rainier area. We inverted P wave arrival times from local earthquakes recorded at permanent and temporary Pacific Northwest Seismograph Network seismographs between 1980 and 1996. We used a method similar to that described by Lees and Crosson [1989], modified to incorporate the parameter separation method for decoupling the hypocenter and velocity problems. In the upper 7 km of the resulting model there is good correlation between velocity anomalies and surface geology. Many focal mechanisms within the St. Helens seismic zone have nodal planes parallel to the epicentral trend as well as to a north-south trending low-velocity trough, leading us to speculate that the trough represents a zone of structural weakness in which a moderate ( M 6.5–7.0) earthquake could occur. In contrast, the western Rainier seismic zone does not correlate in any simple way with anomaly patterns or focal mechanism fault planes, leading us to infer that it is less likely to experience a moderate earthquake. A ∼10 km-wide low-velocity anomaly occurs 5 to 18 km beneath the summit of Mount Rainier, which we interpret to be a signal of a region composed of hot, fractured rock with possible small amounts of melt or fluid. No systematic velocity pattern is observed in association with the southern Washington Cascades conductor. A midcrustal anomaly parallels the Olympic-Wallowa lineament as well as several other geophysical trends, indicating that it may play an important role in regional tectonics.

Journal of Geophysical Research B: Solid Earth↗

USGS2022

Flood lavas on Earth, Io and Mars

Flood lavas are major geological features on all the major rocky planetary bodies. They provide important insight into the dynamics and chemistry of the interior of these bodies. On the Earth, they appear to be associated with major and mass extinction events. It is therefore not surprising that there has been significant research on flood lavas in recent years. Initial models suggested eruption durations of days and volumetric fluxes of order 10 7 m 3 s −1 with flows moving as turbulent floods. However, our understanding of how lava flows can be emplaced under an insulating crust was revolutionized by the observations of actively inflating pahoehoe flows in Hawaii. These new ideas led to the hypothesis that flood lavas were emplaced over many years with eruption rates of the order of 10 4 m 3 s −1 . The field evidence indicates that flood lava flows in the Columbia River Basalts, Deccan Traps, Etendeka lavas, and the Kerguelen Plateau were emplaced as inflated pahoehoe sheet flows. This was reinforced by the observation of active lava flows of ≥100 km length on Io being formed as tube-fed flows fed by moderate eruption rates (10 2 –10 3 m 3 s −1 ). More recently it has been found that some flood lavas are also emplaced in a more rapid manner. New high-resolution images from Mars revealed ‘platy–ridged’ flood lava flows, named after the large rafted plates and ridges formed by compression of the flow top. A search for appropriate terrestrial analogues found an excellent example in Iceland: the 1783–1784 Laki Flow Field. The brecciated Laki flow top consists of pieces of pahoehoe, not aa clinker, leading us to call this ‘rubbly pahoehoe’. Similar flows have been found in the Columbia River Basalts and the Kerguelen Plateau. We hypothesize that these flows form with a thick, insulating, but mobile crust, which is disrupted when surges in the erupted flux are too large to maintain the normal pahoehoe mode of emplacement. Flood lavas emplaced in this manner could have (intermittently) reached effusion rates of the order of 10 6 m 3 s −1 .

Journal of the Geological Society↗