The ice sheet-ocean modeling community is making large strides toward developing coupled models capable of examining the interactions and feedbacks between ice shelves and ocean along the Antarctic margin. We present preliminary results and address some of the challenges that have arisen during the development of a coupled ice sheet-ocean model. The ice sheet model is icepack, a shallow-shelf finite element model written in Python. The ocean model is the Regional Ocean Modelling System (ROMS), a terrain-following vertical (sigma) coordinate model that has been modified to interface with a moving ice shelf. These two models are coupled in an online configuration using the Framework for Ice Sheet Ocean Coupling (FISOC). The use of a model with sigma coordinates for the ocean component introduces a simplification and a complication to modeling a moving ice draft. The sigma coordinate system retains the same number of vertical layers at any depth, eliminating the need to convert grid cells between ice and water, when using a fixed grounding line configuration. However, as the ice shelf draft evolves in time, topographic configurations develop that induce pressure gradient errors in ROMS. We quantify these errors in an idealized set-up with an artificially changing ice draft following the ISOMIP+ geometry. We compare results between an ice draft that is smoothed to meet standard ROMS smoothing criteria (rx0, rx1) and a non-smoothed ice draft. Finally, we present a simple parameterization in a buffer zone near the grounding line that uses interpolated melt rates from the ocean model, allowing us to maintain a steep ice topography in the ice model without inducing pressure gradient errors in the short water column in the ocean model. This model configuration will be applied to Pine Island Glacier and used to examine present and possible future states of the ice sheet-ocean system.

Pine Island Glacier Ice Shelf (PIGIS) is melting rapidly from beneath due to the circulation of relatively warm water under the ice shelf, driven primarily by buoyancy of the meltwater plume. Basal melt rates predicted by ocean models with thermodynamically active ice shelves depend on the representation of environmental characteristics including geometry (grounding line location, ice draft and seabed bathymetry) and ocean hydrographic conditions, and subgrid-scale parameterizations. We developed a relatively high resolution (lateral grid spacing of 0.5 km, 24 terrain following levels) model for the PIGIS vicinity based on the Regional Ocean Modeling System (ROMS). Initial stratification was specified with idealized profiles based on observed hydrographic data seaward of the ice front. Predicted basal melt rate distributions were compared with satellite-derived estimates and stratification beneath PIGIS was compared with Autosub profiles. As in previous studies, we found that the melt rate was strongly dependent on the (specified) depth of the thermocline separating cold surface waters from deep, relatively warm waters, and on the presence of a submarine ridge under the ice shelf that impedes circulation of warm deep water into the back portion of the cavity. Melt rates were sensitive to the model’s subgrid-scale parameterizations. The quadratic drag coefficient, which parameterizes roughness of the ice shelf base, had a substantial effect on the melt rate through its role in the three-equation formulation for ice-ocean buoyancy exchange. Turbulent tracer diffusion, which was parameterized by a constant value or various mixed layer models, played an important role in determining stratification in the cavity. Numerical diffusion became significant in some cases. We conclude that flow of warm water into the inner portion of the PIGIS cavity near the deep grounding line is sensitive to poorly constrained mixing parameterizations, both at the ice base and as a mechanism for allowing inflowing ocean heat to cross the sub-ice-shelf sill. Improved understanding of mixing processes is required as the community moves towards fully coupled ocean/ice-sheet models with evolving ice thickness and grounding lines.

Multi-decadal expansion in the winter maximum sea ice extent (SIE) around Antarctica was interrupted by contraction, beginning in 2016 and continuing into 2019. This unexpected behavior motivates a closer look at factors controlling the position of the outer ice margin.We analyzed sea ice concentration (SIC) estimates derived from passive microwave sensors with differing resolutions (SSM/I, AMSR-E and AMSR2) to identify spatial and temporal statistics of the sea ice edge deVned by 15% SIC. The low-pass Vltered position of the ice edge is similar in different products, with the maximum northward position determined by proximity to the relatively warm waters of the Antarctic Circumpolar Current. Higher resolution SIC products reveal greater spatial detail along the convoluted margin, resulting in a relatively longer sea ice perimeter. Spectral analysis does not identify statistically signiVcant peaks in length scales along the margin; however, visual comparison with geostrophic velocities and sea surface temperature inferred from satellite altimetry suggests that advection of sea ice by mesoscale eddies is an important mechanism for deforming the ice edge in some regions, such as the Bellingshausen Sea. We analyze a high-resolution (dx=5 km), coupled ocean-sea ice model which realistically represents the annual expansion of sea ice to quantify the dynamic and thermodynamic roles of eddies in sea-ice mass balance and SIE. These eddy effects on the sea ice edge are not well represented in coarser-grid ocean reanalysis products such as ECCO-2, motivating an investigation of how to represent eddy/sea-ice interactions in global climate models.