Geology Reports⌕ Search

SEARCH · Geology Reports

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 307 records · Page 17Linked to original sources
USGS2012

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗
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↗
USGS2012

Hydrous carbonates on Mars?: Evidence from Mariner 6/7 infrared spectrometer and ground‐based telescopic spectra

Absorption features at 2.28 and 5.4 μm identified in Mariner 6/7 infrared spectrometer and terrestrial telescopic spectra are consistent with the spectra of hydrous magnesium carbonates such as hydromagnesite and artinite. Spectral characteristics of these hydrous carbonates are different from those of the anhydrous carbonates, as the former do not have the strong spectral features typically associated with anhydrous carbonates such as calcite and siderite. Theoretical mixing indicates that, depending on the type of hydrous carbonate, 10–20 wt % can be incorporated into the regolith without contradicting the spectral observations or the Viking x ray fluorescence chemical analysis. Hydrous carbonates form as weathering products of mafic minerals in the presence of H 2 O and CO 2 , even in the Antarctic. Their formation as evaporite minerals from either original magmas or hydrothermally altered rocks is consistent with the Martian environment, provided liquid water is or has been at least transiently present. On Earth, formation of hydrous Mg carbonates is associated with the production of amorphous iron oxides, which is consistent with both the environment and the inferred surface mineralogy of Mars. These minerals are about 60 wt % H 2 O, CO 3 , and OH; if they are abundant everywhere at the 10% level, then about 6% of the surface weight could be volatiles bound in this type of mineral. Although the stability of hydrous carbonates in a Martian environment is uncertain, there may be kinetic factors inhibiting the dehydration of these minerals, which may persist metastably in the current environment. Although the spectroscopic evidence for anhydrous carbonates is scant, the possible presence of hydrous carbonates provides an appealing mechanism for the existence of carbonates on Mars.

Journal of Geophysical Research↗
USGS2009

The 2006-2007 Kuril Islands great earthquake sequence

[1] The southwestern half of a ∼500 km long seismic gap in the central Kuril Island arc subduction zone experienced two great earthquakes with extensive preshock and aftershock sequences in late 2006 to early 2007. The nature of seismic coupling in the gap had been uncertain due to the limited historical record of prior large events and the presence of distinctive upper plate, trench and outer rise structures relative to adjacent regions along the arc that have experienced repeated great interplate earthquakes in the last few centuries. The intraplate region seaward of the seismic gap had several shallow compressional events during the preceding decades (notably an M S 7.2 event on 16 March 1963), leading to speculation that the interplate fault was seismically coupled. This issue was partly resolved by failure of the shallow portion of the interplate megathrust in an M W = 8.3 thrust event on 15 November 2006. This event ruptured ∼250 km along the seismic gap, just northeast of the great 1963 Kuril Island (M w = 8.5) earthquake rupture zone. Within minutes of the thrust event, intense earthquake activity commenced beneath the outer wall of the trench seaward of the interplate rupture, with the larger events having normal-faulting mechanisms. An unusual double band of interplate and intraplate aftershocks developed. On 13 January 2007, an M W = 8.1 extensional earthquake ruptured within the Pacific plate beneath the seaward edge of the Kuril trench. This event is the third largest normal-faulting earthquake seaward of a subduction zone on record, and its rupture zone extended to at least 33 km depth and paralleled most of the length of the 2006 rupture. The 13 January 2007 event produced stronger shaking in Japan than the larger thrust event, as a consequence of higher short-period energy radiation from the source. The great event aftershock sequences were dominated by the expected faulting geometries; thrust faulting for the 2006 rupture zone, and normal faulting for the 2007 rupture zone. A large intraplate compressional event occurred on 15 January 2009 (M w = 7.4) near 45 km depth, below the rupture zone of the 2007 event and in the vicinity of the 16 March 1963 compressional event. The fault geometry, rupture process and slip distributions of the two great events are estimated using very broadband teleseismic body and surface wave observations. The occurrence of the thrust event in the shallowest portion of the interplate fault in a region with a paucity of large thrust events at greater depths suggests that the event removed most of the slip deficit on this portion of the interplate fault. This great earthquake doublet demonstrates the heightened seismic hazard posed by induced intraplate faulting following large interplate thrust events. Future seismic failure of the remainder of the seismic gap appears viable, with the northeastern region that has also experienced compressional activity seaward of the megathrust warranting particular attention.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

Volcano deformation source parameters estimated from InSAR: Sensitivities to uncertainties in seismic tomography

The eruption cycle of a volcano is controlled in part by the upward migration of magma. The characteristics of the magma flux produce a deformation signature at the Earth's surface. Inverse analyses use geodetic data to estimate strategic controlling parameters that describe the position and pressurization of a magma chamber at depth. The specific distribution of material properties controls how observed surface deformation translates to source parameter estimates. Seismic tomography models describe the spatial distributions of material properties that are necessary for accurate models of volcano deformation. This study investigates how uncertainties in seismic tomography models propagate into variations in the estimates of volcano deformation source parameters inverted from geodetic data. We conduct finite element model-based nonlinear inverse analyses of interferometric synthetic aperture radar (InSAR) data for Okmok volcano, Alaska, as an example. We then analyze the estimated parameters and their uncertainties to characterize the magma chamber. Analyses are performed separately for models simulating a pressurized chamber embedded in a homogeneous domain as well as for a domain having a heterogeneous distribution of material properties according to seismic tomography. The estimated depth of the source is sensitive to the distribution of material properties. The estimated depths for the homogeneous and heterogeneous domains are 2666 ± 42 and 3527 ± 56 m below mean sea level, respectively (99% confidence). A Monte Carlo analysis indicates that uncertainties of the seismic tomography cannot account for this discrepancy at the 99% confidence level. Accounting for the spatial distribution of elastic properties according to seismic tomography significantly improves the fit of the deformation model predictions and significantly influences estimates for parameters that describe the location of a pressurized magma chamber.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

CRUST 5.1: A global crustal model at 5° x 5°

We present a new global model for the Earth's crust based on seismic refraction data published in the period 1948-1995 and a detailed compilation of ice and sediment thickness. An extensive compilation of seismic refraction measurements has been used to determine the crustal structure on continents and their margins. Oceanic crust is modeled with both a standard model for normal oceanic crust, and variants for nonstandard regions, such as oceanic plateaus. Our model (CRUST 5.1) consists of 2592 5° × 5° tiles in which the crust and uppermost mantle are described by eight layers: (1) ice, (2) water, (3) soft sediments, (4) hard sediments, (5) crystalline upper, (6) middle, (7) lower crust, and (8) uppermost mantle. Topography and bathymetry are adopted from a Mandard database (ETOPO-5). Compressional wave velocity in each layer is based on field measurements, and shear wave velocity and density are estimated using recently published empirical V p -V s and V p -density relationships. The crustal model differs from previous models in that (1) the thickness and seismic/density structure of sedimentary basins is accounted for more completely, (2) the velocity structure of unmeasured regions is estimated using statistical averages that are based on a significantly larger database of crustal structure, (3) the compressional wave, shear wave, and density structure have been explicitly specified using newly available constraints from field and laboratory studies. Thus - this global crustal model is based on substantially more data than previous models and differs from them in many important respects. A new map of the thickness of the Earth's crust is presented, and we illustrate the application of this model by using it to provide the crustal correction for surface wave phase velocity maps. Love waves at 40 s are dominantly sensitive to crustal structure, and there is a very close correspondence between observed phase velocities at this period and those predicted by CRUST 5.1. We find that the application of crustal corrections to long-period (167 s) Rayleigh waves significantly increases the variance in the phase velocity maps and strengthens the upper mantle , velocity anomalies beneath stable continental regions. A simple calculation of crustal isostacy indicates significant lateral variations in upper mantle density. The model CRUST 5.1 provides a complete description of the physical properties of the Earth's crust at a scale of 5° × 5° and can be used for a wide range of seismological and nonseismological problems.

Journal of Geophysical Research B: Solid Earth↗
USGS2022

Microgravity change during the 2008-2018 Kı̄lauea summit eruption: Nearly a decade of subsurface mass accumulation

Results from nine microgravity campaigns from Kı̄lauea, Hawaiʻi, spanning most of the volcano's 2008–2018 summit eruption, indicate persistent mass accumulation at shallow levels. A weighted least squares approach is used to recover microgravity results from a network of benchmarks around Kı̄lauea's summit, eliminate instrumental drift, and restore suspected data tares. A total mass of 1.9 × 10 11 kg was determined from these microgravity campaigns to have accumulated below Kı̄lauea Caldera during 2009–2015 at an estimated depth of 1.3 km below sea level. Only a fraction of this mass is reflected in surface deformation, and this is consistent with previously reported discrepancies between subsurface mass accumulation and observed surface deformation. The discrepancy, amongst other independent evidence from gas emissions, seismicity, and continuous gravimetry, indicate densification of magma in the reservoirs below the volcano summit. This densification may have been driven by degassing through the summit vent. It is hypothesized that during the final years of the summit eruption, magma densification resulted in a buildup of pressure in the reservoirs that may have contributed to the lower East Rift Zone outbreak of 2018. The observed mass accumulation beneath Kı̄lauea could not have been detected through other techniques and illustrates the importance of microgravity measurements in volcano monitoring.

Hawai'i↗
USGS1998

Poroelastic rebound along the Landers 1992 earthquake surface rupture

Maps of surface displacement following the 1992 Landers, California, earthquake, generated by interferometric processing of ERS-1 synthetic aperture radar (SAR) images, reveal effects of various postseismic deformation processes along the 1992 surface rupture. The large-scale pattern of the postseismic displacement field includes large lobes, mostly visible on the west side of the fault, comparable in shape with the lobes observed in the coseismic displacement field. This pattern and the steep displacement gradient observed near the Emerson-Camp Rock fault cannot be simply explained by afterslip on deep sections of the 1992 rupture. Models show that horizontal slip occurring on a buried dislocation in a Poisson's material produces a characteristic quadripole pattern in the surface displacement field with several centimeters of vertical motion at distances of 10–20 km from the fault, yet this pattern is not observed in the postseismic interferograms. As previously proposed to explain local strain in the fault step overs [ Peltzer et al ., 1996b], we argue that poroelastic rebound caused by pore fluid flow may also occur over greater distances from the fault, compensating the vertical ground shift produced by fault afterslip. Such a rebound is explained by the gradual change of the crustal rocks' Poisson's ratio value from undrained (coseismic) to drained (postseismic) conditions as pore pressure gradients produced by the earthquake dissipate. Using the Poisson's ratio values of 0.27 and 0.31 for the drained and undrained crustal rocks, respectively, elastic dislocation models show that the combined contributions of afterslip on deep sections of the fault and poroelastic rebound can account for the range change observed in the SAR data and the horizontal displacement measured at Global Positioning System (GPS) sites along a 60-km-long transect across the Emerson fault [ Savage and Svarc , 1997]. Using a detailed surface slip distribution on the Homestead Valley, Kickapoo, and Johnson Valley faults, we modeled the poroelastic rebound in the Homestead Valley pull apart. A Poisson's ratio value of 0.35 for the undrained gouge rocks in the fault zone is required to account for the observed surface uplift in the 3.5 years following the earthquake. This large value implies a seismic velocity ratio V p /V s of 2.1, consistent with the observed low V s values of fault zone guided waves at shallow depth [ Li et al ., 1997]. The SAR data also reveal postseismic creep along shallow patches of the Eureka Peak and Burnt Mountain faults with a characteristic decay time of 0.8 years. Coseismic, dilatant hardening (locking process) followed by post-seismic, pore pressure controlled fault creep provide a plausible mechanism to account for the decay time of the observed slip rate along this section of the fault.

Journal of Geophysical Research B: Solid Earth↗
USGS2019

Coseismic slip and early afterslip of the M6.0 August 24, 2014 South Napa, California, earthquake

We employ strong motion seismograms and static offsets from the Global Positioning System, Interferometric Synthetic Aperture Radar, and other measurements in order to derive a coseismic slip and afterslip model of the M6.0 24 August 2014 South Napa earthquake. This earthquake ruptured an ∼13‐km‐long portion of the West Napa fault with predominantly right‐lateral strike slip. In the kinematic seismic slip inversions, we couple the coseismic slip and afterslip distributions by requiring both distributions to involve right‐lateral strike slip with positive amplitude, with the net static slip being the sum of the two. We consider several candidate fault geometries: a first involving two steeply east dipping fault planes that reach Earth's surface at the western surface trace (STW), where most surface rupture was observed, a second involving a steeply west dipping plane that also reaches Earth's surface at the STW, and a third involving a combination of two variably west dipping planes constrained to pass through the locus of postseismic seismicity located ∼1 km west of the STW. The data are best fit using the model of two east dipping fault planes, with coseismic slip up to ∼1.2 m on a dominant shallow asperity about 10 km north of the hypocenter and on deeper asperities on the southern part of the rupture. Afterslip up to 1 m is concentrated along the southern part of the rupture at depths 5 km, consistent with surface observations of afterslip. Seismic moments associated with coseismic slip and afterslip are 1.13×10 18 N m (Mw 6.00) and 3.64×10 17 N m, respectively.

California↗
USGS2016

Coseismic slip and early afterslip of the 2015 Illapel, Chile, earthquake: Implications for frictional heterogeneity and coastal uplift

Great subduction earthquakes are thought to rupture portions of the megathrust, where interseismic coupling is high and velocity-weakening frictional behavior is dominant, releasing elastic deformation accrued over a seismic cycle. Conversely, postseismic afterslip is assumed to occur primarily in regions of velocity-strengthening frictional characteristics that may correlate with lower interseismic coupling. However, it remains unclear if fixed frictional properties of the subduction interface, coseismic or aftershock-induced stress redistribution, or other factors control the spatial distribution of afterslip. Here we use interferometric synthetic aperture radar and Global Position System observations to map the distribution of coseismic slip of the 2015 M w 8.3 Illapel, Chile, earthquake and afterslip within the first 38 days following the earthquake. We find that afterslip overlaps the coseismic slip area and propagates along-strike into regions of both high and moderate interseismic coupling. The significance of these observations, however, is tempered by the limited resolution of geodetic inversions for both slip and coupling. Additional afterslip imaged deeper on the fault surface bounds a discrete region of deep coseismic slip, and both contribute to net uplift of the Chilean Coastal Cordillera. A simple partitioning of the subduction interface into regions of fixed frictional properties cannot reconcile our geodetic observations. Instead, stress heterogeneities, either preexisting or induced by the earthquake, likely provide the primary control on the afterslip distribution for this subduction zone earthquake. We also explore the occurrence of coseismic and postseismic coastal uplift in this sequence and its implications for recent hypotheses concerning the source of permanent coastal uplift along subduction zones.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Evidence for a shallow pluton beneath the Goat Rocks Wilderness, Washington, from gravity and magnetic data

A gravity and aeromagnetic study was conducted in and adjacent to the Goat Rocks Wilderness area, Washington. This work was done in conjunction with geologic and geochemical studies which were critically important to our interpretation. The Bouguer gravity anomaly map displays a large central positive gravity anomaly and three superimposed smaller positive anomalies. The main anomaly, about 26 by 11 km with an amplitude of about 16 mGal seems to delineate a dense large pluton or sheeted dike complex intruding less dense, older volcanic and sedimentary rocks lying directly below the volcano. The pluton or sheeted dike complex may have fed the Goat Rocks volcano throughout its 1–2 m.y. history. Three small feeders, mapped on the surface as andesite intrusions, may be the cause of three smaller gravity anomalies superimposed on the main anomaly. The magnetic anomaly maps also exhibit a large central anomaly whose source is the Goat Rocks pluton. The magnetic data suggest that the three shallow intrusive cupolas have undergone a significant amount of hydrothermal alteration. In addition, to some extent the magnetic data reflect the terrain of the area. Topographic highs produce magnetic highs, and topographic lows produce magnetic lows. The deviations from this pattern delineate subsurface geology. Magnetic lows may be of particular interest in locating zones of mineralization where magnetic minerals have been altered and have become significantly less magnetic.

Washington↗
USGS2012

Magnetic signals from the core of the earth and secular variation

An oscillating, radial magnetic dipole source was assumed to exist in the core of the earth, 100 km beneath the core-mantle boundary. As an approximation, electromagnetic propagation was assumed in the core in lieu of hydromagnetic propagation, which could not be used because of unknown internal fields.Using Debye potentials, the radial and horizontal components of the surface fields were calculated using various assumed conductivity parameters in the core and in the mantle. It is concluded that most spherical harmonic models of the earth's magnetic field do not include enough terms to describe properly the field of core sources with periods of tens of years. These short-period variations are especially important in describing the secular variation. Because of this, a proper description of the secular variation requires more spherical harmonic terms than are required for the field itself. Inadequate representation of short-period variations in spherical harmonic models may contribute to the rapid deterioration of predictive models. Alternatives to spherical harmonic analysis for secular variation should be investigated; regardless of the method used, a much greater spatial distribution of high quality secular variation data is needed.

Journal of Geophysical Research Solid Earth↗
USGS1996

Cyclic stressing and seismicity at strongly coupled subduction zones

We use the finite element method to analyze stress variations in and near a strongly coupled subduction zone during an earthquake cycle. Deformation is assumed to be uniform along strike (plane strain on a cross section normal to the trench axis), and periodic earthquake slip is imposed consistent with the long-term rate of plate convergence and degree of coupling. Simulations of stress and displacement rate fields represent periodic fluctuations in time superimposed on an average field. The oceanic plate, descending slab, and continental lithosphere are assumed here to respond elastically to these fluctuations, and the remaining mantle under and between plates is assumed to respond as Maxwell viscoelastic. In the first part of the analysis we find that computed stress fluctuations in space and time are generally consistent with observed earthquake mechanism variations with time since a great thrust event. In particular, trench-normal extensional earthquakes tend to occur early in the earthquake cycle toward the outer rise but occur more abundantly late in the cycle in the subducting slab downdip of the main thrust zone. Compressional earthquakes, when they occur at all, have the opposite pattern. Our results suggest also that the actual timing of extensional outer rise events is controlled by the rheology of the shallow aseismic portion of the thrust interface. The second part of the analysis shows the effects of mantle relaxation on the rate of ground surface deformation during the earthquake cycle. Models without relaxation predict a strong overall compressional strain rate in the continental plate above the main thrust zone, with the strain rate constant between mainshocks. However with significant relaxation present, a localized region of unusually low compressional, or even slightly extensional, strain rate develops along the surface of the continental plate above and somewhat inland from the downdip edge of the locked main thrust zone. The low strain rate starts in the middle or late part of the cycle, depending on position. This result suggests that the negligible or small contraction measured on the Shumagin Islands, Alaska, during 1980 to 1991, may not invalidate an interpretation of that region as being a moderately coupled subduction zone. In contrast, mantle relaxation causes only modest temporal nonuniformity of uplift rates in the overriding plate and of extensional stress rates in the subducting plate, even when the Maxwell time is an order of magnitude less than the recurrence interval.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Eolian features in the Western Desert of Egypt and some applications to Mars

Relations of landform types to wind regimes, bedrock composition, sediment supply, and topography are shown by field studies and satellite photographs of the Western Desert of Egypt. This desert, which lies at the core of the largest hyperarid region on earth, provides analogs of Martian wind-formed features. These include sand dunes, alternating light and dark streaks, knob ‘shadows,’ and yardangs. Surface particles have been segregated by wind into deposits (dunes, sand sheets, and light streaks) that can be differentiated by their grain size distributions, surface shapes, and colors. Throughgoing sand of mostly fine to medium grain size is migrating southward in longitudinal dune belts and barchan chains whose long axes lie parallel to the prevailing northerly winds, but topographic variations such as scarps and depressions strongly influence the zones of deposition and dune morphology. Sand from the longitudinal dunes on the plains is commonly redistributed into barchans in the depressions. These barchans are generally simple crescents that are morphologically similar to many of the dunes seen on Viking orbiter pictures of the north polar sand sea on mars. Light streaks are depositional features consisting of dune belts and elongate sheets of coarse to medium sand and granules. Intervening dark streaks are erosional features consisting of strips of desert-varnished bedrock and lag gravel surfaces exposed between the sand deposits. The shape of both light and dark streaks is controlled by wind flow around topographic highs. Dark zones (shadows) in the lee of mountains, hills, and knobs are erosional products from the topographic highs; they change shape only in response to movement of the adjacent lighter-colored sand deposits. Streamlined yardangs carved in crystalline limestone constitute one of the largest yardang fields on earth. Yardangs occur also in sandstone of the Nubian Series and in lacustrine sediments. The variables that affect the patterns of wind erosion and deposition in the Western Desert are topographic effects on wind velocities and directions, resistance of the bedrock, sand supply, and climatic change with time; vegetation is essentially absent and is not a controlling factor.

Journal of Geophysical Research Solid Earth↗
USGS2017

Temporal variation of tectonic tremor activity in southern Taiwan around the 2010 ML6.4 Jiashian Earthquake

Deep tectonic tremor, which is extremely sensitive to small stress variations, could be used to monitor fault-zone processes during large earthquake cycles and aseismic processes before large earthquakes. In this study, we develop an algorithm for the automatic detection and location of tectonic tremor beneath the southern Central Range of Taiwan and examine the spatio-temporal relationship between tremor and the 4 March 2010 ML6.4 Jiashian earthquake, located about 20 km from active tremor sources. We find that tremor in this region has a relatively short duration, short recurrence time, and no consistent correlation with surface GPS data. We find a short-term increase in the tremor rate 19 days before the Jiashian mainshock, and around the time when the tremor rate began to rise, one GPS station recorded a flip in its direction of motion. We hypothesize that tremor is driven by a slow-slip event that preceded the occurrence of the shallower Jiashian mainshock, even though the inferred slip is too small to be observed by all GPS stations. Our study shows that tectonic tremor may reflect stress variation during the pre-nucleation process of a nearby earthquake.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Television observations from Surveyor 3

A total of 6315 pictures were taken by the television camera on Surveyor 3 after the lunar landing. These pictures have provided much new information about the location of the landing site on the moon, the detailed topographic and geologic characteristics of the lunar surface, and the appearance of the earth as seen from the moon, both during eclipse of the sun and during partial direct illumination by the sun. Surveyor 3 landed in a subdued crater slightly more than 200 meters in diameter, which has a low rounded rim and is about 15 meters deep. The spacecraft is situated on the east wall of the crater, about half way between the center of the crater and the rim crest. The spacecraft is inclined 14.7° ± 1.0° toward the west. The selenographic coordinates of the landing site are 2.94°S latitude, 23.34°W longitude, relative to selenodetic control adopted by the Aeronautical Chart and Information Center. Small morphologic elements of the landing site include small craters, linear ridges and troughs, and fragmental debris. The craters and fragmental debris resemble those observed at the Surveyor 1 landing site, both in distribution of shape and in distribution of size. Most of the craters in Surveyor 1 and 3 pictures are inferred to be of impact origin. Their size-frequency distribution corresponds to the distribution that would be produced by repetitive bombardment by meteoroids, a bombardment sufficiently prolonged that the crater population has reached a steady state or has come to equilibrium. Some of the craters observed at the Surveyor 3 landing site are inferred to be of secondary impact origin, and some probably have been formed either by subsidence or by drainage of fragmental debris into cracks or fissures in the subsurface. Fragmental debris at the landing site is inferred to have been derived primarily by the same process of repetitive bombardment that produced the majority of craters. The inferred volumetric size-frequency distribution of fragments, derived from the observed size distribution of fragments on the surface, is similar to the distribution that would be produced by repetitive bombardment of coherent rocks by meteoroids with a mass-frequency distribution like that found from observed meteors and recovered meteorites on earth. Two prominent strewn fields of blocky debris were observed around two craters, 13 and 15 meters across, at the Surveyor 3 site. The 13-meter crater has a sharp raised rim, and the 15-meter crater has a more subdued rounded rim. The blocks associated with the subdued crater have twice as high a mean roundness as the blocks associated with the raised-rim crater, and they are much more deeply buried. The size-frequency distribution function for the fragments in each of the strewn fields of blocks resembles the size-frequency distribution for fragments ejected by impacts in strong rock, such as Meteor Crater, Arizona. Most of the fragments at the surface of the Surveyor 3 landing site are evidently part of a layer of fragmental material of low cohesion that is at least 1 meter thick along the upper parts of the wall of the crater in which Surveyor 3 landed and may be much thicker near the center of the crater. The evidence from the Surveyor 3 pictures suggests that this layer of debris, or regolith, is subject to downslope creep or mass movement. Creep is probably caused by seismic shaking, due mainly to near and far impact events and perhaps due partly to internal lunar seismicity. Disturbances of the lunar surface produced by Surveyor 3, like those produced by Surveyor 1, exposed material at depths of a few centimeters or less that was darker than the material at the surface. The albedo of the fine-grained fragmental debris is probably 20 to 30% lower at depths of only a fraction of a millimeter than it is at the optically observed surface. All coarse fragments protruding above the general level of the surface have a higher albedo than the fine-grained matrix of the surface. These general photometric relationships can be explained if it is assumed that the surfaces of the particles in the shallow lunar subsurface tend to become coated with a dark substance; the term proposed for this hypothetical substance is ‘lunar varnish’ On the protruding surface of blocks and coarse fragments the lunar varnish is scrubbed off by the processes that cause rounding. The exposed surfaces of fine particles on the lunar surface are similarly affected, but, because they are mixed relatively rapidly with particles just beneath the surface, the process is incomplete and the fine-grained material exposed at the lunar surface, therefore, has a lower albedo than blocks and other large fragments. Coating of particles by lunar varnish evidently takes place just beneath the surface. The estimated normal luminance factor (normal albedo) of an undisturbed part of the lunar surface next to footpad 2 of the Surveyor 3 spacecraft is 8.5%. An area of the lunar surface disturbed by the surface sampler has an estimated normal luminance factor of 6.6%, and fine-grained material placed on footpad 2 by the surface sampler has an estimated normal luminance factor of 7.6%. The errors in all these estimates may be as high as 25% because of uncertainties of correction required for light scattered from the camera mirror. Preliminary search for color differences, by color reconstitution methods, revealed no determinable differences in color among various coarse blocks, the fine-grained matrix of the surface, or fine-grained material disturbed by the surface sampler. Surveyor 3 pictures of the eclipse of the sun by the earth revealed a bright region in the refraction halo surrounding the earth, which was correlated with the position of the sun, and a series of bright beads that occurred over regions of the earth largely clear of clouds. Clouds tend to occult the refracted rays of the sun, most of which pass through the lower part of the atmosphere at the limb; the beads occurred in the depressions in the optical silhouette of the earth. Preliminary reduction of the color of the refracted light showed that the brightest region, near the position of the sun, exhibited a correlated color temperature close to 4800°K. The color temperature tended to be lower for light that followed paths of greater atmospheric absorption. Preliminary analysis of Surveyor 3 pictures of the partly illuminated earth revealed colors similar to the colors recorded from orbit by the Mercury and Gemini astronauts.

Journal of Geophysical Research↗
USGS2008

Wind-driven particle mobility on Mars: Insights from Mars Exploration Rover observations at "El Dorado" and surroundings at Gusev Crater

The ripple field known as 'El Dorado' was a unique stop on Spirit's traverse where dust-raising, active mafic sand ripples and larger inactive coarse-grained ripples interact, illuminating several long-standing issues of Martian dust mobility, sand mobility, and the origin of transverse aeolian ridges. Strong regional wind events endured by Spirit caused perceptible migration of ripple crests in deposits SSE of El Dorado, erasure of tracks in sandy areas, and changes to dust mantling the site. Localized thermal vortices swept across El Dorado, leaving paths of reduced dust but without perceptibly damaging nearly cohesionless sandy ripple crests. From orbit, winds responsible for frequently raising clay-sized dust into the atmosphere do not seem to significantly affect dunes composed of (more easily entrained) sand-sized particles, a long-standing paradox. This disparity between dust mobilization and sand mobilization on Mars is due largely to two factors: (1) dust occurs on the surface as fragile, low-density, sand-sized aggregates that are easily entrained and disrupted, compared with clay-sized air fall particles; and (2) induration of regolith is pervasive. Light-toned bed forms investigated at Gusev are coarse-grained ripples, an interpretation we propose for many of the smallest linear, light-toned bed forms of uncertain origin seen in high-resolution orbital images across Mars. On Earth, wind can organize bimodal or poorly sorted loose sediment into coarse-grained ripples. Coarse-grained ripples could be relatively common on Mars because development of durable, well-sorted sediments analogous to terrestrial aeolian quartz sand deposits is restricted by the lack of free quartz and limited hydraulic sediment processing. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2018

Using stereo satellite imagery to account for ablation, entrainment, and compaction in volume calculations for rock avalanches on Glaciers: Application to the 2016 Lamplugh Rock Avalanche in Glacier Bay National Park, Alaska

The use of preevent and postevent digital elevation models (DEMs) to estimate the volume of rock avalanches on glaciers is complicated by ablation of ice before and after the rock avalanche, scour of material during rock avalanche emplacement, and postevent ablation and compaction of the rock avalanche deposit. We present a model to account for these processes in volume estimates of rock avalanches on glaciers. We applied our model by calculating the volume of the 28 June 2016 Lamplugh rock avalanche in Glacier Bay National Park, Alaska. We derived preevent and postevent 2‐m resolution DEMs from WorldView satellite stereo imagery. Using data from DEM differencing, we reconstructed the rock avalanche and adjacent surfaces at the time of occurrence by accounting for elevation changes due to ablation and scour of the ice surface, and postevent deposit changes. We accounted for uncertainties in our DEMs through precise coregistration and an assessment of relative elevation accuracy in bedrock control areas. The rock avalanche initially displaced 51.7 ± 1.5 Mm3 of intact rock and then scoured and entrained 13.2 ± 2.2 Mm3 of snow and ice during emplacement. We calculated the total deposit volume to be 69.9 ± 7.9 Mm3. Volume estimates that did not account for topographic changes due to ablation, scour, and compaction underestimated the deposit volume by 31.0–46.8 Mm3. Our model provides an improved framework for estimating uncertainties affecting rock avalanche volume measurements in glacial environments. These improvements can contribute to advances in the understanding of rock avalanche hazards and dynamics.

Alaska↗