Search USGSSearch

SEARCH · Search USGS

Results for “Climate of the Past”

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 1,189 records · Page 66Linked to original sources

Imprint of regional oceanography on foraminifera of eastern Pacific Coral Reefs

The marginal marine environments of the eastern tropical Pacific (ETP) serve as an ideal natural laboratory to study how oceanographic and climatic variability influence coral-reef ecosystems. Reefs along the Pacific coast of Panamá span a natural gradient of nutrients, pH, and temperature as a result of stronger seasonal upwelling in the Gulf of Panamá relative to the Gulf of Chiriquí. The ecosystems are not only influenced by spatial and seasonal variations in oceanography but are affected by the climatic variability of the El Niño-Southern Oscillation (ENSO). Foraminifera can be robust indicators of ecosystem condition because the composition of their assemblages and the geochemistry of their tests can change rapidly in response to environmental variability. We studied benthic foraminifera in sediment samples collected from 3 m below mean sea level in the Gulf of Panamá and the Gulf of Chiriquí. Temperature loggers deployed from 2016 to 2019 showed that average temperatures were lower and more variable in the Gulf of Panamá due to seasonal upwelling. All sites in both gulfs were dominated by heterotrophic foraminifera, which was likely the result of nutrient enrichment due to upwelling, combined with ENSO effects. However, the Gulf of Chiriquí was characterized by higher abundances of symbiont-bearing foraminifera than the Gulf of Panamá. The orders Miliolida and Rotaliida dominated the foraminiferal assemblages in both gulfs, with Quinqueloculina and Rosalina being the most abundant genera in the two orders, respectively. Miliolids were less abundant in the Gulf of Panamá than in the Gulf of Chiriquí, whereas rotaliid densities were not significantly different between the two gulfs. Lower pH in the Gulf of Panamá as a result of upwelling may have contributed to the lower abundance of miliolids, which secrete tests of high-magnesium calcite. Geochemical analysis of tests of the symbiont-bearing miliolid Sorites marginalis revealed that foraminiferal Mg/Ca ratios were lower in the Gulf of Panamá than in the Gulf of Chiriquí. The offset in foraminiferal Mg/Ca is consistent with the lower mean annual temperature observed in the Gulf of Panamá due to stronger seasonal upwelling. Because the geochemistry and assemblages of foraminifera reflect differences in environmental conditions, they could potentially be used in tandem with coral proxies to reconstruct past environmental change and project the future of coral-reef systems within the ETP.

Gulf of Chiriqui, Gulf of Panama

Mapping burned areas using dense time-series of Landsat data

Complete and accurate burned area data are needed to document patterns of fires, to quantify relationships between the patterns and drivers of fire occurrence, and to assess the impacts of fires on human and natural systems. Unfortunately, in many areas existing fire occurrence datasets are known to be incomplete. Consequently, the need to systematically collect burned area information has been recognized by the United Nations Framework Convention on Climate Change and the Intergovernmental Panel on Climate Change, which have both called for the production of essential climate variables (ECVs), including information about burned area. In this paper, we present an algorithm that identifies burned areas in dense time-series of Landsat data to produce the Landsat Burned Area Essential Climate Variable (BAECV) products. The algorithm uses gradient boosted regression models to generate burn probability surfaces using band values and spectral indices from individual Landsat scenes, lagged reference conditions, and change metrics between the scene and reference predictors. Burn classifications are generated from the burn probability surfaces using pixel-level thresholding in combination with a region growing process. The algorithm can be applied anywhere Landsat and training data are available. For this study, BAECV products were generated for the conterminous United States from 1984 through 2015. These products consist of pixel-level burn probabilities for each Landsat scene, in addition to, annual composites including: the maximum burn probability and a burn classification. We compared the BAECV burn classification products to the existing Global Fire Emissions Database (GFED; 1997–2015) and Monitoring Trends in Burn Severity (MTBS; 1984–2013) data. We found that the BAECV products mapped 36% more burned area than the GFED and 116% more burned area than MTBS. Differences between the BAECV products and the GFED were especially high in the West and East where the BAECV products mapped 32% and 88% more burned area, respectively. However, the BAECV products found less burned area than the GFED in regions with frequent agricultural fires. Compared to the MTBS data, the BAECV products identified 31% more burned area in the West, 312% more in the Great Plains, and 233% more in the East. Most pixels in the MTBS data were detected by the BAECV, regardless of burn severity. The BAECV products document patterns of fire similar to those in the GFED but also showed patterns of fire that are not well characterized by the existing MTBS data. We anticipate the BAECV products will be useful to studies that seek to understand past patterns of fire occurrence, the drivers that created them, and the impacts fires have on natural and human systems.

continental United States

The unprecedented loss of Florida's reef-building corals and the emergence of a novel coral-reef assemblage

Over the last half century, climate change, coral disease, and other anthropogenic disturbances have restructured coral-reef ecosystems on a global scale. The disproportionate loss of once-dominant, reef-building taxa has facilitated relative increases in the abundance of “weedy” or stress-tolerant coral species. Although the recent transformation of coral-reef assemblages is unprecedented on ecological timescales, determining whether modern coral reefs have truly reached a novel ecosystem state requires evaluating the dynamics of reef composition over much longer periods of time. Here, we provide a geologic perspective on the shifting composition of Florida's reefs by reconstructing the millennial-scale spatial and temporal variability in reef assemblages using 59 Holocene reef cores collected throughout the Florida Keys Reef Tract (FKRT). We then compare the relative abundances of reef-building species in the Holocene reef framework to data from contemporary reef surveys to determine how much Florida's modern reef assemblages have diverged from long-term baselines. We show that the composition of Florida's reefs was, until recently, remarkably stable over the last 8000 yr. The same corals that have dominated shallow-water reefs throughout the western Atlantic for hundreds of thousands of years, Acropora palmata , Orbicella spp., and other massive coral taxa, accounted for nearly 90% of Florida's Holocene reef framework. In contrast, the species that now have the highest relative abundances on the FKRT, primarily Porites astreoides and Siderastrea siderea , were rare in the reef framework, suggesting that recent shifts in species assemblages are unprecedented over millennial timescales. Although it may not be possible to return coral reefs to pre-Anthropocene states, our results suggest that coral-reef management focused on the conservation and restoration of the reef-building species of the past, will optimize efforts to preserve coral reefs, and the valuable ecosystem services they provide into the future.

Florida

Threshold effects of flood duration on the vegetation and soils of the Upper Mississippi River floodplain, USA

Most large rivers have experienced major changes in hydrology and land use over the past century, with concomitant effects on sedimentation, nutrient cycling and biodiversity. To restore and/or enhance these ecosystems, managers need to know where their efforts are most likely to succeed under current hydrologic regimes as well as under potential future hydrologic regimes. We therefore examined changes in forest vegetation and soils across a hydrologic gradient, expressed as flood duration during the growing season, for 320 km of the Upper Mississippi River (UMR) floodplain. Soil texture was highly variable but trended toward finer grained sediments and >5% organic matter as flood duration increased from 0% to ~40% of the growing season. Beyond 40%, soil texture was exclusively silt plus clay with >5% organic matter. The diversity of both the understory and overstory tree communities was also highly variable at sites that flooded for <40% of the growing season. However, understory diversity decreased as flood duration increased from 0% to ~25% of the growing season and overstory diversity declined as flood duration increased from 0% to ~40% of the growing season. Diversity estimates for both strata were uniformly low at sites that flooded for longer than ~40% of the growing season. Beyond this point the proportional abundance of Acer saccharinum in the overstory exceeded 70%. Our results suggest that there is a threshold along the elevation gradient of this floodplain, corresponding with flood durations lasting ~40% of the growing season. At lower elevation sites, flooding exerts primary control over forest soils and vegetation, restricting the former to silt plus clay with higher organic matter and the latter to a few highly flood tolerant species. The existence of such thresholds have implications for management of floodplain soil nutrient dynamics and plant diversity under existing hydrologic regimes, more natural hydrologic regimes and more extreme hydrologic regimes that may result from climate change.

Iowa;Minnesota;Wisconsin

The mineral diversity of Jezero crater: Evidence for possible lacustrine carbonates on Mars

Noachian-aged Jezero crater is the only known location on Mars where clear orbital detections of carbonates are found in close proximity to clear fluvio-lacustrine features indicating the past presence of a paleolake; however, it is unclear whether or not the carbonates in Jezero are related to the lacustrine activity. This distinction is critical for evaluating the astrobiological potential of the site, as lacustrine carbonates on Earth are capable of preserving biosignatures at scales that may be detectable by a landed mission like the Mars 2020 rover, which is planned to land in Jezero in February 2021. In this study, we conduct a detailed investigation of the mineralogical and morphological properties of geological units within Jezero crater in order to better constrain the origin of carbonates in the basin and their timing relative to fluvio-lacustrine activity. Using orbital visible/near-infrared hyperspectral images from the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) along with high resolution imagery and digital elevation models, we identify a distinct carbonate-bearing unit, the “Marginal Carbonates,” located along the inner margin of the crater, near the largest inlet valley and the western delta. Based on their strong carbonate signatures, topographic properties, and location in the crater, we propose that this unit may preserve authigenic lacustrine carbonates, precipitated in the near-shore environment of the Jezero paleolake. Comparison to carbonate deposits from terrestrial closed basin lakes suggests that if the Marginal Carbonates are lacustrine in origin, they could preserve macro- and microscopic biosignatures in microbialite rocks like stromatolites, some of which would likely be detectable by Mars 2020. The Marginal Carbonates may represent just one phase of a complex fluvio-lacustrine history in Jezero crater, as we find that the spectral diversity of the fluvio-lacustrine deposits in the crater is consistent with a long-lived lake system cataloging the deposition and erosion of regional geologic units. Thus, Jezero crater may contain a unique record of the evolution of surface environments, climates, and habitability on early Mars.

Icarus

Deglaciation and postglacial treeline fluctuation in the northern San Juan Mountains, Colorado

The San Juan Mountains of southwestern Colorado contain numerous lakes and bogs at and above treeline. In June 1978, Lake Emma, a tarn above present-day treeline, was suddenly drained by the collapse of underground mine workings. This study was initiated because the draining exposed a well-preserved archive of subfossil coniferous wood fragments that provided a unique opportunity to further our understanding of the paleoclimatic history of this region. These paleoclimatic studies-coniferous macrofossil identification in conjunction with radiocarbon dating, deuterium analysis of the dated conifer fragments, as well as pollen and fossil insect analyses-yielded new information regarding Holocene climate and accompanying treeline changes in the northern San Juan Mountains. This report synthesizes previously published reports by the author and other investigators, and unpublished information of the author bearing on late Pleistocene and Holocene treeline and climate in this region. Retreat of the glacier that occupied the upper Animas River valley from its Pinedale terminal position began about 19.4 + or - 1.5 10 Be thousands of years ago and was essentially complete by about 12.3 + or - 1.0 10 Be thousands of years ago. Two sets of late Pleistocene cirque moraines were identified in the northern San Juan Mountains. The older set is widespread and probably correlates with the Younger Dryas (11,000-10,000 radiocarbon years before present; 12,800-11,500 calendar years). The younger set is found only in the Grenadier Range and represents remnant glacier ice lying in well-shaded niches in a mountain range undergoing rapid deglaciation. A snowbank at the northern base of this range appears to be fronted by a Little Ice Age moraine. Soon after deglaciation the average July temperature is estimated to have been about 5&deg;C cooler and timberline about 650 meters lower than at present. However, timberline (and treeline) responded rapidly to the postglacial warming and reached higher-than-present elevations by the early Holocene. A comparison of recently obtained accelerator mass spectrometry radiocarbon ages of coniferous wood fragments from Lake Emma, previously dated by conventional radiocarbon methods during the 1980s, led to a slight modification of previously published ages of Holocene treeline fluctuations. As early as 9,200 radiocarbon years before present (about 10,400 calendar years) and probably to about 5,400 radiocarbon years before present (about 6,200 calendar years), treeline was at least 80 meters higher than at present (about 3,660 meters). Furthermore, a large conifer fragment with a complacent annual ring record suggests that timberline may have been at least 140 meters higher than present (about 3,600 meters) about 8,000 radiocarbon years before present (about 8,900 calendar years). These past elevations of treeline and timberline suggest that growing-season temperatures were at least 0.5&deg;-0.9&deg;C warmer than at present. Deuterium data from the Lake Emma wood samples suggests that the maximum average temperature change from about 9,000 to 5,400 radiocarbon years before present (about 10,150 to 6,200 calendar years) was about 4&deg;C. Owing to these warmer temperatures the summer monsoon circulation, which currently brings a large part of the annual precipitation to the San Juan Mountains, probably was more intense during the early and middle Holocene than it is today. Between about 5,400 and 3,500 radiocarbon years before present (about 6,200 and 3,770 calendar years) it appears that treeline was near its present-day limit. After 3,500 radiocarbon years before present (about 3,770 calendar years), evidence of treeline position is very sparse, suggesting that treeline lay at, or below, its present-day elevation. However, a spruce krummholz fragment from the Lake Emma site provided two radiocarbon ages of about 3,100 radiocarbon years before present (about 3,300 cal

Colorado

Magmatic-tectonic settings of Cenozoic epithermal gold-silver deposits of the Great Basin, western United States

Numerous epithermal gold-silver deposits formed during the past 40 Ma are irregularly distributed across the Great Basin. These deposits formed in six major magmatic-tectonic settings that varied during the complex evolution of the continental margin of western North America: (1) slab rollback–ignimbrite flareup (~45–17 Ma), (2) slab rollback–ancestral Cascade arc (~35 Ma–present), (3) Yellowstone hotspot– bimodal (~16.7–3 Ma), (4) slab window (~16 Ma–present), (5) Basin and Range bimodal extensional (≤ 7 Ma), and (6) amagmatic extensional (≤ 5 Ma). Most large (> 1 Moz gold produced) deposits are Miocene (~20–8 Ma), low-, intermediate-, and high-sulfidation deposits in the southern part of the ancestral Cascade arc; late Miocene post-subduction, low-sulfidation deposits formed over the slab window; and lowsulfidation deposits related to early (16.7–15 Ma) Yellowstone hotspot magmatism formed along the northern Nevada rift and related fracture zones to the west. The world-class Round Mountain low-sulfidation deposit is the only large deposit in ignimbrite flareup rocks despite these rocks constituting the largest eruptive volume of Cenozoic magmas in the Great Basin. Intermediate to silicic composition lava dome complexes are the most common setting for epithermal deposits in the western Great Basin, whereas deposits formed in a wide range of settings and rock types during Yellowstone hotspot activity. With exception of the Round Mountain caldera, the dozens of calderas of the ignimbrite flareup do not host large epithermal deposits. Several young (≤ 5 Ma), “amagmatic” low-sulfidation deposits formed along Basin and Range fault zones in sedimentary rocks that lack proximal magmatic activity. The types and characteristics of epithermal gold-silver deposits in the Great Basin systematically vary with magmatic-tectonic setting and magma composition, and their distribution reflects the combined effects of tectonic setting of magma genesis; magma source, composition and eruptive style; crustal thickness and composition; presence of crustal-scale structural zones; climate; and preservation of deposits.

Great Basin

Late Holocene estuarine-inner shelf interactions; is there evidence of an estuarine retreat path for Tampa Bay, Florida?

The purpose of this study was to determine if and how a large, modern estuarine system, situated in the middle of an ancient carbonate platform, has affected its adjacent inner shelf both in the past during the last, post-glacial sea-level rise and during the present. An additional purpose was to determine if and how this inner shelf seaward of a major estuary differed from the inner shelves located just to the north and south but seaward of barrier-island shorelines. Through side-scan sonar mosaicking, bathymetric studies, and ground-truthing using surface grab samples as well as diver observations, two large submarine sand plains were mapped - one being the modern ebb-tidal delta and the other interpreted to be a relict ebb-tidal delta formed earlier in the Holocene. The most seaward portion of the inner shelf studied consists of a field of lobate, bathymetrically elevated, fine-sand accumulations, which were interpreted to be sediment-starved 3D dunes surrounded by small 2D dunes composed of coarse molluscan shell gravel. Additionally, exposed limestone hardbottoms supporting living benthic communities were found as well. This modern shelf sedimentary environment is situated on a large, buried shelf valley, which extends eastward beneath the modern Tampa Bay estuary. These observations plus the absence of an incised shelf valley having surficial bathymetric expression, and the absence of sand bodies normally associated with back-tracking estuarine systems indicate that there was no cross-shelf estuarine retreat path formed during the last rise in sea level. Instead, the modern Tampa Bay formed within a mid-platform, low-relief depression, which was flooded by rising marine waters late in the Holocene. With continued sea-level rise in the late Holocene, this early embayment was translated eastward or landward to its present position, whereby a larger ebb-tidal delta prograded out onto the inner shelf. Extensive linear sand ridges, common to the inner shelves to the north and south, did not form in this shelf province because it was a low-energy, open embayment lacking the wave climate and nearshore zone necessary to create such sand bodies. The distribution of bedforms on the inner shelf and the absence of seaward-oriented 2D dunes on the modern ebb-tidal delta indicate that the modern estuarine system has had little effect on its adjacent inner shelf. ?? 2003 Elsevier B.V. All rights reserved.

Marine Geology

Our evolving conceptual model of the coastal eutrophication problem

A primary focus of coastal science during the past 3 decades has been the question: How does anthropogenic nutrient enrichment cause change in the structure or function of nearshore coastal ecosystems? This theme of environmental science is recent, so our conceptual model of the coastal eutrophication problem continues to change rapidly. In this review, I suggest that the early (Phase I) conceptual model was strongly influenced by limnologists, who began intense study of lake eutrophication by the 1960s. The Phase I model emphasized changing nutrient input as a signal, and responses to that signal as increased phytoplankton biomass and primary production, decomposition of phytoplankton-derived organic matter, and enhanced depletion of oxygen from bottom waters. Coastal research in recent decades has identified key differences in the responses of lakes and coastal-estuarine ecosystems to nutrient enrichment. The contemporary (Phase II) conceptual model reflects those differences and includes explicit recognition of (1) system-specific attributes that act as a filter to modulate the responses to enrichment (leading to large differences among estuarine-coastal systems in their sensitivity to nutrient enrichment); and (2) a complex suite of direct and indirect responses including linked changes in: water transparency, distribution of vascular plants and biomass of macroalgae, sediment biogeochemistry and nutrient cycling, nutrient ratios and their regulation of phytoplankton community composition, frequency of toxic/harmful algal blooms, habitat quality for metazoans, reproduction/growth/survival of pelagic and benthic invertebrates, and subtle changes such as shifts in the seasonality of ecosystem functions. Each aspect of the Phase II model is illustrated here with examples from coastal ecosystems around the world. In the last section of this review I present one vision of the next (Phase III) stage in the evolution of our conceptual model, organized around 5 questions that will guide coastal science in the early 21st century: (1) How do system-specific attributes constrain or amplify the responses of coastal ecosystems to nutrient enrichment? (2) How does nutrient enrichment interact with other stressors (toxic contaminants, fishing harvest, aquaculture, nonindigenous species, habitat loss, climate change, hydrologic manipulations) to change coastal ecosystems? (3) How are responses to multiple stressors linked? (4) How does human-induced change in the coastal zone impact the Earth system as habitat for humanity and other species? (5) How can a deeper scientific understanding of the coastal eutrophication problem be applied to develop tools for building strategies at ecosystem restoration or rehabilitation?

Marine Ecology Progress Series

Forecasting selenium discharges to the San Francisco Bay-Delta Estuary: Ecological effects of a proposed San Luis Drain extension

Selenium discharges to the San Francisco Bay-Delta Estuary (Bay-Delta) could change significantly if federal and state agencies (1) approve an extension of the San Luis Drain to convey agricultural drainage from the western San Joaquin Valley to the North Bay (Suisun Bay, Carquinez Strait, and San Pablo Bay); (2) allow changes in flow patterns of the lower San Joaquin River and Bay-Delta while using an existing portion of the San Luis Drain to convey agricultural drainage to a tributary of the San Joaquin River; or (3) revise selenium criteria for the protection of aquatic life or issue criteria for the protection of wildlife. Understanding the biotransfer of selenium is essential to evaluating effects of selenium on Bay-Delta ecosystems. Confusion about selenium threats to fish and wildlife stem from (1) monitoring programs that do not address specific protocols necessary for an element that bioaccumulates; and (2) failure to consider the full complexity of the processes that result in selenium toxicity. Past studies show that predators are more at risk from selenium contamination than their prey, making it difficult to use traditional methods to predict risk from environmental concentrations alone. This report presents an approach to conceptualize and model the fate and effects of selenium under various load scenarios from the San Joaquin Valley. For each potential load, progressive forecasts show resulting (1) water-column concentration; (2) speciation; (3) transformation to particulate form; (4) particulate concentration; (5) bioaccumulation by invertebrates; (6) trophic transfer to predators; and (7) effects on those predators. Enough is known to establish a first-order understanding of relevant conditions, biological response, and ecological risks should selenium be discharged directly into the North Bay through a conveyance such as a proposed extension of the San Luis Drain. The approach presented here, the Bay-Delta selenium model, determines the mass, fate, and effects of selenium released to the Bay-Delta through use of (1) historical land-use, drainage, alluvial-fill, and runoff databases; (2) existing knowledge concerning biogeochemical reactions and physiological parameters of selenium (e.g., speciation, partitioning between dissolved and particulate forms, and bivalve assimilation efficiency); and (3) site-specific data mainly from 1986 to 1996 for clams and bottom-feeding fish and birds. Selenium load scenarios consider effluents from North Bay oil refineries and discharges of agricultural drainage from the San Joaquin Valley to enable calculation of (a) a composite freshwater endmember selenium concentration at the head of the estuary; and (b) a selenium concentration at a selected seawater location (Carquinez Strait) as a foundation for modeling. Analysis of selenium effects also takes into account the mode of conveyance for agricultural drainage (i.e., the San Luis Drain or San Joaquin River); and flows of the Sacramento River and San Joaquin River on a seasonal or monthly basis. Load scenarios for San Joaquin Valley mirror predictions made since 1955 of a worsening salt (and by inference, selenium) build-up exacerbated by an arid climate and massive irrigation. The reservoir of selenium in the San Joaquin Valley is sufficient to provide loading at an annual rate of approximately 42,500 pounds of selenium to a Bay-Delta disposal point for 63 to 304 years at the lower range of projections presented here, even if influx of selenium from the California Coast Ranges could be curtailed. Disposal of wastewaters on an annual basis outside of the San Joaquin Valley may slow the degradation of valley resources, but drainage alone cannot alleviate the salt and selenium build-up in the San Joaquin Valley, at least within a century. Load scenarios also show the different proportions of selenium loading to the Bay-Delta. Oil refinery loads from 1986 to 1992 ranged from 8.5 to 20 pounds of selenium per day; with treatment and cleanup, loads decreased to 3.0 pounds of selenium per day in 1999. In contrast, San Joaquin Valley agricultural drainage loads disposed of in a San Luis Drain extension could range from 45 to 117 pounds of selenium per day across a set of historical and future conditions. Components of this valley-wide load include five source subareas (i.e., Grassland, Westlands, Tulare, Kern, and Northern) defined by water and drainage management. Loads vary per subarea mainly because of proximity of the subarea to geologic sources of selenium and irrigation history. Loads from the Sacramento River, depending on flow conditions, range from 0.8 to 10 pounds of selenium per day. Loads from the San Joaquin River vary depending on restoration and flow conditions, which are considered. A consistent picture of ecological risk emerges under modeled selenium discharges from a proposed San Luis Drain extension. The threat to the estuary is greatest during low flow seasons and critically dry years. Where selenium undergoes reactions typical of low flow or longer residence time, highly problematic bioaccumulation in prey (food) is forecast. Surf scoter, greater and lesser scaup, and white sturgeon appear to be most at risk because these Bay-Delta predators feed on deposit and filter-feeding bivalves. Recent findings add Sacramento splittail and Dungeness crab to that list. During the low flow season of critically dry years, forecasted selenium concentrations in water, particulate matter, prey (diet), and predator tissue exceed guidelines with a high certainty of producing adverse effects under the most likely load scenario from a proposed San Luis Drain extension. High flows afford some protection under certain conditions in modeled San Joaquin River scenarios. However, meeting a combined goal of releasing a specific load during maximum flows and keeping selenium concentrations in the river below a certain objective to protect against bioaccumulation may not always be attainable. Management of the San Joaquin River on a constant concentration basis also could create problematic bioaccumulation during a wet year, especially during the low flow season, because high flows translate to high loads that are not always offset by seasonal river inflows. Prior to refinery cleanup, selenium contamination was sufficient to threaten reproduction in key species within the Bay-Delta ecosystems and human health advisories were posted based on selenium concentrations in tissues of diving ducks. During this time, selenium concentrations in the Bay-Delta were well below the most stringent recommended water quality criterion [1 microgram per liter (1 µg/L)]. Enhanced biogeochemical transformations to bioavailable particulate selenium and efficient bioaccumulation by bivalves characterized the system. If these biogeochemical conditions continue to prevail and agricultural selenium sources replace or exceed refinery sources, ecological forecasts suggest the risk of adverse effects will be difficult to eliminate under an out-of-valley resolution to the selenium problem. The Bay-Delta selenium model presented here is a systematic approach for conducting forecasts of the ecological effects from selenium on aquatic food webs. It is a new tool that links and models the major processes leading from loads through consumer organisms to predators. It also is a feasible approach for site-specific analysis and could provide a framework for developing new protective selenium foodweb guidelines and predator criteria. Model components that help ensure understanding ecosystems and the basis of environmental protection are (1) contaminant concentrations and speciation in sources, such as particulate material, that most influence bioavailability; (2) bioaccumulation models that calculate concentrations in diet, specifically in bivalves of the Bay-Delta that act as sensitive indicators of selenium contamination; (3) food-web type that determines what animals are threatened and when; and (4) multiple media concentrations (water, particulate material, and tissue of prey and predators) that, in-combination, determine risk or hazard.

California

A 42 year inference of cloud base height trends in the Luquillo Mountains of northeastern Puerto Rico

The Luquillo Mountains of eastern Puerto Rico are home to the only tropical rainforest managed by the United States Forest Service, with cloud-immersed forests historically occupying the highest elevations. However, within the past 50 yr, studies of the Luquillo cloud forest have suggested an increase in cloud base heights (CBH), although the CBH in the area was not quantified until recently. The present work uses radiosonde observations from nearby San Juan, Puerto Rico, to contextualize the present-day CBH within a 42 yr (1975−2016) proxy record and determine evidence for rising cloud base. Two key questions are addressed: (1) Can theoretical CBH calculations from San Juan provide a reasonable proxy for CBHs in the Luquillo Mountains? (2) Does a significant trend accompany the CBH lifting inferred from recent work in the region? The mean-layer lifted condensation level (MLLCL), a thermodynamic parameter expressing the altitude at which a rising air parcel reaches 100% relative humidity, serves as the proxy. The 42 yr MLLCL time series corroborates both the low CBHs claimed in the 1980s and the higher CBHs documented by recent work. When considering all available radiosonde data, statistically significant increasing CBH trends are detected for all seasons. However, when the record is standardized to correct for progressive vertical resolution improvements to radiosonde observations, recent CBH increases are more modest than initially indicated, and statistically significant increases are only apparent in the late rainfall season.

Puerto Rico

Pleistocene and Holocene landscape development of the South Platte River Corridor, Northeastern Colorado

This report provides a synthesis of geologic mapping and geochronologic research along the South Platte River between the town of Masters and the city of Fort Morgan, northeastern Colorado. This work was undertaken to better understand landscape development along this part of the river corridor. The focus is on times of rapid change within the fluvial system that had a marked effect on the landscape. The study area is susceptible to drought, which destabilizes vegetation and makes the landscape vulnerable to eolian activity. This is reflected in a landscape that is largely covered by eolian sand and lesser amounts of loess. Past glaciation of the river’s headwaters had a major influence on river discharge and sediment supply, as have major flood events particularly on unglaciated tributaries heading on the piedmont. In the mapping area, fluvial deposits of the South Platte River system span the Pliocene and early Pleistocene(?) deposits of Nussbaum Alluvium to present-day deposits of the active channel and floodplain. Results of the study indicate that along this stretch of the South Platte River, the early Pleistocene and first half of the middle Pleistocene were times of net incision, periodically interrupted by episodes of aggradation that resulted in deposition of alluvium that has been correlated to Rocky Flats Alluvium, Verdos Alluvium, and Slocum Alluvium. Net incision between depositional events formed a series of poorly preserved terrace deposits along the valley sides that are now largely covered by eolian deposits. Sometime after about 380 thousand years, the river cut a deep paleovalley into Upper Cretaceous Pierre Shale that was then filled with a thick sequence of inferred Louviers Alluvium (coeval with Bull Lake glaciation). Net aggradation continued during the late Pleistocene, resulting in burial of the Louviers paleovalley with a thick sequence of mainstream and sidestream Broadway Alluvium (coeval with Pinedale glaciation). Subsequent incision during the late Pleistocene–Holocene transition formed the Kersey (Broadway) terrace, whose riser forms a prominent bluff on the south side of the river valley. This episode of incision spanned a very short period and was followed by renewed aggradation that deposited the next-lower terrace alluvium (Kuner terrace alluvium). The Kuner terrace level was probably abandoned sometime around the beginning of the middle Holocene. Low terraces on the valley floor indicate that the river has been primarily cutting and backfilling laterally rather than incising during the late Holocene. Synthesis of geologic mapping and chronologic data generated in this study indicate that the South Platte River in northeastern Colorado likely was highly sensitive to rapidly changing environmental conditions or crossed threshold conditions that triggered rapid geomorphic response during major climate changes associated with the late Pleistocene–Holocene transition. Historical times have been another period marked by rapid incision, reflected by gully incision and headward erosion in tributary valleys draining the north side of the South Platte River. This historical erosion could be related at least in part to extensive construction of irrigation ditches and reservoirs in the late 1800s–early 1900s, which altered drainage paths and groundwater flow and could have amplified natural factors such as climate change or intrinsic geomorphic instabilities within the system.

Colorado

Research, monitoring, and evaluation of emerging issues and measures to recover the Snake River Fall Chinook Salmon ESU: January 2019 - December 2019

The portion of the Snake River fall Chinook salmon Oncorhynchus tshawytscha evolutionary significant unit (ESU) that spawns upstream of Lower Granite Dam transitioned from low to high abundance during 19922019 in association with U.S. Endangered Species Act recovery efforts and other federally mandated actions. This annual report focuses on (1) numeric and habitat use responses by natural- and hatchery-origin spawners, (2) phenotypic and numeric responses by natural-origin juveniles, and (3) use of a small unmanned aerial system (sUAS) to search for fall Chinook salmon redds and carcasses. Spawners have located and used most of the available spawning habitat and that habitat is gradually approaching redd capacity. Timing of spawning and fry emergence have been relatively stable, but effects of density dependence are evident in juvenile life stages. Apparent abundance of juvenile fall Chinook salmon has increased and we noted the following responses: parr dispersal from riverine rearing habitat into Lower Granite Reservoir has become earlier; growth rate (g/d) and dispersal size of parr declined; and passage timing of smolts from the two Snake River reaches has become earlier and downstream movement rate faster. These findings coupled with stock-recruitment analyses presented in this report provide evidence for density-dependence in the Snake River reaches and in Lower Granite Reservoir that was influenced by the expansion of the recovery program. The long-term goal is to use this information in a comprehensive modeling effort to conduct action-effectiveness and uncertainty research and to inform Fish Population, Hydrosystem, Harvest, Hatchery, and Predation and Invasive Species Management Research, Monitoring, and Evaluation (RM&E) progams. In 2019, the U.S. Geological Survey (USGS) shifted survey efforts in the Snake River toward deepwater redd searches and fish collection for parentage-based tagging (PBT) analyses because all unmanned aerial system (UAS) activities were suspended by the Department of the Interior two weeks into the spawning season. We counted 81 deepwater redds at 17 of the 29 sites surveyed. Redd depths averaged 3.6 m. We collected a total of 123 live fall Chinook salmon from 16 unique geographic locations that spanned 55 river kilometers. Forty-six fish were recovered at Eureka Bar (rkm 307.1) and Kirby Creek (rkm 352.0), which accounted for 37% of all collected fish in 2019. Most (73 fish) post-spawned salmon were collected from early to mid-November just after peak spawning. A summary of 2019 PBT results can be found in Appendix A.1. In 2019, we PIT tagged subyearling fall Chinook salmon in both the Snake and Clearwater rivers. In the Snake River, we tagged 410 fish with 8-mm tags, 666 fish with 9-mm tags, and 1,082 fish 12-mm tags. During seining, our recapture rate of previously tagged fish was slightly higher in the lower reach at 11.4% than in the upper reach at 10.6%. In an effort to represent more of the population through tagging, we tagged fish as small as 45 mm with 8-mm tags at one site in the upper reach. This allowed us to increase the number of fish tagged in that reach by 17.7%. An additional 10.6% of collected fish could have been tagged in the lower reach had we used 8-mm tags in that reach. In the Clearwater River, we tagged 2,451 subyearlings and recaptured 260 (10.6%) fish in the river and 66 fish (48 tagged by USGS, 18 tagged by the Nez Perce Tribe) at Lower Granite Dam during October to provide information for growth estimation. Within riverine habitats, growth in both length and mass were higher for fish tagged with 8-mm tags than with 9- and 12-mm tags. Estimated growth in length and mass of subyearlings was higher in Lower Granite Reservoir than in riverine habitats. We adapted existing statistical models used to estimate abundance of steelhead and spring/summer Chinook salmon for fall Chinook salmon passing Lower Granite Dam. Run reconstruction efforts to date at Lower Granite Dam for Snake River fall Chinook salmon have provided estimates of the number of returning adults but with no measure of uncertainty about the estimates. The objective of this study was to estimate the abundance, with uncertainty, of marked (coded-wire tagged CWT or adipose clipped) and unmarked fall Chinook salmon past Lower Granite Dam for return years 20032018. Estimating uncertainty is important for informing the state-space life cycle model (Chapter 5), which incorporates both observation and process uncertainty into parameter estimates. The coefficient of variation (CV) for log- abundance was 1.0% or less in all years, whereas the CV for abundance averaged 4.2% and ranged from 1.4% to 10.4% among years. Over the past five years, we have been developing a two-stage state-space life-cycle model for naturally produced fall Chinook salmon in the Snake River basin. Initial efforts focused on generating juvenile and adult abundance estimates, with estimates of uncertainty, for informing the life-cycle model. In this report we 1) describe the statistical life-cycle model and improvements made since our last report to the Independent Scientific Advisory Board (ISAB), 2) estimate the effects of covariates on key demographic parameters, and 3) use the fitted life- cycle model to simulate population trajectories under hydrosystem actions proposed for the NOAA 2020 Biological Opinion (hereafter, the Proposed Action). Major recent advancements to the model include revised juvenile abundance estimates, the ability to estimate smolt-to-adult return rates (SAR) separately for subyearling and yearling juvenile fall Chinook salmon, and improvements in the observation model for estimating age, sex, and outmigration structure in adult returns. We examined the effect of numerous environmental, hydrosystem, and ocean covariates on key demographic parameters but only a few covariates were significant. For the adult-to-juvenile transition, we found the maximum weekly river flows during the winter egg- incubation period had a significant negative effect on the juvenile outmigrant abundance from that brood year. For subyearling outmigrants, percent spill during the summer had a significant positive effect on SAR and the mean winter PDO (Pacific Decadal Oscillation) had a significant negative effect on SAR. For yearling outmigrants, NPGO (North Pacific Gyre Oscillation) had a significant positive effect on SAR. We used the fitted model to simulate population trajectories under the Proposed Action, and the median 10-year geometric mean abundance was 8,222 female spawners (interquartile range: 2,592 26,714). Overall, the probability of quasi- extinction (probability of falling below 50 female spawners for 4 consecutive years) was low, with only 1.6% of all simulations having a quasi-extinction probability >0.95. Although quasi-extinction probability was low, we did not assess the additional effect of climate change, which would be expected to increase quasi-extinction probability.

Idaho, Oregon, Washington

Anomalously high relief on Denali, Alaska, caused by tectonic, lithologic, and climatic drivers

We assess the growth of anomalously high relief on Denali, located in the Alaska Range, Alaska, and the tallest mountain in North America (6190 masl). Denali is 3000 m taller than most surrounding peaks. It lies inside a 19° restraining bend in the active Denali fault system that is moving at about 7 mm/yr, and thus there is a tectonic and structural driver for ongoing rock uplift. High relief around Denali is also due, in part, to its granitic rock type and low fracture density relative to adjacent metasediments. Here we show that unique climatic conditions at high elevations also contribute to the growth of relief. We examine 10Be concentrations in 1) three new gravel samples between 3500 and 5200 m elevation from sites unaffected by recent glacial erosion, 2) previously published samples from a sidewall of the Kahiltna Glacier from 2400 to 2800 masl, 3) previously published samples from medial moraines along the length of the Kahiltna Glacier from 1400 to ~500 masl, and 4) previously published alluvial samples along the Kahiltna River at an elevation of ~200 masl. These samples constitute a transect extending more than 5000 vertical meters, and the data establish that erosion rates decrease with elevation and contribute to the growth of relief. Erosion rates for the three new high-elevation samples are calculated to 4.6±0.6 mm/ka at 5.2 km, 28.6±3.7 mm/ka at 4.0 km, and 38±5 mm/ka at 3.5 km. Erosion rates at intermediate elevations, on the sidewall of the Kahiltna Glacier, range between 160-327 mm/ka. Along the medial moraines inferred erosion rates range between 140-537 mm/ka, and basin-wide erosion rates calculated from sediments in the river below the glacier range between 450-896 mm/ka. These differences in erosion rates can create relief of 3 km within 1-10 Ma, well within the estimated period of increase in rock uplift and exhumation on Denali over the last ~6 Ma. Meteorological data from 2130 to 5550 masl at 5 sites show it rarely gets above freezing above 4 km elevation, indicating that frost weathering currently plays a diminished role in erosion at high elevations. The immediate implication of this temperature and erosional correlation is an increase in relief. This is the first study to directly measure a significant decrease in erosion rates at high elevations in the relative absence of frost weathering. The results highlight the combined influence of rock type, glacial erosion, and permanent sub-zero temperatures on erosion rates. In combination with active faulting, the data explain the resultant increase in relief along the southern side of the Alaska Range over the past 100 ka.

Alaska

Drought resistance and resilience: The role of soil moisture–plant interactions and legacies in a dryland ecosystem

In many regions of the world, climate change is projected to reduce water availability through changes in the hydrological cycle, including more frequent and intense droughts, as well as seasonal shifts in precipitation. In water-limited ecosystems, such as drylands, lower soil water availability may exceed the adaptive capacity of many organisms, leading to cascading ecological effects during (concurrent effects) and after drought (legacy effects). The magnitude and duration of concurrent and legacy effects depends on drought intensity, duration and timing as well as the resistance and resilience of the ecosystem. Here, we investigated the effects of drought seasonality and plant community composition on two dominant perennial grasses, Achnatherum hymenoides (C 3 photosynthesis) and Pleuraphis jamesii (C 4 photosynthesis), in a dryland ecosystem. The experiment consisted of three precipitation treatments: control (ambient precipitation), cool-season drought (−66% ambient precipitation November–April) and warm-season drought (−66% ambient precipitation May–October), applied in two plant communities (perennial grasses with or without a large shrub, Ephedra viridis ) over a 3-year period. We examined the concurrent and legacy effects of seasonal drought on soil moisture, phenology and biomass. Drought treatments had strong concurrent and legacy effects on soil moisture, which impacted the phenology and biomass of the two grasses. Drought reduced growing season length by delaying green-up (cool-season drought) or advancing senescence (warm-season drought) and reduced biomass for both species. Biomass and phenology legacy effects from drought emerged in the second and third years of the experiment. While we observed differential sensitivity to drought legacies between the two grasses, we found limited evidence that shrub presence had interactive effects with the drought treatment. Synthesis . The results from this study highlight how abiotic and biotic legacies can develop and influence a community's resistance and resilience to subsequent droughts. When the frequency of repeated extreme events, such as recurring seasonal droughts, exceeds the capacity of organisms or ecosystems to recover (i.e. resilience), persistent drought legacies can reduce the resistance to subsequent drought events. Overall, these results highlight how drought legacies are a product of ecological resistance and resilience to past drought and can influence ecosystem vulnerability to future droughts.

Utah

Detection and characterization of coastal tidal wetland change in the northeastern US using Landsat time series

Coastal tidal wetlands are highly altered ecosystems exposed to substantial risk due to widespread and frequent land-use change coupled with sea-level rise, leading to disrupted hydrologic and ecologic functions and ultimately, significant reduction in climate resiliency. Knowing where and when the changes have occurred, and the nature of those changes, is important for coastal communities and natural resource management . Large-scale mapping of coastal tidal wetland changes is extremely difficult due to their inherent dynamic nature. To bridge this gap, we developed an automated algorithm for DEtection and Characterization of cOastal tiDal wEtlands change (DECODE) using dense Landsat time series. DECODE consists of three elements, including spectral break detection, land cover classification and change characterization. DECODE assembles all available Landsat observations and introduces a water level regressor for each pixel to flag the spectral breaks and estimate harmonic time-series models for the divided temporal segments. Each temporal segment is classified (e.g., vegetated wetlands, open water, and others – including unvegetated areas and uplands) based on the phenological characteristics and the synthetic surface reflectance values calculated from the harmonic model coefficients, as well as a generic rule-based classification system. This harmonic model-based approach has the advantage of not needing the acquisition of satellite images at optimal conditions (i.e., low tide status) to avoid underestimating coastal vegetation caused by the tidal fluctuation. At the same time, DECODE can also characterize different kinds of changes including land cover change and condition change (i.e., land cover modification without conversion). We used DECODE to track status of coastal tidal wetlands in the northeastern United States from 1986 to 2020. The overall accuracy of land cover classification and change detection is approximately 95.8% and 99.8%, respectively. The vegetated wetlands and open water were mapped with user's accuracy of 94.6% and 99.0%, and producer's accuracy of 98.1% and 93.5%, respectively. The cover change and condition change were mapped with user's accuracy of 68.0% and 80.0%, and producer's accuracy of 80.5% and 97.1%, respectively. Approximately 3283 km 2 of the coastal landscape within our study area in the northeastern United States changed at least once (12% of the study area), and condition changes were the dominant change type (84.3%). Vegetated coastal tidal wetland decreased consistently (~2.6 km 2 per year) in the past 35 years, largely due to conversion to open water in the context of sea-level rise.

Remote Sensing of Environment

Estimating the potential for submergence for two wetlands in the Mississippi River Delta

We used a combined field and modeling approach to estimate the potential for submergence for one rapidly deteriorating (Bayou Chitigue Marsh) and one apparently stable (Old Oyster Bayou Marsh) saltmarsh wetland in coastal Louisiana, given two eustatic sea level rise scenarios: the current rate (0.15 cm year −1 ); and the central value predicted by the Intergovernmental Panel on Climate Change (48 cm by the year 2100). We also used the model to determine what processes were most critical for maintaining and influencing salt marsh elevation including, mineral matter deposition, organic matter production, shallow subsidence (organic matter decomposition + primary sediment compaction), deep subsidence, and sediment pulsing events (e.g., hurricanes). Eight years of field measurements from feldspar marker horizons and surface elevation tables revealed that the rates of vertical accretion at the Bayou Chitigue Marsh were high (2.26 (0.09) cm yr −1 (mean ± SE)) because the marsh exists at the lower end of the tidal range. The rate of shallow subsidence was also high (2.04 (0.1) cm yr −1 ), resulting in little net elevation gain (0.22 (0.06) cm yr −1 ). In contrast, vertical accretion at the Old Oyster Bayou Marsh, which is 10 cm higher in elevation, was 0.48 (0.09) cm yr −1 . However, there was a net elevation gain of 0.36 (0.08) cm yr −1 because there was no significant shallow subsidence. When these rates of elevation gain were compared to rates of relative sea level rise (deep subsidence plus eustatic sea level rise), both sites showed a net elevation deficit although the Bayou Chitigue site was subsiding at approximately twice the rate of the Old Oyster Bayou site (1.1 cm yr −1 versus 0.49 cm yr −1 respectively). These field data were used to modify, initialize, and calibrate a previously published wetland soil development model that simulates primary production and mineral matter deposition as, feedback functions of elevation. Sensitivity analyses revealed that wetland elevation was most sensitive to changes in the rates of deep subsidence, a model forcing function that is difficult to measure in the field and for which estimates in the literature vary widely. The model also revealed that, given both the current rate of sea level rise and the central value estimate, surface elevation at both sites would fall below mean sea level over the next 100 years. Although these results were in agreement with the field study, they contradicted long term observations that the Old Oyster Bayou site has been in equilibrium with sea level for at least the past 50 years. Further simulations showed that the elevation at the Old Oyster Bayou site could keep pace with current rates of sea level rise if either a lower rate for deep subsidence was used as a forcing function, or if a periodic sediment pulsing function (e.g., from hurricanes) was programmed into the model.

Louisiana

Synoptic analysis of a decade of daily measurements of SO2 emission in the troposphere from volcanoes of the global ground-based Network for Observation of Volcanic and Atmospheric Change

Volcanic plumes are common and far-reaching manifestations of volcanic activity during and between eruptions. Observations of the rate of emission and composition of volcanic plumes are essential to recognize and, in some cases, predict the state of volcanic activity. Measurements of the size and location of the plumes are important to assess the impact of the emission from sporadic or localized events to persistent or widespread processes of climatic and environmental importance. These observations provide information on volatile budgets on Earth, chemical evolution of magmas, and atmospheric circulation and dynamics. Space-based observations during the last decades have given us a global view of Earth's volcanic emission, particularly of sulfur dioxide ( SO 2 ). Although none of the satellite missions were intended to be used for measurement of volcanic gas emission, specially adapted algorithms have produced time-averaged global emission budgets. These have confirmed that tropospheric plumes, produced from persistent degassing of weak sources, dominate the total emission of volcanic SO 2 . Although space-based observations have provided this global insight into some aspects of Earth's volcanism, it still has important limitations. The magnitude and short-term variability of lower-atmosphere emissions, historically less accessible from space, remain largely uncertain. Operational monitoring of volcanic plumes, at scales relevant for adequate surveillance, has been facilitated through the use of ground-based scanning differential optical absorption spectrometer (ScanDOAS) instruments since the beginning of this century, largely due to the coordinated effort of the Network for Observation of Volcanic and Atmospheric Change (NOVAC). In this study, we present a compilation of results of homogenized post-analysis of measurements of SO 2 flux and plume parameters obtained during the period March 2005 to January 2017 of 32 volcanoes in NOVAC. This inventory opens a window into the short-term emission patterns of a diverse set of volcanoes in terms of magma composition, geographical location, magnitude of emission, and style of eruptive activity. We find that passive volcanic degassing is by no means a stationary process in time and that large sub-daily variability is observed in the flux of volcanic gases, which has implications for emission budgets produced using short-term, sporadic observations. The use of a standard evaluation method allows for intercomparison between different volcanoes and between ground- and space-based measurements of the same volcanoes. The emission of several weakly degassing volcanoes, undetected by satellites, is presented for the first time. We also compare our results with those reported in the literature, providing ranges of variability in emission not accessible in the past. The open-access data repository introduced in this article will enable further exploitation of this unique dataset, with a focus on volcanological research, risk assessment, satellite-sensor validation, and improved quantification of the prevalent tropospheric component of global volcanic emission.

Earth System Science Data