Advances and emerging challenges in luminescence geochronology: Insights from the 2024 New World Luminescence Dating Workshop
No abstract available.
Quaternary Research↗
Geology topics
Harrison J. Gray
Publications and source records attributed to Harrison J. Gray.
Late Miocene Colorado River arrival in the Bidahochi basin supports spillover origin of Grand Canyon
The timing and mechanism of the integration of the Colorado River and incision of the Grand Canyon remain among geology’s enduring controversies. A key question is the configuration of the upper Colorado River watershed between 11 and 6 million years ago. In this study, we present new evidence from zircon uranium-lead geochronology for the arrival of distinctive Colorado–Green River sediment in the Bidahochi basin by 6.6 million years ago derived from the Browns Park Formation. This is coeval with an order-of-magnitude increase in depositional rate, an increase in carbonate strontium isotope ( 87 Sr/ 86 Sr) ratios, the appearance of large fish species characteristic of fast-flowing waters, and other sedimentological changes. This evidence is consistent with the Colorado River supplying water and sediment to the Bidahochi basin before spillover integration of the river through the Grand Canyon.
Arizona, Nevada, Utah↗
Spatiotemporal variations in strain release and seismic rupture in multifault systems: An example from Panamint Valley, southeastern California
Geometrically complex, multifault ruptures have been observed in recent, damaging earthquakes in southeastern California, sparking renewed efforts to identify physical conditions that promote or inhibit fault discontinuity-spanning coseismic ruptures. The likelihood of ruptures propagating across fault discontinuities is thought to be partly controlled by fault geometries, rupture direction, and the history of strain release. However, these parameters vary in space and time over multiple earthquake cycles, making it difficult to forecast the likelihood that an earthquake on one fault will trigger rupture on a nearby fault. Here we use tectono-geomorphic mapping of a geometrically complex fault zone in Panamint Valley, southeastern California, to assess spatiotemporal variations of paleo-rupture patterns and geometries of fault discontinuities over multiple earthquake cycles. First, we identify ten generations of late Pleistocene to Holocene alluvium using geomorphic parameters and luminescence dating to constrain ages of alluvium and bracket late Holocene earthquake timing. Then, we quantify slip kinematics using high-resolution structure from motion digital surface models. We find the Panamint Valley transtensional relay (PVTR) hosted four late Holocene earthquakes, bracketed to ~5.8–3.4 ka, ~3.8–2.2 ka, ~2.4–0.6 ka, and ~0.64–0.16 ka, with ~0.6–1.1 m of slip per event, correlative to Mw ≈ 6.7–6.9 earthquakes. Additionally, we find similarities in earthquake timing on the Ash Hill, PVTR, and Panamint Valley faults and similarities in the slip magnitude and slip kinematics between the Ash Hill and PVTR faults, implying that the PVTR may co-rupture with nearby faults. Paleo-rupture patterns indicate that seismogenic strain transfer may occur through the PVTR, along different combinations of fault segments and jump distances, over multiple earthquake cycles. These data highlight the utility of tectono-geomorphic mapping in evaluating paleo-rupture patterns and suggest that the PVTR may act to propagate and/or arrest rupture between the Ash Hill and Panamint Valley faults.
California↗
Evidence for nonlocal sediment transport on hillslopes from fault scarp morphology
Hillslope sediment transport processes such as bioturbation, rainsplash, and granular mechanics occur across the entire planet. Yet, it remains uncertain how these small-scale processes act together to shape landscapes. Longstanding hillslope diffusion theory posits that hillslope processes are spatially limited, whereas new concepts of nonlocal sediment transport argue otherwise. However, each theory produces subtly different, but distinct, predictions for the evolution of fault scarps. We use the topographic change of fault scarps to demonstrate that hillslope processes produce nonlocal sediment transport. Analysis of a global compilation of 340 dated single-earthquake scarp profiles reveals a statistically significant ( p < 0.05 ) relationship between scarp age and scarp asymmetry, here defined as the ratio of imaginary to real components of the Fourier transform of absolute slope. Numerical simulations show that nonlocal models predict this relationship, whereas hillslope diffusion models do not. To further investigate this result, we examined the depositional geometry of a well-exposed colluvial wedge along the Wasatch fault in central Utah, United States. Our quantitative comparison between the exposure and numerical simulations reveals better agreement with the nonlocal model. Nonlocal sediment transport theory appears to best capture the physics of how hillslope processes shape fault scarps, yet hillslope diffusion provides a useful approximation in many cases. As the processes that act on fault scarps are nearly identical to those acting on hillslopes, our results provide evidence supporting nonlocality as a generalized model of hillslope sediment transport.
Utah↗
Timing and geometry of the Chemehuevi Formation reveal a late Pleistocene sediment pulse into the Lower Colorado River
The Chemehuevi Formation is a distinctive 50−150-m-thick wedge-shaped Pleistocene sedimentary unit deposited by the Colorado River. It lines the perimeters of the river’s floodplains and bedrock canyons for more than 600 km between the mouth of the Grand Canyon and the delta region in the Gulf of California. The formation is composed of a basal tan to light-yellowish-brown and pale-orange mud-dominated facies overlain and interbedded by a light-yellow-brown sand-dominated facies. The unit is one of two extensively exposed aggradational packages in the Lower Colorado River corridor, in addition to a series of other smaller alluvial terrace deposits. The Chemehuevi Formation appears to represent the response of a fully integrated Colorado River system to a significant perturbation, in contrast to the Bullhead Alluvium, which is likely a unique result of Pliocene river integration. The aggradation of the Chemehuevi Formation in the Lower Colorado River corridor may be similarly due to a unique event in the Colorado River system, or it may instead be a well-preserved sedimentary sequence recording typical behavior of the Colorado River below the Grand Canyon in the late Pleistocene. As such, multiple causal mechanisms have been proposed, but no study to date has conclusively explained the Chemehuevi Formation. To help resolve its timing, duration, and origin, we applied post-infrared infrared stimulated luminescence, carbonate U-Th series, and zircon sensitive high-resolution ion microprobe U-Th series geochronology to determine the ages of key exposures of the unit over a wide spatial area. These new data demonstrate that the Chemehuevi Formation was deposited ca. 110−90 ka. The depositional ages collectively overlap, suggesting that deposition occurred rapidly relative to the resolution of the geochronometers. The new depositional timing coincides with a shift from glacial to interglacial conditions after the marine isotope stage 5-6 transition. This observation is consistent with a climate-induced sediment pulse as a causal mechanism, yet correlations with similar deposits in the Colorado River headwaters or in neighboring catchments appear elusive. Potentially, climate transitions between glacial and interglacial periods induced a sediment pulse from hillslopes of the Colorado River system that resulted in the Chemehuevi Formation. An alternative or additional explanation is that the Chemehuevi Formation represents release of lava dam−impounded sediment in the Grand Canyon. The surface geometry of the Chemehuevi Formation projects upstream to the approximate location of lava dams, and the largest possible lava dam impoundment (the Upper Prospect dam) is comparable in volume to the formation. The lava dam hypothesis appears to be a possible explanation for the Chemehuevi Formation. However, tying deposition to a specific lava dam or series of lava dams remains challenging due to discrepancies in timing and volume. The combined effects of a series of lava dams may have led to the Chemehuevi Formation, as the last Pleistocene lava dam eruption coincides with the onset of deposition. Alternatively, the formation may result from the combined effects of both regional climate transitions and the lava dams that created a transient reservoir to compound a climate transition−driven sediment pulse. The geochronologic data presented here do not allow us to distinguish between the lava dam or climate transition hypotheses but will need to be reconciled with any future proposed depositional model.
Arizona, California, Nevada↗
Neotectonic mapping of Puerto Rico
Puerto Rico is part of the Puerto Rico-Virgin Islands microplate, along the Caribbean–North American plate boundary between the Puerto Rico trench subduction zone and the Muertos Trough incipient subduction zone. Despite recent seismicity and geodetically constrained deformation of ~3 mm/yr of left-lateral shear across the island, Quaternary fault locations remain largely uncertain. Preservation of recent faulting in the landscape is masked by distributed faulting, high weathering rates of the tropical climate, steep topography, frequent landsliding, and extensive agriculture and urbanization along coastlines and river valleys. We present remote neotectonic mapping of historical aerial imagery from the 1930s-60s and <1-m lidar topography, integrated with field observations, to create an updated active fault map of Puerto Rico. We focus on faults that offset younger geomorphic surfaces, with ages estimated by geologic mapping, OSL and radiocarbon dating, and morphologic interpretations. We present new evidence for Quaternary activity on seven faults, including the South Lajas, Salinas, Punta Montalva, Great Southern Puerto Rico, Cerro Goden, Parguera, and San Marcos faults. We find that active faulting occurs preferentially along the southern and western coasts and does not always spatially coincide with preexisting bedrock faults, possibly reflecting partial reactivation of older faults in the modern strain field.
Puerto Rico↗
The geochronology of White Sands Locality 2 is resolved
Rhode et al. (2024) allege that there are many “unresolved issues” with the geochronology of White Sands National Park (WHSA) Locality 2. They suggest there are substantial age offsets due to hard-water effects in the aquatic plants that were dated and that radiocarbon ages of pollen may be anomalously old due to reworking. In their view, the luminescence ages are likely to be maximum ages because of the probable presence of partially bleached quartz grains, overestimation of water content, and stratigraphic position of the samples. They also assert the ages of the footprint trackways are not as internally consistent as suggested and can be interpreted in various ways. We review each of these issues and show they are without merit, often irrelevant, at odds with first principles, and stem from a lack of firsthand understanding of the studies we conducted at WHSA Locality 2.
New Mexico↗
Paleoseismology of the Sawtooth fault and implications for fault behavior in the epicentral region of the 2020 Mw 6.5 Stanley, Idaho, earthquake
The 2020 moment magnitude ( M w ) 6.5 Stanley, Idaho, earthquake raised questions about the history and extent of complex faulting in the northwestern Centennial Tectonic Belt (CTB) and its relation to the Sawtooth normal fault and Eocene Trans‐Challis fault system (TCFS). To explore faulting in this area, we excavated a paleoseismic trench across the Sawtooth fault along the western margin of the CTB, and compared an early Holocene (9.1 ± 2.1 ka, 1δ ) rupture at the site with lacustrine paleoseismic data and fault mapping in the 2020 epicentral region. We find: (1) a history of partial to full rupture of the Sawtooth fault ( M w 6.8–7.4), (2) that shorter ruptures ( M w ≤ 6.9 ) are likely along distributed and discontinuous faults in the epicentral region, (3) that this complex system that hosted the 2020 earthquake is not directly linked to the Sawtooth fault, (4) that the northeast‐trending TCFS likely plays a role in controlling fault length and rupture continuity for adjacent faults, and (5) that parts of the TCFS may facilitate displacement transfer between normal faults that accommodate crustal extension and rotation. Our results help unravel complex faulting in the CTB and imply that relict structures can help inform regional seismic hazard assessments.
Idaho↗
High-resolution geophysical and geochronological analysis of a relict shoreface deposit offshore central California: Implications for slip rate along the Hosgri fault
The Cross-Hosgri slope is a bathymetric lineament that crosses the main strand of the Hosgri fault offshore Point Estero, central California. Recently collected chirp seismic reflection profiles and sediment cores provide the basis for a reassessment of Cross-Hosgri slope origin and the lateral slip rate of the Hosgri fault based on offset of the lower slope break of the Cross-Hosgri slope. The Cross-Hosgri slope is comprised of two distinct stratigraphic units. The lower unit (unit 1) overlies the post–Last Glacial Maximum transgressive erosion surface and is interpreted as a Younger Dryas (ca. 12.85–11.65 ka) shoreface deposit based on radiocarbon and optically stimulated luminescence (OSL) ages, Bayesian age modeling, seismic facies, sediment texture, sediment infauna, and heavy mineral component. The shoreface was abandoned and partly eroded during rapid sea-level rise from ca. 11.5 to 7 ka. Unit 2 consists of fine sand and silt deposited in a midshelf environment when the rate of sea-level rise slowed between ca. 7 ka and the present. Although unit 2 provides a thin, relatively uniform cover over the lower slope break of the older shoreface, this feature still represents a valuable piercing point, providing a Hosgri fault slip rate of 2.6 ± 0.8 mm/yr. Full-waveform processing of chirp data resulted in significantly higher resolution in coarser-grained strata, which are typically difficult to interpret with more traditional envelope processing. Our novel combination of offshore radiocarbon and OSL dating is the first application to offshore paleoseismic studies, and our results indicate the utility of this approach for future marine neotectonic investigations.
California↗
Paleohydrologic history of pluvial lake San Agustin, New Mexico: Tracking changing effective moisture in southwest North America through the last glacial transition
Paleoclimate records across the Intermountain West region of North America show significant regional variation in timing and magnitude of wet conditions that accompanied the last glacial-interglacial transition. To understand the climate controls on paleohydrologic change, well-dated records are needed across the region. The Plains of San Agustin (New Mexico, USA) is a closed-basin watershed of the American Southwest influenced by both winter westerly and summer North American Monsoon precipitation. The flat valley floors of the Plains contain lake and groundwater discharge deposits that record multiple periods of past wet climate. We present a record of hydroclimate for the past 26,000 years based on radiocarbon, U–Th series, and OSL dating of these deposits and stratigraphic correlation across the three sub-basins of the lake system. We find that two major lake oscillations occurred, coincident with the global Last Glacial Maximum (∼23–18 ka) and with Heinrich Stadial 1 (∼17–14 ka). The LGM lake cycle created a deep lake in the lowermost sub-basin, fed by marsh/lake overflow in the upper sub-basins. The Heinrich Stadial 1 wet interval attained the highest recorded lake level between ∼17.0 and 15.3 ka, merging the lower two sub-basins into one lake. Both lake cycles agree well in timing and magnitude with other lake-based records from the southwestern U.S., supporting previous interpretations that a south-shifted cool season storm track brought a dipole-like pattern of enhanced moisture to the southwest at the expense of drier conditions in the northwest during the LGM and deglaciation. A transition from lake to groundwater discharge conditions followed during the Bølling-Allerød (14.7–12.9 ka) across the study area. Wet meadows prevailed in the lower sub-basin during the Younger Dryas (12.9–11.7 ka), with marsh and open lake conditions in the upper and middle sub-basins, respectively. During the early Holocene, discrete wet intervals are recorded by ages for wet meadow deposits in all sub-basins with centennial-millennial frequency at 9.9, 8.8, 8.2 ka. These events agree well with other Intermountain West records showing wetter-than-present conditions prior to 8 ka. Two additional wet periods, centered at 6.8 and 5.4 ka, occurred during the driest known interval of the middle Holocene, and likely were partly supported in the Plains of San Agustin by summer moisture associated with the peak strength of the North American Monsoon. Our record highlights that both winter and summer moisture support water resources in New Mexico watersheds.
New Mexico↗
Independent age estimates resolve the controversy of ancient human footprints at White Sands
Human footprints at White Sands National Park, New Mexico, USA, reportedly date to between ~23,000 and 21,000 years ago according to radiocarbon dating of seeds from the aquatic plant Ruppia cirrhosa . These ages remain controversial because of potential old carbon reservoir effects that could compromise their accuracy. We present new calibrated 14C ages of terrestrial pollen collected from the same stratigraphic horizons as those of the Ruppia seeds, along with optically stimulated luminescence ages of sediments from within the human footprint–bearing sequence, to evaluate the veracity of the seed ages. The results show that the chronologic framework originally established for the White Sands footprints is robust and reaffirm that humans were present in North America during the Last Glacial Maximum.
New Mexico↗
Construction and modification of debris-flow alluvial fans as captured in the geomorphic and sedimentary record: Examples from the western Sangre de Cristo Mountains, south-central Colorado
Debris-flow alluvial fans are iconic features of dynamic landscapes and are hypothesized to record tectonic and climatic change. Here, we highlight their complex formation and evolution through an exemplary suite of Quaternary debris-flow alluvial fans emanating from the western range front of the Sangre de Cristo Mountains in south-central Colorado, USA. To evaluate the constructive and modifying processes that produce fan form and the associated sedimentary signatures, we applied a combined geomorphologic and sedimentologic approach using sedimentary facies analysis, soils mapping, high-resolution topographic data, and luminescence geochronology to document timing of fan construction and modification. We explored two subsets of fans in the study area: a southern set sourced from the extensively glaciated drainages of the Blanca Peak massif, and a northern set from the unglaciated drainages south of Great Sand Dunes National Park. Both sets of fans have: (1) active and successively abandoned surfaces that show evolving degradation of primary features through modification by secondary processes, (2) associated facies that display distinct characteristics representative of primary depositional and secondary modifying sedimentary processes, and (3) evidence of primary debris flow with subsequent modification by secondary processes. We found that surface geomorphology and facies assemblages in exposed alluvial-fan deposits represent sediment transport processes on both active and abandoned lobes. The link between fan surface morphologies and the sedimentary facies of their deposits provides a basis for an evolutionary process–based interpretation of debris-flow alluvial-fan geomorphology and provides a better understanding of complexities in buried paleosurfaces (intraformational progressive unconformities), surficial deformation, and landform development as recorded in debris-flow fan deposits in the sedimentary record.
Colorado↗
Mid- to late Quaternary geomorphic and paleoseismic event history, Cheraw fault, southeastern Colorado
Despite its subdued expression and isolated location within the Great Plains of southeastern Colorado, the 80 km long Cheraw fault may be one of the most active faults in North America east of the Rocky Mountains. We present geomorphic analyses, geochronology, and paleoseismic trenching data to (1) document the rupture history of the ∼45 km long southwestern section of the Cheraw fault over the past ∼19 ka, and (2) evaluate slip rate changes for the entire fault over the past ∼200 ka. Results from new trenches excavated at the Old Ranch site show evidence of four surface rupture events since ∼19 ka, each with an average vertical displacement of 0.75 ± 0.15 m. An additional event is likely only slightly older than ∼19 ka. Evidence for relatively small displacements at and near the Old Ranch site indicates that most of these earthquakes were magnitude ( M ) 7 or less and likely did not rupture the full length of the Cheraw fault. Since ∼19 ka, the average slip rate is ∼0.16 ± 0.3 mm/yr near the Old Ranch site with an average interevent time of 3–5 ka. New geochronologic data for mid‐ to late‐Quaternary geomorphic surfaces cut by the Cheraw fault imply rapid incision by local Arkansas River tributaries from ∼145 to ∼100 ka. Maximum vertical offsets of 7–9 m for these surfaces indicate that from ∼19 to >200 ka the average slip rate was no greater than ∼0.03 mm/yr. The accelerated slip rate since ∼19 ka indicates a possible response to rapid erosional unloading and/or a limited late Cenozoic, <40 ka, paleoseismic history for the Cheraw fault.
Colorado↗
A geomorphic-process-based cellular automata model of colluvial wedge morphology and stratigraphy
The development of colluvial wedges at the base of fault scarps following normal-faulting earthquakes serves as a sedimentary record of paleoearthquakes and is thus crucial in assessing seismic hazard. Although there is a large body of observations of colluvial wedge development, connecting this knowledge to the physics of sediment transport can open new frontiers in our understanding. To explore theoretical colluvial wedge evolution, we develop a cellular automata model driven by the production and disturbance (e.g., bioturbative reworking) of mobile regolith and fault-scarp collapse. We consider both 90 and 60 ∘ dipping faults and allow the colluvial wedges to develop over 2000 model years. By tracking sediment transport time, velocity, and provenance, we classify cells into analogs for the debris and wash sedimentary facies commonly described in paleoseismic studies. High values of mobile regolith production and disturbance rates produce relatively larger and more wash-facies-dominated wedges, whereas lower values produced relatively smaller, debris-facies-dominated wedges. Higher lateral collapse rates lead to more debris facies relative to wash facies. Many of the modeled colluvial wedges fully developed within 2000 model years after the earthquake, with many being much faster when process rates are high. Finally, for scenarios with the same amount of vertical displacement, differently sized colluvial wedges developed depending on the rates of geomorphic processes and fault dip. A change in these variables, say by environmental change such as precipitation rates, could theoretically result in different colluvial wedge facies assemblages for the same characteristic earthquake rupture scenario. Finally, the stochastic nature of collapse events, when coupled with high disturbance, illustrates that multiple phases of colluvial deposition are theoretically possible for a single earthquake event.
Earth Surface Dynamics↗
How similar was the 1983 Mw 6.9 Borah Peak earthquake rupture to its surface-faulting predecessors along the northern Lost River fault zone (Idaho, USA)?
We excavated trenches at two paleoseismic sites bounding a trans-basin bedrock ridge (the Willow Creek Hills) along the northern Lost River fault zone to explore the uniqueness of the 1983 M w 6.9 Borah Peak earthquake compared to its prehistoric predecessors. At the Sheep Creek site on the southernmost Warm Springs section, two earthquakes occurred at 9.8−14.0 ka (95% confidence) and 6.5−7.1 ka; each had ∼1.9 m of vertical displacement. About 4 km to the southeast, across the Willow Creek Hills, two ruptures at the Arentson Gulch site on the northernmost Thousand Springs section occurred at 9.0−14.7 ka and 6.1−7.5 ka with ∼1.9 m of vertical displacement each. We synthesize these and previous paleoseismic results into a model of five postglacial (<15 ka) ruptures along a ∼65 km reach of the northern Lost River fault zone. Our results show that the Borah Peak earthquake (34 km; 0.9 m mean displacement) was unique compared to previous ruptures that had both longer and shorter rupture lengths (∼25−38 km), more displacement (mean of ∼1.3−1.4 m), and equal or greater magnitude ( M w 6.9−7.1) than that in the 1983 earthquake. These ruptures support a hypothesis of variable rupture length and displacement on the northern Lost River fault zone and show that predecessors to the 1983 rupture have passed unimpeded through the Willow Creek Hills. Our work demonstrates that normal faults are capable of producing variable spatial-temporal patterns of rupture that, together with comparisons of fault geometry and historical rupture length, improve our understanding of fault segmentation and help inform models of earthquake rupture probability.
Idaho↗
Portable optically stimulated luminescence age map of a paleoseismic exposure
The quality and quantity of geochronologic data used to constrain the history of major earthquakes in a region exerts a first-order control on the accuracy of seismic hazard assessments that affect millions of people. However, evaluations of geochronological data are limited by uncertainties related to inherently complex depositional processes that may vary spatially and temporally. To improve confidence in models of earthquake timing, we use a high-density suite of radiocarbon and optically stimulated luminescence (OSL) ages with a grid of 342 portable OSL samples to explore spatiotemporal trends in geochronological data across an exemplary normal fault colluvial wedge exposure. The data reveal a two-dimensional age map of the paleoseismic exposure and demonstrate how vertical and horizontal trends in age relate to dominant sedimentary facies and soil characteristics at the site. Portable OSL data provide critical context for the interpretation of 14 C and OSL ages, show that geochronologic age boundaries between pre- and post-earthquake deposits do not match stratigraphic contacts, and provide the basis for selecting alternate Bayesian models of earthquake timing. Our results demonstrate the potential to use emergent, portable OSL methods to dramatically improve paleoseismic constraints on earthquake timing.
Utah↗
Luminescence sediment tracing reveals the complex dynamics of colluvial wedge formation
Paleoearthquake studies that inform seismic hazard rely on assumptions of sediment transport that remain largely untested. Here, we test a widespread conceptual model and a new numerical model on the formation of colluvial wedges, a key deposit used to constrain the timing of paleoearthquakes. We perform this test by applying luminescence, a sunlight-sensitive sediment tracer, at a field site displaying classic colluvial wedge morphostratigraphy. The model and data comparison reveals complex sediment transport processes beyond the predictions of either conceptual or numerical models, including periods of simultaneous debris and wash facies forming processes, erosion, and reworking. These processes could lead to preservation bias, such as incomplete or overinterpretable paleoearthquake records, given the right environmental conditions. Attention to the site-specific mechanics of fault zone depositional systems, such as via sediment tracing, may buffer against the possible effects of preservation bias on paleoseismic study.
Science↗
Landscape evolution in eastern Chuckwalla Valley, Riverside County, California
This study investigates sedimentary and geomorphic processes in eastern Chuckwalla Valley, Riverside County, California, a region of arid, basin-and-range terrain where extensive solar-energy development is planned. The objectives of this study were to (1) measure local weather parameters and use them to model aeolian sediment-transport potential; (2) identify surface sedimentary characteristics in representative localities; and (3) evaluate long-term landscape evolution rates and processes by analyzing stratigraphy in combination with luminescence geochronology. The new stratigraphic and geochronologic data presented in this report demonstrate the varying local significance of aeolian, alluvial fan, lacustrine (playa), and possibly Colorado River influence over a range of time scales. The dominant sand-transport direction in eastern Chuckwalla Valley is toward the northeast, consistent with the recognized regional west-to-east wind direction. However, occasional strong wind events from the north can transport large quantities of sand southward and temporarily reshape local geomorphic features. Influence of a northwest wind direction is also locally dominant around mountain ranges and controls the modern morphology of the Palen dune field. Modeled sand fluxes are on the order of 10 5 kilograms per meter width per year at the site of weather monitoring, 5 kilometers northwest of the Mule Mountains. Aeolian dunes are locally well developed and actively migrating. Their location and activity are determined largely by sediment supply from playa surfaces and ephemeral stream channels, which also control the dunes’ spatial extent and migration potential; stream channels act as both source and sink for aeolian sediment in this environment. Excavations at five sites along a northwest-to-southeast transect reveal that playa deposits formed around 266–226 thousand years ago south of the McCoy Mountains and immediately north of the present location of Interstate 10. The playa material is overlain by late Pleistocene to Holocene alluvial fan deposits. To the southeast (south of Interstate 10, but north of the Mule Mountains), we identified rapid accumulation of alluvial sediment around the time of the Last Glacial Maximum (23–20 thousand years ago), unconformably overlain by a locally varying assemblage of recent aeolian material or Holocene alluvial fan sediment. We have used stratigraphic characteristics and luminescence ages to calculate accumulation rates for sites in eastern Chuckwalla Valley, and thereby to identify spatial variation in landscape stability over decadal and longer time scales. If future solar-energy development plans are to include natural sand-transport corridors, plans would entail retaining the ability for sand to be transported eastward from the ephemeral stream channels and playas that supply sediment to the dunes, sand sheets, and sand ramps of Chuckwalla Valley, and also to allow for southward transport during episodic strong weather events several times per year. The aeolian sediment-transport corridors are dynamic spatially and temporally, reorganizing on the basis of seasonal changes to wind drift potential. Future landscape stability also will be determined by climate-driven changes to vegetation and thereby to aeolian sediment availability. In a warmer, drier climate, aeolian sediment activity is expected to increase, owing to a decrease in stabilizing vegetation cover and more extreme rain that supplies sediment to ephemeral stream channels and playas from which it is remobilized by wind.
California↗