USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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

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

At least 73 records · Page 4Linked to original sources

Triggered tremor sweet spots in Alaska

To better understand what controls fault slip along plate boundaries, we have exploited the abundance of seismic and geodetic data available from the richly varied tectonic environments composing Alaska. A search for tremor triggered by 11 large earthquakes throughout all of seismically monitored Alaska reveals two tremor “sweet spots”—regions where large-amplitude seismic waves repeatedly triggered tremor between 2006 and 2012. The two sweet spots locate in very different tectonic environments—one just trenchward and between the Aleutian islands of Unalaska and Akutan and the other in central mainland Alaska. The Unalaska/Akutan spot corroborates previous evidence that the region is ripe for tremor, perhaps because it is located where plate-interface frictional properties transition between stick-slip and stably sliding in both the dip direction and laterally. The mainland sweet spot coincides with a region of complex and uncertain plate interactions, and where no slow slip events or major crustal faults have been noted previously. Analyses showed that larger triggering wave amplitudes, and perhaps lower frequencies (<~0.03 Hz), may enhance the probability of triggering tremor. However, neither the maximum amplitude in the time domain or in a particular frequency band, nor the geometric relationship of the wavefield to the tremor source faults alone ensures a high probability of triggering. Triggered tremor at the two sweet spots also does not occur during slow slip events visually detectable in GPS data, although slow slip below the detection threshold may have facilitated tremor triggering.

Alaska

Aftershocks of the 2010 M w 7.2 El Mayor-Cucapah earthquake revealcomplex faulting in the Yuha Desert, California

We detect and precisely locate over 9500 aftershocks that occurred in the Yuha Desert region during a 2 month period following the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake. Events are relocated using a series of absolute and relative relocation procedures that include Hypoinverse, Velest, and hypoDD. Location errors are reduced to ~40 m horizontally and ~120 m vertically.Aftershock locations reveal a complex pattern of faulting with en echelon fault segments trending toward the northwest, approximately parallel to the North American-Pacific plate boundary and en echelon, conjugate features trending to the northeast. The relocated seismicity is highly correlated with published surface mapping of faults that experienced triggered surface slip in response to the EMC main shock. Aftershocks occurred between 2 km and 11 km depths, consistent with previous studies of seismogenic thickness in the region. Three-dimensional analysis reveals individual and intersecting fault planes that are limited in their along-strike length. These fault planes remain distinct structures at depth, indicative of conjugate faulting, and do not appear to coalesce onto a throughgoing fault segment. We observe a complex spatiotemporal migration of aftershocks, with seismicity that jumps between individual fault segments that are active for only a few days to weeks. Aftershock rates are roughly consistent with the expected earthquake production rates of Dieterich (1994). The conjugate pattern of faulting and nonuniform aftershock migration patterns suggest that strain in the Yuha Desert is being accommodated in a complex manner.

California

Maximum magnitude earthquakes induced by fluid injection

Analysis of numerous case histories of earthquake sequences induced by fluid injection at depth reveals that the maximum magnitude appears to be limited according to the total volume of fluid injected. Similarly, the maximum seismic moment seems to have an upper bound proportional to the total volume of injected fluid. Activities involving fluid injection include (1) hydraulic fracturing of shale formations or coal seams to extract gas and oil, (2) disposal of wastewater from these gas and oil activities by injection into deep aquifers, and (3) the development of enhanced geothermal systems by injecting water into hot, low-permeability rock. Of these three operations, wastewater disposal is observed to be associated with the largest earthquakes, with maximum magnitudes sometimes exceeding 5. To estimate the maximum earthquake that could be induced by a given fluid injection project, the rock mass is assumed to be fully saturated, brittle, to respond to injection with a sequence of earthquakes localized to the region weakened by the pore pressure increase of the injection operation and to have a Gutenberg-Richter magnitude distribution with a b value of 1. If these assumptions correctly describe the circumstances of the largest earthquake, then the maximum seismic moment is limited to the volume of injected liquid times the modulus of rigidity. Observations from the available case histories of earthquakes induced by fluid injection are consistent with this bound on seismic moment. In view of the uncertainties in this analysis, however, this should not be regarded as an absolute physical limit.

Journal of Geophysical Research B: Solid Earth

Observations of static Coulomb stress triggering of the November 2011 M5.7 Oklahoma earthquake sequence

In November 2011, a M 5.0 earthquake occurred less than a day before a M 5.7 earthquake near Prague, Oklahoma, which may have promoted failure of the mainshock and thousands of aftershocks along the Wilzetta fault, including a M 5.0 aftershock. The M 5.0 foreshock occurred in close proximity to active fluid injection wells; fluid injection can cause a buildup of pore fluid pressure, decrease the fault strength, and may induce earthquakes. Keranen et al . [ 2013 ] links the M 5.0 foreshock with fluid injection, but the relationship between the foreshock and successive events has not been investigated. Here we examine the role of coseismic Coulomb stress transfer on earthquakes that follow the M 5.0 foreshock, including the M 5.7 mainshock. We resolve the static Coulomb stress change onto the focal mechanism nodal plane that is most consistent with the rupture geometry of the three M ≥ 5.0 earthquakes, as well as specified receiver fault planes that reflect the regional stress orientation. We find that Coulomb stress is increased, e.g., fault failure is promoted, on the nodal planes of ~60% of the events that have focal mechanism solutions, and more specifically, that the M 5.0 foreshock promoted failure on the rupture plane of the M 5.7 mainshock. We test our results over a range of effective coefficient of friction values. Hence, we argue that the M 5.0 foreshock, induced by fluid injection, potentially triggered a cascading failure of earthquakes along the complex Wilzetta fault system.

Oklahoma

Real-time inversions for finite fault slip models and rupture geometry based on high-rate GPS data

We present an inversion strategy capable of using real-time high-rate GPS data to simultaneously solve for a distributed slip model and fault geometry in real time as a rupture unfolds. We employ Bayesian inference to find the optimal fault geometry and the distribution of possible slip models for that geometry using a simple analytical solution. By adopting an analytical Bayesian approach, we can solve this complex inversion problem (including calculating the uncertainties on our results) in real time. Furthermore, since the joint inversion for distributed slip and fault geometry can be computed in real time, the time required to obtain a source model of the earthquake does not depend on the computational cost. Instead, the time required is controlled by the duration of the rupture and the time required for information to propagate from the source to the receivers. We apply our modeling approach, called Bayesian Evidence-based Fault Orientation and Real-time Earthquake Slip, to the 2011 Tohoku-oki earthquake, 2003 Tokachi-oki earthquake, and a simulated Hayward fault earthquake. In all three cases, the inversion recovers the magnitude, spatial distribution of slip, and fault geometry in real time. Since our inversion relies on static offsets estimated from real-time high-rate GPS data, we also present performance tests of various approaches to estimating quasi-static offsets in real time. We find that the raw high-rate time series are the best data to use for determining the moment magnitude of the event, but slightly smoothing the raw time series helps stabilize the inversion for fault geometry.

Journal of Geophysical Research

Crustal earthquake triggering by pre-historic great earthquakes on subduction zone thrusts

Triggering of earthquakes on upper plate faults during and shortly after recent great ( M> 8.0) subduction thrust earthquakes raises concerns about earthquake triggering following Cascadia subduction zone earthquakes. Of particular regard to Cascadia was the previously noted, but only qualitatively identified, clustering of M> ~6.5 crustal earthquakes in the Puget Sound region between about 1200–900 cal yr B.P. and the possibility that this was triggered by a great Cascadia thrust subduction thrust earthquake, and therefore portends future such clusters. We confirm quantitatively the extraordinary nature of the Puget Sound region crustal earthquake clustering between 1200–900 cal yr B.P., at least over the last 16,000. We conclude that this cluster was not triggered by the penultimate, and possibly full-margin, great Cascadia subduction thrust earthquake. However, we also show that the paleoseismic record for Cascadia is consistent with conclusions of our companion study of the global modern record outside Cascadia, that M> 8.6 subduction thrust events have a high probability of triggering at least one or more M> ~6.5 crustal earthquakes.

Cascadia;Puget Sound

Volcanic tremor masks its seismogenic source: Results from a study of noneruptive tremor recorded at Mount St. Helens, Washington

On 2 October 2004, a significant noneruptive tremor episode occurred during the buildup to the 2004–2008 eruption of Mount St. Helens (Washington). This episode was remarkable both because no explosion followed, and because seismicity abruptly stopped following the episode. This sequence motivated us to consider a model for volcanic tremor that does not involve energetic gas release from magma but does involve movement of conduit magma through extension on its way toward the surface. We found that the tremor signal was composed entirely of Love and Rayleigh waves and that its spectral bandwidth increased and decreased with signal amplitude, with broader bandwidth signals containing both higher and lower frequencies. Our modeling results demonstrate that the forces giving rise to this tremor were largely normal to conduit walls, generating hybrid head waves along conduit walls that are coupled to internally reflected waves. Together these form a crucial part of conduit resonance, giving tremor wavefields that are largely a function of waveguide geometry and velocity. We find that the mechanism of tremor generation fundamentally masks the nature of the seismogenic source giving rise to resonance. Thus multiple models can be invoked to explain volcanic tremor, requiring that information from other sources (such as visual observations, geodesy, geology, and gas geochemistry) be used to constrain source models. With concurrent GPS and field data supporting rapid rise of magma, we infer that tremor resulted from drag of nearly solid magma along rough conduit walls as magma was forced toward the surface.

Washington

A continuous record of intereruption velocity change at Mount St. Helens from coda wave interferometry

In September 2004, Mount St. Helens volcano erupted after nearly 18&thinsp;years of quiescence. However, it is unclear from the limited geophysical observations when or if the magma chamber replenished following the 1980&ndash;1986 eruptions in the years before the 2004&ndash;2008 extrusive eruption. We use coda wave interferometry with repeating earthquakes to measure small changes in the velocity structure of Mount St. Helens volcano that might indicate magmatic intrusion. By combining observations of relative velocity changes from many closely located earthquake sources, we solve for a continuous function of velocity changes with time. We find that seasonal effects dominate the relative velocity changes. Seismicity rates and repeating earthquake occurrence also vary seasonally; therefore, velocity changes and seismicity are likely modulated by snow loading, fluid saturation, and/or changes in groundwater level. We estimate hydrologic effects impart stress changes on the order of tens of kilopascals within the upper 4&thinsp;km, resulting in annual velocity variations of 0.5 to 1%. The largest nonseasonal change is a decrease in velocity at the time of the deep M w &thinsp;=&thinsp;6.8 Nisqually earthquake. We find no systematic velocity changes during the most likely times of intrusions, consistent with a lack of observable surface deformation. We conclude that if replenishing intrusions occurred, they did not alter seismic velocities where this technique is sensitive due to either their small size or the finite compressibility of the magma chamber. We interpret the observed velocity changes and shallow seasonal seismicity as a response to small stress changes in a shallow, pressurized system.

Journal of Geophysical Research B: Solid Earth

The 17 May 2012 M4.8 earthquake near Timpson, East Texas: An event possibly triggered by fluid injection

This study summarizes our investigation of the 17 May 2012 MW-RMT4.8 earthquake near Timpson, Texas, the largest earthquake recorded historically in eastern Texas. To identify preshocks and aftershocks of the 17 May event we examined the arrivals recorded at Nacogdoches (NATX) 30 km from the 17 May epicenter, at nearby USArray Transportable Array stations, and at eight temporary stations deployed between 26 May 2012 and mid-2013. At NATX we identified seven preshocks, the earliest occurring in April 2008. Reliably located aftershocks recorded by the temporary stations lie along a 6 km long NW-SE linear trend corresponding to a previously mapped basement fault that extends across the highest-intensity (MMI VII) region of the 17 May main shock. Earthquakes in this sequence are relatively shallow—with focal depths ranging from 1.6 to 4.6 km. Evidence supporting these depths include: hypocentral locations of exceptionally well-recorded aftershocks, S-P intervals at the nearest stations, and comparisons of synthetics and observed seismograms. Within 3 km of the linear trend of aftershock activity there are two Class II injection disposal wells injecting at 1.9 km depth beginning in August 2006 and February 2007, with injection rates averaging 42,750 m3/mo and 15,600 m3/mo, respectively. Several observations support the hypothesis that fluid injection triggered the Timpson sequence: well-located epicenters are situated near a mapped basement fault and near high-volume injection wells, focal depths are at or below the depths of injection, and the earliest preshock (April 2008) occurred after the onset of injection in 2006.

Texas

Three-dimensional seismic velocity structure of Mauna Loa and Kilauea volcanoes in Hawaii from local seismic tomography

We present a new three-dimensional seismic velocity model of the crustal and upper mantle structure for Mauna Loa and Kilauea volcanoes in Hawaii. Our model is derived from the first-arrival times of the compressional and shear waves from about 53,000 events on and near the Island of Hawaii between 1992 and 2009 recorded by the Hawaiian Volcano Observatory stations. The V p model generally agrees with previous studies, showing high-velocity anomalies near the calderas and rift zones and low-velocity anomalies in the fault systems. The most significant difference from previous models is in V p / V s structure. The high- V p and high- V p / V s anomalies below Mauna Loa caldera are interpreted as mafic magmatic cumulates. The observed low- V p and high- V p / V s bodies in the Kaoiki seismic zone between 5 and 15 km depth are attributed to the underlying volcaniclastic sediments. The high- V p and moderate- to low- V p / V s anomalies beneath Kilauea caldera can be explained by a combination of different mafic compositions, likely to be olivine-rich gabbro and dunite. The systematically low- V p and low- V p / V s bodies in the southeast flank of Kilauea may be caused by the presence of volatiles. Another difference between this study and previous ones is the improved V p model resolution in deeper layers, owing to the inclusion of events with large epicentral distances. The new velocity model is used to relocate the seismicity of Mauna Loa and Kilauea for improved absolute locations and ultimately to develop a high-precision earthquake catalog using waveform cross-correlation data.

Hawaii

Suspended particulate layers and internal waves over the southern Monterey Bay continental shelf: an important control on shelf mud belts?

Physical and optical measurements taken over the mud belt on the southern continental shelf of Monterey Bay, California documented the frequent occurrence of suspended particulate matter features, the majority of which were detached from the seafloor, centered 9–33 m above the bed. In fall 2011, an automated profiling mooring and fixed instrumentation, including a thermistor chain and upward-looking acoustic Doppler current profiler, were deployed at 70 m depth for 5 weeks, and from 12 to 16 October a long-range autonomous underwater vehicle performed across-shelf transects. Individual SPM events were uncorrelated with local bed shear stress caused by surface waves and bottom currents. Nearly half of all observed SPM layers occurred during 1 week of the study, 9–16 October 2011, and were advected past the fixed profiling mooring by the onshore phase of semidiurnal internal tide bottom currents. At the start of the 9–16 October period, we observed intense near-bed vertical velocities capable of lifting particulates into the middle of the water column. This “updraft” event appears to have been associated with nonlinear adjustment of high-amplitude internal tides over the mid and outer shelf. These findings suggest that nonlinear internal tidal motions can erode material over the outer shelf and that, once suspended, this SPM can then be transported shoreward to the middle and shallow sections of the mud belt. This represents a fundamental broadening of our understanding of how shelf mud belts may be built up and sustained.

California

Interactions between waves, sediment, and turbulence on a shallow estuarine mudflat

Measurements were collected on a shallow estuarine mudflat in northern San Francisco Bay to examine the physical processes controlling waves, turbulence, sediment resuspension, and their interactions. Tides alone forced weak to moderate currents of 10–30 cm s -1 in depths of 0–3 m, and maintained a background suspension of 30–50 mg L21 of fine sediment. In the presence of wind waves, bottom orbital velocities spanned 20–30 cm s -1 , suspended-sediment concentrations (SSC) at 15 and 30 cm above the bed (cmab) increased by 1–2 orders of magnitude, and vertical gradients in SSC were strong enough to produce turbulence-limiting stratification, with gradient Richardson numbers exceeding 0.25. Simultaneously, turbulent stresses (decomposed from wave motions) increased by an order of magnitude. The apparent contradiction of energetic turbulence in the presence of strong stratification was reconciled by considering the turbulent kinetic energy (TKE) budget: in general, dissipation and buoyancy flux were balanced by local shear production, and each of these terms increased during wave events. The classic wave-current boundary layer model represented the observations qualitatively, but not quantitatively since the velocity profile could not be approximated as logarithmic. Rather, the mean shear was elevated by the Stokes drift return flow and wind-generated surface stress, which diffused sediment upward and limited stratification. Our findings highlight a pathway for waves to supply energy to both the production and destruction of turbulence, and demonstrate that in such shallow depths, TKE and SSC can be elevated over more of the water column than predicted by traditional models.

California

Igneous mineralogy at Bradbury Rise: the first ChemCam campaign at Gale crater

Textural and compositional analyses using ChemCam Remote Micro Imager (RMI) and Laser Induced Breakdown Spectroscopy (LIBS) have been performed on 5 float rocks and coarse gravels along the first 100 meters of the Curiosity traverse at Bradbury rise. ChemCam, the first LIBS instrument sent to another planet, offers the opportunity to assess mineralogic diversity at grain-size scales (∼100μm) and, from this, lithologic diversity. Depth profiling indicates that targets are relatively free of surface coatings. One type of igneous rock is volcanic and includes both aphanitic (Coronation) and porphyritic (Mara) samples. The porphyritic sample shows dark grains that are likely pyroxene megacrysts in a finegrained mesostasis containing andesine needles. Both types have magnesium-poor basaltic compositions and in this respect are similar to the evolved Jake-Matijevic rock [Stolper et al., 2013] analyzed further along the Curiosity traverse both with APXS and ChemCam instruments. The second rock-type encountered is a coarse-grained intrusive rock (Thor Lake) showing equigranular texture with mm size crystals of feldspars and Fe-Ti oxides. Such a rock is not unique at Gale as the surrounding coarse gravels (such as Beaulieu) and the conglomerate Link are dominated by feldspathic (andesine-bytownite) clasts. Finally, alkali feldspar compositions associated with a silica polymorph have been analyzed in fractured filling material of Preble rock and in Stark, a putative pumice or an impact melt. These observations document magmatic diversity at Gale and describe the first fragments of feldspar-rich lithologies (possibly an anorthosite) that may be ancient crust transported from the crater rim and now forming float rocks, coarse gravel or conglomerate clasts.

Journal of Geophysical Research E: Planets

Trace element geochemistry (Li, Ba, Sr, and Rb) using Curiosity's ChemCam: early results for Gale crater from Bradbury Landing Site to Rocknest

The ChemCam instrument package on the Mars rover, Curiosity, provides new capabilities to probe the abundances of certain trace elements in the rocks and soils on Mars using the laser-induced breakdown spectroscopy technique. We focus on detecting and quantifying Li, Ba, Rb, and Sr in targets analyzed during the first 100 sols, from Bradbury Landing Site to Rocknest. Univariate peak area models and multivariate partial least squares models are presented. Li, detected for the first time directly on Mars, is generally low (<15 ppm). The lack of soil enrichment in Li, which is highly fluid mobile, is consistent with limited influx of subsurface waters contributing to the upper soils. Localized enrichments of up to ~60 ppm Li have been observed in several rocks but the host mineral for Li is unclear. Bathurst_Inlet is a fine-grained bedrock unit in which several analysis locations show a decrease in Li and other alkalis with depth, which may imply that the unit has undergone low-level aqueous alteration that has preferentially drawn the alkalis to the surface. Ba (~1000 ppm) was detected in a buried pebble in the Akaitcho sand ripple and it appears to correlate with Si, Al, Na, and K, indicating a possible feldspathic composition. Rb and Sr are in the conglomerate Link at abundances >100 ppm and >1000 ppm, respectively. These analysis locations tend to have high Si and alkali abundances, consistent with a feldspar composition. Together, these trace element observations provide possible evidence of magma differentiation and aqueous alteration.

Journal of Geophysical Research E: Planets

Sands at Gusev Crater, Mars

Processes, environments, and the energy associated with the transport and deposition of sand at Gusev Crater are characterized at the microscopic scale through the comparison of statistical moments for particle size and shape distributions. Bivariate and factor analyses define distinct textural groups at 51 sites along the traverse completed by the Spirit rover as it crossed the plains and went into the Columbia Hills. Fine-to-medium sand is ubiquitous in ripples and wind drifts. Most distributions show excess fine material, consistent with a predominance of wind erosion over the last 3.8 billion years. Negative skewness at West Valley is explained by the removal of fine sand during active erosion, or alternatively, by excess accumulation of coarse sand from a local source. The coarse to very coarse sand particles of ripple armors in the basaltic plains have a unique combination of size and shape. Their distribution display significant changes in their statistical moments within the ~400 m that separate the Columbia Memorial Station from Bonneville Crater. Results are consistent with aeolian and/or impact deposition, while the elongated and rounded shape of the grains forming the ripples, as well as their direction of origin, could point to Ma'adim Vallis as a possible source. For smaller particles on the traverse, our findings confirm that aeolian processes have dominated over impact and other processes to produce sands with the observed size and shape patterns across a spectrum of geologic (e.g., ripples and plains soils) and aerographic settings (e.g., wind shadows).

Journal of Geophysical Research E: Planets

Volcanic glass signatures in spectroscopic survey of newly proposed lunar pyroclastic deposits

Moon Mineralogy Mapper spectroscopic observations are used to assess the mineralogy of five sites that have recently been proposed to include lunar dark mantle deposits (DMDs). Volcanic glasses have, for the first time, clearly been identified at the location of three of the proposed pyroclastic deposits. This is the first time that volcanic glasses have been identified at such a small scale on the lunar surface from remote sensing observations. Deposits at Birt E, Schluter, and Walther A appear to be glassy DMDs. Deposits at Birt E and Schluter show (1) morphological evidence suggesting a likely vent and (2) mineralogical evidence indicative of the presence of volcanic glasses. The Walther A deposits, although they show no morphological evidence of vents, have the spectroscopic characteristics diagnostic of volcanic glasses. The deposits of the Freundlich-Sharonov basin are separated in two areas: (1) the Buys-Ballot deposits lack mineralogical and morphological evidence and thus are found to be associated with mare volcanism not with DMDs and (2) the Anderson crater deposits, which do not exhibit glassy DMD signatures, but they appear to be associated with possible vent structures and so may be classifiable as DMDs. Finally, dark deposits near the crater Kopff are found to be associated with likely mare volcanism and not associated with DMDs. The spectral identification of volcanic glass seen in many of the potential DMDs is a strong indicator of their pyroclastic origin.

Journal of Geophysical Research E: Planets

The rock abrasion record at Gale Crater: Mars Science Laboratory results from Bradbury Landing to Rocknest

Ventifacts, rocks abraded by wind-borne particles, are found in Gale Crater, Mars. In the eastward drive from &ldquo;Bradbury Landing&rdquo; to &ldquo;Rocknest,&rdquo; they account for about half of the float and outcrop seen by Curiosity's cameras. Many are faceted and exhibit abrasion textures found at a range of scales, from submillimeter lineations to centimeter-scale facets, scallops, flutes, and grooves. The drive path geometry in the first 100 sols of the mission emphasized the identification of abrasion facets and textures formed by westerly flow. This upwind direction is inconsistent with predictions based on models and the orientation of regional dunes, suggesting that these ventifact features formed from very rare high-speed winds. The absence of active sand and evidence for deflation in the area indicates that most of the ventifacts are fossil features experiencing little abrasion today.

Journal of Geophysical Research E: Planets

Chemistry and texture of the rocks at Rocknest, Gale Crater: Evidence for sedimentary origin and diagenetic alteration

A suite of eight rocks analyzed by the Curiosity Rover while it was stopped at the Rocknest sand ripple shows the greatest chemical divergence of any potentially sedimentary rocks analyzed in the early part of the mission. Relative to average Martian soil and to the stratigraphically lower units encountered as part of the Yellowknife Bay formation, these rocks are significantly depleted in MgO, with a mean of 1.3 wt %, and high in Fe, averaging over 20 wt % FeOT, with values between 15 and 26&thinsp;wt % FeOT. The variable iron and low magnesium and rock texture make it unlikely that these are igneous rocks. Rock surface textures range from rough to smooth, can be pitted or grooved, and show various degrees of wind erosion. Some rocks display poorly defined layering while others seem to show possible fractures. Narrow vertical voids are present in Rocknest 3, one of the rocks showing the strongest layering. Rocks in the vicinity of Rocknest may have undergone some diagenesis similar to other rocks in the Yellowknife Bay Formation as indicated by the presence of soluble calcium phases. The most reasonable scenario is that fine-grained sediments, potentially a mixture of feldspar-rich rocks from Bradbury Rise and normal Martian soil, were lithified together by an iron-rich cement.

Journal of Geophysical Research E: Planets