Massachusetts Institute of Technology
Stochastic Sea Ice Modeling with the Dynamically Orthogonal Equations
Abstract
dc:description.abstractAccurate numerical models are essential to predict the complex evolution of rapidly changing sea ice conditions and study impacts on climate and navigation. However, sea ice models contain uncertainties associated with initial conditions and forcing (wind, ocean), as well as with parameter values, functional forms of the constitutive relations, and state variables themselves, all of which limit predictive capabilities. Due to the multiple types and scales of sea ice and the complex nonlinear mechanics and high dimensionality of differential equations, efficient ocean and sea ice probabilistic modeling, Bayesian inversion, and machine learning are challenging. In this work, we implement a deterministic 2D viscoplastic sea ice solver and derive and implement new sea ice probabilistic models based on the dynamically orthogonal (DO) equations. We focus on the stochastic two-dimensional sea ice momentum equations with nonlinear viscoplastic constitutive law. We first implement and verify a deterministic 2D viscoplastic sea ice solver. Next, we derive the new stochastic Sea Ice Dynamically Orthogonal equations and develop numerical schemes for their solution. These equations and schemes preserve nonlinearities in the underlying spatiotemporal dynamics and evolve the non-Gaussianity of the statistics. We evaluate and illustrate the new stochastic sea ice modeling and schemes using idealized stochastic test cases. We employ two stochastic test cases with different types of sea ice: ice sheets and frozen ice cover with uncertain initial velocities. We showcase the ability to evolve non-Gaussian statistics and capture complex nonlinear dynamics efficiently. We study the convergence to the physical discretization, and stochastic convergence to the stochastic subspace size and coefficient samples. Finally, we assess and show significant computational and memory efficiency compared to the direct Monte Carlo method.
Degree
thesis:*- Name thesis:degree_name
- Master
- Department dc:contributor.department
- Massachusetts Institute of Technology. Department of Mechanical Engineering
- Grantor dc:publisher
- Massachusetts Institute of Technology
- Year dc:date.issued
- 2023
Author and committee
dc:creator, dc:contributor.*- Author dc:creator
-
- Suresh Babu, Anantha Narayanan
- Advisor dc:contributor.advisor
-
- Lermusiaux, Pierre F.J.
Rights
dc:rights- Statement dc:rights
-
- In Copyright - Educational Use Permitted
- Copyright retained by author(s)
- Licence dc:rights.uri
Identifiers
dc:identifier.*- Handle dc:identifier.uri
- https://hdl.handle.net/1721.1/152699
- OAI identifier oai:identifier
- oai:dspace.mit.edu:1721.1/152699