ALBERT

Description: Journal of Climate, Ahead of Print. 〈br/〉

Description: Ongoing oblique slip at the Pacific–North America plate boundary in the Salton Trough produced the Imperial Valley (California, USA), a seismically active area with deformation distributed across a complex network of exposed and buried faults. To better understand the shallow crustal structure in this region and the connectivity of faults and seismicity lineaments, we used data primarily from the Salton Seismic Imaging Project to construct a three-dimensional P-wave velocity model down to 8 km depth and a velocity profile to 15 km depth, both at 1 km grid spacing. A V P = 5.65–5.85 km/s layer of possibly metamorphosed sediments within, and crystalline basement outside, the valley is locally as thick as 5 km, but is thickest and deepest in fault zones and near seismicity lineaments, suggesting a causative relationship between the low velocities and faulting. Both seismicity lineaments and surface faults control the structural architecture of the western part of the larger wedge-shaped basin, where two deep subbasins are located. We estimate basement depths, and show that high velocities at shallow depths and possible basement highs characterize the geothermal areas.

Description: We applied cosmogenic 26Al/10Be burial dating to sedimentary deposits of the ancestral Colorado River. We compared cosmogenic burial ages of sediments to the age of an independently well-dated overlying basalt flow at one site, and also applied cosmogenic burial dating to sediments with less precise independent age constraints. All dated gravels yielded old ages that suggest several episodes of sediment burial over the past ~5.3 m.y. Comparison of burial ages to the overlying 4.4 Ma basalt yielded good agreement and suggests that under the most favorable conditions, cosmogenic burial dating can extend back 4–5 m.y. In contrast, results from other sites with more broadly independent age constraints highlight the complexities inherent in burial dating; these complexities arise from unknown and complicated burial histories, insufficient shielding, postburial production of cosmogenic isotopes by muons, and unknown initial 26Al/10Be ratios. Nevertheless, and in spite of the large range of burial ages and large uncertainties, we identify samples that provide reasonable burial age constraints on the depositional history of sediment along the lower ancestral Colorado River. These samples suggest possible sediment deposition and burial at ca. 5.3, 4.7, and 3.6 Ma.Our calculated basinwide erosion rate for sediment transported by the modern Colorado River (~187 mm k.y.-1) is higher than the modern erosion rates inferred from the historic sediment load (80–100 mm k.y.-1). In contrast, basinwide paleo-erosion rates calculated from Pliocene sediments are all under 40 mm k.y.-1 The comparatively lower denudation rates calculated for the Pliocene sediment samples are surprising given that the sampled time intervals include significant Pliocene aggradation and may include much incision of the Grand Canyon and its tributaries. This conflict may arise from extensive storage of sediment along the route of the Colorado River, slower paleobedrock erosion, or the inclusion of sediments that were derived preferentially from higher elevations in the watershed.

Description: We present quantitative analyses of recent large rock falls in Yosemite Valley, California, using integrated high-resolution imaging techniques. Rock falls commonly occur from the glacially sculpted granitic walls of Yosemite Valley, modifying this iconic landscape but also posing significant potential hazards and risks. Two large rock falls occurred from the cliff beneath Glacier Point in eastern Yosemite Valley on 7 and 8 October 2008, causing minor injuries and damaging structures in a developed area. We used a combination of gigapixel photography, airborne laser scanning (ALS) data, and ground-based terrestrial laser scanning (TLS) data to characterize the rock-fall detachment surface and adjacent cliff area, quantify the rock-fall volume, evaluate the geologic structure that contributed to failure, and assess the likely failure mode. We merged the ALS and TLS data to resolve the complex, vertical to overhanging topography of the Glacier Point area in three dimensions, and integrated these data with gigapixel photographs to fully image the cliff face in high resolution. Three-dimensional analysis of repeat TLS data reveals that the cumulative failure consisted of a near-planar rock slab with a maximum length of 69.0 m, a mean thickness of 2.1 m, a detachment surface area of 2750 m2, and a volume of 5663 {+/-} 36 m3. Failure occurred along a surface-parallel, vertically oriented sheeting joint in a clear example of granitic exfoliation. Stress concentration at crack tips likely propagated fractures through the partially attached slab, leading to failure. Our results demonstrate the utility of high-resolution imaging techniques for quantifying far-range (〉1 km) rock falls occurring from the largely inaccessible, vertical rock faces of Yosemite Valley, and for providing highly accurate and precise data needed for rock-fall hazard assessment.

Description: The magnitude of late Cenozoic rock uplift of the Sierra Nevada (California, USA) remains unresolved despite more than a century of investigation, with estimates ranging from essentially zero to ~3 km of uplift at the range crest. Two sets of two-dimensional end-member mechanical models bracket how normal faulting along the eastern escarpment of the Sierra Nevada contributed to uplift of the range over a time span of millions of years. The short-term models are based on dislocations in an elastic half-space. The long-term models involve thin elastic beams resting on an inviscid fluid. Both sets of models predict that if the regional topography were entirely a response to faulting along the eastern escarpment, then the bedrock floors immediately east of the range should consistently lie thousands of meters below sea level, instead of thousands of meters above sea level as they generally do. Both sets of analyses indicate that although faulting would lift the range crest, it would drop the rock east of the range-front faults at least as much, and perhaps much more; model results suggest that ~66%–85% of the current escarpment relief stems from subsidence of the grabens east of the Sierra Nevada, with only ~15%–34% resulting from crestal uplift. Our results strongly indicate that range-front faulting in the last 3–10 m.y. uplifted rock at the Sierra Nevada crest by hundreds of meters to as much as 1 km, and that this uplift was superposed on high topography that predated the origin of the eastern escarpment. These conclusions are compatible with diverse geologic observations and measurements.

Description: The 4 April 2010 moment magnitude (M w ) 7.2 El Mayor–Cucapah earthquake revealed the existence of a previously unidentified fault system in Mexico that extends ~120 km from the northern tip of the Gulf of California to the U.S.–Mexico border. The system strikes northwest and is composed of at least seven major faults linked by numerous smaller faults, making this one of the most complex surface ruptures ever documented along the Pacific–North America plate boundary. Rupture propagated bilaterally through three distinct kinematic and geomorphic domains. Southeast of the epicenter, a broad region of distributed fracturing, liquefaction, and discontinuous fault rupture was controlled by a buried, southwest-dipping, dextral-normal fault system that extends ~53 km across the southern Colorado River delta. Northwest of the epicenter, the sense of vertical slip reverses as rupture propagated through multiple strands of an imbricate stack of east-dipping dextral-normal faults that extend ~55 km through the Sierra Cucapah. However, some coseismic slip (10–30 cm) was partitioned onto the west-dipping Laguna Salada fault, which extends parallel to the main rupture and defines the western margin of the Sierra Cucapah. In the northernmost domain, rupture terminates on a series of several north-northeast–striking cross-faults with minor offset (〈8 cm) that cut uplifted and folded sediments of the northern Colorado River delta in the Yuha Desert. In the Sierra Cucapah, primary rupture occurred on four major faults separated by one fault branch and two accommodation zones. The accommodation zones are distributed in a left-stepping en echelon geometry, such that rupture passed systematically to structurally lower faults. The structurally lowest fault that ruptured in this event is inclined as shallowly as ~20°. Net surface offsets in the Sierra Cucapah average ~200 cm, with some reaching 300–400 cm, and rupture kinematics vary greatly along strike. Nonetheless, instantaneous extension directions are consistently oriented ~085° and the dominant slip direction is ~310°, which is slightly (~10°) more westerly than the expected azimuth of relative plate motion, but considerably more oblique to other nearby historical ruptures such as the 1992 Landers earthquake. Complex multifault ruptures are common in the central portion of the Pacific North American plate margin, which is affected by restraining bend tectonics, gravitational potential energy gradients, and the inherently three-dimensional strain of the transtensional and transpressional shear regimes that operate in this region.

Description: Erosion of volcanic ocean islands creates dramatic landscapes, modulates Earth’s carbon cycle, and delivers sediment to coasts and reefs. Because many volcanic islands have large climate gradients and minimal variations in lithology and tectonic history, they are excellent natural laboratories for studying climatic effects on the evolution of topography. Despite concerns that modern sediment fluxes to island coasts may exceed long-term fluxes, little is known about how erosion rates and processes vary across island interiors, how erosion rates are influenced by the strong climate gradients on many islands, and how modern island erosion rates compare to long-term rates. Here, we present new measurements of erosion rates over 5 yr to 5 m.y. timescales on the Hawaiian island of Kaua‘i, across which mean annual precipitation ranges from 0.5 to 9.5 m/yr. Eroded rock volumes from basins across Kaua‘i indicate that million-year-scale erosion rates are correlated with modern mean annual precipitation and range from 8 to 335 t km –2 yr –1 . In Kaua‘i’s Hanalei River basin, 3 He concentrations in detrital olivines imply millennial-scale erosion rates of 〉126 to 〉390 t km –2 yr –1 from olivine-bearing hillslopes, while fluvial suspended sediment fluxes measured from 2004 to 2009 plus estimates of chemical and bed-load fluxes imply basin-averaged erosion rates of 545 ± 128 t km –2 yr –1 . Mapping of landslide scars in satellite imagery of the Hanalei basin from 2004 and 2010 implies landslide-driven erosion rates of 30–47 t km –2 yr –1 . These measurements imply that modern erosion rates in the Hanalei basin are no more than 2.3 ± 0.6 times faster than millennial-scale erosion rates, and, to the extent that modern precipitation patterns resemble long-term patterns, they are consistent with a link between precipitation rates and long-term erosion rates.

Description: Anthropogenic climate change is projected to lead to ocean warming, acidification, deoxygenation, reductions in near-surface nutrients, and changes to primary production, all of which are expected to affect marine ecosystems. Here we assess projections of these drivers of environmental change over the twenty-first century from Earth system models (ESMs) participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6) that were forced under the CMIP6 Shared Socioeconomic Pathways (SSPs). Projections are compared to those from the previous generation (CMIP5) forced under the Representative Concentration Pathways (RCPs). A total of 10 CMIP5 and 13 CMIP6 models are used in the two multi-model ensembles. Under the high-emission scenario SSP5-8.5, the multi-model global mean change (2080–2099 mean values relative to 1870–1899) ± the inter-model SD in sea surface temperature, surface pH, subsurface (100–600 m) oxygen concentration, euphotic (0–100 m) nitrate concentration, and depth-integrated primary production is +3.47±0.78 ∘C, -0.44±0.005, -13.27±5.28, -1.06±0.45 mmol m−3 and -2.99±9.11 %, respectively. Under the low-emission, high-mitigation scenario SSP1-2.6, the corresponding global changes are +1.42±0.32 ∘C, -0.16±0.002, -6.36±2.92, -0.52±0.23 mmol m−3, and -0.56±4.12 %. Projected exposure of the marine ecosystem to these drivers of ocean change depends largely on the extent of future emissions, consistent with previous studies. The ESMs in CMIP6 generally project greater warming, acidification, deoxygenation, and nitrate reductions but lesser primary production declines than those from CMIP5 under comparable radiative forcing. The increased projected ocean warming results from a general increase in the climate sensitivity of CMIP6 models relative to those of CMIP5. This enhanced warming increases upper-ocean stratification in CMIP6 projections, which contributes to greater reductions in upper-ocean nitrate and subsurface oxygen ventilation. The greater surface acidification in CMIP6 is primarily a consequence of the SSPs having higher associated atmospheric CO2 concentrations than their RCP analogues for the same radiative forcing. We find no consistent reduction in inter-model uncertainties, and even an increase in net primary production inter-model uncertainties in CMIP6, as compared to CMIP5.

Description: The coastal ecosystem of the Gulf of Alaska (GOA) is especially vulnerable to the effects of ocean acidification and climate change. Detection of these long-term trends requires a good understanding of the system’s natural state. The GOA is a highly dynamic system that exhibits large inorganic carbon variability on subseasonal to interannual timescales. This variability is poorly understood due to the lack of observations in this expansive and remote region. We developed a new model setup for the GOA that couples the three-dimensional Regional Oceanic Model System (ROMS) and the Carbon, Ocean Biogeochemistry and Lower Trophic (COBALT) ecosystem model. To improve our conceptual understanding of the system, we conducted a hindcast simulation from 1980 to 2013. The model was explicitly forced with temporally and spatially varying coastal freshwater discharges from a high-resolution terrestrial hydrological model, thereby affecting salinity, alkalinity, dissolved inorganic carbon, and nutrient concentrations. This represents a substantial improvement over previous GOA modeling attempts. Here, we evaluate the model on seasonal to interannual timescales using the best available inorganic carbon observations. The model was particularly successful in reproducing observed aragonite oversaturation and undersaturation of near-bottom water in May and September, respectively. The largest deficiency in the model is its inability to adequately simulate springtime surface inorganic carbon chemistry, as it overestimates surface dissolved inorganic carbon, which translates into an underestimation of the surface aragonite saturation state at this time. We also use the model to describe the seasonal cycle and drivers of inorganic carbon parameters along the Seward Line transect in under-sampled months. Model output suggests that the majority of the near-bottom water along the Seward Line is seasonally undersaturated with respect to aragonite between June and January, as a result of upwelling and remineralization. Such an extensive period of reoccurring aragonite undersaturation may be harmful to ocean acidification-sensitive organisms. Furthermore, the influence of freshwater not only decreases the aragonite saturation state in coastal surface waters in summer and fall, but it simultaneously decreases the surface partial pressure of carbon dioxide (pCO2), thereby decoupling the aragonite saturation state from pCO2. The full seasonal cycle and geographic extent of the GOA region is under-sampled, and our model results give new and important insights for months of the year and areas that lack in situ inorganic carbon observations.