USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research”

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 253 records · Page 14Linked to original sources
USGS2016

The role of crystallization-driven exsolution on the sulfur mass balance in volcanic arc magmas

The release of large amounts of sulfur to the stratosphere during explosive eruptions affects the radiative balance in the atmosphere and consequentially impacts climate for up to several years after the event. Quantitative estimations of the processes that control the mass balance of sulfur between melt, crystals, and vapor bubbles is needed to better understand the potential sulfur yield of individual eruption events and the conditions that favor large sulfur outputs to the atmosphere. The processes that control sulfur partitioning in magmas are (1) exsolution of volatiles (dominantly H 2 O) during decompression (first boiling) and during isobaric crystallization (second boiling), (2) the crystallization and breakdown of sulfide or sulfate phases in the magma, and (3) the transport of sulfur-rich vapor (gas influx) from deeper unerupted regions of the magma reservoir. Vapor exsolution and the formation/breakdown of sulfur-rich phases can all be considered as closed-system processes where mass balance arguments are generally easier to constrain, whereas the contribution of sulfur by vapor transport (open system process) is more difficult to quantify. The ubiquitous “excess sulfur” problem, which refers to the much higher sulfur mass released during eruptions than what can be accounted for by amount of sulfur originally dissolved in erupted melt, as estimated from melt inclusion sulfur concentrations (the “petrologic estimate”), reflects the challenges in closing the sulfur mass balance between crystals, melt, and vapor before and during a volcanic eruption. In this work, we try to quantify the relative importance of closed- and open-system processes for silicic arc volcanoes using kinetic models of sulfur partitioning during exsolution. Our calculations show that crystallization-induced exsolution (second boiling) can generate a significant fraction of the excess sulfur observed in crystal-rich arc magmas. This result does not negate the important role of vapor migration in sulfur mass balance but rather points out that second boiling (in situ exsolution) can provide the necessary yield to drive the excess sulfur to the levels observed for crystal-rich systems. In contrast, in crystal-poor systems, magma recharge that releases sulfur-rich bubbles is necessary and most likely the primary contributor to sulfur mass balance. Finally, we apply our model to account for the effect of sulfur partitioning during second boiling and its impact on sulfur released during the Cerro Galan supereruption in Argentina (2.08 Ma) and show the potential importance of second boiling in releasing a large amount of sulfur to the atmosphere during the eruption of large crystal-rich ignimbrites.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Reexamination of the magnitudes for the 1906 and 1922 Chilean earthquakes using Japanese tsunami amplitudes: Implications for source depth constraints

Far-field tsunami records from the Japanese tide gauge network allow the reexamination of the moment magnitudes ( M w ) for the 1906 and 1922 Chilean earthquakes, which to date rely on limited information mainly from seismological observations alone. Tide gauges along the Japanese coast provide extensive records of tsunamis triggered by six great ( M w >8) Chilean earthquakes with instrumentally determined moment magnitudes. These tsunami records are used to explore the dependence of tsunami amplitudes in Japan on the parent earthquake magnitude of Chilean origin. Using the resulting regression parameters together with tide gauge amplitudes measured in Japan we estimate apparent moment magnitudes of M w 8.0–8.2 and M w 8.5–8.6 for the 1906 central and 1922 north-central Chile earthquakes. The large discrepancy of the 1906 magnitude estimated from the tsunami observed in Japan as compared with those previously determined from seismic waves ( M s 8.4) suggests a deeper than average source with reduced tsunami excitation. A deep dislocation along the Chilean megathrust would favor uplift of the coast rather than beneath the sea, giving rise to a smaller tsunami and producing effects consistent with those observed in 1906. The 1922 magnitude inferred from far-field tsunami amplitudes appear to better explain the large extent of damage and the destructive tsunami that were locally observed following the earthquake than the lower seismic magnitudes ( M s 8.3) that were likely affected by the well-known saturation effects. Thus, a repeat of the large 1922 earthquake poses seismic and tsunami hazards in a region identified as a mature seismic gap.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

A new strategy for earthquake focal mechanisms using waveform-correlation-derived relative polarities and cluster analysis: Application to the 2014 Long Valley Caldera earthquake swarm

In microseismicity analyses, reliable focal mechanisms can typically be obtained for only a small subset of located events. We address this limitation here, presenting a framework for determining robust focal mechanisms for entire populations of very small events. To achieve this, we resolve relative P and S wave polarities between pairs of waveforms by using their signed correlation coefficients—a by-product of previously performed precise earthquake relocation. We then use cluster analysis to group events with similar patterns of polarities across the network. Finally, we apply a standard mechanism inversion to the grouped data, using either catalog or correlation-derived P wave polarity data sets. This approach has great potential for enhancing analyses of spatially concentrated microseismicity such as earthquake swarms, mainshock-aftershock sequences, and industrial reservoir stimulation or injection-induced seismic sequences. To demonstrate its utility, we apply this technique to the 2014 Long Valley Caldera earthquake swarm. In our analysis, 85% of the events (7212 out of 8494 located by Shelly et al . [2016]) fall within five well-constrained mechanism clusters, more than 12 times the number with network-determined mechanisms. Of the earthquakes we characterize, 3023 (42%) have magnitudes smaller than 0.0. We find that mechanism variations are strongly associated with corresponding hypocentral structure, yet mechanism heterogeneity also occurs where it cannot be resolved by hypocentral patterns, often confined to small-magnitude events. Small (5–20°) rotations between mechanism orientations and earthquake location trends persist when we apply 3-D velocity models and might reflect a geometry of en echelon, interlinked shear, and dilational faulting.

California↗
USGS2017

Ground-rupturing earthquakes on the northern Big Bend of the San Andreas Fault, California, 800 A.D. to Present

Paleoseismic data on the timing of ground-rupturing earthquakes constrain the recurrence behavior of active faults and can provide insight on the rupture history of a fault if earthquakes dated at neighboring sites overlap in age and are considered correlative. This study presents the evidence and ages for 11 earthquakes that occurred along the Big Bend section of the southern San Andreas Fault at the Frazier Mountain paleoseismic site. The most recent earthquake to rupture the site was the M w 7.7–7.9 Fort Tejon earthquake of 1857. We use over 30 trench excavations to document the structural and sedimentological evolution of a small pull-apart basin that has been repeatedly faulted and folded by ground-rupturing earthquakes. A sedimentation rate of 0.4 cm/yr and abundant organic material for radiocarbon dating contribute to a record that is considered complete since 800 A.D. and includes 10 paleoearthquakes. Earthquakes have ruptured this location on average every ~100 years over the last 1200 years, but individual intervals range from ~22 to 186 years. The coefficient of variation of the length of time between earthquakes (0.7) indicates quasiperiodic behavior, similar to other sites along the southern San Andreas Fault. Comparison with the earthquake chronology at neighboring sites along the fault indicates that only one other 1857-size earthquake could have occurred since 1350 A.D., and since 800 A.D., the Big Bend and Mojave sections have ruptured together at most 50% of the time in M w ≥ 7.3 earthquakes.

California↗
USGS2017

Integrated geophysical characteristics of the 2015 Illapel, Chile, earthquake

On 16 September 2015, a Mw 8.3 earthquake ruptured the subduction zone offshore of Illapel, Chile, generating an aftershock sequence with 14 Mw 6.0–7.0 events. A double source W phase moment tensor inversion consists of a Mw 7.2 subevent and the main Mw 8.2 phase. We determine two slip models for the mainshock, one using teleseismic broadband waveforms and the other using static GPS and InSAR surface displacements, which indicate high slip north of the epicenter and west-northwest of the epicenter near the oceanic trench. These models and slip distributions published in other studies suggest spatial slip uncertainties of ~25 km and have peak slip values that vary by a factor of 2. We relocate aftershock hypocenters using a Bayesian multiple-event relocation algorithm, revealing a cluster of aftershocks under the Chilean coast associated with deep (20–45 km depth) mainshock slip. Less vigorous aftershock activity also occurred near the trench and along strike of the main aftershock region. Most aftershocks are thrust-faulting events, except for normal-faulting events near the trench. Coulomb failure stress change amplitudes and signs are uncertain for aftershocks collocated with deeper mainshock slip; other aftershocks are more clearly associated with loading from the mainshock. These observations reveal a frictionally heterogeneous interface that ruptured in patches at seismogenic depths (associated with many aftershocks) and with homogeneous slip (and few aftershocks) up to the trench. This event likely triggered seismicity separate from the main slip region, including along-strike events on the megathrust and intraplate extensional events.

Illapel↗
USGS2017

Frictional strength of wet and dry montmorillonite

Montmorillonite is a common mineral in fault zones, and its low strength relative to other common gouge minerals is important in many models of fault rheology. However, the coefficient of friction, μ , varies with degree of saturation and is not well constrained in the literature due to the difficulty of establishing fully drained or fully dried states in the laboratory. We measured μ of both saturated and oven-dried montmorillonite at normal stresses up to 700 MPa. Care was taken to shear saturated samples slowly enough to avoid pore fluid overpressure. For saturated samples, μ increased from 0.10 to 0.28 with applied effective normal stress, while for dry samples μ decreased from 0.78 to 0.45. The steady state rate dependence of friction, ( a − b ), was positive, promoting stable sliding. The wide disparity in reported frictional strengths can be attributed to experimental procedures that promote differing degrees of partial saturation or overpressured pore fluid conditions.

Journal of Geophysical Research↗
USGS2017

Sediment gravity flows triggered by remotely generated earthquake waves

Recent great earthquakes and tsunamis around the world have heightened awareness of the inevitability of similar events occurring within the Cascadia Subduction Zone of the Pacific Northwest. We analyzed seafloor temperature, pressure, and seismic signals, and video stills of sediment-enveloped instruments recorded during the 2011–2015 Cascadia Initiative experiment, and seafloor morphology. Our results led us to suggest that thick accretionary prism sediments amplified and extended seismic wave durations from the 11 April 2012 M w 8.6 Indian Ocean earthquake, located more than 13,500 km away. These waves triggered a sequence of small slope failures on the Cascadia margin that led to sediment gravity flows culminating in turbidity currents. Previous studies have related the triggering of sediment-laden gravity flows and turbidite deposition to local earthquakes, but this is the first study in which the originating seismic event is extremely distant (> 10,000 km). The possibility of remotely triggered slope failures that generate sediment-laden gravity flows should be considered in inferences of recurrence intervals of past great Cascadia earthquakes from turbidite sequences. Future similar studies may provide new understanding of submarine slope failures and turbidity currents and the hazards they pose to seafloor infrastructure and tsunami generation in regions both with and without local earthquakes.

Cascadia subduction zone↗
USGS2017

Integrating puffing and explosions in a general scheme for Strombolian-style activity

Strombolian eruptions are among the most common subaerial styles of explosive volcanism worldwide. Distinctive features of each volcano lead to a correspondingly wide range of variations of magnitude and erupted products, but most papers focus on a single type of event at a single volcano. Here, in order to emphasize the common features underlying this diversity of styles, we scrutinize a database from 35 different erupting vents, including 21 thermal infrared videos from Stromboli (Italy), Etna (Italy), Yasur (Vanuatu), and Batu Tara (Indonesia), from puffing, through rapid explosions to normal explosions, with variable ejection parameters and relative abundance of gas, ash, and bombs. Using field observations and high-speed thermal infrared videos processed by a new algorithm, we identify the distinguishing characteristics of each type of activity and how they may relate and interact. In particular, we record that ash-poor normal explosions may be preceded and followed by the onset or the increase of the puffing activity, while ash-rich explosions are emergent, i.e., with inflation of the free surface followed directly by emission of increasingly large gas pockets. Overall, we see that all Strombolian activities form a continuum arising from a common mechanism and are modulated by the combination of two well-established controls: (1) the length of the bursting gas pocket with respect to the vent diameter and (2) the presence and thickness of a high-viscosity layer in the uppermost part of the volcanic conduit.

Journal of Geophysical Research B: 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↗
USGS2016

Processes of multibathyal aragonite undersaturation in the Arctic Ocean

During 3 years of study (2010–2012), the western Arctic Ocean was found to have unique aragonite saturation profiles with up to three distinct aragonite undersaturation zones. This complexity is produced as inflow of Atlantic-derived and Pacific-derived water masses mix with Arctic-derived waters, which are further modified by physiochemical and biological processes. The shallowest aragonite undersaturation zone, from the surface to ∼30 m depth is characterized by relatively low alkalinity and other dissolved ions. Besides local influence of biological processes on aragonite undersaturation of shallow coastal waters, the nature of this zone is consistent with dilution by sea-ice melt and invasion of anthropogenic CO 2 from the atmosphere. A second undersaturated zone at ∼90–220 m depth (salinity ∼31.8–35.4) occurs within the Arctic Halocline and is characterized by elevated p CO 2 and nutrients. The nature of this horizon is consistent with remineralization of organic matter on shallow continental shelves bordering the Canada Basin and the input of the nutrients and CO 2 entrained by currents from the Pacific Inlet. Finally, the deepest aragonite undersaturation zone is at greater than 2000 m depth and is controlled by similar processes as deep aragonite saturation horizons in the Atlantic and Pacific Oceans. The comparatively shallow depth of this deepest aragonite saturation horizon in the Arctic is maintained by relatively low temperatures, and stable chemical composition. Understanding the mechanisms controlling the distribution of these aragonite undersaturation zones, and the time scales over which they operate will be crucial to refine predictive models.

Journal of Geophysical Research: Oceans↗
USGS2017

Behavior of a wave-driven buoyant surface jet on a coral reef

A wave-driven surface buoyant jet exiting a coral reef was studied in order to quantify the amount of water re-entrained over the reef crest. Both moored observations and Lagrangian drifters were used to study the fate of the buoyant jet. To investigate in detail the effects of buoyancy and along-shore flow variations, we developed an idealized numerical model of the system. Consistent with previous work, the ratio of along-shore velocity to jet-velocity and the jet internal Froude number were found to be important determinants of the fate of the jet. In the absence of buoyancy, the entrainment of fluid at the reef crest, creates a significant amount of retention, keeping 60% of water in the reef system. However, when the jet is lighter than the ambient ocean-water, the net effect of buoyancy is to enhance the separation of the jet from shore, leading to a greater export of reef water. Matching observations, our modeling predicts that buoyancy limits retention to 30% of the jet flow for conditions existing on the Moorea reef. Overall, the combination of observations and modeling we present here shows that reef-ocean temperature gradients can play an important role in reef-ocean exchanges.

Journal of Geophysical Research C: Oceans↗
USGS2017

Sediment transport in the presence of large reef bottom roughness

The presence of large bottom roughness, such as that formed by benthic organisms on coral reef flats, has important implications for the size, concentration, and transport of suspended sediment in coastal environments. A 3 week field study was conducted in approximately 1.5 m water depth on the reef flat at Ningaloo Reef, Western Australia, to quantify the cross-reef hydrodynamics and suspended sediment dynamics over the large bottom roughness (∼20–40 cm) at the site. A logarithmic mean current profile consistently developed above the height of the roughness; however, the flow was substantially reduced below the height of the roughness (canopy region). Shear velocities inferred from the logarithmic profile and Reynolds stresses measured at the top of the roughness, which are traditionally used in predictive sediment transport formulations, were similar but much larger than that required to suspend the relatively coarse sediment present at the bed. Importantly, these stresses did not represent the stresses imparted on the sediment measured in suspension and are therefore not relevant to the description of suspended sediment transport in systems with large bottom roughness. Estimates of the bed shear stresses that accounted for the reduced near-bed flow in the presence of large roughness vastly improved the relationship between the predicted and observed grain sizes that were in suspension. Thus, the impact of roughness, not only on the overlying flow but also on bed stresses, must be accounted for to accurately estimate suspended sediment transport in regions with large bottom roughness, a common feature of many shallow coastal ecosystems.

Journal of Geophysical Research C: Oceans↗
USGS2016

Barrier island breach evolution: Alongshore transport and bay-ocean pressure gradient interactions

Physical processes controlling repeated openings and closures of a barrier island breach between a bay and the open ocean are studied using aerial photographs and atmospheric and hydrodynamic observations. The breach site is located on Pea Island along the Outer Banks, separating Pamlico Sound from the Atlantic Ocean. Wind direction was a major control on the pressure gradients between the bay and the ocean to drive flows that initiate or maintain the breach opening. Alongshore sediment flux was found to be a major contributor to breach closure. During the analysis period from 2011 to 2016, three hurricanes had major impacts on the breach. First, Hurricane Irene opened the breach with wind-driven flow from bay to ocean in August 2011. Hurricane Sandy in October 2012 quadrupled the channel width from pressure gradient flows due to water levels that were first higher on the ocean side and then higher on the bay side. The breach closed sometime in Spring 2013, most likely due to an event associated with strong alongshore sediment flux but minimal ocean-bay pressure gradients. Then, in July 2014, Hurricane Arthur briefly opened the breach again from the bay side, in a similar fashion to Irene. In summary, opening and closure of breaches are shown to follow a dynamic and episodic balance between along-channel pressure gradient driven flows and alongshore sediment fluxes.

Pea Island↗
USGS2017

Rip currents and alongshore flows in single channels dredged in the surf zone

To investigate the dynamics of flows near nonuniform bathymetry, single channels (on average 30 m wide and 1.5 m deep) were dredged across the surf zone at five different times, and the subsequent evolution of currents and morphology was observed for a range of wave and tidal conditions. In addition, circulation was simulated with the numerical modeling system COAWST, initialized with the observed incident waves and channel bathymetry, and with an extended set of wave conditions and channel geometries. The simulated flows are consistent with alongshore flows and rip-current circulation patterns observed in the surf zone. Near the offshore-directed flows that develop in the channel, the dominant terms in modeled momentum balances are wave-breaking accelerations, pressure gradients, advection, and the vortex force. The balances vary spatially, and are sensitive to wave conditions and the channel geometry. The observed and modeled maximum offshore-directed flow speeds are correlated with a parameter based on the alongshore gradient in breaking-wave-driven-setup across the nonuniform bathymetry (a function of wave height and angle, water depths in the channel and on the sandbar, and a breaking threshold) and the breaking-wave-driven alongshore flow speed. The offshore-directed flow speed increases with dissipation on the bar and reaches a maximum (when the surf zone is saturated) set by the vertical scale of the bathymetric variability.

Journal of Geophysical Research: Oceans↗
USGS2017

Observations and a linear model of water level in an interconnected inlet-bay system

A system of barrier islands and back-barrier bays occurs along southern Long Island, New York, and in many coastal areas worldwide. Characterizing the bay physical response to water level fluctuations is needed to understand flooding during extreme events and evaluate their relation to geomorphological changes. Offshore sea level is one of the main drivers of water level fluctuations in semienclosed back-barrier bays. We analyzed observed water levels (October 2007 to November 2015) and developed analytical models to better understand bay water level along southern Long Island. An increase (∼0.02 m change in 0.17 m amplitude) in the dominant M 2 tidal amplitude (containing the largest fraction of the variability) was observed in Great South Bay during mid-2014. The observed changes in both tidal amplitude and bay water level transfer from offshore were related to the dredging of nearby inlets and possibly the changing size of a breach across Fire Island caused by Hurricane Sandy (after December 2012). The bay response was independent of the magnitude of the fluctuations (e.g., storms) at a specific frequency. An analytical model that incorporates bay and inlet dimensions reproduced the observed transfer function in Great South Bay and surrounding areas. The model predicts the transfer function in Moriches and Shinnecock bays where long-term observations were not available. The model is a simplified tool to investigate changes in bay water level and enables the evaluation of future conditions and alternative geomorphological settings.

Journal of Geophysical Research C: Oceans↗
USGS2017

Physical response of a back-barrier estuary to a post-tropical cyclone

This paper presents a modeling investigation of the hydrodynamic and sediment transport response of Chincoteague Bay (VA/MD, USA) to Hurricane Sandy using the Coupled Ocean-Atmosphere-Wave-Sediment-Transport (COAWST) modeling system. Several simulation scenarios with different combinations of remote and local forces were conducted to identify the dominant physical processes. While 80% of the water level increase in the bay was due to coastal sea level at the peak of the storm, a rich spatial and temporal variability in water surface slope was induced by local winds and waves. Local wind increased vertical mixing, horizontal exchanges, and flushing through the inlets. Remote waves (swell) enhanced southward flow through wave setup gradients between the inlets, and increased locally generated wave heights. Locally generated waves had a negligible effect on water level but reduced the residual flow up to 70% due to enhanced apparent roughness and breaking-induced forces. Locally generated waves dominated bed shear stress and sediment resuspension in the bay. Sediment transport patterns mirrored the interior coastline shape and generated deposition on inundated areas. The bay served as a source of fine sediment to the inner shelf, and the ocean-facing barrier island accumulated sand from landward-directed overwash. Despite the intensity of the storm forcing, the bathymetric changes in the bay were on the order of centimeters. This work demonstrates the spectrum of responses to storm forcing, and highlights the importance of local and remote processes on back-barrier estuarine function.

Chincoteague Bay↗
USGS2017

Compositional variations in sands of the Bagnold Dunes, Gale Crater, Mars, from visible-shortwave infrared spectroscopy and comparison with ground truth from the Curiosity Rover

During its ascent up Mount Sharp, the Mars Science Laboratory Curiosity rover traversed the Bagnold Dune Field. We model sand modal mineralogy and grain size at four locations near the rover traverse, using orbital shortwave infrared single scattering albedo spectra and a Markov-Chain Monte Carlo implementation of Hapke's radiative transfer theory to fully constrain uncertainties and permitted solutions. These predictions, evaluated against in situ measurements at one site from the Curiosity rover, show that XRD-measured mineralogy of the basaltic sands is within the 95% confidence interval of model predictions. However, predictions are relatively insensitive to grain size and are non-unique, especially when modeling the composition of minerals with solid solutions. We find an overall basaltic mineralogy and show subtle spatial variations in composition in and around the Bagnold dunes, consistent with a mafic enrichment of sands with cumulative transport distance by sorting of olivine, pyroxene, and plagioclase grains during aeolian saltation. Furthermore, the large variations in Fe and Mg abundances (~20 wt%) at the Bagnold Dunes suggest that compositional variability induced by wind sorting may be enhanced by local mixing with proximal sand sources. Our estimates demonstrate a method for orbital quantification of composition with rigorous uncertainty determination and provide key constraints for interpreting in situ measurements of compositional variability within martian aeolian sandstones.

Journal of Geophysical Research E: Planets↗
USGS2016

Geologic history of Martian regolith breccia Northwest Africa 7034: Evidence for hydrothermal activity and lithologic diversity in the Martian crust

The timing and mode of deposition for Martian regolith breccia Northwest Africa (NWA) 7034 were determined by combining petrography, shape analysis, and thermochronology. NWA 7034 is composed of igneous, impact, and brecciated clasts within a thermally annealed submicron matrix of pulverized crustal rocks and devitrified impact/volcanic glass. The brecciated clasts are likely lithified portions of Martian regolith with some evidence of past hydrothermal activity. Represented lithologies are primarily ancient crustal materials with crystallization ages as old as 4.4 Ga. One ancient zircon was hosted by an alkali-rich basalt clast, confirming that alkalic volcanism occurred on Mars very early. NWA 7034 is composed of fragmented particles that do not exhibit evidence of having undergone bed load transport by wind or water. The clast size distribution is similar to terrestrial pyroclastic deposits. We infer that the clasts were deposited by atmospheric rainout subsequent to a pyroclastic eruption(s) and/or impact event(s), although the ancient ages of igneous components favor mobilization by impact(s). Despite ancient components, the breccia has undergone a single pervasive thermal event at 500–800°C, evident by groundmass texture and concordance of ~1.5 Ga dates for bulk rock K-Ar, U-Pb in apatite, and U-Pb in metamict zircons. The 1.5 Ga age is likely a thermal event that coincides with rainout/breccia lithification. We infer that the episodic process of regolith lithification dominated sedimentary processes during the Amazonian Epoch. The absence of pre-Amazonian high-temperature metamorphic events recorded in ancient zircons indicates source domains of static southern highland crust punctuated by episodic impact modification.

Journal of Geophysical Research E: Planets↗