Ocean circulation and sea ice distribution in a finite element global sea ice–ocean model

A newly developed global Finite Element Sea Ice–Ocean Model (FESOM) is presented. The ocean component is based on the Finite Element model of the North Atlantic (FENA) but has been substantially updated and extended. In addition to a faster realization of the numerical code, state-of-the-art parameterizations of subgrid-scale processes have been implemented. A Redi/GM scheme is employed to parameterize the effects of mesoscale eddies on lateral 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 ice–model has been developed and coupled to the ocean component. Sea ice thermodynamics have been derived from the standard AWI sea ice model featuring a prognostic snow layer but neglecting internal heat storage. The dynamic part offers the viscous-plastic and elastic-viscous-plastic rheologies. All model components are discretized on a triangular/tetrahedral grid with a continuous, conforming representation of model variables. The coupled model is run in a global configuration and forced with NCEP daily atmospheric reanalysis data for 1948–2007. Results are analysed with a slight focus on the Southern Hemisphere. Many aspects of sea ice distribution and hydrography are found to be in good agreement with observations. As in most coarse-scale models, Gulf Stream transport is underestimated, but transports of the Kuroshio and the Antarctic Circumpolar Current appear realistic. The seasonal cycles of Arctic and Antarctic sea ice extents and Antarctic sea ice thickness are well captured; long- and short-term variability of ice coverage is found to be reproduced realistically in both hemispheres. The coupled model is now ready to be used in a wide range of applications.     

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

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.     

A global coupled sea ice–ocean model

A new sea ice model, GELATO, was developed at Centre National de Recherches Météorologiques (CNRM) and coupled with OPA global ocean model. The sea ice model includes elastic–viscous–plastic rheology, redistribution of ice floes of different thicknesses, and it also takes into account leads, snow cover and snow ice formation. Climatologies of atmospheric surface parameters are used to perform a 20-year global ocean–sea ice simulation, in order to compute surface heat fluxes from diagnosed sea ice or ocean surface temperature. A surface salinity restoring term is applied only to ocean grid cells with no sea ice to avoid significant surface salinity drifts, but no correction of sea surface temperature is introduced. In the Arctic the use of an ocean model substantially improves the representation of sea ice, and particularly of the ice edge in all seasons, as advection of heat and salt can be more accurately accounted for than in the case of, for example, a sea ice–ocean mixed layer model. In contrast, in the Antarctic, a region where ocean convective processes bear a much stronger influence in shaping sea ice characteristics, a better representation of convection and probably of sea ice (for example, of frazil sea ice, brine rejection) would be needed to improve the simulation of the annual cycle of the sea ice cover. The effect of the inclusion of several ice categories in the sea ice model is assessed by running a sensitivity experiment in which only one category of sea ice is considered, along with leads. In the Arctic, such an experiment clearly shows that a multicategory sea ice model better captures the position of the sea ice edge and yields much more realistic sea ice concentrations in most of the region, which is in agreement with results from Bitz et al. [J. Geophys. Res. 106 (C2) (2001) 2441–2463].     

Reproduction of the large-scale state of water and sea ice in the arctic ocean in 1948–2002: Part I. Numerical model

Calculations were performed using a model of the combined circulation of the Atlantic Ocean (from 20° S), the Arctic Ocean, and the Bering Sea with a resolution of 0.25° by latitude and longitude for 1958–2006. The results are compared with observational data and results obtained by other models. Model estimates were obtained for the evolution of the Atlantic water inflow into the Arctic basin through the Fram Strait and the Barents Sea. Increased transports of Atlantic water inflow into the Arctic basin were found for the first half of the 1990s and 2004–2006. The relation between Atlantic water transports into the Arctic basin and variations in the North Atlantic oscillation is shown. A positive trend of Atlantic water inflow into the Arctic basin through the Fram Strait (0.061 Sv per year) was revealed. The evolution of the freshwater-layer thickness in the Beaufort Circulation (BC) is considered. There are three periods of its increased values combined with the increased anticyclonic vorticity of BC currents: the 1960s, the 1980s, and from 1999 until now. The model estimate for a statistical mean timescale of the cycle of freshwater concentration and sink from the BC is 16 years, which is close to currently existing estimates. The evolution of anticyclonic vorticity of currents leads the variations in the freshwater-layer thickness of the BC by 1.75 years. Since the mid-1970s, there have been long positive trends of both the freshwater-layer thickness and anticyclonic vorticity of currents in the BC. In the same time period, there has been a satellite-registered negative trend in the ice area in the Arctic, which was reproduced by the model.     

Wave–ice interactions in the marginal ice zone. Part 2: Numerical implementation and sensitivity studies along 1D transects of the ocean surface

The theoretical foundation of a wave–ice interaction model is reported in Part 1 of this study. The model incorporates attenuation of ocean surface waves by sea ice floes and the concomitant breaking of the floes by waves that determines the structure of the marginal ice zone (MIZ). A numerical implementation of the method is presented here. Convergence of the numerical method is demonstrated, as temporal and spatial grids are refined. A semi-analytical method, which does not require time-stepping, is also developed to validate the numerical results, when dispersion is neglected. The wave energy lost during ice breakage is parameterized, as part of the numerical method. Sensitivity studies are conducted in relation to the energy loss and also dispersive effects, the choice of the attenuation model, the properties of the wave field, and sea ice properties such as concentration, thickness and breaking strain. Example simulations intended to represent conditions in the Fram Strait in 2007, which exploit reanalyzed wave and ice model data, are shown to conclude the results section. These are compared to estimates of MIZ widths based on a concentration criteria, and obtained from remotely-sensed passive microwave images.     

Impact of ice cover in the Arctic on ocean–atmosphere turbulent heat fluxes

The impact of spatiotemporal variability of the ice-covered area in the Arctic on the value and interannual dynamics of turbulent heat fluxes on the ocean–atmosphere border is considered. An expected inverse dependence of the heat fluxes integrated over the Arctic area and the area of ice is not detected. The largest interannual oscillations of heat fluxes from the ocean to the atmosphere are timed to the varying position of the ice edge and, to a lesser extent, are connected with total area of ice. The role of the marginal ice zone in oceanic heat transfer is analyzed. In particular, it is shown that while moving along the marginal zone from the ice-free surface to the surface with an ice concentration of 0.8, latent and sensible heat fluxes are reduced by a factor of 2.5–3.     

A coupled ice–ocean model for the Greenland,Iceland and Norwegian Seas

Simulations from a coupled ice–ocean model that highlight the importance of synoptic forcing on sea-ice dynamics are described. The ocean model is a non-hydrostatic primitive equation model coupled to a dynamic thermodynamic sea ice model. The ice modelling sensitivity study presented here is part of an ongoing research programme to define the role played by sea ice in the energy balance of the Greenland Sea. The different categories of sea ice found in the subpolar regions are simulated through the use of equations for thin ice, thick ice and the Marginal Ice Zone. A basin scale numerical model of the Greenland, Iceland and Norwegian Seas has a horizontal resolution of 20 km and a vertical grid spacing of 50 m. This resolution is adequate for resolving the mesoscale topographic structures known to control the circulation in this region. The spin-up reproduces the main features of the circulation, including the cyclonic gyres in the Norwegian and Greenland Basins and Iceland Plateau. Topographic steering of the flow is evident. The baroclinic Rossby radius of deformation is between 5 and 10 km so that the model is not eddy-resolving. The coupled ice–ocean model was run for a period of two weeks. The influence of horizontal resolution of the atmospheric model was tested by comparing simulations using six hourly wind fields from the ECMWF with those generated using six hourly fields from a HIRLAM, with horizontal resolutions of 1° and 0.18° respectively. The simulations show reasonable agreement with satellite ice compactness data and data of ice transports across sections at 79°N, 75°N and Denmark Strait.     

Sea-ice retreat in the Sea of Okhotsk and the ice–ocean albedo feedback effect on it

Sea-ice retreat processes are examined in the Sea of Okhotsk. A heat budget analysis in the sea-ice zone shows that net heat flux from the atmosphere at the water surface is about 77 W m⁻² on average in the active ice melt season (April) due to large solar heating, while that at the ice surface is about 12 W m⁻² because of the difference in surface albedo. The temporal variation of the heat input into the upper ocean through the open water fraction corresponds well to that of the latent heat required for ice retreat. These results suggest that heat input into the ice–upper ocean system from the atmosphere mainly occurs at the open water fraction, and this heat input into the upper ocean is an important heat source for ice melting. The decrease in ice area in the active melt season (April) and the geostrophic wind just before the melt season (March) show a correlation: the decrease is large when the offshoreward wind is strong. This relationship can be explained by the following process. Once ice concentration is decreased (increased) by the offshoreward (onshoreward) wind just before 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. This positive feedback is regarded as the ice–ocean albedo feedback, and explains in part the large interannual variability of the ice cover in the ice melt season.     

A three-dimensional ocean general circulation model for mesoscale eddies——Ⅱ Diagnostic analysis

A three-dimensional density field associated with mesoscaie unstable waves generated by the 3-D, primitive-equation model (Wang and Ikeda, 1996) is provided to the quasi-geostrophic pressure tendency and ω-equations, and to the (ageostrophic) Q-vector equation. Diagnostic analyses, analogous to the approaches in meteorology: ω-equation and Q-vector method, are for the first time developed to examine the mesoscaie dynamical processes and mechanisms of the unstable waves propagating in the mid-latitude ocean. The weaknesses and strengths of these two diagnostic approaches are evaluated and compared to the model results. The Q-vector method is then recommended to diagnose the vertical motion associated with the mesoscaie dynamics from a hydrographic CTD (conductivity-temperature-depth) array, while the quasi-geostrophic equations produce some small-scale features (errors) in the diagnosed fields.     

Numerical simulations of wave–current flow in an ocean basin

Wave–current flow is a phenomenon that is present in many practical engineering situations. Over the past several decades, this type of flow has been increasingly investigated under controlled laboratory conditions. This paper presents a numerical study of wave–current flow in the ocean basin of the LabOceano (COPPE/UFRJ). A homogeneous multiphase model based on the RANS equations and the k–ɛ turbulence model implemented in ANSYS-CFX code were used. A cross section of the ocean basin was represented. A regular wave with a height of 0.08 m and a period of 1.80 s (i.e., a wave steepness of H/L = 0.016), propagating on favourable currents, was simulated. The behaviour of the free surface elevation over time and the streamlines along the basin for wave and wave–current flows were presented. The numerical results were compared to the non-viscous theory given by the Rayleigh equation applied to the problem of wave–current interaction. Good agreement was found between the wave length estimated by the numerical results and the analytical solutions, with a deviation of less than 2%.     

A parameterization of ice shelf–ocean interaction for climate models

Model results from a regional model (BRIOS) of the Southern Ocean that includes ice shelf cavities and the interaction between ocean and ice shelves are used to derive a simple parameterization for ice shelf melting and the corresponding fresh water flux in large-scale ocean climate models. The parameterization assumes that the heat loss and fresh water gain due to the ice shelves are proportional to the difference in freezing temperature at the ice shelf edge base and the oceanic temperature on the shelf/slope area of the adjacent ocean as well as an effective area of interaction. This area is proportional to the along-shelf width of ice shelf and an effective cross-shelf distance, which turns out to be rather uniform (5–15 km) for a variety of different ice shelves. The proposed parameterization is easy to implement and valid for a wide range of circumstances. An application of the proposed scheme in a global ice ocean model (CLIO) supports our hypothesis that it can be used successfully and improves both the ocean and sea ice component of the model. This parameterization should also be used in models of the climate system that include a coupling between an ice sheet and an oceanic component.     

Variability of the Bering Sea circulation in the period 1992–2010

Sea surface height anomalies observed by satellites in 1992–2010 are combined with monthly climatologies of temperature and salinity to estimate circulation in the southern Bering Sea. The estimated surface and deep currents are consistent with independent velocity observations by surface drifters and Argo floats parked at 1,000?m. Analysis reveals 1–3-Sv interannual transport variations of the major currents with typical intra-annual variability of 3–7?Sv. On the seasonal scale, the Alaskan Stream transport is well correlated with the Kamchatka (0.81), Near Strait (0.53) and the Bering Slope (0.37) currents. Lagged correlations reveal a gradual increase of the time the lags between the transports of the Alaskan Stream, the Bering Slope Current and the Kamchatka Current, supporting the concept that the Bering Sea basin is ventilated by the waters carried by the Alaskan Stream south of the Aleutian Arc and by the flow through the Near Strait. Correlations of the Bering Sea currents with the Bering Strait transport are dominated by the seasonal cycle. On the interannual time scale, significant negative correlations are diagnosed between the Near Strait transport and the Bering Slope and Alaskan Stream currents. Substantial correlations are also diagnosed between the eddy kinetic energy and Pacific Decadal Oscillation.     

Numerical simulation of the world ocean circulation and its climatic variability for 1948–2007 using the INMOM

The results of simulating global ocean circulation and its interannual variability in 1948–2007 using INM RAS ocean general circulation model INMOM (Institute of Numerical Mathematics Ocean Model) are presented. One of the INMOM versions is also used for the Black Sea dynamics simulation. The CORE datasets were used to set realistic atmospheric forcing. Sea ice area decrease by 2007 was reproduced in the Arctic Ocean that is in good agreement with observations. The interdecadal climatic variability was revealed with significant decrease of Atlantic thermohaline circulation (ATHC) and meridional heat transport (MHT) in North Atlantic (NA) since the late 1990’s. MHT presents decrease of heat transport from NA to the atmosphere since the mid-1990’s. Therefore the negative feedback is revealed in the Earth climate system that leads to reducing of climate warming caused primarily by anthropogenic factor for the last decades. Long-term variability (60 years) of ATHC is revealed as well which influences NA thermal state with 10 year delay. The assumption is argued that this mechanism can make a contribution in the ATHC own long-term variability.     

Can an atmospherically forced ocean model accurately simulate sea surface temperature during ENSO events?

The performance of an atmospherically forced ocean general circulation model (OGCM) in simulating daily and monthly sea surface temperature (SST) is examined during the historical El Niño Southern Oscillation (ENSO) events during the time period 1993–2003. For this purpose, we use the HYbrid Coordinate Ocean Model (HYCOM) configured for the North Pacific north of 20°S at a resolution of ≈9 km. There is no assimilation of (or relaxation to) SST data and no date-specific assimilation of any data type. The ability of the model in simulating temporal variations of SST anomalies is discussed by comparing model results with two satellite-based SST products. The HYCOM simulation gives a basin-averaged monthly mean bias of 0.3 °C and rms difference of 0.6 °C over the North Pacific Ocean during 1993–2003. While the model is able to simulate SST anomalies with mean biases  <0.5 °C  in comparison to observations during most of the ENSO events, limitations in the accuracy of atmospheric forcing (specifically, net short-wave radiation) have some influence on the accuracy of simulations. This is specifically demonstrated during the 1998 transition period from El Niño to La Niña, when a record large SST drop of  ≈7 °C  occurred in the eastern equatorial Pacific Ocean.     

The concentrations of “aluminosilicates” in particulates from some surface waters of the world ocean

The concentrations of aluminium, and of “aluminosilicates”, have been determined in 56 centrifuge-collected total-particulate samples from “open-ocean” and “coastal” surface waters.The average percentages of “aluminosilicates” in the particulates have been tentatively related to the supply of land-derived solids to the oceans. The China Sea, which receives a relatively large input of river-transported solids from the mainland of Asia, has particulates which contain an average of 62% “aluminosilicate” material. The eastern margins of the North Atlantic have a relatively large input of wind-transported solids, and the particulates contain an average of 9.6% “aluminosilicates”. “Open-ocean” areas remote from the land masses have particulates in which, on average, “aluminosilicate” material constitutes ? 5%.     

A three-dimensional ocean general circulation model for mesoscale eddies——Ⅰ. Meander simulation and linear growth rate

Ａｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｏｃｅａｎｇｅｎｅｒａｌｃｉｒｃｕｌａｔｉｏｎｍｏｄｅｌｆｏｒｍｅｓｏｓｃａｌｅｅｄｄｉｅｓ──ⅠＭｅａｎｄｅｒｓｉｍｕｌａｔｉｏｎａｎｄｌｉｎｅａｒｇｒｏｗｔｈｒａｔｅ￥ＷａｎｇＪｉａａｎｄＭｏｔｏＩｋｅｄａ（Ｒ...     

A model of the formation and development of ocean shear wave

Ａｍｏｄｅｌｏｆｔｈｅｆｏｒｍａｔｉｏｎａｎｄｄｅｖｅｌｏｐｍｅｎｔｏｆｏｃｅａｎｓｈｅａｒｗａｖｅ￥ＱｉａｏＦａｎｇｌｉａｎｄＺｈａｎｇＱｉｎｇｈｕａ（ＦｉｒｓｔＩｎｓｔｉｔｕｔｅｏｆＯｃｅａｎｏｇｒａｐｈｙ，ｓｔａｔｅｏｃｅａｎｉｃＡｄｍｉｎｉｓ...     

Biological versus physical processes as drivers of large oscillations of the air–sea CO2 flux in the Antarctic marginal ice zone during summer

The fugacity of CO₂ and abundance of chlorophyll a (Chla) were determined in two long transects from the Polar Front to the Antarctic Continent in austral summer, December 1995–January 1996. Large undersaturations of CO₂ in the surface water were observed coinciding with high Chla content. In the major hydrographic regions the mean air–sea fluxes were found to range from −3 to +7 mmol m⁻² d⁻¹ making these regions act as a sink as well as a source for CO₂. In the total 40-d period, the summation of the several strong source and sink regions revealed an overall modest net source of 0.3 mmol m⁻² d⁻¹, this based on the Wanninkhof (J. Geophys. Res. 97 (1992) 7373) quadratic relationship at in situ windspeed. A simple budget approach was used to quantify the role of phytoplankton blooms in the inorganic carbonate system of the Antarctic seas in a time frame spanning several weeks. The major controlling physical factors such as air–sea flux, Ekman pumping and upwelling are included. Net community production varies between −9 and +7 mmol m⁻² d⁻¹, because of the large oscillations in the dominance of autotrophic (CO₂ fixation) versus heterotrophic (CO₂ respiration) activity. Here the mixed layer depth is the major controlling factor. When integrated over time the gross influx and efflux of CO₂ from air to sea is large, but the net residual air/sea exchange is a modest efflux from sea to atmosphere.