Validation of a global finite element sea ice-ocean model

Results from a global Finite Element Sea iceOcean Model (FESOM) are evaluated using a wide range of observational datasets. FESOMs ocean component is a primitive-equation, hydrostatic ocean model using isopycnic diffusion and a Gent-McWilliams scheme to parameterize the effects of sub-gridscale turbulence on tracer distribution. Vertical mixing and convection are parameterized as a function of the Richardson number and the Monin-Obukhov length. A finite element dynamic-thermodynamic sea icemodel with elastic-viscous-plastic rheology has been developed and coupled to the ocean component. The model features a prognostic snow layer but neglects internal heat storage. All model components are discretized on a triangular/tetrahedral grid with a continuous, conformingrepresentation of model variables. The coupled model has been run in a global configuration and forced with NCEP daily atmospheric reanalysis data for 1948-2007. Results are analysed with a focus on the Southern Hemisphere. While summer ice extent is underestimated in both hemispheres, winter ice extents are in good agreement with satellite data. Southern Ocean sea ice thickness distribution agrees well with ship-based observations and even quantitatively with data from upwards looking sonars (ULS). Sea ice freezing rates have been validated using repeated salinity profiles from Southern Elephant Seals. Gulf Stream transport is underestimated, but transports of the Kuroshio and the Antarctic Circumpolar Current appear realistic. A comparison of numerical tracer studies to observed CFC distribution indicates that bottom layer ventilation occurs on realistic pathways. Global meridional overturning features a strong Antarctic Bottom Water (AABW) cell, while the formation of North Atlantic Deep Water (NADW) appears to be on the weak side. Besides pure model validation, the study also identifies regions and processes that critically require a locally increased horizontal resolution in order to be represented adequately

Southern Ocean warming: Increase in basal melting and grounded ice loss

We apply a global finite element sea ice/ice shelf/ocean model (FESOM) to the Antarctic marginal seas to analyze projections of ice shelf basal melting in a warmer climate. The model is forced with the atmospheric output from two climate models: (1) the Hadley Centre Climate Model (HadCM3) and (2) Max Planck Institute’s ECHAM5/MPI-OM. Results from their 20th-century simulations are used to evaluate the modeled present-day ocean state. Sea-ice coverage is largely realistic in both simulations. Modeled ice shelf basal melt rates compare well with observations in both cases, but are consistently smaller for ECHAM5/MPI-OM. Projections for future ice shelf basal melting are computed using atmospheric output for IPCC scenarios E1 and A1B. While trends in sea ice coverage, ocean heat content, and ice shelf basal melting are small in simulations forced with ECHAM5 data, a substantial shift towards a warmer regime is found in experiments forced with HadCM3 output. A strong sensitivity of basal melting to increased ocean temperatures is found for the ice shelves in the Amundsen Sea. For the cold-water ice shelves in the Ross and Weddell Seas,decreasing convection on the continental shelf in the HadCM3 scenarios leads to an erosion of the continental slope front and to warm water of open ocean origin entering the continental shelf. As this water reaches deep into the Filchner-Ronne Ice Shelf (FRIS) cavity, basal melting increases by a factor of three to six compared to the present value of about 100 Gt/yr. Highest melt rates at the deep FRIS grounding line causes a retreat of > 200km, equivalent to an land ice loss of 110 Gt/yr

Feedbacks between ice and ocean dynamics at the West Antarctic Filchner-Ronne Ice Shelf in future global warming scenarios

The ice flow at the margins of the West Antarctic Ice Sheet is moderated by large ice shelves. Their buttressing effect substantially controls the mass balance of the WAIS and thus its contribution to sea level rise. The stability of these ice shelves results from the balance of mass gain by accumulation and ice flow from the adjacent ice sheet and mass loss by calving and basal melting due to the ocean heat flux. Recent results of ocean circulation models indicate that warm circumpolar water of the Southern Ocean may override the submarine slope front of the Antarctic Continent and boost basal ice shelf melting. In particular, ocean simulations for several of the IPCC's future climate scenarios demonstrate the redirection of a warm coastal current into the Filchner Trough and underneath the Filchner-Ronne Ice Shelf within the next decades. In this study, we couple the finite elements ocean circulation model FESOM and the three-dimensional thermomechanical ice flow model RIMBAY to investigate the complex interactions between ocean and ice dynamics at the Filchner-Ronne Ice Shelf. We focus on the impact of a changing ice shelf cavity on ocean dynamics as well as the feedback of the resulting sub-shelf melting rates on the ice shelf geometry and implications for the dynamics of the adjacent marine-based Westantarctic Ice Sheet. Our simulations reveal the high sensitivity of grounding line migration to ice-ocean interactions within the Filchner-Ronne Ice Shelf and emphasize the importance of coupled model studies for realistic assessments of the Antarctic mass balance in future global warming scenarios

Impact of coastal polynyas on dense shelf water formation in the Weddell Sea

Dense shelf water is an essential ingredient to the formation of Antarctic Bottom Water (AABW). It is formed on the continental shelves surrounding Antarctica, when freezing rates are sufficiently high to push ocean salinity to values of 34.65 and higher. Coastal polynyas, where the ice is driven away from the coastline, maintain the highest freezing rates in Antarctic winter. Since theWeddell Sea is considered the most productive source region of AABW, we investigate the dense water formation on the continental shelves of the southwestern Weddell Sea, with a focus on the role of coastal polynyas, using the Finite Element Sea ice-Ocean Model (FESOM), a primitive-equation, hydrostatic ocean model coupled with a dynamic-thermodynamic sea ice model. The horizontal resolution of the global, unstructured mesh is up to 3 km at the southwestern Weddell Sea coastline; in vertical direction the mesh features 37 depth levels (resolution increases toward the surface). The model was initialized on 01/01/1980 with data from the Polar Hydrographic Climatology and forced with NCEP/NCAR Reanalysis data. The 20-year period 1990-2009 is used for analysis. Our results indicate that in an average year, the polynya freezing rates of 9 cm d--1 (corresponding to a salt input of 2.5 kg m--2d--1) cause a seasonal variation in salinity of 0.3 psu under the Ronne polynya and result in the production of 5.10-4 km-3 dense shelf water, which leaves the continental shelf (outlined by the 700 m isobath in this study) at a long-term mean volume flux of 5.2 Sv. Some of this water contributes to the formation of Weddell Sea Deep/BottomWater, but a large fraction is diluted by mixing with ambient water and leaves the Weddell Sea at intermediate levels

The role of ice streams in a coupled ice flow-ocean modeling approach at the Filchner-Ronne Ice Shelf, Antarctica

The ice flow at the margins of the Antarctic Ice Sheet (AIS) is moderated by large ice shelves. Their buttressing effect substantially controls the mass balance of the AIS and thus its contribution to sea level rise. Recent results of ocean circulation models indicate that warm circumpolar water of the Southern Ocean may override the continental slope front and boost basal ice shelf melting. In particular, simulations demonstrate the redirection of a warm coastal current into the Filchner Trough and underneath the Filchner-Ronne Ice Shelf (FRIS) within the next decades. The increase of water temperature in the sub-shelf cavity is estimated to dramatically raise the basal shelf melting. Coupled simulations with a finite elements ocean model and a three-dimensional thermomechanical ice flow model reveal that the consequent thinning of the FRIS would lead to an extensive grounding line retreat associated with a vast mass loss of the AIS. In this subsequent study, we aim for an enhanced understanding of the complex feedbacks between ocean circulation and ice dynamics of the grounded AIS. Therefor, we focus on the ice streams which are draining into the FRIS and dominating the mass transport from grounded to floating ice. For a better representation of these fast-flowing ice streams we expand the above ice flow model by the incorporation of local processes at the ice base. There, sediment deformation and lubrication by subglacial hydrology locally allow high basal sliding rates and thus create the precondition for the development of ice streams. Based on satellite-observed ice surface velocity patterns we identify such areas with low basal drag and parametrize the ice flow model accordingly. As a result, the modeled ice flow patterns will depict velocity and locations of observed ice streams in the catchment of the FRIS more realistically. We present first results of this advanced ice-flow modeling approach, anticipating an even larger response of the AIS on increased sub-shelf melting rates in future coupled simulations

The Effect of Flow at Maud Rise on the Sea Ice Cover - Numerical Experiments

The role of seamounts in the formation and evolution of sea ice isinvestigated in a series of numerical experiments with a coupled seaice-ocean model. Bottom topography, stratification and forcing areconfigured for the Maud Rise region in the Weddell Sea. The specificflow regime that develops at the seamount as the combined response tosteady and tidal forcing consists of free and trapped waves and aTaylor column, which is caused by mean flow and tidal flowrectification. The enhanced variability through tidal motion inparticular is capable of modifying the mixed layer above the seamountenough to delay and reduce sea ice formation throughout the winter.The induced sea ice anomaly spreads and moves westward and affects anarea of several 100~000 km$^{2}$. Process studies reveal the complexinteraction between wind, steady and periodic ocean currents: allthree are required in the process of generation of the sea ice andmixed layer anomalies (mainly through tidal flow), their detachmentfrom the topography (caused by steady oceanic flow), and the westwardtranslation of the sea ice anomaly (driven by the time-mean wind)

Simulation of coastal polynyas in the western Weddell Sea

Coastal polynyas play a prominent role in the formation and modification of water masses in the polar oceans. A coastal polynya is usually kept open mechanically, primarily by winds, and the ocean surface is at freezing point. Thus a major fraction of the annual ice production of the high-latitude oceans occurs in polynyas and hence the duration and extent of their appearance has a substantial effect on bottom water formation. In the western Weddell Sea, recurring coastal polynyas are formed in front of the Filchner-Ronne Ice Shelf and in the area of the decayed Larsen A/B Ice Shelf. Simulations to study polynya formation and their impact on ice production and bottom water formation in the western Weddell Sea were performed with the Finite Element Sea ice-Ocean Model (FESOM) of Alfred-Wegener-Institute (AWI). FESOM is a fully coupled system of a primitive-equation, hydrostatic ocean model and a dynamic-thermodynamic sea ice model. The simulations were conducted on a global grid with a resolution varying between roughly 300 km in tropical latitudes and <5 km along the coast of the southwestern Weddell Sea. In vertical direction, the grid uses terrain-following coordinates. The model results give insight into the mechanisms governing the formation of transient and persistent polynyas and their influence on ice production and deep water formation. Water mass formation and ice export rates are quantified and compared to observation-based estimates