Katherine R. Barnhart

An Exhumation History of Continents over Billion-Year Time Scales

Continental Thermocouple The patchy presence of billions-of-years-old continental crust indicates a complex coupling between the buoyant forces keeping the lithosphere floating on the mantle and the persistent erosional forces gradually wearing the crust away. Measuring long-term rates of exhumation—the creation of new rock surfaces due to erosion—can reveal how the crust is thermally coupled to the underlying mantle, but techniques to do so have often only been able to resolve a limited temperature range across narrow slices of geologic time. Blackburn et al. (p. 73) used uranium-lead thermochronology, which is sensitive to the much higher temperatures representative of lower crustal depths, to construct a long-term quantitative model of exhumation and erosion for North America. Thermochronology indicates a balance between low erosion rates and slow thermal cooling in old continental crust. The continental lithosphere contains the oldest and most stable structures on Earth, where fragments of ancient material have eluded destruction by tectonic and surface processes operating over billions of years. Although present-day erosion of these remnants is slow, a record of how they have uplifted, eroded, and cooled over Earth’s history can provide insight into the physical properties of the continents and the forces operating to exhume them over geologic time. We constructed a continuous record of ancient lithosphere cooling with the use of uranium-lead (U-Pb) thermochronology on volcanically exhumed lower crustal fragments. Combining these measurements with thermal and Pb-diffusion models constrains the range of possible erosion histories. Measured U-Pb data are consistent with extremely low erosion rates persisting over time scales approaching the age of the continents themselves.

Modeling erosion of ice‐rich permafrost bluffs along the Alaskan Beaufort Sea coast

The Arctic climate is changing, inducing accelerating retreat of ice-rich permafrost coastal bluffs. Along Alaskas Beaufort Sea coast, erosion rates have increased roughly threefold from 6.8 to 19 m yr−1 since 1955 while the sea ice-free season has increased roughly twofold from 45 to 100 days since 1979. We develop a numerical model of bluff retreat to assess the relative roles of the length of sea ice-free season, sea level, water temperature, nearshore wavefield, and permafrost temperature in controlling erosion rates in this setting. The model captures the processes of erosion observed in short-term monitoring experiments along the Beaufort Sea coast, including evolution of melt notches, topple of ice wedge-bounded blocks, and degradation of these blocks. Model results agree with time-lapse imagery of bluff evolution and time series of ocean-based instrumentation. Erosion is highly episodic with 40% of erosion is accomplished during less than 5% of the sea ice-free season. Among the formulations of the submarine erosion rate we assessed, we advocate those that employ both water temperature and nearshore wavefield. As high water levels are a prerequisite for erosion, any future changes that increase the frequency with which water levels exceed the base of the bluffs will increase rates of coastal erosion. The certain increases in sea level and potential changes in storminess will both contribute to this effect. As water temperature also influences erosion rates, any further expansion of the sea ice-free season into the midsummer period of greatest insolation is likely to result in an additional increase in coastal retreat rates.

Deep crustal xenoliths from central Montana, USA: Implications for the timing and mechanisms of high-velocity lower crust formation

Integration of petrologic, chronologic and petrophysical xenolith data with geophysical observations can offer fundamental insights into understanding the evolution of continental crust. We present the results of a deep crustal xenolith study from the northern Rocky Mountain region of the western U.S., where seismic experiments reveal an anomalously thick (10–30 km), high seismic velocity (compressional body wave, Vp > 7.0 km/s) lower crustal layer, herein referred to as the 7.x layer. Xenoliths exhumed by Eocene minettes from the Bearpaw Mountains of central Montana, within the Great Falls tectonic zone, include mafic and intermediate garnet granulites, mafic hornblende eclogite, and felsic granulites. Calculated pressures of 0.6–1.5 GPa are consistent with derivation from 23–54 km depths. Samples record diverse and commonly polymetamorphic pressure-temperature histories including prograde burial and episodes of decompression. Samples with barometrically determined depths consistent with residence within the seismically defined 7.x layer have calculated bulk P-wave velocities of 6.9–7.8 km/s, indicating heterogeneity in the layer. Shallower samples have markedly slower velocities consistent with seismic models. New monazite total U-Th-Pb data and a variety of additional published geochronology indicate a prolonged and episodic metamorphic history, beginning with protolith ages as old as Archean and followed by metamorphic and deep crustal fluid-flow events ca. 2.1 Ga, 1.8–1.7 Ga, and 1.5–1.3 Ga. We suggest that the 7.x layer in this region owes its character to a variety of processes, including magmatic underplating and intraplating, associated with multiple tectonic events from the Neoarchean to the Mesoproterozoic.

Quantifying the stratigraphic completeness of delta shoreline trajectories

Understanding the incomplete nature of the stratigraphic record is fundamental for interpreting stratigraphic sequences. Methods for quantifying stratigraphic completeness for one-dimensional stratigraphic columns, defined as the proportion of time intervals of some length that contain stratigraphy, are commonplace; however, quantitative assessments of completeness in higher dimensions are lacking. Here we present a metric for defining stratigraphic completeness of two-dimensional shoreline trajectories using topset-foreset rollover positions in dip-parallel sections and describe the preservation potential of a shoreline trajectory derived from the geometry of the delta surface profile and the kinematics of the geomorphic shoreline trajectory. Two end-member forward models are required to fully constrain the preservation potential of the shoreline dependent on whether or not a topset is eroded during base level fall. A laboratory fan-delta was constructed under nonsteady boundary conditions, and one-dimensional stratigraphic column and two-dimensional shoreline completeness curves were calculated. Results are consistent with the hypothesis derived from conservation of sediment mass that completeness over all timescales should increase given increasing dimensions of analysis. Stratigraphic trajectories and completeness curves determined from forward models using experimental geomorphic trajectories compare well to values from transects when subsampled to the equivalent stratigraphic resolution as observed in the actual preserved sequence. The concept of stratigraphic completeness applied to two-dimensional trajectory analysis and the end-member forward models presented here provide novel tools for a conceptual understanding of the nature of stratigraphic preservation at basin scales.

Decompression during Late Proterozoic Al2SiO5 Triple-Point Metamorphism at Cerro Colorado, New Mexico

An outstanding problem in understanding the late Proterozoic tectonic assembly of the southwest is identifying the tectonic setting associated with regional metamorphism at 1.4 Ga. Both isobaric heating and cooling, and counterclockwise looping PT paths are proposed for this time. We present a study of the Proterozoic metamorphic and deformation history of the Cerro Colorado area, southern Tusas Mountains, New Mexico, which shows that the metamorphism in this area records near-isothermal decompression from 6 to 4 kbar at ca. 1.4 Ga. We do not see evidence for isobaric heating at this time. Decompression from peak pressures is recorded by the reaction , with a negative slope in PT space; the reaction , which is nearly horizontal in PT space; and partial to total pseudomorphing of kyanite by sillimanite during the main phase of deformation. The clearest reaction texture indicating decompression near peak metamorphic temperature is the replacement of garnet by clots of sillimanite, which are surrounded by halos of biotite. The sillimanite clots, most without relict garnet in the cores and with highly variable aspect ratios, are aligned. They define a lineation that formed with the dominant foliation. An inverted metamorphic gradient is locally defined by sillimanite-garnet schists (625°C) structurally above staurolite-garnet schists (550°C) and implies ductile thrusting during the main phase of deformation. The exhumation that led to the recorded decompression was likely in response to crustal thickening due to ductile thrusting and subsequent denudation.

Lithology: A Landlab submodule for spatially variable rock properties

1 University of Colorado at Boulder, Department of Geological Sciences 2 University of Colorado at Boulder, Cooperative Institute for Research in Environmental Sciences 3 University of Colorado at Boulder, Community Surface Dynamics Modeling System Integration Facility 4 University of Colorado at Boulder, Institute for Arctic and Alpine Research 5 Tulane University, Department of Earth and Environmental Sciences