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 19 records

Implementation of Interactive Sea Ice Radiative Transfer into the GISS GCM and Its Impact on Sea Ice Simulation

To improve the representation of solar radiation in sea ice in the Goddard Institute for Space Studies (GISS) climate model (ModelE), we have implemented an advanced radiative transfer (RT) scheme to interactively calculate the surface albedo and transmittance of the ice layer and the solar absorption within. The ice physical properties (temperature, salinity, and density), which are simulated by the thermodynamic module in ModelE, are directly used as the input of the RT scheme to determine the inherent optical properties (IOPs) of brine pockets and air bubbles in the ice layers. This provides a more sophisticated, accurate, and complete treatment of sea ice RT for ModelE. The sea ice RT model is also validated with measured spectral albedo and transmittance for different ice types. Compared with ModelE simulations with current radiative parameterization for sea ice, the new RT implementation will slightly increase the ice surface albedo in the summer in the Arctic and greatly increase the radiative transmittance. Consequently, the ice surface melting is reduced and the base melting is increased but the difference in the simulated ice area is small. The effect of this new RT implementation on the Antarctic ice is negligible. The observed decrease in Arctic sea ice extent is well captured by the ModelE simulation.

Sea ice↗

Arctic Sea Ice Simulation in the PlioMIP Ensemble

Eight general circulation models have simulated the mid-Pliocene warm period (mid-Pliocene, 3.264 to 3.025 Ma) as part of the Pliocene Modelling Intercomparison Project (PlioMIP). Here, we analyse and compare their simulation of Arctic sea ice for both the pre-industrial period and the mid-Pliocene. Mid-Pliocene sea ice thickness and extent is reduced, and the model spread of extent is more than twice the pre-industrial spread in some summer months. Half of the PlioMIP models simulate ice-free conditions in the mid-Pliocene. This spread amongst the ensemble is in line with the uncertainties amongst proxy reconstructions for mid-Pliocene sea ice extent. Correlations between mid-Pliocene Arctic temperatures and sea ice extents are almost twice as strong as the equivalent correlations for the pre-industrial simulations. The need for more comprehensive sea ice proxy data is highlighted, in order to better compare model performances.

simulation↗

The Effects of Snow Depth Forcing on Southern Ocean Sea Ice Simulations

The spatial and temporal distribution of snow on sea ice is an important factor for sea ice and climate models. First, it acts as an efficient insulator between the ocean and the atmosphere, and second, snow is a source of fresh water for altering the already weak Southern Ocean stratification. For the Antarctic, where the ice thickness is relatively thin, snow can impact the ice thickness in two ways: a) As mentioned above snow on sea ice reduces the ocean-atmosphere heat flux and thus reduces freezing at the base of the ice flows; b) a heavy snow load can suppress the ice below sea level which causes flooding and, with subsequent freezing, a thickening of the sea ice (snow-to-ice conversion). In this paper, we compare different snow fall paramterizations (incl. the incorporation of satellite-derived snow depth) and study the effect on the sea ice using a sea ice model.

Powel, Dylan C.↗

Sea ice simulations based on fields generated by the GLAS GCM

The GLAS General Circulation Model (GCM) was applied to the four-month simulation of the thermodynamic part of the Parkinson-Washington sea ice model using atmospheric boundary conditions. The sea ice thickness and distribution were predicted for the Jan. 1-Apr. 30 period using the GCM-fields of solar and infrared radiation, specific humidity and air temperature at the surface, and snow accumulation; the sensible heat and evaporative surface fluxes were consistent with the ground temperatures produced by the ice model and the air temperatures determined by the atmospheric concept. It was concluded that the Parkinson-Washington sea ice model results in acceptable ice concentrations and thicknesses when used with GLAS GCM for the Jan.-Apr. period suggesting the feasibility of fully coupled ice-atmosphere simulations with these two approaches.

Parkinson, C. L.↗

Evaluation of Arctic Sea Ice Thickness Simulated by AOMIP Models

We compare results from six AOMIP model simulations with estimates of sea ice thickness obtained from ICESat, moored and submarine-based upward looking sensors, airborne electromagnetic measurements and drill holes. Our goal is to find patterns of model performance to guide model improvement. The satellite data is pan-arctic from 2004-2008, ice-draft data is from moored instruments in Fram Strait, the Greenland Sea and the Beaufort Sea from 1992-2008 and from submarines from 1975-2000. The drill hole data are from the Laptev and East Siberian marginal seas from 1982-1986 and from coastal stations from 1998-2009. While there are important caveats when comparing modeled results with measurements from different platforms and time periods such as these, the models agree well with moored ULS data. In general, the AOMIP models underestimate the thickness of measured ice thicker than about 2 m and overestimate thickness of ice thinner than 2 m. The simulated results are poor over the fast ice and marginal seas of the Siberian shelves. Averaging over all observational data sets, the better correlations and smaller differences from observed thickness are from the ECCO2 and UW models.

Johnson, Mark↗

Evaluating the Dominant Components of Warming in Pliocene Climate Simulations

The Pliocene Model Intercomparison Project (PlioMIP) is the first coordinated climate model comparison for a warmer palaeoclimate with atmospheric CO2 significantly higher than pre-industrial concentrations. The simulations of the mid-Pliocene warm period show global warming of between 1.8 and 3.6 C above pre-industrial surface air temperatures, with significant polar amplification. Here we perform energy balance calculations on all eight of the coupled ocean-atmosphere simulations within PlioMIP Experiment 2 to evaluate the causes of the increased temperatures and differences between the models. In the tropics simulated warming is dominated by greenhouse gas increases, with the cloud component of planetary albedo enhancing the warming in most of the models, but by widely varying amounts. The responses to mid-Pliocene climate forcing in the Northern Hemisphere midlatitudes are substantially different between the climate models, with the only consistent response being a warming due to increased greenhouse gases. In the high latitudes all the energy balance components become important, but the dominant warming influence comes from the clear sky albedo, only partially offset by the increases in the cooling impact of cloud albedo. This demonstrates the importance of specified ice sheet and high latitude vegetation boundary conditions and simulated sea ice and snow albedo feedbacks. The largest components in the overall uncertainty are associated with clouds in the tropics and polar clear sky albedo, particularly in sea ice regions. These simulations show that albedo feedbacks, particularly those of sea ice and ice sheets, provide the most significant enhancements to high latitude warming in the Pliocene.

climate models↗

Evaluation of Arctic Sea Ice Thickness Simulated by Arctic Ocean Model Intercomparison Project Models

Six Arctic Ocean Model Intercomparison Project model simulations are compared with estimates of sea ice thickness derived from pan-Arctic satellite freeboard measurements (2004-2008); airborne electromagnetic measurements (2001-2009); ice draft data from moored instruments in Fram Strait, the Greenland Sea, and the Beaufort Sea (1992-2008) and from submarines (1975-2000); and drill hole data from the Arctic basin, Laptev, and East Siberian marginal seas (1982-1986) and coastal stations (1998-2009). Despite an assessment of six models that differ in numerical methods, resolution, domain, forcing, and boundary conditions, the models generally overestimate the thickness of measured ice thinner than approximately 2 mand underestimate the thickness of ice measured thicker than about approximately 2m. In the regions of flat immobile landfast ice (shallow Siberian Seas with depths less than 25-30 m), the models generally overestimate both the total observed sea ice thickness and rates of September and October ice growth from observations by more than 4 times and more than one standard deviation, respectively. The models do not reproduce conditions of fast ice formation and growth. Instead, the modeled fast ice is replaced with pack ice which drifts, generating ridges of increasing ice thickness, in addition to thermodynamic ice growth. Considering all observational data sets, the better correlations and smaller differences from observations are from the Estimating the Circulation and Climate of the Ocean, Phase II and Pan-Arctic Ice Ocean Modeling and Assimilation System models.

sea ice↗

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗

Arctic Sea ice decay simulated for a CO2-induced temperature rise

A large scale numerical time-dependent model of sea ice that takes into account the heat fluxes in and out of the ice, the seasonal occurrence of snow, and ice motions was used in an experiment to determine the response of the Arctic Ocean ice pack to a warming of the atmosphere. The degree of warming specified is that expected for a doubling of atmospheric carbon dioxide with its associated greenhouse effect, a condition that could occur before the middle of the next century. The results of three 5-year simulations with a warmer atmosphere and varied boundary conditions were: (1) that in the face of a 5 K surface atmospheric temperature increase the ice pack disappeared completely in August and September but reformed in the central Arctic Ocean in mid fall; (2) that the simulations were moderately dependence on assumptions concerning cloud cover; and (3) that even when atmospheric temperature increases of 6-9 K were combined with an order-of-magnitude increase in the upward heat flux from the ocean, the ice still appeared in winter. It should be noted that a year-round ice-free Arctic Ocean has apparently not existed for a million years or more.

Heat flux↗

Improving the Simulation of Sea Ice Lead Conditions and Turbulent Fluxes Using RGPS Products and Merged RADARSAT, AVHRR and MODIS Data

The importance of sea ice leads in the ice-ocean-atmosphere system lies in the fact that each of the boxes in the 'surface processes' interface in this diagram is closely linked to lead conditions. For example, heat, moisture and salt exchange between the Ocean and atmosphere within the ice pack occur nearly entirely through leads. The shear, divergence and convergence associated with lead formation and closure alter surface and basal roughness and topography, which in turn affects momentum transfer in the atmosphere and ocean boundary layers, and modifies the accumulation of snow on the ice surface, which then affects heat conduction and summertime albedo. In addition to providing openings for loss of heat and moisture fluxes to the atmosphere, leads absorb solar energy, which is used to melt ice and is transmitting to the underlying ocean. Given that leads dominate the ice-ocean interface in this manner, then it stands to reason that focusing on lead treatments within models can identify performance limitations of models and yield routes for significant improvements.

Maslanik, James A.↗

Simulating Snow Over Sea Ice In Climate Models

We have evaluated two methods of simulating the seasonal cycle of snow over sea ice in and around the Arctic: The NCAR global climate model CCM3, with its standard snow hydrology, and the snow pack model SNTHERM, forced with hourly atmospheric output from CCM3. A new dataset providing dates for the onset of snow melt over Arctic sea ice provides a means for assessing basin-wide how well the models simulate melt onset, but contains no information on how long it then takes for all the snow to melt. Use of data from the SHEBA site provides very detailed information on the behavior of the snow before and during the melt season, but only for a very limited area. Russian drift data provide climatological data on the seasonal cycle of snow water equivalent and snow density, over multi-year sea ice in the central Arctic basin. These datasets are used to compare the two modeling methods, and to see if use of the more physically-realistic SNTHERM provides any significant improvements. Conclusions obtained so far include: 1. Both CCM3 and CCM3/SNTHERM do a good job overall of matching the onset of snow melt dataset; although CCM3/SNTHERM consistently trends to underestimate the date and CCM3 to overestimate it. 2. SHEBA and ice drift data for the Arctic show that CCM3/ SNTHERM does a better job than CCM3 at simulating the total melt period. 3. Ice drift snow density and accumulation data suggest that while providing superior results, CCM3/SNTHERM may still suffer from overly vigorous melting. 4. Both the large-scale atmospheric forcing and snow pack physical processes are important in proper simulation of the snow seasonal cycle. Ongoing work includes further diagnosis of CCM3/SNTHERM, use of more observational datasets, especially from marginal seas in the pan-Arctic, and full coupling of SNTHERM into CCM3 (work to date has all been off-line simulations).

Arnold, James E.↗

Upper Arctic Ocean Properties and Relationships with Sea Ice in CMIP6 Historical Simulations

While current-generation CMIP and OMIP models have clear biases in their upper Arctic Ocean hydrography, it is less clear how these biases impact the models' ability to simulate the observed Arctic sea ice mean state and trends. In this study we seek to quantify cross-relationship between sea ice and ocean states in CMIP6 historical simulations and identify common model behaviors. Multi-model mean (MMM) simulations exhibit accelerated changes in the ice and ocean system since the late 20th century. Underlying the MMM is strong inter-model variation in the simulated ice and ocean mean states and their temporal variability including trends. Despite such inter-model differences, all models show a similar ratio between sea ice reduction and upper ocean warming such that models with higher ocean warming also have higher SIE reduction and vice versa. Our results also highlight the urgent needs of reliable Arctic Ocean observations or data products in order to better contextualize modeling results.

Upper Arctic Ocean↗

Assessment of Arctic and Antarctic Sea Ice Predictability in CMIP5 Decadal Hindcasts

This paper examines the ability of coupled global climate models to predict decadal variability of Arctic and Antarctic sea ice. We analyze decadal hindcasts/predictions of 11 Coupled Model Intercomparison Project Phase 5 (CMIP5) models. Decadal hindcasts exhibit a large multimodel spread in the simulated sea ice extent, with some models deviating significantly from the observations as the predicted ice extent quickly drifts away from the initial constraint. The anomaly correlation analysis between the decadal hindcast and observed sea ice suggests that in the Arctic, for most models, the areas showing significant predictive skill become broader associated with increasing lead times. This area expansion is largely because nearly all the models are capable of predicting the observed decreasing Arctic sea ice cover. Sea ice extent in the North Pacific has better predictive skill than that in the North Atlantic (particularly at a lead time of 3-7 years), but there is a reemerging predictive skill in the North Atlantic at a lead time of 6-8 years. In contrast to the Arctic, Antarctic sea ice decadal hindcasts do not show broad predictive skill at any timescales, and there is no obvious improvement linking the areal extent of significant predictive skill to lead time increase. This might be because nearly all the models predict a retreating Antarctic sea ice cover, opposite to the observations. For the Arctic, the predictive skill of the multi-model ensemble mean outperforms most models and the persistence prediction at longer timescales, which is not the case for the Antarctic. Overall, for the Arctic, initialized decadal hindcasts show improved predictive skill compared to uninitialized simulations, although this improvement is not present in the Antarctic.

sea ice↗

High-Resolution Global Ocean Simulation in the SWOT ERA: Introducing ECCO Llc4320v2

The Surface Water and Ocean Topography (SWOT) mission is opening up new opportunities to evaluate and drive improvement in km-scale global ocean and atmosphere-ocean models via its simultaneous swath observations of sea-surface height, wind speed, and significant wave height down to 2 km scales globally. Motivated by this, we have initiated a collaborative effort at NASA Ames to run a new global ocean-sea-ice simulation in MITgcm at a nominal 1/48-deg resolution of the year 2023 that enables direct comparisons with SWOT observations. New features relative to a prior simulation at a similar resolution of the year 2012 include much more realistic tidal forcing, thermodynamic coupling in ice-shelf cavities, high-frequency hourly atmospheric forcing, much improved river runoff, double the vertical resolution in the deep ocean (173 vs 90 levels), updated bathymetry including improvements in finescale bathymetry derived from SWOT, and updated vertical mixing scheme parameters that improve the equatorial Pacific circulation.

ECCO↗

Synergistic and Collaborative Development Strategies for FV3 Powered Next Generation Unified Global Modeling System

The GFDL (Geophysical Fluid Dynamics Laboratory, NOAA - National Oceanic and Atmospheric Administration) Finite-­Volume Cubed-Sphere Dynamical Core (FV3) is a scalable and flexible dynamical core capable of both hydrostatic and non-hydrostatic atmospheric simulations. FV3 has been chosen as the dynamical core for the Next Generation Global Prediction System project (NGGPS), designed to upgrade the current NCEP (National Centers for Environmental Prediction, NOAA) operational Global Forecast System (GFS) to run as a unified, fully-coupled system in NOAA's Environmental Modeling System infrastructure. FV3 dynamic core has a long history of serving as the main engine for global atmospheric models at various Government and Academic Research Laboratories including NOAA: GFDL Climate Modeling Suite (AM4 (Atmosphere Modeling 4), CM4 (Climate Model 4) ESM4 (Earth System Model 4), Hiram (HIgh Resolution Atmospheric Model)); NASA: GMAO (Global Modeling and Assimilation Office) Goddard Earth Observing System Model (GEOS); and NCAR (National Center for Atmospheric Research) Community Earth System Model (CESM). The three primary stakeholders in FV3 (GFDL, GMAO, EMC (Environmental Modeling Center - NOAA)) have embarked on synergistic and collaborative strategies with focus on the advancement of non-hydrostratic dynamic core, physics, chemistry, and data assimilation efforts, leveraging collective strengths of each of the partnering agencies. This talk presents the ongoing plans for model advancements at GFDL, GMAO and EMC, with special emphasis on contributions towards developing community based unified global modeling system for operational and research applications at respective organizations. Future plans include extending the collaborations to the development of earth system components including GMAO's aerosol chemistry models (GOCART/MAM (Goddard Chemistry Aerosol Radiation and Transport / Modal Aerosol Model)), Land Information System (LIS), and advanced data assimilation techniques; and GFDL's Modular Ocean Model (MOM) and Sea Ice Simulator (SIS) for transition to operations at NCEP.

Putman, William↗

1985 International Geoscience and Remote Sensing Symposium (IGARSS '85), University of Massachusetts, Amherst, October 7-9, 1985, Proceedings

Papers are presented on an EM subsurface radar based on the transient field radiated by a wire antenna; the microwave dielectric, structural, and salinity properties of simulated sea ice; the extraction of sea-ice data from satellite SAR imagery, and the probing of thick vegetation canopies with a field microwave scatterometer. Also discussed are the bidirectional reflectance modeling of a conifer forest canopy, a microwave dielectric model for aggregated soils, and the estimation of soil hydraulic parameters with passive microwave data.

Carver, K. R.↗

The Impact of Sea Ice Concentration Accuracies on Climate Model Simulations with the GISS GCM

The Goddard Institute for Space Studies global climate model (GISS GCM) is used to examine the sensitivity of the simulated climate to sea ice concentration specifications in the type of simulation done in the Atmospheric Modeling Intercomparison Project (AMIP), with specified oceanic boundary conditions. Results show that sea ice concentration uncertainties of +/- 7% can affect simulated regional temperatures by more than 6 C, and biases in sea ice concentrations of +7% and -7% alter simulated annually averaged global surface air temperatures by -0.10 C and +0.17 C, respectively, over those in the control simulation. The resulting 0.27 C difference in simulated annual global surface air temperatures is reduced by a third, to 0.18 C, when considering instead biases of +4% and -4%. More broadly, least-squares fits through the temperature results of 17 simulations with ice concentration input changes ranging from increases of 50% versus the control simulation to decreases of 50% yield a yearly average global impact of 0.0107 C warming for every 1% ice concentration decrease, i.e., 1.07 C warming for the full +50% to -50% range. Regionally and on a monthly average basis, the differences can be far greater, especially in the polar regions, where wintertime contrasts between the +50% and -50% cases can exceed 30 C. However, few statistically significant effects are found outside the polar latitudes, and temperature effects over the non-polar oceans tend to be under 1 C, due in part to the specification of an unvarying annual cycle of sea surface temperatures. The +/- 7% and 14% results provide bounds on the impact (on GISS GCM simulations making use of satellite data) of satellite-derived ice concentration inaccuracies, +/- 7% being the current estimated average accuracy of satellite retrievals and +/- 4% being the anticipated improved average accuracy for upcoming satellite instruments. Results show that the impact on simulated temperatures of imposed ice concentration changes is least in summer, encouragingly the same season in which the satellite accuracies are thought to be worst. Hence the impact of satellite inaccuracies is probably less than the use of an annually averaged satellite inaccuracy would suggest.

Parkinson, Claire L.↗

Numerical Simulation of North Atlantic Sea Ice Variability, 1951 - 1980

A two-level dynamic-thermodynamic sea ice model is used to simulate the growth, drift and decay of sea ice in the Northern Hemisphere during a 30-year period, 1951 to 1980. The model is run with a daily timestep on a 222 km grid and is forced by interanually varying fields of geostrophic wind and temperature-derived thermodynamic fluxes. The objective is a quantitative description of large-scale sea ice variability in terms of the dynamic and thermodynamic processes responsible for the fluctuations, especially in the North Atlantic where sea ice represents a substantial input of fresh water. The fields of ice velocity and thickness contain strong seasonal as well as interannual variability. The mean drift pattern results in thicknesses of 4 to 5 m offshore of northern Canada and Greenland, while winter thicknesses of approximately 2 m are typical of Alaskan. Eurasian and East Greenland coastal waters. The 30-year mean fields are characterized by ecessive ice in the North Atlantic during winter and by a summer retreat that is more rapid than observed.

Walsh, J. E.↗