USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Science”

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 757 records · Page 42Linked to original sources
USGS2018

A snow density dataset for improving surface boundary conditions in Greenland ice sheet firn modeling

The surface snow density of glaciers and ice sheets is of fundamental importance in converting volume to mass in both altimetry and surface mass balance studies, yet it is often poorly constrained. Site-specific surface snow densities are typically derived from empirical relations based on temperature and wind speed. These parameterizations commonly calculate the average density of the top meter of snow, thereby systematically overestimating snow density at the actual surface. Therefore, constraining surface snow density to the top 0.1 m can improve boundary conditions in high-resolution firn-evolution modeling. We have compiled an extensive dataset of 200 point measurements of surface snow density from firn cores and snow pits on the Greenland ice sheet. We find that surface snow density within 0.1 m of the surface has an average value of 315 kg m −3 with a standard deviation of 44 kg m −3 , and has an insignificant annual air temperature dependency. We demonstrate that two widely-used surface snow density parameterizations dependent on temperature systematically overestimate surface snow density over the Greenland ice sheet by 17–19%, and that using a constant density of 315 kg m −3 may give superior results when applied in surface mass budget modeling.

Frontiers in Earth Science↗
USGS2018

Analysis of the Alaska Volcano Observatory’s response time to volcanic explosions-1989 to 2016

A major goal of volcano monitoring is the rapid identification of volcanic explosions and subsequent warning of associated hazards. Between 1988 and 2016 the Alaska Volcano Observatory (AVO) responded to at least 54 separate volcanic eruptions. During this period, AVO's monitoring program relied principally on seismic and satellite remote sensing data, supplemented with geodetic, gas, and visual observations to track volcanic unrest. In this study we focus on AVO's response time, or the time required for AVO to (1) identify seismic signals associated with large ash-producing volcanic explosions and (2) initiate public warnings. We restrict this analysis to volcanoes monitored by a local seismic network and explosive in character. We focus on the 1989–90 eruption of Redoubt Volcano (VEI 3), the 1992 eruption of Mount Spurr (VEI 4), the 1999 eruption of Shishaldin Volcano (VEI 3), the 2006 eruption of Augustine Volcano (VEI 3) and the 2016 eruption of Pavlof Volcano (VEI 2) as detailed records of the timing of formal warnings are preserved. These eruption sequences allow us to evaluate AVO's response time under a number of monitoring scenarios, including both expected (those with recognized precursory unrest) and surprise eruptions (those without identified precursory unrest) as well as individual and repetitive sequences of explosive events. Recorded response time ranges from ~1 to 86 min. The shorter response times (~1–13 min) were achieved during sequences of explosive events at Redoubt (1989–90), Spurr (1992) and Augustine (2006). The longer response times (31– 86 min) are recorded for unexpected or surprise explosions such as Spurr (August 18, 1992) and Pavlof (2016) and the only or first explosions in an eruptive sequence such as Shishaldin (1999) and Augustine (2006).

Alaska↗
USGS2018

Alaska Volcano Observatory alert and forecasting timeliness: 1989–2017

The Alaska Volcano Observatory (AVO) monitors volcanoes in Alaska and issues notifications and warnings of volcanic unrest and eruption. We evaluate the timeliness and accuracy of eruption forecasts for 53 eruptions at 20 volcanoes, beginning with Mount Redoubt's 1989–1990 eruption. Successful forecasts are defined as those where AVO issued a formal warning before eruption onset. These warning notifications are now part of AVO's Aviation Color Code and Volcanic Alert Level. This analysis considers only the start of an eruption, although many eruptions have multiple phases of activity. For the 21 eruptions at volcanoes with functioning local seismic networks, AVO has high forecasting success at volcanoes with: >15 years repose intervals and magmatic eruptions (4 out of 4, 100%); or larger eruptions (Volcanic Explosivity Index (VEI) 3 or greater; 6 out of 10, 60%). Therefore, AVO successfully forecast all four monitored, longer-repose period, VEI 3+ eruptions: Redoubt 1989–1990 and 2009, Spurr 1992, and Augustine 2005–2006. For volcanoes with functioning seismic monitoring networks, success rates are lower for: volcanoes with shorter repose periods (3 out of 16, 19%); more mafic compositions (3 out of 18, 17%); or smaller eruption size (VEI 2 or less, 1 out of 11, 9%). These eruptions (Okmok, Pavlof, Veniaminof, and Shishaldin) often lack detectable precursory signals. For 32 eruptions at volcanoes without functioning local seismic networks, the forecasting success rate is much lower (2, 6%; Kasatochi 2008 and Shishaldin 2014). For remote volcanoes where the main hazard is to aviation, rapid detection is a goal in the absence of in situ monitoring. Eruption detection has improved in recent years, shown by a decrease in the time between eruption onset and notification. Even limited seismic monitoring can detect precursory activity at volcanoes with certain characteristics (intermediate composition, longer repose times, larger eruptions), but difficulty persists in detecting subtle precursory activity at frequently active volcanoes with more mafic compositions. This suggests that volcano-specific characteristics should be considered when designing monitoring programs and evaluating forecasting success. More proximally-located sensors and data types are likely needed to forecast eruptive activity at frequently-active, more mafic volcanoes that generally produce smaller eruptions.

Alaska↗
USGS2018

Prevalence of seismic rate anomalies preceding volcanic eruptions in Alaska

Seismic rate increases often precede eruptions at volcanoes worldwide. However, many eruptions occur without such precursors. Additionally, identifying seismic rate increases near volcanoes with high levels of background seismicity is non-trivial and many periods of elevated seismicity occur without ensuing eruptions, limiting their usefulness for forecasting in some cases. Although these issues are commonly known, efforts to quantify them are limited. In this study, we consistently apply a common statistical tool, the β-statistic, to seismically monitored eruptions in Alaska of various styles to determine the overall prevalence of seismic rate anomalies immediately preceding eruptions. We find that 6 out of 20 (30%) eruptions have statistically significant precursory seismic rate increases. Of these 6 eruptions, 3 of them occur at volcanoes with relatively felsic compositions, repose periods >15 years, and VEI ≥ 3. Overall, our results confirm that seismic rate increases are common prior to larger eruptions at long dormant, “closed-system” volcanoes, but uncommon preceding smaller eruptions at more frequently active, “open-system” volcanoes with more mafic magmas. We also explore the rate of other anomalies not precursory to eruptions and investigate their origins. Some of these non-eruptive anomalies can be explained by aftershocks of regional seismic events, magmatic activity that did not lead to eruption, or unrest at other nearby volcanoes. Some open-system volcanoes have high non-eruptive anomaly rates and low pre-eruptive anomaly rates and are thus not amenable to forecasting based on earthquake catalogs. In this study, we find that 31% of anomalies lead to eruption. With continued calibration at more volcanoes, the β-statistic that we apply may be used more broadly to analyze future periods of seismic unrest at other volcanoes, properly placing such episodes into the context of the long-term background rate. These results may be useful for informing future eruption forecasts around the world, and the statistical tool may aid volcano observatories in identifying future seismic rate anomalies under changing network conditions.

Alaska↗
USGS2018

Short-term forecasting and detection of explosions during the 2016–2017 eruption of Bogoslof volcano, Alaska

We describe a multidisciplinary approach to forecast, rapidly detect, and characterize explosive events during the 2016–2017 eruption of Bogoslof volcano, a back-arc shallow submarine volcano in Alaska’s Aleutian arc. The eruptive sequence began in December 2016 and included about 70 discrete explosive events. Because the volcano has no local monitoring stations, we used distant stations on the nearest volcanoes, Okmok (54 km) and Makushin (72 km), combined with regional infrasound sensors and lightning detection from the Worldwide Lightning Location Network (WWLLN). Pre-eruptive seismicity was detected for 12 events during the first half of the eruption; for all other events co-eruptive signals allowed for detection only. Monitoring of activity used a combination of scheduled checks combined with automated alarms. Alarms triggered on real-time data included real-time seismic amplitude measurement (RSAM); infrasound from several arrays, the closest being on Okmok; and lightning strokes detected from WWLLN within a 20-km radius of the volcano. During periods of unrest, a multidisciplinary response team of four people fulfilled specific roles to evaluate geophysical and remote-sensing data, run event-specific ash-cloud dispersion models, ensure interagency coordination, and develop and distribute of formalized warning products. Using this approach, for events that produced ash clouds ≥7.5 km above sea level, Alaska Volcano Observatory (AVO) called emergency response partners 15 min, and issued written notices 30 min, after event onset (mean times). Factors that affect timeliness of written warnings include event size and number of data streams available; bigger events and more data both decrease uncertainty and allow for faster warnings. In remote areas where airborne ash is the primary hazard, the approach used at Bogoslof is an effective strategy for hazard mitigation.

Alaska↗
USGS2018

Mount St. Helens retrospective: Lessons learned since 1980 and remaining challenges

Since awakening from a 123-year repose in 1980, Mount St. Helens has provided an opportunity to study changes in crustal magma storage at an active arc volcano—a process of fundamental importance to eruption forecasting and hazards mitigation. There has been considerable progress, but important questions remain unanswered. Was the 1980 eruption triggered by an injection of magma into an upper crustal reservoir? If so, when? How did magma rise into the edifice without producing detectable seismicity deeper than ∼2.5 km or measurable surface deformation beyond the volcano’s north flank? Would precursory activity have been recognized earlier if current monitoring techniques had been available? Despite substantial improvements in monitoring capability, similar questions remain after the dome-forming eruption of 2004–2008. Did additional magma accumulate in the reservoir between the end of the 1980–1986 eruption and the start of the 2004–2008 eruption? If so, when? What is the significance of a relative lull in seismicity and surface deformation for several years prior to the 2004–2008 eruption onset? How did magma reach the surface without producing seismicity deeper than ∼2 km or measurable deformation more than a few hundred meters from the vent? Has the reservoir been replenished since the eruption ended, and is it now primed for the next eruption? What additional precursors, if any, should be expected? This paper addresses these questions, explores possible answers, and identifies unresolved issues in need of additional study. The 1980–1986 and 2004–2008 eruptions could have resulted from second boiling during crystallization of magma long-resident in an upper crustal reservoir, rather than from injection of fresh magma from below. If reservoir pressurization and magma ascent were slow enough, resulting strain might have been accommodated by viscoelastic deformation, without appreciable seismicity or surface deformation, until rising magma entered a brittle regime within 2–2.5 km of the surface. Given the remarkably gas-poor nature of the 2004–2008 dome lava, future eruptive activity might require a relatively long period of quiescence and reservoir pressurization or a large injection of fresh magma—an event that arguably has not occurred since the Kalama eruptive period (C.E. 1479–1720).

Washington↗
USGS2018

Measuring SO2 emission rates at Kīlauea Volcano, Hawaii, using an array of upward-looking UV spectrometers, 2014-2017

Retrieving accurate volcanic sulfur dioxide (SO 2 ) gas emission rates is important for a variety of purposes. It is an indicator of shallow subsurface magma, and thus may signal impending eruption or unrest. SO 2 emission rates are significant for accurately assessing climate impact, and providing context for assessing environmental, agricultural, and human health effects during volcanic eruptions. The U.S. Geological Survey Hawaiian Volcano Observatory uses an array of ten fixed, upward-looking ultraviolet spectrometer systems to measure SO 2 emission rates at 10-s sample intervals from the Kīlauea summit. We present Kīlauea SO 2 emission rates from the volcano’s summit and middle East Rift Zone during 2014–2017 and discuss the major sources of error for these measurements. Due to the wide range of SO 2 emissions encountered at the summit vent, we used a variable wavelength spectral analysis range to accurately quantify both high and low SO 2 column densities. We compare measured emission rates from the fixed spectrometer array to independent road and helicopter-based traverse measurements and evaluate the magnitudes and sources of uncertainties for each method. To address the challenge of obtaining accurate plume speed measurements, we examine ground-based wind-speed, plume speed tracking via spectrometer, and SO 2 camera derived plume speeds. Our analysis shows that: (1) the summit array column densities calculated using a dual fit window, are within -6 to +22% of results obtained with a variety of other conventional and experimental retrieval methods; (2) emission rates calculated from the summit array located ∼3 km downwind provide the best, practical estimate of summit SO 2 release under normal trade wind conditions; (3) ground-based anemometer wind speeds are 22% less than plume speeds determined by cross-correlation of plume features; (4) our best estimate of average Kīlauea SO 2 release for 2014–2017 is 5100 t/d, which is comparable to the space-based OMI emissions of 5518 t/d; and (5) short-term variability of SO 2 emissions reflects Kīlauea lava lake dynamics.

Hawaii↗
USGS2018

Best practices for elevation-based assessments of sea-level rise and coastal flooding exposure

Elevation data are critical for assessments of sea-level rise (SLR) and coastal flooding exposure. Previous research has demonstrated that the quality of data used in elevation-based assessments must be well understood and applied to properly model potential impacts. The cumulative vertical uncertainty of the input elevation data substantially controls the minimum increments of SLR and the minimum planning horizons that can be effectively used in assessments. For regional, continental, or global assessments, several digital elevation models (DEMs) are available for the required topographic information to project potential impacts of increased coastal water levels, whether a simple inundation model is used or a more complex process-based or probabilistic model is employed. When properly characterized, the vertical accuracy of the DEM can be used to report assessment results with the uncertainty stated in terms of a specific confidence level or likelihood category. An accuracy evaluation has been conducted of global DEMs to quantify their inherent vertical uncertainty to demonstrate how accuracy information should be considered when planning and implementing a SLR or coastal flooding assessment. The evaluation approach includes comparison of the DEMs with high-accuracy geodetic control points as the independent reference data over a variety of coastal relief settings. The global DEMs evaluated include SRTM, ASTER GDEM, ALOS World 3D, TanDEM-X, NASADEM, and MERIT. High-resolution, high-accuracy DEM sources, such as airborne lidar and stereo imagery, are also included to give context to the results from the global DEMs. The accuracy characterization results show that current global DEMs are not adequate for high confidence mapping of exposure to fine increments (<1 m) of SLR or with shorter planning horizons (<100 years) and thus they should not be used for such mapping, but they are suitable for general delineation of low elevation coastal zones. In addition to the best practice of rigorous accounting for vertical uncertainty, other recommended procedures are presented for delineation of different types of impact areas (marine and groundwater inundation) and use of regional relative SLR scenarios. The requirement remains for a freely available, high-accuracy, high-resolution global elevation model that supports quantitative SLR and coastal inundation assessments at high confidence levels.

Frontiers in Earth Science↗
USGS2019

Isotopic and petrologic investigation, and a thermomechanical model of genesis of large-volume rhyolites in arc environments: Karymshina Volcanic Complex, Kamchatka, Russia

The Kamchatka Peninsula of eastern Russia is currently one of the most volcanically active areas on Earth where a combination of >8 cm/yr subduction convergence rate and thick continental crust generates large silicic magma chambers, reflected by abundant large calderas and caldera complexes. This study examines the largest center of silicic 4-0.5 Ma Karymshina Volcanic Complex, which includes the 25 × 15 km Karymshina caldera, the largest in Kamchatka. A series of rhyolitic tuff eruptions at 4 Ma were followed by the main eruption at 1.78 Ma and produced an estimated 800 km 3 of rhyolitic ignimbrites followed by high-silica rhyolitic post-caldera extrusions. The postcaldera domes trace the 1.78 Ma right fracture and form a continuous compositional series with ignimbrites. We here present results of a geologic, petrologic, and isotopic study of the Karymshina eruptive complex, and present new Ar-Ar ages, and isotopic values of rocks for the oldest pre- 1.78 Ma caldera ignimbrites and intrusions, which include a diversity of compositions from basalts to rhyolites. Temporal trends in δ 18 O, 87 Sr/ 86 Sr, and 144 Nd/ 143 Nd indicate values comparable to neighboring volcanoes, increase in homogeneity, and temporal increase in mantle-derived Sr and Nd with increasing differentiation over the last 4 million years. Data are consistent with a batholithic scale magma chamber formed by primarily fractional crystallization of mantle derived composition and assimilation of Cretaceous and younger crust, driven by basaltic volcanism and mantle delaminations. All rocks have 35–45% quartz, plagioclase, biotite, and amphibole phenocrysts. Rhyolite-MELTS crystallization models favor shallow (2 kbar) differentiation conditions and varying quantities of assimilated amphibolite partial melt and hydrothermally-altered silicic rock. Thermomechanical modeling with a typical 0.001 km 3 /yr eruption rate of hydrous basalt into a 38 km Kamchatkan arc crust produces two magma bodies, one near the Moho and the other engulfing the entire section of upper crust. Rising basalts are trapped in the lower portion of an upper crustal magma body, which exists in a partially molten to solid state. Differentiation products of basalt periodically mix with the resident magma diluting its crustal isotopic signatures. At the end of the magmatism crust is thickened by 8 km. Thermomechanical modeling show that the most likely way to generate large spikes of rhyolitic magmatism is through delamination of cumulates and mantle lithosphere after many millions of years of crustal thickening. The paper also presents a chemical dataset for Pacific ashes from ODDP 882 and 883 and compares them to Karymshina ignimbrites and two other Pleistocene calderas studied by us in earlier works.

Kamchatka↗
USGS2019

Insights into the mechanisms of phreatic eruptions from continuous high frequency volcanic gas monitoring: Rincón de la Vieja volcano, Costa Rica

Understanding the trigger mechanisms of phreatic eruptions is key to mitigating the effects of these hazardous but poorly forecastable volcanic events. It has recently been established that high-rate volcanic gas observations are potentially very suitable to identifying the source processes driving phreatic eruptions, and to eventually detecting precursory changes prior to individual phreatic blasts. In February-May 2017, we deployed a Multi-GAS instrument to continuously monitor gas concentrations in the crater lake plume of Rincón de la Vieja, a remote and poorly monitored active volcano in Costa Rica, site of frequent phreatic/phreatomagmatic eruptions. Forty-two phreatic/phreatomagmatic eruptions were seismically recorded during our investigated period, 9 of which were also recorded for gas by the Multi-GAS. To the best of our knowledge, these represent the first instrumentally measured gas compositions during individual phreatic/phreatomagmatic explosions at an active volcano. Our results show that during background quiescent degassing the Rincón de la Vieja crater lake plume was characterized by high CO 2 /SO 2 ratios of 64 ± 59 and H 2 S/SO 2 ratios of 0.57 ± 0.20. This composition is interpreted as reflecting hydrothermal (re)processing of magma-sourced gas in the sub-limnic environment. Phreatic blasts were recorded by the Multi-GAS as brief (1–2 min long) pulses of elevated gas mixing ratios (up to ~52 ppmv SO 2 and >3,000 ppmv CO 2 ), or more than an order of magnitude higher than during background degassing (~1 ppmv SO 2 and ~450 ppmv CO 2 ). During the phreatic eruption(s), the H 2 S/SO 2 ratio was systematically lower (<0.18) than during background degassing, but the CO 2 /SO 2 ratio remained high (and variable), ranging from 37 to 390. These S-poor compositions for the eruptive gas imply extensive processing of the source magmatic gas during pre-eruptive hydrothermal storage, likely by deposition of native S and/or sulfate. Our gas results are thus overall consistent with a mechanism of phreatic eruptions triggered by accumulation of magmatic-hydrothermal gases beneath a hydrothermal seal. We claim that real-time Multi-GAS monitoring is urgently needed at other crater lake-hosting volcanoes (e.g., Ruapehu, Aso), where phreatic eruptions may similarly be preceded by phases of reduced S degassing at the surface.

Rincón de la Vieja Volcano↗
USGS2019

Space-based imaging radar studies of U.S. volcanoes

The arrival of space-based imaging radar as a revolutionary land-surface mapping and monitoring tool little more than a quarter century ago enabled a spate of innovative volcano research worldwide. Soon after launch of European Space Agency’s ERS-1 spacecraft in 1991, the U.S. Geological Survey began SAR and InSAR studies of volcanoes in the Aleutian and Cascades arcs, in Hawai’i, and elsewhere in the western U.S. including the Yellowstone and Long Valley calderas. This paper summarizes results of that effort and presents new findings concerning: (1) prevalence of volcano deformation in the Aleutian and Cascade arcs; (2) surface-change detection and hazard assessment during eruptions at Aleutian and Hawaiian volcanoes; (3) geodetic imaging of magma storage and transport systems in Hawai’i; and (4) deformation sources and processes at the Yellowstone and Long Valley calderas. Surface deformation caused by a variety of processes is common in arc settings and could easily escape detection without systematic InSAR surveillance. Space-based SAR imaging of active lava flows and domes in remote or heavily vegetated settings, including during periods of bad weather and darkness, extends land-based monitoring capabilities and improves hazards assessments. At Kīlauea Volcano, comprehensive SAR and InSAR observations identify multiple magma storage zones beneath the summit area and along the East Rift Zone, and illuminate magma transport pathways. The same approach at Yellowstone tracks the ascent of magmatic volatiles from a mid-crustal intrusion to shallow depth and relates that process to increased hydrothermal activity at the surface. Together with recent and planned launches of highly capable imaging-radar satellites, these findings support an optimistic outlook for near-real time surveillance of volcanoes at global scale in the coming decade.

Frontiers in Earth Science↗
USGS2019

Size distributions of Arctic waterbodies reveal consistent relations in their statistical moments in space and time

Arctic lowlands are characterized by large numbers of small waterbodies, which are known to affect surface energy budgets and the global carbon cycle. Statistical analysis of their size distributions has been hindered by the shortage of observations at sufficiently high spatial resolutions. This situation has now changed with the high-resolution (<5 m) circum-Arctic Permafrost Region Pond and Lake (PeRL) database recently becoming available. We have used this database to make the first consistent, high-resolution estimation of Arctic waterbody size distributions, with surface areas ranging from 0.0001 km 2 (100 m 2 ) to 1 km 2 . We found that the size distributions varied greatly across the thirty study regions investigated and that there was no single universal size distribution function (including power-law distribution functions) appropriate across all of the study regions. We did, however, find close relationships between the statistical moments (mean, variance, and skewness) of the waterbody size distributions from different study regions. Specifically, we found that the spatial variance increased linearly with mean waterbody size ( R 2 = 0.97, p < 2.2e-16) and that the skewness decreased approximately hyperbolically. We have demonstrated that these relationships (1) hold across the 30 Arctic study regions covering a variety of (bio)climatic and permafrost zones, (2) hold over time in two of these study regions for which multi-decadal satellite imagery is available, and (3) can be reproduced by simulating rising water levels in a high-resolution digital elevation model. The consistent spatial and temporal relationships between the statistical moments of the waterbody size distributions underscore the dominance of topographic controls in lowland permafrost areas. These results provide motivation for further analyses of the factors involved in waterbody development and spatial distribution and for investigations into the possibility of using statistical moments to predict future hydrologic dynamics in the Arctic.

Alaska↗
USGS2019

An assessment of plant species differences on cellulose oxygen isotopes from two Kenai Peninsula, Alaska peatlands: Implications for hydroclimatic reconstructions

Peat cores are valuable archives of past environmental change because they accumulate plant organic matter over millennia. While studies have primarily focused on physical, ecological, and some biogeochemical proxies, cores from peatlands have increasingly been used to interpret hydroclimatic change using stable isotope analyses of cellulose preserved in plant remains. Previous studies indicate that the stable oxygen isotope compositions (δ 18 O) preserved in alpha cellulose extracted from specific plant macrofossils reflect the δ 18 O values of past peatland water and thereby provide information on long-term changes in hydrology in response to climate. Oxygen isotope analyses of peat cellulose (δ 18 O cellulose ) have been successfully developed from peat cores that accumulate the same species for millennia. However, to fully exploit the potential of this proxy in species-diverse fens, studies are needed that account for the isotopic variations caused by changes in dominant species composition. This study assesses variation in δ 18 O values among peatland plant species and how they relate to environmental waters in two fens informally named Horse Trail and Goldfin, located on the leeward (dry) and windward (wet) side, respectively, of the climatic gradient across the Kenai Peninsula, Alaska. Environmental water δ 18 O values at both fens reflect unmodified δ 18 O values of mean annual precipitation, although at Goldfin standing pools were slightly influenced by evaporation. Modern plant [mosses and Carex spp. (sedges)] δ 18 O cellulose values indicate that all Carex spp. are higher (~2.5‰) than those of mosses, likely driven by their vascular structure and ecophysiological difference from non-vascular mosses. Moss δ 18 O cellulose values within each peatland are similar among the species, and differences appear related to evaporation effects on environmental waters within hummocks and hollows. The plant taxa-environmental water δ 18 O differences are applied to the previously determined Horse Trail Fen untreated bulk δ 18 O record. Results include significant changes to inferred millennial-to-centennial scale hydroclimatic trends where dominant taxa shift from moss to Carex spp., indicating that modern calibration datasets are necessary for interpreting stable isotopes from fens, containing a mix of vascular and nonvascular plants. Accounting for isotopic offsets through macrofossil analysis and modern plant-water isotope measurements opens new opportunities for hydroclimatic reconstructions from fen peatlands.

Alaska↗
USGS2019

Late-Glacial paleoecology of the Middle Susitna Valley, Alaska: Environmental context for human dispersal

We present here the results of multi-proxy analyses (sediment geochemistry, diatoms, and pollen) from sediment cores collected at four lakes in the middle Susitna Valley, Alaska. These lakes form a transect from the tundra to the boreal forest. The retrieved cores span from ~12,000 cal yr BP to the present, with age control provided by radiometric dates and tephra deposits, some of which are newly identified. Results indicate that deglaciation occurred before 12,000 cal yr BP and that by that time, the lakes were deep, productive, and surrounded by shrub tundra. The lake with the highest sampling resolution indicates a brief climatic reversal ~11,500 cal yr BP with decreased diatom-inferred lake level and lowered lake productivity, and reduced shrub presence. During the early to middle Holocene, all of the sedimentary records provide evidence of climatic amelioration with tree expansion and productive lakes. A middle to late Holocene climatic deterioration with reduced trees and a shallower, less productive lake is also indicated. In addition, the prominent Watana tephra at ~4,000 cal yr BP likely reduced lake productivity and affected the vegetation. Even though the region was relatively productive soon after deglaciation, people did not occupy the region until ~11,000 cal yr BP, about 1000 years later, and then only sparsely. By the middle and late Holocene, the region was more densely populated and this shift in human occupancy presumably reflects changes in resource abundance, especially caribou. Whether the Watana ashfall influenced caribou abundance and thus people, is still under investigation, but given the tephra’s effect on vegetation and lake productivity, it seems likely.

Alaska↗
USGS2019

Holocene thermokarst lake dynamics in northern Interior Alaska: The interplay of climate, fire, and subsurface hydrology

The current state of permafrost in Alaska and meaningful expectations for its future evolution are informed by long-term perspectives of previous permafrost degradation. Thermokarst processes in permafrost landscapes often lead to widespread lake formation and the spatial and temporal evolution of thermokarst lake landscapes reflects the combined effects of climate, ground conditions, vegetation, and fire. This study provides detailed analyses of thermokarst lake sediments of Holocene age from the southern loess uplands of the Yukon Flats; including bathymetry and sediment core analyses across a water depth transect. The sediment core results, dated by radiocarbon and 210Pb, indicate the onset of finely laminated lacustrine sedimentation between ~10,000 and 9,000 cal yr BP following basin development through inferred thermokarst processes. Thermokarst expansion to modern shoreline configurations continued until ~5000 cal yr BP, which may have been influenced by increased fire. Between ~5000 and 2000 cal yr BP, the preservation of fine laminations at intermediate and deep-water depths indicate higher lake levels than present. At that time, the lake likely overflowed into an over-deepened gully system that is no longer occupied by perennial streams. By ~2000 cal yr BP, massive sedimentation at intermediate water depths indicates that lake levels lowered, which is interpreted to reflect a response to drier conditions based on correspondence with Yukon Flats regional fire and local paleoclimate reconstructions. Consideration of additional contributing mechanisms include the possible influence of catastrophic lake drainages on downgradient base flow levels that may have enhanced subsurface water loss, although this mechanism is untested. The overall consistency between the millennial lake level trends documented here with regional paleoclimate trends indicates that after lakes formed, their size and depth has likely been affected directly by North Pacific atmospheric circulation changes and indirectly through evolution of permafrost, ground ice and sub-surface hydrology. As the first detailed study of Holocene thermokarst basin expansion, stabilization and subsequent climate-driven lake level variations in a loess upland, results provide a framework for future investigations of paleoclimatic signals from similar lake systems that characterize large regions of Alaska and Siberia.

Alaska↗
USGS2019

A stratigraphic approach to inferring depositional ages from detrital geochronology data

With the increasing use of detrital geochronology data for provenance analyses, we have also developed new constraints on the age of otherwise undateable sedimentary deposits. Because a deposit can be no older than its youngest mineral constituent, the youngest defensible detrital mineral age defines the maximum depositional age of the sampled bed. Defining the youngest `defensible' age in the face of uncertainty (e.g., analytical and geological uncertainty, or sample contamination) is challenging. The current standard practice of finding multiple detrital minerals with indistinguishable ages provides confidence that a given age is not an artifact; however, we show how requiring this overlap reduces the probability of identifying the true youngest component age. Barring unusual complications, the principle of superposition dictates that sedimentary deposits must get younger upsection. This fundamental constraint can be incoporated into the analysis of depositional ages in sedimentary sections through the use of Bayesian statistics, allowing for the inference of bounded estimates of true depositional ages and uncertainties from detrital geochronology so long as some minimum age constraints are present. We present two approaches for constructing a Bayesian model of deposit ages, first solving directly for the ages of deposits with the prior constraint that the ages of units must obey stratigraphic ordering, and second describing the evolution of ages with a curve that represents the sediment accumulation rate. Using synthetic examples we highlight how this method preforms in less-than-ideal circumstances. In an example from the Magallanes Basin of Patagonia, we demonstrate how introducing other age information from the stratigraphic section (e.g., fossil assemblages or radiometric dates) and formalizing the stratigraphic context of samples provides additional constraints on and information regarding depositional ages or derived quantities (e.g., sediment accumulation rates) compared to isolated analysis of individual samples.

Frontiers in Earth Science↗
USGS2019

Stratigraphic occurrences of sub-polar planktic foraminifera in pleistocene sediments on the Lomonosov Ridge, Arctic Ocean

Turborotalita quinqueloba is a species of planktic foraminifera commonly found in the sub-polar North Atlantic along the pathway of Atlantic waters in the Nordic seas and sometimes even in the Arctic Ocean, although its occurrence there remains poorly understood. Existing data show that T. quinqueloba is scarce in Holocene sediments from the central Arctic but abundance levels increase in sediments from the last interglacial period [Marine isotope stage (MIS) 5, 71–120 ka] in cores off the northern coast of Greenland and the southern Mendeleev Ridge. Turborotalita also occurs in earlier Pleistocene interglacials in these regions, with a unique and widespread occurrence of the less known Turborotalita egelida morphotype, proposed as a biostratigraphic marker for MIS 11 (474–374 ka). Here we present results from six new sediment cores, extending from the central to western Lomonosov Ridge, that show a consistent Pleistocene stratigraphy over 575 km. Preliminary semi-quantitative assessments of planktic foraminifer abundance and assemblage composition in two of these records (LOMROG12-7PC and AO16-5PC) reveal two distinct stratigraphic horizons containing Turborotalita in MIS 5. Earlier occurrences in Pleistocene interglacials are recognized, but contain significantly fewer specimens and do not appear to be stratigraphically coeval in the studied sequences. In all instances, the Turborotalita specimens resemble the typical T. quinqueloba morphotype but are smaller (63–125 μm), smooth-walled and lack the final thickened calcite layer common to adults of the species. These results extend the geographical range for T. quinqueloba in MIS 5 sediments of the Arctic Ocean and provide compelling evidence for recurrent invasions during Pleistocene interglacials.

Frontiers in Earth Science↗
USGS2019

Late Quaternary vegetation development following deglaciation of northwestern Alexander Archipelago, Alaska

The Cordilleran Ice sheet covered most of southeastern Alaska during the Last Glacial Interval (LGI: Marine Isotope Stage 2). Ice began to recede from western Alexander Archipelago ~17,000 + 700 yr BP. In this study pollen analysis and radiocarbon dating of three sediment cores were used to reconstruct, for the first time, the postglacial development of vegetation of the northwestern Alexander Archipelago during the past ~15,240 cal yr. Hummingbird Lake (HL), on southwestern Baranof Island, yielded a sediment core with one of the longest dated records from southeastern Alaska. The earliest part of the HL pollen record (~15,240-14,040 yr BP) indicates that the earliest vegetation was pine (Pinus contorta subsp. contorta) parkland with willows (Salix), heaths (Ericaceae), sedges (Cyperaceae), grasses (Poaceae), herbs and ferns. Starting at ~14,040 yr BP, alder (Alnus) rapidly colonized the area as pine populations declined. By 11,400 yr BP, Sitka spruce (Picea sitchensis) colonized the area, and soon became the dominant conifer. Mountain hemlock (Tsuga mertensiana) also colonized the area by ~11,400 yr BP, followed by western hemlock (Tsuga heterophylla) at ~10,200 yr BP. By ~9200 yr BP, western hemlock had become the dominant species in the area. During the late Holocene yellow cedar (Chamaecyparis nootkatensis) became established. Two marine sediment cores were also analyzed for pollen, with the oldest core from Lower Sitka Sound, between Kruzof and Baranof Islands. The lower part of the core consists of interlayered tephras and freshwater lake muds that are estimated to be ~13,150 to 14,000 yr BP. Pollen evidence indicates that the early postglacial vegetation around Sitka Sound was pine parkland with alders and abundant ferns. Damage to vegetation around Sitka Sound by volcanic eruptions is suggested by abrupt, large shifts in alder and pine pollen, and fern spores in samples adjacent to tephra layers. A marine sediment core from Slocum Arm, a fjord on the western coast of Chichagof Island, has a basal age of ~10,000 yr BP. The pollen record is similar to the Holocene pollen record at Hummingbird Lake. The sequence of vegetation changes interpreted from the three northwestern Alexander Archipelago pollen records are similar to those from other well-dated sites in southeastern Alaska, although chronologies differ between sites.

Alaska↗