Search USGS⌕ Search

SEARCH · Search USGS

Results for “JGR Solid Earth”

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.

310 records · Page 8Linked to original sources
USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗
USGS2026

Dynamic drainage reorganization in Eastern Tibet: Insights from the Yangtze River first bend

The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10 5 –10 6 -year, rather than remaining in a long-term stationary configuration on ∼10 7 -year timescales.

eastern Tibetan Plateau, first bend of the Yangtze↗
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗
USGS2026

Multi-site evaluation of a postfire debris-flow runout forecast method

Postfire debris flows pose a hazard to human life, property, and infrastructure when they travel from steep source areas to urbanized alluvial fans or other developed areas. Existing methods for rapid (<1 week) postfire debris-flow hazard assessment document the increase in the likelihood and size of debris flows as the magnitude of high-intensity rain necessary to initiate debris flows increases but do not indicate the extent of downstream debris-flow runout. Although many models for the simulation of debris-flow motion are available, there is no established approach for using these models to delineate locations susceptible to postfire debris-flow runout that (a) is feasible to use at the spatial scale of an entire fire; (b) is appropriate for runout onto unconfined areas; (c) reproduces observed relations between runout and rainfall intensity; and (d) characterizes inherent uncertainty in runout, even without spatiotemporally variable rainfall. We propose and evaluate a method for generating postfire debris-flow runout hazard maps that has all the above qualities. Selection of case studies prioritized events triggered by a range of rainfall intensities, locations within and outside of southern California, and observed runout onto unconfined topography. Qualitative and quantitative assessment of performance for four events indicate that simulation results broadly match observations albeit with some discrepancies at a scale larger than structure or land parcel level (approximately 20-m by 20-m). The method may be used to identify potentially hazardous areas immediately following a fire and to provide approximate runout forecasts when a storm is imminent.

Arizona, California↗
USGS2025

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska↗
USGS2025

USGS addresses needs for lithium calibration and quality control materials for pLIBS analysis

Lithium (Li) is a globally important commodity used for energy storage, national defense, human health, and advanced technologies. Lithium resource development requires identifying deposits with elevated concentrations and optimal mineralogy, typically associated with select clays and pegmatites. Lithium is a light, highly reactive alkali metal with low atomic mass that is difficult to detect and quantify using conventional portable geochemical techniques such as X-ray fluorescence (XRF). However, portable laser-induced breakdown spectroscopy (pLIBS) is a powerful analytical technique for lithium exploration due to its ability to analyze solids quickly with minimal preparation. The expanded utility of pLIBS is hampered by the lack of matrix-matched calibration and quality control (QC) materials. The United States Geological Survey (USGS) has developed in-house lithium calibration and QC materials for lithium in clay and pegmatite matrices to address this limitation. We present the workflow and implementation of a custom-built matrix specific calibration on a SciAps Z-300 pLIBS, using proprietary Profile Builder software. The implementation of the custom calibration and quality control standards enables us to collect semiquantitative results directly from the pLIBS while in the field. Ultimately, this calibration has improved confidence in sample selection and collection in the field, providing more efficient site characterization.

Conference Paper↗
USGS2026

Investigating the influence of climate and volcanic surface aging on fluvial erosion: A case study of Réunion Island, Indian Ocean

Precipitation is one of the dominant drivers of landscape erosion and evolution; however, the effects of typical rainfall compared with less frequent, high-magnitude precipitation events on erosion remain unclear. Volcanic islands are ideal locations to study such phenomena due to their simple geometries, nontectonic construction, and strong spatiotemporal rainfall gradients. However, spatial variation in surface age, created during their construction, often complicates their degradation histories by introducing temporal changes in erosion rates as drainage networks develop. Réunion Island (western Indian Ocean) presents a clear example of this, with an east–west gradient in both surface age and mean annual precipitation, as well as infrequent cyclones that alter the background rainfall pattern. In this study, we analyze the effects of surface age, average rainfall, and rainfall variability on basin development and fluvial erosion across the island. We calculate basin-averaged values of basin morphology, age, precipitation, river discharge, eroded volumes, and erosion rates, and use these to analyze the dominant drivers of landscape evolution through a series of correlation analyses. Our results indicate a temporal dependence on the influence of precipitation, with young surfaces being dominantly eroded by high-rainfall events and older surfaces eroded by mean annual rainfall patterns. Furthermore, we show that drainage development of shield volcanoes follows similar trends to other volcano types, and suggest that surface permeability and groundwater structure are important controls on runoff-driven erosion on shield volcanoes. These results add new components to the question of how precipitation impacts erosion.

Réunion Island↗
USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗
USGS2026

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York↗
USGS2025

Relationship of atmospheric nitrogen deposition to soil nitrogen cycling along an elevation gradient in the Colorado Front Range

Microbial processing of atmospheric nitrogen (N) deposition regulates the retention and mobilization of N in soils, with important implications for water quality. Understanding the links between N deposition, microbial communities, N transformations, and water quality is critical as N deposition shifts toward reduced N and remains persistently high in many regions. Here, we investigated these connections along an elevation transect in the Colorado Front Range. Although rates of N deposition and pools of extractable N increased down the elevation transect, soil microbial communities and N transformation rates did not follow clear elevational patterns. The subalpine microbial community was distinct, corresponding to a high C:N ratio and low pH, while the microbial communities at the lower elevation sites were all very similar. Net nitrification, mineralization, and nitrification potential rates were highest at the Plains (1,700 m) and Montane (2,527 m) sites, suggesting that these ecosystems mobilize N. In contrast, the net immobilization of N observed at the Foothills (1,978 m) and Subalpine (3,015 m) sites suggests that these ecosystems retain N deposition. The contrast in N transformation rates between the plains and foothills, both of which receive elevated N deposition, may be due to spatial heterogeneity not captured in this study and warrants further investigation. Stream N concentrations from the subalpine to the foothills were consistently low, indicating that these soils are currently able to process and retain N deposition, but this may be disrupted if drought, wildfire, or land-use change alter the ability of the soils to retain N.

Colorado↗
USGS2024

Seismicity zoning at Coso geothermal field and stress changes from fluid production and migration

The Coso geothermal field is a major geothermal power production site in the western United States. It has been observed that low-magnitude seismic events ( M < 3.71) are unevenly distributed in three distinct zones, namely, nearfield (<3 km), midfield (3–6 km), and farfield (> 6 km) from the Coso geothermal plant. These zones exhibit distinct changes in earthquake location before and during geothermal production episodes that began in 1986. After 1986, the midfield region of the main flank experiences a significantly lower seismicity rate than the surrounding areas before production episodes. During 2014–2019, the farfield earthquakes cluster in the eastern and western parts of the greater Coso area, which is discernably different from how those pre-production earthquake events were distributed along the conjugate NW-SE and SW-NW trending structures across the main flank. Here, we analyze the stage of stress with finite-element-based poroelastic simulations to illustrate how the spatiotemporal evolution of the seismicity is associated with the pattern of stress perturbations caused by fluid migration amid the operations of geothermal power plants. Generally, ∼70% of co-production seismicity is found in zones of increased Coulomb stress between 2014 and 2019 at >99% confidence. Meanwhile, the midfield zone of seismic paucity overlaps with the zone of decreasing pore-fluid pressure. Overall, the results provide a physical explanation of how decadal geothermal operations at Coso have perturbed stress-field changes and contributed to the evolving characteristic seismic pattern, shedding insights into assessing the seismic hazard in other geothermal settings.

California↗
USGS2025

New constraints on location and timing of the Great Lakes tectonic zone, central Upper Peninsula, Michigan, USA

The Great Lakes tectonic zone (GLTZ) forms the boundary between the Wawa–Abitibi and Minnesota River Valley subprovinces within the Archean Superior Province. The GLTZ is concealed for all of its 1100 km length, except for a segment in the central Upper Peninsula of Michigan. There, it is exposed as a northwest-striking mylonite zone along a 11 km segment, extending to the onlap of Paleozoic rocks to the east. Farther east, its location has been unknown. Here, we use aeromagnetic and gravity data to develop interpretations of the expression of the GLTZ and to define its extent under cover. Aeromagnetic gradients over the mylonite zone are interpreted to be produced by structurally juxtaposed rocks with varying magnetizations. Gravity data show a regional gradient along the GLTZ, produced by the juxtaposition of a dense greenstone belt on the north against lower-density gneisses and granites on the south. The GLTZ is interpreted to extend ∼55 km under cover to the east. The GLTZ is terminated on the east by the buried eastern arm of the ca. 1100 Ma Midcontinent Rift. An undeformed granitic dike that cuts the mylonitic foliation produces a U–Pb apatite age of 2523 ± 33 Ma, implying no major post-Archean shearing occurred, and is at odds with previous interpretations of major Proterozoic reactivation. A granite intrusion in the Minnesota River Valley subprovince produces a Pb–Pb zircon age of 2606.9 ± 3.6/7.4 Ma. This suggests that magmatism related to the Sacred Heart orogeny, previously known in Minnesota, extended to Michigan.

Michigan↗
USGS2026

Geophysical architecture and geochronology of the Neoarchean Mentor anorthosite intrusive complex, northwestern Minnesota: Largest anorthosite complex of the Superior Province?

The Neoarchean Mentor anorthosite intrusive complex (MAIC) lies within the Wawa subprovince in northwestern Minnesota, in a region where the Wawa, Quetico, and Wabigoon subprovinces are juxtaposed in close proximity. Archean rocks are entirely concealed, and interpretations are developed from aeromagnetic, gravity, and borehole samples. The MAIC includes both anorthosite (dense, weakly magnetized) and oxide-rich gabbro (strongly magnetized) lithologies. Anorthosite is coarse-grained to megacrystic, intensely altered, and locally brecciated. Pervasive epidote alteration enhanced the density of the anorthosite via introduction of higher-density mineral assemblages, explaining why the MAIC produces a significant gravity high. Oxide-rich gabbro forms a border phase of the MAIC and has potential for vanadium, chromium, and titanium mineralization, and produces a strong aeromagnetic high. The MAIC is interpreted to extend over an area of 640 km 2 , making it the largest known anorthosite complex of the Superior Province, as measured by preserved areal extent. Modeling indicates the MAIC extends more than 6 km into the subsurface. A new Pb–Pb zircon age of 2737.2 ± 4.5 Ma is interpretated as the crystallization age of anorthosite within the MAIC, showing that the MAIC formed well before the ca. 2690 Ma Shebandowanian orogeny, and raising new questions about correlations with other parts of the Wawa subprovince. A low-density batholith, here informally called the Fertile batholith, is interpreted to intrude the southern part of the MAIC. A new Pb–Pb zircon age of 2701.1 ± 6 Ma is interpreted as the magmatic age of the Fertile batholith.

Minnesota↗
USGS2026

Channel morphology and large wood control postfire debris-flow erosion and deposition

Runoff-generated debris flows are a known response to wildfire, and accurately predicting the volume of these debris flows is important for estimating the magnitude of downstream hazards. Prior data collection efforts have focused on debris-flow volume measurements at catchment outlets, but few studies have considered how erosion and deposition modulate the volume of sediment arriving at catchment outlets. This study takes advantage of a high-resolution dataset to document the factors that control the total debris-flow volume reaching the catchment outlet during a fatal postfire debris flow. Using pre- and post-event airborne lidar, satellite imagery and field mapping, we found that a postfire debris flow in the Black Hollow catchment in northern Colorado eroded 136,000 ± 30,000 m 3 and redeposited 27,000 ± 7,500 m 3 in the main channel. Most of the in-channel deposition (52% by volume) occurred where a confined channel reach transitioned to an unconfined channel reach downstream, allowing the flow to widen and deposit material. Wood jams played multiple roles in the debris-flow dynamics, both nucleating deposition (25% of the deposit volume was stored behind wood jams) and exacerbating erosion (50% of the total erosion occurred downstream from a wood dam break). The remaining deposition occurred due to spatial changes in channel slope as well as deposition observed at newly formed channel bars. Using these data in this study, we identified topographic and vegetation metrics that can be used (pre-event) to anticipate where deposition may occur in channels prior to a debris flow.

Colorado↗
USGS2025

Effect of land cover type on 3D deformation recovery from synthetically deformed high-resolution satellite optical imagery

The limits of detection for earthquake surface deformation in the spatial domain have improved with advances in remote sensing imagery data availability, resolution, and analysis. Sub-pixel correlation and digital elevation model (DEM) differencing from sub-meter, earthquake-spanning satellite optical imagery has enhanced surface rupture mapping and deformation measurements. However, knowledge of measurement accuracy and uncertainty is limited. To address this, we construct orthophotos and digital elevation models (DEMs) from repeat high resolution (∼0.5 m) satellite optical imagery along two sections of the Garlock fault, California with clear fault geomorphology and differing land cover. We deform later sets of DEMs and images with synthetic earthquakes containing both diffuse and discrete horizontal and vertical displacements. Sub-pixel image correlation and DEM differencing demonstrate how vegetation degrades recovered displacement accuracy. In barren land cover, horizontal displacements are detectable to an expected ∼1/10th-pixel size. With shrubs, trees, and grass, detectable displacements increase to >1/2-pixel size, and filtering results by correlation score and using elevation values as input rather than image values improves accuracy. Vertical displacement detection thresholds remain lower in vegetation, at >1-pixel size. Higher slope angles degrade displacement recovery, worsened by vegetation. Diminishing seasonal separation improves accuracy over vegetated regions, though not to the level achieved in barren environments. These results will inform research and operational efforts on the utility of high resolution satellite optical imagery for detecting deformation in varied land cover. Furthermore, they reveal where alternative measurements, such as from LiDAR or radar interferometry, are required to mitigate the effects of vegetation and capture fine-scale crustal deformation.

California↗
USGS2026

Recycling radiogenic osmium by crustal foundering in subduction zones: Evidence from pyroxenite xenoliths in the north Andean arc

Removal of mafic-ultramafic lower crust (e.g., via delamination) is fundamental to making andesitic continental crust, yet direct evidence of this process remains elusive. A unique suite of garnet clinopyroxenite and hornblendite (arclogite) xenoliths from the northern volcanic zone of the Andes, erupted to the surface in the mid-Pleistocene, have bulk-rock geochemistry and osmium (Os) isotopic compositions that demonstrate their crustal affinity, yet have equilibration pressures and temperatures below the arc Mohorovičić discontinuity (~53 km) and some as deep as within the sub-arc mantle wedge (~105 km). Garnet websterites from the same xenolith suite, sourced from depths approaching the Wadati-Benioff zone (~140 km), have elemental and isotopic compositions indicative of a mantle origin, likely formed as products of peridotite-liquid reactions above the subducting slab. Variability in bulk-rock 143 Nd/ 144 Nd values and garnet oxygen isotope ratios ( 𝛿 18 O) for these samples is attributed to minor assimilation of subducted components and/or older crustal material, but assimilation alone cannot account for the highly radiogenic Os isotopic compositions. These results provide direct petrologic evidence for modern arclogite formation and foundering in the Andean orogen, the archetypal active continental subduction system. Elevated 187 Os/ 188 Os values of Mercaderes arclogites imply that recycling of gravitationally unstable arc cumulates during continental crust formation introduces highly radiogenic Os into the convective mantle, which has implications for mass fluxes across the crust-mantle boundary as well as the sources and evolution of mantle heterogeneity as seen in orogenic peridotite massifs and in the sources of oceanic basalts

Earth and Planetary Science Letters↗
USGS2026

Widespread landslide activity in an extreme wet season and implications for regional sediment management, eastern San Francisco Bay area, California

Watershed sediment production is expected to increase in a warmer future with more extreme rain, with cascading effects throughout drainage and sediment-transport networks. This study investigated landscape-scale sediment movement in the eastern San Francisco Bay area, California, USA, during the extreme 2016–2017 wet season that brought major rainfall, landslides, and flooding. Mapping 8,928 landslides across a 1,050-km 2 study area revealed new sediment yield of 510–956 t/km 2 , equivalent to denudation of 193–361 mm/ky. These results correspond closely to long-term denudation rates in the northern and central California Coast Ranges, indicating that mass wasting in very wet years dominates long-term sediment mobilization. However, due to long residence times in drainage networks, the 2017 landslides contributed at most ∼1%–2% of the estimated locally derived fluvial sediment transport to San Francisco Bay. Although the amount of sediment mobilized did not threaten municipal water supplies, small rangeland impoundments in this mixed-use landscape lost storage capacity to new sedimentation. Considering regional sediment supply and demand, even the exceptionally large sediment delivery in an extreme wet year cannot meet the need for sediment to accrete tidal wetlands in the bay. To keep pace with rising sea levels, this abnormally high terrestrial sediment input would need to occur in 50 of the next 75 years, an unlikely occurrence due to the prevalence of recent drought years. Shoreline protection and restoration in the bay would need additional sources of sediment, such as through management of dredged sediment through beneficial-reuse programs.

California↗
USGS2024

Intense alteration on early Mars revealed by high-aluminum rocks at Jezero Crater

The NASA Perseverance rover discovered light-toned float rocks scattered across the surface of Jezero crater that are particularly rich in alumina ( ~ 35 wt% Al 2 O 3 ) and depleted in other major elements (except silica). These unique float rocks have heterogeneous mineralogy ranging from kaolinite/halloysite-bearing in hydrated samples, to spinel-bearing in dehydrated samples also containing a dehydrated Al-rich phase. Here we describe SuperCam and Mastcam-Z observations of the float rocks, including the first in situ identification of kaolinite or halloysite on another planet, and dehydrated phases including spinel and apparent partially dehydroxylated kaolinite. The presence of spinel in these samples is likely detrital in origin, surviving kaolinitization, pointing to an ultramafic origin. However, the association of low hydration with increased Al 2 O 3 abundances suggests heating-induced dehydration which could have occurred during the lithification or impact excavation of these rocks. Given the orbital context of kaolinite-bearing megabreccia in the Jezero crater rim, we propose an origin for these rocks involving intense aqueous alteration of the parent material, followed by dehydration/lithification potentially through impact processes, and dispersion into Jezero crater through flood or impact-related processes.

Communications Earth & Environment↗