skip to main content

This content will become publicly available on January 1, 2023

Title: A statistical framework for integrating nonparametric proxy distributions into geological reconstructions of relative sea level
Abstract. Robust, proxy-based reconstructions of relative sea-level (RSL) change are critical to distinguishing the processes that drive spatial and temporal sea-level variability. The relationships between individual proxies and RSL can be complex and are often poorly represented by traditional methods that assume Gaussian likelihood distributions. We develop a new statistical framework to estimate past RSL change based on nonparametric, empirical modern distributions of proxies in relation to RSL, applying the framework to corals and mangroves as an illustrative example. We validate our model by comparing its skill in reconstructing RSL and rates of change to two previous RSL models using synthetic time-series datasets based on Holocene sea-level data from South Florida. The new framework results in lower bias, better model fit, and greater accuracy and precision than the two previous RSL models. We also perform sensitivity tests using sea-level scenarios based on two periods of interest – meltwater pulses (MWPs) and the Holocene – to analyze the sensitivity of the statistical reconstructions to the quantity and precision of proxy data; we define high-precision indicators, such as mangroves and the reef-crest coral Acropora palmata, with 2σ vertical uncertainties within ± 3 m and lower-precision indicators, such as Orbicella spp., with 2σ vertical uncertainties more » within ± 10 m. For reconstructing rapid rates of change in RSL of up to ∼ 40 m kyr−1, such as those that may have characterized MWPs during deglacial periods, we find that employing the nonparametric model with 5 to 10 high-precision data points per kiloyear enables us to constrain rates to within ± 3 m kyr−1 (1σ). For reconstructing RSL with rates of up to ∼ 15 m kyr−1, as observed during the Holocene, we conclude that employing the model with 5 to 10 high-precision (or a combination of high- and low-precision) data points per kiloyear enables precise estimates of RSL within ±∼ 2 m (2σ) and accurate RSL reconstructions with errors ≲ 0.7 m. Employing the nonparametric model with only lower-precision indicators also produces fairly accurate estimates of RSL with errors ≲1.50 m, although with less precision, only constraining RSL to ±∼ 3–4 m (2σ). Although the model performs better than previous models in terms of bias, model fit, accuracy, and precision, it is computationally expensive to run because it requires inverting large matrices for every sample. The new model also provides minimal gains over similar models when a large quantity of high-precision data are available. Therefore, we recommend incorporating the nonparametric likelihood distributions when no other information (e.g., reef facies or epibionts indicative of shallow-water environments to refine coral elevational uncertainties) or no high-precision data are available at a location or during a given time period of interest. « less
; ; ; ;
Award ID(s):
2002437 2041325 1702740
Publication Date:
Journal Name:
Advances in Statistical Climatology, Meteorology and Oceanography
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Recent ice-mass loss driven by warming along the Antarctic Peninsula has resulted in rapid changes in uplift rates across the region. Are such events only a function of recent warming? If not, does the Earth response to such events last long enough to be preserved in Holocene records of relative sea level (RSL), and thus have a bearing on global-scale glacial isostatic adjustment (GIA) models (e.g. ICE-6G)? Answering such questions in Antarctica is hindered by the scarcity of RSL reconstructions within the region. Here, a new RSL reconstruction for Antarctica is presented based on beach ridges from Joinville Islandmore »on the Antarctic Peninsula. We find that RSL has fallen 4.9 ± 0.58 m over the past 3100 yr, and that the island experienced a significant increase in the rate of RSL fall from 1540 ± 125 cal. (calibrated) yr B.P. to 1320 ± 125 cal. yr B.P. This increase in the rate of RSL fall is likely due to the viscoelastic response of the solid Earth to terrestrial ice-mass loss from the Antarctic Peninsula, similar to the Earth response experienced after ice-mass loss following acceleration of glaciers behind the collapsed Larsen B ice shelf in 2002 C.E. Additionally, slower rates of beach-ridge progradation from 695 ± 190 cal. yr B.P. to 235 ± 175 cal. yr B.P. potentially reflect erosion of beach ridges from a RSL rise induced by a local glacial advance. The rapid response of the Earth to minor ice-mass changes recorded in the RSL record further supports recent assertions of a more responsive Earth to glacial unloading and at time scales relevant for GIA of Holocene and Pleistocene sea levels. Thus, current continental and global GIA models may not accurately capture the ice-mass changes of the Antarctic ice sheets at decadal and centennial time scales.« less
  2. ABSTRACT Salt-marsh foraminifera are sea-level proxies used to quantitatively reconstruct Holocene paleo-marsh elevations (PME) and subsequently relative sea level (RSL). The reliability of these reconstructions is partly dependent upon counting enough foraminifera to accurately characterize assemblages, while counting fewer tests allows more samples to be processed. We test the influence of count size on PME reconstructions by repeatedly subsampling foraminiferal assemblages preserved in a core of salt-marsh peat (from Newfoundland, Canada) with unusually large counts (up to 1595). Application of a single, weighted-averaging transfer function developed from a regional-scale modern training set to these ecologically-plausible simulated assemblages generated PME reconstructionsmore »at count sizes of 10–700. Reconstructed PMEs stabilize at counts sizes greater than ∼50 and counts exceeding ∼250 tests show little return for the additional time invested. The absence of some rare taxa in low counts is unlikely to markedly influence results from weighted-averaging transfer functions. Subsampling of modern foraminifera indicates that cross-validated transfer function performance shows only modest improvement when more than ∼40 foraminifera are counted. Studies seeking to understand multi-meter and millennial scale RSL trends should count more than ∼50 tests. The precision sought by studies aiming to resolve decimeter- and decadal-scale RSL variability is best achieved with counts greater than ∼75. In most studies seeking to reconstruct PME, effort is more productively allocated by counting relatively fewer foraminifera in more core samples than in counting large numbers of individuals. Target count sizes of 100–300 in existing studies are likely conservative and robust. Given the low diversity of salt-marsh foraminiferal assemblages, our results are likely applicable throughout and beyond northeastern North America.« less
  3. The West Antarctic Ice Sheet (WAIS) is largely marine based and thus highly sensitive to both climatic and oceanographic changes. Therefore, the WAIS has likely had a very dynamic history over the last several million years. A complete collapse of the WAIS would result in a global sea level rise of 3.3–4.3 m, yet the world’s scientific community is not able to predict its future behavior. Moreover, knowledge about past behavior of the WAIS is poor, in particular during geological times with climatic conditions similar to those expected for the near and distant future. Reconstructions and quantifications of partial ormore »complete WAIS collapses in the past are urgently needed for constraining and testing ice sheet models that aim to predict future WAIS behavior and the potential contribution of the WAIS to global sea level rise. Large uncertainties exist regarding the chronology, extent, rates, and spatial and temporal variability of past advances and retreats of the WAIS across the continental shelves. These uncertainties largely result from the fundamental lack of data from drill cores recovered proximal to the WAIS. The continental shelf and rise of the Amundsen Sea are prime targets for drilling because the records are expected to yield archives of pure WAIS dynamics unaffected by other ice sheets and the WAIS sector draining into the Amundsen Sea Embayment (ASE) currently experiences the largest ice loss in Antarctica (Paolo et al., 2015). We propose a series of drill sites for the ASE shelf where seismic data reveal seaward-dipping sedimentary sequences that span from the preglacial depositional phase to the most recent glacial periods. Our strategy is to drill a transect from the oldest sequences close to the bedrock/basin boundary at the middle–inner shelf transition to the youngest sequences on the outer shelf in the eastern ASE. If the eastern ASE is inaccessible due to sea ice cover, a similar transect of sites can be drilled on the western ASE. The core transect will provide a detailed history of the glacial cycles in the Amundsen Sea region and allow comparison to the glacial history from the Ross Sea sector. In addition, deep-water sites on the continental rise of the Amundsen Sea are selected for recovering continuous records of glacially transported sediments and detailed archives of climatic and oceanographic changes throughout glacial–interglacial cycles. We will apply a broad suite of analytical techniques, including multiproxy analyses, to address our objectives of reconstructing the onset of glaciation in the greenhouse to icehouse transition, processes of dynamic ice sheet behavior during the Neogene and Quaternary, and ocean conditions associated with the glacial cycles. The five principal objectives of Expedition 379 are as follows: 1. To reconstruct the glacial history of West Antarctica from the Paleogene to recent times and the dynamic behavior of the WAIS during the Neogene and Quaternary, especially possible partial or full WAIS collapses, and the WAIS contribution to past sea level changes. Emphasis is placed in particular on studying the response of the WAIS at times when the pCO2 in Earth’s atmosphere exceeded 400 ppm and atmospheric and oceanic temperatures were higher than at present. 2. To correlate the WAIS-proximal records of ice sheet dynamics in the Amundsen Sea with global records of ice volume changes and proxy records for air and seawater temperatures. 3. To study the relationship between incursions of warm Circumpolar Deep Water (CDW) onto the continental shelf of the Amundsen Sea Embayment and the stability of marine-based ice sheet margins under warm water conditions. 4. To reconstruct the processes of major WAIS advances onto the middle and outer shelf that are likely to have occurred since the middle Miocene and compare their timing and processes to those of other Antarctic continental shelves. 5. To identify the timing of the first ice sheet expansion onto the continental shelf of the ASE and its possible relationship to the uplift of Marie Byrd Land.« less
  4. Accelerating erosion of the Alaska Beaufort Sea coast is increasing inputs of organic matter from land to the Arctic Ocean, and improved estimates of organic matter stocks in eroding coastal permafrost are needed to assess their mobilization rates under contemporary conditions. We collected three permafrost cores (4.5–7.5 m long) along a geomorphic gradient near Drew Point, Alaska, where recent erosion rates average 17.2 m year −1 . Down-core patterns indicate that organic-rich soils and lacustrine sediments (12–45% total organic carbon; TOC) in the active layer and upper permafrost accumulated during the Holocene. Deeper permafrost (below 3 m elevation) mainly consists of Late Pleistocene marinemore »sediments with lower organic matter content (∼1% TOC), lower C:N ratios, and higher δ 13 C values. Radiocarbon-based estimates of organic carbon accumulation rates were 11.3 ± 3.6 g TOC m −2  year −1 during the Holocene and 0.5 ± 0.1 g TOC m −2  year −1 during the Late Pleistocene (12–38 kyr BP). Within relict marine sediments, porewater salinities increased with depth. Elevated salinity near sea level (∼20–37 in thawed samples) inhibited freezing despite year-round temperatures below 0°C. We used organic matter stock estimates from the cores in combination with remote sensing time-series data to estimate carbon fluxes for a 9 km stretch of coastline near Drew Point. Erosional fluxes of TOC averaged 1,369 kg C m −1  year −1 during the 21st century (2002–2018), nearly doubling the average flux of the previous half-century (1955–2002). Our estimate of the 21st century erosional TOC flux year −1 from this 9 km coastline (12,318 metric tons C year −1 ) is similar to the annual TOC flux from the Kuparuk River, which drains a 8,107 km 2 area east of Drew Point and ranks as the third largest river on the North Slope of Alaska. Total nitrogen fluxes via coastal erosion at Drew Point were also quantified, and were similar to those from the Kuparuk River. This study emphasizes that coastal erosion represents a significant pathway for carbon and nitrogen trapped in permafrost to enter modern biogeochemical cycles, where it may fuel food webs and greenhouse gas emissions in the marine environment.« less
  5. Abstract. Numerical simulations of the Greenland Ice Sheet (GrIS) over geologictimescales can greatly improve our knowledge of the critical factors drivingGrIS demise during climatically warm periods, which has clear relevance forbetter predicting GrIS behavior over the upcoming centuries. To assess thefidelity of these modeling efforts, however, observational constraints ofpast ice sheet change are needed. Across southwestern Greenland, geologicrecords detail Holocene ice retreat across both terrestrial-based and marine-terminating environments, providing an ideal opportunity to rigorouslybenchmark model simulations against geologic reconstructions of ice sheetchange. Here, we present regional ice sheet modeling results using theIce-sheet and Sea-level System Model (ISSM) of Holocene icemore »sheet historyacross an extensive fjord region in southwestern Greenland covering thelandscape around the Kangiata Nunaata Sermia (KNS) glacier and extendingoutward along the 200 km Nuup Kangerula (Godthåbsfjord). Oursimulations, forced by reconstructions of Holocene climate and recentlyimplemented calving laws, assess the sensitivity of ice retreat across theKNS region to atmospheric and oceanic forcing. Our simulations reveal thatthe geologically reconstructed ice retreat across the terrestrial landscapein the study area was likely driven by fluctuations in surface mass balancein response to Early Holocene warming – and was likely not influencedsignificantly by the response of adjacent outlet glaciers to calving andocean-induced melting. The impact of ice calving within fjords, however,plays a significant role by enhancing ice discharge at the terminus, leadingto interior thinning up to the ice divide that is consistent withreconstructed magnitudes of Early Holocene ice thinning. Our results,benchmarked against geologic constraints of past ice-margin change, suggestthat while calving did not strongly influence Holocene ice-margin migrationacross terrestrial portions of the KNS forefield, it strongly impactedregional mass loss. While these results imply that the implementation andresolution of ice calving in paleo-ice-flow models is important towardsmaking more robust estimations of past ice mass change, they also illustratethe importance these processes have on contemporary and future long-term icemass change across similar fjord-dominated regions of the GrIS.« less