Ocean Sea Ice Model Research Articles

Arctic sea‐ice export as a mechanism for cold climate events during the last deglaciation

ABSTRACTRecently it was proposed that a sudden drainage of accumulated Arctic sea ice could have caused cold climate events during the last deglaciation. To explore this mechanism, we performed numerical experiments with an atmosphere–ocean–sea ice model. In these experiments, the impact of a large flush of Arctic sea ice was compared to a reference glacial state of the climate. In our results, the sea ice flush produces a major surface freshening of the North Atlantic Ocean and a 17% weakening of the strength of the Atlantic Meridional Overturning Circulation. Together with an increase in surface albedo, this weakening in ocean circulation leads to a cooling over the North Atlantic, extending to the downwind continents. Compared to our reference glacial state, the cooling reaches 5°C and lasts about 80 years. This climate anomaly is similar in magnitude and duration to relatively short cooling events during the last deglaciation, such as the Older Dryas, the Inter‐Allerød cold period and the Preboreal Oscillation. We thus conclude that the sea‐ice flush mechanism is consistent with the occurrence of such cooling events in the North Atlantic region. However, longer cooling events such as the Younger Dryas would require additional mechanisms.

Impact of ocean vertical-mixing parameterization on Arctic sea ice and upper-ocean properties using the NEMO-SI3 model

Abstract. We evaluate the vertical turbulent-kinetic-energy (TKE) mixing scheme of the NEMO-SI3 ocean–sea-ice model in sea-ice-covered regions of the Arctic Ocean. Specifically, we assess the parameters involved in TKE mixed-layer-penetration (MLP) parameterization. This ad hoc parameterization aims to capture processes that impact the ocean surface boundary layer, such as near-inertial oscillations, ocean swells, and waves, which are often not well represented in the default TKE scheme. We evaluate this parameterization for the first time in three regions of the Arctic Ocean: the Makarov, Eurasian, and Canada basins. We demonstrate the strong effect of the scaling parameter that accounts for the presence of sea ice. Our results confirm that TKE MLP must be scaled down below sea ice to avoid unrealistically deep mixed layers. The other parameters evaluated are the percentage of energy penetrating below the mixed layer and the length scale of its decay with depth. All these parameters affect mixed-layer depth and its seasonal cycle, surface temperature, and salinity, as well as underlying stratification. Shallow mixed layers are associated with stronger stratification and fresh surface anomalies, and deeper mixed layers correspond to weaker stratification and salty surface anomalies. Notably, we observe significant impacts on sea-ice thickness across the Arctic Ocean in two scenarios: when the scaling parameter due to sea ice is absent and when the TKE mixed-layer-penetration process vanishes. In the former case, we observe an increase of several meters in mixed-layer depth, along with a reduction in sea-ice thickness ranging from 30 to 40 cm, reflecting the impact of stronger mixing. Conversely, in the latter case, we notice that a shallower mixed layer is accompanied by a moderate increase in sea-ice thickness, ranging from 10 to 20 cm, as expected from weaker mixing. Additionally, interannual variability suggests that experiments incorporating a scaling parameter based on sea-ice concentration display an increased mixed-layer depth during periods of reduced sea ice, which is consistent with observed trends. These findings underscore the influence of enhanced ocean mixing, through specific parameterizations, on the physical properties of the upper ocean and sea ice.

Split-explicit external mode solver in the finite volume sea ice–ocean model FESOM2

Abstract. A novel split-explicit (SE) external mode solver for the Finite volumE Sea ice–Ocean Model (FESOM2) and its sub-versions is presented. It is compared with the semi-implicit (SI) solver currently used in FESOM2. The split-explicit solver for the external mode utilizes a dissipative asynchronous (forward–backward) time-stepping scheme that does not require additional filtering during the barotropic sub-cycling. Its implementation with arbitrary Lagrangian–Eulerian vertical coordinates like Z star (z*) and Z tilde (z̃) is explored. The comparisons are performed through multiple test cases involving idealized and realistic global simulations. The SE solver demonstrates lower phase errors and dissipation, while maintaining a simulated mean ocean state very similar to the SI solver. The SE solver is also shown to possess better run-time performance and parallel scalability across all workloads tested.

The role of atmospheric conditions in the Antarctic sea ice extent summer minima

Abstract. Understanding the variability of Antarctic sea ice is still a challenge. After decades of modest growth, an unprecedented minimum in the sea ice extent (SIE) was registered in summer 2017, and, following years of anomalously low SIE, a new record was established in early 2022. These two memorable minima have received great attention as single cases, but a comprehensive analysis of summer SIE minima is currently lacking. Indeed, other similar events are present in the observational record, although they are minor compared to the most recent ones, and a full analysis of all summer SIE minima is essential to separate potential common drivers from event-specific dynamics in order to ultimately improve our understanding of the Antarctic sea ice and climate variability. In this work, we examine sea ice and atmospheric conditions during and before all summer SIE minima over the satellite period up to 2022. We use observations and reanalysis data and compare our main findings with results from an ocean–sea ice model (NEMO–LIM) driven by prescribed atmospheric fields from ERA5. Examining SIE and sea ice concentration (SIC) anomalies, we find that the main contributors to the summer minima are the Ross and Weddell sectors. However, the two regions play different roles, and the variability of the Ross Sea explains most of the minima, with typical negative SIE anomalies about twice as large as the ones in the Weddell Sea. Furthermore, the distribution of SIC anomalies is also different: in the Weddell Sea, they exhibit a dipolar structure, with increased SIC next to the continent and decreased SIC at the sea ice margin, while the Ross Sea displays a more homogenous decrease. We also examine the role of wintertime sea ice conditions before the summer SIE minima and find mixed results depending on the period: the winter conditions are relevant in the most recent events, after 2017, but they are marginal for previous years. Next, we consider the influence of the atmosphere on the SIE minima, which is shown to play a major role: after analyzing the anomalous atmospheric circulation during the preceding spring, we find that different large-scale anomalies can lead to similar regional prevailing winds that drive the summer minima. Specifically, the SIE minima are generally associated with dominant northwesterly anomalous winds in the Weddell Sea, while a southwesterly anomalous flow prevails in the Ross Sea. Finally, we investigate the relative contribution of dynamic (e.g., ice transport) and thermodynamic (e.g., local melting) processes to the summer minima. Our results indicate that the exceptional sea ice loss in both the Ross and Weddell sectors is dominated at the large scale by thermodynamic processes, while dynamics are also present but with a minor role.

Bathymetry-constrained ocean geostrophic currents play a key role in shaping the sea ice circulation in the Canada Basin, Arctic Ocean

Abstract Satellite-based observations and a pan-Arctic coupled sea ice-ocean model are utilized to study the effect of ocean geostrophic currents on large-scale sea ice circulation in the Canada Basin, Arctic Ocean. We find that surface winds primarily drive sea ice drifts in the west-east direction, while the geostrophic currents in the Beaufort Gyre promote north-south ice drifts. Wind fluctuations can create variable ice drifts, yet geostrophic currents respond more slowly due to their larger vertical scale, serving as a slowly-evolving conveyor belt for maintaining the anticyclonic ice circulation. It is further demonstrated that the bathymetry can regulate the movement of sea ice via constraining the expansion of ocean circulation. This mechanism is indirect in the sense that the ice is far from the seafloor. Our research underscores the necessity of considering the bathymetry-constrained geostrophic currents in understanding Arctic sea ice dynamics. With the rapid retreat of Arctic sea ice, the multi-scale interactions between ice drifts and ocean currents may have significant implications for the Arctic ecosystem, climate, and shipping corridors.

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

Dynamical reconstruction of the upper-ocean state in the central Arctic during the winter period of the MOSAiC expedition

Abstract. This paper presents a methodological tool for dynamic reconstruction of the state of the ocean, based, as an example, on observations from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) experiment. The data used in this study were collected in the Amundsen Basin between October 2019 and January 2020. Analysing observational data to assess tracer field and upper-ocean dynamics is highly challenging when measurement platforms drift with the ice pack due to continuous drift speed and direction changes. We have equipped the new version of the coastal branch of the global Finite-volumE sea ice–Ocean Model (FESOM-C) with a nudging method. Model nudging was carried out assuming a quasi-steady state. Overall, the model can reproduce the lateral and vertical structure of the temperature, salinity, and density fields, which allows for projecting dynamically consistent features of these fields onto a regular grid. We identify two separate depth ranges of enhanced eddy kinetic energy located around two maxima in buoyancy frequency: the depth of the upper halocline and the depth of the warm (modified) Atlantic Water. Simulations reveal a notable decrease in surface layer salinity and density in the Amundsen Basin towards the north but no significant gradient from east to west. However, we find a mixed-layer deepening from east to west, with a 0.084 m km−1 gradient at 0.6 m km−1 standard deviation, compared to a weak deepening from south to north. The model resolves several stationary eddies in the warm Atlantic Water and provides insights into the associated dynamics. The model output can be used to further analyse the thermohaline structure and related dynamics associated with mesoscale and submesoscale processes in the central Arctic, such as estimates of heat fluxes or mass transport. The developed nudging method can be utilized to incorporate observational data from a diverse set of instruments and for further analysis of data from the MOSAiC expedition.

How Much Can Riverine Biogeochemical Fluxes Affect the Arctic Ocean Acidification?

AbstractArctic rivers carry not only large amounts of freshwater but also biogeochemical materials into the ocean and play important roles in Arctic ocean acidification (OA). This study quantitatively evaluated the effects of riverine biogeochemical fluxes (R‐BGC; carbon and nutrients) on the Arctic marine carbonate system and OA using multi‐decadal experiments (1979–2018) with a pan‐Arctic sea ice–ocean model. Improved initial and lateral boundary conditions of carbonate properties, observation‐based riverine biogeochemical data, and land model‐based interannually varying riverine freshwater discharge were adopted to enable more realistic experiments. The model simulated negative trends in aragonite saturation state (Ω) and pH in most regions of Arctic Ocean regardless of R‐BGC. The increased riverine freshwater promoted more OA through the higher dilution effect. Compared to the experiment with riverine discharge of only freshwater, the inclusion of R‐BGC caused positive anomalies in Ω and pH (by ∼0.14 and ∼0.03 in central basins, and by ∼0.15 and ∼0.06 in shelf seas, respectively). In the central basins, these anomalies were caused mostly by carbon (total alkalinity and dissolved inorganic carbon) of the R‐BGC. In the shelf seas, nutrient (nitrate and silicate) fluxes also contributed ∼14% and ∼32% of the anomalies owing to the enhanced primary production and a corresponding reduction in seawater pCO2. R‐BGC mitigated OA (ΔΩ = −1.53 × 10−3 year−1 and ΔpH = −0.56 × 10−3 year−1) in regions where riverine freshwater was accumulated (i.e., the Canada Basin, Chukchi Cap, Eurasian Basin, and East Siberian Sea). This study stressed the importance of including R‐BGC for OA model projection.

Advanced Sea Ice Modeling for Short-Term Forecasting for Alaska’s Coasts

Abstract In Alaska’s coastal environment, accurate information of sea ice conditions is desired by operational forecasters, emergency managers, and responders. Complicated interactions among atmosphere, waves, ocean circulation, and sea ice collectively impact the ice conditions, intensity of storm surges, and flooding, making accurate predictions challenging. A collaborative work to build the Alaska Coastal Ocean Forecast System established an integrated storm surge, wave, and sea ice model system for the coasts of Alaska, where the verified model components are linked using the Earth System Modeling Framework and the National Unified Operational Prediction Capability. We present the verification of the sea ice model component based on the Los Alamos Sea Ice Model, version 6. The regional, high-resolution (3 km) configuration of the model was forced by operational atmospheric and ocean model outputs. Extensive numerical experiments were conducted from December 2018 to August 2020 to verify the model’s capability to represent detailed nearshore and offshore sea ice behavior, including landfast ice, ice thickness, and evolution of air–ice drag coefficient. Comparisons of the hindcast simulations with the observations of ice extent presented the model’s comparable performance with the Global Ocean Forecast System 3.1 (GOFS3.1). The model’s skill in reproducing landfast ice area significantly outperformed GOFS3.1. Comparison of the modeled sea ice freeboard with the Ice, Cloud, and Land Elevation Satellite-2 product showed a mean bias of −4.6 cm. Daily 5-day forecast simulations for October 2020–August 2021 presented the model’s promising performance for future implementation in the coupled model system. Significance Statement Accurate sea ice information along Alaska’s coasts is desired by the communities for preparedness of hazardous events, such as storm surges and flooding. However, such information, in particular predicted conditions, remains to be a gap. This study presents the verification of the state-of-art sea ice model for Alaska’s coasts for future use in the more comprehensive coupled model system where ocean circulation, wave, and sea ice models are integrated. The model demonstrates comparable performance with the existing operational ocean–ice coupled model product in reproducing overall sea ice extent and significantly outperformed it in reproducing landfast ice cover. Comparison with the novel satellite product presented the model’s ability to capture sea ice freeboard in the stable ice season.

Eastern Boundary Upwelling Systems in Ocean–Sea Ice Simulations Forced by CORE and JRA55-do: Mean State and Variability at the Surface

Abstract In this paper we summarize improvements in climate model simulation of eastern boundary upwelling systems (EBUS) when changing the forcing dataset from the Coordinated Ocean-Ice Reference Experiments (CORE; ∼2° winds) to the higher-resolution Japanese 55-year Atmospheric Reanalysis for driving ocean–sea ice models (JRA55-do, ∼0.5°) and also due to refining ocean grid spacing from 1° to 0.1°. The focus is on sea surface temperature (SST), a key variable for climate studies, and which is typically too warm in climate model representation of EBUS. The change in forcing leads to a better-defined atmospheric low-level coastal jet, leading to more equatorward ocean flow and coastal upwelling, both in turn acting to reduce SST over the upwelling regions off the west coast of North America, Peru, and Chile. The refinement of ocean resolution then leads to narrower and stronger alongshore ocean flow and coastal upwelling, and the emergence of strong across-shore temperature gradients not seen with the coarse ocean model. Off northwest Africa the SST bias mainly improves with ocean resolution but not with forcing, while in the Benguela, JRA55-do with high-resolution ocean leads to lower SST but a substantial bias relative to observations remains. Reasons for the Benguela bias are discussed in the context of companion regional ocean model simulations. Finally, we address to what extent improvements in mean state lead to changes to the monthly to interannual variability. It is found that large-scale SST variability in EBUS on monthly and longer time scales is largely governed by teleconnections from climate modes and less sensitive to model resolution and forcing than the mean state.

Assimilation of satellite swaths versus daily means of sea ice concentration in a regional coupled ocean–sea ice model

Abstract. Operational forecasting systems routinely assimilate daily means of sea ice concentration (SIC) from microwave radiometers in order to improve the accuracy of the forecasts. However, the temporal and spatial averaging of the individual satellite swaths into daily means of SIC entails two main drawbacks: (i) the spatial resolution of the original product is blurred (especially critical in periods with strong sub-daily sea ice movement), and (ii) the sub-daily frequency of passive microwave observations in the Arctic are not used, providing less temporal resolution in the data assimilation (DA) analysis and, therefore, in the forecast. Within the SIRANO (Sea Ice Retrievals and data Assimilation in NOrway) project, we investigate how challenges (i) and (ii) can be avoided by assimilating individual satellite swaths (level 3 uncollated) instead of daily means (level 3) of SIC. To do so, we use a regional configuration of the Barents Sea (2.5 km grid) based on the Regional Ocean Modeling System (ROMS) and the Los Alamos Sea Ice Model (CICE) together with the ensemble Kalman filter (EnKF) as the DA system. The assimilation of individual swaths significantly improves the EnKF analysis of SIC compared to the assimilation of daily means; the mean absolute difference (MAD) shows a 10 % improvement at the end of the assimilation period and a 7 % improvement at the end of the 7 d forecast period. This improvement is caused by better exploitation of the information provided by the SIC swath data, in terms of both spatial and temporal variance, compared to the case when the swaths are combined to form a daily mean before assimilation.

CD-type discretization for sea ice dynamics in FESOM version 2

Abstract. Two recently proposed variants of CD-type discretizations of sea ice dynamics on triangular meshes are implemented in the Finite-VolumE Sea ice–Ocean Model (FESOM version 2). The implementations use the finite element method in spherical geometry with longitude–latitude coordinates. Both are based on the edge-based sea ice velocity vectors but differ in the basis functions used to represent the velocities. The first one uses nonconforming linear (Crouzeix–Raviart) basis functions, and the second one uses continuous linear basis functions on sub-triangles obtained by splitting parent triangles into four smaller triangles. Test simulations are run to show how the performance of the new discretizations compares with the A-grid discretization using linear basis functions. Both CD discretizations are found to simulate a finer structure of linear kinematic features (LKFs). Both show some sensitivity to the representation of scalar fields (sea ice concentration and thickness). Cell-based scalars lead to a finer LKF structure for the first CD discretization, but the vertex-based scalars may be advantageous in the second case.

Mechanisms of Offshore Solid and Liquid Freshwater Flux from the East Greenland Current

Abstract The mechanisms that control the export of freshwater from the East Greenland Current, in both liquid and solid form, are explored using an idealized numerical model and scaling theory. A regional, coupled ocean–sea ice model is applied to a series of calculations in which key parameters are varied and the scaling theory is used to interpret the model results. The offshore ice flux, occurring in late winter, is driven primarily by internal stresses and is most sensitive to the thickness of sea ice on the shelf coming out of Fram Strait and the strength of alongshore winds over the shelf. The offshore liquid freshwater flux is achieved by eddy fluxes in late summer while there is an onshore liquid freshwater flux in winter due to the ice–ocean stress, resulting in only weak annual mean flux. The scaling theory identifies the key nondimensional parameters that control the behavior and reproduces the general parameter dependence found in the numerical model. Climate models predict that winds will increase and ice export from the Arctic will decrease in the future, both of which will lead to a decrease in the offshore flux of sea ice, while the influence on liquid freshwater may increase or decrease, depending on the relative changes in the onshore Ekman transport and offshore eddy fluxes. Additional processes that have not been considered here, such as more complex topography and synoptic wind events, may also contribute to cross-shelf exchange. Significance Statement The purpose of this study is to provide a basic understanding of what controls the flux of sea ice and low-salinity water from the East Greenland shelf into the interior of the Greenland and Iceland Seas. This is a potentially important process since it has been shown that sufficient freshening of the surface waters in the interior of the Nordic seas can inhibit deep convection and the associated air–sea heat flux and water mass transformation. A combination of idealized computer models and basic theory indicates that the fluxes of liquid and solid freshwater are controlled by different mechanisms and occur at different times of the year. Accurate representation in climate models will require representation of small-scale processes such as mesoscale eddies and gradients of ice thickness across the shelf.

Eddy activity in the Arctic Ocean projected to surge in a warming world

Ocean eddies play a critical role in climate and marine life. In the rapidly warming Arctic, little is known about how ocean eddy activity will change because existing climate models cannot resolve Arctic Ocean mesoscale eddies. Here, by employing a next-generation global sea ice–ocean model with kilometre-scale horizontal resolution in the Arctic, we find a surge of eddy kinetic energy in the upper Arctic Ocean, tripling on average in a four-degree-warmer world. The driving mechanism behind this surge is an increase in eddy generation due to enhanced baroclinic instability. Despite the decline of sea ice, eddy killing (a process in which eddies are dampened by sea ice and winds) will not weaken in its annual mean effect in the considered warming scenario. Our study suggests the importance of adequately representing Arctic eddy activity in climate models for understanding the impacts of its increase on climate and ecosystems.

The Southern Ocean Freshwater Input from Antarctica (SOFIA) Initiative: scientific objectives and experimental design

Abstract. As the climate warms, the grounded ice sheet and floating ice shelves surrounding Antarctica are melting and releasing additional freshwater into the Southern Ocean. Nonetheless, almost all existing coupled climate models have fixed ice sheets and lack the physics required to represent the dominant sources of Antarctic melt. These missing ice dynamics represent a key uncertainty that is typically unaccounted for in current global climate change projections. Previous modelling studies that have imposed additional Antarctic meltwater have demonstrated regional impacts on Southern Ocean stratification, circulation, and sea ice, as well as remote changes in atmospheric circulation, tropical precipitation, and global temperature. However, these previous studies have used widely varying rates of freshwater forcing, have been conducted using different climate models and configurations, and have reached differing conclusions on the magnitude of meltwater–climate feedbacks. The Southern Ocean Freshwater Input from Antarctica (SOFIA) initiative brings together a team of scientists to quantify the climate system response to Antarctic meltwater input along with key aspects of the uncertainty. In this paper, we summarize the state of knowledge on meltwater discharge from the Antarctic ice sheet and ice shelves to the Southern Ocean and explain the scientific objectives of our initiative. We propose a series of coupled and ocean–sea ice model experiments, including idealized meltwater experiments, historical experiments with observationally consistent meltwater input, and future scenarios driven by meltwater inputs derived from stand-alone ice sheet models. Through coordinating a multi-model ensemble of simulations using a common experimental design, open data archiving, and facilitating scientific collaboration, SOFIA aims to move the community toward better constraining our understanding of the climate system response to Antarctic melt.

Improvements in September Arctic Sea Ice Predictions Via Assimilation of Summer CryoSat‐2 Sea Ice Thickness Observations

AbstractBecause of a spring predictability barrier, the seasonal forecast skill of Arctic summer sea ice is limited by the availability of melt‐season sea ice thickness (SIT) observations. The first year‐round SIT observations, retrieved from CryoSat‐2 from 2011 to 2020, are assimilated into the GFDL ocean–sea ice model. The model's SIT anomaly field is brought into significantly better agreement with the observations, particularly in the Central Arctic. Although the short observational period makes forecast assessment challenging, we find that the addition of May–August SIT assimilation improves September local sea ice concentration (SIC) and extent forecasts similarly to SIC‐only assimilation. Although most regional forecasts are improved by SIT assimilation, the Chukchi Sea forecasts are degraded. This degradation is likely due to the introduction of negative correlations between September SIC and earlier SIT introduced by SIT assimilation, contrary to the increased correlations found in other regions.

Response of the Arctic sea ice–ocean system to meltwater perturbations based on a one-dimensional model study

Abstract. A one-dimensional coupled sea ice–ocean model is used to investigate how the Arctic Ocean stratification and sea ice respond to changes in meltwater. In the control experiments, the model is capable of accurately simulating seasonal changes in the upper-ocean stratification structure compared with observations, and the results suggest that ocean stratification is important for ice thickness development during the freezing season. The sensitivity experiments reveal the following: (1) a decrease in meltwater release weakens ocean stratification and creates a deeper, higher-salinity mixed layer. (2) Meltwater reduced ice melting by 17 % by strengthening ocean stratification. (3) The impact of meltwater released during the previous melting season on ice growth in winter depends on the strength of stratification. After removing all the meltwater during the summer, ice formation in areas with strong stratification increased by 12 % during the winter, while it decreased by 43 % in areas with weak stratification. (4) In some areas of the Nansen Basin where stratification is nearly absent, the warm Atlantic Water can reach the ice directly in early spring, leading to early melting of the sea ice in winter if all meltwater is removed from the model. These findings contribute to our understanding of the complex interactions between ocean stratification, meltwater and sea ice growth and have important implications for climate models and future change prediction in the Arctic.

Basin-scale tracer replacement timescales in a one-degree global OGCM

This study investigates basin-scale tracer replacement timescales of the two polar oceans and the Atlantic, Indian, and Pacific Oceans using a one-degree global ocean-sea ice model that represents oceans under the largest Antarctic ice shelves, the Filchner-Ronne and Ross Ice Shelf (FRIS and RIS). After a long spin-up with present-day surface boundary conditions, we confirm that the model has a typical representation of wind-driven and thermohaline circulations in one-degree ocean models. We use virtual passive tracers placed in the five oceans and examine the behavior of the passive tracers to estimate the tracer replacement timescales and pathways of the basin-scale ocean waters. Replacement timescales in the polar oceans (114 years for the Southern Ocean and 109 years for the Arctic Ocean) are found to be shorter than those in the three oceans (217 years for the Atlantic Ocean, 163 years for the Indian Ocean, and 338 years for the Pacific Ocean). The Southern Ocean tracer has two clear pathways to the Northern Hemisphere: the surface route in the Atlantic Ocean and the bottom route in the Pacific and Indian Oceans. This surface route is a rapid conduit to transport the Southern Ocean signal to the North Atlantic and Arctic Oceans. The Atlantic Ocean tracer is transported to both polar regions along the North Atlantic Current and the Antarctic Circumpolar Current (ACC). The tracer experiments clearly demonstrate that Atlantic Meridional Overturning Circulation (AMOC) plays a vital role in transporting the water masses in the Atlantic and Arctic Oceans to the Southern Ocean. The southward flow of the AMOC at the intermediate depths carries the northern waters to the ACC region, and then the water spreads over the Southern Ocean along the eastward-flowing ACC. The decay timescales of water in the ice-shelf cavities exposed to the water outside the Southern Ocean are estimated to be approximately 150 years for both the FIRS and RIS. The decay timescales in the Antarctic coastal region are short at the surface and long in the deep layers, with a noticeable reduction in the areas where ACC flows southward toward the Antarctic continent.

Sensitivity of Antarctic Shelf Waters and Abyssal Overturning to Local Winds

Abstract Winds around the Antarctic continental margin are known to exert a strong control on the local ocean stratification and circulation. However, past work has largely focused on the ocean response to changing winds in limited regional sectors and the circumpolar dynamical response to polar wind change remains uncertain. In this work, we use a high-resolution global ocean–sea ice model to investigate how dense shelf water formation and the temperature of continental shelf waters respond to changes in the zonal and meridional components of the polar surface winds. Increasing the zonal easterly wind component drives an enhanced southward Ekman transport in the surface layer, raising sea level over the continental shelf and deepening coastal isopycnals. The downward isopycnal movement cools the continental shelf, as colder surface waters replace warmer waters below. However, in this model the zonal easterly winds do not impact the strength of the abyssal overturning circulation, in contrast to past idealized model studies. Instead, increasing the meridional wind speed strengthens the abyssal overturning circulation via a sea ice advection mechanism. Enhanced offshore meridional wind speed increases the northward export of sea ice, resulting in decreased sea ice thickness over the continental shelf. The reduction in sea ice coverage leads to increased air–sea heat loss, sea ice formation, brine rejection, dense shelf water formation, and abyssal overturning circulation. Increasing the meridional winds causes warming at depth over most of the continental shelf, due to a heat advection feedback associated with the enhanced overturning circulation.

Ocean biogeochemistry in the coupled ocean–sea ice–biogeochemistry model FESOM2.1–REcoM3

Abstract. The cycling of carbon in the oceans is affected by feedbacks driven by changes in climate and atmospheric CO2. Understanding these feedbacks is therefore an important prerequisite for projecting future climate. Marine biogeochemistry models are a useful tool but, as with any model, are a simplification and need to be continually improved. In this study, we coupled the Finite-volumE Sea ice–Ocean Model (FESOM2.1) to the Regulated Ecosystem Model version 3 (REcoM3). FESOM2.1 is an update of the Finite-Element Sea ice–Ocean Model (FESOM1.4) and operates on unstructured meshes. Unlike standard structured-mesh ocean models, the mesh flexibility allows for a realistic representation of small-scale dynamics in key regions at an affordable computational cost. Compared to the previous coupled model version of FESOM1.4–REcoM2, the model FESOM2.1–REcoM3 utilizes a new dynamical core, based on a finite-volume discretization instead of finite elements, and retains central parts of the biogeochemistry model. As a new feature, carbonate chemistry, including water vapour correction, is computed by mocsy 2.0. Moreover, REcoM3 has an extended food web that includes macrozooplankton and fast-sinking detritus. Dissolved oxygen is also added as a new tracer. In this study, we assess the ocean and biogeochemical state simulated with FESOM2.1–REcoM3 in a global set-up at relatively low spatial resolution forced with JRA55-do (Tsujino et al., 2018) atmospheric reanalysis. The focus is on the recent period (1958–2021) to assess how well the model can be used for present-day and future climate change scenarios on decadal to centennial timescales. A bias in the global ocean–atmosphere preindustrial CO2 flux present in the previous model version (FESOM1.4–REcoM2) could be significantly reduced. In addition, the computational efficiency is 2–3 times higher than that of FESOM1.4–REcoM2. Overall, it is found that FESOM2.1–REcoM3 is a skilful tool for ocean biogeochemical modelling applications.