Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sea ice simulator”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2015

Is Ice-Rafted Sediment in a North Pole Marine Record Evidence for Perennial Sea-ice Cover?

Ice-rafted sediments of Eurasian and North American origin are found consistently in the upper part (13 Ma BP to present) of the Arctic Coring Expedition (ACEX) ocean core from the Lomonosov Ridge, near the North Pole (approximately 88 degrees N). Based on modern sea-ice drift trajectories and speeds, this has been taken as evidence of the presence of a perennial sea-ice cover in the Arctic Ocean from the middle Miocene onwards. However, other high latitude land and marine records indicate a long-term trend towards cooling broken by periods of extensive warming suggestive of a seasonally ice-free Arctic between the Miocene and the present. We use a coupled sea-ice slab-ocean model including sediment transport tracers to map the spatial distribution of ice-rafted deposits in the Arctic Ocean. We use 6 hourly wind forcing and surface heat fluxes for two different climates: one with a perennial sea-ice cover similar to that of the present day and one with seasonally ice-free conditions, similar to that simulated in future projections. Model results confirm that in the present-day climate, sea ice takes more than 1 year to transport sediment from all its peripheral seas to the North Pole. However, in a warmer climate, sea-ice speeds are significantly faster (for the same wind forcing) and can deposit sediments of Laptev, East Siberian and perhaps also Beaufort Sea origin at the North Pole. This is primarily because of the fact that sea-ice interactions are much weaker with a thinner ice cover and there is less resistance to drift. We conclude that the presence of ice-rafted sediment of Eurasian and North American origin at the North Pole does not imply a perennial sea-ice cover in the Arctic Ocean, reconciling the ACEX ocean core data with other land and marine records.

Palaeoclimate↗
NASA NTRS1996

Projected Impact of Climate Change on the Water and Salt Budgets of the Arctic Ocean by a Global Climate Model

The annual flux of freshwater into the Arctic Ocean by the atmosphere and rivers is balanced by the export of sea ice and oceanic freshwater. Two 150-year simulations of a global climate model are used to examine how this balance might change if atmospheric greenhouse gases (GHGs) increase. Relative to the control, the last 50-year period of the GHG experiment indicates that the total inflow of water from the atmosphere and rivers increases by 10% primarily due to an increase in river discharge, the annual sea-ice export decreases by about half, the oceanic liquid water export increases, salinity decreases, sea-ice cover decreases, and the total mass and sea-surface height of the Arctic Ocean increase. The closed, compact, and multi-phased nature of the hydrologic cycle in the Arctic Ocean makes it an ideal test of water budgets that could be included in model intercomparisons.

Miller, James R.↗
NASA NTRS2007

A New Sea Ice Model: Pancakes to Metrics

We are developing, testing, and validating a new sea ice dynamics model that treats the ice cover as an elastic/decohesive material in the permanent pack and includes the correct frazil/pancake behavior in the marginal zone. Two salient features of present ice dynamics models are that they do not: 1) reproduce the oriented fracture patterns of openings and closings in the pack ice, and 2) accurately model the effects of frazil/pancake ice formation in the ice margin. These poorly modeled areas account for a substantial portion of the ice growth, turbulent heat flux to the atmosphere, salt flux to the ocean, and energy dissipation due to slippage, ridging, and rafting, in the Arctic. Existing sea ice models have shown limited success in predicting the degree to which a lead will open for prescribed or observed forcing conditions. An important aspect of the new model we are developing is that the existence of cracks, along with their orientation, opening, and closing, is predicted. To put this effort in perspective a short history of ice dynamics modeling and data collection is presented. The RGPS data set is used to validate the model. As part of the testing and validation of the model, we are working on a new metric for comparing linear features (leads and ridges) in the data and model to be used in data assimilation for this model. The model framework is presented as well as some results showing the creation and development of leads in a simulation of ice dynamics in the Beaufort Sea. Other presentations by the authors will show other results from this effort.

Coon, M.↗
NASA NTRS2012

The Distribution of Snow Black Carbon observed in the Arctic and Compared to the GISS-PUCCINI Model

In this study, we evaluate the ability of the latest NASA GISS composition-climate model, GISS-E2- PUCCINI, to simulate the spatial distribution of snow BC (sBC) in the Arctic relative to present-day observations. Radiative forcing due to BC deposition onto Arctic snow and sea ice is also estimated. Two sets of model simulations are analyzed, where meteorology is linearly relaxed towards National Centers for Environmental Prediction (NCEP) and towards NASA Modern Era Reanalysis for Research and Applications (MERRA) reanalyses. Results indicate that the modeled concentrations of sBC are comparable with presentday observations in and around the Arctic Ocean, except for apparent underestimation at a few sites in the Russian Arctic. That said, the model has some biases in its simulated spatial distribution of BC deposition to the Arctic. The simulations from the two model runs are roughly equal, indicating that discrepancies between model and observations come from other sources. Underestimation of biomass burning emissions in Northern Eurasia may be the main cause of the low biases in the Russian Arctic. Comparisons of modeled aerosol BC (aBC) with long-term surface observations at Barrow, Alert, Zeppelin and Nord stations show significant underestimation in winter and spring concentrations in the Arctic (most significant in Alaska), although the simulated seasonality of aBC has been greatly improved relative to earlier model versions. This is consistent with simulated biases in vertical profiles of aBC, with underestimation in the lower and middle troposphere but overestimation in the upper troposphere and lower stratosphere, suggesting that the wet removal processes in the current model may be too weak or that vertical transport is too rapid, although the simulated BC lifetime seems reasonable. The combination of observations and modeling provides a comprehensive distribution of sBC over the Arctic. On the basis of this distribution, we estimate the decrease in snow and sea ice albedo and the resulting radiative forcing. We suggest that the albedo reduction due to BC deposition presents significant space-time variations, with highest mean reductions of 1.25% in the Russian Arctic, which are much larger than those in other Arctic regions (0.39% to 0.64 %). The averaged value over the Arctic north of 66degN is 0.4-0.6% during spring, leading to regional surface radiative forcings of 0.7, 1.1 and 1.0Wm(exp−2) in spring 2007, 2008 and 2009, respectively.

Dou, T.↗
NASA NTRS2013

Sensitivity of the Atmospheric Response to Warm Pool El Nino Events to Modeled SSTs and Future Climate Forcings

Warm pool El Nino (WPEN) events are characterized by positive sea surface temperature (SST) anomalies in the central equatorial Pacific. Under present-day climate conditions, WPEN events generate poleward propagating wavetrains and enhance midlatitude planetary wave activity, weakening the stratospheric polar vortices. The late 21st century extratropical atmospheric response to WPEN events is investigated using the Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM), version 2. GEOSCCM simulations are forced by projected late 21st century concentrations of greenhouse gases (GHGs) and ozone-depleting substances (ODSs) and by SSTs and sea ice concentrations from an existing ocean-atmosphere simulation. Despite known ocean-atmosphere model biases, the prescribed SST fields represent a best estimate of the structure of late 21st century WPEN events. The future Arctic vortex response is qualitatively similar to that observed in recent decades but is weaker in late winter. This response reflects the weaker SST forcing in the Nino 3.4 region and subsequently weaker Northern Hemisphere tropospheric teleconnections. The Antarctic stratosphere does not respond to WPEN events in a future climate, reflecting a change in tropospheric teleconnections: The meridional wavetrain weakens while a more zonal wavetrain originates near Australia. Sensitivity simulations show that a strong poleward wavetrain response to WPEN requires a strengthening and southeastward extension of the South Pacific Convergence Zone; this feature is not captured by the late 21st century modeled SSTs. Expected future increases in GHGs and decreases in ODSs do not affect the polar stratospheric responses to WPEN.

stratosphere↗
DOE OSTI2022

Robust but weak winter atmospheric circulation response to future Arctic sea ice loss

The possibility that Arctic sea ice loss weakens mid-latitude westerlies, promoting more severe cold winters, has sparked more than a decade of scientific debate, with apparent support from observations but inconclusive modelling evidence. Here we show that sixteen models contributing to the Polar Amplification Model Intercomparison Project simulate a weakening of mid-latitude westerlies in response to projected Arctic sea ice loss. We develop an emergent constraint based on eddy feedback, which is 1.2 to 3 times too weak in the models, suggesting that the real-world weakening lies towards the higher end of the model simulations. Still, the modelled response to Arctic sea ice loss is weak: the North Atlantic Oscillation response is similar in magnitude and offsets the projected response to increased greenhouse gases, but would only account for around 10% of variations in individual years. We further find that relationships between Arctic sea ice and atmospheric circulation have weakened recently in observations and are no longer inconsistent with those in models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

File Specification for M2AMIP Products: Version Number - 1.0

The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) is an atmospheric reanalysis computed with the Goddard Earth Observing (EOS) System, Version 5.12.4 (GEOS) data assimilation system (Gelaro et al., 2017). To supplement the reanalysis, the GEOS General Circulation Model (GCM) used in MERRA-2 has been used to generate a 10-member ensemble of simulations, configured following the convention of the Atmospheric Model Intercomparison Project (AMIP; Gates et al., 1992). Each ensemble member was initialized using meteorological fields from a different date in November 1979. The AMIP simulations used the sea-surface temperature (SST) and sea-ice boundary conditions that were used in MERRA-2 (Bosilovich et al., 2016). This 10-member ensemble of AMIP simulations, denoted M2AMIP, is available for download in a group of self-describing files, which are documented in this office note. All data collections are provided on the same horizontal grid as MERRA-2. This grid has 576 points in the longitudinal direction and 361 points in the latitudinal direction, corresponding to a resolution of 0.625 degrees by 0.5 degrees. Although data collections are available at this grid, all fields are computed on a cubed-sphere grid with an approximate resolution of 50 km by 50 km and are then spatially interpolated to the latitude-longitude grid. There are no changes in the vertical grids used: variables are provided on either the native vertical grid of 72 model layers, or interpolated to 42 standard pressure levels. Unlike MERRA, no data collections are available at the vertical layer edges. More details on the grid are provided in Section 4. MERRA-2 introduced observation-based precipitation forcing for the land surface parameterization and the corresponding variable PRECTOTCORR in the MERRA-2 FLX (surface turbulent fluxes and related quantities) and LFO (land-surface forcing) collections (see Section 6; Reichle et al., 2017). While this variable is still available for M2AMIP, there was no observation-based forcing, making the value identical to the model derived precipitation, PRECTOT. Similarly, without data assimilation, the values for the analysis increments, D*DTANA, in the tendency and vertically integrated file collections are zero. The M2AMIP data are available for download online through the NASA Center for Climate Simulation (NCCS) DataPortal (https://portal.nccs.nasa.gov/datashare/gmao_m2amip/). Data are arranged in subdirectories based on ensemble member, followed by year and month. Control files that are compatible with the Grid Analysis and Display System (GrADS) are available in the ctl_daily and ctl_monthly directories for the hourly, three hourly, and monthly mean data. Control files for the monthly mean diurnal cycle can be found in the ctl_diurnal subdirectory within the directory for each individual ensemble member.

M2AMIP↗
NASA NTRS2011

The Impact of Warm Pool El Nino Events on Antarctic Ozone

Warm pool El Nino (WPEN) events are characterized by positive sea surface temperature (SST) anomalies in the central equatorial Pacific in austral spring and summer. Previous work found an enhancement in planetary wave activity in the South Pacific in austral spring, and a warming of 3-5 K in the Antarctic lower stratosphere during austral summer, in WPEN events as compared with ENSO neutral. In this presentation, we show that weakening of the Antarctic vortex during WPEN affects the structure and magnitude of high-latitude total ozone. We use total ozone data from TOMS and OMI, as well as station data from Argentina and Antarctica, to identify shifts in the longitudinal location of the springtime ozone minimum from its climatological position. In addition, we examine the sensitivity of the WPEN-related ozone response to the phase of the quasi-biennial oscillation (QBO). We then compare the observed response to WPEN events with Goddard Earth Observing System chemistry-climate model, version 2 (GEOS V2 CCM) simulations. Two, 50-year time-slice simulations are forced by annually repeating SST and sea ice climatologies, one set representing observed WPEN events and the second set representing neutral ENSO events, in a present-day climate. By comparing the two simulations, we isolate the impact of WPEN events on lower stratospheric ozone, and furthermore, examine the sensitivity of the WPEN ozone response to the phase of the QBO.

Hurwitz, Margaret M.↗
NASA NTRS1993

An Arctic source for the Great Salinity Anomaly - A simulation of the Arctic ice-ocean system for 1955-1975

The paper employs a fully prognostic Arctic ice-ocean model to study the interannual variability of sea ice during the period 1955-1975 and to explain the large variability of the ice extent in the Greenland and Iceland seas during the late 1960s. The model is used to test the contention of Aagaard and Carmack (1989) that the Great Salinity Anomaly (GSA) was a consequence of the anomalously large ice export in 1968. The high-latitude ice-ocean circulation changes due to wind field changes are explored. The ice export event of 1968 was the largest in the simulation, being about twice as large as the average and corresponding to 1600 cu km of excess fresh water. The simulations suggest that, besides the above average ice export to the Greenland Sea, there was also fresh water export to support the larger than average ice cover. The model results show the origin of the GSA to be in the Arctic, and support the view that the Arctic may play an active role in climate change.

Hakkinen, Sirpa↗
NASA NTRS2017

Evolving Understanding of Antarctic Ice-Sheet Physics and Ambiguity in Probabilistic Sea-Level Projections

Mechanisms such as ice-shelf hydrofracturing and ice-cliff collapse may rapidly increase discharge from marine-based ice sheets. Here, we link a probabilistic framework for sea-level projections to a small ensemble of Antarctic ice-sheet (AIS) simulations incorporating these physical processes to explore their influence on global-mean sea-level (GMSL) and relative sea-level (RSL). We compare the new projections to past results using expert assessment and structured expert elicitation about AIS changes. Under high greenhouse gas emissions (Representative Concentration Pathway [RCP] 8.5), median projected 21st century GMSL rise increases from 79 to 146 cm. Without protective measures, revised median RSL projections would by 2100 submerge land currently home to 153 million people, an increase of 44 million. The use of a physical model, rather than simple parameterizations assuming constant acceleration of ice loss, increases forcing sensitivity: overlap between the central 90% of simulations for 2100 for RCP 8.5 (93-243 cm) and RCP 2.6 (26-98 cm) is minimal. By 2300, the gap between median GMSL estimates for RCP 8.5 and RCP 2.6 reaches >10 m, with median RSL projections for RCP 8.5 jeopardizing land now occupied by 950 million people (versus 167 million for RCP 2.6). The minimal correlation between the contribution of AIS to GMSL by 2050 and that in 2100 and beyond implies current sea-level observations cannot exclude future extreme outcomes. The sensitivity of post-2050 projections to deeply uncertain physics highlights the need for robust decision and adaptive management frameworks.

Kopp, Robert E.↗
DOE OSTI2021

AMOC stability and diverging response to Arctic sea ice decline in two climate models

This study compares the impacts of Arctic sea ice decline on the Atlantic meridional overturning circulation (AMOC) in two configurations of the Community Earth System Model with different horizontal resolution. In a suite of model experiments, we impose radiative imbalance at the ice surface, replicating a loss of sea ice cover comparable to that observed during 1979–2014, and we find dramatic differences in the AMOC response between the two models. In the lower-resolution configuration, the AMOC weakens by about one-third over the first 100 years, approaching a new quasi-equilibrium. By contrast, in the higher-resolution configuration, the AMOC weakens by ~10% during the first 20–30 years followed by a full recovery driven by invigorated deep water formation in the Labrador Sea and adjacent regions. We investigate these differences using a diagnostic AMOC stability indicator, which reflects the AMOC freshwater transport in and out of the basin and hence the strength of the basin-scale salt-advection feedback. This indicator suggests that the AMOC in the lower-resolution model is less stable and more sensitive to surface perturbations, as confirmed by hosing experiments mimicking Arctic freshening due to sea ice decline. Differences between the models’ mean states, including the Atlantic Ocean mean surface freshwater fluxes, control the differences in AMOC stability. Our results demonstrate that the AMOC stability indicator is indeed useful for evaluating AMOC sensitivity to perturbations. We emphasize that, despite the differences in the long-term adjustment, both models simulate a multidecadal AMOC weakening caused by Arctic sea ice decline, relevant to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

More frequent atmospheric rivers slow the seasonal recovery of Arctic sea ice

In recent decades, Arctic sea ice coverage experienced a drastic decline in winter, when sea ice is expected to recover following the melting season. It is unclear to what extent atmospheric processes such as atmospheric rivers (ARs), intense corridors of moisture transport, contribute to this reduced recovery of sea ice. Here, in this study, using observations and climate model simulations we find a robust frequency increase in ARs in early winter over the Barents-Kara Seas and the central Arctic (ABK) between 1979 and 2021. The moisture carried by more frequent ARs has intensified surface downward longwave radiation and rainfall, caused stronger melting of thin, fragile ice cover, and slowed the seasonal recovery of sea ice, accounting for 34% of the sea ice cover decline in ABK. A series of model ensemble experiments suggests that, in addition to a uniform AR increase in response to anthropogenic warming, tropical Pacific variability also contributes to the observed Arctic AR changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Ocean and Sea Ice and their Interactions around Greenland and the West Antarctic Peninsula in Forced Fine-Resolution Global Simulations

Anomalously warm ocean waters have been implicated as the cause of accelerated ablation along the margin of the Greenland ice sheet, as well as the cause of increasing basal melt rate, mass loss, and grounding line retreat of many Antarctic ice shelves. The overarching science objective of this project is to investigate how the interplay of regional processes and decadal changes in local and remote forcing impact the delivery and end member composition of waters over the continental shelves of Greenland and Antarctica. To address this goal, an existing global eddy-active ocean general circulation model (OGCM), coupled to a thermodynamic/dynamic ice model and forced with atmospheric reanalysis fluxes from 1948-2009, was analyzed to understand these processes around the continental margins of Greenland and Antarctica. This simulation did not contain a representation of ice sheet melt or resolve mesoscale eddies over these continental margins, so two new coupled OGCM/sea-ice models were constructed and analyzed to address the roles of small scale features and processes (< 5km) and freshwater flux surpluses.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Components of the ice age circulation

The effects of ice age boundary conditions on atmospheric dynamics and regional climate patterns are investigated using four GCM simulations. Particular consideration is given to sea surface temperature-sea ice distribution, the appearance of land ice, and the increased elevation of land ice. It is observed that the ice-age sea surface temperature stabilizes the atmosphere over the oceans, increases the frequency of storm tracking through central North America, and amplifies transient eddy energy without increasing baroclinic generation. It is detected that low-elevation ice generates low pressure over eastern North America and southern Europe in winter, while increasing cloud cover and cooling the land in summer. Elevation of the ice sheets cools the land in winter, further intensifies storms off northeastern North America, induces subsidence warming downstream of the European ice sheets in summer, and increases the transient and stationary eddy energy through increased baroclinicity.

Rind, D.↗
NASA NTRS1994

The Area Coverage of Geophysical Fields as a Function of Sensor Field-of View

In many remote sensing studies of geophysical fields such as clouds, land cover, or sea ice characteristics, the fractional area coverage of the field in an image is estimated as the proportion of pixels that have the characteristic of interest (i.e., are part of the field) as determined by some thresholding operation. The effect of sensor field-of-view on this estimate is examined by modeling the unknown distribution of subpixel area fraction with the beta distribution, whose two parameters depend upon the true fractional area coverage, the pixel size, and the spatial structure of the geophysical field. Since it is often not possible to relate digital number, reflectance, or temperature to subpixel area fraction, the statistical models described are used to determine the effect of pixel size and thresholding operations on the estimate of area fraction for hypothetical geophysical fields. Examples are given for simulated cumuliform clouds and linear openings in sea ice, whose spatial structures are described by an exponential autocovariance function. It is shown that the rate and direction of change in total area fraction with changing pixel size depends on the true area fraction, the spatial structure, and the thresholding operation used.

Key, Jeffrey R.↗
NASA NTRS2022

Seasonal Tidewater Glacier Terminus Oscillations Bias Multi-Decadal Projections of Ice Mass Change

Numerical, process-based simulations of tidewater glacier evolution are necessary to project future sea-level change under various climate scenarios. Previous work has shown that nonlinearities in tidewater glacier and ice stream dynamics can lead to biases in simulated ice mass change in the presence of noisy forcings. Ice sheet modeling projections that are used in the IPCC Assessment Report 6 (AR6) utilize atmospheric and oceanic forcings at annual temporal resolution, omitting any higher frequency forcings. Here, we quantify the effect of seasonal (<1 year) tidewater glacier terminus oscillations on decadal-scale (30 years) mass change. We use an idealized geometry to mimic realistic tidewater glacier geometries and investigate the impact of the magnitude of seasonal oscillations, bed slope at the glacier terminus, and basal friction law. We find that omitting seasonal terminus motion results in biased mass change projections, with up to an 39% underestimate or 25% overestimate of mass loss when seasonality is neglected, depending on bed slope and magnitude of oscillations. The bias is most sensitive to the magnitude of the seasonal terminus oscillations and exhibits very little sensitivity to choice of friction law. Although omitting seasonal terminus motion may bias mass change projections of individual glaciers, further work is needed to investigate the impact on ice sheet mass change projections. In order to achieve this, seasonality in atmospheric and oceanic forcings must be adequately represented and observations of seasonal terminus positions and tidewater glacier thickness changes must be acquired to evaluate numerical models. Note: Zip folder contains all videos listed in the supplemental document and can be downloaded for viewing.

D. Felikson↗