Search NASA⌕ Search

SEARCH · Search NASA

Results for “ice models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Validation of sea ice models using an uncertainty-based distance metric for multiple model variables: NEW METRIC FOR SEA ICE MODEL VALIDATION

Here, we implement a variance-based distance metric (D n ) to objectively assess skill of sea ice models when multiple output variables or uncertainties in both model predictions and observations need to be considered. The metric compares observations and model data pairs on common spatial and temporal grids improving upon highly aggregated metrics (e.g., total sea ice extent or volume) by capturing the spatial character of model skill. The D n metric is a gamma-distributed statistic that is more general than the χ 2 statistic commonly used to assess model fit, which requires the assumption that the model is unbiased and can only incorporate observational error in the analysis. The D n statistic does not assume that the model is unbiased, and allows the incorporation of multiple observational data sets for the same variable and simultaneously for different variables, along with different types of variances that can characterize uncertainties in both observations and the model. This approach represents a step to establish a systematic framework for probabilistic validation of sea ice models. The methodology is also useful for model tuning by using the D n metric as a cost function and incorporating model parametric uncertainty as part of a scheme to optimize model functionality. We apply this approach to evaluate different configurations of the standalone Los Alamos sea ice model (CICE) encompassing the parametric uncertainty in the model, and to find new sets of model configurations that produce better agreement than previous configurations between model and observational estimates of sea ice concentration and thickness.

58 GEOSCIENCES↗

The CMIP6 Sea-Ice Model Intercomparison Project (SIMIP): Understanding sea ice through climate-model simulations

A better understanding of the role of sea ice for the changing climate of our planet is the central aim of the diagnostic Coupled Model Intercomparison Project 6 (CMIP6)-endorsed Sea-Ice Model Intercomparison Project (SIMIP). To reach this aim, SIMIP requests sea-ice-related variables from climate-model simulations that allow for a better understanding and, ultimately, improvement of biases and errors in sea-ice simulations with large-scale climate models. This then allows us to better understand to what degree CMIP6 model simulations relate to reality, thus improving our confidence in answering sea-ice-related questions based on these simulations. Furthermore, the SIMIP protocol provides a standard for sea-ice model output that will streamline and hence simplify the analysis of the simulated sea-ice evolution in research projects independent of CMIP. To reach its aims, SIMIP provides a structured list of model output that allows for an examination of the three main budgets that govern the evolution of sea ice, namely the heat budget, the momentum budget, and the mass budget. Furthermore, we explain the aims of SIMIP in more detail and outline how its design allows us to answer some of the most pressing questions that sea ice still poses to the international climate-research community.

54 ENVIRONMENTAL SCIENCES↗

A New Discrete Element Sea-Ice Model for Earth System Modeling

Sea ice forms a frozen crust of sea water oating in high-latitude oceans. It is a critical component of the Earth system because its formation helps to drive the global thermohaline circulation, and its seasonal waxing and waning in the high north and Southern Ocean signi cantly affects planetary albedo. Usually 4{6% of Earth's marine surface is covered by sea ice at any one time, which limits the exchange of heat, momentum, and mass between the atmosphere and ocean in the polar realms. Snow accumulates on sea ice and inhibits its vertical growth, increases its albedo, and contributes to pooled water in melt ponds that darken the Arctic ice surface in the spring. Ice extent and volume are subject to strong seasonal, inter-annual and hemispheric variations, and climatic trends, which Earth System Models (ESMs) are challenged to simulate accurately (Stroeve et al., 2012; Stocker et al., 2013). This is because there are strong coupled feedbacks across the atmosphere-ice-ocean boundary layers, including the ice-albedo feedback, whereby a reduced ice cover leads to increased upper ocean heating, further enhancing sea-ice melt and reducing incident solar radiation re ected back into the atmosphere (Perovich et al., 2008). A reduction in perennial Arctic sea-ice during the satellite era has been implicated in mid-latitude weather changes, including over North America (Overland et al., 2015). Meanwhile, most ESMs have been unable to simulate observed inter-annual variability and trends in Antarctic sea-ice extent during the same period (Gagne et al., 2014).

54 ENVIRONMENTAL SCIENCES↗

CICE, The Los Alamos Sea Ice Model

The Los Alamos sea ice model (CICE) is the result of an effort to develop a computationally efficient sea ice component for a fully coupled atmosphere–land–ocean–ice global climate model. It was originally designed to be compatible with the Parallel Ocean Program (POP), an ocean circulation model developed at Los Alamos National Laboratory for use on massively parallel computers. CICE has several interacting components: a vertical thermodynamic model that computes local growth rates of snow and ice due to vertical conductive, radiative and turbulent fluxes, along with snowfall; an elastic-viscous-plastic model of ice dynamics, which predicts the velocity field of the ice pack based on a model of the material strength of the ice; an incremental remapping transport model that describes horizontal advection of the areal concentration, ice and snow volume and other state variables; and a ridging parameterization that transfers ice among thickness categories based on energetic balances and rates of strain. It also includes a biogeochemical model that describes evolution of the ice ecosystem. The CICE sea ice model is used for climate research as one component of complex global earth system models that include atmosphere, land, ocean and biogeochemistry components. It is also used for operational sea ice forecasting in the polar regions and in numerical weather prediction models.

Hunke, Elizabeth↗

CICE, The Los Alamos Sea Ice Model

The Los Alamos sea ice model (CICE) is the result of an effort to develop a computationally efficient sea ice component for a fully coupled atmosphere–land–ocean–ice global climate model. It was originally designed to be compatible with the Parallel Ocean Program (POP), an ocean circulation model developed at Los Alamos National Laboratory for use on massively parallel computers. CICE has several interacting components: a vertical thermodynamic model that computes local growth rates of snow and ice due to vertical conductive, radiative and turbulent fluxes, along with snowfall; an elastic-viscous-plastic model of ice dynamics, which predicts the velocity field of the ice pack based on a model of the material strength of the ice; an incremental remapping transport model that describes horizontal advection of the areal concentration, ice and snow volume and other state variables; and a ridging parameterization that transfers ice among thickness categories based on energetic balances and rates of strain. It also includes a biogeochemical model that describes evolution of the ice ecosystem. The CICE sea ice model is used for climate research as one component of complex global earth system models that include atmosphere, land, ocean and biogeochemistry components. It is also used for operational sea ice forecasting in the polar regions and in numerical weather prediction models.

Hunke, Elizabeth↗

Quality control for community-based sea-ice model development

A new collaborative organization for sea-ice model development, the CICE Consortium, has devised quality control procedures to maintain the integrity of its numerical codes' physical representations, enabling broad participation from the scientific community in the Consortium's open software development environment. Using output from five coupled and uncoupled configurations of the Los Alamos Sea Ice Model, CICE, we formulate quality control methods that exploit common statistical properties of sea-ice thickness, and test for significant changes in model results in a computationally efficient manner. New additions and changes to CICE are graded into four categories, ranging from bit-for-bit amendments to significant, answer-changing upgrades. These modifications are assessed using criteria that account for the high level of autocorrelation in sea-ice time series, along with a quadratic skill metric that searches for hemispheric changes in model answers across an array of different CICE configurations. These metrics also provide objective guidance for assessing new physical representations and code functionality.

58 GEOSCIENCES↗

Quantifying uncertainty and sensitivity in sea ice models

The Los Alamos Sea Ice model has a number of input parameters for which accurate values are not always well established. We conduct a variance-based sensitivity analysis of hemispheric sea ice properties to 39 input parameters. The method accounts for non-linear and non-additive effects in the model.

54 ENVIRONMENTAL SCIENCES↗

Deviation from the dipole-ice model in the spinel spin-ice candidate MgEr 2 Se 4

In spin-ice research, small variations in structure or interactions drive a multitude of different behaviors, yet the collection of known materials relies heavily on the “227” pyrochlore structure. Here, we present thermodynamic, structural, and inelastic-neutron-scattering data on a spin-ice material, MgEr 2 Se 4 , which contributes to the relatively underexplored family of rare-earth spinel chalcogenides. X-ray and neutron diffraction confirm a normal spinel structure and place Er 3+ moments on an ideal pyrochlore sublattice. Measurement of crystal electric field excitations with inelastic neutron scattering confirms that the moments have perfect Ising character, and further identifies the ground-state Kramers doublet as having dipolar-octupolar form with a significant multipolar character. Heat capacity and magnetic neutron diffuse scattering have icelike features, but are inconsistent with Monte Carlo simulations of the nearest-neighbor and next-nearest-neighbor dipolar spin-ice (DSI) models. A significant remnant entropy is observed as T → 0 K, but again falls short of the full Pauling expectation for DSI, unless significant disorder is added. We show that these observations are fully in line with what has been recently reported for CdEr 2 Se 4 , and point to the importance of quantum fluctuations in these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Should Sea-Ice Modeling Tools Designed for Climate Research Be Used for Short-Term Forecasting?

In theory, the same sea-ice models could be used for both research and operations, but in practice, differences in scientific and software requirements and computational and human resources complicate the matter. Although sea-ice modeling tools developed for climate studies and other research applications produce output of interest to operational forecast users, such as ice motion, convergence, and internal ice pressure, the relevant spatial and temporal scales may not be sufficiently resolved. For instance, sea-ice research codes are typically run with horizontal resolution of more than 3 km, while mariners need information on scales less than 300 m. Certain sea-ice processes and coupled feedbacks that are critical to simulating the Earth system may not be relevant on these scales; and therefore, the most important model upgrades for improving sea-ice predictions might be made in the atmosphere and ocean components of coupled models or in their coupling mechanisms, rather than in the sea-ice model itself. This paper discusses some of the challenges in applying sea-ice modeling tools developed for research purposes for operational forecasting on short time scales, and highlights promising new directions in sea-ice modeling.

54 ENVIRONMENTAL SCIENCES↗

Hybrid (PDE+ML) models in the context of land ice modeling

Focal Area(s): Primary focal areas are: predictive modeling through the use of AI-derived model components; advanced methods including network design/optimization/deep learning. Science Challenge: The atmospheric, ocean and ice dynamics components of the Energy Exascale Earth System Model (E3SM) are governed by Partial Differential Equations (PDEs) and significant efforts have been made during the last decades to develop such computational models. Here we propose to fundamentally improve these PDE-based codes by enhancing them with Machine Learning (ML) sub-models for complex, poorly understood physical processes in the context of ice sheet modeling. We propose to train these models with a novel approach that allows the assimilation of the different sources of data available (direct/indirect observations and possibly simulation data), improving on existing simplified models. We also highlight computational challenges originating from the coexistence of PDE-based and ML-based models.

54 ENVIRONMENTAL SCIENCES↗

A modified thermodynamic sea ice model and its application

This paper presents a modified thermodynamic sea ice model based on the Winton’s three-layer model framework. Several improvements focused on the vertical thermodynamics have been made in the model. Results from a series of one-dimensional experiments show that equilibrium ice thickness in our modified model is increased by about 45 cm when compared with the original Winton’s model. Sensitivity tests indicate that all modifications mentioned above contribute to this change of ice thickness, among which the increase of ice layer has the most significant effect followed by the sea ice conductivity parameterization. The modified thermodynamic sea ice model is then coupled into the sea ice component SIS in GFDL MOM4 to exert long-term integration experiments forced by CORE2.0 data. Results indicate that the modified model shows significant improvement in simulating the Arctic sea ice, including an increase in both the sea ice thickness and volume over the whole Arctic region. This confirms the above founding from 1-D simulations.

54 ENVIRONMENTAL SCIENCES↗

Simulating Sea-Ice Deformation in Viscous-Plastic Sea-Ice Models With CD-Grids

Linear kinematic features (LKFs) are found everywhere in the Arctic sea-ice cover. They are strongly localized deformations often associated with the formation of leads and pressure ridges. In viscous-plastic (VP) sea-ice models, the simulation of LKFs depends on several factors such as the grid resolution, the numerical solver convergence, and the placement of the variables on the mesh. In this study, we compare two recently proposed discretization with a CD-grid placement with respect to their ability to reproduce LKFs. The first (CD1) is based on a nonconforming finite element discretization, whereas the second (CD2) uses a conforming subgrid discretization. To analyze their resolution properties, we evaluate runs from different models (e.g., FESOM, MPAS) on a benchmark problem using quadrilateral, hexagonal and triangular meshes. Our findings show that the CD1 setup simulates more deformation structure than the CD2 setup. This highlights the importance of the type of spatial discretization for the simulation of LKFs. Due to the higher number of degrees of freedom, both CD-grids resolve more LKFs than traditional A, B, and C-grids at fixed mesh level. This is an advantage of the CD-grid approach, as high spatial mesh resolution is needed in VP sea-ice models to simulate LKFs.

54 ENVIRONMENTAL SCIENCES↗

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ~2 km, as opposed to ~0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ~2 km (arguably 4 km) when the connectivity is low or moderate and ~1 km (arguably 2 km) when the connectivity is strong.

54 ENVIRONMENTAL SCIENCES↗

Impacts of a mushy‐layer thermodynamic approach in global sea‐ice simulations using the CICE sea‐ice model

Abstract We perform global simulations of the Los Alamos sea‐ice model, CICE, with a new thermodynamics component that has a fully prognostic, variable bulk salinity vertical profile based on mushy layer physics. The processes of gravity drainage, melt‐water flushing and snow‐ice formation are parameterized to allow the bulk salinity to evolve with time. We analyze the seasonal and spatial variation of sea‐ice bulk salinity, area, volume and thickness and compare these quantities to simulations using the previous thermodynamic component. Adjusting one of the gravity drainage parameters, we find good agreement between simulation results and fieldwork ice‐core observations of sea‐ice bulk salinity. As expected, bulk salinity is highest during periods of ice growth and lowest after periods of ice melt. In the northern hemisphere the new thermodynamics component produces thicker ice than the previous thermodynamics component. Of the nine major differences between the two models, differences in how salinities are calculated and how melt‐pond flushing is parameterized are the principal causes of this thickness difference. Thickness differences are smaller in the southern hemisphere than in the northern hemisphere since a greater fraction of ice melts, and differences cannot accumulate year‐on‐year. Model differences in how ice thickness changes and snow‐ice formation are calculated are the most important causes of the different thickness between the two thermodynamic components in the southern hemisphere. The melt‐pond area and volume are found to be highly sensitive to a parameter choice controlling drainage through macroscopic holes in the ice, in both hemispheres.

Turner, Adrian K.↗

Uncertainty quantification and global sensitivity analysis of the Los Alamos sea ice model

Changes in the high-latitude climate system have the potential to affect global climate through feedbacks with the atmosphere and connections with midlatitudes. Sea ice and climate models used to understand these changes have uncertainties that need to be characterized and quantified. We present a quantitative way to assess uncertainty in complex computer models, which is a new approach in the analysis of sea ice models. We characterize parametric uncertainty in the Los Alamos sea ice model (CICE) in a standalone configuration and quantify the sensitivity of sea ice area, extent, and volume with respect to uncertainty in 39 individual model parameters. Unlike common sensitivity analyses conducted in previous studies where parameters are varied one at a time, this study uses a global variance-based approach in which Sobol' sequences are used to efficiently sample the full 39-dimensional parameter space. We implement a fast emulator of the sea ice model whose predictions of sea ice extent, area, and volume are used to compute the Sobol' sensitivity indices of the 39 parameters. Main effects and interactions among the most influential parameters are also estimated by a nonparametric regression technique based on generalized additive models. A ranking based on the sensitivity indices indicates that model predictions are most sensitive to snow parameters such as snow conductivity and grain size, and the drainage of melt ponds. Lastly, it is recommended that research be prioritized toward more accurately determining these most influential parameter values by observational studies or by improving parameterizations in the sea ice model.

54 ENVIRONMENTAL SCIENCES↗

Overestimated Eddy Kinetic Energy in the Eddy‐Rich Regions Simulated by Eddy‐Resolving Global Ocean–Sea Ice Models

Abstract The performance of eddy‐resolving global ocean–sea ice models in simulating mesoscale eddies is evaluated using six eddy‐resolving experiments forced by different atmospheric reanalysis products. Interestingly, eddy‐resolving ocean general circulation models (OGCMs) tend to simulate more (less) energetic eddy‐rich (eddy‐poor) regions with a smaller (larger) spatial extent than satellite observation, which finally shows that larger (smaller) mesoscale energy intensity (EI) is simulated in the eddy‐rich (eddy‐poor) regions. Quantitatively, there is an approximately 27%–60% overestimation of EI in the eddy‐rich regions, which are mainly located in the Kuroshio–Oyashio Extension, the Gulf Stream, and the Antarctic Circumpolar Currents regions, although the global mean EI is underestimated by 25%–45%. Apparently, the eddy kinetic energy in the eddy‐poor region is underestimated. Further analyses based on coherent mesoscale eddy properties show that the overestimation in the eddy‐rich regions is mainly attributed to mesoscale eddies’ intensity and is more prominent when mesoscale eddies are in their growth stage.

54 ENVIRONMENTAL SCIENCES↗