Sensitivity studies with a coupled ice-ocean model of the marginal ice zone

NASA Technical Reports Server (NTRS)

Roed, L. P.

1983-01-01

An analytical coupled ice-ocean model is considered which is forced by a specified wind stress acting on the open ocean as well as the ice. The analysis supports the conjecture that the upwelling dynamics at ice edges can be understood by means of a simple analytical model. In similarity with coastal problems it is shown that the ice edge upwelling is determined by the net mass flux at the boundaries of the considered region. The model is used to study the sensitivity of the upwelling dynamics in the marginal ice zone to variation in the controlling parameters. These parameters consist of combinations of the drag coefficients used in the parameterization of the stresses on the three interfaces atmosphere-ice, atmosphere-ocean, and ice-ocean. The response is shown to be sensitive to variations in these parameters in that one set of parameters may give upwelling while a slightly different set of parameters may give downwelling.

Modeling ocean wave propagation under sea ice covers

NASA Astrophysics Data System (ADS)

Zhao, Xin; Shen, Hayley H.; Cheng, Sukun

2015-02-01

Operational ocean wave models need to work globally, yet current ocean wave models can only treat ice-covered regions crudely. The purpose of this paper is to provide a brief overview of ice effects on wave propagation and different research methodology used in studying these effects. Based on its proximity to land or sea, sea ice can be classified as: landfast ice zone, shear zone, and the marginal ice zone. All ice covers attenuate wave energy. Only long swells can penetrate deep into an ice cover. Being closest to open water, wave propagation in the marginal ice zone is the most complex to model. The physical appearance of sea ice in the marginal ice zone varies. Grease ice, pancake ice, brash ice, floe aggregates, and continuous ice sheet may be found in this zone at different times and locations. These types of ice are formed under different thermal-mechanical forcing. There are three classic models that describe wave propagation through an idealized ice cover: mass loading, thin elastic plate, and viscous layer models. From physical arguments we may conjecture that mass loading model is suitable for disjoint aggregates of ice floes much smaller than the wavelength, thin elastic plate model is suitable for a continuous ice sheet, and the viscous layer model is suitable for grease ice. For different sea ice types we may need different wave ice interaction models. A recently proposed viscoelastic model is able to synthesize all three classic models into one. Under suitable limiting conditions it converges to the three previous models. The complete theoretical framework for evaluating wave propagation through various ice covers need to be implemented in the operational ocean wave models. In this review, we introduce the sea ice types, previous wave ice interaction models, wave attenuation mechanisms, the methods to calculate wave reflection and transmission between different ice covers, and the effect of ice floe breaking on shaping the sea ice morphology

Modeling seasonality of ice and ocean carbon production in the Arctic

NASA Astrophysics Data System (ADS)

Jin, M.; Deal, C. M.; Ji, R.

2011-12-01

In the Arctic Ocean, both phytoplankton and sea ice algae are important contributors to the primary production and the arctic food web. Copepod in the arctic regions have developed their feeding habit depending on the timing between the ice algal bloom and the subsequent phytoplankton bloom. A mismatch of the timing due to climate changes could have dramatic consequences on the food web as shown by some regional observations. In this study, a global coupled ice-ocean-ecosystem model was used to assess the seasonality of the ice algal and phytoplankton blooms in the arctic. The ice-ocean ecosystem modules are fully coupled in the physical model POP-CICE (Parallel Ocean Program- Los Alamos Sea Ice Model). The model results are compared with various observations. The modeled ice and ocean carbon production were analyzed by regions and their linkage to the physical environment changes (such as changes of ice concentration and water temperature, and light intensity etc.) between low- and high-ice years.

Slush Fund: Modeling the Multiphase Physics of Oceanic Ices

NASA Astrophysics Data System (ADS)

Buffo, J.; Schmidt, B. E.

2016-12-01

The prevalence of ice interacting with an ocean, both on Earth and throughout the solar system, and its crucial role as the mediator of exchange between the hydrosphere below and atmosphere above, have made quantifying the thermodynamic, chemical, and physical properties of the ice highly desirable. While direct observations of these quantities exist, their scarcity increases with the difficulty of obtainment; the basal surfaces of terrestrial ice shelves remain largely unexplored and the icy interiors of moons like Europa and Enceladus have never been directly observed. Our understanding of these entities thus relies on numerical simulation, and the efficacy of their incorporation into larger systems models is dependent on the accuracy of these initial simulations. One characteristic of seawater, likely shared by the oceans of icy moons, is that it is a solution. As such, when it is frozen a majority of the solute is rejected from the forming ice, concentrating in interstitial pockets and channels, producing a two-component reactive porous media known as a mushy layer. The multiphase nature of this layer affects the evolution and dynamics of the overlying ice mass. Additionally ice can form in the water column and accrete onto the basal surface of these ice masses via buoyancy driven sedimentation as frazil or platelet ice. Numerical models hoping to accurately represent ice-ocean interactions should include the multiphase behavior of these two phenomena. While models of sea ice have begun to incorporate multiphase physics into their capabilities, no models of ice shelves/shells explicitly account for the two-phase behavior of the ice-ocean interface. Here we present a 1D multiphase model of floating oceanic ice that includes parameterizations of both density driven advection within the `mushy layer' and buoyancy driven sedimentation. The model is validated against contemporary sea ice models and observational data. Environmental stresses such as supercooling and

Decadal-Scale Response of the Antarctic Ice sheet to a Warming Ocean using the POPSICLES Coupled Ice Sheet-Ocean model

NASA Astrophysics Data System (ADS)

Martin, D. F.; Asay-Davis, X.; Cornford, S. L.; Price, S. F.; Ng, E. G.; Collins, W.

2015-12-01

We present POPSICLES simulation results covering the full Antarctic Ice Sheet and the Southern Ocean spanning the period from 1990 to 2010. We use the CORE v. 2 interannual forcing data to force the ocean model. Simulations are performed at 0.1o(~5 km) ocean resolution with adaptive ice sheet resolution as fine as 500 m to adequately resolve the grounding line dynamics. We discuss the effect of improved ocean mixing and subshelf bathymetry (vs. the standard Bedmap2 bathymetry) on the behavior of the coupled system, comparing time-averaged melt rates below a number of major ice shelves with those reported in the literature. We also present seasonal variability and decadal melting trends from several Antarctic regions, along with the response of the ice shelves and the consequent dynamic response of the grounded ice sheet.POPSICLES couples the POP2x ocean model, a modified version of the Parallel Ocean Program, and the BISICLES ice-sheet model. POP2x includes sub-ice-shelf circulation using partial top cells and the commonly used three-equation boundary layer physics. Standalone POP2x output compares well with standard ice-ocean test cases (e.g., ISOMIP) and other continental-scale simulations and melt-rate observations. BISICLES makes use of adaptive mesh refinement and a 1st-order accurate momentum balance similar to the L1L2 model of Schoof and Hindmarsh to accurately model regions of dynamic complexity, such as ice streams, outlet glaciers, and grounding lines. Results of BISICLES simulations have compared favorably to comparable simulations with a Stokes momentum balance in both idealized tests (MISMIP-3d) and realistic configurations.The figure shows the BISICLES-computed vertically-integrated grounded ice velocity field 5 years into a 20-year coupled full-continent Antarctic-Southern-Ocean simulation. Submarine melt rates are painted onto the surface of the floating ice shelves. Grounding lines are shown in green.

Planetary Ice-Oceans: Numerical Modeling Study of Ice-Shell Growth in Convecting Two-Phase Systems

NASA Astrophysics Data System (ADS)

Allu Peddinti, Divya; McNamara, Allen

2017-04-01

Several icy bodies in the Solar system such as the icy moons Europa and Enceladus exhibit signs of subsurface oceans underneath an ice-shell. For Europa, the geologically young surface, the presence of surface features and the aligned surface chemistry pose interesting questions about formation of the ice-shell and its interaction with the ocean below. This also ties in with its astrobiological potential and implications for similar ice-ocean systems elsewhere in the cosmos. The overall thickness of the H2O layer on Europa is estimated to be 100-150 km while the thickness of the ice-shell is debated. Additionally, Europa is subject to tidal heating due to interaction with Jupiter's immense gravity field. It is of interest to understand how the ice-shell thickness varies in the presence of tidal internal heating and the localization of heating in different regions of the ice-shell. Thus this study aims to determine the effect of tidal internal heating on the growth rate of the ice-shell over time. We perform geodynamic modeling of the ice-ocean system in order to understand how the ice-shell thickness changes with time. The convection code employs the ice Ih-water phase diagram in order to model the two-phase convecting ice-ocean system. All the models begin from an initial warm thick ocean that cools from the top. The numerical experiments analyze three cases: case 1 with no tidal internal heating in the system, case 2 with constant tidal internal heating in the ice and case 3 with viscosity-dependent tidal internal heating in the ice. We track the ice-shell thickness as a function of time as the system cools. Modeling results so far have identified that the shell growth rate changes substantially at a point in time that coincides with a change in the planform of ice-convection cells. Additionally, the velocity vs depth plots indicate a shift from a conduction dominant to a convection dominant ice regime. We compare the three different cases to provide a

An ice-ocean coupled model for the Northern Hemisphere

NASA Technical Reports Server (NTRS)

Cheng, Abe; Preller, Ruth

1992-01-01

The Hibler ice model has been modified and adapted to a domain that includes most of the sea ice-covered areas in the Northern Hemisphere. This model, joined with the Cox ocean model, is developed as an enhancement to the U.S. Navy's sea ice forecasting, PIPS, and is termed PIPS2.0. Generally, the modeled ice edge is consistent with the Navy-NOAA Joint Ice Center weekly analysis, and the modeled ice thickness distribution agrees with submarine sonar data in the central Arctic basin.

One hundred years of Arctic ice cover variations as simulated by a one-dimensional, ice-ocean model

NASA Astrophysics Data System (ADS)

Hakkinen, S.; Mellor, G. L.

1990-09-01

A one-dimensional ice-ocean model consisting of a second moment, turbulent closure, mixed layer model and a three-layer snow-ice model has been applied to the simulation of Arctic ice mass and mixed layer properties. The results for the climatological seasonal cycle are discussed first and include the salt and heat balance in the upper ocean. The coupled model is then applied to the period 1880-1985, using the surface air temperature fluctuations from Hansen et al. (1983) and from Wigley et al. (1981). The analysis of the simulated large variations of the Arctic ice mass during this period (with similar changes in the mixed layer salinity) shows that the variability in the summer melt determines to a high degree the variability in the average ice thickness. The annual oceanic heat flux from the deep ocean and the maximum freezing rate and associated nearly constant minimum surface salinity flux did not vary significantly interannually. This also implies that the oceanic influence on the Arctic ice mass is minimal for the range of atmospheric variability tested.

Calibration of sea ice dynamic parameters in an ocean-sea ice model using an ensemble Kalman filter

NASA Astrophysics Data System (ADS)

Massonnet, F.; Goosse, H.; Fichefet, T.; Counillon, F.

2014-07-01

The choice of parameter values is crucial in the course of sea ice model development, since parameters largely affect the modeled mean sea ice state. Manual tuning of parameters will soon become impractical, as sea ice models will likely include more parameters to calibrate, leading to an exponential increase of the number of possible combinations to test. Objective and automatic methods for parameter calibration are thus progressively called on to replace the traditional heuristic, "trial-and-error" recipes. Here a method for calibration of parameters based on the ensemble Kalman filter is implemented, tested and validated in the ocean-sea ice model NEMO-LIM3. Three dynamic parameters are calibrated: the ice strength parameter P*, the ocean-sea ice drag parameter Cw, and the atmosphere-sea ice drag parameter Ca. In twin, perfect-model experiments, the default parameter values are retrieved within 1 year of simulation. Using 2007-2012 real sea ice drift data, the calibration of the ice strength parameter P* and the oceanic drag parameter Cw improves clearly the Arctic sea ice drift properties. It is found that the estimation of the atmospheric drag Ca is not necessary if P* and Cw are already estimated. The large reduction in the sea ice speed bias with calibrated parameters comes with a slight overestimation of the winter sea ice areal export through Fram Strait and a slight improvement in the sea ice thickness distribution. Overall, the estimation of parameters with the ensemble Kalman filter represents an encouraging alternative to manual tuning for ocean-sea ice models.

Observed and Modeled Trends in Southern Ocean Sea Ice

NASA Technical Reports Server (NTRS)

Parkinson, Claire L.

2003-01-01

Conceptual models and global climate model (GCM) simulations have both indicated the likelihood of an enhanced sensitivity to climate change in the polar regions, derived from the positive feedbacks brought about by the polar abundance of snow and ice surfaces. Some models further indicate that the changes in the polar regions can have a significant impact globally. For instance, 37% of the temperature sensitivity to a doubling of atmospheric CO2 in simulations with the GCM of the Goddard Institute for Space Studies (GISS) is attributable exclusively to inclusion of sea ice variations in the model calculations. Both sea ice thickness and sea ice extent decrease markedly in the doubled CO, case, thereby allowing the ice feedbacks to occur. Stand-alone sea ice models have shown Southern Ocean hemispherically averaged winter ice-edge retreats of 1.4 deg latitude for each 1 K increase in atmospheric temperatures. Observations, however, show a much more varied Southern Ocean ice cover, both spatially and temporally, than many of the modeled expectations. In fact, the satellite passive-microwave record of Southern Ocean sea ice since late 1978 has revealed overall increases rather than decreases in ice extents, with ice extent trends on the order of 11,000 sq km/year. When broken down spatially, the positive trends are strongest in the Ross Sea, while the trends are negative in the Bellingshausen/Amundsen Seas. Greater spatial detail can be obtained by examining trends in the length of the sea ice season, and those trends show a coherent picture of shortening sea ice seasons throughout almost the entire Bellingshausen and Amundsen Seas to the west of the Antarctic Peninsula and in the far western Weddell Sea immediately to the east of the Peninsula, with lengthening sea ice seasons around much of the rest of the continent. This pattern corresponds well with the spatial pattern of temperature trends, as the Peninsula region is the one region in the Antarctic with a strong

Geodynamic Modeling of Planetary Ice-Oceans: Evolution of Ice-Shell Thickness in Convecting Two-Phase Systems

NASA Astrophysics Data System (ADS)

Allu Peddinti, D.; McNamara, A. K.

2016-12-01

Along with the newly unveiled icy surface of Pluto, several icy planetary bodies show indications of an active surface perhaps underlain by liquid oceans of some size. This augments the interest to explore the evolution of an ice-ocean system and its surface implications. The geologically young surface of the Jovian moon Europa lends much speculation to variations in ice-shell thickness over time. Along with the observed surface features, it suggests the possibility of episodic convection and conduction within the ice-shell as it evolved. What factors would control the growth of the ice-shell as it forms? If and how would those factors determine the thickness of the ice-shell and consequently the heat transfer? Would parameters such as tidal heating or initial temperature affect how the ice-shell grows and to what significance? We perform numerical experiments using geodynamical models of the two-phase ice-water system to study the evolution of planetary ice-oceans such as that of Europa. The models evolve self-consistently from an initial liquid ocean as it cools with time. The effects of presence, absence and magnitude of tidal heating on ice-shell thickness are studied in different models. The vigor of convection changes as the ice-shell continues to thicken. Initial modeling results track changes in the growth rate of the ice-shell as the vigor of the convection changes. The magnitude and temporal location of the rate change varies with different properties of tidal heating and values of initial temperature. A comparative study of models is presented to demonstrate how as the ice-shell is forming, its growth rate and convection are affected by processes such as tidal heating.

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

NASA Technical Reports Server (NTRS)

Johnson, Mark; Proshuntinsky, Andrew; Aksenov, Yevgeny; Nguyen, An T.; Lindsay, Ron; Haas, Christian; Zhang, Jinlun; Diansky, Nikolay; Kwok, Ron; Maslowski, Wieslaw;

Modelling of Sea Ice Thermodynamics and Biogeochemistry during the N-ICE2015 Expedition in the Arctic Ocean

NASA Astrophysics Data System (ADS)

Meyer, A.; Duarte, P.; Mork Olsen, L.; Kauko, H.; Assmy, P.; Rösel, A.; Itkin, P.; Hudson, S. R.; Granskog, M. A.; Gerland, S.; Sundfjord, A.; Steen, H.; Jeffery, N.; Hunke, E. C.; Elliott, S.; Turner, A. K.

2016-12-01

Changes in the sea ice regime of the Arctic Ocean over the last decades from a thick perennial multiyear ice to a first year ice have been well documented. These changes in the sea ice regime will affect feedback mechanisms between the sea ice, atmosphere and ocean. Here we evaluate the performance of the Los Alamos Sea Ice Model (CICE), a state of the art sea ice model, to predict sea ice physical and biogeochemical properties at time scales of a few weeks. We also identify the most problematic prognostic variables and what is necessary to improve their forecast. The availability of a complete data set of forcing collected during the Norwegian Young sea Ice (N-ICE-2015) expedition north of Svalbard opens the possibility to properly test CICE. Oceanographic, atmospheric, sea ice, snow, and biological data were collected above, on, and below the ice using R/V Lance as the base for the ice camps that were drifting south towards the Fram Strait. Over six months, four different drifts took place, from the Nansen Basin, through the marginal ice zone, to the open ocean. Obtained results from the model show a good performance regarding ice thickness, salinity and temperature. Nutrients and sea ice algae are however not modelled as accurately. We hypothesize that improvements in biogeochemical modeling may be achieved by complementing brine drainage with a diffusion parameterization and biogeochemical modeling with the introduction of an explicit formulation to forecast chlorophyll and regulate photosynthetic efficiency.

Wave-Ice interaction in the Marginal Ice Zone: Toward a Wave-Ocean-Ice Coupled Modeling System

DTIC Science & Technology

2015-09-30

MIZ using WW3 (3 frequency bins, ice retreat in August and ice advance in October); Blue (solid): Based on observations near Antarctica by Meylan...1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Wave- Ice interaction in the Marginal Ice Zone: Toward a...Wave-Ocean- Ice Coupled Modeling System W. E. Rogers Naval Research Laboratory, Code 7322 Stennis Space Center, MS 39529 phone: (228) 688-4727

Investigation of land ice-ocean interaction with a fully coupled ice-ocean model: 2. Sensitivity to external forcings

NASA Astrophysics Data System (ADS)

Goldberg, D. N.; Little, C. M.; Sergienko, O. V.; Gnanadesikan, A.; Hallberg, R.; Oppenheimer, M.

2012-06-01

A coupled ice stream-ice shelf-ocean cavity model is used to assess the sensitivity of the coupled system to far-field ocean temperatures, varying from 0.0 to 1.8°C, as well as sensitivity to the parameters controlling grounded ice flow. A response to warming is seen in grounding line retreat and grounded ice loss that cannot be inferred from the response of integrated melt rates alone. This is due to concentrated thinning at the ice shelf lateral margin, and to processes that contribute to this thinning. Parameters controlling the flow of grounded ice have a strong influence on the response to sub-ice shelf melting, but this influence is not seen until several years after an initial perturbation in temperatures. The simulated melt rates are on the order of that observed for Pine Island Glacier in the 1990s. However, retreat rates are much slower, possibly due to unrepresented bedrock features.

Adaptation of an unstructured-mesh, finite-element ocean model to the simulation of ocean circulation beneath ice shelves

NASA Astrophysics Data System (ADS)

Kimura, Satoshi; Candy, Adam S.; Holland, Paul R.; Piggott, Matthew D.; Jenkins, Adrian

2013-07-01

Several different classes of ocean model are capable of representing floating glacial ice shelves. We describe the incorporation of ice shelves into Fluidity-ICOM, a nonhydrostatic finite-element ocean model with the capacity to utilize meshes that are unstructured and adaptive in three dimensions. This geometric flexibility offers several advantages over previous approaches. The model represents melting and freezing on all ice-shelf surfaces including vertical faces, treats the ice shelf topography as continuous rather than stepped, and does not require any smoothing of the ice topography or any of the additional parameterisations of the ocean mixed layer used in isopycnal or z-coordinate models. The model can also represent a water column that decreases to zero thickness at the 'grounding line', where the floating ice shelf is joined to its tributary ice streams. The model is applied to idealised ice-shelf geometries in order to demonstrate these capabilities. In these simple experiments, arbitrarily coarsening the mesh outside the ice-shelf cavity has little effect on the ice-shelf melt rate, while the mesh resolution within the cavity is found to be highly influential. Smoothing the vertical ice front results in faster flow along the smoothed ice front, allowing greater exchange with the ocean than in simulations with a realistic ice front. A vanishing water-column thickness at the grounding line has little effect in the simulations studied. We also investigate the response of ice shelf basal melting to variations in deep water temperature in the presence of salt stratification.

Representing grounding line migration in synchronous coupling between a marine ice sheet model and a z-coordinate ocean model

NASA Astrophysics Data System (ADS)

Goldberg, D. N.; Snow, K.; Holland, P.; Jordan, J. R.; Campin, J.-M.; Heimbach, P.; Arthern, R.; Jenkins, A.

2018-05-01

Synchronous coupling is developed between an ice sheet model and a z-coordinate ocean model (the MITgcm). A previously-developed scheme to allow continuous vertical movement of the ice-ocean interface of a floating ice shelf ("vertical coupling") is built upon to allow continuous movement of the grounding line, or point of floatation of the ice sheet ("horizontal coupling"). Horizontal coupling is implemented through the maintenance of a thin layer of ocean ( ∼ 1 m) under grounded ice, which is inflated into the real ocean as the ice ungrounds. This is accomplished through a modification of the ocean model's nonlinear free surface evolution in a manner akin to a hydrological model in the presence of steep bathymetry. The coupled model is applied to a number of idealized geometries and shown to successfully represent ocean-forced marine ice sheet retreat while maintaining a continuous ocean circulation.

Submesoscale Sea Ice-Ocean Interactions in Marginal Ice Zones

NASA Astrophysics Data System (ADS)

Manucharyan, Georgy E.; Thompson, Andrew F.

2017-12-01

Signatures of ocean eddies, fronts, and filaments are commonly observed within marginal ice zones (MIZs) from satellite images of sea ice concentration, and in situ observations via ice-tethered profilers or underice gliders. However, localized and intermittent sea ice heating and advection by ocean eddies are currently not accounted for in climate models and may contribute to their biases and errors in sea ice forecasts. Here, we explore mechanical sea ice interactions with underlying submesoscale ocean turbulence. We demonstrate that the release of potential energy stored in meltwater fronts can lead to energetic submesoscale motions along MIZs with spatial scales O(10 km) and Rossby numbers O(1). In low-wind conditions, cyclonic eddies and filaments efficiently trap the sea ice and advect it over warmer surface ocean waters where it can effectively melt. The horizontal eddy diffusivity of sea ice mass and heat across the MIZ can reach O(200 m2 s-1). Submesoscale ocean variability also induces large vertical velocities (order 10 m d-1) that can bring relatively warm subsurface waters into the mixed layer. The ocean-sea ice heat fluxes are localized over cyclonic eddies and filaments reaching about 100 W m-2. We speculate that these submesoscale-driven intermittent fluxes of heat and sea ice can contribute to the seasonal evolution of MIZs. With the continuing global warming and sea ice thickness reduction in the Arctic Ocean, submesoscale sea ice-ocean processes are expected to become increasingly prominent.

Results from ISOMIP+ and MISOMIP1, two interrelated marine ice sheet and ocean model intercomparison projects

NASA Astrophysics Data System (ADS)

Asay-Davis, X.; Galton-Fenzi, B.; Gwyther, D.; Jourdain, N.; Martin, D. F.; Nakayama, Y.; Seroussi, H. L.

2016-12-01

MISMIP+ (the third Marine Ice Sheet MIP), ISOMIP+ (the second Ice Shelf-Ocean MIP) and MISOMIP1 (the first Marine Ice Sheet-Ocean MIP) prescribe a set of idealized experiments for marine ice-sheet models, ocean models with ice-shelf cavities, and coupled ice sheet-ocean models, respectively. Here, we present results from ISOMIP+ and MISOMIP1 experiments using several ocean-only and coupled ice sheet-ocean models. Among the ocean models, we show that differences in model behavior are significant enough that similar results can only be achieved by tuning model parameters (the heat- and salt-transfer coefficients across the sub-ice-shelf boundary layer) for each model. This tuning is constrained by a desired mean melt rate in quasi-steady state under specified forcing conditions, akin to tuning the models to match observed melt rates. We compare the evolution of ocean temperature transects, melt rate, friction velocity and thermal driving between ocean models for the five ISOMIP+ experiments (Ocean0-4), which have prescribed ice-shelf topography. We find that melt patterns differ between models based on the relative importance of overturning strength and vertical mixing of temperature even when the models have been tuned to achieve similar melt rates near the grounding line. For the two MISOMIP1 experiments (IceOcean1 without dynamic calving and IceOcean2 with a simple calving parameterization), we compare temperature transects, melt rate, ice-shelf topography and grounded area across models and for several model configurations. Consistent with preliminary results from MISMIP+, we find that for a given coupled model, the use of a Coulomb-limited basal friction parameterization below grounded ice and the application of dynamic calving both significantly increase the rate of grounding-line retreat, whereas the rate of retreat appears to be less sensitive to the ice stress approximation (shallow-shelf approximation, higher-order, etc.). We show that models with similar

Submesoscale sea ice-ocean interactions in marginal ice zones

NASA Astrophysics Data System (ADS)

Thompson, A. F.; Manucharyan, G.

2017-12-01

Signatures of ocean eddies, fronts and filaments are commonly observed within the marginal ice zones (MIZ) from satellite images of sea ice concentration, in situ observations via ice-tethered profilers or under-ice gliders. Localized and intermittent sea ice heating and advection by ocean eddies are currently not accounted for in climate models and may contribute to their biases and errors in sea ice forecasts. Here, we explore mechanical sea ice interactions with underlying submesoscale ocean turbulence via a suite of numerical simulations. We demonstrate that the release of potential energy stored in meltwater fronts can lead to energetic submesoscale motions along MIZs with sizes O(10 km) and Rossby numbers O(1). In low-wind conditions, cyclonic eddies and filaments efficiently trap the sea ice and advect it over warmer surface ocean waters where it can effectively melt. The horizontal eddy diffusivity of sea ice mass and heat across the MIZ can reach O(200 m2 s-1). Submesoscale ocean variability also induces large vertical velocities (order of 10 m day-1) that can bring relatively warm subsurface waters into the mixed layer. The ocean-sea ice heat fluxes are localized over cyclonic eddies and filaments reaching about 100 W m-2. We speculate that these submesoscale-driven intermittent fluxes of heat and sea ice can potentially contribute to the seasonal evolution of MIZs. With continuing global warming and sea ice thickness reduction in the Arctic Ocean, as well as the large expanse of thin sea ice in the Southern Ocean, submesoscale sea ice-ocean processes are expected to play a significant role in the climate system.

Improved Upper Ocean/Sea Ice Modeling in the GISS GCM for Investigating Climate Change

NASA Technical Reports Server (NTRS)

1997-01-01

This project built on our previous results in which we highlighted the importance of sea ice in overall climate sensitivity by determining that for both warming and cooling climates, when sea ice was not allowed to change, climate sensitivity was reduced by 35-40%. We also modified the Goddard Institute for Space Studies (GISS) 8 deg x lO deg atmospheric General Circulation Model (GCM) to include an upper-ocean/sea-ice model involving the Semtner three-layer ice/snow thermodynamic model, the Price et al. (1986) ocean mixed layer model and a general upper ocean vertical advection/diffusion scheme for maintaining and fluxing properties across the pycnocline. This effort, in addition to improving the sea ice representation in the AGCM, revealed a number of sensitive components of the sea ice/ocean system. For example, the ability to flux heat through the ice/snow properly is critical in order to resolve the surface temperature properly, since small errors in this lead to unrestrained climate drift. The present project, summarized in this report, had as its objectives: (1) introducing a series of sea ice and ocean improvements aimed at overcoming remaining weaknesses in the GCM sea ice/ocean representation, and (2) performing a series of sensitivity experiments designed to evaluate the climate sensitivity of the revised model to both Antarctic and Arctic sea ice, determine the sensitivity of the climate response to initial ice distribution, and investigate the transient response to doubling CO2.

Processes driving sea ice variability in the Bering Sea in an eddying ocean/sea ice model: Mean seasonal cycle

NASA Astrophysics Data System (ADS)

Li, Linghan; McClean, Julie L.; Miller, Arthur J.; Eisenman, Ian; Hendershott, Myrl C.; Papadopoulos, Caroline A.

2014-12-01

The seasonal cycle of sea ice variability in the Bering Sea, together with the thermodynamic and dynamic processes that control it, are examined in a fine resolution (1/10°) global coupled ocean/sea-ice model configured in the Community Earth System Model (CESM) framework. The ocean/sea-ice model consists of the Los Alamos National Laboratory Parallel Ocean Program (POP) and the Los Alamos Sea Ice Model (CICE). The model was forced with time-varying reanalysis atmospheric forcing for the time period 1970-1989. This study focuses on the time period 1980-1989. The simulated seasonal-mean fields of sea ice concentration strongly resemble satellite-derived observations, as quantified by root-mean-square errors and pattern correlation coefficients. The sea ice energy budget reveals that the seasonal thermodynamic ice volume changes are dominated by the surface energy flux between the atmosphere and the ice in the northern region and by heat flux from the ocean to the ice along the southern ice edge, especially on the western side. The sea ice force balance analysis shows that sea ice motion is largely associated with wind stress. The force due to divergence of the internal ice stress tensor is large near the land boundaries in the north, and it is small in the central and southern ice-covered region. During winter, which dominates the annual mean, it is found that the simulated sea ice was mainly formed in the northern Bering Sea, with the maximum ice growth rate occurring along the coast due to cold air from northerly winds and ice motion away from the coast. South of St Lawrence Island, winds drive the model sea ice southwestward from the north to the southwestern part of the ice-covered region. Along the ice edge in the western Bering Sea, model sea ice is melted by warm ocean water, which is carried by the simulated Bering Slope Current flowing to the northwest, resulting in the S-shaped asymmetric ice edge. In spring and fall, similar thermodynamic and dynamic

A coupled dynamic-thermodynamic model of an ice-ocean system in the marginal ice zone

NASA Technical Reports Server (NTRS)

Hakkinen, Sirpa

1987-01-01

Thermodynamics are incorporated into a coupled ice-ocean model in order to investigate wind-driven ice-ocean processes in the marginal zone. Upswelling at the ice edge which is generated by the difference in the ice-air and air-water surface stresses is found to give rise to a strong entrainment by drawing the pycnocline closer to the surface. Entrainment is shown to be negligible outside the areas affected by the ice edge upswelling. If cooling at the top is included in the model, the heat and salt exchanges are further enhanced in the upswelling areas. It is noted that new ice formation occurs in the region not affected by ice edge upswelling, and it is suggested that the high-salinity mixed layer regions (with a scale of a few Rossby radii of deformation) will overturn due to cooling, possibly contributing to the formation of deep water.

Sensitivity of a climatologically-driven sea ice model to the ocean heat flux

NASA Technical Reports Server (NTRS)

Parkinson, C. L.; Good, M. R.

1982-01-01

Ocean heat flux sensitivity was studied on a numerical model of sea ice covering the Weddell Sea region of the southern ocean. The model is driven by mean monthly climatological atmospheric variables. For each model run, the ocean heat flux is uniform in both space and time. Ocean heat fluxes below 20 W m to the minus 2 power do not provide sufficient energy to allow the ice to melt to its summertime thicknesses and concentrations by the end of the 14 month simulation, whereas ocean heat fluxes of 30 W m to the minus 2 power and above result in too much ice melt, producing the almost total disappearance of ice in the Weddell Sea by the end of the 14 months. These results are dependent on the atmospheric forcing fields.

Ocean-ice interaction in the marginal ice zone

NASA Technical Reports Server (NTRS)

Liu, Antony K.; Peng, Chich Y.

1994-01-01

Ocean ice interaction processes in the Marginal Ice Zone (MIZ) by wind, waves, and mesoscale features, such as upwelling and eddies, are studied using ERS-1 Synthetic Aperture Radar (SAR) images and ocean ice interaction model. A sequence of SAR images of the Chukchi Sea MIZ with three days interval are studied for ice edge advance/retreat. Simultaneous current measurements from the northeast Chukchi Sea as well as the Barrow wind record are used to interpret the MIZ dynamics.

Improved Upper Ocean/Sea Ice Modeling in the GISS GCM for Investigating Climate Change

NASA Technical Reports Server (NTRS)

1998-01-01

This project built on our previous results in which we highlighted the importance of sea ice in overall climate sensitivity by determining that for both warming and cooling climates, when sea ice was not allowed to change, climate sensitivity was reduced by 35-40%. We also modified the GISS 8 deg x lO deg atmospheric GCM to include an upper-ocean/sea-ice model involving the Semtner three-layer ice/snow thermodynamic model, the Price et al. (1986) ocean mixed layer model and a general upper ocean vertical advection/diffusion scheme for maintaining and fluxing properties across the pycnocline. This effort, in addition to improving the sea ice representation in the AGCM, revealed a number of sensitive components of the sea ice/ocean system. For example, the ability to flux heat through the ice/snow properly is critical in order to resolve the surface temperature properly, since small errors in this lead to unrestrained climate drift. The present project, summarized in this report, had as its objectives: (1) introducing a series of sea ice and ocean improvements aimed at overcoming remaining weaknesses in the GCM sea ice/ocean representation, and (2) performing a series of sensitivity experiments designed to evaluate the climate sensitivity of the revised model to both Antarctic and Arctic sea ice, determine the sensitivity of the climate response to initial ice distribution, and investigate the transient response to doubling CO2.

Effects of Model Resolution and Ocean Mixing on Forced Ice-Ocean Physical and Biogeochemical Simulations Using Global and Regional System Models

NASA Astrophysics Data System (ADS)

Jin, Meibing; Deal, Clara; Maslowski, Wieslaw; Matrai, Patricia; Roberts, Andrew; Osinski, Robert; Lee, Younjoo J.; Frants, Marina; Elliott, Scott; Jeffery, Nicole; Hunke, Elizabeth; Wang, Shanlin

2018-01-01

The current coarse-resolution global Community Earth System Model (CESM) can reproduce major and large-scale patterns but is still missing some key biogeochemical features in the Arctic Ocean, e.g., low surface nutrients in the Canada Basin. We incorporated the CESM Version 1 ocean biogeochemical code into the Regional Arctic System Model (RASM) and coupled it with a sea-ice algal module to investigate model limitations. Four ice-ocean hindcast cases are compared with various observations: two in a global 1° (40˜60 km in the Arctic) grid: G1deg and G1deg-OLD with/without new sea-ice processes incorporated; two on RASM's 1/12° (˜9 km) grid R9km and R9km-NB with/without a subgrid scale brine rejection parameterization which improves ocean vertical mixing under sea ice. Higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice extent, ice thickness, and ice algae. In the Bering Sea shelf, only higher resolution contributed to lower model errors in salinity, nitrate (NO3), and chlorophyll-a (Chl-a). In the Arctic Basin, model errors in mixed layer depth (MLD) were reduced 36% by brine rejection parameterization, 20% by new sea-ice processes, and 6% by higher resolution. The NO3 concentration biases were caused by both MLD bias and coarse resolution, because of excessive horizontal mixing of high NO3 from the Chukchi Sea into the Canada Basin in coarse resolution models. R9km showed improvements over G1deg on NO3, but not on Chl-a, likely due to light limitation under snow and ice cover in the Arctic Basin.

Modeling South Pacific Ice-Ocean Interactions in the Global Climate System

NASA Technical Reports Server (NTRS)

Holland, David M.; Jenkins, Adrian; Jacobs, Stanley S.

2001-01-01

The objective of this project has been to improve the modeling of interactions between large Antarctic ice shelves and adjacent regions of the Southern Ocean. Our larger goal is to gain a better understanding of the extent to which the ocean controls ice shelf attrition, thereby influencing the size and dynamics of the Antarctic Ice Sheet. Melting and freezing under ice shelves also impacts seawater properties, regional upwelling and sinking and the larger-scale ocean circulation. Modifying an isopycnal coordinate general circulation model for use in sub-ice shelf cavities, we found that the abrupt change in water column thickness at an ice shelf front does not form a strong barrier to buoyancy-driven circulation across the front. Outflow along the ice shelf base, driven by melting of the thickest ice, is balanced by deep inflow. Substantial effort was focused on the Filchner-Ronne cavity, where other models have been applied and time-series records are available from instruments suspended beneath the ice. A model comparison indicated that observed changes in the production of High Salinity Shelf Water could have a major impact on circulation within the cavity. This water propagates into the cavity with an asymmetric seasonal signal that has similar phasing and shape in the model and observations, and can be related to winter production at the sea surface. Even remote parts of the sub-ice shelf cavity are impacted by external forcing on sub-annual time scales. This shows that cavity circulations and products, and therefore cavity shape, will respond to interannual variability in sea ice production and longer-term climate change. The isopycnal model gives generally lower net melt rates than have been obtained from other models and oceanographic data, perhaps due to its boundary layer formulation, or the lack of tidal forcing. Work continues on a manuscript describing the Ross cavity results.

Design of the MISMIP+, ISOMIP+, and MISOMIP ice-sheet, ocean, and coupled ice sheet-ocean intercomparison projects

NASA Astrophysics Data System (ADS)

Asay-Davis, Xylar; Cornford, Stephen; Martin, Daniel; Gudmundsson, Hilmar; Holland, David; Holland, Denise

2015-04-01

The MISMIP and MISMIP3D marine ice sheet model intercomparison exercises have become popular benchmarks, and several modeling groups have used them to show how their models compare to both analytical results and other models. Similarly, the ISOMIP (Ice Shelf-Ocean Model Intercomparison Project) experiments have acted as a proving ground for ocean models with sub-ice-shelf cavities.As coupled ice sheet-ocean models become available, an updated set of benchmark experiments is needed. To this end, we propose sequel experiments, MISMIP+ and ISOMIP+, with an end goal of coupling the two in a third intercomparison exercise, MISOMIP (the Marine Ice Sheet-Ocean Model Intercomparison Project). Like MISMIP3D, the MISMIP+ experiments take place in an idealized, three-dimensional setting and compare full 3D (Stokes) and reduced, hydrostatic models. Unlike the earlier exercises, the primary focus will be the response of models to sub-shelf melting. The chosen configuration features an ice shelf that experiences substantial lateral shear and buttresses the upstream ice, and so is well suited to melting experiments. Differences between the steady states of each model are minor compared to the response to melt-rate perturbations, reflecting typical real-world applications where parameters are chosen so that the initial states of all models tend to match observations. The three ISOMIP+ experiments have been designed to to make use of the same bedrock topography as MISMIP+ and using ice-shelf geometries from MISMIP+ results produced by the BISICLES ice-sheet model. The first two experiments use static ice-shelf geometries to simulate the evolution of ocean dynamics and resulting melt rates to a quasi-steady state when far-field forcing changes in either from cold to warm or from warm to cold states. The third experiment prescribes 200 years of dynamic ice-shelf geometry (with both retreating and advancing ice) based on a BISICLES simulation along with similar flips between warm and

Ocean sea-ice modelling in the Southern Ocean around Indian Antarctic stations

NASA Astrophysics Data System (ADS)

Kumar, Anurag; Dwivedi, Suneet; Rajak, D. Ram

2017-07-01

An eddy-resolving coupled ocean sea-ice modelling is carried out in the Southern Ocean region (9°-78°E; 51°-71°S) using the MITgcm. The model domain incorporates the Indian Antarctic stations, Maitri (11.7{°}E; 70.7{°}S) and Bharati (76.1{°}E; 69.4{°}S). The realistic simulation of the surface variables, namely, sea surface temperature (SST), sea surface salinity (SSS), surface currents, sea ice concentration (SIC) and sea ice thickness (SIT) is presented for the period of 1997-2012. The horizontal resolution of the model varies between 6 and 10 km. The highest vertical resolution of 5 m is taken near the surface, which gradually increases with increasing depths. The seasonal variability of the SST, SSS, SIC and currents is compared with the available observations in the region of study. It is found that the SIC of the model domain is increasing at a rate of 0.09% per month (nearly 1% per year), whereas, the SIC near Maitri and Bharati regions is increasing at a rate of 0.14 and 0.03% per month, respectively. The variability of the drift of the sea-ice is also estimated over the period of simulation. It is also found that the sea ice volume of the region increases at the rate of 0.0004 km3 per month (nearly 0.005 km3 per year). Further, it is revealed that the accumulation of sea ice around Bharati station is more as compared to Maitri station.

Ice-sheet response to oceanic forcing.

PubMed

Joughin, Ian; Alley, Richard B; Holland, David M

2012-11-30

The ice sheets of Greenland and Antarctica are losing ice at accelerating rates, much of which is a response to oceanic forcing, especially of the floating ice shelves. Recent observations establish a clear correspondence between the increased delivery of oceanic heat to the ice-sheet margin and increased ice loss. In Antarctica, most of these processes are reasonably well understood but have not been rigorously quantified. In Greenland, an understanding of the processes by which warmer ocean temperatures drive the observed retreat remains elusive. Experiments designed to identify the relevant processes are confounded by the logistical difficulties of instrumenting ice-choked fjords with actively calving glaciers. For both ice sheets, multiple challenges remain before the fully coupled ice-ocean-atmosphere models needed for rigorous sea-level projection are available.

Wave effects on ocean-ice interaction in the marginal ice zone

NASA Technical Reports Server (NTRS)

Liu, Antony K.; Hakkinen, Sirpa; Peng, Chih Y.

1993-01-01

The effects of wave train on ice-ocean interaction in the marginal ice zone are studied through numerical modeling. A coupled two-dimensional ice-ocean model has been developed to include wave effects and wind stress for the predictions of ice edge dynamics. The sea ice model is coupled to the reduced-gravity ocean model through interfacial stresses. The main dynamic balance in the ice momentum is between water-ice stress, wind stress, and wave radiation stresses. By considering the exchange of momentum between waves and ice pack through radiation stress for decaying waves, a parametric study of the effects of wave stress and wind stress on ice edge dynamics has been performed. The numerical results show significant effects from wave action. The ice edge is sharper, and ice edge meanders form in the marginal ice zone owing to forcing by wave action and refraction of swell system after a couple of days. Upwelling at the ice edge and eddy formation can be enhanced by the nonlinear effects of wave action; wave action sharpens the ice edge and can produce ice meandering, which enhances local Ekman pumping and pycnocline anomalies. The resulting ice concentration, pycnocline changes, and flow velocity field are shown to be consistent with previous observations.

Present-day Circum-Antarctic Simulations using the POPSICLES Coupled Ice Sheet-Ocean Model

NASA Astrophysics Data System (ADS)

Asay-Davis, X.; Martin, D. F.; Price, S. F.; Maltrud, M. E.; Collins, W.

2014-12-01

We present POPSICLES simulation results covering the full Antarctic Ice Sheet and the Southern Ocean spanning the period 1990 to 2010. Simulations are performed at 0.1o (~5 km) ocean resolution and with adaptive ice-sheet model resolution as fine as 500 m. We compare time-averaged melt rates below a number of major ice shelves with those reported by Rignot et al. (2013) as well as other recent studies. We also present seasonal variability and decadal trends in submarine melting from several Antarctic regions. Finally, we explore the influence on basal melting and system dynamics resulting from two different choices of climate forcing: a "normal-year" climatology and the CORE v. 2 forcing data (Large and Yeager 2008).POPSICLES couples the POP2x ocean model, a modified version of the Parallel Ocean Program (Smith and Gent, 2002), and the BISICLES ice-sheet model (Cornford et al., 2012). POP2x includes sub-ice-shelf circulation using partial top cells (Losch, 2008) and boundary layer physics following Holland and Jenkins (1999), Jenkins (2001), and Jenkins et al. (2010). Standalone POP2x output compares well with standard ice-ocean test cases (e.g., ISOMIP; Losch, 2008) and other continental-scale simulations and melt-rate observations (Kimura et al., 2013; Rignot et al., 2013). BISICLES makes use of adaptive mesh refinement and a 1st-order accurate momentum balance similar to the L1L2 model of Schoof and Hindmarsh (2009) to accurately model regions of dynamic complexity, such as ice streams, outlet glaciers, and grounding lines. Results of BISICLES simulations have compared favorably to comparable simulations with a Stokes momentum balance in both idealized tests (MISMIP-3D; Pattyn et al., 2013) and realistic configurations (Favier et al. 2014).A companion presentation, "Response of the Antarctic Ice Sheet to ocean forcing using the POPSICLES coupled ice sheet-ocean model" in session C024 covers the ice-sheet response to these melt rates in the coupled simulation

Operational use of high-resolution sst in a coupled sea ice-ocean model

NASA Astrophysics Data System (ADS)

Albretsen, A.

2003-04-01

A high-latitude, near real time, sea surface temperature (SST) product with 10 km resolution is developed at the Norwegian Meteorological Institute (met.no) through the EUMETSAT project OSI-SAF (Ocean and Sea Ice Satellite Application Facility). The product covers the Atlantic Ocean from 50N to 90N and is produced twice daily. A digitized SST and sea ice map is produced manually once a week at the Ice Mapping Service at met.no using all available information from the previous week. This map is the basis for a daily SST analysis, in which the most recent OSI-SAF SST products are successively overlaid. The resulting SST analysis field is then used in a simple data assimilation scheme in a coupled ice-ocean model to perform daily 10 days forecasts of ocean and sea ice variables. Also, the associated OSI-SAF sea ice concentration product, built from different polar orbiting satellites, is assimilated into the sea ice model. Preliminary estimates of impact on forecast skill and error statistics will be presented.

Evidence for ice-ocean albedo feedback in the Arctic Ocean shifting to a seasonal ice zone.

PubMed

Kashiwase, Haruhiko; Ohshima, Kay I; Nihashi, Sohey; Eicken, Hajo

2017-08-15

Ice-albedo feedback due to the albedo contrast between water and ice is a major factor in seasonal sea ice retreat, and has received increasing attention with the Arctic Ocean shifting to a seasonal ice cover. However, quantitative evaluation of such feedbacks is still insufficient. Here we provide quantitative evidence that heat input through the open water fraction is the primary driver of seasonal and interannual variations in Arctic sea ice retreat. Analyses of satellite data (1979-2014) and a simplified ice-upper ocean coupled model reveal that divergent ice motion in the early melt season triggers large-scale feedback which subsequently amplifies summer sea ice anomalies. The magnitude of divergence controlling the feedback has doubled since 2000 due to a more mobile ice cover, which can partly explain the recent drastic ice reduction in the Arctic Ocean.

High-frequency and meso-scale winter sea-ice variability in the Southern Ocean in a high-resolution global ocean model

NASA Astrophysics Data System (ADS)

Stössel, Achim; von Storch, Jin-Song; Notz, Dirk; Haak, Helmuth; Gerdes, Rüdiger

2018-03-01

This study is on high-frequency temporal variability (HFV) and meso-scale spatial variability (MSV) of winter sea-ice drift in the Southern Ocean simulated with a global high-resolution (0.1°) sea ice-ocean model. Hourly model output is used to distinguish MSV characteristics via patterns of mean kinetic energy (MKE) and turbulent kinetic energy (TKE) of ice drift, surface currents, and wind stress, and HFV characteristics via time series of raw variables and correlations. We find that (1) along the ice edge, the MSV of ice drift coincides with that of surface currents, in particular such due to ocean eddies; (2) along the coast, the MKE of ice drift is substantially larger than its TKE and coincides with the MKE of wind stress; (3) in the interior of the ice pack, the TKE of ice drift is larger than its MKE, mostly following the TKE pattern of wind stress; (4) the HFV of ice drift is dominated by weather events, and, in the absence of tidal currents, locally and to a much smaller degree by inertial oscillations; (5) along the ice edge, the curl of the ice drift is highly correlated with that of surface currents, mostly reflecting the impact of ocean eddies. Where ocean eddies occur and the ice is relatively thin, ice velocity is characterized by enhanced relative vorticity, largely matching that of surface currents. Along the ice edge, ocean eddies produce distinct ice filaments, the realism of which is largely confirmed by high-resolution satellite passive-microwave data.

Ice Shelf-Ocean Interactions Near Ice Rises and Ice Rumples

NASA Astrophysics Data System (ADS)

Lange, M. A.; Rückamp, M.; Kleiner, T.

2013-12-01

The stability of ice shelves depends on the existence of embayments and is largely influenced by ice rises and ice rumples, which act as 'pinning-points' for ice shelf movement. Of additional critical importance are interactions between ice shelves and the water masses underlying them in ice shelf cavities, particularly melting and refreezing processes. The present study aims to elucidate the role of ice rises and ice rumples in the context of climate change impacts on Antarctic ice shelves. However, due to their smaller spatial extent, ice rumples react more sensitively to climate change than ice rises. Different forcings are at work and need to be considered separately as well as synergistically. In order to address these issues, we have decided to deal with the following three issues explicitly: oceanographic-, cryospheric and general topics. In so doing, we paid particular attention to possible interrelationships and feedbacks in a coupled ice-shelf-ocean system. With regard to oceanographic issues, we have applied the ocean circulation model ROMBAX to ocean water masses adjacent to and underneath a number of idealized ice shelf configurations: wide and narrow as well as laterally restrained and unrestrained ice shelves. Simulations were performed with and without small ice rises located close to the calving front. For larger configurations, the impact of the ice rises on melt rates at the ice shelf base is negligible, while for smaller configurations net melting rates at the ice-shelf base differ by a factor of up to eight depending on whether ice rises are considered or not. We employed the thermo-coupled ice flow model TIM-FD3 to simulate the effects of several ice rises and one ice rumple on the dynamics of ice shelf flow. We considered the complete un-grounding of the ice shelf in order to investigate the effect of pinning points of different characteristics (interior or near calving front, small and medium sized) on the resulting flow and stress fields

Coupled ice-ocean dynamics in the marginal ice zones Upwelling/downwelling and eddy generation

NASA Technical Reports Server (NTRS)

Hakkinen, S.

1986-01-01

This study is aimed at modeling mesoscale processes such as upwelling/downwelling and ice edge eddies in the marginal ice zones. A two-dimensional coupled ice-ocean model is used for the study. The ice model is coupled to the reduced gravity ocean model through interfacial stresses. The parameters of the ocean model were chosen so that the dynamics would be nonlinear. The model was tested by studying the dynamics of upwelling. Wings parallel to the ice edge with the ice on the right produce upwelling because the air-ice momentum flux is much greater than air-ocean momentum flux; thus the Ekman transport is greater than the ice than in the open water. The stability of the upwelling and downwelling jets is discussed. The downwelling jet is found to be far more unstable than the upwelling jet because the upwelling jet is stabilized by the divergence. The constant wind field exerted on a varying ice cover will generate vorticity leading to enhanced upwelling/downwelling regions, i.e., wind-forced vortices. Steepening and strengthening of vortices are provided by the nonlinear terms. When forcing is time-varying, the advection terms will also redistribute the vorticity. The wind reversals will separate the vortices from the ice edge, so that the upwelling enhancements are pushed to the open ocean and the downwelling enhancements are pushed underneath the ice.

Transient sensitivities of sea ice export through the Canadian Arctic Archipelago inferred from a coupled ocean/sea-ice adjoint model

NASA Astrophysics Data System (ADS)

Heimbach, P.; Losch, M.; Menemenlis, D.; Campin, J.; Hill, C.

2008-12-01

The sensitivity of sea-ice export through the Canadian Arctic Archipelago (CAA), measured in terms of its solid freshwater export through Lancaster Sound, to changes in various elements of the ocean and sea-ice state, and to elements of the atmospheric forcing fields through time and space is assessed by means of a coupled ocean/sea-ice adjoint model. The adjoint model furnishes full spatial sensitivity maps (also known as Lagrange multipliers) of the export metric to a variety of model variables at any chosen point in time, providing the unique capability to quantify major drivers of sea-ice export variability. The underlying model is the MIT ocean general circulation model (MITgcm), which is coupled to a Hibler-type dynamic/thermodynamic sea-ice model. The configuration is based on the Arctic face of the ECCO3 high-resolution cubed-sphere model, but coarsened to 36-km horizontal grid spacing. The adjoint of the coupled system has been derived by means of automatic differentiation using the software tool TAF. Finite perturbation simulations are performed to check the information provided by the adjoint. The sea-ice model's performance in the presence of narrow straits is assessed with different sea-ice lateral boundary conditions. The adjoint sensitivity clearly exposes the role of the model trajectory and the transient nature of the problem. The complex interplay between forcing, dynamics, and boundary condition is demonstrated in the comparison between the different calculations. The study is a step towards fully coupled adjoint-based ocean/sea-ice state estimation at basin to global scales as part of the ECCO efforts.

Sea ice thermohaline dynamics and biogeochemistry in the Arctic Ocean: Empirical and model results

DOE Office of Scientific and Technical Information (OSTI.GOV)

Duarte, Pedro; Meyer, Amelie; Olsen, Lasse M.

Here, large changes in the sea ice regime of the Arctic Ocean have occurred over the last decades justifying the development of models to forecast sea ice physics and biogeochemistry. The main goal of this study is to evaluate the performance of the Los Alamos Sea Ice Model (CICE) to simulate physical and biogeochemical properties at time scales of a few weeks and to use the model to analyze ice algal bloom dynamics in different types of ice. Ocean and atmospheric forcing data and observations of the evolution of the sea ice properties collected from 18 April to 4 Junemore » 2015, during the Norwegian young sea ICE expedition, were used to test the CICE model. Our results show the following: (i) model performance is reasonable for sea ice thickness and bulk salinity; good for vertically resolved temperature, vertically averaged Chl a concentrations, and standing stocks; and poor for vertically resolved Chl a concentrations. (ii) Improving current knowledge about nutrient exchanges, ice algal recruitment, and motion is critical to improve sea ice biogeochemical modeling. (iii) Ice algae may bloom despite some degree of basal melting. (iv) Ice algal motility driven by gradients in limiting factors is a plausible mechanism to explain their vertical distribution. (v) Different ice algal bloom and net primary production (NPP) patterns were identified in the ice types studied, suggesting that ice algal maximal growth rates will increase, while sea ice vertically integrated NPP and biomass will decrease as a result of the predictable increase in the area covered by refrozen leads in the Arctic Ocean.« less

Sea ice thermohaline dynamics and biogeochemistry in the Arctic Ocean: Empirical and model results

DOE PAGES

Duarte, Pedro; Meyer, Amelie; Olsen, Lasse M.; ...

2017-06-08

Here, large changes in the sea ice regime of the Arctic Ocean have occurred over the last decades justifying the development of models to forecast sea ice physics and biogeochemistry. The main goal of this study is to evaluate the performance of the Los Alamos Sea Ice Model (CICE) to simulate physical and biogeochemical properties at time scales of a few weeks and to use the model to analyze ice algal bloom dynamics in different types of ice. Ocean and atmospheric forcing data and observations of the evolution of the sea ice properties collected from 18 April to 4 Junemore » 2015, during the Norwegian young sea ICE expedition, were used to test the CICE model. Our results show the following: (i) model performance is reasonable for sea ice thickness and bulk salinity; good for vertically resolved temperature, vertically averaged Chl a concentrations, and standing stocks; and poor for vertically resolved Chl a concentrations. (ii) Improving current knowledge about nutrient exchanges, ice algal recruitment, and motion is critical to improve sea ice biogeochemical modeling. (iii) Ice algae may bloom despite some degree of basal melting. (iv) Ice algal motility driven by gradients in limiting factors is a plausible mechanism to explain their vertical distribution. (v) Different ice algal bloom and net primary production (NPP) patterns were identified in the ice types studied, suggesting that ice algal maximal growth rates will increase, while sea ice vertically integrated NPP and biomass will decrease as a result of the predictable increase in the area covered by refrozen leads in the Arctic Ocean.« less

Sea ice thermohaline dynamics and biogeochemistry in the Arctic Ocean: Empirical and model results

NASA Astrophysics Data System (ADS)

Duarte, Pedro; Meyer, Amelie; Olsen, Lasse M.; Kauko, Hanna M.; Assmy, Philipp; Rösel, Anja; Itkin, Polona; Hudson, Stephen R.; Granskog, Mats A.; Gerland, Sebastian; Sundfjord, Arild; Steen, Harald; Hop, Haakon; Cohen, Lana; Peterson, Algot K.; Jeffery, Nicole; Elliott, Scott M.; Hunke, Elizabeth C.; Turner, Adrian K.

2017-07-01

Large changes in the sea ice regime of the Arctic Ocean have occurred over the last decades justifying the development of models to forecast sea ice physics and biogeochemistry. The main goal of this study is to evaluate the performance of the Los Alamos Sea Ice Model (CICE) to simulate physical and biogeochemical properties at time scales of a few weeks and to use the model to analyze ice algal bloom dynamics in different types of ice. Ocean and atmospheric forcing data and observations of the evolution of the sea ice properties collected from 18 April to 4 June 2015, during the Norwegian young sea ICE expedition, were used to test the CICE model. Our results show the following: (i) model performance is reasonable for sea ice thickness and bulk salinity; good for vertically resolved temperature, vertically averaged Chl a concentrations, and standing stocks; and poor for vertically resolved Chl a concentrations. (ii) Improving current knowledge about nutrient exchanges, ice algal recruitment, and motion is critical to improve sea ice biogeochemical modeling. (iii) Ice algae may bloom despite some degree of basal melting. (iv) Ice algal motility driven by gradients in limiting factors is a plausible mechanism to explain their vertical distribution. (v) Different ice algal bloom and net primary production (NPP) patterns were identified in the ice types studied, suggesting that ice algal maximal growth rates will increase, while sea ice vertically integrated NPP and biomass will decrease as a result of the predictable increase in the area covered by refrozen leads in the Arctic Ocean.

Ocean Tide Influences on the Antarctic and Greenland Ice Sheets

NASA Astrophysics Data System (ADS)

Padman, Laurie; Siegfried, Matthew R.; Fricker, Helen A.

2018-03-01

Ocean tides are the main source of high-frequency variability in the vertical and horizontal motion of ice sheets near their marine margins. Floating ice shelves, which occupy about three quarters of the perimeter of Antarctica and the termini of four outlet glaciers in northern Greenland, rise and fall in synchrony with the ocean tide. Lateral motion of floating and grounded portions of ice sheets near their marine margins can also include a tidal component. These tide-induced signals provide insight into the processes by which the oceans can affect ice sheet mass balance and dynamics. In this review, we summarize in situ and satellite-based measurements of the tidal response of ice shelves and grounded ice, and spatial variability of ocean tide heights and currents around the ice sheets. We review sensitivity of tide heights and currents as ocean geometry responds to variations in sea level, ice shelf thickness, and ice sheet mass and extent. We then describe coupled ice-ocean models and analytical glacier models that quantify the effect of ocean tides on lower-frequency ice sheet mass loss and motion. We suggest new observations and model developments to improve the representation of tides in coupled models that are used to predict future ice sheet mass loss and the associated contribution to sea level change. The most critical need is for new data to improve maps of bathymetry, ice shelf draft, spatial variability of the drag coefficient at the ice-ocean interface, and higher-resolution models with improved representation of tidal energy sinks.

Biogeochemical Coupling between Ocean and Sea Ice

NASA Astrophysics Data System (ADS)

Wang, S.; Jeffery, N.; Maltrud, M. E.; Elliott, S.; Wolfe, J.

2016-12-01

Biogeochemical processes in ocean and sea ice are tightly coupled at high latitudes. Ongoing changes in Arctic and Antarctic sea ice domain likely influence the coupled system, not only through physical fields but also biogeochemical properties. Investigating the system and its changes requires representation of ocean and sea ice biogeochemical cycles, as well as their coupling in Earth System Models. Our work is based on ACME-HiLAT, a new offshoot of the Community Earth System Model (CESM), including a comprehensive representation of marine ecosystems in the form of the Biogeochemical Elemental Cycling Module (BEC). A full vertical column sea ice biogeochemical module has recently been incorporated into the sea ice component. We have further introduced code modifications to couple key growth-limiting nutrients (N, Si, Fe), dissolved and particulate organic matter, and phytoplankton classes that are important in polar regions between ocean and sea ice. The coupling of ocean and sea ice biology-chemistry will enable representation of key processes such as the release of important climate active constituents or seeding algae from melting sea ice into surface waters. Sensitivity tests suggest sea ice and ocean biogeochemical coupling influences phytoplankton competition, biological production, and the CO2 flux. Sea ice algal seeding plays an important role in determining phytoplankton composition of Arctic early spring blooms, since different groups show various responses to the seeding biomass. Iron coupling leads to increased phytoplankton biomass in the Southern Ocean, which also affects carbon uptake via the biological pump. The coupling of macronutrients and organic matter may have weaker influences on the marine ecosystem. Our developments will allow climate scientists to investigate the fully coupled responses of the sea ice-ocean BGC system to physical changes in polar climate.

Influence of coupling on atmosphere, sea ice and ocean regional models in the Ross Sea sector, Antarctica

NASA Astrophysics Data System (ADS)

Jourdain, Nicolas C.; Mathiot, Pierre; Gallée, Hubert; Barnier, Bernard

2011-04-01

Air-sea ice-ocean interactions in the Ross Sea sector form dense waters that feed the global thermohaline circulation. In this paper, we develop the new limited-area ocean-sea ice-atmosphere coupled model TANGO to simulate the Ross Sea sector. TANGO is built up by coupling the atmospheric limited-area model MAR to a regional configuration of the ocean-sea ice model NEMO. A method is then developed to identify the mechanisms by which local coupling affects the simulations. TANGO is shown to simulate realistic sea ice properties and atmospheric surface temperatures. These skills are mostly related to the skills of the stand alone atmospheric and oceanic models used to build TANGO. Nonetheless, air temperatures over ocean and winter sea ice thickness are found to be slightly improved in coupled simulations as compared to standard stand alone ones. Local atmosphere ocean feedbacks over the open ocean are found to significantly influence ocean temperature and salinity. In a stand alone ocean configuration, the dry and cold air produces an ocean cooling through sensible and latent heat loss. In a coupled configuration, the atmosphere is in turn moistened and warmed by the ocean; sensible and latent heat loss is therefore reduced as compared to the stand alone simulations. The atmosphere is found to be less sensitive to local feedbacks than the ocean. Effects of local feedbacks are increased in the coastal area because of the presence of sea ice. It is suggested that slow heat conduction within sea ice could amplify the feedbacks. These local feedbacks result in less sea ice production in polynyas in coupled mode, with a subsequent reduction in deep water formation.

Efficient Flowline Simulations of Ice Shelf-Ocean Interactions: Sensitivity Studies with a Fully Coupled Model

NASA Technical Reports Server (NTRS)

Walker, Ryan Thomas; Holland, David; Parizek, Byron R.; Alley, Richard B.; Nowicki, Sophie M. J.; Jenkins, Adrian

2013-01-01

Thermodynamic flowline and plume models for the ice shelf-ocean system simplify the ice and ocean dynamics sufficiently to allow extensive exploration of parameters affecting ice-sheet stability while including key physical processes. Comparison between geophysically and laboratory-based treatments of ice-ocean interface thermodynamics shows reasonable agreement between calculated melt rates, except where steep basal slopes and relatively high ocean temperatures are present. Results are especially sensitive to the poorly known drag coefficient, highlighting the need for additional field experiments to constrain its value. These experiments also suggest that if the ice-ocean interface near the grounding line is steeper than some threshold, further steepening of the slope may drive higher entrainment that limits buoyancy, slowing the plume and reducing melting; if confirmed, this will provide a stabilizing feedback on ice sheets under some circumstances.

Ocean-ice interaction in the marginal ice zone using synthetic aperture radar imagery

NASA Technical Reports Server (NTRS)

Liu, Antony K.; Peng, Chich Y.; Weingartner, Thomas J.

1994-01-01

Ocean-ice interaction processes in the marginal ice zone (MIZ) by wind, waves, and mesoscale features, such as up/downwelling and eddies are studied using Earth Remote-Sensing Satellite (ERS) 1 synthetic aperture radar (SAR) images and an ocean-ice interaction model. A sequence of seven SAR images of the MIZ in the Chukchi Sea with 3 or 6 days interval are investigated for ice edge advance/retreat. Simultaneous current measurements from the northeast Chukchi Sea, as well as the Barrow wind record, are used to interpret the MIZ dynamics. SAR spectra of waves in ice and ocean waves in the Bering and Chukchi Sea are compared for the study of wave propagation and dominant SAR imaging mechanism. By using the SAR-observed ice edge configuration and wind and wave field in the Chukchi Sea as inputs, a numerical simulation has been performed with the ocean-ice interaction model. After 3 days of wind and wave forcing the resulting ice edge configuration, eddy formation, and flow velocity field are shown to be consistent with SAR observations.

Short-term sea ice forecasts with the RASM-ESRL coupled model: A testbed for improving simulations of ocean-ice-atmosphere interactions in the marginal ice zone

NASA Astrophysics Data System (ADS)

Solomon, A.; Cox, C. J.; Hughes, M.; Intrieri, J. M.; Persson, O. P. G.

2015-12-01

The dramatic decrease of Arctic sea-ice has led to a new Arctic sea-ice paradigm and to increased commercial activity in the Arctic Ocean. NOAA's mission to provide accurate and timely sea-ice forecasts, as explicitly outlined in the National Ocean Policy and the U.S. National Strategy for the Arctic Region, needs significant improvement across a range of time scales to improve safety for human activity. Unfortunately, the sea-ice evolution in the new Arctic involves the interaction of numerous physical processes in the atmosphere, ice, and ocean, some of which are not yet understood. These include atmospheric forcing of sea-ice movement through stress and stress deformation; atmospheric forcing of sea-ice melt and formation through energy fluxes; and ocean forcing of the atmosphere through new regions of seasonal heat release. Many of these interactions involve emerging complex processes that first need to be understood and then incorporated into forecast models in order to realize the goal of useful sea-ice forecasting. The underlying hypothesis for this study is that errors in simulations of "fast" atmospheric processes significantly impact the forecast of seasonal sea-ice retreat in summer and its advance in autumn in the marginal ice zone (MIZ). We therefore focus on short-term (0-20 day) ice-floe movement, the freeze-up and melt-back processes in the MIZ, and the role of storms in modulating stress and heat fluxes. This study uses a coupled ocean-atmosphere-seaice forecast model as a testbed to investigate; whether ocean-sea ice-atmosphere coupling improves forecasts on subseasonal time scales, where systematic biases develop due to inadequate parameterizations (focusing on mixed-phase clouds and surface fluxes), how increased atmospheric resolution of synoptic features improves the forecasts, and how initialization of sea ice area and thickness and snow depth impacts the skill of the forecasts. Simulations are validated with measurements at pan-Arctic land

Dynamics of coupled ice-ocean system in the marginal ice zone: Study of the mesoscale processes and of constitutive equations for sea ice

NASA Technical Reports Server (NTRS)

Hakkinen, S.

1984-01-01

This study is aimed at the modelling of mesoscale processed such as up/downwelling and ice edge eddies in the marginal ice zones. A 2-dimensional coupled ice-ocean model is used for the study. The ice model is coupled to the reduced gravity ocean model (f-plane) through interfacial stresses. The constitutive equations of the sea ice are formulated on the basis of the Reiner-Rivlin theory. The internal ice stresses are important only at high ice concentrations (90-100%), otherwise the ice motion is essentially free drift, where the air-ice stress is balanced by the ice-water stress. The model was tested by studying the upwelling dynamics. Winds parallel to the ice edge with the ice on the right produce upwilling because the air-ice momentum flux is much greater that air-ocean momentum flux, and thus the Ekman transport is bigger under the ice than in the open water. The upwelling simulation was extended to include temporally varying forcing, which was chosen to vary sinusoidally with a 4 day period. This forcing resembles successive cyclone passings. In the model with a thin oceanic upper layer, ice bands were formed.

Role of ice-ocean interaction on glacier instability: Results from numerical modelling applied to Petermann Glacier

NASA Astrophysics Data System (ADS)

Nick, Faezeh M.; Hubbard, Alun; van der Veen, Kees; Vieli, Andreas

2010-05-01

Calving of icebergs and bottom melting from ice shelves accounts for roughly half the ice transferred from the Greenland Ice Sheet into the surrounding ocean, and virtually all of the ice loss from the Antarctic Ice Sheet. Petermann Glacier (north Greenland) with its 16 km wide and 80 km long floating tongue, experiences massive bottom melting. We apply a numerical ice flow model using a physically-based calving criterion based on crevasse depth to investigate the contribution of processes such as bottom melting, sea ice or sikkusak disintegration, surface run off and iceberg calving to the mass balance and instability of Petermann Glacier and its ice shelf. Our modelling study provides insights into the role of ice-ocean interaction, and on how to incorporate calving in ice sheet models, improving our ability to predict future ice sheet change.

Astronomical Ice: The Effects of Treating Ice as a Porous Media on the Dynamics and Evolution of Extraterrestrial Ice-Ocean Environments

NASA Astrophysics Data System (ADS)

Buffo, J.; Schmidt, B. E.

2015-12-01

With the prevalence of water and ice rich environments in the solar system, and likely the universe, becoming more apparent, understanding the evolutionary dynamics and physical processes of such locales is of great importance. Piqued interest arises from the understanding that the persistence of all known life depends on the presence of liquid water. As in situ investigation is currently infeasible, accurate numerical modeling is the best technique to demystify these environments. We will discuss an evolving model of ice-ocean interaction aimed at realistically describing the behavior of the ice-ocean interface by treating basal ice as a porous media, and its possible implications on the formation of astrobiological niches. Treating ice as a porous media drastically affects the thermodynamic properties it exhibits. Thus inclusion of this phenomenon is critical in accurately representing the dynamics and evolution of all ice-ocean environments. This model utilizes equations that describe the dynamics of sea ice when it is treated as a porous media (Hunke et. al. 2011), coupled with a basal melt and accretion model (Holland and Jenkins 1999). Combined, these two models produce the most accurate description of the processes occurring at the base of terrestrial sea ice and ice shelves, capable of resolving variations within the ice due to environmental pressures. While these models were designed for application to terrestrial environments, the physics occurring at any ice-water interface is identical, and these models can be used to represent the evolution of a variety of icy astronomical bodies. As terrestrial ice shelves provide a close analog to planetary ice-ocean environments, we truth test the models validity against observations of ice shelves. We apply this model to the ice-ocean interface of the icy Galilean moon Europa. We include profiles of temperature, salinity, solid fraction, and Darcy velocity, as well as temporally and spatially varying melt and

The Influence of Ice-Ocean Interactions on Europa's Overturning Circulation

NASA Astrophysics Data System (ADS)

Zhu, P.; Manucharyan, G. E.; Thompson, A. F.; Goodman, J. C.; Vance, S.

2016-12-01

Jupiter's moon Europa appears to have a global liquid ocean, which is located beneath an ice shell that covers the moon's entire surface. Linking ocean dynamics and ice-ocean interactions is crucial to understanding observed surface features on Europa as well as other satellite measurements. Ocean properties and circulation may also provide clues as to whether the moon has the potential to support extraterrestrial life through chemical transport governed by ice-ocean interactions. Previous studies have identified a Hadley cell-like overturning circulation extending from the equator to mid latitudes. However, these model simulations do not consider ice-ocean interactions. In this study, our goal is to investigate how the ocean circulation may be affected by ice. We study two ice-related processes by building idealized models. One process is horizontal convection driven by an equator-to-pole buoyancy difference due to latitudinal ice transport at the ocean surface, which is found to be much weaker than the convective overturning circulation. The second process we consider is the freshwater layer formed by ice melting at the equator. A strong buoyancy contrast between the freshwater layer and the underlying water suppresses convection and turbulent mixing, which may modify the surface heat flux from the ocean to the bottom of the ice. We find that the salinity of the ocean below the freshwater layer tends to be homogeneous both vertically and horizontally with the presence of an overturning circulation. Critical values of circulation strength constrain the freshwater layer depth, and this relationship is sensitive to the average salinity of the ocean. Further coupling of temperature and salinity of the ice and the ocean that includes mutual influences between the surface heat flux and the freshwater layer may provide additional insights into the ice-ocean feedback, and its influence on the latitudinal difference of heat transport.

Amundsen Sea simulation with optimized ocean, sea ice, and thermodynamic ice shelf model parameters

NASA Astrophysics Data System (ADS)

Nakayama, Y.; Menemenlis, D.; Schodlok, M.; Heimbach, P.; Nguyen, A. T.; Rignot, E. J.

2016-12-01

Ice shelves and glaciers of the West Antarctic Ice Sheet are thinning and melting rapidly in the Amundsen Sea (AS). This is thought to be caused by warm Circumpolar Deep Water (CDW) that intrudes via submarine glacial troughs located at the continental shelf break. Recent studies, however, point out that the depth of thermocline, or thickness of Winter Water (WW, potential temperature below -1 °C located above CDW) is critical in determining the melt rate, especially for the Pine Island Glacier (PIG). For example, the basal melt rate of PIG, which decreased by 50% during summer 2012, has been attributed to thickening of WW. Despite the possible importance of WW thickness on ice shelf melting, previous modeling studies in this region have focused primarily on CDW intrusion and have evaluated numerical simulations based on bottom or deep CDW properties. As a result, none of these models have shown a good representation of WW for the AS. In this study, we adjust a small number of model parameters in a regional Amundsen and Bellingshausen Seas configuration of the Massachusetts Institute of Technology general circulation model (MITgcm) to better fit the available observations during the 2007-2010 period. We choose this time period because summer observations during these years show small interannual variability in the eastern AS. As a result of adjustments, our model shows significantly better match with observations than previous modeling studies, especially for WW. Since density of sea water depends largely on salinity at low temperature, this is crucial for assessing the impact of WW on PIG melt rate. In addition, we conduct several sensitivity studies, showing the impact of surface heat loss on the thickness and properties of WW. We also discuss some preliminary results pertaining to further optimization using the adjoint method. Our work is a first step toward improved representation of ice-shelf ocean interactions in the ECCO (Estimating the Circulation and

Role of ice-ocean interaction on glacier instability: Results from numerical modeling applied to Petermann Glacier (Invited)

NASA Astrophysics Data System (ADS)

Nick, F.; Hubbard, A.; Vieli, A.; van der Veen, C. J.; Box, J. E.; Bates, R.; Luckman, A. J.

2009-12-01

Calving of icebergs and bottom melting from ice shelves accounts for roughly half the ice transferred from the Greenland Ice Sheet into the surrounding ocean, and virtually all of the ice loss from the Antarctic Ice Sheet. Petermann Glacier (north Greenland) with its 16 km wide and 80 km long floating tongue, experiences massive bottom melting. We apply a numerical ice flow model using a physically-based calving criterion based on crevasse depth to investigate the contribution of processes such as bottom melting, sea ice or sikkusak disintegration, surface run off and iceberg calving to the mass balance and instability of Petermann Glacier and its ice shelf. Our modeling study provides insights into the role of ice-ocean interaction, and on how to incorporate calving in ice sheet models, improving our ability to predict future ice sheet change.

Spurious sea ice formation caused by oscillatory ocean tracer advection schemes

NASA Astrophysics Data System (ADS)

Naughten, Kaitlin A.; Galton-Fenzi, Benjamin K.; Meissner, Katrin J.; England, Matthew H.; Brassington, Gary B.; Colberg, Frank; Hattermann, Tore; Debernard, Jens B.

2017-08-01

Tracer advection schemes used by ocean models are susceptible to artificial oscillations: a form of numerical error whereby the advected field alternates between overshooting and undershooting the exact solution, producing false extrema. Here we show that these oscillations have undesirable interactions with a coupled sea ice model. When oscillations cause the near-surface ocean temperature to fall below the freezing point, sea ice forms for no reason other than numerical error. This spurious sea ice formation has significant and wide-ranging impacts on Southern Ocean simulations, including the disappearance of coastal polynyas, stratification of the water column, erosion of Winter Water, and upwelling of warm Circumpolar Deep Water. This significantly limits the model's suitability for coupled ocean-ice and climate studies. Using the terrain-following-coordinate ocean model ROMS (Regional Ocean Modelling System) coupled to the sea ice model CICE (Community Ice CodE) on a circumpolar Antarctic domain, we compare the performance of three different tracer advection schemes, as well as two levels of parameterised diffusion and the addition of flux limiters to prevent numerical oscillations. The upwind third-order advection scheme performs better than the centered fourth-order and Akima fourth-order advection schemes, with far fewer incidents of spurious sea ice formation. The latter two schemes are less problematic with higher parameterised diffusion, although some supercooling artifacts persist. Spurious supercooling was eliminated by adding flux limiters to the upwind third-order scheme. We present this comparison as evidence of the problematic nature of oscillatory advection schemes in sea ice formation regions, and urge other ocean/sea-ice modellers to exercise caution when using such schemes.

Modeling the interplay between sea ice formation and the oceanic mixed layer: Limitations of simple brine rejection parameterizations

NASA Astrophysics Data System (ADS)

Barthélemy, Antoine; Fichefet, Thierry; Goosse, Hugues; Madec, Gurvan

2015-02-01

The subtle interplay between sea ice formation and ocean vertical mixing is hardly represented in current large-scale models designed for climate studies. Convective mixing caused by the brine release when ice forms is likely to prevail in leads and thin ice areas, while it occurs in models at the much larger horizontal grid cell scale. Subgrid-scale parameterizations have hence been developed to mimic the effects of small-scale convection using a vertical distribution of the salt rejected by sea ice within the mixed layer, instead of releasing it in the top ocean layer. Such a brine rejection parameterization is included in the global ocean-sea ice model NEMO-LIM3. Impacts on the simulated mixed layers and ocean temperature and salinity profiles, along with feedbacks on the sea ice cover, are then investigated in both hemispheres. The changes are overall relatively weak, except for mixed layer depths, which are in general excessively reduced compared to observation-based estimates. While potential model biases prevent a definitive attribution of this vertical mixing underestimation to the brine rejection parameterization, it is unlikely that the latter can be applied in all conditions. In that case, salt rejections do not play any role in mixed layer deepening, which is unrealistic. Applying the parameterization only for low ice-ocean relative velocities improves model results, but introduces additional parameters that are not well constrained by observations.

Modelling the interplay between sea ice formation and the oceanic mixed layer: limitations of simple brine rejection parameterizations

NASA Astrophysics Data System (ADS)

Barthélemy, Antoine; Fichefet, Thierry; Goosse, Hugues; Madec, Gurvan

2015-04-01

The subtle interplay between sea ice formation and ocean vertical mixing is hardly represented in current large-scale models designed for climate studies. Convective mixing caused by the brine release when ice forms is likely to prevail in leads and thin ice areas, while it occurs in models at the much larger horizontal grid cell scale. Subgrid-scale parameterizations have hence been developed to mimic the effects of small-scale convection using a vertical distribution of the salt rejected by sea ice within the mixed layer, instead of releasing it in the top ocean layer. Such a brine rejection parameterization is included in the global ocean--sea ice model NEMO-LIM3. Impacts on the simulated mixed layers and ocean temperature and salinity profiles, along with feedbacks on the sea ice cover, are then investigated in both hemispheres. The changes are overall relatively weak, except for mixed layer depths, which are in general excessively reduced compared to observation-based estimates. While potential model biases prevent a definitive attribution of this vertical mixing underestimation to the brine rejection parameterization, it is unlikely that the latter can be applied in all conditions. In that case, salt rejections do not play any role in mixed layer deepening, which is unrealistic. Applying the parameterization only for low ice--ocean relative velocities improves model results, but introduces additional parameters that are not well constrained by observations.

Under-ice melt ponds and the oceanic mixed layer

NASA Astrophysics Data System (ADS)

Flocco, D.; Smith, N.; Feltham, D. L.

2017-12-01

Under-ice melt ponds are pools of freshwater beneath the Arctic sea ice that form when melt from the surface of the sea ice percolates down through the porous sea ice. Through double diffusion, a sheet of ice can form at the interface between the ocean and the under-ice melt pond, completely isolating the pond from the mixed layer below and forming a false bottom to the sea ice. As such, they insulate the sea ice from the ocean below. It has been estimated that these ponds could cover between 5 and 40 % of the base of the Arctic sea ice, and so could have a notable impact on the mass balance of the sea ice. We have developed a one-dimensional model to calculate the thickness and thermodynamic properties of a slab of sea ice, an under-ice melt pond, and a false bottom, as these layers evolve. Through carrying out sensitivity studies, we have identified a number of interesting ways that under-ice melt ponds affect the ice above them and the rate of basal ablation. We found that they result in thicker sea ice above them, due to their insulation of the ice, and have found a possible positive feedback cycle in which less ice will be gained due to under-ice melt ponds as the Arctic becomes warmer. More recently, we have coupled this model to a simple Kraus-Turner type model of the oceanic mixed layer to investigate how these ponds affect the ocean water beneath them. Through altering basal ablation rates and ice thickness, they change the fresh water and salt fluxes into the mixed layer, as well as incoming radiation. Multi-year simulations have, in particular, shown how these effects work on longer time-scales.

North Greenland's Ice Shelves and Ocean Warming

NASA Astrophysics Data System (ADS)

Muenchow, A.; Schauer, U.; Padman, L.; Melling, H.; Fricker, H. A.

2014-12-01

Rapid disintegration of ice shelves (the floating extensions of marine-terminating glaciers) can lead to increasing ice discharge, thinning upstream ice sheets, rising sea level. Pine Island Glacier, Antarctica, and Jacobshavn Isbrae, Greenland, provide prominent examples of these processes which evolve at decadal time scales. We here focus on three glacier systems north of 78 N in Greenland, each of which discharges more than 10 Gt per year of ice and had an extensive ice shelf a decade ago; Petermann Gletscher (PG), Niogshalvfjerdsfjorden (79N), and Zachariae Isstrom (ZI). We summarize and discuss direct observations of ocean and glacier properties for these systems as they have evolved in the northwest (PG) and northeast (79N and ZI) of Greenland over the last two decades. We use a combination of modern and historical snapshots of ocean temperature and salinity (PG, 79N, ZI), moored observations in Nares Strait (PG), and snapshots of temperature and velocity fields on the broad continental shelf off northeast Greenland (79N, ZI) collected between 1993 and 2014. Ocean warming adjacent to PG has been small relative to the ocean warming adjacent to 79N and ZI; however, ZI lost its entire ice shelf during the last decade while 79N, less than 70 km to the north of ZI, remained stable. In contrast, PG has thinned by about 10 m/y just prior to shedding two ice islands representing almost half its ice shelf area or a fifth by volume. At PG advective ice flux divergence explains about half of the dominantly basal melting while response to non-steady external forcing explains the other half. The observations at PG,79N, and ZI suggest that remotely sensed ambient surface ocean temperatures are poor proxies to explain ice shelf thinning and retreat. We posit that local dynamics of the subsurface ocean heat flux matters most. Ocean heat must first be delivered over the sill into the fjord and then within the ice shelf cavity to the base of the shelf near the grounding line

Arctic Ice-Ocean Coupling and Gyre Equilibration Observed With Remote Sensing

NASA Astrophysics Data System (ADS)

Dewey, Sarah; Morison, James; Kwok, Ronald; Dickinson, Suzanne; Morison, David; Andersen, Roger

2018-02-01

Model and observational evidence has shown that ocean current speeds in the Beaufort Gyre have increased and recently stabilized. Because these currents rival ice drift speeds, we examine the potential for the Beaufort Gyre's shift from a system in which the wind drives the ice and the ice drives a passive ocean to one in which the ocean often, in the absence of high winds, drives the ice. The resultant stress exerted on the ocean by the ice and the resultant Ekman pumping are reversed, without any change in average wind stress curl. Through these curl reversals, the ice-ocean stress provides a key feedback in Beaufort Gyre stabilization. This manuscript constitutes one of the first observational studies of ice-ocean stress inclusive of geostrophic ocean currents, by making use of recently available remote sensing data.

Observational Evidence of a Hemispheric-wide Ice-ocean Albedo Feedback Effect on Antarctic Sea-ice Decay

NASA Technical Reports Server (NTRS)

Nihashi, Sohey; Cavalieri, Donald J.

2007-01-01

The effect of ice-ocean albedo feedback (a kind of ice-albedo feedback) on sea-ice decay is demonstrated over the Antarctic sea-ice zone from an analysis of satellite-derived hemispheric sea ice concentration and European Centre for Medium-Range Weather Forecasts (ERA-40) atmospheric data for the period 1979-2001. Sea ice concentration in December (time of most active melt) correlates better with the meridional component of the wind-forced ice drift (MID) in November (beginning of the melt season) than the MID in December. This 1 month lagged correlation is observed in most of the Antarctic sea-ice covered ocean. Daily time series of ice , concentration show that the ice concentration anomaly increases toward the time of maximum sea-ice melt. These findings can be explained by the following positive feedback effect: once ice concentration decreases (increases) at the beginning of the melt season, solar heating of the upper ocean through the increased (decreased) open water fraction is enhanced (reduced), leading to (suppressing) a further decrease in ice concentration by the oceanic heat. Results obtained fi-om a simple ice-ocean coupled model also support our interpretation of the observational results. This positive feedback mechanism explains in part the large interannual variability of the sea-ice cover in summer.

The Sensitivity of the Greenland Ice Sheet to Glacial-Interglacial Oceanic Forcing

NASA Astrophysics Data System (ADS)

Tabone, I.; Blasco Navarro, J.; Robinson, A.; Alvarez-Solas, J.; Montoya, M.

2017-12-01

Up to now, the scientific community has mainly focused on the sensitivity of the Greenland Ice Sheet (GrIS) to atmospheric variations. However, several studies suggest that the enhanced ice mass loss experienced by the GrIS in the past decades is directly connected to the increasing North Atlantic temperatures. Melting of GrIS outlet glaciers triggers grounding-line retreat increasing ice discharge into the ocean. This new evidence leads to consider the ocean as a relevant driver to be taken into account when modeling the evolution of the GrIS. The ice-ocean interaction is a primary factor controling not only the likely future retreat of GrIS outlet glaciers, or the huge ice loss in past warming climates, but also, and more strongly, the past GrIS glacial expansion. The latter assumption is supported by reconstructions which propose the GrIS to be fully marine-based during glacials, and thus more exposed to the influence of the ocean. Here, for the first time, we investigate the response of the GrIS to past oceanic changes using a three-dimensional hybrid ice-sheet/ice-shelf model, which combines the Shallow Ice Approximation (SIA) for slow grounded ice sheets and the Shallow Shelf Approximation (SSA) in ice shelves and ice streams. The model accounts for a time-dependent parametrisation of the marine basal melting rate, which is used to reproduce past oceanic variations. In this work simulations of the last two glacial cycles are performed. Our results show that the GrIS is very sensitive to the ocean-triggered submarine melting (freezing). Mild oceanic temperature variations lead to a rapid retreat (expansion) of the GrIS margins, which, inducing a dynamic adjustment of the grounded ice sheet, drive the evolution of the whole ice sheet. Our results strongly suggest the need to consider the ocean as an active forcing in paleo ice sheet models.

A New Ice-sheet / Ocean Interaction Model for Greenland Fjords using High-Order Discontinuous Galerkin Methods

NASA Astrophysics Data System (ADS)

Kopera, M. A.; Maslowski, W.; Giraldo, F.

2015-12-01

One of the key outstanding challenges in modeling of climate change and sea-level rise is the ice-sheet/ocean interaction in narrow, elongated and geometrically complicated fjords around Greenland. To address this challenge we propose a new approach, a separate fjord model using discontinuous Galerkin (DG) methods, or FDG. The goal of this project is to build a separate, high-resolution module for use in Earth System Models (ESMs) to realistically represent the fjord bathymetry, coastlines, exchanges with the outside ocean, circulation and fine-scale processes occurring within the fjord and interactions at the ice shelf interface. FDG is currently at the first stage of development. The DG method provides FDG with high-order accuracy as well as geometrical flexibility, including the capacity to handle non-conforming adaptive mesh refinement to resolve the processes occurring near the ice-sheet/ocean interface without introducing prohibitive computational costs. Another benefit of this method is its excellent performance on multi- and many-core architectures, which allows for utilizing modern high performance computing systems for high-resolution simulations. The non-hydrostatic model of the incompressible Navier-Stokes equation will account for the stationary ice-shelf with sub-shelf ocean interaction, basal melting and subglacial meltwater influx and with boundary conditions at the surface to account for floating sea ice. The boundary conditions will be provided to FDG via a flux coupler to emulate the integration with an ESM. Initially, FDG will be tested for the Sermilik Fjord settings, using real bathymetry, boundary and initial conditions, and evaluated against available observations and other model results for this fjord. The overarching goal of the project is to be able to resolve the ice-sheet/ocean interactions around the entire coast of Greenland and two-way coupling with regional and global climate models such as the Regional Arctic System Model (RASM

The frequency response of a coupled ice sheet-ice shelf-ocean system to climate forcing variability

NASA Astrophysics Data System (ADS)

Goldberg, D.; Snow, K.; Jordan, J. R.; Holland, P.; Arthern, R. J.

2017-12-01

Changes at the West Antarctic ice-ocean boundary in recent decades has triggered significant increases in the regions contribution to global sea-level rise, coincident with large scale, and in some cases potentially unstable, grounding line retreat. Much of the induced change is thought to be driven by fluctuations in the oceanic heat available at the ice-ocean boundary, transported on-shelf via warm Circumpolar Deep Water (CDW). However, the processes in which ocean heat drives ice-sheet loss remains poorly understood, with observational studies routinely hindered by the extreme environment notorious to the Antarctic region. In this study we apply a novel synchronous coupled ice-ocean model, developed within the MITgcm, and are thus able to provide detailed insight into the impacts of short time scale (interannual to decadal) climate variability and feedbacks within the ice-ocean system. Feedbacks and response are assessed in an idealised ice-sheet/ocean-cavity configuration in which the far field ocean condition is adjusted to emulate periodic climate variability patterns. We reveal a non-linear response of the ice-sheet to periodic variations in thermocline depth. These non-linearities illustrate the heightened sensitivity of fast flowing ice-shelves to periodic perturbations in heat fluxes occurring at interannual and decadal time scales. The results thus highlight how small perturbations in variable climate forcing, like that of ENSO, may trigger large changes in ice-sheet response.

Research Spotlight: No tipping point for Arctic Ocean ice

NASA Astrophysics Data System (ADS)

Schultz, Colin

2011-03-01

Declines in the summer sea ice extent have led to concerns within the scientific community that the Arctic Ocean may be nearing a tipping point, beyond which the sea ice cap could not recover. In such a scenario, greenhouse gases in the atmosphere trap outgoing radiation, and as the Sun beats down 24 hours a day during the Arctic summer, temperatures rise and melt what remains of the polar sea ice cap. The Arctic Ocean, now less reflective, would absorb more of the Sun’s warmth, a feedback loop that would keep the ocean ice free. However, new research by Tietsche et al. suggests that even if the Arctic Ocean sees an ice-free summer, it would not lead to catastrophic runaway ice melt. The researchers, using a general circulation model of the global ocean and the atmosphere, found that Arctic sea ice recovers within 2 years of an imposed ice-free summer to the conditions dictated by general climate conditions during that time. Furthermore, they found that this quick recovery occurs whether the ice-free summer is triggered in 2000 or in 2060, when global temperatures are predicted to be 2°C warmer. (Geophysical Research Letters, doi:10.1029/2010GL045698, 2011)

Downscaled ice-ocean simulations for the Chukchi and Eastern Siberian Seas from an oceanic re-analysis product

NASA Astrophysics Data System (ADS)

Fujisaki-Manome, A.; Wang, J.

2016-12-01

Arctic summer sea ice has been declining at the rate that is much faster than any climate models predict. While the accelerated sea ice melting in the recent few decades could be attributed to several mechanisms such as the Arctic temperature amplification and the ice-albedo feedback, this does not necessarily explain why climate models underestimate the observed rate of summer sea ice loss. Clearly, an improved understanding is needed in what processes could be missed in climate models and could play roles in unprecedented loss of sea ice. This study evaluates contributions of sub-mesoscale processes in the ice edge (i.e. the boundary region between open water and ice covered area), which include eddies, ice bands, and the vertical mixing associated with ice bands, to the melting of sea ice and how they explain the underestimation of sea ice loss in the current state-of-art climate models. The focus area is in the pacific side of the Arctic Ocean. First, several oceanic re-analysis products including NCEP-Climate Forecast System Reanalysis (CFSR) and Modern-Era Retrospective Analysis for Research and Applications (MERRA) are evaluated in comparison with the in-situ observations from the Russian-American Long-term Census of the Arctic (RUSALCA) project. Second, the downscaled ice-ocean simulations are conducted for the Chukchi and East Siberian Seas with initial and open boundary conditions provided from a selected oceanic re-analysis product.

Investigating role of ice-ocean interaction on glacier dynamic: Results from numerical modeling applied to Petermann Glacier

NASA Astrophysics Data System (ADS)

Nick, F. M.; van der Veen, C. J.; Vieli, A.; Pattyn, F.; Hubbard, A.; Box, J. E.

2010-12-01

Calving of icebergs and bottom melting from ice shelves accounts for roughly half the ice transferred from the Greenland Ice Sheet into the surrounding ocean, and virtually all of the ice loss from the Antarctic Ice Sheet. Petermann Glacier (north Greenland) with its ~17 km wide and ~ 60 km long floating ice-shelf is experiencing high rates of bottom melting. The recent partial disintegration of its shelf (in August 2010) presents a natural experiment to investigate the dynamic response of the ice sheet to its shelf retreat. We apply a numerical ice flow model using a physically-based calving criterion based on crevasse depth to investigate the contribution of processes such as shelf disintegration, bottom melting, sea ice or sikkusak disintegration and surface run off to the mass balance of Petermann Glacier and assess its stability. Our modeling study provides insights into the role of ice-ocean interaction, and on response of Petermann Glacier to its recent massive ice loss.

Influences of Ocean Thermohaline Stratification on Arctic Sea Ice

NASA Astrophysics Data System (ADS)

Toole, J. M.; Timmermans, M.-L.; Perovich, D. K.; Krishfield, R. A.; Proshutinsky, A.; Richter-Menge, J. A.

2009-04-01

The Arctic Ocean's surface mixed layer constitutes the dynamical and thermodynamical link between the sea ice and the underlying waters. Wind stress, acting directly on the surface mixed layer or via wind-forced ice motion, produce surface currents that can in turn drive deep ocean flow. Mixed layer temperature is intimately related to basal sea ice growth and melting. Heat fluxes into or out of the surface mixed layer can occur at both its upper and lower interfaces: the former via air-sea exchange at leads and conduction through the ice, the latter via turbulent mixing and entrainment at the layer base. Variations in Arctic Ocean mixed layer properties are documented based on more than 16,000 temperature and salinity profiles acquired by Ice-Tethered Profilers since summer 2004 and analyzed in conjunction with sea ice observations from Ice Mass Balance Buoys and atmospheric heat flux estimates. Guidance interpreting the observations is provided by a one-dimensional ocean mixed layer model. The study focuses attention on the very strong density stratification about the mixed layer base in the Arctic that, in regions of sea ice melting, is increasing with time. The intense stratification greatly impedes mixed layer deepening by vertical convection and shear mixing, and thus limits the flux of deep ocean heat to the surface that could influence sea ice growth/decay. Consistent with previous work, this study demonstrates that the Arctic sea ice is most sensitive to changes in ocean mixed layer heat resulting from fluxes across its upper (air-sea and/or ice-water) interface.

A Tale of Two Forcings: Present-Day Coupled Antarctic Ice-sheet/Southern Ocean dynamics using the POPSICLES model.

NASA Astrophysics Data System (ADS)

Martin, Daniel; Asay-Davis, Xylar; Cornford, Stephen; Price, Stephen; Ng, Esmond; Collins, William

2015-04-01

We present POPSICLES simulation results covering the full Antarctic Ice Sheet and the Southern Ocean spanning the period 1990 to 2010 resulting from two different choices of climate forcing: a 'normal-year' climatology and the CORE v. 2 interannual forcing data (Large and Yeager 2008). Simulations are performed at 0.1o (~5 km) ocean resolution and adaptive ice sheet resolution as fine as 500 m. We compare time-averaged melt rates below a number of major ice shelves with those reported by Rignot et al. (2013) as well as other recent studies. We also present seasonal variability and decadal melting trends from several Antarctic regions, along with the response of the ice shelves and consequent dynamics of the grounded ice sheet. POPSICLES couples the POP2x ocean model, a modified version of the Parallel Ocean Program (Smith and Gent, 2002), and the BISICLES ice-sheet model (Cornford et al., 2012). POP2x includes sub-ice-shelf circulation using partial top cells (Losch, 2008) and boundary layer physics following Holland and Jenkins (1999), Jenkins (2001), and Jenkins et al. (2010). Standalone POP2x output compares well with standard ice-ocean test cases (e.g., ISOMIP; Losch, 2008) and other continental-scale simulations and melt-rate observations (Kimura et al., 2013; Rignot et al., 2013). BISICLES makes use of adaptive mesh refinement and a 1st-order accurate momentum balance similar to the L1L2 model of Schoof and Hindmarsh (2009) to accurately model regions of dynamic complexity, such as ice streams, outlet glaciers, and grounding lines. Results of BISICLES simulations have compared favorably to comparable simulations with a Stokes momentum balance in both idealized tests (MISMIP-3d; Pattyn et al., 2013) and realistic configurations (Favier et al. 2014).

Geophysical Investigations of Habitability in Ice-Covered Ocean Worlds

NASA Astrophysics Data System (ADS)

Vance, Steven D.; Panning, Mark P.; Stähler, Simon; Cammarano, Fabio; Bills, Bruce G.; Tobie, Gabriel; Kamata, Shunichi; Kedar, Sharon; Sotin, Christophe; Pike, William T.; Lorenz, Ralph; Huang, Hsin-Hua; Jackson, Jennifer M.; Banerdt, Bruce

2018-01-01

Geophysical measurements can reveal the structures and thermal states of icy ocean worlds. The interior density, temperature, sound speed, and electrical conductivity thus characterize their habitability. We explore the variability and correlation of these parameters using 1-D internal structure models. We invoke thermodynamic consistency using available thermodynamics of aqueous MgSO4, NaCl (as seawater), and NH3; pure water ice phases I, II, III, V, and VI; silicates; and any metallic core that may be present. Model results suggest, for Europa, that combinations of geophysical parameters might be used to distinguish an oxidized ocean dominated by MgSO4 from a more reduced ocean dominated by NaCl. In contrast with Jupiter's icy ocean moons, Titan and Enceladus have low-density rocky interiors, with minimal or no metallic core. The low-density rocky core of Enceladus may comprise hydrated minerals or anhydrous minerals with high porosity. Cassini gravity data for Titan indicate a high tidal potential Love number (k2>0.6), which requires a dense internal ocean (ρocean>1,200 kg m-3) and icy lithosphere thinner than 100 km. In that case, Titan may have little or no high-pressure ice, or a surprisingly deep water-rock interface more than 500 km below the surface, covered only by ice VI. Ganymede's water-rock interface is the deepest among known ocean worlds, at around 800 km. Its ocean may contain multiple phases of high-pressure ice, which will become buoyant if the ocean is sufficiently salty. Callisto's interior structure may be intermediate to those of Titan and Europa, with a water-rock interface 250 km below the surface covered by ice V but not ice VI.

Simulations of coupled, Antarctic ice-ocean evolution using POP2x and BISICLES (Invited)

NASA Astrophysics Data System (ADS)

Price, S. F.; Asay-Davis, X.; Martin, D. F.; Maltrud, M. E.; Hoffman, M. J.

2013-12-01

We present initial results from Antarctic, ice-ocean coupled simulations using large-scale ocean circulation and land ice evolution models. The ocean model, POP2x is a modified version of POP, a fully eddying, global-scale ocean model (Smith and Gent, 2002). POP2x allows for circulation beneath ice shelf cavities using the method of partial top cells (Losch, 2008). Boundary layer physics, which control fresh water and salt exchange at the ice-ocean interface, are implemented following Holland and Jenkins (1999), Jenkins (1999), and Jenkins et al. (2010). Standalone POP2x output compares well with standard ice-ocean test cases (e.g., ISOMIP; Losch, 2008; Kimura et al., 2013) and with results from other idealized ice-ocean coupling test cases (e.g., Goldberg et al., 2012). The land ice model, BISICLES (Cornford et al., 2012), includes a 1st-order accurate momentum balance (L1L2) and uses block structured, adaptive-mesh refinement to more accurately model regions of dynamic complexity, such as ice streams, outlet glaciers, and grounding lines. For idealized test cases focused on marine-ice sheet dynamics, BISICLES output compares very favorably relative to simulations based on the full, nonlinear Stokes momentum balance (MISMIP-3d; Pattyn et al., 2013). Here, we present large-scale (southern ocean) simulations using POP2x with fixed ice shelf geometries, which are used to obtain and validate modeled submarine melt rates against observations. These melt rates are, in turn, used to force evolution of the BISICLES model. An offline-coupling scheme, which we compare with the ice-ocean coupling work of Goldberg et al. (2012), is then used to sequentially update the sub-shelf cavity geometry seen by POP2x.

Modeling ice front Dynamics of Northwest Greenland in response to ocean thermal forcing, using ISSM and OMG data

NASA Astrophysics Data System (ADS)

Morlighem, M.; Bondzio, J. H.; Seroussi, H. L.; Wood, M.; Rignot, E. J.

2016-12-01

Glacier-front dynamics is an important control on Greenland's ice mass balance. Warmer ocean waters trigger ice-front retreats of marine-terminating glaciers, and the corresponding loss in resistive stress leads to glacier acceleration and thinning. Here, we quantify the sensitivity and vulnerability of marine-terminating glaciers along the Northwest coast of Greenland (from 73°N to 7°N) to ocean-induced melt using the Ice Sheet System Model (ISSM) and bathymetry data collected by NASA's Occreans Melting Greenland (OMG). We first combine OMG bathymetry data with ice velocity from satellites and ice thickness from airborne radars using a mass conservation approach on land to produce ice thickness and bed elevation mapping across the ice-ocean boundary that are more precise and reliable than ever before. Using this new map, we then develop a plan-view model of this region that includes a level set based moving boundary capability, a parameterized ocean-induced melt and a calving law based on a Von Mises criterion. We find that some glaciers, such as Dietrichson Gletscher or Alison Gletscher, are sensitive to small increases in ocean-induced melt, while others, such as Steenstrup Gletscher or Qeqertarsuup Sermia, are very difficult to destabilize, even with a quadrupling of the melt. Under the most intense melt experiment of 12 m/day in the summer, we find that Hayes Gletscher retreats by more than 50 km inland into a deep trough and its velocity increases by a factor of 10 over only 15 years. The model suggests that ice-ocean interactions are the triggering mechanism of glacier retreat, but the bed controls its magnitude. This work was performed at the University of California Irvine under a contract with the National Aeronautics and Space Administration, Cryospheric Sciences Program, grant NNX15AD55G.

IGLOO: an Intermediate Complexity Framework to Simulate Greenland Ice-Ocean Interactions

NASA Astrophysics Data System (ADS)

Perrette, M.; Calov, R.; Beckmann, J.; Alexander, D.; Beyer, S.; Ganopolski, A.

2017-12-01

The Greenland ice-sheet is a major contributor to current and future sea level rise associated to climate warming. It is widely believed that over a century time scale, surface melting is the main driver of Greenland ice volume change, in contrast to melting by the ocean. It is due to relatively warmer air and less ice area exposed to melting by ocean water compared to Antarctica, its southern, larger twin. Yet most modeling studies do not have adequate grid resolution to represent fine-scale outlet glaciers and fjords at the margin of the ice sheet, where ice-ocean interaction occurs, and must use rather crude parameterizations to represent this process. Additionally, the ice-sheet area grounded below sea level has been reassessed upwards in the most recent estimates of bedrock elevation under the Greenland ice sheet, revealing a larger potential for marine-mediated melting than previously thought. In this work, we develop an original approach to estimate potential Greenland ice sheet contribution to sea level rise from ocean melting, in an intermediate complexity framework, IGLOO. We use a medium-resolution (5km) ice-sheet model coupled interactively to a number of 1-D flowline models for the individual outlet glaciers. We propose a semi-objective methodology to derive 1-D glacier geometries from 2-D Greenland datasets, as well as preliminary results of coupled ice-sheet-glaciers simulations with IGLOO.

Strong and highly variable push of ocean waves on Southern Ocean sea ice.

PubMed

Stopa, Justin E; Sutherland, Peter; Ardhuin, Fabrice

2018-06-05

Sea ice in the Southern Ocean has expanded over most of the past 20 y, but the decline in sea ice since 2016 has taken experts by surprise. This recent evolution highlights the poor performance of numerical models for predicting extent and thickness, which is due to our poor understanding of ice dynamics. Ocean waves are known to play an important role in ice break-up and formation. In addition, as ocean waves decay, they cause a stress that pushes the ice in the direction of wave propagation. This wave stress could not previously be quantified due to insufficient observations at large scales. Sentinel-1 synthetic aperture radars (SARs) provide high-resolution imagery from which wave height is measured year round encompassing Antarctica since 2014. Our estimates give an average wave stress that is comparable to the average wind stress acting over 50 km of sea ice. We further reveal highly variable half-decay distances ranging from 400 m to 700 km, and wave stresses from 0.01 to 1 Pa. We expect that this variability is related to ice properties and possibly different floe sizes and ice thicknesses. A strong feedback of waves on sea ice, via break-up and rafting, may be the cause of highly variable sea-ice properties.

The Met Office Coupled Atmosphere/Land/Ocean/Sea-Ice Data Assimilation System

NASA Astrophysics Data System (ADS)

Lea, Daniel; Mirouze, Isabelle; Martin, Matthew; Hines, Adrian; Guiavarch, Catherine; Shelly, Ann

2014-05-01

The Met Office has developed a weakly-coupled data assimilation (DA) system using the global coupled model HADGEM3 (Hadley Centre Global Environment Model, version 3). This model combines the atmospheric model UM (Unified Model) at 60 km horizontal resolution on 85 vertical levels, the ocean model NEMO (Nucleus for European Modeling of the Ocean) at 25 km (at the equator) horizontal resolution on 75 vertical levels, and the sea-ice model CICE at the same resolution as NEMO. The atmosphere and the ocean/sea-ice fields are coupled every 1-hour using the OASIS coupler. The coupled model is corrected using two separate 6-hour window data assimilation systems: a 4D-Var for the atmosphere with associated soil moisture content nudging and snow analysis schemes on the one hand, and a 3D-Var FGAT for the ocean and sea-ice on the other hand. The background information in the DA systems comes from a previous 6-hour forecast of the coupled model. To show the impact of coupled DA, one-month experiments have been carried out, including 1) a full atmosphere/land/ocean/sea-ice coupled DA run, 2) an atmosphere-only run forced by OSTIA SSTs and sea-ice with atmosphere and land DA, and 3) an ocean-only run forced by atmospheric fields from run 2 with ocean and sea-ice DA. In addition, 5-day forecast runs, started twice a day, have been produced from initial conditions generated by either run 1 or a combination of runs 2 and 3. The different results have been compared to each other and, whenever possible, to other references such as the Met Office atmosphere and ocean operational analyses or the OSTIA data. These all show the coupled DA system functioning well. Evidence of imbalances and initialisation shocks has also been looked for.

The sources of Antarctic bottom water in a global ice ocean model

NASA Astrophysics Data System (ADS)

Goosse, Hugues; Campin, Jean-Michel; Tartinville, Benoı̂t

Two mechanisms contribute to the formation of Antarctic bottom water (AABW). The first, and probably the most important, is initiated by the brine released on the Antarctic continental shelf during ice formation which is responsible for an increase in salinity. After mixing with ambient water at the shelf break, this salty and dense water sinks along the shelf slope and invades the deepest part of the global ocean. For the second one, the increase of surface water density is due to strong cooling at the ocean-atmosphere interface, together with a contribution from brine release. This induces deep convection and the renewal of deep waters. The relative importance of these two mechanisms is investigated in a global coupled ice-ocean model. Chlorofluorocarbon (CFC) concentrations simulated by the model compare favourably with observations, suggesting a reasonable deep water ventilation in the Southern Ocean, except close to Antarctica where concentrations are too high. Two artificial passive tracers released at surface on the Antarctic continental shelf and in the open-ocean allow to show clearly that the two mechanisms contribute significantly to the renewal of AABW in the model. This indicates that open-ocean convection is overestimated in our simulation. Additional experiments show that the amount of AABW production due to the export of dense shelf waters is quite sensitive to the parameterisation of the effect of downsloping and meso-scale eddies. Nevertheless, shelf waters always contribute significantly to deep water renewal. Besides, increasing the P.R. Gent, J.C. McWilliams [Journal of Physical Oceanography 20 (1990) 150-155] thickness diffusion can nearly suppress the AABW formation by open-ocean convection.

Ocean interactions with the base of Amery Ice Shelf, Antarctica

NASA Technical Reports Server (NTRS)

Hellmer, Hartmut H.; Jacobs, Stanley S.

1992-01-01

Using a two-dimensional ocean themohaline circulation model, we varied the cavity shape beneath Amery Ice Shelf in an attempt to reproduce the 150-m-thick marine ice layer observed at the 'G1' ice core site. Most simulations caused melting rates which decrease the ice thickness by as much as 400 m between grounding line and G1, but produce only minor accumulation at the ice core site and closer to the ice front. Changes in the sea floor and ice topographies revealed a high sensitivity of the basal mass balance to water column thickness near the grounding line, to submarine sills, and to discontinuities in ice thickness. Model results showed temperature/salinity gradients similar to observations from beneath other ice shelves where ice is melting into seawater. Modeled outflow characteristics at the ice front are in general agreement with oceanographic data from Prydz Bay. We concur with Morgan's inference that the G1 core may have been taken in a basal crevasse filled with marine ice. This ice is formed from water cooled by ocean/ice shelf interactions along the interior ice shelf base.

The evolution of a coupled ice shelf-ocean system under different climate states

NASA Astrophysics Data System (ADS)

Grosfeld, Klaus; Sandhäger, Henner

2004-07-01

Based on a new approach for coupled applications of an ice shelf model and an ocean general circulation model, we investigate the evolution of an ice shelf-ocean system and its sensitivity to changed climatic boundary conditions. Combining established 3D models into a coupled model system enabled us to study the reaction and feedbacks of each component to changes at their interface, the ice shelf base. After calculating the dynamics for prescribed initial ice shelf and bathymetric geometries, the basal mass balance determines the system evolution. In order to explore possible developments for given boundary conditions, an idealized geometry has been chosen, reflecting basic features of the Filchner-Ronne Ice Shelf, Antarctica. The model system is found to be especially sensitive in regions where high ablation or accretion rates occur. Ice Shelf Water formation as well as the build up of a marine ice body, resulting from accretion of marine ice, is simulated, indicating strong interaction processes. To improve consistency between modeled and observed ice shelf behavior, we incorporate the typical cycle of steady ice front advance and sudden retreat due to tabular iceberg calving in our time-dependent simulations. Our basic hypothesis is that iceberg break off is associated with abrupt crack propagation along elongated anomalies of the inherent stress field of the ice body. This new concept yields glaciologically plausible results and represents an auspicious basis for the development of a thorough calving criterion. Experiments under different climatic conditions (ocean warming of 0.2 and 0.5 °C and doubled surface accumulation rates) show the coupled model system to be sensitive especially to ocean warming. Increased basal melt rates of 100% for the 0.5 °C ocean warming scenario and an asymmetric development of ice shelf thicknesses suggest a high vulnerability of ice shelf regions, which represent pivotal areas between the Antarctic Ice Sheet and the Southern

Liquid Water Oceans in Ice Giants

NASA Technical Reports Server (NTRS)

Wiktorowicz, Sloane J.; Ingersoll, Andrew P.

2007-01-01

Aptly named, ice giants such as Uranus and Neptune contain significant amounts of water. While this water cannot be present near the cloud tops, it must be abundant in the deep interior. We investigate the likelihood of a liquid water ocean existing in the hydrogen-rich region between the cloud tops and deep interior. Starting from an assumed temperature at a given upper tropospheric pressure (the photosphere), we follow a moist adiabat downward. The mixing ratio of water to hydrogen in the gas phase is small in the photosphere and increases with depth. The mixing ratio in the condensed phase is near unity in the photosphere and decreases with depth; this gives two possible outcomes. If at some pressure level the mixing ratio of water in the gas phase is equal to that in the deep interior, then that level is the cloud base. The gas below the cloud base has constant mixing ratio. Alternately, if the mixing ratio of water in the condensed phase reaches that in the deep interior, then the surface of a liquid ocean will occur. Below this ocean surface, the mixing ratio of water will be constant. A cloud base occurs when the photospheric temperature is high. For a family of ice giants with different photospheric temperatures, the cooler ice giants will have warmer cloud bases. For an ice giant with a cool enough photospheric temperature, the cloud base will exist at the critical temperature. For still cooler ice giants, ocean surfaces will result. A high mixing ratio of water in the deep interior favors a liquid ocean. We find that Neptune is both too warm (photospheric temperature too high) and too dry (mixing ratio of water in the deep interior too low) for liquid oceans to exist at present. To have a liquid ocean, Neptune s deep interior water to gas ratio would have to be higher than current models allow, and the density at 19 kbar would have to be approx. equal to 0.8 g/cu cm. Such a high density is inconsistent with gravitational data obtained during the Voyager

Outlet Glacier-Ice Shelf-Ocean Interactions: Is the Tail Wagging the Dog?

NASA Astrophysics Data System (ADS)

Parizek, B. R.; Walker, R. T.; Rinehart, S. K.

2009-12-01

While the massive interior regions of the Antarctic and Greenland Ice Sheets are presently ``resting quietly", the lower elevations of many outlet glaciers are experiencing dramatic adjustments due to changes in ice dynamics and/or surface mass balance. Oceanic and/or atmospheric forcing in these marginal regions often leads to mass deficits for entire outlet basins. Therefore, coupling the wagging tail of ice-ocean interactions with the vast ice-sheet reservoirs is imperative for accurate assessments of future sea-level rise. To study ice-ocean dynamic processes, we couple an ocean-plume model that simulates ice-shelf basal melting rates based on temperature and salinity profiles combined with plume dynamics associated with the geometry of the ice-shelf cavity (following Jenkins, 1991 and Holland and Jenkins, 1999) with a two-dimensional, isothermal model of outlet glacier-ice shelf flow (as used in Alley et al., 2007; Walker et al., 2008; Parizek et al., in review). Depending on the assigned temperature and salinity profiles, the ocean model can simulate both water-mass end-members: either cold High Salinity Shelf Water (HSSW) or relatively warm Circumpolar Deep Water (CDW), as well as between-member conditions. Notably, the coupled system exhibits sensitivity to the initial conditions. In particular, melting concentrated near the grounding line has the greatest effect in forcing grounding-line retreat. Retreat is further enhanced by a positive feedback between the ocean and ice, as the focused melt near the grounding line leads to an increase in the local slope of the basal ice, thereby enhancing buoyancy-driven plume flow and subsequent melt rates.

Century/millennium internal climate oscillations in an ocean-atmosphere-continental ice sheet model

NASA Technical Reports Server (NTRS)

Birchfield, Edward G.; Wang, Huaxiao; Rich, Jonathan J.

1994-01-01

We demonstrate in a simple climate model that there exist nonlinear feedbacks between the atmosphere, ocean, and ice sheets capable of producing century/millennium timescale internal oscillations resembling those seen in the paleoclimate record. Feedbacks involve meridional heat and salt transports in the North Atlantic, surface ocean freshwater fluxes associated with melting and growing continental ice sheets in the northen hemisphere and with Atlantic to Pacific water vapor transport. The positive feedback between the production of North Atlantic Deep Water (NADW) and the meridional salt transport by the Atlantic thermohaline circulation tends to destabilize the climate system, while the negative feedback between the freshwater flux, either to or from the continental ice sheets, and meridional heat flux to the high-latitude North Atlantic, accomplished by the thermohaline circulation, stabilizes the system. The thermohaline circulation plays a central role in both positive and negative feedbacks because of its transport of both heat and salt. Because of asymmetries between the growth and melt phases the oscillations are, in general, accompanied by a growing or decreasing ice volume over each cycle, which in the model is reflected by increasing or decreasing mean salinity.

Strong control of Southern Ocean cloud reflectivity by ice-nucleating particles

NASA Astrophysics Data System (ADS)

Vergara-Temprado, Jesús; Miltenberger, Annette K.; Furtado, Kalli; Grosvenor, Daniel P.; Shipway, Ben J.; Hill, Adrian A.; Wilkinson, Jonathan M.; Field, Paul R.; Murray, Benjamin J.; Carslaw, Ken S.

2018-03-01

Large biases in climate model simulations of cloud radiative properties over the Southern Ocean cause large errors in modeled sea surface temperatures, atmospheric circulation, and climate sensitivity. Here, we combine cloud-resolving model simulations with estimates of the concentration of ice-nucleating particles in this region to show that our simulated Southern Ocean clouds reflect far more radiation than predicted by global models, in agreement with satellite observations. Specifically, we show that the clouds that are most sensitive to the concentration of ice-nucleating particles are low-level mixed-phase clouds in the cold sectors of extratropical cyclones, which have previously been identified as a main contributor to the Southern Ocean radiation bias. The very low ice-nucleating particle concentrations that prevail over the Southern Ocean strongly suppress cloud droplet freezing, reduce precipitation, and enhance cloud reflectivity. The results help explain why a strong radiation bias occurs mainly in this remote region away from major sources of ice-nucleating particles. The results present a substantial challenge to climate models to be able to simulate realistic ice-nucleating particle concentrations and their effects under specific meteorological conditions.

Strong control of Southern Ocean cloud reflectivity by ice-nucleating particles.

PubMed

Vergara-Temprado, Jesús; Miltenberger, Annette K; Furtado, Kalli; Grosvenor, Daniel P; Shipway, Ben J; Hill, Adrian A; Wilkinson, Jonathan M; Field, Paul R; Murray, Benjamin J; Carslaw, Ken S

2018-03-13

Large biases in climate model simulations of cloud radiative properties over the Southern Ocean cause large errors in modeled sea surface temperatures, atmospheric circulation, and climate sensitivity. Here, we combine cloud-resolving model simulations with estimates of the concentration of ice-nucleating particles in this region to show that our simulated Southern Ocean clouds reflect far more radiation than predicted by global models, in agreement with satellite observations. Specifically, we show that the clouds that are most sensitive to the concentration of ice-nucleating particles are low-level mixed-phase clouds in the cold sectors of extratropical cyclones, which have previously been identified as a main contributor to the Southern Ocean radiation bias. The very low ice-nucleating particle concentrations that prevail over the Southern Ocean strongly suppress cloud droplet freezing, reduce precipitation, and enhance cloud reflectivity. The results help explain why a strong radiation bias occurs mainly in this remote region away from major sources of ice-nucleating particles. The results present a substantial challenge to climate models to be able to simulate realistic ice-nucleating particle concentrations and their effects under specific meteorological conditions. Copyright © 2018 the Author(s). Published by PNAS.

Strong control of Southern Ocean cloud reflectivity by ice-nucleating particles

PubMed Central

Miltenberger, Annette K.; Furtado, Kalli; Grosvenor, Daniel P.; Shipway, Ben J.; Hill, Adrian A.; Wilkinson, Jonathan M.; Field, Paul R.

2018-01-01

Large biases in climate model simulations of cloud radiative properties over the Southern Ocean cause large errors in modeled sea surface temperatures, atmospheric circulation, and climate sensitivity. Here, we combine cloud-resolving model simulations with estimates of the concentration of ice-nucleating particles in this region to show that our simulated Southern Ocean clouds reflect far more radiation than predicted by global models, in agreement with satellite observations. Specifically, we show that the clouds that are most sensitive to the concentration of ice-nucleating particles are low-level mixed-phase clouds in the cold sectors of extratropical cyclones, which have previously been identified as a main contributor to the Southern Ocean radiation bias. The very low ice-nucleating particle concentrations that prevail over the Southern Ocean strongly suppress cloud droplet freezing, reduce precipitation, and enhance cloud reflectivity. The results help explain why a strong radiation bias occurs mainly in this remote region away from major sources of ice-nucleating particles. The results present a substantial challenge to climate models to be able to simulate realistic ice-nucleating particle concentrations and their effects under specific meteorological conditions. PMID:29490918

Sensitivity of an Antarctic Ice Sheet Model to Sub-Ice-Shelf Melting

NASA Astrophysics Data System (ADS)

Lipscomb, W. H.; Leguy, G.; Urban, N. M.; Berdahl, M.

2017-12-01

Theory and observations suggest that marine-based sectors of the Antarctic ice sheet could retreat rapidly under ocean warming and increased melting beneath ice shelves. Numerical models of marine ice sheets vary widely in sensitivity, depending on grid resolution and the parameterization of key processes (e.g., calving and hydrofracture). Here we study the sensitivity of the Antarctic ice sheet to ocean warming and sub-shelf melting in standalone simulations of the Community Ice Sheet Model (CISM). Melt rates either are prescribed based on observations and high-resolution ocean model output, or are derived from a plume model forced by idealized ocean temperature profiles. In CISM, we vary the model resolution (between 1 and 8 km), Stokes approximation (shallow-shelf, depth-integrated higher-order, or 3D higher-order) and calving scheme to create an ensemble of plausible responses to sub-shelf melting. This work supports a broader goal of building statistical and reduced models that can translate large-scale Earth-system model projections to changes in Antarctic ocean temperatures and ice sheet discharge, thus better quantifying uncertainty in Antarctic-sourced sea-level rise.

On the assimilation of absolute geodetic dynamic topography in a global ocean model: impact on the deep ocean state

NASA Astrophysics Data System (ADS)

Androsov, Alexey; Nerger, Lars; Schnur, Reiner; Schröter, Jens; Albertella, Alberta; Rummel, Reiner; Savcenko, Roman; Bosch, Wolfgang; Skachko, Sergey; Danilov, Sergey

2018-05-01

General ocean circulation models are not perfect. Forced with observed atmospheric fluxes they gradually drift away from measured distributions of temperature and salinity. We suggest data assimilation of absolute dynamical ocean topography (DOT) observed from space geodetic missions as an option to reduce these differences. Sea surface information of DOT is transferred into the deep ocean by defining the analysed ocean state as a weighted average of an ensemble of fully consistent model solutions using an error-subspace ensemble Kalman filter technique. Success of the technique is demonstrated by assimilation into a global configuration of the ocean circulation model FESOM over 1 year. The dynamic ocean topography data are obtained from a combination of multi-satellite altimetry and geoid measurements. The assimilation result is assessed using independent temperature and salinity analysis derived from profiling buoys of the AGRO float data set. The largest impact of the assimilation occurs at the first few analysis steps where both the model ocean topography and the steric height (i.e. temperature and salinity) are improved. The continued data assimilation over 1 year further improves the model state gradually. Deep ocean fields quickly adjust in a sustained manner: A model forecast initialized from the model state estimated by the data assimilation after only 1 month shows that improvements induced by the data assimilation remain in the model state for a long time. Even after 11 months, the modelled ocean topography and temperature fields show smaller errors than the model forecast without any data assimilation.

Early ice retreat and ocean warming may induce copepod biogeographic boundary shifts in the Arctic Ocean

NASA Astrophysics Data System (ADS)

Feng, Zhixuan; Ji, Rubao; Campbell, Robert G.; Ashjian, Carin J.; Zhang, Jinlun

2016-08-01

Early ice retreat and ocean warming are changing various facets of the Arctic marine ecosystem, including the biogeographic distribution of marine organisms. Here an endemic copepod species, Calanus glacialis, was used as a model organism, to understand how and why Arctic marine environmental changes may induce biogeographic boundary shifts. A copepod individual-based model was coupled to an ice-ocean-ecosystem model to simulate temperature- and food-dependent copepod life history development. Numerical experiments were conducted for two contrasting years: a relatively cold and normal sea ice year (2001) and a well-known warm year with early ice retreat (2007). Model results agreed with commonly known biogeographic distributions of C. glacialis, which is a shelf/slope species and cannot colonize the vast majority of the central Arctic basins. Individuals along the northern boundaries of this species' distribution were most susceptible to reproduction timing and early food availability (released sea ice algae). In the Beaufort, Chukchi, East Siberian, and Laptev Seas where severe ocean warming and loss of sea ice occurred in summer 2007, relatively early ice retreat, elevated ocean temperature (about 1-2°C higher than 2001), increased phytoplankton food, and prolonged growth season created favorable conditions for C. glacialis development and caused a remarkable poleward expansion of its distribution. From a pan-Arctic perspective, despite the great heterogeneity in the temperature and food regimes, common biogeographic zones were identified from model simulations, thus allowing a better characterization of habitats and prediction of potential future biogeographic boundary shifts.

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

NASA Technical Reports Server (NTRS)

Powel, Dylan C.; Markus, Thorsten; Stoessel, Achim

2003-01-01

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.

Retention of ice-associated amphipods: possible consequences for an ice-free Arctic Ocean.

PubMed

Berge, J; Varpe, O; Moline, M A; Wold, A; Renaud, P E; Daase, M; Falk-Petersen, S

2012-12-23

Recent studies predict that the Arctic Ocean will have ice-free summers within the next 30 years. This poses a significant challenge for the marine organisms associated with the Arctic sea ice, such as marine mammals and, not least, the ice-associated crustaceans generally considered to spend their entire life on the underside of the Arctic sea ice. Based upon unique samples collected within the Arctic Ocean during the polar night, we provide a new conceptual understanding of an intimate connection between these under-ice crustaceans and the deep Arctic Ocean currents. We suggest that downwards vertical migrations, followed by polewards transport in deep ocean currents, are an adaptive trait of ice fauna that both increases survival during ice-free periods of the year and enables re-colonization of sea ice when they ascend within the Arctic Ocean. From an evolutionary perspective, this may have been an adaptation allowing success in a seasonally ice-covered Arctic. Our findings may ultimately change the perception of ice fauna as a biota imminently threatened by the predicted disappearance of perennial sea ice.

Antarctic icebergs melt over the Southern Ocean : Climatology and impact on sea ice

NASA Astrophysics Data System (ADS)

Merino, Nacho; Le Sommer, Julien; Durand, Gael; Jourdain, Nicolas C.; Madec, Gurvan; Mathiot, Pierre; Tournadre, Jean

2016-08-01

Recent increase in Antarctic freshwater release to the Southern Ocean is suggested to contribute to change in water masses and sea ice. However, climate models differ in their representation of the freshwater sources. Recent improvements in altimetry-based detection of small icebergs and in estimates of the mass loss of Antarctica may help better constrain the values of Antarctic freshwater releases. We propose a model-based seasonal climatology of iceberg melt over the Southern Ocean using state-of-the-art observed glaciological estimates of the Antarctic mass loss. An improved version of a Lagrangian iceberg model is coupled with a global, eddy-permitting ocean/sea ice model and compared to small icebergs observations. Iceberg melt increases sea ice cover, about 10% in annual mean sea ice volume, and decreases sea surface temperature over most of the Southern Ocean, but with distinctive regional patterns. Our results underline the importance of improving the representation of Antarctic freshwater sources. This can be achieved by forcing ocean/sea ice models with a climatological iceberg fresh-water flux.

RTopo-2: A global high-resolution dataset of ice sheet topography, ice shelf cavity geometry and ocean bathymetry

NASA Astrophysics Data System (ADS)

Timmermann, Ralph; Schaffer, Janin

2016-04-01

The RTopo-1 data set of Antarctic ice sheet/shelf geometry and global ocean bathymetry has proven useful not only for modelling studies of ice-ocean interaction in the southern hemisphere. Following the spirit of this data set, we introduce a new product (RTopo-2) that contains consistent maps of global ocean bathymetry, upper and lower ice surface topographies for Greenland and Antarctica, and global surface height on a spherical grid with now 30 arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO_2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the International Bathymetric Chart of the Southern Ocean (IBCSO) version 1. To achieve a good representation of the fjord and shelf bathymetry around the Greenland continent, we corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Helheim Glacier assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model including all available multibeam survey data for the region. Radar data for ice surface and ice base topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities. The data set is available in full and in regional subsets in NetCDF format from the PANGAEA database.

Evolution of Planetary Ice-Ocean Systems: Effects of Salinity

NASA Astrophysics Data System (ADS)

Allu Peddinti, D.; McNamara, A. K.

2015-12-01

Planetary oceanography is enjoying renewed attention thanks to not only the detection of several exoplanetary ocean worlds but also due to the expanding family of ocean worlds within our own star system. Our solar system is now believed to host about nine ocean worlds including Earth, some dwarf planets and few moons of Jupiter and Saturn. Amongst them, Europa, like Earth is thought to have an ice Ih-liquid water system. However, the thickness of the Europan ice-ocean system is much larger than that of the Earth. The evolution of this system would determine the individual thicknesses of the ice shell and the ocean. In turn, these thicknesses can alter the course of evolution of the system. In a pure H2O system, the thickness of the ice shell would govern if heat loss occurs entirely by conduction or if the shell begins to convect as it attains a threshold thickness. This switch between conduction-convection regimes could determine the longevity of the subsurface ocean and hence define the astrobiological potential of the planetary body at any given time. In reality, however, the system is not pure water ice. The detected induced magnetic field infers a saline ocean layer. Salts are expected to act as an anti-freeze allowing a subsurface ocean to persist over long periods but the amount of salts would determine the extent of that effect. In our current study, we use geodynamic models to examine the effect of salinity on the evolution of ice-ocean system. An initial ocean with different salinities is allowed to evolve. The effect of salinity on thickness of the two layers at any time is examined. We also track how salinity controls the switch between conductive-convective modes. The study shows that for a given time period, larger salinities can maintain a thick vigorously convecting ocean while the smaller salinities behave similar to a pure H2O system leading to a thick convecting ice-shell. A range of salinities identified can potentially predict the current state

A laboratory examination of the three-equation model of ice-ocean interactions

NASA Astrophysics Data System (ADS)

McConnochie, Craig; Kerr, Ross

2017-11-01

Numerical models of ice-ocean interactions are typically unable to resolve the transport of heat and salt to the ice face. As such, models rely upon parameterizations that have not been properly validated by data. Recent laboratory experiments of ice-saltwater interactions allow us to test the standard parameterization of heat and salt transport to ice faces - the `three equation model'. We find a significant disagreement in the dependence of the melt rate on the fluid velocity. The three-equation model predicts that the melt rate is proportional to the fluid velocity while the experimental results typically show that the melt rate is independent of the fluid velocity. By considering a theoretical analysis of the boundary layer next to a melting ice face we suggest a resolution to this disagreement. We show that the three-equation model assumes that the thickness of the diffusive sublayer is set by a shear instability. However, at low flow velocities, the sublayer is instead set by a convective instability. This distinction leads to a threshold velocity of approximately 4 cm/s at geophysically relevant conditions, above which the form of the parameterization should be valid. In contrast, at flow speeds below 4 cm/s, the three-equation model will underestimate the melt rate. ARC DP120102772.

Ocean Wave-to-Ice Energy Transfer Determined from Seafloor Pressure and Ice Shelf Seismic Observations

NASA Astrophysics Data System (ADS)

Chen, Z.; Bromirski, P. D.; Gerstoft, P.; Stephen, R. A.; Wiens, D.; Aster, R. C.; Nyblade, A.

2017-12-01

Ice shelves play an important role in buttressing land ice from reaching the sea, thus restraining the rate of sea level rise. Long-period gravity wave impacts excite vibrations in ice shelves that may trigger tabular iceberg calving and/or ice shelf collapse events. Three kinds of seismic plate waves were continuously observed by broadband seismic arrays on the Ross Ice Shelf (RIS) and on the Pine Island Glacier (PIG) ice shelf: (1) flexural-gravity waves, (2) flexural waves, and (3) extensional Lamb waves, suggesting that all West Antarctic ice shelves are subjected to similar gravity wave excitation. Ocean gravity wave heights were estimated from pressure perturbations recorded by an ocean bottom differential pressure gauge at the RIS front, water depth 741 m, about 8 km north of an on-ice seismic station that is 2 km from the shelf front. Combining the plate wave spectrum, the frequency-dependent energy transmission and reflection at the ice-water interface were determined. In addition, Young's modulus and Poisson's ratio of the RIS are estimated from the plate wave motions, and compared with the widely used values. Quantifying these ice shelf parameters from observations will improve modeling of ice shelf response to ocean forcing, and ice shelf evolution.

Antarctic Climate Variability: Covariance of Ozone and Sea Ice in Atmosphere - Ocean Coupled Model Simulations

NASA Astrophysics Data System (ADS)

Jrrar, Amna; Abraham, N. Luke; Pyle, John A.; Holland, David

2014-05-01

Changes in sea ice significantly modulate climate change because of its high reflective and insulating nature. While Arctic Sea Ice Extent (SIE) shows a negative trend. Antarctic SIE shows a weak but positive trend, estimated at 0.127 x 106 km2 per decade. The trend results from large regional cancellations, more ice in the Weddell and the Ross seas, and less ice in the Amundsen - Bellingshausen seas. A number of studies had demonstrated that stratospheric ozone depletion has had a major impact on the atmospheric circulation, causing a positive trend in the Southern Annular Mode (SAM), which has been linked to the observed positive trend in autumn sea ice in the Ross Sea. However, other modelling studies show that models forced with prescribed ozone hole simulate decreased sea ice in all regions comparative to a control run. A recent study has also shown that stratospheric ozone recovery will mitigate Antarctic sea ice loss. To verify this assumed relationship, it is important first to investigate the covariance between ozone's natural (dynamical) variability and Antarctic sea ice distribution in pre-industrial climate, to estimate the trend due to natural variability. We investigate the relationship between anomalous Antarctic ozone years and the subsequent changes in Antarctic sea ice distribution in a multidecadal control simulation using the AO-UMUKCA model. The model has a horizontal resolution of 3.75 X 2.5 degrees in longitude and latitude; and 60 hybrid height levels in the vertical, from the surface up to a height of 84 km. The ocean component is the NEMO ocean model on the ORCA2 tripolar grid, and the sea ice model is CICE. We evaluate the model's performance in terms of sea ice distribution, and we calculate sea ice extent trends for composites of anomalously low versus anomalously high SH polar ozone column. We apply EOF analysis to the seasonal anomalies of sea ice concentration, MSLP, and Z 500, and identify the leading climate modes controlling the

Sea-level response to ice sheet evolution: An ocean perspective

NASA Technical Reports Server (NTRS)

Jacobs, Stanley S.

1991-01-01

The ocean's influence upon and response to Antarctic ice sheet changes is considered in relation to sea level rise over recent and future decades. Assuming present day ice fronts are in approximate equilibrium, a preliminary budget for the ice sheet is estimated from accumulation vs. iceberg calving and the basal melting that occurs beneath floating ice shelves. Iceberg calving is derived from the volume of large bergs identified and tracked by the Navy/NOAA Joint Ice Center and from shipboard observations. Basal melting exceeds 600 cu km/yr and is concentrated near the ice fronts and ice shelf grounding lines. An apparent negative mass balance for the Antarctic ice sheet may result from an anomalous calving rate during the past decade, but there are large uncertainties associated with all components of the ice budget. The results from general circulation models are noted in the context of projected precipitation increases and ocean temperature changes on and near the continent. An ocean research program that could help refine budget estimates is consistent with goals of the West Antarctic Ice Sheet Initiative.

The interaction between sea ice and salinity-dominated ocean circulation: implications for halocline stability and rapid changes of sea-ice cover

NASA Astrophysics Data System (ADS)

Jensen, M. F.; Nilsson, J.; Nisancioglu, K. H.

2016-02-01

In this study, we develop a simple conceptual model to examine how interactions between sea ice and oceanic heat and freshwater transports affect the stability of an upper-ocean halocline in a semi-enclosed basin. The model represents a sea-ice covered and salinity stratified ocean, and consists of a sea-ice component and a two-layer ocean; a cold, fresh surface layer above a warmer, more saline layer. The sea-ice thickness depends on the atmospheric energy fluxes as well as the ocean heat flux. We introduce a thickness-dependent sea-ice export. Whether sea ice stabilizes or destabilizes against a freshwater perturbation is shown to depend on the representation of the vertical mixing. In a system where the vertical diffusivity is constant, the sea ice acts as a positive feedback on a freshwater perturbation. If the vertical diffusivity is derived from a constant mixing energy constraint, the sea ice acts as a negative feedback. However, both representations lead to a circulation that breaks down when the freshwater input at the surface is small. As a consequence, we get rapid changes in sea ice. In addition to low freshwater forcing, increasing deep-ocean temperatures promote instability and the disappearance of sea ice. Generally, the unstable state is reached before the vertical density difference disappears, and small changes in temperature and freshwater inputs can provoke abrupt changes in sea ice.

Experimental design for three interrelated marine ice sheet and ocean model intercomparison projects: MISMIP v. 3 (MISMIP +), ISOMIP v. 2 (ISOMIP +) and MISOMIP v. 1 (MISOMIP1)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Asay-Davis, Xylar S.; Cornford, Stephen L.; Durand, Gaël

Coupled ice sheet-ocean models capable of simulating moving grounding lines are just becoming available. Such models have a broad range of potential applications in studying the dynamics of marine ice sheets and tidewater glaciers, from process studies to future projections of ice mass loss and sea level rise. The Marine Ice Sheet-Ocean Model Intercomparison Project (MISOMIP) is a community effort aimed at designing and coordinating a series of model intercomparison projects (MIPs) for model evaluation in idealized setups, model verification based on observations, and future projections for key regions of the West Antarctic Ice Sheet (WAIS). Here we describe computationalmore » experiments constituting three interrelated MIPs for marine ice sheet models and regional ocean circulation models incorporating ice shelf cavities. These consist of ice sheet experiments under the Marine Ice Sheet MIP third phase (MISMIP+), ocean experiments under the Ice Shelf-Ocean MIP second phase (ISOMIP+) and coupled ice sheet-ocean experiments under the MISOMIP first phase (MISOMIP1). All three MIPs use a shared domain with idealized bedrock topography and forcing, allowing the coupled simulations (MISOMIP1) to be compared directly to the individual component simulations (MISMIP+ and ISOMIP+). The experiments, which have qualitative similarities to Pine Island Glacier Ice Shelf and the adjacent region of the Amundsen Sea, are designed to explore the effects of changes in ocean conditions, specifically the temperature at depth, on basal melting and ice dynamics. In future work, differences between model results will form the basis for the evaluation of the participating models.« less

Tidal Impacts on Oceanographic and Sea-ice Processes in the Southern Ocean

NASA Astrophysics Data System (ADS)

Padman, L.; Muench, R. D.; Howard, S.; Mueller, R.

2008-12-01

We review recent field and modeling results that demonstrate the importance of tides in establishing the oceanographic and sea-ice conditions in the boundary regions of the Southern Ocean. The tidal component dominates the total oceanic kinetic energy throughout much of the circum-Antarctic seas. This domination is especially pronounced over the continental slope and shelf including the sub-ice-shelf cavities. Tides provide most of the energy that forces diapycnal mixing under ice shelves and thereby contributes to basal melting. The resulting Ice Shelf Water is a significant component of the Antarctic Bottom Water (AABW) filling much of the deep global ocean. Tides exert significant divergent forcing on sea ice along glacial ice fronts and coastal regions, contributing to creation and maintenance of the coastal polynyas where much of the High Salinity Shelf Water component of AABW is formed. Additional tidally forced ice divergence along the shelf break and upper slope significantly impacts area-averaged ice growth and upper-ocean salinity. Tidally forced cross- slope advection, and mixing by the benthic stress associated with tidal currents along the shelf break and upper slope, strongly influence the paths, volume fluxes and hydrographic properties of benthic outflows of dense water leaving the continental shelf. These outflows provide primary source waters for the AABW. These results confirm that general ocean circulation and coupled ocean/ice/atmosphere climate models must incorporate the impacts of tides.

MIZEX. A Program for Mesoscale Air-Ice-Ocean Interaction Experiments in Arctic Marginal Ice Zones. III. Modeling the Marginal Ice Zone,

DTIC Science & Technology

1984-04-01

Ii TS C]r.I2 TAB 0] Unzanro’ unoed 0 justi fica ~r: 0 April 1984 vs - ASValabilitY Codes lvyall and/or U.S. Army Cold Regions Research and Engineering...coupled model. Fig. 1. Annual average simulated velocity fields. 3 192 Aloka 190 / 902 190+ WOO S’,. o Ice OnlY Mod" D"’, 55*w F~tth Yea’ Ice Ocean Model...A more precise delinga- inflow boundary conditions. 12 4- a. [ o ll ii traspert 00 0 0- 0e a I " i i , - - I I 1161 63 15 67 69 Ti 73 75 77 1980 *= 4h

Simulating Heinrich events in a coupled atmosphere-ocean-ice sheet model

NASA Astrophysics Data System (ADS)

Mikolajewicz, Uwe; Ziemen, Florian

2016-04-01

Heinrich events are among the most prominent events of long-term climate variability recorded in proxies across the northern hemisphere. They are the archetype of ice sheet - climate interactions on millennial time scales. Nevertheless, the exact mechanisms that cause Heinrich events are still under discussion, and their climatic consequences are far from being fully understood. We contribute to answering the open questions by studying Heinrich events in a coupled ice sheet model (ISM) atmosphere-ocean-vegetation general circulation model (AOVGCM) framework, where this variability occurs as part of the model generated internal variability without the need to prescribe external perturbations, as was the standard approach in almost all model studies so far. The setup consists of a northern hemisphere setup of the modified Parallel Ice Sheet Model (mPISM) coupled to the global coarse resolution AOVGCM ECHAM5/MPIOM/LPJ. The simulations used for this analysis were an ensemble covering substantial parts of the late Glacial forced with transient insolation and prescribed atmospheric greenhouse gas concentrations. The modeled Heinrich events show a marked influence of the ice discharge on the Atlantic circulation and heat transport, but none of the Heinrich events during the Glacial did show a complete collapse of the North Atlantic meridional overturning circulation. The simulated main consequences of the Heinrich events are a freshening and cooling over the North Atlantic and a drying over northern Europe.

Modelled ocean changes at the Plio-Pleistocene transition driven by Antarctic ice advance

PubMed Central

Hill, Daniel J.; Bolton, Kevin P.; Haywood, Alan M.

2017-01-01

The Earth underwent a major transition from the warm climates of the Pliocene to the Pleistocene ice ages between 3.2 and 2.6 million years ago. The intensification of Northern Hemisphere Glaciation is the most obvious result of the Plio-Pleistocene transition. However, recent data show that the ocean also underwent a significant change, with the convergence of deep water mass properties in the North Pacific and North Atlantic Ocean. Here we show that the lack of coastal ice in the Pacific sector of Antarctica leads to major reductions in Pacific Ocean overturning and the loss of the modern North Pacific Deep Water (NPDW) mass in climate models of the warmest periods of the Pliocene. These results potentially explain the convergence of global deep water mass properties at the Plio-Pleistocene transition, as Circumpolar Deep Water (CDW) became the common source. PMID:28252023

Sea-level response to abrupt ocean warming of Antarctic ice shelves

NASA Astrophysics Data System (ADS)

Pattyn, Frank

2016-04-01

Antarctica's contribution to global sea-level rise increases steadily. A fundamental question remains whether the ice discharge will lead to marine ice sheet instability (MISI) and collapse of certain sectors of the ice sheet or whether ice loss will increase linearly with the warming trends. Therefore, we employ a newly developed ice sheet model of the Antarctic ice sheet, called f.ETISh (fast Elementary Thermomechanical Ice Sheet model) to simulate ice sheet response to abrupt perturbations in ocean and atmospheric temperature. The f.ETISh model is a vertically integrated hybrid (SSA/SIA) ice sheet model including ice shelves. Although vertically integrated, thermomechanical coupling is ensured through a simplified representation of ice sheet thermodynamics based on an analytical solution of the vertical temperature profile, including strain heating and horizontal advection. The marine boundary is represented by a flux condition either coherent with power-law basal sliding (Pollard & Deconto (2012) based on Schoof (2007)) or according to Coulomb basal friction (Tsai et al., 2015), both taking into account ice-shelf buttressing. Model initialization is based on optimization of the basal friction field. Besides the traditional MISMIP tests, new tests with respect to MISI in plan-view models have been devised. The model is forced with stepwise ocean and atmosphere temperature perturbations. The former is based on a parametrised sub-shelf melt (limited to ice shelves), while the latter is based on present-day mass balance/surface temperature and corrected for elevation changes. Surface melting is introduced using a PDD model. Results show a general linear response in mass loss to ocean warming. Nonlinear response due to MISI occurs under specific conditions and is highly sensitive to the basal conditions near the grounding line, governed by both the initial conditions and the basal sliding/deformation model. The Coulomb friction model leads to significantly higher

Bathymetric and oceanic controls on Abbot Ice Shelf thickness and stability

NASA Astrophysics Data System (ADS)

Cochran, J. R.; Jacobs, S. S.; Tinto, K. J.; Bell, R. E.

2014-05-01

Ice shelves play key roles in stabilizing Antarctica's ice sheets, maintaining its high albedo and returning freshwater to the Southern Ocean. Improved data sets of ice shelf draft and underlying bathymetry are important for assessing ocean-ice interactions and modeling ice response to climate change. The long, narrow Abbot Ice Shelf south of Thurston Island produces a large volume of meltwater, but is close to being in overall mass balance. Here we invert NASA Operation IceBridge (OIB) airborne gravity data over the Abbot region to obtain sub-ice bathymetry, and combine OIB elevation and ice thickness measurements to estimate ice draft. A series of asymmetric fault-bounded basins formed during rifting of Zealandia from Antarctica underlie the Abbot Ice Shelf west of 94° W and the Cosgrove Ice Shelf to the south. Sub-ice water column depths along OIB flight lines are sufficiently deep to allow warm deep and thermocline waters observed near the western Abbot ice front to circulate through much of the ice shelf cavity. An average ice shelf draft of ~200 m, 15% less than the Bedmap2 compilation, coincides with the summer transition between the ocean surface mixed layer and upper thermocline. Thick ice streams feeding the Abbot cross relatively stable grounding lines and are rapidly thinned by the warmest inflow. While the ice shelf is presently in equilibrium, the overall correspondence between draft distribution and thermocline depth indicates sensitivity to changes in characteristics of the ocean surface and deep waters.

Observed platelet ice distributions in Antarctic sea ice: An index for ocean-ice shelf heat flux

NASA Astrophysics Data System (ADS)

Langhorne, P. J.; Hughes, K. G.; Gough, A. J.; Smith, I. J.; Williams, M. J. M.; Robinson, N. J.; Stevens, C. L.; Rack, W.; Price, D.; Leonard, G. H.; Mahoney, A. R.; Haas, C.; Haskell, T. G.

2015-07-01

Antarctic sea ice that has been affected by supercooled Ice Shelf Water (ISW) has a unique crystallographic structure and is called platelet ice. In this paper we synthesize platelet ice observations to construct a continent-wide map of the winter presence of ISW at the ocean surface. The observations demonstrate that, in some regions of coastal Antarctica, supercooled ISW drives a negative oceanic heat flux of -30 Wm-2 that persists for several months during winter, significantly affecting sea ice thickness. In other regions, particularly where the thinning of ice shelves is believed to be greatest, platelet ice is not observed. Our new data set includes the longest ice-ocean record for Antarctica, which dates back to 1902 near the McMurdo Ice Shelf. These historical data indicate that, over the past 100 years, any change in the volume of very cold surface outflow from this ice shelf is less than the uncertainties in the measurements.

Observations of Recent Arctic Sea Ice Volume Loss and Its Impact on Ocean-Atmosphere Energy Exchange and Ice Production

NASA Technical Reports Server (NTRS)

Kurtz, N. T.; Markus, T.; Farrell, S. L.; Worthen, D. L.; Boisvert, L. N.

2011-01-01

Using recently developed techniques we estimate snow and sea ice thickness distributions for the Arctic basin through the combination of freeboard data from the Ice, Cloud, and land Elevation Satellite (ICESat) and a snow depth model. These data are used with meteorological data and a thermodynamic sea ice model to calculate ocean-atmosphere heat exchange and ice volume production during the 2003-2008 fall and winter seasons. The calculated heat fluxes and ice growth rates are in agreement with previous observations over multiyear ice. In this study, we calculate heat fluxes and ice growth rates for the full distribution of ice thicknesses covering the Arctic basin and determine the impact of ice thickness change on the calculated values. Thinning of the sea ice is observed which greatly increases the 2005-2007 fall period ocean-atmosphere heat fluxes compared to those observed in 2003. Although there was also a decline in sea ice thickness for the winter periods, the winter time heat flux was found to be less impacted by the observed changes in ice thickness. A large increase in the net Arctic ocean-atmosphere heat output is also observed in the fall periods due to changes in the areal coverage of sea ice. The anomalously low sea ice coverage in 2007 led to a net ocean-atmosphere heat output approximately 3 times greater than was observed in previous years and suggests that sea ice losses are now playing a role in increasing surface air temperatures in the Arctic.

The interaction between sea ice and salinity-dominated ocean circulation: implications for halocline stability and rapid changes of sea ice cover

NASA Astrophysics Data System (ADS)

Jensen, Mari F.; Nilsson, Johan; Nisancioglu, Kerim H.

2016-11-01

Changes in the sea ice cover of the Nordic Seas have been proposed to play a key role for the dramatic temperature excursions associated with the Dansgaard-Oeschger events during the last glacial. In this study, we develop a simple conceptual model to examine how interactions between sea ice and oceanic heat and freshwater transports affect the stability of an upper-ocean halocline in a semi-enclosed basin. The model represents a sea ice covered and salinity stratified Nordic Seas, and consists of a sea ice component and a two-layer ocean. The sea ice thickness depends on the atmospheric energy fluxes as well as the ocean heat flux. We introduce a thickness-dependent sea ice export. Whether sea ice stabilizes or destabilizes against a freshwater perturbation is shown to depend on the representation of the diapycnal flow. In a system where the diapycnal flow increases with density differences, the sea ice acts as a positive feedback on a freshwater perturbation. If the diapycnal flow decreases with density differences, the sea ice acts as a negative feedback. However, both representations lead to a circulation that breaks down when the freshwater input at the surface is small. As a consequence, we get rapid changes in sea ice. In addition to low freshwater forcing, increasing deep-ocean temperatures promote instability and the disappearance of sea ice. Generally, the unstable state is reached before the vertical density difference disappears, and the temperature of the deep ocean do not need to increase as much as previously thought to provoke abrupt changes in sea ice.

Two new ways of mapping sea ice thickness using ocean waves

NASA Astrophysics Data System (ADS)

Wadhams, P.

2010-12-01

TWO NEW METHODS OF MAPPING SEA ICE THICKNESS USING OCEAN WAVES. P. Wadhams (1,2), Martin Doble (1,2) and F. Parmiggiani (3) (1) Dept. of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 0WA, UK. (2) Laboratoire d’Océanographie de Villefranche, Université Pierre et Marie Curie, 06234 Villefranche-sur-Mer, France (2) ISAC-CNR, Bologna, Italy Two new methods of mapping ice thickness have been recently developed and tested, both making use of the dispersion relation of ocean waves in ice of radically different types. In frazil-pancake ice, a young ice type in which cakes less than 5 m across float in a suspension of individual ice crystals, the propagation of waves has been successfully modelled by treating the ice layer as a highly viscous fluid. The model predicts a shortening of wavelengths within the ice. Two-dimensional Fourier analysis of successive SAR subscenes to track the directional spectrum of a wave field as it enters an ice edge shows that waves do indeed shorten within the ice, and the change has been successfully used to predict the thickness of the frazil-pancake layer. Concurrent shipborne sampling in the Antarctic has shown that the method is accurate, and we now propose its use throughout the important frazil-pancake regimes in the world ocean (Antarctic circumpolar ice edge zone, Greenland Sea, Bering Sea and others). A radically different type of dispersion occurs when ocean waves enter the continuous icefields of the central Arctic, when they couple with the elastic ice cover to propagate as a flexural-gravity wave. A two-axis tiltmeter array has been used to measure the resulting change in the dispersion relation for long ocean swell (15-30 s) originating from storms in the Greenland Sea. The dispersion relation is slightly different from swell in the open ocean, so if two such arrays are placed a substantial distance (100s of km) apart and used to observe the changing wave period of arrivals from a given

CICE, The Los Alamos Sea Ice Model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hunke, Elizabeth; Lipscomb, William; Jones, Philip

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 themore » 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.« less

Projecting Antarctic ice discharge using response functions from SeaRISE ice-sheet models

NASA Astrophysics Data System (ADS)

Levermann, A.; Winkelmann, R.; Nowicki, S.; Fastook, J. L.; Frieler, K.; Greve, R.; Hellmer, H. H.; Martin, M. A.; Meinshausen, M.; Mengel, M.; Payne, A. J.; Pollard, D.; Sato, T.; Timmermann, R.; Wang, W. L.; Bindschadler, R. A.

2014-08-01

The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we attempt to estimate the uncertainty range of future ice discharge from Antarctica by combining uncertainty in the climatic forcing, the oceanic response and the ice-sheet model response. The uncertainty in the global mean temperature increase is obtained from historically constrained emulations with the MAGICC-6.0 (Model for the Assessment of Greenhouse gas Induced Climate Change) model. The oceanic forcing is derived from scaling of the subsurface with the atmospheric warming from 19 comprehensive climate models of the Coupled Model Intercomparison Project (CMIP-5) and two ocean models from the EU-project Ice2Sea. The dynamic ice-sheet response is derived from linear response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. The resulting uncertainty range for the historic Antarctic contribution to global sea-level rise from 1992 to 2011 agrees with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time-delayed warming of the oceanic subsurface compared to the surface air temperature. The median of the additional ice loss for the 21st century is computed to 0.07 m (66% range: 0.02-0.14 m; 90% range: 0.0-0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 (Representative Concentration Pathway) scenario and 0.09 m (66% range: 0.04-0.21 m; 90% range: 0.01-0.37 m) for the strongest RCP-8.5. Assuming no time delay between the atmospheric warming and the oceanic subsurface, these

The transient response of ice-shelf melting to ocean change

NASA Astrophysics Data System (ADS)

Holland, P.

2017-12-01

Idealised modelling studies show that the melting of ice shelves varies as a quadratic function of ocean temperature. This means that warm-water ice shelves have higher melt rates and are also more sensitive to ocean warming. However, this result is the equilibrium response, derived from a set of ice—ocean simulations subjected to a fixed ocean forcing and run until steady. This study considers instead the transient response of melting, using unsteady simulations subjected to forcing conditions that are oscillated in time with a range of periods. The results show that when the ocean forcing is varied slowly, the melt rates follow the equililbrium response. However, for rapid ocean change melting deviates from the equilibrium response in interesting ways. The residence time of water in the sub-ice cavity offers a critical timescale. When the forcing varies slowly (period of oscillation >> residence time), the cavity is fully-flushed with forcing anomalies at all stages of the cycle and melting follows the equilibrium response. When the forcing varies rapidly (period ≤ residence time), multiple cold and warm anomalies coexist in the cavity, cancelling each other in the spatial mean and thus inducing a relatively steady melt rate. This implies that all ice shelves have a maximum frequency of ocean variability that can be manifested in melting. The results also show that ice shelves forced by warm water have high melt rates, high equilibrium sensitivity, and short residence times, hence a short timescale over which the equilibrium sensitivity is manifest. The most rapid melting adjustment is induced by warm anomalies that are also saline. Thus, ice shelves in the Amundsen and Bellingshausen seas, Antarctica, are highly sensitive to ocean change.

The ocean mixed layer under Southern Ocean sea-ice: Seasonal cycle and forcing

NASA Astrophysics Data System (ADS)

Pellichero, Violaine; Sallée, Jean-Baptiste; Schmidtko, Sunke; Roquet, Fabien; Charrassin, Jean-Benoît

2017-02-01

The oceanic mixed layer is the gateway for the exchanges between the atmosphere and the ocean; in this layer, all hydrographic ocean properties are set for months to millennia. A vast area of the Southern Ocean is seasonally capped by sea-ice, which alters the characteristics of the ocean mixed layer. The interaction between the ocean mixed layer and sea-ice plays a key role for water mass transformation, the carbon cycle, sea-ice dynamics, and ultimately for the climate as a whole. However, the structure and characteristics of the under-ice mixed layer are poorly understood due to the sparseness of in situ observations and measurements. In this study, we combine distinct sources of observations to overcome this lack in our understanding of the polar regions. Working with elephant seal-derived, ship-based, and Argo float observations, we describe the seasonal cycle of the ocean mixed-layer characteristics and stability of the ocean mixed layer over the Southern Ocean and specifically under sea-ice. Mixed-layer heat and freshwater budgets are used to investigate the main forcing mechanisms of the mixed-layer seasonal cycle. The seasonal variability of sea surface salinity and temperature are primarily driven by surface processes, dominated by sea-ice freshwater flux for the salt budget and by air-sea flux for the heat budget. Ekman advection, vertical diffusivity, and vertical entrainment play only secondary roles. Our results suggest that changes in regional sea-ice distribution and annual duration, as currently observed, widely affect the buoyancy budget of the underlying mixed layer, and impact large-scale water mass formation and transformation with far reaching consequences for ocean ventilation.

The Met Office Coupled Atmosphere/Land/Ocean/Sea-Ice Data Assimilation System

NASA Astrophysics Data System (ADS)

Lea, Daniel; Mirouze, Isabelle; King, Robert; Martin, Matthew; Hines, Adrian

2015-04-01

The Met Office has developed a weakly-coupled data assimilation (DA) system using the global coupled model HadGEM3 (Hadley Centre Global Environment Model, version 3). At present the analysis from separate ocean and atmosphere DA systems are combined to produced coupled forecasts. The aim of coupled DA is to produce a more consistent analysis for coupled forecasts which may lead to less initialisation shock and improved forecast performance. The HadGEM3 coupled model combines the atmospheric model UM (Unified Model) at 60 km horizontal resolution on 85 vertical levels, the ocean model NEMO (Nucleus for European Modelling of the Ocean) at 25 km (at the equator) horizontal resolution on 75 vertical levels, and the sea-ice model CICE at the same resolution as NEMO. The atmosphere and the ocean/sea-ice fields are coupled every 1-hour using the OASIS coupler. The coupled model is corrected using two separate 6-hour window data assimilation systems: a 4D-Var for the atmosphere with associated soil moisture content nudging and snow analysis schemes on the one hand, and a 3D-Var FGAT for the ocean and sea-ice on the other hand. The background information in the DA systems comes from a previous 6-hour forecast of the coupled model. To isolate the impact of the coupled DA, 13-month experiments have been carried out, including 1) a full atmosphere/land/ocean/sea-ice coupled DA run, 2) an atmosphere-only run forced by OSTIA SSTs and sea-ice with atmosphere and land DA, and 3) an ocean-only run forced by atmospheric fields from run 2 with ocean and sea-ice DA. In addition, 5-day and 10-day forecast runs, have been produced from initial conditions generated by either run 1 or a combination of runs 2 and 3. The different results have been compared to each other and, whenever possible, to other references such as the Met Office atmosphere and ocean operational analyses or the OSTIA SST data. The performance of the coupled DA is similar to the existing separate ocean and atmosphere

Southern Ocean frontal structure and sea-ice formation rates revealed by elephant seals

PubMed Central

Charrassin, J.-B.; Hindell, M.; Rintoul, S. R.; Roquet, F.; Sokolov, S.; Biuw, M.; Costa, D.; Boehme, L.; Lovell, P.; Coleman, R.; Timmermann, R.; Meijers, A.; Meredith, M.; Park, Y.-H.; Bailleul, F.; Goebel, M.; Tremblay, Y.; Bost, C.-A.; McMahon, C. R.; Field, I. C.; Fedak, M. A.; Guinet, C.

2008-01-01

Polar regions are particularly sensitive to climate change, with the potential for significant feedbacks between ocean circulation, sea ice, and the ocean carbon cycle. However, the difficulty in obtaining in situ data means that our ability to detect and interpret change is very limited, especially in the Southern Ocean, where the ocean beneath the sea ice remains almost entirely unobserved and the rate of sea-ice formation is poorly known. Here, we show that southern elephant seals (Mirounga leonina) equipped with oceanographic sensors can measure ocean structure and water mass changes in regions and seasons rarely observed with traditional oceanographic platforms. In particular, seals provided a 30-fold increase in hydrographic profiles from the sea-ice zone, allowing the major fronts to be mapped south of 60°S and sea-ice formation rates to be inferred from changes in upper ocean salinity. Sea-ice production rates peaked in early winter (April–May) during the rapid northward expansion of the pack ice and declined by a factor of 2 to 3 between May and August, in agreement with a three-dimensional coupled ocean–sea-ice model. By measuring the high-latitude ocean during winter, elephant seals fill a “blind spot” in our sampling coverage, enabling the establishment of a truly global ocean-observing system. PMID:18695241

The effects of ocean circulation on ocean-ice interaction and potential feedbacks in an idealized shelf cavity

NASA Astrophysics Data System (ADS)

Bishop, S. P.; Thompson, A. F.; Schodlok, M.

2016-02-01

The West Antarctic ice sheet is melting at unprecedented rates, which will impact global sea level rise. The ocean may be playing the dominant role in this ice melt through the upwelling of warm and salty Circumpolar Deep Water (CDW) in regions such as Pine Island Glacier (PIG). There is evidence that the Antarctic Slope Front at the continental shelf constrains shoreward transport of CDW by mesoscale eddies. However, little is known about the ocean-ice interaction and potential feedbacks that take place once this water is advected into ice shelf cavities. In this talk we use MITgcm to simulate an idealized setup of the PIG ice shelf cavity, similar to the setup in De Rydt et al. 2014, to understand the effects of ocean circulation and potential feedbacks of ice-shelf melt on the ocean circulation. To do this we run the model in two different configurations with and without a wind-driven current at the northern edge of the ice shelf and annually updating the geometry of the ice shelf based on the parameterized ice-shelf melt. Eddy heat and potential vorticity fluxes are diagnosed and presented for each of the simulations and compared with control simulations where the ice-shelf cavity is not modified. Results show high ice shelf melt during the first year with maximum values in excess of 60 meters near the grounding line, but settle to tens of meters during the following years.

Wave-Ice and Air-Ice-Ocean Interaction During the Chukchi Sea Ice Edge Advance

DTIC Science & Technology

2014-09-30

During cruise CU-B UAF UW Airborne expendable Ice Buoy (AXIB) Ahead, at and inside ice edge Surface meteorology T, SLP ~1 year CU-B UW...Balance (IMB) buoys Inside ice edge w/ >50cm thickness Ice mass balance T in snow-ice-ocean, T, SLP at surface ~1 year WHOI CRREL (SeaState DRI

Evolution of ocean-induced ice melt beneath Zachariæ Isstrøm, Northeast Greenland combining observations and an ocean general circulation model from 1978 to present

NASA Astrophysics Data System (ADS)

Cai, C.; Rignot, E. J.; Menemenlis, D.; Millan, R.; Bjørk, A. A.; Khan, S. A.; Charolais, A.

2017-12-01

Zachariæ Isstrøm, a major ice stream in northeast Greenland, lost a large fraction of its ice shelf during the last decade. We study the evolution of subaqueous melting of its floating section from 1978 to present. The ice shelf melt rate depends on thermal forcing from warm, salty, subsurface ocean waters of Atlantic origin (AW), the mixing of AW with fresh, buoyant subglacial discharge at the calving margin, and the shape of the sub-ice-shelf cavity. Subglacial discharge doubled as a result of enhanced ice sheet runoff caused by warmer air temperatures. Ocean thermal forcing has increased due to enhanced advection of AW. Using an Eulerian method, MEaSUREs ice velocity, Operation IceBridge (OIB) ice thickness, and RACMO2.3 surface balance data, we evaluate the ice shelf melt rate in 1978, 1999 and 2010. The melt rate doubled from 1999 to 2010. Using a Lagrangian method with World View imagery, we map the melt rate in detail from 2011 to 2016. We compare the results with 2D simulations from the Massachusetts Institute of Technology general circulation model (MITgcm), at a high spatial resolution (20-m horizontal and 40-m vertical grid spacing), using OIB ice thickness and sub-ice-shelf cavity for years 1978, 1996, 2010 and 2011, combined with in-situ ocean temperature/salinity data from Ocean Melting Greenland (OMG) 2017. We find that winter melt rates are 2 3 times smaller than summer rates and melt rates increase by one order magnitude during the transition from ice shelf termination to near-vertical calving wall termination. As the last remaining bits of floating ice shelf disappear, ice-ocean interaction will therefore play an increasing role in driving the glacier retreat into its marine-based basin. This work was performed under a contract with NASA Cryosphere Program at UC Irvine and Caltech's Jet Propulsion Laboratory.

Dynamic ocean topography from CryoSat-2: examining recent changes in ice-ocean stress and advancing a theory for Beaufort Gyre stabilization

NASA Astrophysics Data System (ADS)

Dewey, S.; Morison, J.; Kwok, R.; Dickinson, S.; Morison, D.; Andersen, R.

2017-12-01

Model and sparse observational evidence has shown the ocean current speed in the Beaufort Gyre to have increased and recently stabilized. However, full-basin altimetric observations of dynamic ocean topography (DOT) and ocean surface currents have yet to be applied to the dynamics of gyre stabilization. DOT fields from retracked CryoSat-2 retrievals in Arctic Ocean leads have enabled us to calculate 2-month average ocean geostrophic currents. These currents are crucial to accurately computing ice-ocean stress, especially because they have accelerated so that their speed rivals that of the overlying sea ice. Given these observations, we can shift our view of the Beaufort Gyre as a system in which the wind drives the ice and the ice drives a passive ocean to a system with the following feedback: After initial input of energy by wind, ice velocity decreases due to water drag and internal ice stress and the ocean drives the ice, reversing Ekman pumping and decelerating the gyre. This reversal changes the system from a persistently convergent regime to one in which freshwater is released from the gyre and doming of the gyre decreases, without any change in long-term average wind stress curl. Through these processes, the ice-ocean stress provides a key feedback in Beaufort Gyre stabilization.

Stress and deformation characteristics of sea ice in a high resolution numerical sea ice model.

NASA Astrophysics Data System (ADS)

Heorton, Harry; Feltham, Daniel; Tsamados, Michel

2017-04-01

The drift and deformation of sea ice floating on the polar oceans is due to the applied wind and ocean currents. The deformations of sea ice over ocean basin length scales have observable patterns; cracks and leads in satellite images and within the velocity fields generated from floe tracking. In a climate sea ice model the deformation of sea ice over ocean basin length scales is modelled using a rheology that represents the relationship between stresses and deformation within the sea ice cover. Here we investigate the link between observable deformation characteristics and the underlying internal sea ice stresses and force balance using the Los Alamos numerical sea ice climate model. In order to mimic laboratory experiments on the deformation of small cubes of sea ice we have developed an idealised square domain that tests the model response at spatial resolutions of up to 500m. We use the Elastic Anisotropic Plastic and Elastic Viscous Plastic rheologies, comparing their stability over varying resolutions and time scales. Sea ice within the domain is forced by idealised winds in order to compare the confinement of wind stresses and internal sea ice stresses. We document the characteristic deformation patterns of convergent, divergent and rotating stress states.

Sensitivity of Totten Glacier Ice Shelf extent and grounding line to oceanic forcing

NASA Astrophysics Data System (ADS)

Pelle, T.; Morlighem, M.; Choi, Y.

2017-12-01

Totten Glacier is a major outlet glacier of the East Antarctic Ice Sheet and has been shown to be vulnerable to ocean-induced melt in both its past and present states. The intrusion of warm, circumpolar deep water beneath the Totten Glacier Ice Shelf (TGIS) has been observed to accelerate ice shelf thinning and promote iceberg calving, a primary mechanism of mass discharge from Totten. As such, accurately simulating TGIS's ice front dynamics is crucial to the predictive capabilities of ice sheet models in this region. Here, we study the TGIS using the Ice Sheet System Model (ISSM) and test the applicability of three calving laws: Crevasse Formation calving, Eigen calving, and Tensile Stress calving. We simulate the evolution of Totten Glacier through 2100 under enhanced oceanic forcing in order to investigate both future changes in ice front dynamics and possible thresholds of instability. In addition, we artificially retreat Totten's ice front position and allow the model to proceed dynamically in order to analyze the response of the glacier to calving events. Our analyses show that Tensile Stress calving most accurately reproduces Totten Glacier's observed ice front position. Furthermore, unstable grounding line retreat is projected when Totten is simulated under stronger oceanic thermal forcing scenarios and when the calving front is significantly retreated.

Duality of Ross Ice Shelf systems: crustal boundary, ice sheet processes and ocean circulation from ROSETTA-Ice surveys

NASA Astrophysics Data System (ADS)

Tinto, K. J.; Siddoway, C. S.; Padman, L.; Fricker, H. A.; Das, I.; Porter, D. F.; Springer, S. R.; Siegfried, M. R.; Caratori Tontini, F.; Bell, R. E.

2017-12-01

Bathymetry beneath Antarctic ice shelves controls sub-ice-shelf ocean circulation and has a major influence on the stability and dynamics of the ice sheets. Beneath the Ross Ice Shelf, the sea-floor bathymetry is a product of both tectonics and glacial processes, and is influenced by the processes it controls. New aerogeophysical surveys have revealed a fundamental crustal boundary bisecting the Ross Ice Shelf and imparting a duality to the Ross Ice Shelf systems, encompassing bathymetry, ocean circulation and ice flow history. The ROSETTA-Ice surveys were designed to increase the resolution of Ross Ice Shelf mapping from the 55 km RIGGS survey of the 1970s to a 10 km survey grid, flown over three years from New York Air National Guard LC130s. Radar, LiDAR, gravity and magnetic instruments provide a top to bottom profile of the ice shelf and the underlying seafloor, with 20 km resolution achieved in the first two survey seasons (2015 and 2016). ALAMO ocean-profiling floats deployed in the 2016 season are measuring the temperature and salinity of water entering and exiting the sub-ice water cavity. A significant east-west contrast in the character of the magnetic and gravity fields reveals that the lithospheric boundary between East and West Antarctica exists not at the base of the Transantarctic Mountains (TAM), as previously thought, but 300 km further east. The newly-identified boundary spatially coincides with the southward extension of the Central High, a rib of shallow basement identified in the Ross Sea. The East Antarctic side is characterized by lower amplitude magnetic anomalies and denser TAM-type lithosphere compared to the West Antarctic side. The crustal structure imparts a fundamental duality on the overlying ice and ocean, with deeper bathymetry and thinner ice on the East Antarctic side creating a larger sub-ice cavity for ocean circulation. The West Antarctic side has a shallower seabed, more restricted ocean access and a more complex history of

The ocean mixed layer under Southern Ocean sea-ice: seasonal cycle and forcing.

NASA Astrophysics Data System (ADS)

Violaine, P.; Sallee, J. B.; Schmidtko, S.; Roquet, F.; Charrassin, J. B.

2016-02-01

The mixed-layer at the surface of the ocean is the gateway for all exchanges between air and sea. A vast area of the Southern Ocean is however seasonally capped by sea-ice, which alters this gateway and the characteristic the ocean mixed-layer. The interaction between the ocean mixed-layer and sea-ice plays a key role for water-mass formation and circulation, carbon cycle, sea-ice dynamics, and ultimately for the climate as a whole. However, the structure and characteristics of the mixed layer, as well as the processes responsible for its evolution, are poorly understood due to the lack of in-situ observations and measurements. We urgently need to better understand the forcing and the characteristics of the ocean mixed-layer under sea-ice if we are to understand and predict the world's climate. In this study, we combine a range of distinct sources of observation to overcome this lack in our understanding of the Polar Regions. Working on Elephant Seal-derived data as well as ship-based observations and Argo float data, we describe the seasonal cycle of the characteristics and stability of the ocean mixed layer over the entire Southern Ocean (South of 40°S), and specifically under sea-ice. Mixed-layer budgets of heat and freshwater are used to investigate the main forcings of the mixed-layer seasonal cycle. The seasonal variability of sea surface salinity and temperature are primarily driven by surface processes, dominated by sea-ice freshwater flux for the salt budget, and by air-sea flux for the heat budget. Ekman advection, vertical diffusivity and vertical entrainment play only secondary role.Our results suggest that changes in regional sea-ice distribution or sea-ice seasonal cycle duration, as currently observed, would widely affect the buoyancy budget of the underlying mixed-layer, and impacts large-scale water-mass formation and transformation.

Ecology of southern ocean pack ice.

PubMed

Brierley, Andrew S; Thomas, David N

2002-01-01

Around Antarctica the annual five-fold growth and decay of sea ice is the most prominent physical process and has a profound impact on marine life there. In winter the pack ice canopy extends to cover almost 20 million square kilometres--some 8% of the southern hemisphere and an area larger than the Antarctic continent itself (13.2 million square kilometres)--and is one of the largest, most dynamic ecosystems on earth. Biological activity is associated with all physical components of the sea-ice system: the sea-ice surface; the internal sea-ice matrix and brine channel system; the underside of sea ice and the waters in the vicinity of sea ice that are modified by the presence of sea ice. Microbial and microalgal communities proliferate on and within sea ice and are grazed by a wide range of proto- and macrozooplankton that inhabit the sea ice in large concentrations. Grazing organisms also exploit biogenic material released from the sea ice at ice break-up or melt. Although rates of primary production in the underlying water column are often low because of shading by sea-ice cover, sea ice itself forms a substratum that provides standing stocks of bacteria, algae and grazers significantly higher than those in ice-free areas. Decay of sea ice in summer releases particulate and dissolved organic matter to the water column, playing a major role in biogeochemical cycling as well as seeding water column phytoplankton blooms. Numerous zooplankton species graze sea-ice algae, benefiting additionally because the overlying sea-ice ceiling provides a refuge from surface predators. Sea ice is an important nursery habitat for Antarctic krill, the pivotal species in the Southern Ocean marine ecosystem. Some deep-water fish migrate to shallow depths beneath sea ice to exploit the elevated concentrations of some zooplankton there. The increased secondary production associated with pack ice and the sea-ice edge is exploited by many higher predators, with seals, seabirds and whales

Modelling sea ice dynamics

NASA Astrophysics Data System (ADS)

Murawski, Jens; Kleine, Eckhard

2017-04-01

Sea ice remains one of the frontiers of ocean modelling and is of vital importance for the correct forecasts of the northern oceans. At large scale, it is commonly considered a continuous medium whose dynamics is modelled in terms of continuum mechanics. Its specifics are a matter of constitutive behaviour which may be characterised as rigid-plastic. The new developed sea ice dynamic module bases on general principles and follows a systematic approach to the problem. Both drift field and stress field are modelled by a variational property. Rigidity is treated by Lagrangian relaxation. Thus one is led to a sensible numerical method. Modelling fast ice remains to be a challenge. It is understood that ridging and the formation of grounded ice keels plays a role in the process. The ice dynamic model includes a parameterisation of the stress associated with grounded ice keels. Shear against the grounded bottom contact might lead to plastic deformation and the loss of integrity. The numerical scheme involves a potentially large system of linear equations which is solved by pre-conditioned iteration. The entire algorithm consists of several components which result from decomposing the problem. The algorithm has been implemented and tested in practice.

Ganymede's internal structure including thermodynamics of magnesium sulfate oceans in contact with ice

NASA Astrophysics Data System (ADS)

Vance, Steve; Bouffard, Mathieu; Choukroun, Mathieu; Sotin, Christophe

2014-06-01

The large icy moons of Jupiter contain vast quantities of liquid water, a key ingredient for life. Ganymede and Callisto are weaker candidates for habitability than Europa, in part because of the model-based assumption that high-pressure ice layers cover their seafloors and prevent significant water-rock interaction. Water-rock interactions may occur, however, if heating at the rock-ice interface melts the high pressure ice. Highly saline fluids would be gravitationally stable, and might accumulate under the ice due to upward migration, refreezing, and fractionation of salt from less concentrated liquids. To assess the influence of salinity on Ganymede's internal structure, we use available phase-equilibrium data to calculate activity coefficients and predict the freezing of water ice in the presence of aqueous magnesium sulfate. We couple this new equation of state with thermal profiles in Ganymede's interior-employing recently published thermodynamic data for the aqueous phase-to estimate the thicknesses of layers of ice I, III, V, and VI. We compute core and silicate mantle radii consistent with available constraints on Ganymede's mass and gravitational moment of inertia. Mantle radii range from 800 to 900 km for the values of salt and heat flux considered here (4-44 mW m-2 and 0 to 10 wt% MgSO4). Ocean concentrations with salinity higher than 10 wt% have little high pressure ice. Even in a Ganymede ocean that is mostly liquid, achieving such high ocean salinity is permissible for the range of likely S/Si ratios. However, elevated salinity requires a smaller silicate mantle radius to satisfy mass and moment-of-inertia constraints, so ice VI is always present in Ganymede's ocean. For lower values of heat flux, oceans with salinity as low as 3 wt% can co-exist with ice III. Available experimental data indicate that ice phases III and VI become buoyant for salinity higher than 5 wt% and 10 wt%, respectively. Similar behavior probably occurs for ice V at salinities

Contributions of the atmosphere-land and ocean-sea ice model components to the tropical Atlantic SST bias in CESM1

NASA Astrophysics Data System (ADS)

Song, Z.; Lee, S. K.; Wang, C.; Kirtman, B. P.; Qiao, F.

2016-02-01

In order to identify and quantify intrinsic errors in the atmosphere-land and ocean-sea ice model components of the Community Earth System Model version 1 (CESM1) and their contributions to the tropical Atlantic sea surface temperature (SST) bias in CESM1, we propose a new method of diagnosis and apply it to a set of CESM1 simulations. Our analyses of the model simulations indicate that both the atmosphere-land and ocean-sea ice model components of CESM1 contain large errors in the tropical Atlantic. When the two model components are fully coupled, the intrinsic errors in the two components emerge quickly within a year with strong seasonality in their growth rates. In particular, the ocean-sea ice model contributes significantly in forcing the eastern equatorial Atlantic warm SST bias in early boreal summer. Further analysis shows that the upper thermocline water underneath the eastern equatorial Atlantic surface mixed layer is too warm in a stand-alone ocean-sea ice simulation of CESM1 forced with observed surface flux fields, suggesting that the mixed layer cooling associated with the entrainment of upper thermocline water is too weak in early boreal summer. Therefore, although we acknowledge the potential importance of the westerly wind bias in the western equatorial Atlantic and the low-level stratus cloud bias in the southeastern tropical Atlantic, both of which originate from the atmosphere-land model, we emphasize here that solving those problems in the atmosphere-land model alone does not resolve the equatorial Atlantic warm bias in CESM1.

Wave-Ice and Air-Ice-Ocean Interaction During the Chukchi Sea Ice Edge Advance

DTIC Science & Technology

2015-09-30

1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Wave -Ice and Air-Ice-Ocean Interaction During the...Chukchi Sea in the late summer have potentially changed the impact of fall storms by creating wave fields in the vicinity of the advancing ice edge. A...first) wave -ice interaction field experiment that adequately documents the relationship of a growing pancake ice cover with a time and space varying

Using Satellite-derived Ice Concentration to Represent Antarctic Coastal Polynyas in Ocean Climate Models

NASA Technical Reports Server (NTRS)

Stoessel, Achim; Markus, Thorsten

2003-01-01

The focus of this paper is on the representation of Antarctic coastal polynyas in global ice-ocean general circulation models (OGCMs), in particular their local, regional, and high-frequency behavior. This is verified with the aid of daily ice concentration derived from satellite passive microwave data using the NASATeam 2 (NT2) and the bootstrap (BS) algorithms. Large systematic regional and temporal discrepancies arise, some of which are related to the type of convection parameterization used in the model. An attempt is made to improve the fresh-water flux associated with melting and freezing in Antarctic coastal polynyas by ingesting (assimilating) satellite ice concentration where it comes to determining the thermodynamics of the open-water fraction of a model grid cell. Since the NT2 coastal open-water fraction (polynyas) tends to be less extensive than the simulated one in the decisive season and region, assimilating NT2 coastal ice concentration yields overall reduced net freezing rates, smaller formation rates of Antarctic Bottom Water, and a stronger southward flow of North Atlantic Deep Water across 30 S. Enhanced net freezing rates occur regionally when NT2 coastal ice concentration is assimilated, concomitant with a more realistic ice thickness distribution and accumulation of High-Salinity Shelf Water. Assimilating BS rather than NT2 coastal ice concentration, the differences to the non-assimilated simulation are generally smaller and of opposite sign. This suggests that the model reproduces coastal ice concentration in closer agreement with the BS data than with the NT2 data, while more realistic features emerge when NT2 data are assimilated.

Regular network model for the sea ice-albedo feedback in the Arctic.

PubMed

Müller-Stoffels, Marc; Wackerbauer, Renate

2011-03-01

The Arctic Ocean and sea ice form a feedback system that plays an important role in the global climate. The complexity of highly parameterized global circulation (climate) models makes it very difficult to assess feedback processes in climate without the concurrent use of simple models where the physics is understood. We introduce a two-dimensional energy-based regular network model to investigate feedback processes in an Arctic ice-ocean layer. The model includes the nonlinear aspect of the ice-water phase transition, a nonlinear diffusive energy transport within a heterogeneous ice-ocean lattice, and spatiotemporal atmospheric and oceanic forcing at the surfaces. First results for a horizontally homogeneous ice-ocean layer show bistability and related hysteresis between perennial ice and perennial open water for varying atmospheric heat influx. Seasonal ice cover exists as a transient phenomenon. We also find that ocean heat fluxes are more efficient than atmospheric heat fluxes to melt Arctic sea ice.

Collaborative Project. Understanding the effects of tides and eddies on the ocean dynamics, sea ice cover and decadal/centennial climate prediction using the Regional Arctic Climate Model (RACM)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hutchings, Jennifer; Joseph, Renu

2013-09-14

The goal of this project is to develop an eddy resolving ocean model (POP) with tides coupled to a sea ice model (CICE) within the Regional Arctic System Model (RASM) to investigate the importance of ocean tides and mesoscale eddies in arctic climate simulations and quantify biases associated with these processes and how their relative contribution may improve decadal to centennial arctic climate predictions. Ocean, sea ice and coupled arctic climate response to these small scale processes will be evaluated with regard to their influence on mass, momentum and property exchange between oceans, shelf-basin, ice-ocean, and ocean-atmosphere. The project willmore » facilitate the future routine inclusion of polar tides and eddies in Earth System Models when computing power allows. As such, the proposed research addresses the science in support of the BER’s Climate and Environmental Sciences Division Long Term Measure as it will improve the ocean and sea ice model components as well as the fully coupled RASM and Community Earth System Model (CESM) and it will make them more accurate and computationally efficient.« less

Atmosphere-Ice-Ocean-Ecosystem Processes in a Thinner Arctic Sea Ice Regime: The Norwegian Young Sea ICE (N-ICE2015) Expedition

NASA Astrophysics Data System (ADS)

Granskog, Mats A.; Fer, Ilker; Rinke, Annette; Steen, Harald

2018-03-01

Arctic sea ice has been in rapid decline the last decade and the Norwegian young sea ICE (N-ICE2015) expedition sought to investigate key processes in a thin Arctic sea ice regime, with emphasis on atmosphere-snow-ice-ocean dynamics and sea ice associated ecosystem. The main findings from a half-year long campaign are collected into this special section spanning the Journal of Geophysical Research: Atmospheres, Journal of Geophysical Research: Oceans, and Journal of Geophysical Research: Biogeosciences and provide a basis for a better understanding of processes in a thin sea ice regime in the high Arctic. All data from the campaign are made freely available to the research community.

How does ice sheet loading affect ocean flow around Antarctica?

NASA Astrophysics Data System (ADS)

Dijkstra, H. A.; Rugenstein, M. A.; Stocchi, P.; von der Heydt, A. S.

2012-12-01

Interactions and dynamical feedbacks between ocean circulation, heat and atmospheric moisture transport, ice sheet evolution, and Glacial Isostatic Adjustment (GIA) are overlooked issues in paleoclimatology. Here we will present first results on how ocean flows were possibly affected by the glaciation of Antarctica across the Eocene-Oligocene Transition (~ 34 Ma) through GIA and bathymetry variations. GIA-induced gravitationally self-consistent bathymetry variations are determined by solving the Sea Level Equation (SLE), which describes the time dependent shape of (i) the solid Earth and (ii) the equipotential surface of gravity. Since the ocean circulation equations are defined relative to the equipotential surface of gravity, only bathymetry variations can influence ocean flows, although the sea surface slope will also change through time due to gravitational attraction. We use the Hallberg Isopycnal Model under late Eocene conditions to calculate equilibrium ocean flows in a domain in which the bathymetry evolves under ice loading according to the SLE. The bathymetric effects of the glaciation of Antarctica lead to substantial spatial changes in ocean flows, and close to the coast, the flow even reverses direction. Volume transports through the Drake Passage and Tasman Seaway adjust to the new bathymetry. The results indicate that GIA-induced ocean flow variations alone may have had an impact on sedimentation and erosion patterns, the repositioning of fronts, ocean heat transport and grounding line and ice sheet stability.

Ocean Disposal of Man-Made Ice Piers

EPA Pesticide Factsheets

The National Science Foundation is permitted to ocean dump man-made ice piers from its base at McMurdo Sound in Antarctica under a MPRSA general permit. Information is provided about ice piers and impacts of ice pier disposal.

Year-Round Pack Ice in the Weddell Sea, Antarctica: Response and Sensitivity to Atmospheric and Oceanic Forcing

NASA Technical Reports Server (NTRS)

Geiger, Cathleen A.; Ackley, Stephen F.; Hibler, William D., III

1997-01-01

Using a dynamic-thermodynamic numerical sea-ice model, external oceanic and atmospheric forcings on sea ice in the Weddell Sea are examined to identify physical processes associated with the seasonal cycle of pack ice, and to identify further the parameters that coupled models need to consider in predicting the response of the pack ice to climate and ocean-circulation changes. In agreement with earlier studies, the primary influence on the winter ice-edge maximum extent is air temperature. Ocean heat flux has more impact on the minimum-ice-edge extent and in reducing pack-ice thickness, especially in the eastern-Weddell Sea. Low relative humidity enhances ice growth in thin ice and open-water regions, producing a more realistic ice edge along the coastal areas of the western-Weddell Sea where dry continental air has an impact. The modeled extent of the Weddell summer pack is equally sensitive to ocean heat flux and atmospheric relative humidity variations with the more dynamic responses being from the atmosphere. Since the atmospheric regime in the eastern Weddell is dominated by marine intrusions from lower latitudes, with high humidity already, it is unlikely that either the moisture trans- port could be further raised or that it could be significantly lowered because of its distance from the continent (the lower humidity source). Ocean heat-transport variability is shown to lead to overall ice thinning in the model response and is a known feature of the actual system, as evidenced by the occurrence of the Weddell Polynya in the mid 1970s.

The sensitivity of the Greenland Ice Sheet to glacial-interglacial oceanic forcing

NASA Astrophysics Data System (ADS)

Tabone, Ilaria; Blasco, Javier; Robinson, Alexander; Alvarez-Solas, Jorge; Montoya, Marisa

2018-04-01

Observations suggest that during the last decades the Greenland Ice Sheet (GrIS) has experienced a gradually accelerating mass loss, in part due to the observed speed-up of several of Greenland's marine-terminating glaciers. Recent studies directly attribute this to warming North Atlantic temperatures, which have triggered melting of the outlet glaciers of the GrIS, grounding-line retreat and enhanced ice discharge into the ocean, contributing to an acceleration of sea-level rise. Reconstructions suggest that the influence of the ocean has been of primary importance in the past as well. This was the case not only in interglacial periods, when warmer climates led to a rapid retreat of the GrIS to land above sea level, but also in glacial periods, when the GrIS expanded as far as the continental shelf break and was thus more directly exposed to oceanic changes. However, the GrIS response to palaeo-oceanic variations has yet to be investigated in detail from a mechanistic modelling perspective. In this work, the evolution of the GrIS over the past two glacial cycles is studied using a three-dimensional hybrid ice-sheet-shelf model. We assess the effect of the variation of oceanic temperatures on the GrIS evolution on glacial-interglacial timescales through changes in submarine melting. The results show a very high sensitivity of the GrIS to changing oceanic conditions. Oceanic forcing is found to be a primary driver of GrIS expansion in glacial times and of retreat in interglacial periods. If switched off, palaeo-atmospheric variations alone are not able to yield a reliable glacial configuration of the GrIS. This work therefore suggests that considering the ocean as an active forcing should become standard practice in palaeo-ice-sheet modelling.

Operational coupled atmosphere - ocean - ice forecast system for the Gulf of St. Lawrence, Canada

NASA Astrophysics Data System (ADS)

Faucher, M.; Roy, F.; Desjardins, S.; Fogarty, C.; Pellerin, P.; Ritchie, H.; Denis, B.

2009-09-01

A fully interactive coupled atmosphere-ocean-ice forecasting system for the Gulf of St. Lawrence (GSL) has been running in experimental mode at the Canadian Meteorological Centre (CMC) for the last two winter seasons. The goal of this project is to provide more accurate weather and sea ice forecasts over the GSL and adjacent coastal areas by including atmosphere-oceanice interactions in the CMC operational forecast system using a formal coupling strategy between two independent modeling components. The atmospheric component is the Canadian operational GEM model (Côté et al. 1998) and the oceanic component is the ocean-ice model for the Gulf of St. Lawrence developed at the Maurice Lamontagne Institute (IML) (Saucier et al. 2003, 2004). The coupling between those two models is achieved by exchanging surface fluxes and variables through MPI communication. The re-gridding of the variables is done with a package developed at the Recherche en Prevision Numerique centre (RPN, Canada). Coupled atmosphere - ocean - ice forecasts are issued once a day based on 00GMT data. Results for the past two years have demonstrated that the coupled system produces improved forecasts in and around the GSL during all seasons, proving that atmosphere-ocean-ice interactions are indeed important even for short-term Canadian weather forecasts. This has important implications for other coupled modeling and data assimilation partnerships that are in progress involving EC, the Department of Fisheries and Oceans (DFO) and the National Defense (DND). Following this experimental phase, it is anticipated that this GSL system will be the first fully interactive coupled system to be implemented at CMC.

Landfast Sea Ice Breakouts: Stabilizing Ice Features, Oceanic and Atmospheric Forcing at Barrow, Alaska

NASA Astrophysics Data System (ADS)

Jones, J.; Eicken, H.; Mahoney, A. R.; MV, R.; Kambhamettu, C.; Fukamachi, Y.; Ohshima, K. I.; George, C.

2016-12-01

Landfast sea ice is an important seasonal feature along most Arctic coastlines, such as that of the Chukchi Sea near Barrow, Alaska. Its stability throughout the ice season is determined by many factors but grounded pressure ridges are the primary stabilizing component. Landfast ice breakouts occur when these grounded ridges fail or unground, and previously stationary ice detaches from the coast and drifts away. Using ground-based radar imagery from a coastal ice and ocean observatory at Barrow, we have developed a method to estimate the extent of grounded ridges by tracking ice motion and deformation over the course of winter and have derived ice keel depth and potential for grounding from cumulative convergent ice motion. Estimates of landfast ice grounding strength have been compared to the atmospheric and oceanic stresses acting on the landfast ice before and during breakout events to determine prevailing causes for the failure of stabilizing features. Applying this approach to two case studies in 2008 and 2010, we conclude that a combination of atmospheric and oceanic stresses may have caused the breakouts analyzed in this study, with the latter as the dominant force. Preconditioning (as weakening) of grounded ridges by sea level variations may facilitate failure of the ice sheet leading to breakout events.

Landfast sea ice breakouts: Stabilizing ice features, oceanic and atmospheric forcing at Barrow, Alaska

NASA Astrophysics Data System (ADS)

Jones, Joshua; Eicken, Hajo; Mahoney, Andrew; MV, Rohith; Kambhamettu, Chandra; Fukamachi, Yasushi; Ohshima, Kay I.; George, J. Craig

2016-09-01

Landfast sea ice is an important seasonal feature along most Arctic coastlines, such as that of the Chukchi Sea near Barrow, Alaska. Its stability throughout the ice season is determined by many factors but grounded pressure ridges are the primary stabilizing component. Landfast ice breakouts occur when these grounded ridges fail or unground, and previously stationary ice detaches from the coast and drifts away. Using ground-based radar imagery from a coastal ice and ocean observatory at Barrow, we have developed a method to estimate the extent of grounded ridges by tracking ice motion and deformation over the course of winter and have derived ice keel depth and potential for grounding from cumulative convergent ice motion. Estimates of landfast ice grounding strength have been compared to the atmospheric and oceanic stresses acting on the landfast ice before and during breakout events to determine prevailing causes for the failure of stabilizing features. Applying this approach to two case studies in 2008 and 2010, we conclude that a combination of atmospheric and oceanic stresses may have caused the breakouts analyzed in this study, with the latter as the dominant force. Preconditioning (as weakening) of grounded ridges by sea level variations may facilitate failure of the ice sheet leading to breakout events.

Exploring the sensitivity of global ocean circulation to future ice loss from Antarctica

DOE Office of Scientific and Technical Information (OSTI.GOV)

Condron, Alan

The sensitivity of the global ocean circulation and climate to large increases in iceberg calving and meltwater discharges from the Antarctic Ice Sheet (AIS) are rarely studied and poorly understood. The requirement to investigate this topic is heightened by growing evidence that the West Antarctic Ice Sheet (WAIS) is vulnerable to rapid retreat and collapse on multidecadal-to-centennial timescales. Observations collected over the last 30 years indicate that the WAIS is now losing mass at an accelerated and that a collapse may have already begun in the Amundsen Sea sector. In addition, some recent future model simulations of the AIS showmore » the potential for rapid ice sheet retreat in the next 50 – 300 years. Such a collapse would be associated with the discharge of enormous volumes of ice and meltwater to the Southern Ocean. This project funds PI Condron to begin assessing the sensitivity of the global ocean circulation to projected increases in meltwater discharge and iceberg calving from the AIS for the next 50 – 100 years. A series of climate model simulations will determine changes in ocean circulation and temperature at the ice sheet grounding line, the role of mesoscale ocean eddies in mixing and transporting freshwater away from the continent to deep water formation regions, and the likely impact on the northward transport of heat to Europe and North America.« less

Modeling the basal melting and marine ice accretion of the Amery Ice Shelf

NASA Astrophysics Data System (ADS)

Galton-Fenzi, B. K.; Hunter, J. R.; Coleman, R.; Marsland, S. J.; Warner, R. C.

2012-09-01

The basal mass balance of the Amery Ice Shelf (AIS) in East Antarctica is investigated using a numerical ocean model. The main improvements of this model over previous studies are the inclusion of frazil formation and dynamics, tides and the use of the latest estimate of the sub-ice-shelf cavity geometry. The model produces a net basal melt rate of 45.6 Gt year-1 (0.74 m ice year-1) which is in good agreement with reviewed observations. The melting at the base of the ice shelf is primarily due to interaction with High Salinity Shelf Water created from the surface sea-ice formation in winter. The temperature difference between the coldest waters created in the open ocean and the in situ freezing point of ocean water in contact with the deepest part of the AIS drives a melt rate that can exceed 30 m of ice year-1. The inclusion of frazil dynamics is shown to be important for both melting and marine ice accretion (refreezing). Frazil initially forms in the supercooled water layer adjacent to the base of the ice shelf. The net accretion of marine ice is 5.3 Gt year-1, comprised of 3.7 Gt year-1 of frazil accretion and 1.6 Gt year-1 of direct basal refreezing.

Ocean stratification reduces melt rates at the grounding zone of the Ross Ice Shelf

NASA Astrophysics Data System (ADS)

Begeman, C. B.; Tulaczyk, S. M.; Marsh, O.; Mikucki, J.; Stanton, T. P.; Hodson, T. O.; Siegfried, M. R.; Powell, R. D.; Christianson, K. A.; King, M. A.

2017-12-01

Ocean-driven melting of ice shelves is often invoked as the primary mechanism for triggering ice loss from Antarctica. However, due to the difficulty in accessing the sub-ice-shelf ocean cavity, the relationship between ice-shelf melt rates and ocean conditions is poorly understood, particularly near the transition from grounded to floating ice, known as the grounding zone. Here we present the first borehole oceanographic observations from the grounding zone of Antarctica's largest ice shelf. Contrary to predictions that tidal currents near grounding zones should mix the water column, driving high ice-shelf melt rates, we find a stratified sub-ice-shelf water column. The vertical salinity gradient dominates stratification over a weakly unstable vertical temperature gradient; thus, stratification takes the form of a double-diffusive staircase. These conditions limit vertical heat fluxes and lead to low melt rates in the ice-shelf grounding zone. While modern grounding zone melt rates may presently be overestimated in models that assume efficient tidal mixing, the high sensitivity of double-diffusive staircases to ocean freshening and warming suggests future melt rates may be underestimated, biasing projections of global sea-level rise.

The frequency and extent of sub-ice phytoplankton blooms in the Arctic Ocean

PubMed Central

Horvat, Christopher; Jones, David Rees; Iams, Sarah; Schroeder, David; Flocco, Daniela; Feltham, Daniel

2017-01-01

In July 2011, the observation of a massive phytoplankton bloom underneath a sea ice–covered region of the Chukchi Sea shifted the scientific consensus that regions of the Arctic Ocean covered by sea ice were inhospitable to photosynthetic life. Although the impact of widespread phytoplankton blooms under sea ice on Arctic Ocean ecology and carbon fixation is potentially marked, the prevalence of these events in the modern Arctic and in the recent past is, to date, unknown. We investigate the timing, frequency, and evolution of these events over the past 30 years. Although sea ice strongly attenuates solar radiation, it has thinned significantly over the past 30 years. The thinner summertime Arctic sea ice is increasingly covered in melt ponds, which permit more light penetration than bare or snow-covered ice. Our model results indicate that the recent thinning of Arctic sea ice is the main cause of a marked increase in the prevalence of light conditions conducive to sub-ice blooms. We find that as little as 20 years ago, the conditions required for sub-ice blooms may have been uncommon, but their frequency has increased to the point that nearly 30% of the ice-covered Arctic Ocean in July permits sub-ice blooms. Recent climate change may have markedly altered the ecology of the Arctic Ocean. PMID:28435859

Sea Ice Biogeochemistry: A Guide for Modellers

PubMed Central

Tedesco, Letizia; Vichi, Marcello

2014-01-01

Sea ice is a fundamental component of the climate system and plays a key role in polar trophic food webs. Nonetheless sea ice biogeochemical dynamics at large temporal and spatial scales are still rarely described. Numerical models may potentially contribute integrating among sparse observations, but available models of sea ice biogeochemistry are still scarce, whether their relevance for properly describing the current and future state of the polar oceans has been recently addressed. A general methodology to develop a sea ice biogeochemical model is presented, deriving it from an existing validated model application by extension of generic pelagic biogeochemistry model parameterizations. The described methodology is flexible and considers different levels of ecosystem complexity and vertical representation, while adopting a strategy of coupling that ensures mass conservation. We show how to apply this methodology step by step by building an intermediate complexity model from a published realistic application and applying it to analyze theoretically a typical season of first-year sea ice in the Arctic, the one currently needing the most urgent understanding. The aim is to (1) introduce sea ice biogeochemistry and address its relevance to ocean modelers of polar regions, supporting them in adding a new sea ice component to their modelling framework for a more adequate representation of the sea ice-covered ocean ecosystem as a whole, and (2) extend our knowledge on the relevant controlling factors of sea ice algal production, showing that beyond the light and nutrient availability, the duration of the sea ice season may play a key-role shaping the algal production during the on going and upcoming projected changes. PMID:24586604

Using Laboratory Experiments to Improve Ice-Ocean Parameterizations

NASA Astrophysics Data System (ADS)

McConnochie, C. D.; Kerr, R. C.

2017-12-01

Numerical models of ice-ocean interactions are typically unable to resolve the transport of heat and salt to the ice face. Instead, models rely upon parameterizations that have not been sufficiently validated by observations. Recent laboratory experiments of ice-saltwater interactions allow us to test the standard parameterization of heat and salt transport to ice faces - the three-equation model. The three-equation model predicts that the melt rate is proportional to the fluid velocity while the experimental results typically show that the melt rate is independent of the fluid velocity. By considering an analysis of the boundary layer that forms next to a melting ice face, we suggest a resolution to this disagreement. We show that the three-equation model makes the implicit assumption that the thickness of the diffusive sublayer next to the ice is set by a shear instability. However, at low flow velocities, the sublayer is instead set by a convective instability. This distinction leads to a threshold velocity of approximately 4 cm/s at geophysically relevant conditions, above which the form of the parameterization should be valid. In contrast, at flow speeds below 4 cm/s, the three-equation model will underestimate the melt rate. By incorporating such a minimum velocity into the three-equation model, predictions made by numerical simulations could be easily improved.

Sea Ice on the Southern Ocean

NASA Technical Reports Server (NTRS)

Jacobs, Stanley S.

1998-01-01

Year-round satellite records of sea ice distribution now extend over more than two decades, providing a valuable tool to investigate related characteristics and circulations in the Southern Ocean. We have studied a variety of features indicative of oceanic and atmospheric interactions with Antarctic sea ice. In the Amundsen & Bellingshausen Seas, sea ice extent was found to have decreased by approximately 20% from 1973 through the early 1990's. This change coincided with and probably contributed to recently warmer surface conditions on the west side of the Antarctic Peninsula, where air temperatures have increased by approximately 0.5 C/decade since the mid-1940's. The sea ice decline included multiyear cycles of several years in length superimposed on high interannual variability. The retreat was strongest in summer, and would have lowered the regional mean ice thickness, with attendant impacts upon vertical heat flux and the formation of snow ice and brine. The cause of the regional warming and loss of sea ice is believed to be linked to large-scale circulation changes in the atmosphere and ocean. At the eastern end of the Weddell Gyre, the Cosmonaut Polyna revealed greater activity since 1986, a recurrence pattern during recent winters and two possible modes of formation. Persistence in polynya location was noted off Cape Ann, where the coastal current can interact more strongly with the Antarctic Circumpolar Current. As a result of vorticity conservation, locally enhanced upwelling brings warmer deep water into the mixed layer, causing divergence and melting. In the Ross Sea, ice extent fluctuates over periods of several years, with summer minima and winter maxima roughly in phase. This leads to large interannual cycles of sea ice range, which correlate positively with meridinal winds, regional air temperatures and subsequent shelf water salinities. Deep shelf waters display considerable interannual variability, but have freshened by approximately 0.03/decade

Determination of a Critical Sea Ice Thickness Threshold for the Central Arctic Ocean

NASA Astrophysics Data System (ADS)

Ford, V.; Frauenfeld, O. W.; Nowotarski, C. J.

2017-12-01

While sea ice extent is readily measurable from satellite observations and can be used to assess the overall survivability of the Arctic sea ice pack, determining the spatial variability of sea ice thickness remains a challenge. Turbulent and conductive heat fluxes are extremely sensitive to ice thickness but are dominated by the sensible heat flux, with energy exchange expected to increase with thinner ice cover. Fluxes over open water are strongest and have the greatest influence on the atmosphere, while fluxes over thick sea ice are minimal as heat conduction from the ocean through thick ice cannot reach the atmosphere. We know that turbulent energy fluxes are strongest over open ocean, but is there a "critical thickness of ice" where fluxes are considered non-negligible? Through polar-optimized Weather Research and Forecasting model simulations, this study assesses how the wintertime Arctic surface boundary layer, via sensible heat flux exchange and surface air temperature, responds to sea ice thinning. The region immediately north of Franz Josef Land is characterized by a thickness gradient where sea ice transitions from the thickest multi-year ice to the very thin marginal ice seas. This provides an ideal location to simulate how the diminishing Arctic sea ice interacts with a warming atmosphere. Scenarios include both fixed sea surface temperature domains for idealized thickness variability, and fixed ice fields to detect changes in the ocean-ice-atmosphere energy exchange. Results indicate that a critical thickness threshold exists below 1 meter. The threshold is between 0.4-1 meters thinner than the critical thickness for melt season survival - the difference between first year and multi-year ice. Turbulent heat fluxes and surface air temperature increase as sea ice thickness transitions from perennial ice to seasonal ice. While models predict a sea ice free Arctic at the end of the warm season in future decades, sea ice will continue to transform

Ice, Ocean and Atmosphere Interactions in the Arctic Marginal Ice Zone

DTIC Science & Technology

2015-09-30

the northward retreat of the ice edge. Through the long-term measurement of the key oceanic, atmospheric, and sea ice processes that...began to move southward towards the Alaskan coast. In 2104 the anomalous areas of ice retreat were the region north of Alaska...and Siberia. (see figures below). This is not uncommon as these regions have seen the greatest retreat in sea ice. See http://nsidc.org

Arctic sea ice in the global eddy-permitting ocean reanalysis ORAP5

NASA Astrophysics Data System (ADS)

Tietsche, Steffen; Balmaseda, Magdalena A.; Zuo, Hao; Mogensen, Kristian

2017-08-01

We discuss the state of Arctic sea ice in the global eddy-permitting ocean reanalysis Ocean ReAnalysis Pilot 5 (ORAP5). Among other innovations, ORAP5 now assimilates observations of sea ice concentration using a univariate 3DVar-FGAT scheme. We focus on the period 1993-2012 and emphasize the evaluation of model performance with respect to recent observations of sea ice thickness. We find that sea ice concentration in ORAP5 is close to assimilated observations, with root mean square analysis residuals of less than 5 % in most regions. However, larger discrepancies exist for the Labrador Sea and east of Greenland during winter owing to biases in the free-running model. Sea ice thickness is evaluated against three different observational data sets that have sufficient spatial and temporal coverage: ICESat, IceBridge and SMOSIce. Large-scale features like the gradient between the thickest ice in the Canadian Arctic and thinner ice in the Siberian Arctic are simulated well by ORAP5. However, some biases remain. Of special note is the model's tendency to accumulate too thick ice in the Beaufort Gyre. The root mean square error of ORAP5 sea ice thickness with respect to ICESat observations is 1.0 m, which is on par with the well-established PIOMAS model sea ice reconstruction. Interannual variability and trend of sea ice volume in ORAP5 also compare well with PIOMAS and ICESat estimates. We conclude that, notwithstanding a relatively simple sea ice data assimilation scheme, the overall state of Arctic sea ice in ORAP5 is in good agreement with observations and will provide useful initial conditions for predictions.

Nudging the Arctic Ocean to quantify Arctic sea ice feedbacks

NASA Astrophysics Data System (ADS)

Dekker, Evelien; Severijns, Camiel; Bintanja, Richard

2017-04-01

It is well-established that the Arctic is warming 2 to 3 time faster than rest of the planet. One of the great uncertainties in climate research is related to what extent sea ice feedbacks amplify this (seasonally varying) Arctic warming. Earlier studies have analyzed existing climate model output using correlations and energy budget considerations in order to quantify sea ice feedbacks through indirect methods. From these analyses it is regularly inferred that sea ice likely plays an important role, but details remain obscure. Here we will take a different and a more direct approach: we will keep the sea ice constant in a sensitivity simulation, using a state-of -the-art climate model (EC-Earth), applying a technique that has never been attempted before. This experimental technique involves nudging the temperature and salinity of the ocean surface (and possibly some layers below to maintain the vertical structure and mixing) to a predefined prescribed state. When strongly nudged to existing (seasonally-varying) sea surface temperatures, ocean salinity and temperature, we force the sea ice to remain in regions/seasons where it is located in the prescribed state, despite the changing climate. Once we obtain fixed' sea ice, we will run a future scenario, for instance 2 x CO2 with and without prescribed sea ice, with the difference between these runs providing a measure as to what extent sea ice contributes to Arctic warming, including the seasonal and geographical imprint of the effects.

Modelling the Climate - Greenland Ice Sheet Interaction in the Coupled Ice-sheet/Climate Model EC-EARTH - PISM

NASA Astrophysics Data System (ADS)

Yang, S.; Madsen, M. S.; Rodehacke, C. B.; Svendsen, S. H.; Adalgeirsdottir, G.

2014-12-01

Recent observations show that the Greenland ice sheet (GrIS) has been losing mass with an increasing speed during the past decades. Predicting the GrIS changes and their climate consequences relies on the understanding of the interaction of the GrIS with the climate system on both global and local scales, and requires climate model systems with an explicit and physically consistent ice sheet module. A fully coupled global climate model with a dynamical ice sheet model for the GrIS has recently been developed. The model system, EC-EARTH - PISM, consists of the EC-EARTH, an atmosphere, ocean and sea ice model system, and the Parallel Ice Sheet Model (PISM). The coupling of PISM includes a modified surface physical parameterization in EC-EARTH adapted to the land ice surface over glaciated regions in Greenland. The PISM ice sheet model is forced with the surface mass balance (SMB) directly computed inside the EC-EARTH atmospheric module and accounting for the precipitation, the surface evaporation, and the melting of snow and ice over land ice. PISM returns the simulated basal melt, ice discharge and ice cover (extent and thickness) as boundary conditions to EC-EARTH. This coupled system is mass and energy conserving without being constrained by any anomaly correction or flux adjustment, and hence is suitable for investigation of ice sheet - climate feedbacks. Three multi-century experiments for warm climate scenarios under (1) the RCP85 climate forcing, (2) an abrupt 4xCO2 and (3) an idealized 1% per year CO2 increase are performed using the coupled model system. The experiments are compared with their counterparts of the standard CMIP5 simulations (without the interactive ice sheet) to evaluate the performance of the coupled system and to quantify the GrIS feedbacks. In particular, the evolution of the Greenland ice sheet under the warm climate and its impacts on the climate system are investigated. Freshwater fluxes from the Greenland ice sheet melt to the Arctic

The CONCEPTS Global Ice-Ocean Prediction System: Establishing an Environmental Prediction Capability in Canada

NASA Astrophysics Data System (ADS)

Pellerin, Pierre; Smith, Gregory; Testut, Charles-Emmanuel; Surcel Colan, Dorina; Roy, Francois; Reszka, Mateusz; Dupont, Frederic; Lemieux, Jean-Francois; Beaudoin, Christiane; He, Zhongjie; Belanger, Jean-Marc; Deacu, Daniel; Lu, Yimin; Buehner, Mark; Davidson, Fraser; Ritchie, Harold; Lu, Youyu; Drevillon, Marie; Tranchant, Benoit; Garric, Gilles

2015-04-01

Here we describe a new system implemented recently at the Canadian Meteorological Centre (CMC) entitled the Global Ice Ocean Prediction System (GIOPS). GIOPS provides ice and ocean analyses and 10 day forecasts daily at 00GMT on a global 1/4° resolution grid. GIOPS includes a full multivariate ocean data assimilation system that combines satellite observations of sea level anomaly and sea surface temperature (SST) together with in situ observations of temperature and salinity. In situ observations are obtained from a variety of sources including: the Argo network of autonomous profiling floats, moorings, ships of opportunity, marine mammals and research cruises. Ocean analyses are blended with sea ice analyses produced by the Global Ice Analysis System.. GIOPS has been developed as part of the Canadian Operational Network of Coupled Environmental PredicTion Systems (CONCEPTS) tri-departmental initiative between Environment Canada, Fisheries and Oceans Canada and National Defense. The development of GIOPS was made through a partnership with Mercator-Océan, a French operational oceanography group. Mercator-Océan provided the ocean data assimilation code and assistance with the system implementation. GIOPS has undergone a rigorous evaluation of the analysis, trial and forecast fields demonstrating its capacity to provide high-quality products in a robust and reliable framework. In particular, SST and ice concentration forecasts demonstrate a clear benefit with respect to persistence. These results support the use of GIOPS products within other CMC operational systems, and more generally, as part of a Government of Canada marine core service. Impact of a two-way coupling between the GEM atmospheric model and NEMO-CICE ocean-ice model will also be presented.

Ocean Profile Measurements During the Seasonal Ice Zone Reconnaissance Surveys Ocean Profiles

DTIC Science & Technology

2017-01-01

repeated ocean, ice, and atmospheric measurements across the Beaufort-Chukchi sea seasonal sea ice zone (SIZ) utilizing US Coast Guard Arctic Domain...contributing to the rapid decline in summer ice extent that has occurred in recent years. The SIZ is the region between maximum winter sea ice extent and...minimum summer sea ice extent. As such, it contains the full range of positions of the marginal ice zone (MIZ) where sea ice interacts with open water

The effect of sea ice on the solar energy budget in the astmosphere-sea ice-ocean system: A model study

NASA Technical Reports Server (NTRS)

Jin, Z.; Stamnes, Knut; Weeks, W. F.; Tsay, Si-Chee

1994-01-01

A coupled one-dimensional multilayer and multistream radiative transfer model has been developed and applied to the study of radiative interactions in the atmosphere, sea ice, and ocean system. The consistent solution of the radiative transfer equation in this coupled system automatically takes into account the refraction and reflection at the air-ice interface and allows flexibility in choice of stream numbers. The solar radiation spectrum (0.25 micron-4.0 micron) is divided into 24 spectral bands to account adequately for gaseous absorption in the atmosphere. The effects of ice property changes, including salinity and density variations, as well as of melt ponds and snow cover variations over the ice on the solar energy distribution in the entire system have been studied quantitatively. The results show that for bare ice it is the scattering, determined by air bubbles and brine pockets, in just a few centimeters of the top layer of ice that plays the most important role in the solar energy absorption and partitioning in the entire system. Ice thickness is important to the energy distribution only when the ice is thin, while the absorption in the atmosphere is not sensitive to ice thickness exceeds about 70 cm. The presence of clouds moderates all the sensitivities of the absorptive amounts in each layer to the variations in the ice properties and ice thickness. Comparisons with observational spectral albedo values for two simple ice types are also presented.

Ice-Ocean Interactions to the North-West of Greenland: Glaciers, Straits, Ice Bridges, and the Rossby Radius (Invited)

NASA Astrophysics Data System (ADS)

Muenchow, A.; Falkner, K. K.; Melling, H.; Johnson, H. L.; Huntley, H. S.; Ryan, P.; Friends Of Petermann

2010-12-01

Petermann Glacier at 81 N latitude is a major outlet glacier adjacent to Nares Strait. It terminates in a long (70 km), narrow (16 km) and thin (50 m) floating tongue and has a grounding line more than 500 m below sea level. A calving event in 2010 reduced the floating area by 25% and produced a single 240 km2 ice island currently moving south in Nares Strait where it will likely interact with island to potentially create a temporary polynya in Nares Strait. The 2010 calving from Petermann Glacier contributes <10% to its mass balance as more than 80% is lost due to basal melting by the ocean. Hence the largely unexplored physics at the ice-ocean interface determine how a changing climate impacts this outlet glacier. Conducting exploratory surveys inside Petermann Fjord in 2003, 2007, and 2009, we find a 1100 m deep fjord connected to Nares Strait via a sill at 350-450 m depth. The fjord receives about 3 times the amount of heat required for the basal melt rates. Furthermore, limited data and analytical modeling suggests a 3-dimensional circulation over the upper 300-m of the water column with a coastally trapped buoyant outflow. We integrate these findings with more complete oceanic time series data from an array moored in Nares Strait from 2003 through 2009 near 80.5 N. In the past Nares Strait and Petermann Fjord were covered by land fast sea ice during the 9-10 month long winter season. Archeological and remotely sensed records indicate that an ice bridge formed regularly at the southern end of Nares Strait creating the North-Water polynya near 79 N latitude. Since 2006 this ice bridge has largely failed to form, leading, perhaps, to the occasional formation of a secondary ice bridge 300 km to the north where Nares Strait connects to the Arctic Ocean. However, this ice bridge appears to form for shorter periods only. Consequently Arctic sea ice can now exit the Arctic in winter via pathways to the west of Greenland all year. We speculate that this changed ocean

Possible Sea Ice Impacts on Oceanic Deep Convection

NASA Technical Reports Server (NTRS)

Parkinson, C. L.

1984-01-01

Many regions of the world ocean known or suspected to have deep convection are sea-ice covered for at least a portion of the annual cycle. As this suggests that sea ice might have some impact on generating or maintaining this phenomenon, several mechanisms by which sea ice could exert an influence are presented in the following paragraphs. Sea ice formation could be a direct causal factor in deep convection by providing the surface density increase necessary to initiate the convective overturning. As sea ice forms, either by ice accretion or by in situ ice formation in open water or in lead areas between ice floes, salt is rejected to the underlying water. This increases the water salinity, thereby increasing water density in the mixed layer under the ice. A sufficient increase in density will lead to mixing with deeper waters, and perhaps to deep convection or even bottom water formation. Observations are needed to establish whether this process is actually occurring; it is most likely in regions with extensive ice formation and a relatively unstable oceanic density structure.

Slush Fund: The Multiphase Nature of Oceanic Ices and Its Role in Shaping Europa's Icy Shell

NASA Astrophysics Data System (ADS)

Buffo, J.; Schmidt, B. E.; Huber, C.

2017-12-01

The role of Europa's ice shell in mediating ocean-surface interaction, constraining potential habitability of the underlying hydrosphere, and dictating the surface morphology of the moon is discussed extensively in the literature, yet the dynamics and characteristics of the shell itself remain largely unconstrained. Some of the largest unknowns arise from underrepresented physics and varying a priori assumptions built into the current ice shell models. Here we modify and apply a validated one-dimensional reactive transport model designed to simulate the formation and evolution of terrestrial sea ice to the Europa environment. The top-down freezing of sea ice due to conductive heat loss to the atmosphere is akin to the formation of the Jovian moon's outer ice shell, albeit on a different temporal and spatial scale. Nevertheless, the microscale physics that govern the formation of sea ice on Earth (heterogenous solidification leading to brine pockets and channels, multiphase reactive transport phenomena, gravity drainage) likely operate in a similar manner at the ice-ocean interface of Europa, dictating the thermal, chemical, and mechanical properties of the ice shell. Simulations of the European ice-ocean interface at different stages during the ice shell's evolution are interpolated to produce vertical profiles of temperature, salinity, solid fraction, and eutectic points throughout the entire shell. Additionally, the model is coupled to the equilibrium chemistry package FREZCHEM to investigate the impact a diverse range of putative European ocean chemistries has on ice shell properties. This method removes the need for a priori assumptions of impurity entrainment rates and ice shell properties, thus providing a first principles constraint on the stratigraphic characteristics of a simulated European ice shell. These insights have the potential to improve existing estimates for the onset of solid state convection, melt lens formation due to eutectic melting, ice

A two-dimensional ocean model for long-term climatic simulations: Stability and coupling to atmospheric and sea ice models

NASA Astrophysics Data System (ADS)

Harvey, L. D. Danny

1992-06-01

A two-dimensional (latitude-depth) deep ocean model is presented which is coupled to a sea ice model and an Energy Balance Climate Model (EBCM), the latter having land-sea and surface-air resolution. The processes which occur in the ocean model are thermohaline overturning driven by the horizontal density gradient, shallow wind-driven overturning cells, convective overturning, and vertical and horizontal diffusion of heat and salt. The density field is determined from the temperature and salinity fields using a nonlinear equation of state. Mixed layer salinity is affected by evaporation, precipitation, runoff from continents, and sea ice freezing and melting, as well as by advective, convective, and diffusive exchanges with the deep ocean. The ocean model is first tested in an uncoupled mode, in which hemispherically symmetric mixed layer temperature and salinity, or salinity flux, are specified as upper boundary conditions. An experiment performed with previous models is repeated in which a mixed layer salinity perturbation is introduced in the polar half of one hemisphere after switching from a fixed salinity to a fixed salinity flux boundary condition. For small values of the vertical diffusion coefficient KV, the model undergoes self-sustained oscillations with a period of about 1500 years. With larger values of KV, the model locks into either an asymmetric mode with a single overturning cell spanning both hemispheres, or a symmetric quiescent state with downwelling near the equator, upwelling at high latitudes, and a warm deep ocean (depending on the value of KV). When the ocean model is forced with observed mixed layer temperature and salinity, no oscillations occur. The model successfully simulates the very weak meridional overturning and strong Antarctic Circumpolar Current at the latitudes of the Drake Passage. The coupled EBCM-deep ocean model displays internal oscillations with a period of 3000 years if the ocean fraction is uniform with latitude and KV

A global, high-resolution data set of ice sheet topography, cavity geometry, and ocean bathymetry

NASA Astrophysics Data System (ADS)

Schaffer, Janin; Timmermann, Ralph; Arndt, Jan Erik; Savstrup Kristensen, Steen; Mayer, Christoph; Morlighem, Mathieu; Steinhage, Daniel

2016-10-01

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies, and global surface height on a spherical grid with now 30 arcsec grid spacing. For this new data set, called RTopo-2, we used the General Bathymetric Chart of the Oceans (GEBCO_2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the International Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves, and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry surrounding the Greenland continent. We modified data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ, and Sermilik Fjord, assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off Northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79° N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centres of Technical University of Denmark (DTU), Operation Icebridge (NASA

Implications of sea-ice biogeochemistry for oceanic production and emissions of dimethyl sulfide in the Arctic

NASA Astrophysics Data System (ADS)

Hayashida, Hakase; Steiner, Nadja; Monahan, Adam; Galindo, Virginie; Lizotte, Martine; Levasseur, Maurice

2017-06-01

Sea ice represents an additional oceanic source of the climatically active gas dimethyl sulfide (DMS) for the Arctic atmosphere. To what extent this source contributes to the dynamics of summertime Arctic clouds is, however, not known due to scarcity of field measurements. In this study, we developed a coupled sea ice-ocean ecosystem-sulfur cycle model to investigate the potential impact of bottom-ice DMS and its precursor dimethylsulfoniopropionate (DMSP) on the oceanic production and emissions of DMS in the Arctic. The results of the 1-D model simulation were compared with field data collected during May and June of 2010 in Resolute Passage. Our results reproduced the accumulation of DMS and DMSP in the bottom ice during the development of an ice algal bloom. The release of these sulfur species took place predominantly during the earlier phase of the melt period, resulting in an increase of DMS and DMSP in the underlying water column prior to the onset of an under-ice phytoplankton bloom. Production and removal rates of processes considered in the model are analyzed to identify the processes dominating the budgets of DMS and DMSP both in the bottom ice and the underlying water column. When openings in the ice were taken into account, the simulated sea-air DMS flux during the melt period was dominated by episodic spikes of up to 8.1 µmol m-2 d-1. Further model simulations were conducted to assess the effects of the incorporation of sea-ice biogeochemistry on DMS production and emissions, as well as the sensitivity of our results to changes of uncertain model parameters of the sea-ice sulfur cycle. The results highlight the importance of taking into account both the sea-ice sulfur cycle and ecosystem in the flux estimates of oceanic DMS near the ice margins and identify key uncertainties in processes and rates that should be better constrained by new observations.

Enceladus's ice shell thickness and ocean depth from gravity, topography, and libration measurements

NASA Astrophysics Data System (ADS)

Trinh, A.; Rivoldini, A.; Beuthe, M.; Rekier, J.; Baland, R. M.; Van Hoolst, T.

2017-12-01

One of Cassini's major achievements is the discovery of a global ocean a few kilometres beneath Enceladus's south polar terrain. Here we infer the thickness of Enceladus's ice shell and ocean from Cassini's observations using our latest models of isostatic compensation, shell libration, and ocean dynamics.

The emergence of modern sea ice cover in the Arctic Ocean.

PubMed

Knies, Jochen; Cabedo-Sanz, Patricia; Belt, Simon T; Baranwal, Soma; Fietz, Susanne; Rosell-Melé, Antoni

2014-11-28

Arctic sea ice coverage is shrinking in response to global climate change and summer ice-free conditions in the Arctic Ocean are predicted by the end of the century. The validity of this prediction could potentially be tested through the reconstruction of the climate of the Pliocene epoch (5.33-2.58 million years ago), an analogue of a future warmer Earth. Here we show that, in the Eurasian sector of the Arctic Ocean, ice-free conditions prevailed in the early Pliocene until sea ice expanded from the central Arctic Ocean for the first time ca. 4 million years ago. Amplified by a rise in topography in several regions of the Arctic and enhanced freshening of the Arctic Ocean, sea ice expanded progressively in response to positive ice-albedo feedback mechanisms. Sea ice reached its modern winter maximum extension for the first time during the culmination of the Northern Hemisphere glaciation, ca. 2.6 million years ago.

Arctic Ocean sea ice cover during the penultimate glacial and the last interglacial.

PubMed

Stein, Ruediger; Fahl, Kirsten; Gierz, Paul; Niessen, Frank; Lohmann, Gerrit

2017-08-29

Coinciding with global warming, Arctic sea ice has rapidly decreased during the last four decades and climate scenarios suggest that sea ice may completely disappear during summer within the next about 50-100 years. Here we produce Arctic sea ice biomarker proxy records for the penultimate glacial (Marine Isotope Stage 6) and the subsequent last interglacial (Marine Isotope Stage 5e). The latter is a time interval when the high latitudes were significantly warmer than today. We document that even under such warmer climate conditions, sea ice existed in the central Arctic Ocean during summer, whereas sea ice was significantly reduced along the Barents Sea continental margin influenced by Atlantic Water inflow. Our proxy reconstruction of the last interglacial sea ice cover is supported by climate simulations, although some proxy data/model inconsistencies still exist. During late Marine Isotope Stage 6, polynya-type conditions occurred off the major ice sheets along the northern Barents and East Siberian continental margins, contradicting a giant Marine Isotope Stage 6 ice shelf that covered the entire Arctic Ocean.Coinciding with global warming, Arctic sea ice has rapidly decreased during the last four decades. Here, using biomarker records, the authors show that permanent sea ice was still present in the central Arctic Ocean during the last interglacial, when high latitudes were warmer than present.

Oceans Melting Greenland (OMG): 2017 Observations and the First Look at Yearly Ocean/Ice Changes

NASA Astrophysics Data System (ADS)

Willis, J. K.; Rignot, E. J.; Fenty, I. G.; Khazendar, A.; Moller, D.; Tinto, K. J.; Morison, J.; Schodlok, M.; Thompson, A. F.; Fukumori, I.; Holland, D.; Forsberg, R.; Jakobsson, M.; Dinardo, S. J.

2017-12-01

Oceans Melting Greenland (OMG) is an airborne NASA Mission to investigate the role of the oceans in ice loss around the margins of the Greenland Ice Sheet. A five-year campaign, OMG will directly measure ocean warming and glacier retreat around all of Greenland. By relating these two, we will explore one of the most pressing open questions about how climate change drives sea level rise: How quickly are the warming oceans melting the Greenland Ice Sheet from the edges? This year, OMG collected its second set of both elevation maps of marine terminating glaciers and ocean temperature and salinity profiles around all of Greenland. This give us our first look at year-to-year changes in both ice volume at the margins, as well as the volume and extent of warm, salty Atlantic water present on the continental shelf. In addition, we will compare recent data in east Greenland waters with historical ocean observations that suggest a long-term warming trend there. Finally, we will briefly review the multi-beam sonar and airborne gravity campaigns—both of which were completed last year—and the dramatic improvement they had on bathymetry maps over the continental shelf around Greenland.

Fast Response of the Tropics to an Abrupt Loss of Arctic Sea Ice via Ocean Dynamics

NASA Astrophysics Data System (ADS)

Wang, Kun; Deser, Clara; Sun, Lantao; Tomas, Robert A.

2018-05-01

The role of ocean dynamics in the transient adjustment of the coupled climate system to an abrupt loss of Arctic sea ice is investigated using experiments with Community Climate System Model version 4 in two configurations: a thermodynamic slab mixed layer ocean and a full-depth ocean that includes both dynamics and thermodynamics. Ocean dynamics produce a distinct sea surface temperature warming maximum in the eastern equatorial Pacific, accompanied by an equatorward intensification of the Intertropical Convergence Zone and Hadley Circulation. These tropical responses are established within 25 years of ice loss and contrast markedly with the quasi-steady antisymmetric coupled response in the slab-ocean configuration. A heat budget analysis reveals the importance of anomalous vertical advection tied to a monotonic temperature increase below 200 m for the equatorial sea surface temperature warming maximum in the fully coupled model. Ocean dynamics also rapidly modify the midlatitude atmospheric response to sea ice loss.

ARCTIC SEA ICE EXTENT AND DRIFT, MODELED AS A VISCOUS FLUID.

USGS Publications Warehouse

Ling, Chi-Hai; Parkinson, Claire L.

1986-01-01

A dynamic/thermodynamic numerical model of sea ice has been used to calculate the yearly cycle of sea ice thicknesses, concentrations, and velocities in the Arctic Ocean and surrounding seas. The model combines the formulations of two previous models, taking the thermodynamics and momentum equations from the model of Parkinson and Washington and adding the constitutive equation and equation of state from the model of Ling, Rasmussen, and Campbell. Simulated annually averaged ice drift vectors compare well with observed ice drift from the Arctic Ocean Buoy Program.

Ice-free summers predominant in the late Miocene central Arctic Ocean - New insights from a proxy-modeling approach

NASA Astrophysics Data System (ADS)

Stein, Ruediger; Fahl, Kirsten; Schreck, Michael; Knorr, Gregor; Forwick, Matthias; Lohmann, Gerrit; Niessen, Frank

2016-04-01

During Polarstern Expedition PS87/2014, we discovered multiple submarine landslides over a distance of >350 km along Lomonosov Ridge between about 81°N and 84°N (Stein, 2015). The load and erosional behaviour of an extended ice sheet/shelf that probably occurred during major Quaternary glaciations, may have caused physical conditions that triggered these landslides and major down-slope transport of sediments at this part of Lomonosov Ridge (Stein et al., 2016 and further references therein). The removal of younger sediments from steep headwalls has led to exhumation of Miocene to early Quaternary sediments close to the seafloor, allowing the retrieval of such old sediments by gravity coring and multi-proxy studies of theses sediments. Within one of these studies (Stein et al., 2016), we used for the first time the sea-ice biomarker IP25 (for background of approach see Belt et al., 2007; Müller et al., 2009, 2011) together with alkenone-based sea-surface temperatures (SST) to reconstruct upper Miocene Arctic Ocean sea-ice and SST conditions. The presence of IP25 as proxy for spring sea-ice cover and alkenone-based relatively warm summer SST of >4 °C support a seasonal sea-ice cover with an ice-free summer season being dominant during (most of) the late Miocene central Arctic Ocean. A comparison of our proxy data with Miocene climate simulations seems to favour either relatively high late Miocene atmospheric CO2 concentrations and/or an overly weak sensitivity of the model to simulate the magnitude of high-latitude warming in a warmer than modern climate. References: Belt, S.T., Massé, G., Rowland, S.J., Poulin, M., Michel, and C., LeBlanc, B., 2007. A novel chemical fossil of palaeo sea ice: IP25, Organic Geochemistry 38, 16-27. Müller, J., Massé, G., Stein, R., and Belt, S., 2009. Extreme variations in sea ice cover for Fram Strait during the past 30 ka. Nature Geoscience, DOI: 10.1038/NGEO665. Müller, J., Wagner, A., Fahl, K., Stein, R., Prange, M., and

Insights Into Ice-Ocean Interactions on Earth and Europa

NASA Astrophysics Data System (ADS)

Lawrence, J.; Schmidt, B. E.; Winslow, L.; Doran, P. T.; Kim, S.; Walker, C. C.; Buffo, J.; Skidmore, M. L.; Soderlund, K. M.; Blankenship, D. D.; Bramall, N. E.; Johnson, A.; Rack, F. R.; Stone, W.; Kimball, P.; Clark, E.

2016-12-01

Europa and Earth appear to be drastically different worlds, yet below their icy crusts the two likely share similar oceanic conditions including temperatures, pressures (relatively), and salinity. Earth's ice shelves provide an important analog for the physiochemical, and potentially microbial, characteristics of icy worlds. NASA's ASTEP program funded Sub-Ice Marine and PLanetary-analog Ecosystems (SIMPLE) to help address the fundamental processes occurring at ice ocean interfaces, the extent and limitations of life in sub-ice environments, and how environmental properties and biological communities interact. The relationships between currents, temperature, and salinity with physical processes such as melt, freeze, and marine ice accretion at the basal surfaces of ice shelves influence habitability yet are poorly understood even on Earth. Resultant processes such as the inclusion of ocean-derived material in ice shelves and the transport of biotics from the interface towards the surface via ablation, convection, and diapirism also have important astrobiological implications for Europa.Here, we present results from CTD and imaging data gathered at multiple locations beneath the McMurdo Ice Shelf (MIS) to highlight how the ice and ocean interact in a Europan analog environment. Over the course of three years, the SIMPLE team observed heterogeneity in the water column and basal ice beneath the MIS. During the recent 2015 field season we deployed ARTEMIS, an AUV capable of characterizing the interface over multiple kilometer missions, and conducted daily CTD casts to 480 m (bottom depth 529 m) in November adjacent to the terminus of the MIS to capture temporal variation in the water column. These casts show the presence of transient water masses related to the tidal period, each containing a single or double temperature minimum (down to -1.97 °C from -1.93 °C) between 60 to 150 m depth. Further comparisons between years and sampling locations demonstrate the

Response of a comprehensive climate model to a broad range of external forcings: relevance for deep ocean ventilation and the development of late Cenozoic ice ages

NASA Astrophysics Data System (ADS)

Galbraith, Eric; de Lavergne, Casimir

2018-03-01

Over the past few million years, the Earth descended from the relatively warm and stable climate of the Pliocene into the increasingly dramatic ice age cycles of the Pleistocene. The influences of orbital forcing and atmospheric CO2 on land-based ice sheets have long been considered as the key drivers of the ice ages, but less attention has been paid to their direct influences on the circulation of the deep ocean. Here we provide a broad view on the influences of CO2, orbital forcing and ice sheet size according to a comprehensive Earth system model, by integrating the model to equilibrium under 40 different combinations of the three external forcings. We find that the volume contribution of Antarctic (AABW) vs. North Atlantic (NADW) waters to the deep ocean varies widely among the simulations, and can be predicted from the difference between the surface densities at AABW and NADW deep water formation sites. Minima of both the AABW-NADW density difference and the AABW volume occur near interglacial CO2 (270-400 ppm). At low CO2, abundant formation and northward export of sea ice in the Southern Ocean contributes to very salty and dense Antarctic waters that dominate the global deep ocean. Furthermore, when the Earth is cold, low obliquity (i.e. a reduced tilt of Earth's rotational axis) enhances the Antarctic water volume by expanding sea ice further. At high CO2, AABW dominance is favoured due to relatively warm subpolar North Atlantic waters, with more dependence on precession. Meanwhile, a large Laurentide ice sheet steers atmospheric circulation as to strengthen the Atlantic Meridional Overturning Circulation, but cools the Southern Ocean remotely, enhancing Antarctic sea ice export and leading to very salty and expanded AABW. Together, these results suggest that a `sweet spot' of low CO2, low obliquity and relatively small ice sheets would have poised the AMOC for interruption, promoting Dansgaard-Oeschger-type abrupt change. The deep ocean temperature and

An Assessment of Southern Ocean Water Masses and Sea Ice During 1988-2007 in a Suite of Interannual CORE-II Simulations

NASA Technical Reports Server (NTRS)

Downes, Stephanie M.; Farneti, Riccardo; Uotila, Petteri; Griffies, Stephen M.; Marsland, Simon J.; Bailey, David; Behrens, Erik; Bentsen, Mats; Bi, Daohua; Biastoch, Arne;

The influence of meridional ice transport on Europa's ocean stratification and heat content

NASA Astrophysics Data System (ADS)

Zhu, Peiyun; Manucharyan, Georgy E.; Thompson, Andrew F.; Goodman, Jason C.; Vance, Steven D.

2017-06-01

Jupiter's moon Europa likely hosts a saltwater ocean beneath its icy surface. Geothermal heating and rotating convection in the ocean may drive a global overturning circulation that redistributes heat vertically and meridionally, preferentially warming the ice shell at the equator. Here we assess the previously unconstrained influence of ocean-ice coupling on Europa's ocean stratification and heat transport. We demonstrate that a relatively fresh layer can form at the ice-ocean interface due to a meridional ice transport forced by the differential ice shell heating between the equator and the poles. We provide analytical and numerical solutions for the layer's characteristics, highlighting their sensitivity to critical ocean parameters. For a weakly turbulent and highly saline ocean, a strong buoyancy gradient at the base of the freshwater layer can suppress vertical tracer exchange with the deeper ocean. As a result, the freshwater layer permits relatively warm deep ocean temperatures.

The influence of meridional ice transport on Europa's ocean stratification and heat content

NASA Astrophysics Data System (ADS)

Zhu, P.; Manucharyan, G.; Thompson, A. F.; Goodman, J. C.; Vance, S.

2017-12-01

Jupiter's moon Europa likely hosts a saltwater ocean beneath its icy surface. Geothermal heating and rotating convection in the ocean may drive a global overturning circulation that redistributes heat vertically and meridionally, preferentially warming the ice shell at the equator. Here we assess thepreviously unconstrained influence of ocean-ice coupling on Europa's ocean stratification and heat transport. We demonstrate that a relatively fresh layer can form at the ice-ocean interface due to a meridional ice transport forced by the differential ice shell heating between the equator and the poles. We provide analytical and numerical solutions for the layer's characteristics, highlighting their sensitivity to critical ocean parameters. For a weakly turbulent and highly saline ocean, a strong buoyancy gradient at the base of the freshwater layer can suppress vertical tracer exchange with the deeper ocean. As a result, the freshwater layer permits relatively warm deep ocean temperatures.

Modelling and parameterizing the influence of tides on ice-shelf melt rates

NASA Astrophysics Data System (ADS)

Jourdain, N.; Molines, J. M.; Le Sommer, J.; Mathiot, P.; de Lavergne, C.; Gurvan, M.; Durand, G.

2017-12-01

Significant Antarctic ice sheet thinning is observed in several sectors of Antarctica, in particular in the Amundsen Sea sector, where warm circumpolar deep waters affect basal melting. The later has the potential to trigger marine ice sheet instabilities, with an associated potential for rapid sea level rise. It is therefore crucial to simulate and understand the processes associated with ice-shelf melt rates. In particular, the absence of tides representation in ocean models remains a caveat of numerous ocean hindcasts and climate projections. In the Amundsen Sea, tides are relatively weak and the melt-induced circulation is stronger than the tidal circulation. Using a regional 1/12° ocean model of the Amundsen Sea, we nonetheless find that tides can increase melt rates by up to 36% in some ice-shelf cavities. Among the processes that can possibly affect melt rates, the most important is an increased exchange at the ice/ocean interface resulting from the presence of strong tidal currents along the ice drafts. Approximately a third of this effect is compensated by a decrease in thermal forcing along the ice draft, which is related to an enhanced vertical mixing in the ocean interior in presence of tides. Parameterizing the effect of tides is an alternative to the representation of explicit tides in an ocean model, and has the advantage not to require any filtering of ocean model outputs. We therefore explore different ways to parameterize the effects of tides on ice shelf melt. First, we compare several methods to impose tidal velocities along the ice draft. We show that getting a realistic spatial distribution of tidal velocities in important, and can be deduced from the barotropic velocities of a tide model. Then, we explore several aspects of parameterized tidal mixing to reproduce the tide-induced decrease in thermal forcing along the ice drafts.

Oceans, Ice Shells, and Life on Europa

NASA Technical Reports Server (NTRS)

Schenk, Paul

2002-01-01

The four large satellites of Jupiter are famous for their planet-like diversity and complexity, but none more so than ice-covered Europa. Since the provocative Voyager images of Europa in 1979, evidence has been mounting that a vast liquid water ocean may lurk beneath the moon's icy surface. Europa has since been the target of increasing and sometimes reckless speculation regarding the possibility that giant squid and other creatures may be swimming its purported cold, dark ocean. No wonder Europa tops everyone's list for future exploration in the outer solar system (after the very first reconnaissance of Pluto and the Kuiper belt, of course). Europa may be the smallest of the Galilean moons (so-called because they were discovered by Galileo Galilei in the early 17th century) but more than makes up for its diminutive size with a crazed, alien landscape. The surface is covered with ridges hundreds of meters high, domes tens of kilometers across, and large areas of broken and disrupted crust called chaos. Some of the geologic features seen on Europa resemble ice rafts floating in polar seas here on Earth-reinforcing the idea that an ice shell is floating over an ocean on this Moon-size satellite. However, such features do not prove that an ocean exists or ever did. Warm ice is unusually soft and will flow under its own weight. If the ice shell is thick enough, the warm bottom of the shell will flow, as do terrestrial glaciers. This could produce all the observed surface features on Europa through a variety of processes, the most important of which is convection. (Convection is the vertical overturn of a layer due to heating or density differences-think of porridge or sauce boiling on the stove.) Rising blobs from the base of the crust would then create the oval domes dotting Europa's surface. The strongest evidence for a hidden ocean beneath Europa's surface comes from the Galileo spacecraft's onboard magnetometer, which detected fluctuations in Jupiter's magnetic

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

DOE Office of Scientific and Technical Information (OSTI.GOV)

Turner, Adrian Keith

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 pooledmore » 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).« less

Ice Floe Breaking in Contemporary Third Generation Operational Wave Models

NASA Astrophysics Data System (ADS)

Sévigny, C.; Baudry, J.; Gauthier, J. C.; Dumont, D.

2016-02-01

The dynamical zone observed at the edge of the consolidated ice area where are found the wave-fractured floes (i.e. marginal ice zone or MIZ) has become an important topic in ocean modeling. As both operational and climate ocean models now seek to reproduce the complex atmosphere-ice-ocean system with realistic coupling processes, many theoretical and numerical studies have focused on understanding and modeling this zone. Few attempts have been made to embed wave-ice interactions specific to the MIZ within a two-dimensional model, giving the possibility to calculate both the attenuation of surface waves by sea ice and the concomitant breaking of the sea ice-cover into smaller floes. One of the first challenges consists in improving the parameterization of wave-ice dynamics in contemporary third generation operational wave models. A simple waves-in-ice model (WIM) similar to the one proposed by Williams et al. (2013a,b) was implemented in WAVEWATCH III. This WIM considers ice floes as floating elastic plates and predicts the dimensionless attenuation coefficient by the use of a lookup-table-based, wave scattering scheme. As in Dumont et al. (2011), the different frequencies are treated individually and floe breaking occurs for a particular frequency when the expected wave amplitude exceeds the allowed strain amplitude, which considers ice floes properties and wavelength in ice field. The model is here further refined and tested in idealized two-dimensional cases, giving preliminary results of the performance and sensitivity of the parameterization to initial wave and ice conditions. The effects of the wave-ice coupling over the incident wave spectrum are analyzed as well as the resulting floe size distribution. The model gives prognostic values of the lateral extent of the marginal ice zone with maximum ice floe diameter that progressively increases with distance from the ice edge.

Free oscillations in a climate model with ice-sheet dynamics

NASA Technical Reports Server (NTRS)

Kallen, E.; Crafoord, C.; Ghil, M.

1979-01-01

A study of stable periodic solutions to a simple nonlinear model of the ocean-atmosphere-ice system is presented. The model has two dependent variables: ocean-atmosphere temperature and latitudinal extent of the ice cover. No explicit dependence on latitude is considered in the model. Hence all variables depend only on time and the model consists of a coupled set of nonlinear ordinary differential equations. The globally averaged ocean-atmosphere temperature in the model is governed by the radiation balance. The reflectivity to incoming solar radiation, i.e., the planetary albedo, includes separate contributions from sea ice and from continental ice sheets. The major physical mechanisms active in the model are (1) albedo-temperature feedback, (2) continental ice-sheet dynamics and (3) precipitation-rate variations. The model has three-equilibrium solutions, two of which are linearly unstable, while one is linearly stable. For some choices of parameters, the stability picture changes and sustained, finite-amplitude oscillations obtain around the previously stable equilibrium solution. The physical interpretation of these oscillations points to the possibility of internal mechanisms playing a role in glaciation cycles.

Energizing the Discussion of Ice-Ocean World Habitability

NASA Astrophysics Data System (ADS)

Schmidt, B.

2014-04-01

The outer solar system boasts a wide range of worlds with oceans - moons orbiting the gas giants as well as putative ocean worlds in the Kuiper Belt. These objects span sizes from a few hundred kilometers to larger than Mercury. How do we understand these bodies as a class as well as evaluate the habitability of individual environments? Recognizing that there is more to habitability than a set of ingredients, "Follow the Energy' has become an important mantra. Earth's biosphere is strongly coupled to its geologic activity that maintains a sort of stable chemical disequilibria that is employed by life. From this perspective, we can think of geologic activity as a planetary proxy for energy, setting up redox environments of which life can take advantage. With this as a backdrop, we will explore two of the most intriguing bodies: Europa and Enceladus. With an icy outer shell hiding a global ocean, Europa (r=1565 km) exists in a dynamic environment, where immense tides from Jupiter potentially power an active deeper interior. Intense irradiation and impacts bathe the top of the ice shell. These processes are sources of energy that could sustain a biosphere. In the past few decades the debate about habitability of Europa has been focused strongly on the thickness of its ice shell. However an arguably more critical question is: how does the ice shell really work? Galileo data indicated that Europa has undergone recent resurfacing, and implied that near-surface water was likely involved. Now the detection of potential water ice plumes, subduction-like behavior as well as shallow subsurface "lakes" within the past few years implies that rapid ice shell recycling could create a conveyor belt between the ice and ocean. Mediated by processes at the ice-ocean interface, exchange between Europa's surface and subsurface could allow ocean material to one day be detected or sampled by spacecraft. At least at this level, Europa passes the energy test. But the question remains: is

On the impact of ice-ocean interaction on Greenland glaciers versus calving speed.

NASA Astrophysics Data System (ADS)

Rignot, E. J.; Menemenlis, D.; Morlighem, M.; Wood, M.; Millan, R.; Mouginot, J.; An, L.

2016-12-01

Glacier retreat from frontal ablation is a delicate balance between subaqueous melt, calving processes and bed geometry. Here, we model subaqueous melt from a large number of Greenland tidewater glaciers using generalized 3D, high resolution simulations of ice melt from the MITgcm ocean model constrained by subglacial melt from RACMO2.3 and ISSM, ocean temperature from ECCO2-4km Arctic, and bed topography from OMG and MC for 1992-2015. The results are analyzed in combination with ice-front retreat and glacier speed from Landsat and imaging radar data since the 1990s. We find that subaqueous melt is 2-3 times greater in summer than in winter and doubled in magnitude since the 1990s because of enhanced ice sheet runoff and warmer ocean temperature. Glaciers that retreated rapidly are characterized by subaqueous melt rates comparable to their calving speed and favorable bed geometry. Glaciers dominated by calving processes are in contrast more resilient to thermal forcing from the ocean, especially in the presence of stabilizing geometry. The study highlights the fundamental importance of calving processes in controlling glacier retreat in Greenland.

Ice sheet-ocean interactions and sea level change

NASA Astrophysics Data System (ADS)

Heimbach, Patrick

2014-03-01

Mass loss from the Greenland and Antarctic ice sheets has increased rapidly since the mid-1990s. Their combined loss now accounts for about one-third of global sea level rise. In Greenland, a growing body of evidence points to the marine margins of these glaciers as the region from which this dynamic response originated. Similarly, ice streams in West Antarctica that feed vast floating ice shelves have exhibited large decadal changes. We review observational evidence and present physical mechanisms that might explain the observed changes, in particular in the context of ice sheet-ocean interactions. Processes involve cover 7 orders of magnitudes of scales, ranging from mm boundary-layer processes to basin-scale coupled atmosphere-ocean variability. We discuss observational needs to fill the gap in our mechanistic understanding.

Modeling tabular icebergs submerged in the ocean

NASA Astrophysics Data System (ADS)

Stern, A. A.; Adcroft, A.; Sergienko, O.; Marques, G.

2017-08-01

Large tabular icebergs calved from Antarctic ice shelves have long lifetimes (due to their large size), during which they drift across large distances, altering ambient ocean circulation, bottom-water formation, sea-ice formation, and biological primary productivity in the icebergs' vicinity. However, despite their importance, the current generation of ocean circulation models usually do not represent large tabular icebergs. In this study, we develop a novel framework to model large tabular icebergs submerged in the ocean. In this framework, tabular icebergs are represented by pressure-exerting Lagrangian elements that drift in the ocean. The elements are held together and interact with each other via bonds. A breaking of these bonds allows the model to emulate calving events (i.e., detachment of a tabular iceberg from an ice shelf) and tabular icebergs breaking up into smaller pieces. Idealized simulations of a calving tabular iceberg, its drift, and its breakup demonstrate capabilities of the developed framework.

The effect of changing wind forcing on Antarctic ice shelf melting in high-resolution, global sea ice-ocean simulations with the Accelerated Climate Model for Energy (ACME)

NASA Astrophysics Data System (ADS)

Asay-Davis, Xylar; Price, Stephen; Petersen, Mark; Wolfe, Jonathan

2017-04-01

The capability for simulating sub-ice shelf circulation and submarine melting and freezing has recently been added to the U.S. Department of Energy's Accelerated Climate Model for Energy (ACME). With this new capability, we use an eddy permitting ocean model to conduct two sets of simulations in the spirit of Spence et al. (GRL, 41, 2014), who demonstrate increased warm water upwelling along the Antarctic coast in response to poleward shifting and strengthening of Southern Ocean westerly winds. These characteristics, symptomatic of a positive Southern Annular Mode (SAM), are projected to continue into the 21st century under anthropogenic climate change (Fyfe et al., J. Clim., 20, 2007). In our first simulation, we force the climate model using the standard CORE interannual forcing dataset (Large and Yeager; Clim. Dyn., 33, 2009). In our second simulation, we force our climate model using an altered version of CORE interannual forcing, based on the latter half of the full time series, which we take as a proxy for a future climate state biased towards a positive SAM. We compare ocean model states and sub-ice shelf melt rates with observations, exploring sources of model biases as well as the effects of the two forcing scenarios.

Evolution of Summer Ocean Mixed Layer Heat Content and Ocean/Ice Fluxes in the Arctic Ocean During the Last Decade

NASA Astrophysics Data System (ADS)

Stanton, T. P.; Shaw, W. J.

2014-12-01

Since 2002, a series of 28 Autonomous Ocean Flux Buoys have been deployed in the Beaufort Sea and from the North Pole Environmental Observatory. These long-term ice-deployed instrument systems primarily measure vertical turbulent fluxes of heat, salt and momentum at a depth of 2 - 6 m below the ocean/ice interface, while concurrently measuring current profile every 2m down to approximately 40-50m depth, within the seasonal pycnocline. Additional sensors have been added to measure local ice melt rates acoustically, and finescale thermal structure from the eddy correlation flux sensor up into the ice to resolve summer near-surface heating. The AOFB buoys have typically been co-located with Ice Tethered Profilers, that measure the upper ocean T/S structure and ice mass balance instruments. Comparisons of near-surface heat fluxes, heat content and vertical structure over the last decade will be made for buoys in the Beaufort Sea and Transpolar Drift between the North Pole and Spitzbergen. The effects of enhanced basal melting from ice/albedo feedbacks can be clearly seen in the low ice concentration summer conditions found more recently in the Beaufort Sea, while there are less pronounced effects of enhanced summer surface heating in the higher ice concentrations still found in the transpolar drift.

Modelling wave-induced sea ice break-up in the marginal ice zone

NASA Astrophysics Data System (ADS)

Montiel, F.; Squire, V. A.

2017-10-01

A model of ice floe break-up under ocean wave forcing in the marginal ice zone (MIZ) is proposed to investigate how floe size distribution (FSD) evolves under repeated wave break-up events. A three-dimensional linear model of ocean wave scattering by a finite array of compliant circular ice floes is coupled to a flexural failure model, which breaks a floe into two floes provided the two-dimensional stress field satisfies a break-up criterion. A closed-feedback loop algorithm is devised, which (i) solves the wave-scattering problem for a given FSD under time-harmonic plane wave forcing, (ii) computes the stress field in all the floes, (iii) fractures the floes satisfying the break-up criterion, and (iv) generates an updated FSD, initializing the geometry for the next iteration of the loop. The FSD after 50 break-up events is unimodal and near normal, or bimodal, suggesting waves alone do not govern the power law observed in some field studies. Multiple scattering is found to enhance break-up for long waves and thin ice, but to reduce break-up for short waves and thick ice. A break-up front marches forward in the latter regime, as wave-induced fracture weakens the ice cover, allowing waves to travel deeper into the MIZ.

The implementation of sea ice model on a regional high-resolution scale

NASA Astrophysics Data System (ADS)

Prasad, Siva; Zakharov, Igor; Bobby, Pradeep; McGuire, Peter

2015-09-01

The availability of high-resolution atmospheric/ocean forecast models, satellite data and access to high-performance computing clusters have provided capability to build high-resolution models for regional ice condition simulation. The paper describes the implementation of the Los Alamos sea ice model (CICE) on a regional scale at high resolution. The advantage of the model is its ability to include oceanographic parameters (e.g., currents) to provide accurate results. The sea ice simulation was performed over Baffin Bay and the Labrador Sea to retrieve important parameters such as ice concentration, thickness, ridging, and drift. Two different forcing models, one with low resolution and another with a high resolution, were used for the estimation of sensitivity of model results. Sea ice behavior over 7 years was simulated to analyze ice formation, melting, and conditions in the region. Validation was based on comparing model results with remote sensing data. The simulated ice concentration correlated well with Advanced Microwave Scanning Radiometer for EOS (AMSR-E) and Ocean and Sea Ice Satellite Application Facility (OSI-SAF) data. Visual comparison of ice thickness trends estimated from the Soil Moisture and Ocean Salinity satellite (SMOS) agreed with the simulation for year 2010-2011.

Modeling Wave-Ice Interactions in the Marginal Ice Zone

NASA Astrophysics Data System (ADS)

Orzech, Mark; Shi, Fengyan; Bateman, Sam; Veeramony, Jay; Calantoni, Joe

2015-04-01

., F. Shi, and J.T. Kirby (2012). Shock-capturing non-hydrostatic model for fully dispersive surface wave processes. Ocean Modelling 43-44, 22-35. >Wadhams P., V. Squire, J.A. Ewing, and R.W. Pascal (1986). The effect of the marginal ice zone on the directional wave spectrum of the ocean. J. Phys. Oceanog., 16(2), 358-376.

Testing a common ice-ocean parameterization with laboratory experiments

NASA Astrophysics Data System (ADS)

McConnochie, C. D.; Kerr, R. C.

2017-07-01

Numerical models of ice-ocean interactions typically rely upon a parameterization for the transport of heat and salt to the ice face that has not been satisfactorily validated by observational or experimental data. We compare laboratory experiments of ice-saltwater interactions to a common numerical parameterization and find a significant disagreement in the dependence of the melt rate on the fluid velocity. We suggest a resolution to this disagreement based on a theoretical analysis of the boundary layer next to a vertical heated plate, which results in a threshold fluid velocity of approximately 4 cm/s at driving temperatures between 0.5 and 4°C, above which the form of the parameterization should be valid.

Monitoring Subsurface Ice-Ocean Processes Using Underwater Acoustics in the Ross Sea

NASA Astrophysics Data System (ADS)

Haxel, J. H.; Dziak, R. P.; Matsumoto, H.; Lee, W. S.; Yun, S.

2016-12-01

The Ross Sea is a dynamic area of ice-ocean interaction, where a large component of the Southern Ocean's sea ice formation occurs within regional polynyas in addition to the destructive processes happening at the seaward boundary of the Ross Ice Shelf. Recent studies show the sea-ice season has been lengthening and the sea ice extent has been growing with more persistent and larger regional polynyas. These trends have important implications for the Ross Sea ecosystem with polynyas supporting high rates of primary productivity in the area. Monitoring trends in sea ice and ice shelf dynamics in the Southern Ocean has relied heavily on satellite imagery and remote sensing methods despite a significant portion of these physical processes occurring beneath the ocean surface. In January 2014, an ocean bottom hydrophone (OBH) was moored on the seafloor in the polynya area of Terra Nova Bay in the northwest region of the Ross Sea, north of the Drygalski Ice Tongue. The OBH recorded a year long record of the underwater low frequency acoustic spectrum up to 500 Hz from January 29 until it was recovered the following December 17, 2014. The acoustic records reveal a complex annual history of ice generated signals with over 50,000 detected events. These ice generated events related to collisions and cracking provide important insight for the timing and intensity of the ice-ocean dynamics happening below the sea surface as the polynya grows and expands and the nearby Drygalski ice tongue flows into Terra Nova Bay. Additionally, high concentrations of baleen whale vocalizations in frequencies ranging from 200-400 Hz from September - December suggest a strong seasonal presence of whales in this ecologically important polynya region.

The Role of Late Summer Melt Pond Water Layers in the Ocean Mixed Layer on Enhancing Ice/Ocean Albedo Feedbacks in the Arctic

NASA Astrophysics Data System (ADS)

Stanton, T. P.; Shaw, W. J.

2016-02-01

Drainage of surface melt pond water into the top of the ocean mixed layer is seen widely in the Arctic ice pack in later summer (for example Gallaher et al 2015). Under calm conditions, this fresh water forms a thin, stratified layer immediately below the ice which is dynamically decoupled from the thicker, underlying seasonal mixed layer by the density difference between the two layers. The ephemeral surface layer is significantly warmer than the underlying ocean water owing to the higher freezing temperature of the fresh melt water. How the presence of this warm ephemeral layer enhances basal melt rate and speeds the destruction of the floes is investigated. High resolution timeseries measurements of T/S profiles in the 2m of the ocean immediately below the ice, and eddy-correlation fluxes of heat, salt and momentum 2.5m below the ice were made from an Autonomous Ocean Flux Buoy over a 2 month interval in later summer of 2015 as a component of the ONR Marginal Ice Zone project. The stratification and turbulent forcing observations are used with a 1 D turbulence closure model to understand how momentum and incoming radiative energy are stored and redistributed within the ephemeral layer. Under low wind forcing conditions both turbulent mixing energy and the water with high departure from freezing are trapped in the ephemeral layer by the strong density gradient at the base of the layer, resulting in rapid basal melting. This case is contrasted with model runs where the ephemeral layer heat is allowed to mix across the seasonal mixed layer, which results in slower basal melt rates. Consequently, the salinity-trapped warm ephemeral layer results in the formation of more open water earlier in the summer season, in turn resulting in increased cumulative heating of the ocean mixed layer, enhancing ice/ocean albedo feedbacks.

Modelling sea ice formation in the Terra Nova Bay polynya

NASA Astrophysics Data System (ADS)

Sansiviero, M.; Morales Maqueda, M. Á.; Fusco, G.; Aulicino, G.; Flocco, D.; Budillon, G.

2017-02-01

Antarctic sea ice is constantly exported from the shore by strong near surface winds that open leads and large polynyas in the pack ice. The latter, known as wind-driven polynyas, are responsible for significant water mass modification due to the high salt flux into the ocean associated with enhanced ice growth. In this article, we focus on the wind-driven Terra Nova Bay (TNB) polynya, in the western Ross Sea. Brine rejected during sea ice formation processes that occur in the TNB polynya densifies the water column leading to the formation of the most characteristic water mass of the Ross Sea, the High Salinity Shelf Water (HSSW). This water mass, in turn, takes part in the formation of Antarctic Bottom Water (AABW), the densest water mass of the world ocean, which plays a major role in the global meridional overturning circulation, thus affecting the global climate system. A simple coupled sea ice-ocean model has been developed to simulate the seasonal cycle of sea ice formation and export within a polynya. The sea ice model accounts for both thermal and mechanical ice processes. The oceanic circulation is described by a one-and-a-half layer, reduced gravity model. The domain resolution is 1 km × 1 km, which is sufficient to represent the salient features of the coastline geometry, notably the Drygalski Ice Tongue. The model is forced by a combination of Era Interim reanalysis and in-situ data from automatic weather stations, and also by a climatological oceanic dataset developed from in situ hydrographic observations. The sensitivity of the polynya to the atmospheric forcing is well reproduced by the model when atmospheric in situ measurements are combined with reanalysis data. Merging the two datasets allows us to capture in detail the strength and the spatial distribution of the katabatic winds that often drive the opening of the polynya. The model resolves fairly accurately the sea ice drift and sea ice production rates in the TNB polynya, leading to

Interannual variability of monthly Southern Ocean sea ice distributions

NASA Technical Reports Server (NTRS)

Parkinson, Claire L.

1992-01-01

The interannual variability of the Southern-Ocean sea-ice distributions was mapped and analyzed using data from Nimbus-5 ESMR and Nimbus-7 SMMR, collected from 1973 to 1987. The set of 12 monthly maps obtained reveals many details on spatial variability that are unobtainable from time series of ice extents. These maps can be used as baseline maps for comparisons against future Southern Ocean sea ice distributions. The maps are supplemented by more detailed maps of the frequency of ice coverage, presented in this paper for one month within each of the four seasons, and by the breakdown of these results to the periods covered individually by each of the two passive-microwave imagers.

Dynamics of the last glacial maximum Antarctic ice-sheet and its response to ocean forcing

PubMed Central

Golledge, Nicholas R.; Fogwill, Christopher J.; Mackintosh, Andrew N.; Buckley, Kevin M.

2012-01-01

Retreat of the Last Glacial Maximum (LGM) Antarctic ice sheet is thought to have been initiated by changes in ocean heat and eustatic sea level propagated from the Northern Hemisphere (NH) as northern ice sheets melted under rising atmospheric temperatures. The extent to which spatial variability in ice dynamics may have modulated the resultant pattern and timing of decay of the Antarctic ice sheet has so far received little attention, however, despite the growing recognition that dynamic effects account for a sizeable proportion of mass-balance changes observed in modern ice sheets. Here we use a 5-km resolution whole-continent numerical ice-sheet model to assess whether differences in the mechanisms governing ice sheet flow could account for discrepancies between geochronological studies in different parts of the continent. We first simulate the geometry and flow characteristics of an equilibrium LGM ice sheet, using pan-Antarctic terrestrial and marine geological data for constraint, then perturb the system with sea level and ocean heat flux increases to investigate ice-sheet vulnerability. Our results identify that fast-flowing glaciers in the eastern Weddell Sea, the Amundsen Sea, central Ross Sea, and in the Amery Trough respond most rapidly to ocean forcings, in agreement with empirical data. Most significantly, we find that although ocean warming and sea-level rise bring about mainly localized glacier acceleration, concomitant drawdown of ice from neighboring areas leads to widespread thinning of entire glacier catchments—a discovery that has important ramifications for the dynamic changes presently being observed in modern ice sheets. PMID:22988078

Dynamics of the last glacial maximum Antarctic ice-sheet and its response to ocean forcing.

PubMed

Golledge, Nicholas R; Fogwill, Christopher J; Mackintosh, Andrew N; Buckley, Kevin M

2012-10-02

Retreat of the Last Glacial Maximum (LGM) Antarctic ice sheet is thought to have been initiated by changes in ocean heat and eustatic sea level propagated from the Northern Hemisphere (NH) as northern ice sheets melted under rising atmospheric temperatures. The extent to which spatial variability in ice dynamics may have modulated the resultant pattern and timing of decay of the Antarctic ice sheet has so far received little attention, however, despite the growing recognition that dynamic effects account for a sizeable proportion of mass-balance changes observed in modern ice sheets. Here we use a 5-km resolution whole-continent numerical ice-sheet model to assess whether differences in the mechanisms governing ice sheet flow could account for discrepancies between geochronological studies in different parts of the continent. We first simulate the geometry and flow characteristics of an equilibrium LGM ice sheet, using pan-Antarctic terrestrial and marine geological data for constraint, then perturb the system with sea level and ocean heat flux increases to investigate ice-sheet vulnerability. Our results identify that fast-flowing glaciers in the eastern Weddell Sea, the Amundsen Sea, central Ross Sea, and in the Amery Trough respond most rapidly to ocean forcings, in agreement with empirical data. Most significantly, we find that although ocean warming and sea-level rise bring about mainly localized glacier acceleration, concomitant drawdown of ice from neighboring areas leads to widespread thinning of entire glacier catchments-a discovery that has important ramifications for the dynamic changes presently being observed in modern ice sheets.

Quaternary Sea-ice history in the Arctic Ocean based on a new Ostracode sea-ice proxy

USGS Publications Warehouse

Cronin, T. M.; Gemery, L.; Briggs, W.M.; Jakobsson, M.; Polyak, L.; Brouwers, E.M.

2010-01-01

Paleo-sea-ice history in the Arctic Ocean was reconstructed using the sea-ice dwelling ostracode Acetabulastoma arcticum from late Quaternary sediments from the Mendeleyev, Lomonosov, and Gakkel Ridges, the Morris Jesup Rise and the Yermak Plateau. Results suggest intermittently high levels of perennial sea ice in the central Arctic Ocean during Marine Isotope Stage (MIS) 3 (25-45 ka), minimal sea ice during the last deglacial (16-11 ka) and early Holocene thermal maximum (11-5 ka) and increasing sea ice during the mid-to-late Holocene (5-0 ka). Sediment core records from the Iceland and Rockall Plateaus show that perennial sea ice existed in these regions only during glacial intervals MIS 2, 4, and 6. These results show that sea ice exhibits complex temporal and spatial variability during different climatic regimes and that the development of modern perennial sea ice may be a relatively recent phenomenon. ?? 2010.

Current State and Recent Changes in the Arctic Ocean from the HYCOM-NCODA Global Ocean and Sea Ice Prediction System

NASA Astrophysics Data System (ADS)

Dukhovskoy, D. S.; Chassignet, E. P.; Hogan, P. J.; Metzger, E. J.; Posey, P.; Smedstad, O. M.; Stefanova, L. B.; Wallcraft, A. J.

2016-12-01

The great potential of numerical models to provide a high-resolution continuous picture of the environmental characteristics of the Arctic system is related to the problem of reliability and accuracy of the simulations. Recent Arctic Ocean model intercomparison projects have identified substantial disagreements in water mass distribution and circulation among the models over the last two decades. In situ and satellite observations cannot yield enough continuous in time and space information to interpret the observed changes in the Arctic system. Observations combined with Arctic Ocean models via data assimilation provide perhaps the most complete knowledge about the state of the Arctic system. We use outputs from the US Navy Global Ocean Forecast System (20-year reanalysis + analysis) to investigate several hypotheses that have been put forward regarding the current state and recent changes in the Arctic Ocean. The system is based on the 0.08-degree HYbrid Coordinate Ocean Model (HYCOM) and can be run with two-way coupling to the Los Alamos Community Ice CodE (CICE) or with an energy-loan ice model. Observations are assimilated by the Navy Coupled Ocean Data Assimilation (NCODA) algorithm. HYCOM temperature and salinity fields are shown to be in good agreement with observational data in the Arctic and North Atlantic. The model reproduces changes in the freshwater budget in the Arctic as reported in other studies. The modeled freshwater fluxes between the Arctic Ocean and the North Atlantic are analyzed to document and discuss the interaction between the two regions over the last two decades.

Reassessment of ice-age cooling of the tropical ocean and atmosphere

USGS Publications Warehouse

Hostetler, S.W.; Mix, A.C.

1999-01-01

The CLIMAP project's reconstruction of past sea surface temperature inferred limited ice-age cooling in the tropical oceans. This conclusion has been controversial, however, because of the greater cooling indicated by other terrestrial and ocean proxy data. A new faunal sea surface temperature reconstruction, calibrated using the variation of foraminiferal species through time, better represents ice-age faunal assemblages and so reveals greater cooling than CLIMAP in the equatorial current systems of the eastern Pacific and tropical Atlantic oceans. Here we explore the climatic implications of this revised sea surface temperature field for the Last Glacial Maximum using an atmospheric general circulation model. Relative to model results obtained using CLIMAP sea surface temperatures, the cooler equatorial oceans modify seasonal air temperatures by 1-2??C or more across parts of South America, Africa and southeast Asia and cause attendant changes in regional moisture patterns. In our simulation of the Last Glacial Maximum, the Amazon lowlands, for example, are cooler and drier, whereas the Andean highlands are cooler and wetter than the control simulation. Our results may help to resolve some of the apparent disagreements between oceanic and continental proxy climate data. Moreover, they suggest a wind-related mechanism for enhancing the export of water vapour from the Atlantic to the Indo-Pacific oceans, which may link variations in deep-water production and high-latitude climate changes to equatorial sea surface temperatures.

The force balance of sea ice in a numerical model of the Arctic Ocean

NASA Astrophysics Data System (ADS)

Steele, Michael; Zhang, Jinlun; Rothrock, Drew; Stern, Harry

1997-09-01

The balance of forces in the sea ice model of Hibler [1979] is examined. The model predicts that internal stress gradients are an important force in much of the Arctic Ocean except in summer, when they are significant only off the northern coasts of Greenland and the Canadian Archipelago. A partition of the internal stress gradient between the pressure gradient and the viscous terms reveals that both are significant, although they operate on very different timescales. The acceleration term is generally negligible, while the sum of Coriolis plus sea surface tilt is small. Thus the seasonal average force balance in fall, winter, and spring is mostly between three terms of roughly equal magnitudes: air drag, water drag, and internal stress gradients. This is also true for the monthly average force balance. However, we find that there is a transition around the weekly timescale and that on a daily basis the force balance at a particular location and time is often between only two terms: either between air drag and water drag or between air drag and internal stress gradients. The model is in agreement with the observations of Thorndike and Colony [1982] in that the correlation between geostrophic wind forcing and the model's ice velocity field is high. This result is discussed in the context of the force balance; we show that the presence of significant internal stress gradients does not preclude high wind-ice correlation. A breakdown of the internal stress gradient into component parts reveals that the shear viscous force is far from negligible, which casts strong doubt on the theoretical validity of the cavitating fluid approximation (in which this component is neglected). Finally, the role of ice pressure is examined by varying the parameter P*. We find a strong sensitivity in terms of the force balance, as well as ice thickness and velocity.

Towards an Ice-Free Arctic Ocean in Summertime

NASA Astrophysics Data System (ADS)

Gascard, Jean Claude

2014-05-01

Dividing the Arctic Ocean in two parts, the so-called Atlantic versus the Pacific sector, two distinct modes of variability appear for characterizing the Arctic sea-ice extent from 70°N up to 80°N in both sectors. The Atlantic sector seasonal sea-ice extent is characterized by a longer time scale than the Pacific sector with a break up melting season starting in May and reaching a peak in June-July, one month earlier than the Pacific sector of the Arctic Ocean revealing a faster time evolution and a larger spatial amplitude than the Atlantic sector. During recent years like 2007, sea-ice extent with sea-ice concentration above 15% retreated from 4 millions km2 to about 1 million km2 in the Arctic Pacific sector between 70° and 80°N except for 2012 when most of sea-ice melted away in this region. That explained most of the differences between the two extreme years 2007 and 2012. In the Atlantic sector, Arctic sea-ice retreated from 2 millions km2 to nearly 0 during recent years including 2007 and 2012. The Atlantic inflow North of Svalbard and Franz Josef Land is more likely responsible for a northward retreat of the ice edge in that region. The important factor is not only that the Arctic summer sea-ice minimum extent decreased by 3 or 4 millions km2 over the past 10 years but also that the melting period was steadily increasing by one to two days every year during that period. An important factor concerns the strength of the freezing that can be quantified in terms of Freezing Degree Days FDD accumulated during the winter-spring season and the strength of the melting (MDD) that can be accumulated during the summer season. FDD and MDD have been calculated for the past 30 years all over the Arctic Ocean using ERA Interim Reanalysis surface temperature at 2m height in the atmosphere. It is clear that FDD decreased significantly by more than 2000 FDD between 1980 and 2012 which is equivalent to the sensible heat flux corresponding to more than a meter of sea-ice

Decadal predictions of Southern Ocean sea ice : testing different initialization methods with an Earth-system Model of Intermediate Complexity

NASA Astrophysics Data System (ADS)

Zunz, Violette; Goosse, Hugues; Dubinkina, Svetlana

2013-04-01

The sea ice extent in the Southern Ocean has increased since 1979 but the causes of this expansion have not been firmly identified. In particular, the contribution of internal variability and external forcing to this positive trend has not been fully established. In this region, the lack of observations and the overestimation of internal variability of the sea ice by contemporary General Circulation Models (GCMs) make it difficult to understand the behaviour of the sea ice. Nevertheless, if its evolution is governed by the internal variability of the system and if this internal variability is in some way predictable, a suitable initialization method should lead to simulations results that better fit the reality. Current GCMs decadal predictions are generally initialized through a nudging towards some observed fields. This relatively simple method does not seem to be appropriated to the initialization of sea ice in the Southern Ocean. The present study aims at identifying an initialization method that could improve the quality of the predictions of Southern Ocean sea ice at decadal timescales. We use LOVECLIM, an Earth-system Model of Intermediate Complexity that allows us to perform, within a reasonable computational time, the large amount of simulations required to test systematically different initialization procedures. These involve three data assimilation methods: a nudging, a particle filter and an efficient particle filter. In a first step, simulations are performed in an idealized framework, i.e. data from a reference simulation of LOVECLIM are used instead of observations, herein after called pseudo-observations. In this configuration, the internal variability of the model obviously agrees with the one of the pseudo-observations. This allows us to get rid of the issues related to the overestimation of the internal variability by models compared to the observed one. This way, we can work out a suitable methodology to assess the efficiency of the

Eddy-resolving simulations of the Fimbul Ice Shelf cavity circulation: Basal melting and exchange with open ocean

NASA Astrophysics Data System (ADS)

Hattermann, T.; Smedsrud, L. H.; Nøst, O. A.; Lilly, J. M.; Galton-Fenzi, B. K.

2014-10-01

Melting at the base of floating ice shelves is a dominant term in the overall Antarctic mass budget. This study applies a high-resolution regional ice shelf/ocean model, constrained by observations, to (i) quantify present basal mass loss at the Fimbul Ice Shelf (FIS); and (ii) investigate the oceanic mechanisms that govern the heat supply to ice shelves in the Eastern Weddell Sea. The simulations confirm the low melt rates suggested by observations and show that melting is primarily determined by the depth of the coastal thermocline, regulating deep ocean heat fluxes towards the ice. Furthermore, the uneven distribution of ice shelf area at different depths modulates the melting response to oceanic forcing, causing the existence of two distinct states of melting at the FIS. In the simulated present-day state, only small amounts of Modified Warm Deep Water enter the continental shelf, and ocean temperatures beneath the ice are close to the surface freezing point. The basal mass loss in this so-called state of "shallow melting" is mainly controlled by the seasonal inflow of solar-heated surface water affecting large areas of shallow ice in the upper part of the cavity. This is in contrast to a state of "deep melting", in which the thermocline rises above the shelf break depth, establishing a continuous inflow of Warm Deep Water towards the deep ice. The transition between the two states is found to be determined by a complex response of the Antarctic Slope Front overturning circulation to varying climate forcings. A proper representation of these frontal dynamics in climate models will therefore be crucial when assessing the evolution of ice shelf basal melting along this sector of Antarctica.

Ocean wave generation by collapsing ice shelves

NASA Astrophysics Data System (ADS)

Macayeal, D. R.; Bassis, J. N.; Okal, E. A.; Aster, R. C.; Cathles, L. M.

2008-12-01

The 28-29 February, 2008, break-up of the Wilkins Ice Shelf, Antarctica, exemplifies the now-familiar, yet largely unexplained pattern of explosive ice-shelf break-up. While environmental warming is a likely ultimate cause of explosive break-up, several key aspects of their short-term behavior need to be explained: (1) The abrupt, near-simultaneous onset of iceberg calving across long spans of the ice front margin; (2) High outward drift velocity (about 0.3 m/s) of a leading phalanx of tabular icebergs that originate from the seaward edge of the intact ice shelf prior to break-up; (3) Rapid coverage of the ocean surface in the wake of this leading phalanx by small, capsized and dismembered tabular icebergs; (4) Extremely large gravitational potential energy release rates, e.g., up to 3 × 1010 W; (5) Lack of proximal iceberg-calving triggers that control the timing of break-up onset and that maintain the high break-up calving rates through to the conclusion of the event. Motivated by seismic records obtained from icebergs and the Ross Ice Shelf that show hundreds of micro- tsunamis emanating from near the ice shelf front, we re-examine the basic dynamic features of ice- shelf/ocean-wave interaction and, in particular, examine the possibility that collapsing ice shelves themselves are a source of waves that stimulate the disintegration process. We propose that ice-shelf generated surface-gravity waves associated with initial calving at an arbitrary seed location produce stress perturbations capable of triggering the onset of calving on the entire ice front. Waves generated by parting detachment rifts, iceberg capsize and break-up act next to stimulate an inverted submarine landslide (ice- slide) process, where gravitational potential energy released by upward movement of buoyant ice is radiated as surface gravity waves in the wake of the advancing phalanx of tabular icebergs. We conclude by describing how field research and remote sensing can be used to test the

Upper Ocean Evolution Across the Beaufort Sea Marginal Ice Zone

NASA Astrophysics Data System (ADS)

Lee, C.; Rainville, L.; Gobat, J. I.; Perry, M. J.; Freitag, L. E.; Webster, S.

2016-12-01

The observed reduction of Arctic summertime sea ice extent and expansion of the marginal ice zone (MIZ) have profound impacts on the balance of processes controlling sea ice evolution, including the introduction of several positive feedback mechanisms that may act to accelerate melting. Examples of such feedbacks include increased upper ocean warming though absorption of solar radiation, elevated internal wave energy and mixing that may entrain heat stored in subsurface watermasses (e.g., the relatively warm Pacific Summer and Atlantic waters), and elevated surface wave energy that acts to deform and fracture sea ice. Spatial and temporal variability in ice properties and open water fraction impact these processes. To investigate how upper ocean structure varies with changing ice cover, how the balance of processes shift as a function of ice fraction and distance from open water, and how these processes impact sea ice evolution, a network of autonomous platforms sampled the atmosphere-ice-ocean system in the Beaufort, beginning in spring, well before the start of melt, and ending with the autumn freeze-up. Four long-endurance autonomous Seagliders occupied sections that extended from open water, through the marginal ice zone, deep into the pack during summer 2014 in the Beaufort Sea. Gliders penetrated up to 200 km into the ice pack, under complete ice cover for up to 10 consecutive days. Sections reveal strong fronts where cold, ice-covered waters meet waters that have been exposed to solar warming, and O(10 km) scale eddies near the ice edge. In the pack, Pacific Summer Water and a deep chlorophyll maximum form distinct layers at roughly 60 m and 80 m, respectively, which become increasingly diffuse late in the season as they progress through the MIZ and into open water. Stratification just above the Pacific Summer Water rapidly weakens near the ice edge and temperature variance increases, likely due to mixing or energetic vertical exchange associated with strong

A conceptual model of oceanic heat transport in the Snowball Earth scenario

NASA Astrophysics Data System (ADS)

Comeau, Darin; Kurtze, Douglas A.; Restrepo, Juan M.

2016-12-01

Geologic evidence suggests that the Earth may have been completely covered in ice in the distant past, a state known as Snowball Earth. This is still the subject of controversy, and has been the focus of modeling work from low-dimensional models up to state-of-the-art general circulation models. In our present global climate, the ocean plays a large role in redistributing heat from the equatorial regions to high latitudes, and as an important part of the global heat budget, its role in the initiation a Snowball Earth, and the subsequent climate, is of great interest. To better understand the role of oceanic heat transport in the initiation of Snowball Earth, and the resulting global ice covered climate state, the goal of this inquiry is twofold: we wish to propose the least complex model that can capture the Snowball Earth scenario as well as the present-day climate with partial ice cover, and we want to determine the relative importance of oceanic heat transport. To do this, we develop a simple model, incorporating thermohaline dynamics from traditional box ocean models, a radiative balance from energy balance models, and the more contemporary "sea glacier" model to account for viscous flow effects of extremely thick sea ice. The resulting model, consisting of dynamic ocean and ice components, is able to reproduce both Snowball Earth and present-day conditions through reasonable changes in forcing parameters. We find that including or neglecting oceanic heat transport may lead to vastly different global climate states, and also that the parameterization of under-ice heat transfer in the ice-ocean coupling plays a key role in the resulting global climate state, demonstrating the regulatory effect of dynamic ocean heat transport.

Evidence for an ice shelf covering the central Arctic Ocean during the penultimate glaciation

USGS Publications Warehouse

Jakobsson, Martin; Nilsson, Johan; Anderson, Leif G.; Backman, Jan; Bjork, Goran; Cronin, Thomas M.; Kirchner, Nina; Koshurnikov, Andrey; Mayer, Larry; Noormets, Riko; O'Regan, Matthew; Stranne, Christian; Ananiev, Roman; Macho, Natalia Barrientos; Cherniykh, Dennis; Coxall, Helen; Eriksson, Bjorn; Floden, Tom; Gemery, Laura; Gustafsson, Orjan; Jerram, Kevin; Johansson, Carina; Khortov, Alexey; Mohammad, Rezwan; Semiletov, Igor

2016-01-01

The hypothesis of a km-thick ice shelf covering the entire Arctic Ocean during peak glacial conditions was proposed nearly half a century ago. Floating ice shelves preserve few direct traces after their disappearance, making reconstructions difficult. Seafloor imprints of ice shelves should, however, exist where ice grounded along their flow paths. Here we present new evidence of ice-shelf groundings on bathymetric highs in the central Arctic Ocean, resurrecting the concept of an ice shelf extending over the entire central Arctic Ocean during at least one previous ice age. New and previously mapped glacial landforms together reveal flow of a spatially coherent, in some regions >1-km thick, central Arctic Ocean ice shelf dated to marine isotope stage 6 (~140 ka). Bathymetric highs were likely critical in the ice-shelf development by acting as pinning points where stabilizing ice rises formed, thereby providing sufficient back stress to allow ice shelf thickening.

Evidence for an ice shelf covering the central Arctic Ocean during the penultimate glaciation

PubMed Central

Jakobsson, Martin; Nilsson, Johan; Anderson, Leif; Backman, Jan; Björk, Göran; Cronin, Thomas M.; Kirchner, Nina; Koshurnikov, Andrey; Mayer, Larry; Noormets, Riko; O'Regan, Matthew; Stranne, Christian; Ananiev, Roman; Barrientos Macho, Natalia; Cherniykh, Denis; Coxall, Helen; Eriksson, Björn; Flodén, Tom; Gemery, Laura; Gustafsson, Örjan; Jerram, Kevin; Johansson, Carina; Khortov, Alexey; Mohammad, Rezwan; Semiletov, Igor

2016-01-01

The hypothesis of a km-thick ice shelf covering the entire Arctic Ocean during peak glacial conditions was proposed nearly half a century ago. Floating ice shelves preserve few direct traces after their disappearance, making reconstructions difficult. Seafloor imprints of ice shelves should, however, exist where ice grounded along their flow paths. Here we present new evidence of ice-shelf groundings on bathymetric highs in the central Arctic Ocean, resurrecting the concept of an ice shelf extending over the entire central Arctic Ocean during at least one previous ice age. New and previously mapped glacial landforms together reveal flow of a spatially coherent, in some regions >1-km thick, central Arctic Ocean ice shelf dated to marine isotope stage 6 (∼140 ka). Bathymetric highs were likely critical in the ice-shelf development by acting as pinning points where stabilizing ice rises formed, thereby providing sufficient back stress to allow ice shelf thickening. PMID:26778247

Sea-ice transport driving Southern Ocean salinity and its recent trends.

PubMed

Haumann, F Alexander; Gruber, Nicolas; Münnich, Matthias; Frenger, Ivy; Kern, Stefan

2016-09-01

Recent salinity changes in the Southern Ocean are among the most prominent signals of climate change in the global ocean, yet their underlying causes have not been firmly established. Here we propose that trends in the northward transport of Antarctic sea ice are a major contributor to these changes. Using satellite observations supplemented by sea-ice reconstructions, we estimate that wind-driven northward freshwater transport by sea ice increased by 20 ± 10 per cent between 1982 and 2008. The strongest and most robust increase occurred in the Pacific sector, coinciding with the largest observed salinity changes. We estimate that the additional freshwater for the entire northern sea-ice edge entails a freshening rate of -0.02 ± 0.01 grams per kilogram per decade in the surface and intermediate waters of the open ocean, similar to the observed freshening. The enhanced rejection of salt near the coast of Antarctica associated with stronger sea-ice export counteracts the freshening of both continental shelf and newly formed bottom waters due to increases in glacial meltwater. Although the data sources underlying our results have substantial uncertainties, regional analyses and independent data from an atmospheric reanalysis support our conclusions. Our finding that northward sea-ice freshwater transport is also a key determinant of the mean salinity distribution in the Southern Ocean further underpins the importance of the sea-ice-induced freshwater flux. Through its influence on the density structure of the ocean, this process has critical consequences for the global climate by affecting the exchange of heat, carbon and nutrients between the deep ocean and surface waters.

Modelling wave-induced sea ice break-up in the marginal ice zone

PubMed Central

Squire, V. A.

2017-01-01

A model of ice floe break-up under ocean wave forcing in the marginal ice zone (MIZ) is proposed to investigate how floe size distribution (FSD) evolves under repeated wave break-up events. A three-dimensional linear model of ocean wave scattering by a finite array of compliant circular ice floes is coupled to a flexural failure model, which breaks a floe into two floes provided the two-dimensional stress field satisfies a break-up criterion. A closed-feedback loop algorithm is devised, which (i) solves the wave-scattering problem for a given FSD under time-harmonic plane wave forcing, (ii) computes the stress field in all the floes, (iii) fractures the floes satisfying the break-up criterion, and (iv) generates an updated FSD, initializing the geometry for the next iteration of the loop. The FSD after 50 break-up events is unimodal and near normal, or bimodal, suggesting waves alone do not govern the power law observed in some field studies. Multiple scattering is found to enhance break-up for long waves and thin ice, but to reduce break-up for short waves and thick ice. A break-up front marches forward in the latter regime, as wave-induced fracture weakens the ice cover, allowing waves to travel deeper into the MIZ. PMID:29118659

Modelling wave-induced sea ice break-up in the marginal ice zone.

PubMed

Montiel, F; Squire, V A

2017-10-01

A model of ice floe break-up under ocean wave forcing in the marginal ice zone (MIZ) is proposed to investigate how floe size distribution (FSD) evolves under repeated wave break-up events. A three-dimensional linear model of ocean wave scattering by a finite array of compliant circular ice floes is coupled to a flexural failure model, which breaks a floe into two floes provided the two-dimensional stress field satisfies a break-up criterion. A closed-feedback loop algorithm is devised, which (i) solves the wave-scattering problem for a given FSD under time-harmonic plane wave forcing, (ii) computes the stress field in all the floes, (iii) fractures the floes satisfying the break-up criterion, and (iv) generates an updated FSD, initializing the geometry for the next iteration of the loop. The FSD after 50 break-up events is unimodal and near normal, or bimodal, suggesting waves alone do not govern the power law observed in some field studies. Multiple scattering is found to enhance break-up for long waves and thin ice, but to reduce break-up for short waves and thick ice. A break-up front marches forward in the latter regime, as wave-induced fracture weakens the ice cover, allowing waves to travel deeper into the MIZ.

Sensitivity study of a dynamic thermodynamic sea ice model

NASA Astrophysics Data System (ADS)

Holland, David M.; Mysak, Lawrence A.; Manak, Davinder K.; Oberhuber, Josef M.

1993-02-01

A numerical simulation of the seasonal sea ice cover in the Arctic Ocean and the Greenland, Iceland, and Norwegian seas is presented. The sea ice model is extracted from Oberhuber's (1990) coupled sea ice-mixed layer-isopycnal general circulation model and is written in spherical coordinates. The advantage of such a model over previous sea ice models is that it can be easily coupled to either global atmospheric or ocean general circulation models written in spherical coordinates. In this model, the thermodynamics are a modification of that of Parkinson and Washington (1979), while the dynamics use the full Hibler (1979) viscous-plastic rheology. Monthly thermodynamic and dynamic forcing fields for the atmosphere and ocean are specified. The simulations of the seasonal cycle of ice thickness, compactness, and velocity, for a control set of parameters, compare favorably with the known seasonal characteristics of these fields. A sensitivity study of the control simulation of the seasonal sea ice cover is presented. The sensitivity runs are carried out under three different themes, namely, numerical conditions, parameter values, and physical processes. This last theme refers to experiments in which physical processes are either newly added or completely removed from the model. Approximately 80 sensitivity runs have been performed in which a change from the control run environment has been implemented. Comparisons have been made between the control run and a particular sensitivity run based on time series of the seasonal cycle of the domain-averaged ice thickness, compactness, areal coverage, and kinetic energy. In addition, spatially varying fields of ice thickness, compactness, velocity, and surface temperature for each season are presented for selected experiments. A brief description and discussion of the more interesting experiments are presented. The simulation of the seasonal cycle of Arctic sea ice cover is shown to be robust.

Stationary Solutions of A One-dimensional Thermodynamic Radiative Sea Ice Model

NASA Astrophysics Data System (ADS)

Taylor, P. D.; Feltham, D. L.

A one-dimensional thermodynamic model of sea ice is coupled to a two-stream radi- ation model and the stationary (time-independent) solutions analysed. The stationary model represents the state of the sea ice subjected to persistent or slowly varying forc- ing. Two physically realisable stationary solutions (real and positive ice thickness) occur for a large range of positive oceanic heat flux ( 20,Wm-2). The two station- ary solutions are due to the two-stream radiation model, which allows radiation to be reflected at the ice-ocean interface. Thick ice ( 1,m) only absorbs radiation near its surface, whereas thin ice ( 0.1,m) absorbs radiation across its entire depth. The two stationary solutions are caused by these two different radiative regimes. The results of this analysis have relevance to the interpretation and implementation of thermody- namic models of sea ice and the interpretation of thickness data.

[Spectral features analysis of sea ice in the Arctic Ocean].

PubMed

Ke, Chang-qing; Xie, Hong-jie; Lei, Rui-bo; Li, Qun; Sun, Bo

2012-04-01

Sea ice in the Arctic Ocean plays an important role in the global climate change, and its quick change and impact are the scientists' focus all over the world. The spectra of different kinds of sea ice were measured with portable ASD FieldSpec 3 spectrometer during the long-term ice station of the 4th Chinese national Arctic Expedition in 2010, and the spectral features were analyzed systematically. The results indicated that the reflectance of sea ice covered by snow is the highest one, naked sea ice the second, and melted sea ice the lowest. Peak and valley characteristics of spectrum curves of sea ice covered by thick snow, thin snow, wet snow and snow crystal are very significant, and the reflectance basically decreases with the wavelength increasing. The rules of reflectance change with wavelength of natural sea ice, white ice and blue ice are basically same, the reflectance of them is medium, and that of grey ice is far lower than natural sea ice, white ice and blue ice. It is very significant for scientific research to analyze the spectral features of sea ice in the Arctic Ocean and to implement the quantitative remote sensing of sea ice, and to further analyze its response to the global warming.

Exploration of the Climate Change Frontier in Polar Regions at the Land Ice-Ocean Boundary.

NASA Astrophysics Data System (ADS)

Rignot, E. J.

2014-12-01

Ice sheets are the largest contributors to sea level rise at present, and responsible for the largest uncertainty in sea level projections. Ice sheets raised sea level 5 m per century 13.5 kyr ago during one period of rapid change. Leading regions for future rapid changes include the marine-based, retrograde bed parts of Greenland (north center and east), West Antarctica (Amundsen Sea), and East Antarctica (Filchner basin and Wilkes Land). Fast changes require an increase in ice melt from a warmer ocean and an increase in iceberg calving. Our understanding of both processes remains limited due to a lack of basic observations. Understanding ocean forcing requires observations on the continental shelf, along bays and glacial fjords and at ice-ocean boundaries, beneath kilometers of ice (Antarctica) or at near-vertical calving cliffs (Greenland), of ocean temperature and sea floor bathymetry. Where such observations exist, the sea floor is much deeper than anticipated because of the carving of deep channels by multiple glacier advances. Warm subsurface waters penetrate throughout the Amundsen Sea Embayment of West Antarctica, the southeast and probably the entire west coasts of Greenland. In Greenland, discharge of subglacial water from surface runoff at the glacier grounding line increases ice melting by the ocean even if the ocean temperature remains the same. Near ice-ocean boundaries, satellite observations are challenged, airborne observations and field surveys are limited, so advanced robotic techniques for cold, deep, remote environments are ultimately required in combination with advanced numerical modeling techniques. Until such technological advances take place and advanced networks are put in place, it is critical to conduct boat surveys, install moorings, and conduct extensive airborne campaigns (for instance, gravity-derived bathymetry and air-dropped CTDs), some of which is already taking place. In the meantime, projections of ice sheet evolution in a

The seasonal and inter-annual variability of sea-ice, ocean circulation and marine ecosystems in the Barents Sea: model results against satellite data

NASA Astrophysics Data System (ADS)

Dvornikov, Anton; Sein, Dmitry; Ryabchenko, Vladimir; Gorchakov, Victor; Pugalova, Svetlana

2015-04-01

This study is aimed at modelling the seasonal and inter-annual variability of sea-ice, ocean circulation and marine ecosystems in the Barents Sea in the modern period. Adequate description of marine ecosystems in the ice-covered seas crucially depends on the accuracy in determining of thicknesses of ice and snow on the sea surface which control penetrating photosynthetically active radiation under the ice. One of the few models of ice able to adequately reproduce the dynamics of sea ice is the sea ice model HELMI [1], containing 7 different categories of ice. This model has been imbedded into the Princeton Ocean Model. With this coupled model 2 runs for the period 1998-2007 were performed under different atmospheric forcing prescribed from NCEP/NCAR and ERA-40 archives. For prescribing conditions at the open boundary, all the necessary information about the horizontal velocity, level, temperature and salinity of the water, ice thickness and compactness was taken from the results of the global ocean general circulation model of the Max Planck Institute for Meteorology (Hamburg, Germany) MPIOM [2]. The resulting solution with NCEP forcing with a high accuracy simulates the seasonal and inter-annual variability of sea surface temperature (SST) estimated from MODIS data. The maximum difference between the calculated and satellite-derived SSTs (averaged over 4 selected areas of the Barents Sea) during the period 2000-2007 does not exceed 1.5 °C. Seasonal and inter-annual variations in the area of ice cover are also in good agreement with satellite-derived estimates. Pelagic ecosystem model developed in [3] has been coupled into the above hydrodynamic model and used to calculate the changes in the characteristics of marine ecosystems under NCEP forcing. Preliminarily the ecosystem model has been improved by introducing a parameterization of detritus deposition on the bottom and through the selection of optimal parameters for photosynthesis and zooplankton grazing

Sea-ice evaluation of NEMO-Nordic 1.0: a NEMO-LIM3.6-based ocean-sea-ice model setup for the North Sea and Baltic Sea

NASA Astrophysics Data System (ADS)

Pemberton, Per; Löptien, Ulrike; Hordoir, Robinson; Höglund, Anders; Schimanke, Semjon; Axell, Lars; Haapala, Jari

2017-08-01

The Baltic Sea is a seasonally ice-covered marginal sea in northern Europe with intense wintertime ship traffic and a sensitive ecosystem. Understanding and modeling the evolution of the sea-ice pack is important for climate effect studies and forecasting purposes. Here we present and evaluate the sea-ice component of a new NEMO-LIM3.6-based ocean-sea-ice setup for the North Sea and Baltic Sea region (NEMO-Nordic). The setup includes a new depth-based fast-ice parametrization for the Baltic Sea. The evaluation focuses on long-term statistics, from a 45-year long hindcast, although short-term daily performance is also briefly evaluated. We show that NEMO-Nordic is well suited for simulating the mean sea-ice extent, concentration, and thickness as compared to the best available observational data set. The variability of the annual maximum Baltic Sea ice extent is well in line with the observations, but the 1961-2006 trend is underestimated. Capturing the correct ice thickness distribution is more challenging. Based on the simulated ice thickness distribution we estimate the undeformed and deformed ice thickness and concentration in the Baltic Sea, which compares reasonably well with observations.

Influence of Sea Ice on the Thermohaline Circulation in the Arctic-North Atlantic Ocean

NASA Technical Reports Server (NTRS)

Mauritzen, Cecilie; Haekkinen, Sirpa

1997-01-01

A fully prognostic coupled ocean-ice model is used to study the sensitivity of the overturning cell of the Arctic-North-Atlantic system to sea ice forcing. The strength of the thermohaline cell will be shown to depend on the amount of sea ice transported from the Arctic to the Greenland Sea and further to the subpolar gyre. The model produces a 2-3 Sv increase of the meridional circulation cell at 25N (at the simulation year 15) corresponding to a decrease of 800 cu km in the sea ice export from the Arctic. Previous modeling studies suggest that interannual and decadal variability in sea ice export of this magnitude is realistic, implying that sea ice induced variability in the overturning cell can reach 5-6 Sv from peak to peak.

Accuracy Assessment of Recent Global Ocean Tide Models around Antarctica

NASA Astrophysics Data System (ADS)

Lei, J.; Li, F.; Zhang, S.; Ke, H.; Zhang, Q.; Li, W.

2017-09-01

Due to the coverage limitation of T/P-series altimeters, the lack of bathymetric data under large ice shelves, and the inaccurate definitions of coastlines and grounding lines, the accuracy of ocean tide models around Antarctica is poorer than those in deep oceans. Using tidal measurements from tide gauges, gravimetric data and GPS records, the accuracy of seven state-of-the-art global ocean tide models (DTU10, EOT11a, GOT4.8, FES2012, FES2014, HAMTIDE12, TPXO8) is assessed, as well as the most widely-used conventional model FES2004. Four regions (Antarctic Peninsula region, Amery ice shelf region, Filchner-Ronne ice shelf region and Ross ice shelf region) are separately reported. The standard deviations of eight main constituents between the selected models are large in polar regions, especially under the big ice shelves, suggesting that the uncertainty in these regions remain large. Comparisons with in situ tidal measurements show that the most accurate model is TPXO8, and all models show worst performance in Weddell sea and Filchner-Ronne ice shelf regions. The accuracy of tidal predictions around Antarctica is gradually improving.

Impact of increasing antarctic glacial freshwater release on regional sea-ice cover in the Southern Ocean

NASA Astrophysics Data System (ADS)

Merino, Nacho; Jourdain, Nicolas C.; Le Sommer, Julien; Goosse, Hugues; Mathiot, Pierre; Durand, Gael

2018-01-01

The sensitivity of Antarctic sea-ice to increasing glacial freshwater release into the Southern Ocean is studied in a series of 31-year ocean/sea-ice/iceberg model simulations. Glaciological estimates of ice-shelf melting and iceberg calving are used to better constrain the spatial distribution and magnitude of freshwater forcing around Antarctica. Two scenarios of glacial freshwater forcing have been designed to account for a decadal perturbation in glacial freshwater release to the Southern Ocean. For the first time, this perturbation explicitly takes into consideration the spatial distribution of changes in the volume of Antarctic ice shelves, which is found to be a key component of changes in freshwater release. In addition, glacial freshwater-induced changes in sea ice are compared to typical changes induced by the decadal evolution of atmospheric states. Our results show that, in general, the increase in glacial freshwater release increases Antarctic sea ice extent. But the response is opposite in some regions like the coastal Amundsen Sea, implying that distinct physical mechanisms are involved in the response. We also show that changes in freshwater forcing may induce large changes in sea-ice thickness, explaining about one half of the total change due to the combination of atmospheric and freshwater changes. The regional contrasts in our results suggest a need for improving the representation of freshwater sources and their evolution in climate models.

Drivers of Antarctic sea-ice expansion and Southern Ocean surface cooling over the past four decades

NASA Astrophysics Data System (ADS)

Purich, Ariaan; England, Matthew

2017-04-01

Despite global warming, total Antarctic sea-ice coverage has increased overall during the past four decades. In contrast, the majority of CMIP5 models simulate a decline. In addition, Southern Ocean surface waters have largely cooled, in stark contrast to almost all historical CMIP5 simulations. Subantarctic Surface Waters have cooled and freshened while waters to the north of the Antarctic Circumpolar Current have warmed and increased in salinity. It remains unclear as to what extent the cooling and Antarctic sea-ice expansion is due to natural variability versus anthropogenic forcing; due for example to changes in the Southern Annular Mode (SAM). It is also unclear what the respective role of surface buoyancy fluxes is compared to internal ocean circulation changes, and what the implications are for longer-term climate change in the region. In this presentation we will outline three distinct drivers of recent Southern Ocean surface trends that have each made a significant contribution to regional cooling: (1) wind-driven surface cooling and sea-ice expansion due to shifted westerly winds, (2) teleconnections of decadal variability from the tropical Pacific, and (3) surface cooling and ice expansion due to large-scale Southern Ocean freshening, most likely driven by SAM-related precipitation trends over the open ocean. We will also outline the main reasons why climate models for the most part miss these Southern Ocean cooling trends, despite capturing overall trends in the SAM.

Measurements to Fill Knowledge Gaps on Ice Nucleating Particle Sources over Oceans

NASA Astrophysics Data System (ADS)

DeMott, P. J.; Hill, T. C.; Ruppel, M. J.; Prather, K. A.; Collins, D. B.; Axson, J. L.; Lee, T.; Hwang, C. Y.; Sullivan, R. C.; McMeeking, G. R.; Mason, R.; Bertram, A. K.; Mayol-Bracero, O. L.; Lewis, E. R.

2013-12-01

Measurements of the temperature spectrum of ice nucleating particle concentrations by two methods in recent specialized laboratory sea spray studies and field campaigns in the Northern Hemisphere will be discussed and compared with historical data from over Southern Oceans. In general, new measurements of the condensation/immersion freezing activation spectra of realistically-generated laboratory sea spray particles (by wave generation or plunging water bubble production) are consistent with previous measurements made over oceans. The number concentrations of ice nuclei tend to be lower than are measured over land regions, at least for modestly supercooled cloud conditions. Certain but complex connections of ice nucleating particle production to ocean microbiological processes affecting the chemical composition of the sea surface microlayer are seen, but the nature of the ice nucleating units of particles remains to be identified. Associations of ice nucleating particle concentrations with heterotrophic bacterial concentrations were noted in some experiments, while correlation with chlorophyll-a concentration in seawater was clearly identified in laboratory simulations of phytoplankton blooms. These data may ultimately serve as the basis for parameterization development for ice initiation in numerical model simulations of mixed-phase clouds. Atmospheric measurements have been made at island sites, via aircraft, and from ship-based filter collections in the Northern Hemisphere. The immersion freezing spectra of these particles are similar to those found in recent laboratory studies and historical measurements, but show the expected natural variability by location. The majority of particles detected thus far as ice nuclei from sea spray and in marine air show minimal or episodic/variable direct participation of biological ice nucleating organisms on the basis of sensitivity to high temperatures (95°C). However, assembled measurements are still sparse, the nuclei

Characteristics of the ocean simulations in the Max Planck Institute Ocean Model (MPIOM) the ocean component of the MPI-Earth system model

NASA Astrophysics Data System (ADS)

Jungclaus, J. H.; Fischer, N.; Haak, H.; Lohmann, K.; Marotzke, J.; Matei, D.; Mikolajewicz, U.; Notz, D.; von Storch, J. S.

2013-06-01

MPI-ESM is a new version of the global Earth system model developed at the Max Planck Institute for Meteorology. This paper describes the ocean state and circulation as well as basic aspects of variability in simulations contributing to the fifth phase of the Coupled Model Intercomparison Project (CMIP5). The performance of the ocean/sea-ice model MPIOM, coupled to a new version of the atmosphere model ECHAM6 and modules for land surface and ocean biogeochemistry, is assessed for two model versions with different grid resolution in the ocean. The low-resolution configuration has a nominal resolution of 1.5°, whereas the higher resolution version features a quasiuniform, eddy-permitting global resolution of 0.4°. The paper focuses on important oceanic features, such as surface temperature and salinity, water mass distribution, large-scale circulation, and heat and freshwater transports. In general, these integral quantities are simulated well in comparison with observational estimates, and improvements in comparison with the predecessor system are documented; for example, for tropical variability and sea ice representation. Introducing an eddy-permitting grid configuration in the ocean leads to improvements, in particular, in the representation of interior water mass properties in the Atlantic and in the representation of important ocean currents, such as the Agulhas and Equatorial current systems. In general, however, there are more similarities than differences between the two grid configurations, and several shortcomings, known from earlier versions of the coupled model, prevail.

The effect of sudden ice sheet melt on ocean circulation and surface climate

NASA Astrophysics Data System (ADS)

Ivanovic, R. F.; Gregoire, L. J.; Wickert, A. D.; Valdes, P. J.; Burke, A.

2017-12-01

Collapse of ice sheets can cause significant sea-level rise and widespread climate change. Around 14.6 thousand years ago, global mean sea level rose by 15 m in less than 350 years during an event known as Meltwater Pulse 1a. Ice sheet modelling and sea-level fingerprinting has suggested that approximately half of this 50 mm yr-1 sea level rise may have come from a North American ice Saddle Collapse that drained into the Arctic and Atlantic Oceans. However, dating uncertainties make it difficult to determine the sequence of events and their drivers, leaving many fundamental questions. For example, was melting from the northern ice sheets responsible for the Older-Dryas or other global-scale cooling events, or did a contribution from Antarctica counteract the climatic effects? What was the role of the abrupt Bølling Warming? And how were all these signals linked to changes in Atlantic Ocean overturning circulation?To address these questions, we examined the effect of the North American ice Saddle Collapse using a high resolution network drainage model coupled to an atmosphere-ocean-vegetation General Circulation Model. Here, we present the quantitative routing estimates of the consequent meltwater discharge and its impact on climate. We also tested a suite of more idealised meltwater forcing scenarios to examine the global influence of Arctic versus Antarctic ice melt. The results show that 50% of the Saddle Collapse meltwater pulse was routed via the Mackenzie River into the Arctic Ocean, and 50% was discharged directly into the Atlantic/Gulf of Mexico. This meltwater flux, equivalent to a total of 7.3 m of sea-level rise, caused a strong (6 Sv) weakening of Atlantic Meridional Overturning Circulation (AMOC) and widespread Northern Hemisphere cooling of 1-5 °C. The greatest cooling is in the Arctic (5-10 °C in the winter), but there is also significant winter warming over eastern North America (1-3 °C). We propose that this robust submillennial mechanism was

Strain in shore fast ice due to incoming ocean waves and swell

NASA Astrophysics Data System (ADS)

Fox, Colin; Squire, Vernon A.

1991-03-01

Using a development from the theoretical model presented by Fox and Squire (1990), this paper investigates the strain field generated in shore fast ice by normally incident ocean waves and swell. After a brief description of the model and its convergence, normalized absolute strain (relative to a 1-m incident wave) is found as a function of distance from the ice edge for various wave periods, ice thicknesses, and water depths. The squared transfer function, giving the relative ability of incident waves of different periods to generate strain in the ice, is calculated, and its consequences are discussed. The ice is then forced with a Pierson-Moskowitz spectrum, and the consequent strain spectra are plotted as a function of penetration into the ice sheet. Finally, rms strain, computed as the incoherent sum of the strains resulting from energy in the open water spectrum, is found. The results have implications to the breakup of shore fast ice and hence to the floe size distribution of the marginal ice zone.

Arctic Sea Ice Classification and Mapping for Surface Albedo Parameterization in Sea Ice Modeling

NASA Astrophysics Data System (ADS)

Nghiem, S. V.; Clemente-Colón, P.; Perovich, D. K.; Polashenski, C.; Simpson, W. R.; Rigor, I. G.; Woods, J. E.; Nguyen, D. T.; Neumann, G.

2016-12-01

A regime shift of Arctic sea ice from predominantly perennial sea ice (multi-year ice or MYI) to seasonal sea ice (first-year ice or FYI) has occurred in recent decades. This shift has profoundly altered the proportional composition of different sea ice classes and the surface albedo distribution pertaining to each sea ice class. Such changes impacts physical, chemical, and biological processes in the Arctic atmosphere-ice-ocean system. The drastic changes upset the traditional geophysical representation of surface albedo of the Arctic sea ice cover in current models. A critical science issue is that these profound changes must be rigorously and systematically observed and characterized to enable a transformative re-parameterization of key model inputs, such as ice surface albedo, to ice-ocean-atmosphere climate modeling in order to obtain re-analyses that accurately reproduce Arctic changes and also to improve sea ice and weather forecast models. Addressing this challenge is a strategy identified by the National Research Council study on "Seasonal to Decadal Predictions of Arctic Sea Ice - Challenges and Strategies" to replicate the new Arctic reality. We review results of albedo characteristics associated with different sea ice classes such as FYI and MYI. Then we demonstrate the capability for sea ice classification and mapping using algorithms developed by the Jet Propulsion Laboratory and by the U.S. National Ice Center for use with multi-sourced satellite radar data at L, C, and Ku bands. Results obtained with independent algorithms for different radar frequencies consistently identify sea ice classes and thereby cross-verify the sea ice classification methods. Moreover, field observations obtained from buoy webcams and along an extensive trek across Elson Lagoon and a sector of the Beaufort Sea during the BRomine, Ozone, and Mercury EXperiment (BROMEX) in March 2012 are used to validate satellite products of sea ice classes. This research enables the mapping

Sea ice, extremophiles and life on extra-terrestrial ocean worlds

NASA Astrophysics Data System (ADS)

Martin, Andrew; McMinn, Andrew

2018-01-01

The primary aim of this review is to highlight that sea-ice microbes would be capable of occupying ice-associated biological niches on Europa and Enceladus. These moons are compelling targets for astrobiological exploration because of the inferred presence of subsurface oceans that have persisted over geological timescales. Although potentially hostile to life in general, Europa and Enceladus may still harbour biologically permissive domains associated with the ice, ocean and seafloor environments. However, validating sources of free energy is challenging, as is qualifying possible metabolic processes or ecosystem dynamics. Here, the capacity for biological adaptation exhibited by microorganisms that inhabit sea ice is reviewed. These ecosystems are among the most relevant Earth-based analogues for considering life on ocean worlds because microorganisms must adapt to multiple physicochemical extremes. In future, these organisms will likely play a significant role in defining the constraints on habitability beyond Earth and developing a mechanistic framework that contrasts the limits of Earth's biosphere with extra-terrestrial environments of interest.

Mounting evidence for intense ocean interaction with the Pine Island Glacier Ice Shelf

NASA Astrophysics Data System (ADS)

Bindschadler, R.; Holland, D.; Vaughan, D.; Vornberger, P.

2008-12-01

The spatial signature of thinning and acceleration of the Pine Island Glacier has led to the inference that these changes originate at the seaward end of the glacier, possibly within or under the ice shelf (Payne et al., 2004; Shepherd et al., 2004). We present new analyses resulting from both new and archived satellite imagery of the ice shelf that supports this inference and provides new insights into strong seasonal and intra- annual characters of ocean-ice shelf interaction. Strong longitudinal variations in both thickness and surface elevation measured by British Antarctic Survey airborne radars (Vaughan et al., 2006) have wavelengths that correspond roughly to the annual motion of the ice shelf. These could be caused by seasonal variations in flow speed, but such variations of flow speed have never been reported and are not seen in the most recent continuous GPS observations of the ice shelf. We suggest that these strong variations in ice thickness, as large as 200 meters in an average thickness of 600 meters, are caused by seasonal variations in the properties of the water circulating underneath the ice shelf. One likely explanation is that the dominant water mass reaching the deepest parts of the ice shelf alternates between cold High Salinity Shelf Water in the winter and warm Circumpolar Deep Water in the summer. Evidence for recent strengthening of the sub- shelf circulation is the sudden occurrence of three persistent polynyas immediately adjacent to the ice front. These are located in precisely the locations expected from modeled sub-shelf circulation (Payne et al., 2007). This mode was never observed in any satellite imagery prior to the 1999-2000 austral summer (data of 7 summers since 1973 were available), but has occurred in 7 of the 9 summers since and persists throughout the summer. Payne, A.J., A. Vieli, A.P. Shepherd, D.J. Wingham and E. Rignot, 2004. Recent dramatic thinning of largest West Antarctic ice stream triggered by oceans, Geophysical

Observations and modeling of ocean-induced melt beneath Petermann Glacier Ice Shelf in northwestern Greenland

NASA Astrophysics Data System (ADS)

Cai, Cilan; Rignot, Eric; Menemenlis, Dimitris; Nakayama, Yoshihiro

2017-08-01

We update observationally based estimates of subaqueous melt, Qm, beneath Petermann Glacier Ice Shelf (PGIS), Greenland, and model its sensitivity to oceanic thermal forcing, TF, and subglacial runoff, Qsg, using the Massachusetts Institute of Technology general circulation model (MITgcm), in a two-dimensional domain, with 20 m vertical and 40 m horizontal resolution at the grounding line. We adjust the drag coefficient to match the observationally based Qm. With the inclusion of Qsg, the maximum melt rate (Qmmax) is 2 times larger in summer and 1/3 larger annually than in winter. Qmmax increases above linear with TF and below linear with Qsg. We estimate that Qmmax increased by 24% (+8.1 m/yr) beneath PGIS from the 1990s to the 2000s from a 0.21°C warming in ocean temperature and a doubling in Qsg, hence contributing to its thinning. If the PGIS is removed, we estimate that the modeled melt rate near the grounding line will increase 13-16 times.

Modeling Vertical Structure and Heat Transport within the Oceans of Ice-covered Worlds (Invited)

NASA Astrophysics Data System (ADS)

Goodman, J. C.

2010-12-01

Indirect observational evidence provides a strong case for liquid oceans beneath the icy crust of Europa and several other frozen moons in the outer solar system. However, little is known about the fluid circulation within these exotic oceans. As a first step toward understanding circulations driven by buoyancy (rather than mechanical forcing from tides), one must understand the typical vertical structure of temperature, salinity, and thus density within the ocean. Following a common approach from terrestrial oceanography, I have built a "single column convection model" for icy world oceans, which describes the density structure of the ocean as a function of depth only: horizontal variations are ignored. On Earth, this approach is of limited utility, because of the strong influence of horizontal wind-driven currents and sea-surface temperature gradients set in concert with the overlying atmosphere. Neither of these confounding issues is present in an icy world's ocean. In the model, mixing of fluid properties via overturning convection is modeled as a strong diffusive process which only acts when the ocean is vertically unstable. "Double diffusive" processes (salt fingering and diffusive layering) are included: these are mixing processes resulting from the unequal molecular diffusivities of heat and salt. Other important processes, such as heating on adiabatic compression, and freshwater fluxes from melting overlying ice, are also included. As a simple test case, I considered an ocean of Europa-like depth (~100 km) and gravity, heated from the seafloor. To simplify matters, I specified an equation of state appropriate to terrestrial seawater, and a simple isothermal ocean as an initial condition. As expected, convection gradually penetrates upward, warming the ocean to an adiabatic, unstratified equilibrium density profile on a timescale of 50 kyr if 4.5 TW of heat are emitted by the silicate interior; the same result is achieved in proportionally more/less time

Large-Ensemble modeling of past and future variations of the Antarctic Ice Sheet with a coupled ice-Earth-sea level model

NASA Astrophysics Data System (ADS)

Pollard, David; DeConto, Robert; Gomez, Natalya

2016-04-01

To date, most modeling of the Antarctic Ice Sheet's response to future warming has been calibrated using recent and modern observations. As an alternate approach, we apply a hybrid 3-D ice sheet-shelf model to the last deglacial retreat of Antarctica, making use of geologic data of the last ~20,000 years to test the model against the large-scale variations during this period. The ice model is coupled to a global Earth-sea level model to improve modeling of the bedrock response and to capture ocean-ice gravitational interactions. Following several recent ice-sheet studies, we use Large Ensemble (LE) statistical methods, performing sets of 625 runs from 30,000 years to present with systematically varying model parameters. Objective scores for each run are calculated using modern data and past reconstructed grounding lines, relative sea level records, cosmogenic elevation-age data and uplift rates. The LE results are analyzed to calibrate 4 particularly uncertain model parameters that concern marginal ice processes and interaction with the ocean. LE's are extended into the future with climates following RCP scenarios. An additional scoring criterion tests the model's ability to reproduce estimated sea-level high stands in the warm mid-Pliocene, for which drastic retreat mechanisms of hydrofracturing and ice-cliff failure are needed in the model. The LE analysis provides future sea-level-rise envelopes with well-defined parametric uncertainty bounds. Sensitivities of future LE results to Pliocene sea-level estimates, coupling to the Earth-sea level model, and vertical profiles of Earth properties, will be presented.

Sea-ice deformation in a coupled ocean-sea-ice model and in satellite remote sensing data

NASA Astrophysics Data System (ADS)

Spreen, Gunnar; Kwok, Ron; Menemenlis, Dimitris; Nguyen, An T.

2017-07-01

A realistic representation of sea-ice deformation in models is important for accurate simulation of the sea-ice mass balance. Simulated sea-ice deformation from numerical simulations with 4.5, 9, and 18 km horizontal grid spacing and a viscous-plastic (VP) sea-ice rheology are compared with synthetic aperture radar (SAR) satellite observations (RGPS, RADARSAT Geophysical Processor System) for the time period 1996-2008. All three simulations can reproduce the large-scale ice deformation patterns, but small-scale sea-ice deformations and linear kinematic features (LKFs) are not adequately reproduced. The mean sea-ice total deformation rate is about 40 % lower in all model solutions than in the satellite observations, especially in the seasonal sea-ice zone. A decrease in model grid spacing, however, produces a higher density and more localized ice deformation features. The 4.5 km simulation produces some linear kinematic features, but not with the right frequency. The dependence on length scale and probability density functions (PDFs) of absolute divergence and shear for all three model solutions show a power-law scaling behavior similar to RGPS observations, contrary to what was found in some previous studies. Overall, the 4.5 km simulation produces the most realistic divergence, vorticity, and shear when compared with RGPS data. This study provides an evaluation of high and coarse-resolution viscous-plastic sea-ice simulations based on spatial distribution, time series, and power-law scaling metrics.

Evidence for ice-free summers in the late Miocene central Arctic Ocean

PubMed Central

Stein, Ruediger; Fahl, Kirsten; Schreck, Michael; Knorr, Gregor; Niessen, Frank; Forwick, Matthias; Gebhardt, Catalina; Jensen, Laura; Kaminski, Michael; Kopf, Achim; Matthiessen, Jens; Jokat, Wilfried; Lohmann, Gerrit

2016-01-01

Although the permanently to seasonally ice-covered Arctic Ocean is a unique and sensitive component in the Earth's climate system, the knowledge of its long-term climate history remains very limited due to the restricted number of pre-Quaternary sedimentary records. During Polarstern Expedition PS87/2014, we discovered multiple submarine landslides along Lomonosov Ridge. Removal of younger sediments from steep headwalls has led to exhumation of Miocene sediments close to the seafloor. Here we document the presence of IP25 as a proxy for spring sea-ice cover and alkenone-based summer sea-surface temperatures >4 °C that support a seasonal sea-ice cover with an ice-free summer season being predominant during the late Miocene in the central Arctic Ocean. A comparison of our proxy data with Miocene climate simulations seems to favour either relatively high late Miocene atmospheric CO2 concentrations and/or a weak sensitivity of the model to simulate the magnitude of high-latitude warming in a warmer than modern climate. PMID:27041737

Ocean as the main driver of Antarctic ice sheet retreat during the Holocene

NASA Astrophysics Data System (ADS)

Crosta, Xavier; Crespin, Julien; Swingedouw, Didier; Marti, Olivier; Masson-Delmotte, Valérie; Etourneau, Johan; Goosse, Hugues; Braconnot, Pascale; Yam, Ruth; Brailovski, Irena; Shemesh, Aldo

2018-07-01

Ocean-driven basal melting has been shown to be the main ablation process responsible for the recession of many Antarctic ice shelves and marine-terminating glaciers over the last decades. However, much less is known about the drivers of ice shelf melt prior to the short instrumental era. Based on diatom oxygen isotope (δ18Odiatom; a proxy for glacial ice discharge in solid or liquid form) records from western Antarctic Peninsula (West Antarctica) and Adélie Land (East Antarctica), higher ocean temperatures were suggested to have been the main driver of enhanced ice melt during the Early-to-Mid Holocene while atmosphere temperatures were proposed to have been the main driver during the Late Holocene. Here, we present a new Holocene δ18Odiatom record from Prydz Bay, East Antarctica, also suggesting an increase in glacial ice discharge since 4500 years before present ( 4.5 kyr BP) as previously observed in Antarctic Peninsula and Adélie Land. Similar results from three different regions around Antarctica thus suggest common driving mechanisms. Combining marine and ice core records along with new transient accelerated simulations from the IPSL-CM5A-LR climate model, we rule out changes in air temperatures during the last 4.5 kyr as the main driver of enhanced glacial ice discharge. Conversely, our simulations evidence the potential for significant warmer subsurface waters in the Southern Ocean during the last 6 kyr in response to enhanced summer insolation south of 60°S and enhanced upwelling of Circumpolar Deep Water towards the Antarctic shelf. We conclude that ice front and basal melting may have played a dominant role in glacial discharge during the Late Holocene.

Aragonite undersaturation in the Arctic Ocean: effects of ocean acidification and sea ice melt.

PubMed

Yamamoto-Kawai, Michiyo; McLaughlin, Fiona A; Carmack, Eddy C; Nishino, Shigeto; Shimada, Koji

2009-11-20

The increase in anthropogenic carbon dioxide emissions and attendant increase in ocean acidification and sea ice melt act together to decrease the saturation state of calcium carbonate in the Canada Basin of the Arctic Ocean. In 2008, surface waters were undersaturated with respect to aragonite, a relatively soluble form of calcium carbonate found in plankton and invertebrates. Undersaturation was found to be a direct consequence of the recent extensive melting of sea ice in the Canada Basin. In addition, the retreat of the ice edge well past the shelf-break has produced conditions favorable to enhanced upwelling of subsurface, aragonite-undersaturated water onto the Arctic continental shelf. Undersaturation will affect both planktonic and benthic calcifying biota and therefore the composition of the Arctic ecosystem.

Sea ice and oceanic processes on the Ross Sea continental shelf

NASA Astrophysics Data System (ADS)

Jacobs, S. S.; Comiso, J. C.

1989-12-01

We have investigated the spatial and temporal variability of Antarctic sea ice concentrations on the Ross Sea continental shelf, in relation to oceanic and atmospheric forcing. Sea ice data were derived from Nimbus 7 scanning multichannel microwave radiometer (SMMR) brightness temperatures from 1979-1986. Ice cover over the shelf was persistently lower than above the adjacent deep ocean, averaging 86% during winter with little month-to-month or interannual variability. The large spring Ross Sea polynya on the western shelf results in a longer period of summer insolation, greater surface layer heat storage, and later ice formation in that region the following autumn. Newly identified Pennell and Ross Passage polynyas near the continental shelf break appear to be maintained in part by divergence above a submarine bank and by upwelling of warmer water near the slope front. Warmer subsurface water enters the shelf region year-round and will retard ice growth and enhance heat flux to the atmosphere when entrained in the strong winter vertical circulation. Temperatures at 125-m depth on a mooring near the Ross Ice Shelf during July 1984 averaged 0.15°C above freezing, sufficient to support a vertical heat flux above 100 W/m2. Monthly average subsurface ocean temperatures along the Ross Ice Shelf lag the air temperature cycle and begin to rise several weeks before spring ice breakout. The coarse SMMR resolution and dynamic ice shelf coastlines can compromise the use of microwave sea ice data near continental boundaries.

Ice-tethered measurement platforms in the Arctic Ocean: a contribution by the FRAM infrastructure program

NASA Astrophysics Data System (ADS)

Hoppmann, Mario; Nicolaus, Marcel; Rabe, Benjamin; Wenzhöfer, Frank; Katlein, Christian; Scholz, Daniel

2016-04-01

additional data acquired by FRAM buoys facilitate the validation of model results and remote sensing products, play an important role in understanding the linkages between the atmosphere, sea ice and upper ocean, and help assess the physical, biological and biogeochemical states of the future Arctic Ocean. Here we present our recent work and future plans, but are also aiming for additional collaborations, especially on technical developments, scientific questions and deployment logistics.

Spatial and Temporal Variability of Surface Energy Fluxes During Autumn Ice Advance: Observations and Model Validation

NASA Astrophysics Data System (ADS)

Persson, O. P. G.; Blomquist, B.; Grachev, A. A.; Guest, P. S.; Stammerjohn, S. E.; Solomon, A.; Cox, C. J.; Capotondi, A.; Fairall, C. W.; Intrieri, J. M.

2016-12-01

From Oct 4 to Nov 5, 2015, the Office of Naval Research - sponsored Sea State cruise in the Beaufort Sea with the new National Science Foundation R/V Sikuliaq obtained extensive in-situ and remote sensing observations of the lower troposphere, the advancing sea ice, wave state, and upper ocean conditions. In addition, a coupled atmosphere, sea ice, upper-ocean model, based on the RASM model, was run at NOAA/PSD in a hindcast mode for this same time period, providing a 10-day simulation of the atmosphere/ice/ocean evolution. Surface energy fluxes quantitatively represent the air-ice, air-ocean, and ice-ocean interaction processes, determining the cooling (warming) rate of the upper ocean and the growth (melting) rate of sea ice. These fluxes also impact the stratification of the lower troposphere and the upper ocean. In this presentation, both direct and indirect measurements of the energy fluxes during Sea State will be used to explore the spatial and temporal variability of these fluxes and the impacts of this variability on the upper ocean, ice, and lower atmosphere during the autumn ice advance. Analyses have suggested that these fluxes are impacted by atmospheric synoptic evolution, proximity to existing ice, ice-relative wind direction, ice thickness and snow depth. In turn, these fluxes impact upper-ocean heat loss and timing of ice formation, as well as stability in the lower troposphere and upper ocean, and hence heat transport to the free troposphere and ocean mixed-layer. Therefore, the atmospheric structure over the advancing first-year ice differs from that over the nearby open water. Finally, these observational analyses will be used to provide a preliminary validation of the spatial and temporal variability of the surface energy fluxes and the associated lower-tropospheric and upper-ocean structures in the simulations.

Challenges in validating model results for first year ice

NASA Astrophysics Data System (ADS)

Melsom, Arne; Eastwood, Steinar; Xie, Jiping; Aaboe, Signe; Bertino, Laurent

2017-04-01

In order to assess the quality of model results for the distribution of first year ice, a comparison with a product based on observations from satellite-borne instruments has been performed. Such a comparison is not straightforward due to the contrasting algorithms that are used in the model product and the remote sensing product. The implementation of the validation is discussed in light of the differences between this set of products, and validation results are presented. The model product is the daily updated 10-day forecast from the Arctic Monitoring and Forecasting Centre in CMEMS. The forecasts are produced with the assimilative ocean prediction system TOPAZ. Presently, observations of sea ice concentration and sea ice drift are introduced in the assimilation step, but data for sea ice thickness and ice age (or roughness) are not included. The model computes the age of the ice by recording and updating the time passed after ice formation as sea ice grows and deteriorates as it is advected inside the model domain. Ice that is younger than 365 days is classified as first year ice. The fraction of first-year ice is recorded as a tracer in each grid cell. The Ocean and Sea Ice Thematic Assembly Centre in CMEMS redistributes a daily product from the EUMETSAT OSI SAF of gridded sea ice conditions which include "ice type", a representation of the separation of regions between those infested by first year ice, and those infested by multi-year ice. The ice type is parameterized based on data for the gradient ratio GR(19,37) from SSMIS observations, and from the ASCAT backscatter parameter. This product also includes information on ambiguity in the processing of the remote sensing data, and the product's confidence level, which have a strong seasonal dependency.

Modeling of subaqueous melting in Petermann Fjord, Northwestern Greenland using an ocean general circulation model

NASA Astrophysics Data System (ADS)

Cai, C.; Rignot, E. J.; Xu, Y.; An, L.; Tinto, K. J.; van den Broeke, M. R.

2014-12-01

Basal melting of the floating tongue of Petermann Glacier, in northwestern Greenland is by far the largest process of mass ablation. Melting of the floating tongue is controlled by the buoyancy of the melt water plume, the pressure-dependence of the melting point of sea ice, and the mixing of warm subsurface water with fresh buoyant subglacial discharge. In prior simulations of this melting process, the role of subglacial discharge has been neglected because in similar configurations (floating ice shelves) in the Antarctic, surface runoff is negligible; this is however not true in Greenland. Here, we use the Mass Institute of Technology general circulation model (MITgcm) at a high spatial resolution (10 m x 10 m) to simulate the melting process of the ice shelf in 2-D. The model is constrained by ice shelf bathymetry and ice thickness (refined model in the immediate vicinity of the grounding line) from NASA Operation IceBridge (2011), ocean temperature/salinity data from Johnson et al. (2011), ocean tide height and current from the Arctic Ocean Tidal Inverse Model (AOTIM-5) by Padman and Erofeeva (2004) and subglacial discharge at the grounding line calculated by the hydrostatic potential of the ice from estimated products of the Regional Atmospheric Climate Model (RACMO) of Royal Netherlands Meteorological Institute (KNMI). We compare the results obtained in winter (no runoff) with summer, and the sensitivity of the results to thermal forcing from the ocean, and to the variation of tide height and current, and to the magnitude of subglacial runoff. We conclude on the impact of the ocean and surface melting on the melting regime of the floating ice tongue of Petermann. The basal melt rate increases ~20% with summer surface runoff. This work is performed under a contract with NASA Cryosphere Program.

An Imminent Revolution in Modeling Interactions of Ice Sheets With Climate

NASA Astrophysics Data System (ADS)

Hughes, T.

2008-12-01

Modeling continental ice sheets was inaugurated by meteorologists William Budd and Uwe Radok, with mathematician Richard Jenssen, in 1971. Their model calculated the thermal and mechanical regime using measured surface accumulation rates, temperatures, and elevations, and bed topography. This top-down approach delivered a basal thermal regime of temperatures or melting rates for an assumed basal geothermal heat flux. When Philippe Huybrechts and others incorporated time, largely unknownpast surface conditions had a major effect on present basal thermal conditions. This approach produced ice-sheet models with only a slow response to external forcing, whereas the glacial geological record and climate records from ice and ocean cores show that ice sheets can have rapid changes in size and shape independent of external forcing. These top-down models were wholly inadequate for reconstructing former ice sheets at the LGM for CLIMAP in 1981. Ice-sheet areas,elevations, and volumes provided the albedo, surface topography, and sea-surface area as input to climate models. A bottom-up model based on dated glacial geology was developed to provide the areal extent and basal thermal regime of ice sheets at the LGM. Basal thermal conditions determined ice-bed coupling and therefore the elevation of ice sheets. High convex ice surfaces for slow sheet flow lower about 20 percent when a frozen bed becomes thawed. As further basal melting drowns bedrock bumps that "pin" basal ice, the ice surface becomes concave in fast stream flow that ends as low floating ice shelves at marine ice margins. A revolution in modeling interactions between glaciation, climate, and sea level is driven by new Greenland and Antarctic data from Earth-orbiting satellites, airborne and surface traverses, and deep drilling. We anticipate continuous data acquisition of surface albedo, accumulation/ablation rates, elevations, velocities, and temperatures over a whole ice sheet, mapping basal thermal conditions

Insights into ice-ocean interactions and fjord circulation from fjord sea surface temperatures at the Petermann Glacier, Greenland

NASA Astrophysics Data System (ADS)

Snow, T.; Shepherd, B.; Abdalati, W.; Scambos, T. A.

2016-12-01

Dynamic processes at marine-terminating outlet glaciers are responsible for over one-third of Greenland Ice Sheet (GIS) mass loss. Enhanced intrusion of warm ocean waters at the termini of these glaciers has contributed to elevated rates of ice thinning and terminus retreat over the last two decades. In situ oceanographic measurements and modeling studies show that basal melting of glaciers and subglacial discharge can cause buoyant plumes of water to rise to the fjord surface and influence fjord circulation characteristics. The temperature of these surface waters holds clues about ice-ocean interactions and small-scale circulation features along the glacier terminus that could contribute to outlet glacier mass loss, but the magnitude and duration of temperature variability remains uncertain. Satellite remote sensing has proven very effectiver for acquiring sea surface temperatuer (SST) data from these remote regions on a long-term, consistent basis and shows promise for identifying temperature anomalies at the ice front. However, these data sets have not been widely utilized to date. Here, we use satellite-derived sea surface temperatures to identify fjord surface outflow characteristics from 2000 to present at the Petermann Glacier, which drains 4% of the GIS and is experiencing 80% of its mass loss from basal melt. We find a general SST warming trend that coincides with early sea ice breakup and precedes two major calving events and ice speedup that began in 2010. Persistent SST anomalies along the terminus provide evidence of warm outflow that is consistent with buoyant plume model predictions. However, the anomalies are not evident early in the time series, suggesting that ocean inflow and ice-ocean interactions have experienced a regime shift since 2000. Our results provide valuable insight into fjord circulation patterns and the forcing mechanisms that contribute to terminus retreat. Comparing our results to ongoing modeling experiments, time series from

Providing Real-time Sea Ice Modeling Support to the U.S. Coast Guard

NASA Astrophysics Data System (ADS)

Allard, Richard; Dykes, James; Hebert, David; Posey, Pamela; Rogers, Erick; Wallcraft, Alan; Phelps, Michael; Smedstad, Ole Martin; Wang, Shouping; Geiszler, Dan

2016-04-01

The Naval Research Laboratory (NRL) supported the U.S. Coast Guard Research Development Center (RDC) through a demonstration project during the summer and autumn of 2015. Specifically, a modeling system composed of a mesoscale atmospheric model, regional sea ice model, and regional wave model were loosely coupled to provide real-time 72-hr forecasts of environmental conditions for the Beaufort/Chukchi Seas. The system components included a 2-km regional Community Ice CodE (CICE) sea ice model, 15-km Coupled Ocean Atmosphere Mesoscale Prediction System (COAMPS) atmospheric model, and a 5-km regional WAVEWATCH III wave model. The wave model utilized modeled sea ice concentration fields to incorporate the effects of sea ice on waves. The other modeling components assimilated atmosphere, ocean, and ice observations available from satellite and in situ sources. The modeling system generated daily 72-hr forecasts of synoptic weather (including visibility), ice drift, ice thickness, ice concentration and ice strength for missions within the economic exclusion zone off the coast of Alaska and a transit to the North Pole in support of the National Science Foundation GEOTRACES cruise. Model forecasts graphics were shared on a common web page with selected graphical products made available via ftp for bandwidth limited users. Model ice thickness and ice drift show very good agreement compared with Cold Regions Research and Engineering Laboratory (CRREL) Ice Mass Balance buoys. This demonstration served as a precursor to a fully coupled atmosphere-ocean-wave-ice modeling system under development. National Ice Center (NIC) analysts used these model data products (CICE and COAMPS) along with other existing model and satellite data to produce the predicted 48-hr position of the ice edge. The NIC served as a liaison with the RDC and NRL to provide feedback on the model predictions. This evaluation provides a baseline analysis of the current models for future comparison studies

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

NASA Technical Reports Server (NTRS)

Hakkinen, Sirpa

1993-01-01

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.

Geophysics of an Oceanic Ice Shell on Snowball Earth

NASA Technical Reports Server (NTRS)

Gaidos, E. J.

2000-01-01

Kirschvink proposed Precambrian low-latitude glaciation could result in an albedo-driven catastrophic runaway to a "Snowball Earth" state in which pack ice up to 1 km thick covered the world ocean. The geophysical state of an ice crust on a Snowball Earth is examined.

Winter ocean-ice interactions under thin sea ice observed by IAOOS platforms during N-ICE2015: Salty surface mixed layer and active basal melt

NASA Astrophysics Data System (ADS)

Koenig, Zoé; Provost, Christine; Villacieros-Robineau, Nicolas; Sennéchael, Nathalie; Meyer, Amelie

2016-10-01

IAOOS (Ice Atmosphere Arctic Ocean Observing System) platforms, measuring physical parameters at the atmosphere-snow-ice-ocean interface deployed as part of the N-ICE2015 campaign, provide new insights on winter conditions North of Svalbard. The three regions crossed during the drifts, the Nansen Basin, the Sofia Deep, and the Svalbard northern continental slope featured distinct hydrographic properties and ice-ocean exchanges. In the Nansen Basin, the quiescent warm layer was capped by a stepped halocline (60 and 110 m) and a deep thermocline (110 m). Ice was forming and the winter mixed layer salinity was larger by ˜0.1 g/kg than previously observed. Over the Svalbard continental slope, the Atlantic Water (AW) was very shallow (20 m from the surface) and extended offshore from the 500 m isobath by a distance of about 70 km, sank along the slope (40 m from the surface) and probably shed eddies into the Sofia Deep. In the Sofia Deep, relatively warm waters of Atlantic origin extended from 90 m downward. Resulting from different pathways, these waters had a wide range of hydrographic characteristics. Sea-ice melt was widespread over the Svalbard continental slope and ocean-to-ice heat fluxes reached values of 400 W m-2 (mean of ˜150 W m-2 over the continental slope). Sea-ice melt events were associated with near 12 h fluctuations in the mixed-layer temperature and salinity corresponding to the periodicity of tides and near-inertial waves potentially generated by winter storms, large barotropic tides over steep topography, and/or geostrophic adjustments.

Sea ice and oceanic processes on the Ross Sea continental shelf

NASA Technical Reports Server (NTRS)

Jacobs, S. S.; Comiso, J. C.

1989-01-01

The spatial and temporal variability of Antarctic sea ice concentrations on the Ross Sea continental shelf have been investigated in relation to oceanic and atmospheric forcing. Sea ice data were derived from Nimbus 7 scanning multichannel microwave radiometer (SMMR) brightness temperatures from 1979-1986. Ice cover over the shelf was persistently lower than above the adjacent deep ocean, averaging 86 percent during winter with little month-to-month of interannual variability. The large spring Ross Sea polynya on the western shelf results in a longer period of summer insolation, greater surface layer heat storage, and later ice formation in that region the following autumn.

Modeling of Antarctic Sea Ice in a General Circulation Model.

NASA Astrophysics Data System (ADS)

Wu, Xingren; Simmonds, Ian; Budd, W. F.

1997-04-01

A dynamic-thermodynamic sea ice model is developed and coupled with the Melbourne University general circulation model to simulate the seasonal cycle of the Antarctic sea ice distribution. The model is efficient, rapid to compute, and useful for a range of climate studies. The thermodynamic part of the sea ice model is similar to that developed by Parkinson and Washington, the dynamics contain a simplified ice rheology that resists compression. The thermodynamics is based on energy conservation at the top surface of the ice/snow, the ice/water interface, and the open water area to determine the ice formation, accretion, and ablation. A lead parameterization is introduced with an effective partitioning scheme for freezing between and under the ice floes. The dynamic calculation determines the motion of ice, which is forced with the atmospheric wind, taking account of ice resistance and rafting. The simulated sea ice distribution compares reasonably well with observations. The seasonal cycle of ice extent is well simulated in phase as well as in magnitude. Simulated sea ice thickness and concentration are also in good agreement with observations over most regions and serve to indicate the importance of advection and ocean drift in the determination of the sea ice distribution.

Modeling of Antarctic sea ice in a general circulation model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wu, Xingren; Budd, W.F.; Simmonds, I.

1997-04-01

A dynamic-thermodynamic sea ice model is developed and coupled with the Melbourne University general circulation model to simulate the seasonal cycle of the Antarctic sea ice distributions The model is efficient, rapid to compute, and useful for a range of climate studies. The thermodynamic part of the sea ice model is similar to that developed by Parkinson and Washington, the dynamics contain a simplified ice rheology that resists compression. The thermodynamics is based on energy conservation at the top surface of the ice/snow, the ice/water interface, and the open water area to determine the ice formation, accretion, and ablation. Amore » lead parameterization is introduced with an effective partitioning scheme for freezing between and under the ice floes. The dynamic calculation determines the motion of ice, which is forced with the atmospheric wind, taking account of ice resistance and rafting. The simulated sea ice distribution compares reasonably well with observations. The seasonal cycle of ice extent is well simulated in phase as well as in magnitude. Simulated sea ice thickness and concentration are also in good agreement with observations over most regions and serve to indicate the importance of advection and ocean drift in the determination of the sea ice distribution. 64 refs., 15 figs., 2 tabs.« less

Subaqueous melting in Zachariae Isstrom, Northeast Greenland combining observations and an ocean general circulation model

NASA Astrophysics Data System (ADS)

Cai, C.; Rignot, E. J.; Menemenlis, D.; Nakayama, Y.

2016-12-01

Zachariae Isstrom, a major ice stream in northeast Greenland, has lost its entire ice shelf in the past decade. Here, we study the evolution of subaqueous melting of its floating section during the transition. Observations show that the rate of ice shelf melting has doubled during 1999-2010 and is twice higher than that maintaining the ice shelf in a steady state. The ice shelf melt rate depends on the thermal forcing from warm, saline, subsurface ocean water of Atlantic origin (AW), and on the mixing of AW with fresh buoyant subglacial discharge. Subglacial discharge has increased as result of enhanced ice sheet runoff driven by warmer air temperature; ocean thermal forcing has increased due to enhanced advection of AW. Here, we employ the Massachusetts Institute of Technology general circulation model (MITgcm) at a high spatial resolution to simulate the melting process in 3-D. The model is constrained by ice thickness from mass conservation, oceanic bathymetry inverted from gravity data by NASA Operation IceBridge and NASA Ocean Melting Greenland missions, in-situ ocean temperature/salinity data, ocean tide height and current from the Arctic Ocean Tidal Inverse Model (AOTIM-5) and reconstructed seasonal subglacial discharge from the Regional Atmospheric Climate Model (RACMO2). We compare the results in winter (small runoff but not negligible) with summer (maximum runoff) at two different stages with (prior to 2012) and without the ice shelf (after 2012) to subaqueous melt rates deduced from remote sensing observations. We show that ice melting by the ocean has increased by one order of magnitude as a result of the transition from ice shelf terminating to near-vertical calving front terminating. We also find that subglacial discharge has a significant impact on ice shelf melt rates in Greenland. We conclude on the impact of ocean warming and air temperature warming on the melting regime of the ice margin of Zachariae Isstrom, Greenland. This work was performed