Finite Element Sea-ice Ocean Model Research Articles

Simulated Signatures of Greenland Melting in the North Atlantic: A Model Comparison With Argo Floats, Satellite Observations, and Ocean Reanalysis

AbstractIncreased Greenland ice sheet melting has an impact on global mean and regional sea level rise and the ocean circulation. In this study, we explore whether Greenland melting signatures found in ocean model simulations are visible in observations from radar altimetry, satellite gravimetry and Argo floats. We have included Greenland freshwater flux (GF) in the global Finite‐Element‐Sea ice‐Ocean Model (FESOM) for the years 1993–2016. The reference run is computed by excluding Greenland freshwater input. These experiments are performed on a low resolution (ca. 24 km) and a high resolution (ca. 6 km) eddy‐permitting mesh. For comparison with the model experiments, we use different observational data, such as Argo floats, satellite observations, and reanalyses. We find that surface GF maps into signatures in temperature and salinity down to about 100 m in the surroundings of Greenland. The simulated melting signatures are particularly visible in steric heights in Baffin Bay and Davis Strait. Here, we find an improvement of the mean square error of up to 30% when including GF. For the Nordic part of the Nordic Seas, however, we find no improvement when including GF. We compare steric heights with reanalysis data and a new setup of the inversion method from gravimetric and altimetric satellite data. We cannot confirm that the GF signatures on variables such as temperature and salinity are visible in the observations on the time scales considered. However, we find that increased model resolution often causes larger improvements than occur due to including the simulated melting effect.

Impacts of strong wind events on sea ice and water mass properties in Antarctic coastal polynyas

Strong offshore wind events (SOWEs) occur frequently near the Antarctic coast during austral winter. These wind events are typically associated with passage of synoptic- or meso-scale cyclones, which interact with the katabatic wind field and affect sea ice and oceanic processes in coastal polynyas. Based on numerical simulations from the coupled Finite Element Sea-ice Ocean Model (FESOM) driven by the CORE-II forcing, two coastal polynyas along the East Antarctica coast––the Prydz Bay Polynya and the Shackleton Polynya are selected to examine the response of sea ice and oceanic properties to SOWEs. In these polynyas, the southern or western flanks of cyclones play a crucial role in increasing the offshore winds depending on the local topography. Case studies for both polynyas show that during SOWEs, when the wind speed is 2–3 times higher than normal values, the offshore component of sea ice velocity can increase by 3–4 times. Sea ice concentration can decrease by 20–40%, and sea ice production can increase up to two to four folds. SOWEs increase surface salinity variability and mixed layer depth, and such effects may persist for 5–10 days. Formation of high salinity shelf water (HSSW) is detected in the coastal regions from surface to 800 m after 10–15 days of the SOWEs, while the HSSW features in deep layers exhibit weak response on the synoptic time scale. HSSW formation averaged over winter is notably greater in years with longer duration of SOWEs.

Properties and dynamics of mesoscale eddies in Fram Strait from a comparison between two high-resolution ocean–sea ice models

Abstract. Fram Strait, the deepest gateway to the Arctic Ocean, is strongly influenced by eddy dynamics. Here we analyse the output from two eddy-resolving models (ROMS – Regional Ocean Modeling System; FESOM – Finite-Element Sea-ice Ocean Model) with around 1 km mesh resolution in Fram Strait, with a focus on their representation of eddy properties and dynamics. A comparison with mooring observations shows that both models reasonably simulate hydrography and eddy kinetic energy. Despite differences in model formulation, they show relatively similar eddy properties. The eddies have a mean radius of 4.9 and 5.6 km in ROMS and FESOM, respectively, with slightly more cyclones (ROMS: 54 %, FESOM: 55 %) than anticyclones. The mean lifetime of detected eddies is relatively short in both simulations (ROMS: 10 d, FESOM: 11 d), and the mean travel distance is 35 km in both models. More anticyclones are trapped in deep depressions or move toward deep locations. The two models show comparable spatial patterns of baroclinic and barotropic instability. ROMS has relatively stronger eddy intensity and baroclinic instability, possibly due to its smaller grid size, while FESOM has stronger eddy kinetic energy in the West Spitsbergen Current. Overall, the relatively good agreement between the two models strengthens our confidence in their ability to realistically represent the Fram Strait ocean dynamics and also highlights the need for very high mesh resolution.

Eddies in the North Greenland Sea and Fram Strait From Satellite Altimetry, SAR and High‐Resolution Model Data

AbstractIn this study, we investigate eddy dynamics in the northern Greenland Sea and the Fram Strait using AVISO altimetry, spaceborne synthetic aperture radar (SAR), and Finite Element Sea ice‐Ocean Model (FESOM) high‐resolution numerical model data. In the region, eddies are thought to play an important role in the redistribution of the warmer and saltier Atlantic Water between the Arctic Ocean and the areas of deep convection in the central Greenland Sea. We found that eddies detected in AVISO and in SAR form two complementary data sets of large mesoscale eddies (with typical radii of 30–50 km) and of small mesoscale/submesoscale eddies (with typical radii of 1–5 km and not exceeding 30 km), respectively. For large mesoscale eddies, the number of cyclones and anticyclones is approximately the same, while for submesoscale eddies, cyclones are strongly dominating. The limitations and possible biases in each of the data sets are discussed and cross‐validated against FESOM results. It is noted that the most energetic eddies are concentrated along the major currents and in the northern part of the region. Eddy translations follow the mean currents in their overall cyclonic circulation around the northern Greenland Sea. A convergence of the eddies toward the Nordbukta area is detected. On seasonal time scale, a higher number of more intense mesoscale eddies is observed during winter, associated with a quasi‐simultaneous intensification of the mean currents. Model results also show an increase in the number of small eddies in spring‐early summer attributed to the decay of large eddies, while in late autumn, the opposite tendency suggests eddy merger.

Low-frequency sea level variability and impact of recent sea ice decline on the sea level trend in the Arctic Ocean from a high-resolution simulation

The Arctic Ocean is undergoing significant changes, with rapid sea ice decline, unprecedented freshwater accumulation, and pronounced regional sea level rise. In this paper, we analyzed the sea level variation in the Arctic Ocean based on a global simulation with 4.5-km resolution in the Arctic Ocean using the multi-resolution Finite Element Sea Ice-Ocean Model (FESOM). The simulation reasonably reproduces both the main spatial features of the sea surface height (SSH) and its temporal evolution in the Arctic Ocean in comparison with tide gauge and satellite data. Using the model results, we investigated the low-frequency variability of the Arctic SSH. Both the first two dominant modes of the annual-mean SSH evolution in the Arctic Ocean present decadal variability and can be mainly attributed to the variability of the halosteric height, thus the freshwater content. The first mode can be explained by the Arctic Oscillation (AO). The AO-related atmospheric circulation drives the accumulation and release of freshwater in the Arctic deep basin and the consequent ocean mass change over the continental shelf, leading to the antiphase changes in SSH between the shelf seas and the deep basin. The second mode shows an antiphase oscillation between the two Arctic deep basins, the Amerasian and Eurasian Basins, which is driven by the Arctic dipole anomaly (DA). The DA-related wind anomaly causes a spatial redistribution of freshwater between the two basins, leading to the antiphase SSH changes. By using a dedicated sensitivity simulation in which the recent sea ice decline is eliminated, we find that the sea ice decline contributed considerably to the observed sea level rise in the Amerasian Basin in the recent decades. Although the sea ice decline did not change the mean SSH averaged over the Arctic Ocean, it significantly changed the spatial pattern of the SSH trend. Our finding indicates that both the wind regime and ongoing sea ice decline should be considered to better understand and predict the changes in regional sea level in the Arctic Ocean.

Geostrophic currents in the northern Nordic Seas from a combination of multi-mission satellite altimetry and ocean modeling

Abstract. A deeper knowledge about geostrophic ocean surface currents in the northern Nordic Seas supports the understanding of ocean dynamics in an area affected by sea ice and rapidly changing environmental conditions. Monitoring these areas by satellite altimetry results in a fragmented and irregularly distributed data sampling and prevents the computation of homogeneous and highly resolved spatio-temporal datasets. In order to overcome this problem, an ocean model is used to fill in data when altimetry observations are missing. The present study provides a novel dataset based on a combination of along-track satellite-altimetry-derived dynamic ocean topography (DOT) elevations and simulated differential water heights (DWHs) from the Finite Element Sea ice Ocean Model (FESOM) version 1.4. This innovative dataset differs from classical assimilation methods because it substitutes altimetry data with the model output when altimetry fails or is not available. The combination approach is mainly based on a principal component analysis (PCA) after reducing both quantities by their constant and seasonal signals. In the main step, the most-dominant spatial patterns of the modeled differential water heights as provided by the PCA are linked with the temporal variability in the estimated DOT from altimetry by performing a principal component synthesis (PCS). After the combination, the annual signal obtained by altimetry and a constant offset are re-added in order to reference the final data product to the altimetry height level. Surface currents are computed by applying the geostrophic flow equations to the combined topography. The resulting final product is characterized by the spatial resolution of the ocean model around 1 km and the temporal variability in the altimetry along-track derived DOT heights. The combined DOT is compared to an independent DOT product, resulting in a positive correlation of about 80 %, to provide more detailed information about short periodic and finer spatial structures. The derived geostrophic velocity components are evaluated by in situ surface drifter observations. Summarizing all drifter observations in equally sized bins and comparing the velocity components shows good agreement in spatial patterns, magnitude and flow direction. Mean differences of 0.004 m s−1 in the zonal and 0.02 m s−1 in the meridional component are observed. A direct pointwise comparison between the combined geostrophic velocity components interpolated onto the drifter locations indicates that about 94 % of all residuals are smaller than 0.15 m s−1. The dataset is able to provide surface circulation information within the sea ice area and can be used to support a deeper comprehension of ocean currents in the northern Nordic Seas affected by rapid environmental changes in the 1995–2012 time period. The data are available at https://doi.org/10.1594/PANGAEA.900691 (Müller et al., 2019).

Effects of High Resolution and Spinup Time on Modeled North Atlantic Circulation

AbstractThe influence of a high horizontal resolution (5–15 km) on the general circulation and hydrography in the North Atlantic is investigated using the Finite Element Sea Ice–Ocean Model (FESOM). We find a stronger shift of the upper-ocean circulation and water mass properties during the model spinup in the high-resolution model version compared to the low-resolution (~1°) control run. In quasi equilibrium, the high-resolution model is able to reduce typical low-resolution model biases. Especially, it exhibits a weaker salinification of the North Atlantic subpolar gyre and a reduced mixed layer depth in the Labrador Sea. However, during the spinup adjustment, we see that initially improved high-resolution features partially reduce over time: the strength of the Atlantic overturning and the path of the North Atlantic Current are not maintained, and hence hydrographic biases known from low-resolution ocean models return in the high-resolution quasi-equilibrium state. We identify long baroclinic Rossby waves as a potential cause for the strong upper-ocean adjustment of the high-resolution model and conclude that a high horizontal resolution improves the state of the modeled ocean but the model integration length should be chosen carefully.

Dynamic ocean topography of the northern Nordic seas: a comparison between satellite altimetry and ocean modeling

Abstract. The dynamic ocean topography (DOT) of the polar seas can be described by satellite altimetry sea surface height observations combined with geoid information as well as by ocean models. The altimetry observations are characterized by an irregular sampling and seasonal sea ice coverage complicating reliable DOT estimations. Models display various spatiotemporal resolutions but are limited to their computational and mathematical context and introduced forcing models. In the present paper, ALES+ retracked altimetry ranges and derived along-track DOT heights of ESA's Envisat and water heights of the Finite Element Sea Ice-Ocean Model (FESOM) are compared to investigate similarities and discrepancies. The goal of the present paper is to identify to what extent pattern and variability of the northern Nordic seas derived from measurements and model agree with each other, respectively. The study period covers the years 2003–2009. An assessment analysis regarding seasonal DOT variabilities shows good agreement and confirms the dominant impact of the annual signal in both datasets. A comparison based on estimated regional annual signal components shows 2–3 times stronger amplitudes of the observations but good agreement of the phase. Reducing both datasets by constant offsets and the annual signal reveals small regional residuals and highly correlated DOT time series (Pearson linear correlation coefficient of at least 0.67). The highest correlations can be found in areas that are ice-free and affected by ocean currents. However, differences are visible in sea-ice-covered shelf regions. Furthermore, remaining constant artificial elevations in the observational data can be attributed to an insufficient representation of the used geoid. In general, the comparison results in good agreement between simulated and altimetry-based descriptions of the DOT in the northern Nordic seas.

Improving the Upper‐Ocean Temperature in an Ocean Climate Model (FESOM 1.4): Shortwave Penetration Versus Mixing Induced by Nonbreaking Surface Waves

AbstractAs the first mature global ocean general circulation model based on unstructured‐mesh methods, the multiresolution Finite Element Sea ice‐Ocean Model (FESOM) has shown great capability in reconstructing the ocean and sea ice in both standalone and coupled simulations at a relatively low computational cost. Parameterizations of some important processes, including the vertical mixing induced by surface waves, however, are still missing, contributing to temperature biases in the upper ocean. In this work we incorporate the vertical mixing induced by nonbreaking surface waves derived from a wave model into FESOM and compare its effect with that of shortwave penetration, another key process to vertically redistribute the heat in the upper ocean. Numerical experiments with and without the shortwave penetration scheme and the nonbreaking surface wave mixing reveal that both processes ameliorate the simulation of upper‐ocean temperature in middle and low latitudes mainly on the summer hemisphere. The role of nonbreaking surface waves is more pronounced in decreasing the mean cold biases at 50 m (by 1.0 °C, in comparison to 0.5 °C achieved by applying shortwave penetration). We conclude that the incorporation of mixing induced by nonbreaking surface waves into FESOM is practically very helpful and suggest that it needs to be considered in other ocean climate models as well.

Recent Sea Ice Decline Did Not Significantly Increase the Total Liquid Freshwater Content of the Arctic Ocean

Abstract The freshwater stored in the Arctic Ocean is an important component of the global climate system. Currently the Arctic liquid freshwater content (FWC) has reached a record high since the beginning of the last century. In this study we use numerical simulations to investigate the impact of sea ice decline on the Arctic liquid FWC and its spatial distribution. The global unstructured-mesh ocean general circulation model Finite Element Sea Ice–Ocean Model (FESOM) with 4.5-km horizontal resolution in the Arctic region is applied. The simulations show that sea ice decline increases the FWC by freshening the ocean through sea ice meltwater and modifies upper ocean circulation at the same time. The two effects together significantly increase the freshwater stored in the Amerasian basin and reduce its amount in the Eurasian basin. The salinification of the upper Eurasian basin is mainly caused by the reduction in the proportion of Pacific Water and the increase in that of Atlantic Water (AW). Consequently, the sea ice decline did not significantly contribute to the observed rapid increase in the Arctic total liquid FWC. However, the changes in the Arctic freshwater spatial distribution indicate that the influence of sea ice decline on the ocean environment is remarkable. Sea ice decline increases the amount of Barents Sea branch AW in the upper Arctic Ocean, thus reducing its supply to the deeper Arctic layers. This study suggests that all the dynamical processes sensitive to sea ice decline should be taken into account when understanding and predicting Arctic changes.

Arctic Ocean biogeochemistry in the high resolution FESOM 1.4-REcoM2 model

In the rapidly changing Arctic Ocean, marine primary productivity is tightly linked to the balance between light and nutrient limitation. To capture this balance in ocean general circulation biogeochemical models (OGCBMs), a good representation of the physics is important due to the tight bio-physical coupling in the Arctic. Using a horizontal model resolution of a few kilometers makes it possible to resolve an increasing number of small scale processes, that otherwise need to be parameterized in OGCBMs. Such high resolution is, however, commonly not possible due to computational constrains. Utilizing an unstructured mesh approach, we have run the finite element sea-ice ocean model (FESOM 1.4) coupled to the biogeochemical model REcoM2 in a global configuration with an Arctic-wide resolution of 4.5 km. This resolution is so far unprecedented for a global biogeochemical setup, and here we present an analysis of the mean state of the model. FESOM-REcoM2’s integrated Arctic net primary production (NPP) averages 445 Tg C yr−1 for the years 2011 to 2015, a value that is in the middle of the range compared to estimates from the literature. Most production takes place in the inflow regions of the Nordic and Chukchi Seas, and 32% is associated with the sea ice zone, the latter including the marginal ice zone and below-ice productivity. Light limits production to some degree at all latitudes north of 60°N, with growth becoming nutrient limited following the initial spring bloom in most places. The model reproduces the relatively low surface concentration of dissolved inorganic nitrogen (DIN) that has been observed in the central Arctic Ocean, as well as the low surface DIN concentration towards the end of the growth season further to the south, thereby also capturing widespread subsurface chlorophyll maxima (SCM). The SCMs are found in the whole Arctic except for the areas where sea ice concentration is high the whole year. They have a duration of two weeks to five months. The balance between nutrient and light limitations, both in the vertical and horizontal direction, highlights that decreased light limitation in a future ice free Arctic Ocean will not necessarily induce an increase in NPP due to increasing nutrient limitation, and that further studies of the role of the SCMs are required.

Intercomparison of Antarctic ice-shelf, ocean, and sea-ice interactions simulated by MetROMS-iceshelf and FESOM 1.4

Abstract. An increasing number of Southern Ocean models now include Antarctic ice-shelf cavities, and simulate thermodynamics at the ice-shelf/ocean interface. This adds another level of complexity to Southern Ocean simulations, as ice shelves interact directly with the ocean and indirectly with sea ice. Here, we present the first model intercomparison and evaluation of present-day ocean/sea-ice/ice-shelf interactions, as simulated by two models: a circumpolar Antarctic configuration of MetROMS (ROMS: Regional Ocean Modelling System coupled to CICE: Community Ice CodE) and the global model FESOM (Finite Element Sea-ice Ocean Model), where the latter is run at two different levels of horizontal resolution. From a circumpolar Antarctic perspective, we compare and evaluate simulated ice-shelf basal melting and sub-ice-shelf circulation, as well as sea-ice properties and Southern Ocean water mass characteristics as they influence the sub-ice-shelf processes. Despite their differing numerical methods, the two models produce broadly similar results and share similar biases in many cases. Both models reproduce many key features of observations but struggle to reproduce others, such as the high melt rates observed in the small warm-cavity ice shelves of the Amundsen and Bellingshausen seas. Several differences in model design show a particular influence on the simulations. For example, FESOM's greater topographic smoothing can alter the geometry of some ice-shelf cavities enough to affect their melt rates; this improves at higher resolution, since less smoothing is required. In the interior Southern Ocean, the vertical coordinate system affects the degree of water mass erosion due to spurious diapycnal mixing, with MetROMS' terrain-following coordinate leading to more erosion than FESOM's z coordinate. Finally, increased horizontal resolution in FESOM leads to higher basal melt rates for small ice shelves, through a combination of stronger circulation and small-scale intrusions of warm water from offshore.

A 4.5 km resolution Arctic Ocean simulation with the global multi-resolution model FESOM 1.4

Abstract. In the framework of developing a global modeling system which can facilitate modeling studies on Arctic Ocean and high- to midlatitude linkage, we evaluate the Arctic Ocean simulated by the multi-resolution Finite Element Sea ice-Ocean Model (FESOM). To explore the value of using high horizontal resolution for Arctic Ocean modeling, we use two global meshes differing in the horizontal resolution only in the Arctic Ocean (24 km vs. 4.5 km). The high resolution significantly improves the model's representation of the Arctic Ocean. The most pronounced improvement is in the Arctic intermediate layer, in terms of both Atlantic Water (AW) mean state and variability. The deepening and thickening bias of the AW layer, a common issue found in coarse-resolution simulations, is significantly alleviated by using higher resolution. The topographic steering of the AW is stronger and the seasonal and interannual temperature variability along the ocean bottom topography is enhanced in the high-resolution simulation. The high resolution also improves the ocean surface circulation, mainly through a better representation of the narrow straits in the Canadian Arctic Archipelago (CAA). The representation of CAA throughflow not only influences the release of water masses through the other gateways but also the circulation pathways inside the Arctic Ocean. However, the mean state and variability of Arctic freshwater content and the variability of freshwater transport through the Arctic gateways appear not to be very sensitive to the increase in resolution employed here. By highlighting the issues that are independent of model resolution, we address that other efforts including the improvement of parameterizations are still required.

Eddy‐Resolving Simulation of the Atlantic Water Circulation in the Fram Strait With Focus on the Seasonal Cycle

AbstractEddy driven recirculation of Atlantic Water (AW) in the Fram Strait modifies the amount of heat that reaches the Arctic Ocean, but is difficult to constrain in ocean models due to very small Rossby radius there. In this study, we explore the effect of resolved eddies on the AW circulation in a locally eddy‐resolving simulation of the global Finite‐Element‐Sea ice‐Ocean Model (FESOM) integrated for the years 2000–2009, by focusing on the seasonal cycle. An eddy‐permitting simulation serves as a control run. Our results suggest that resolving local eddy dynamics is critical to realistically simulate ocean dynamics in the Fram Strait. Strong eddy activity simulated by the eddy‐resolving model, with peak in winter and lower values in summer, is comparable in magnitude and seasonal cycle to observations from a long‐term mooring array, whereas the eddy‐permitting simulation underestimates the observed magnitude. Furthermore, a strong cold bias in the central Fram Strait present in the eddy‐permitting simulation is reduced due to resolved eddy dynamics, and AW transport into the Arctic Ocean is increased with possible implications for the Arctic Ocean heat budget. Given the good agreement between the eddy‐resolving model and measurements, it can help filling gaps that point‐wise observations inevitably leave. For example, the path of the West Spitsbergen Current offshore branch, measured in the winter months by the mooring array, is shown to continue cyclonically around the Molloy Deep in the model, representing the major AW recirculation branch in this season.

Response to Filchner–Ronne Ice Shelf cavity warming in a coupled ocean–ice sheet model – Part 1: The ocean perspective

Abstract. The Regional Antarctic ice and Global Ocean (RAnGO) model has been developed to study the interaction between the world ocean and the Antarctic ice sheet. The coupled model is based on a global implementation of the Finite Element Sea-ice Ocean Model (FESOM) with a mesh refinement in the Southern Ocean, particularly in its marginal seas and in the sub-ice-shelf cavities. The cryosphere is represented by a regional setup of the ice flow model RIMBAY comprising the Filchner–Ronne Ice Shelf and the grounded ice in its catchment area up to the ice divides. At the base of the RIMBAY ice shelf, melt rates from FESOM's ice-shelf component are supplied. RIMBAY returns ice thickness and the position of the grounding line. The ocean model uses a pre-computed mesh to allow for an easy adjustment of the model domain to a varying cavity geometry. RAnGO simulations with a 20th-century climate forcing yield realistic basal melt rates and a quasi-stable grounding line position close to the presently observed state. In a centennial-scale warm-water-inflow scenario, the model suggests a substantial thinning of the ice shelf and a local retreat of the grounding line. The potentially negative feedback from ice-shelf thinning through a rising in situ freezing temperature is more than outweighed by the increasing water column thickness in the deepest parts of the cavity. Compared to a control simulation with fixed ice-shelf geometry, the coupled model thus yields a slightly stronger increase in ice-shelf basal melt rates.

On the ventilation of Bransfield Strait deep basins

The deep basins of the Bransfield Strait (BS) are ventilated by Weddell Sea (WS) waters from different origins. Depending on the source and density, these water masses follow different routes across the complex topography near the tip of the Antarctic Peninsula and thus into the Bransfield Strait abyss.Using a global setup of the Finite Element Sea-ice Ocean Model (FESOM) we show that the WS waters found at the western WS continental shelf break have a higher influence on the short period variability of BS bottom waters than the waters present over the continental shelf.Adding passive tracers to the glacial melt water (GMW) from two different origins, Larsen Ice Shelf (LIS) and Filchner-Ronne Ice Shelf (FRIS), we show that the GMW from FRIS has a larger influence on BS bottom waters than the GMW from LIS. FRIS GMW has a higher concentration in the BS eastern basin, while LIS GMW is more abundant in the BS central basin. This duality mainly leads to the difference between BS central and eastern basins seen on the observations. This is a novel result and we believe is a significant contribution to the understanding of the BS-WS circulation and interactions.

Atlantic Water in the Nordic Seas: Locally eddy-permitting ocean simulation in a global setup

Warm and salty Atlantic Water is transported by the Norwegian Atlantic Current through the Nordic Seas. A fraction of it enters the Arctic Ocean where it contributes significantly to its heat budget. Resolving the complex circulation structure in the Nordic Seas, in particular eddies, presents a numerical challenge in ocean models. Here, we present a hindcast experiment for the years 1969–2009 with a global configuration of the Finite Element Sea-ice Ocean Model, employing a high-resolution mesh in the Nordic Seas and Arctic Ocean (4.5 km). We show that substantial improvements can be achieved in the circulation structure, hydrography and eddy kinetic energy in the Nordic Seas compared with a coarse-resolution reference run. A better represented Norwegian Atlantic Front Current (NwAFC) in the high-resolution setup leads to a reduction of a strong negative temperature bias in the eastern Nordic Seas. The Atlantic Water inflow through the Iceland-Faroe Ridge is found to be very sensitive to mesh resolution, and high resolution is required to adequately represent this inflow and the downstream NwAFC. With increased mesh resolution, the simulated ocean temperature is significantly improved at Barents Sea Opening (BSO), and the Atlantic Water volume transport in Fram Strait becomes much closer to observations in terms of both magnitude and variability. By using passive tracers, the origins of water masses at BSO and Fram Strait are identified. Our study also indicates that eddy-resolving meshes are required to further improve the representation of dynamical processes in the region, in particular in Fram Strait.

Regional sea level change in response to ice mass loss in Greenland, the West Antarctic and Alaska

AbstractBesides the warming of the ocean, sea level is mainly rising due to land ice mass loss of the major ice sheets in Greenland, the West Antarctic, and the Alaskan Glaciers. However, it is not clear yet how these land ice mass losses influence regional sea level. Here, we use the global Finite Element Sea‐ice Ocean Model (FESOM) to simulate sea surface height (SSH) changes caused by these ice mass losses and combine it with the passive ocean response to varying surface loading using the sea level equation. We prescribe rates of fresh water inflow, not only around Greenland, but also around the West Antarctic Ice Sheet and the mountain glaciers in Alaska with approximately present‐day amplitudes of 200, 100, and 50 Gt/yr, respectively. Perturbations in sea level and in freshwater distribution with respect to a reference simulation are computed for each source separately and in their combination. The ocean mass change shows an almost globally uniform behavior. In the North Atlantic and Arctic Ocean, mass is redistributed toward coastal regions. Steric sea level change varies locally in the order of several centimeters on advective timescales of decades. Steric effects to local sea level differ significantly in different coastal locations, e.g., at North American coastal regions the steric effects may have the same order of magnitude as the mass driven effect, whereas at the European coast, steric effects remain small during the simulation period.

Finite-Element Sea Ice Model (FESIM), version 2

Abstract. The Finite-Element Sea Ice Model (FESIM), used as a component of the Finite-Element Sea ice Ocean Model, is presented. Version 2 includes the elastic-viscous-plastic (EVP) and viscous-plastic (VP) solvers and employs a flux corrected transport algorithm to advect the ice and snow mean thicknesses and concentration. The EVP part also includes a modified approach proposed recently by Bouillon et al. (2013), which is characterized by an improved stability compared to the standard EVP approach. The model is formulated on unstructured triangular meshes. It assumes a collocated placement of ice velocities, mean thicknesses and concentration at mesh vertices, and relies on piecewise-linear (P1) continuous elements. Simple tests for the modified EVP and VP solvers are presented to show that they may produce very close results provided the number of iterations is sufficiently high.

Atmospheric forcing of coastal polynyas in the south-western Weddell Sea

AbstractThe development of coastal polynyas, areas of enhanced heat flux and sea ice production strongly depend on atmospheric conditions. In Antarctica, measurements are scarce and models are essential for the investigation of polynyas. A robust quantification of polynya exchange processes in simulations relies on a realistic representation of atmospheric conditions in the forcing dataset. The sensitivity of simulated coastal polynyas in the south-western Weddell Sea to the atmospheric forcing is investigated with the Finite-Element Sea ice-Ocean Model (FESOM) using daily NCEP/NCAR reanalysis data (NCEP), 6 hourly Global Model Europe (GME) data and two different hourly datasets from the high-resolution Consortium for Small-Scale Modelling (COSMO) model. Results are compared for April to August in 2007–09. The two coarse-scale datasets often produce the extremes of the data range, while the finer-scale forcings yield results closer to the median. The GME experiment features the strongest winds and, therefore, the greatest polynya activity, especially over the eastern continental shelf. This results in higher volume and export of High Salinity Shelf Water than in the NCEP and COSMO runs. The largest discrepancies between simulations occur for 2008, probably due to differing representations of the ENSO pattern at high southern latitudes. The results suggest that the large-scale wind field is of primary importance for polynya development.