Paleoclimatic implications of glacial fluctuations in the Sierra Nevada del Cocuy, northern Andes, Colombia, during the Lateglacial and Holocene
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Not Available
Vegetation tolerance to drought depends on an array of site-specific environmental and plant physiological factors. This tolerance is poorly understood for many forest types despite its importance for predicting and managing vegetation stress. We analyzed the relationships between precipitation variability and forest die-off in California's Sierra Nevada and introduce a new measure of drought tolerance that emphasizes plant access to subsurface moisture buffers. We applied this metric to California's severe 2012-2015 drought, and show that it predicted the patterns of tree mortality. We then examined future climate scenarios, and found that the probability of droughts that lead to widespread die-off increases threefold by the end of the 21st century. Our analysis shows that tree mortality in the Sierra Nevada will likely accelerate in the coming decades and that forests in the Central and Northern Sierra Nevada that largely escaped mortality in 2012-2015 are vulnerable to die-off.
This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Ton Tonzi Ranch. This is the FLUXNET version of the carbon flux data for the site US-Ton Tonzi Ranch produced by applying the standard ONEFlux (1F) software. Site Description - The Tonzi Ranch is located on the lower foothills of the Sierra Nevada Mountains. The site is classified as an oak savanna woodland (Davis et al, 2016) and is situated on privately owned land. This site is representative of the oak savanna ecosystem that rings the Great Central Valley of California and inhabits the lower reaches of the Sierra Nevada foothills. This ecosystem is a natural laboratory for studying how ecosystem scale carbon and water fluxes respond to environmental change. The temperatures range between 0 and 40 C, the autumn, winter and spring are wet and cool, while the summer growing season is hot, dry and rainless. The overstory is dominated by blue oak trees (40% of total vegetation) with intermittent grey pine trees (3 trees/ha). Understory species include a variety of grasses and herbs, including purple false brome, smooth cat's ear, and rose clover. These two distinctive layers operate in and out of phase from one another. Growing season of the understory is confined to the wet season only, typically from October to early May. In contrast, the deciduous blue oak trees are dormant and leafless during the rainy winter months. The trees tend to leaf out in March and reach maximum LAI in April. During the rainless summers, the grasses in the understory die and the leaves experience soil moisture deficit stress References Davis, F., et al. (2016). "Oak woodlands." Ecosystems of California: 509–529.
In recent decades, there has been a significant increase in annual area burned in California's Sierra Nevada mountains. This rise in fire activity has prompted the need to understand how historical forest management practices affect fuel composition and emissions. Here we examined the total carbon (TC) concentration and radiocarbon abundance (Δ 14 C) of particulate matter (PM) emitted by the KNP Complex Fire, which occurred during California's 2021 wildfire season and affected several groves of giant sequoia trees in the southern Sierra Nevada. During a 26 h sampling period, we measured concentrations of fine airborne PM (PM 2.5 ), as well as dry air mole fractions of carbon monoxide (CO) and methane (CH 4 ), using a ground-based mobile laboratory. We also collected filter samples of PM 2.5 for analysis of TC concentration and Δ 14 C. High correlation among PM 2.5 , CO, and CH 4 time series confirmed that our PM 2.5 measurements captured variability in wildfire emissions. Using a Keeling plot approach, we determined that the mean Δ 14 C of PM 2.5 was 111.6 ± 7.7‰ (n = 12), which was considerably enriched relative to atmospheric carbon dioxide in the northern hemisphere in 2021 (-3.2 ± 1.4‰). Combining these Δ 14 C data with a steady-state one-box ecosystem model, we estimated that the mean age of fuels combusted in the KNP Complex Fire was 40 years, with a range of 29–57 years. These results provide evidence for emissions originating from woody biomass, larger-diameter fine fuels, and coarse woody debris that have accumulated over multiple decades. This is consistent with independent field observations that indicate high fire intensity contributed to widespread giant sequoia mortality. With the expanded use of prescribed fires planned over the next decade in California to mitigate wildfire impacts, our measurement approach has the potential to provide regionally-integrated estimates of the effectiveness of fuel treatment programs.
Abstract The serial occurrence of atmospheric rivers (ARs) along the US West Coast can lead to prolonged and exacerbated hydrologic impacts, threatening flood‐control and water‐supply infrastructure due to soil saturation and diminished recovery time between storms. Here a statistical approach for quantifying subseasonal temporal clustering among extreme events is applied to a 41‐year (1979–2019) wintertime AR catalog across the western United States (US). Observed AR occurrence, compared against a randomly distributed AR timeseries with the same average event density, reveals temporal clustering at a greater‐than‐random rate across the western US with a distinct geographical pattern. Compared to the Pacific Northwest, significant AR clusters over the northern Coastal Range of California and Sierra Nevada are more frequent and occur over longer time periods. Clusters along the California Coastal Range typically persist for 2 weeks, are composed of 4–5 ARs per cluster, and account for over 85% of total AR occurrence. Across the northwest Coast‐Cascade Ranges, clusters account for ∼50% of total AR occurrence, typically last 8–10 days, and contain 3–4 individual AR events. Based on precipitation data from a high‐resolution dynamical downscaling of reanalysis, the fractions of total and extreme hourly precipitation attributable to AR clusters are largest along the northern California coast and in the Sierra Nevada. Interannual variability among clusters highlights their importance for determining whether a particular water year is anomalously wet or dry. The mechanisms behind this unusual clustering are unclear and require further research.
The overarching objective of this LDRD project was to determine the vulnerability of deep soil organic carbon (SOC) to warming in the Sierra Nevada region. Deep soils (>30 cm) store more than 70% of global SOC, and increased SOC decomposition and CO 2 emissions caused by warming are potentially large climate change feedbacks. According to the Intergovernmental Panel on Climate Change, temperatures are expected to increase by 4°C by the year 2100, warming the land and underlying soil, and making understanding of how warming will influence deep SOC storage and persistence critical to projecting the land carbon sink. However, uncertainty remains in our process-level understanding and ability to quantify how projected warming will impact the stability of carbon in deep soils. We investigated warming effects on deep (up to 16 m) SOC stability across climate, vegetation, and soil mineralogy gradients in California. We sampled soils from the surface to bedrock (down to 16 meters) at four sites representing vastly different ecosystems across the Northern and Southern Sierra Nevada mountains. This study quantified the response of SOC concentration, distribution, and vulnerability to warming using a soil incubation experiment to warm the whole soil profile and radiocarbon and stable isotopes to assess which pools are vulnerable to loss under warming. Our results refine our understanding of the terrestrial carbon cycle by revealing the vulnerability of deep SOC to future changes in climate and how minerology may influence that vulnerability.
High-resolution chronologies are crucial for paleoenvironmental reconstructions, and are particularly challenging for lacustrine records of terrestrial paleoclimate. Accelerator Mass Spectrometry (AMS) radiocarbon measurement of terrestrial macrofossils is the most common technique for building age models for lake sediment cores, but relies on the presence of terrestrial macrofossils in sediments. In the absence of sufficient macrofossils, pollen concentrates represent a valuable source of dates for building high-resolution chronologies. However, pollen isolation and dating may present several challenges, as has been reported by different authors in previous work over the last few decades. Here we present an improved method for extracting, purifying and radiocarbon-dating pollen concentrates using flow cytometry to improve the extraction efficiency and the purity of the pollen concentrates. Overall, the nature of the sediments and the abundance of the pollen represent major considerations in obtaining enough pollen grains and, consequently, enough carbon to be dated. Further, the complete separation of pollen from other forms of organic matter is required to ensure the accuracy of the dates. We apply the method to surface samples and sediment cores recovered from two contrasting lake basins on the eastern side of the Sierra Nevada (California), and describe the variations that may be used to optimize pollen preparation from a variety of sediments.
This study focuses on quantifying future anthropogenic changes in surface runoff associated with extreme precipitation in California's Sierra Nevada. The method involves driving a land surface model with output from a high resolution regional atmospheric simulation of the most extreme atmospheric rivers (ARs). AR events were selected from an ensemble of global climate model simulations of historical and late 21st century climate under the “high-emission” RCP8.5 scenario. Average precipitation during the future ARs increases by ~25% but a much lower proportion falls as snow. The resulting future runoff increase is dramatic—nearly 50%, reflecting both the precipitation increase and simultaneous conversion of snow to rain. The “double whammy” impact on runoff is largest in the 2,000–2,500 m elevation band, where the snowfall loss and precipitation increase are both especially large. Overall, this huge increase in runoff during the most extreme AR events could present major flood control challenges for the region.
Mountain ecosystems typically serve as carbon (C) sinks. However, studies also suggest that they could be C sources due to climate warming, drought and insect-related mortality, wildfires, and management actions. We applied the Lund-Potsdam-Jena General Ecosystem Simulator (LPJ-GUESS), a process-based dynamic vegetation-ecosystem model, to investigate the role of ecosystem management in C storage under Mediterranean climate over the 21 st century. We modified LPJ-GUESS to include implementing mechanical thinning by vegetation size classes, components, and types along with a new mechanistic fire-occurrence model that accounts for wind speed and lightning ignition. Simulations show that mechanical thinning or prescribed fire performed 5-20 years in advance of a high-severity wildfire reduced direct wildfire C emissions by 38-66 %. Our results also show that long-term management actions repeated every 5-20 years, including thinning relatively small trees (diameters up to 7 inches or ∼178 mm), can maintain stable C levels in the forest and lower dead-fuel amounts. We found that, although prescribed fire mitigated wildfire severity, ecosystem C storage from reduced wildfire emissions can be outweighed by the added emissions from the prescribed fire themselves. Thinning plus removing and sequestering the thinned biomass can ensure that forests act as net C sinks through the end of the 21 st century. However, addition of prescribed fire is needed to reduce understory and lower the projected extent of high-severity wildfire. Achieving the competing goals of reducing wildfire and making the Sierra Nevada long-term C sink can be advanced through carefully coordinated thinning, sequestration of thinned biomass, and prescribed fire.
The impact of shifting disturbance regimes on soil carbon (C) storage is a key uncertainty in global change research. Wildfires in coniferous forests are becoming more frequent in many regions, potentially causing large C emissions. Repeated low-intensity prescribed fires can mitigate wildfire severity, but repeated combustion may decrease soil C unless compensatory responses stabilize soil organic matter. Here, we tested how 30 years of decadal prescribed burning affected C and nitrogen (N) in plants, detritus, and soils in coniferous forests in the Sierra Nevada mountains, USA. Tree basal area and litter stocks were resilient to fire, but fire reduced forest floor C by 77% (-36.4 Mg C/ha). In mineral soils, fire reduced C that was free from minerals by 41% (-4.4 Mg C/ha) but not C associated with minerals, and only in depths ≤ 5 cm. Fire also transformed the properties of remaining mineral soil organic matter by increasing the proportion of C in a pyrogenic form (from 3.2% to 7.5%) and associated with minerals (from 46% to 58%), suggesting the remaining soil C is more resistant to decomposition. Laboratory assays illustrated that fire reduced microbial CO2 respiration rates by 55% and the activity of eight extracellular enzymes that degrade cellulosic and aromatic compounds by 40–66%. Lower decomposition was correlated with lower inorganic N (-49%), especially ammonium, suggesting N availability is coupled with decomposition. The relative increase in forms of soil organic matter that are resistant to decay or stabilized onto mineral surfaces, and the associated decline in decomposition suggest that low-intensity fires may promote mineral soil C storage in pools with long mean residence times in coniferous forests.
The California-Nevada 1997 New Year’s flood was an atmospheric river (AR)-driven rain-on-snow (RoS) event and remains the costliest in their history. The joint occurrence of saturated soils, rainfall, and snowmelt generated inundation throughout northern California-Nevada. Although AR RoS events are projected to occur more frequently with climate change, the warming sensitivity of their flood drivers across scales remains understudied. We leverage the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) to recreate the 1997 New Year’s flood with horizontal grid spacings of 3.5 km across California, with forecast lead times of up to 4 days, and across six warming levels ranging from pre-industrial conditions to +3.5° C. We describe the sensitivity of the flood drivers to warming including AR duration and intensity, precipitation phase, intensity and efficiency, snowpack mass and energy changes, and runoff efficiency. Our findings indicate current levels of climate change negligibly influence the flood drivers. At warming levels ≥ 1.7° C, AR hazard potential increases, snowpack nonlinearly decreases, antecedent soil moisture decreases (except where the snowline retreats), and runoff decreases (except in the southern Sierra Nevada where antecedent snowpack persists). Storm total precipitation increases, but at rates below warming-induced increases in saturation-specific humidity. Warming intensifies short-duration, high-intensity rainfall, particularly where snowfall-to-rainfall transitions occur. This study highlights the nonlinear tradeoffs in 21st-century RoS flood hazards with warming and provides water management and infrastructure investment adaptation considerations.
The 1997 New Year’s flood was among the most devastating floods in the Truckee River watershed located in western Nevada. This event resulted from complex interactions of flood drivers, such as extreme precipitation, wet antecedent watershed conditions, warm temperatures and rapid snowmelt. We leveraged simulated forcings from the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) and a process-based hydrological model to recreate the 1997 New Year’s flood for the Truckee River watershed across four climate warming levels ranging from the current temperatures to + 4° C. For each scenario, we conducted ensemble simulations with the same forcing but with 100 different seasonal watershed antecedent conditions, which were randomly sampled from long-term hydrological simulations. The results show that the 1997 New Year’s flood can be reproduced or exceeded consistently only when the antecedent watershed conditions are wet, specifically when streamflows are above the 75th percentile of the climatological value. There is negligible change in ensemble mean peakflows for Truckee River near Reno; however, there are increases of 18% and 14% under the warming levels of + 3° C and + 4° C, respectively. The increases in peakflows under future climate warming are attributed to wetter antecedent watershed conditions and enhanced snowmelt. Furthermore, the largest increases in peakflows occur at small, high-elevation headwater basins along the Sierra Nevada crest. This study highlights that changes in extreme flood events will result from the complex interplay of multiple flood drivers. It also demonstrates the potential of storyline approaches to analyze future realizations of these extreme events under different climate scenarios.
Previous studies have documented how forests influence snow at fine spatial scales, but none have documented the influence that existing forest-snow variability has on streamflow. Here to test how much forest-controlled snow variability influences streamflow, a tiling parameterization based on classifications from high-resolution (1–3 m) vegetation maps was incorporated into the Distributed Hydrology and Soil Vegetation Model (DHSVM). Within each grid cell (90–150 m), the tiling parameterization simulated forest-snow variability with four independently evolving snowpacks. Each tile had unique radiation conditions to represent conditions underneath the canopy, in exposed areas, and along north- and south-facing forest edges. This tiled parameterization was used to test where and when detailed forest-snow modelling should be considered further and where and when the impacts are too small to be worth the effort. To test this, tiled model simulations of streamflow were compared to non-tiled model simulations in the Sierra Nevada, CA, the Jemez Mountains, NM, and the Eastern Cascades, WA. In Tuolumne, CA, the tiled model simulated little difference in grid cell average SWE, and late-season streamflow decreased by only 3%–4% compared to the non-tiled model. In Jemez, NM, the tiled model decreased late-season streamflow by 18% due to increased sublimation. In Chiwawa, WA, the tiled model increased late-season streamflow by 15% due to high shortwave radiation attenuation and less longwave radiation enhancement from the forest. Furthermore, within the Chiwawa, a substantial silvicultural practice was synthetically implemented to increase the north-facing edge's fractional area. This silvicultural experiment, which used the same fractional forest area in all simulations increased late-season streamflow by 35% compared to tiled model simulations that did not represent forest edges. In conclusion, representing forest-SWE variability had an effect on late-season streamflow in some watersheds but not in others based on the fractional area of the forest edges, forest characteristics, and climate conditions.
Many regions worldwide have experienced increasing wildfire activity in recent years and climate changes are predicted to result in more frequent and severe fires. Reconstruction of past fire activity offers paleoenvironmental context for modern and future burning. Pyrogenic polycyclic aromatic hydrocarbons (PAHs) have been increasingly used as a molecular biomarker for fire occurrence in the paleorecord and offer opportunity for nuanced reconstructions of fire characteristics. A suite of PAHs are produced during combustion, and the emission amount and assemblage is influenced by many variables including fuel type, fire temperature, and oxygen availability. Despite recent advances in understanding the controls and taphonomy of these biomass burning markers, the spatial scale of this proxy is unknown. In this paper, measurements of PAH fluxes preserved in a lake sediment archive from the Sierra Nevada, California were compared with a historical geographic information system dataset of area burned up to 150 km distance from the lake to determine the spatial scales for which these biomarkers are reliable proxies of burning. Comparisons of PAH fluxes with charcoal accumulation rates in the same sediments suggest that pyrogenic particulate transport modulates low to mid-molecular weight PAHs via adsorption. Overall, the results indicate that PAH records integrate a combination of spatial signals of area burned and measurement of individual PAHs may enable cross-scale paleofire reconstructions.
Abstract We compare historical and end‐of‐century temperature and precipitation patterns over California from one dynamically downscaled simulation using the Weather Research and Forecast (WRF) model and two simulations statistically downscaled using Localized Constructed Analogs (LOCA). We uniquely separate causes of differences between dynamically and statistically based future climate projections into differences in historical climate (gridded observations versus regional climate model output) and differences in how these downscaling techniques explicitly handle future climate changes (numerical modeling versus analogs). In these methods, solutions between different downscaling techniques differ more in the future compared to the historical period. Changes projected by LOCA are insensitive to the choice of driving data. Only through dynamical downscaling can we simulate physically consistent regional springtime warming patterns across the Sierra Nevada, while the statistical simulations inherit an unphysical signal from their parent Global Climate Model (GCM) or gridded data. The results of our study clarify why these different techniques produce different outcomes and may also provide guidance on which downscaled products to use for certain impact analyses in California and perhaps other Mediterranean regimes.
Understanding potential response of forest carbon (C) and nutrient storage to warming is important for climate mitigation policies. Unfortunately, those responses are difficult to predict in seasonally dry forests, in part, because ecosystem processes are highly sensitive to both changes in temperature and precipitation. We investigated how warming might alter stocks of C, nitrogen (N), and phosphorus (P) in vegetation and the entire regolith (soil + weathered bedrock or “saprock”) using a space-for-time substitution along a bioclimatic gradient in the Sierra Nevada, California. The pine-oak and mixed-conifer forests between 1,160–2,015 m elevation have more optimal climates (not too dry or hot) for ecosystem productivity, soil weathering, and cycling of essential elements than the oak savannah (405 m) and subalpine forest (2,700 m). We found decreases in overstory vegetation nutrient stocks with decreasing elevation because of enhanced water limitation and greater occurrence of disturbances. Stocks of C, N, and P in the entire regolith peaked at the pine-oak and mixed-conifer forests across the bioclimatic gradient, driven by thicker regolith profiles and greater nutrient input rates. These observations suggest long-term warming will decrease ecosystem nutrient storage at the warmer, transitional pine-oak zone, but will increase nutrient storage at the colder, subalpine zone. Assuming steady-state conditions, we found the mean residence time of ecosystem C decreased with projected rising air temperatures and increased following a major drought event across the bioclimatic gradient. Our study emphasizes potentially elevation-dependent changes in nutrient storage and C persistence with warming in seasonally dry forests.
Well-dated lacustrine records are essential to establish the timing and drivers of regional hydroclimate change. Searles Basin, California, records the depositional history of a fluctuating saline-alkaline lake in the terminal basin of the Owens River system draining the eastern Sierra Nevada. Here, we establish a U-Th chronology for the ~76-m-long SLAPP-SLRS17 core collected in 2017 based on dating of evaporite minerals. Ninety-eight dated samples comprising nine different minerals were evaluated based on stratigraphic, mineralogic, textural, chemical, and reproducibility criteria. After the application of these criteria, a total of 37 dated samples remained as constraints for the age model. A lack of dateable minerals between 145 and 110 ka left the age model unconstrained over the penultimate glacial termination (Termination II). We thus established a tie point between plant wax δD values in the core and a nearby speleothem δ 18 O record at the beginning of the Last Interglacial. We construct a Bayesian age model allowing stratigraphy to inform sedimentation rate inflections. We find that the >210 ka SLAPP-SRLS17 record contains five major units that correspond with prior work. The new dating is broadly consistent with previous efforts but provides more precise age estimates and enables a detailed evaluation of evaporite depositional history. We also offer a substantial revision of the age of the Bottom Mud-Mixed Layer contact, shifting it from ~130 ka to 178 ± 3 ka. The new U-Th chronology documents the timing of mud and salt layers and lays the foundation for climate reconstructions.