skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: The Potential of the Ensemble Kalman Filter to Improve Glacier Modeling
Using a simplified two-stage ice sheet model, we explore the potential of statistical data assimilation methods to improve predictions of glacier melt, which has significant implications for reducing uncertainty in projections of sea level rise. Through twin experiments utilizing artificial data, we find that the ensemble Kalman filter improves simulations of glacier evolution initialized with incorrect initial conditions and parameters, providing us with better predictions of future glacier melt. We explore the number of observations necessary to produce an accurate model run. We also explore optimal observation assimilation schemes, and determine that deviations from the true glacier response that stem from having few data points in the pre-satellite era can be corrected with modern observation data. Our results show that statistical data assimilation methods have great potential to improve complex glacier models using real-world observations.  more » « less
Award ID(s):
2051019
PAR ID:
10542783
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Springer
Date Published:
Journal Name:
La Matematica
Volume:
3
Issue:
3
ISSN:
2730-9657
Page Range / eLocation ID:
1085 to 1102
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al (, The Cryosphere)
    Abstract. The discovery of Antarctica's deepest subglacial troughbeneath the Denman Glacier, combined with high rates of basal melt at thegrounding line, has caused significant concern over its vulnerability toretreat. Recent attention has therefore been focusing on understanding thecontrols driving Denman Glacier's dynamic evolution. Here we consider theShackleton system, comprised of the Shackleton Ice Shelf, Denman Glacier,and the adjacent Scott, Northcliff, Roscoe and Apfel glaciers, about whichalmost nothing is known. We widen the context of previously observed dynamicchanges in the Denman Glacier to the wider region of the Shackleton system,with a multi-decadal time frame and an improved biannual temporal frequencyof observations in the last 7 years (2015–2022). We integrate newsatellite observations of ice structure and airborne radar data with changesin ice front position and ice flow velocities to investigate changes in thesystem. Over the 60-year period of observation we find significant riftpropagation on the Shackleton Ice Shelf and Scott Glacier and notablestructural changes in the floating shear margins between the ice shelf andthe outlet glaciers, as well as features indicative of ice with elevatedsalt concentration and brine infiltration in regions of the system. Over theperiod 2017–2022 we observe a significant increase in ice flow speed (up to50 %) on the floating part of Scott Glacier, coincident with small-scalecalving and rift propagation close to the ice front. We do not observe anyseasonal variation or significant change in ice flow speed across the restof the Shackleton system. Given the potential vulnerability of the system toaccelerating retreat into the overdeepened, potentially sediment-filledbedrock trough, an improved understanding of the glaciological,oceanographic and geological conditions in the Shackleton system arerequired to improve the certainty of numerical model predictions, and weidentify a number of priorities for future research. With access to theseremote coastal regions a major challenge, coordinated internationallycollaborative efforts are required to quantify how much the Shackletonregion is likely to contribute to sea level rise in the coming centuries. 
    more » « less
  2. ; ; (, Geochemistry, Geophysics, Geosystems)
    Abstract This study is motivated by the observed variability in trace element isotopic and chemical compositions of primitive (Si52 wt %) basalts in southwest North America (SWNA) during the Cenozoic transition from subduction to extension. Specifically, we focus on processes that may explain the enigmatic observation that in some localities, basalts with low Ta/Th, consistent with parental melts in a subduction setting, have signatures consistent with continental lithospheric mantle (CLM). In locations with the oldest CLM (Proterozoic and Archean), Cenozoic basalts with low Ta/Th have well below zero. We model channelized magma transport through the CLM using simple 1D transport models to explore the extent to which diffusive and reactive mass exchange can modify Nd isotopic compositions via open system melt‐wallrock interactions. For geologically reasonable channel spacings and volume fractions, we quantify the reactive assimilation rates required for incoming melt with a different than the wall‐rock to undergo a substantial isotopic shift during 10 km channelized melt transport. In the presence of grain boundaries, enhanced diffusion between melt‐rich channels and melt‐poor surrounding rock contributes to isotopic equilibration, however this effect is not enough to explain observations; our models suggest a significant contribution from reactive assimilation of wall‐rock. Additionally our models support the idea that the observed covariability in Ta/Th and in Cenozoic basalts cannot be attributed to transport alone and must also reflect the transition from subduction‐related to extension‐related parental melts in SWNA. 
    more » « less
  3. ; ; ; (, Chaos: An Interdisciplinary Journal of Nonlinear Science)
    The reconstruction of electrical excitation patterns through the unobserved depth of the tissue is essential to realizing the potential of computational models in cardiac medicine. We have utilized experimental optical-mapping recordings of cardiac electrical excitation on the epicardial and endocardial surfaces of a canine ventricle as observations directing a local ensemble transform Kalman filter data assimilation scheme. We demonstrate that the inclusion of explicit information about the stimulation protocol can marginally improve the confidence of the ensemble reconstruction and the reliability of the assimilation over time. Likewise, we consider the efficacy of stochastic modeling additions to the assimilation scheme in the context of experimentally derived observation sets. Approximation error is addressed at both the observation and modeling stages through the uncertainty of observations and the specification of the model used in the assimilation ensemble. We find that perturbative modifications to the observations have marginal to deleterious effects on the accuracy and robustness of the state reconstruction. Furthermore, we find that incorporating additional information from the observations into the model itself (in the case of stimulus and stochastic currents) has a marginal improvement on the reconstruction accuracy over a fully autonomous model, while complicating the model itself and thus introducing potential for new types of model errors. That the inclusion of explicit modeling information has negligible to negative effects on the reconstruction implies the need for new avenues for optimization of data assimilation schemes applied to cardiac electrical excitation. 
    more » « less
  4. ; ; ; ; (, Nature Geoscience)
    Abstract Feedbacks between ice melt, glacier flow and ocean circulation can rapidly accelerate ice loss at tidewater glaciers and alter projections of sea-level rise. At the core of these projections is a model for ice melt that neglects the fact that glacier ice contains pressurized bubbles of air due to its formation from compressed snow. Current model estimates can underpredict glacier melt at termini outside the region influenced by the subglacial discharge plume by a factor of 10–100 compared with observations. Here we use laboratory-scale experiments and theoretical arguments to show that the bursting of pressurized bubbles from glacier ice could be a source of this discrepancy. These bubbles eject air into the seawater, delivering additional buoyancy and impulses of turbulent kinetic energy to the boundary layer, accelerating ice melt. We show that real glacier ice melts 2.25 times faster than clear bubble-free ice when driven by natural convection in a laboratory setting. We extend these results to the geophysical scale to show how bubble dynamics contribute to ice melt from tidewater glaciers. Consequently, these results could increase the accuracy of modelled predictions of ice loss to better constrain sea-level rise projections globally. 
    more » « less
  5. ; ; (, Monthly Weather Review)
    Abstract For data assimilation to provide faithful state estimates for dynamical models, specifications of observation uncertainty need to be as accurate as possible. Innovation-based methods based on Desroziers diagnostics, are commonly used to estimate observation uncertainty, but such methods can depend greatly on the prescribed background uncertainty. For ensemble data assimilation, this uncertainty comes from statistics calculated from ensemble forecasts, which require inflation and localization to address under sampling. In this work, we use an ensemble Kalman filter (EnKF) with a low-dimensional Lorenz model to investigate the interplay between the Desroziers method and inflation. Two inflation techniques are used for this purpose: 1) a rigorously tuned fixed multiplicative scheme and 2) an adaptive state-space scheme. We document how inaccuracies in observation uncertainty affect errors in EnKF posteriors and study the combined impacts of misspecified initial observation uncertainty, sampling error, and model error on Desroziers estimates. We find that whether observation uncertainty is over- or underestimated greatly affects the stability of data assimilation and the accuracy of Desroziers estimates and that preference should be given to initial overestimates. Inline estimates of Desroziers tend to remove the dependence between ensemble spread–skill and the initially prescribed observation error. In addition, we find that the inclusion of model error introduces spurious correlations in observation uncertainty estimates. Further, we note that the adaptive inflation scheme is less robust than fixed inflation at mitigating multiple sources of error. Last, sampling error strongly exacerbates existing sources of error and greatly degrades EnKF estimates, which translates into biased Desroziers estimates of observation error covariance. Significance StatementTo generate accurate predictions of various components of the Earth system, numerical models require an accurate specification of state variables at our current time. This step adopts a probabilistic consideration of our current state estimate versus information provided from environmental measurements of the true state. Various strategies exist for estimating uncertainty in observations within this framework, but are sensitive to a host of assumptions, which are investigated in this study. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email