Slab ocean model
Contents
General description of the slab ocean model
The slab ocean model is primarily based on the work of Codron (2012). This model uses the same horizontal grid as the GCM and is composed of two layers. The first layer (50 m depth) represents the surface mixed layer, where the exchanges with the atmosphere take place. The second layer (150 m depth) represents the deep ocean. Each oceanic grid point of the model is defined by the sea ice fraction & thickness and snow thickness.
In general for a slab ocean, the temperature at the ocean surface evolves as per the following equation: $$\rho c_p H \frac{\partial T}{\partial t} = F_{surf} + Q_{flux}$$
The evolution of temperature is driven by surface fluxes ($$F_{surf} = F_{rad} + F_{sens} + F_{lat}$$) arising from the atmosphere model and a prescribed / parameterised flux ($$Q_{flux}$$) that can represent the effect of geothermal heating / ocean circulation (depends on the kind of planet that is being simulated). Both, the surface fluxes and the conductive fluxes coming from within interact with snow or ice that can be present at the ocean surface. Over each grid point, the sea ice is represented by a uniform layer of depth $$H$$ and fractional area $$f$$. It may be covered by a layer of snow, which has its own heat capacity. The temperature at the bottom of the ice layer is always equal to the freezing temperature of sea water $$T_o$$.
In its current form in the Generic-PCM, the slab ocean model handles two key processes, (a) Sea ice formation, and (b) Heat transport by ocean circulation:
Sea ice formation
Sea ice forms when the ocean temperature falls below the freezing point of water ($$T_0=–1.8°C$$, the freezing point of Earth's sea water; function of salinity) and melts when its temperature rises above freezing. The changes in ice extent and thickness are driven by energy conservation in such a way that the ocean temperature is kept at $$T_0$$ as long as ice is present. Specifically, if the temperature of the surface layer instantaneously falls below its freezing point, its temperature is set back to $$T_0$$, and the resulting energy difference is used to create ice. A framework of rules determines how energy differences can make ice thickness and coverage evolve. Conversely, if the ocean temperature becomes positive when ice is already present, its temperature is set back to $$T_0$$, and the resulting energy difference is used to melt ice.
The same arguments as previously stated contribute to the melting of snow if it has been precipitated onto the ice. Further, the surface albedo can be that of snow (if any) or that of bare ice:
$$A = A^\rm{max}_\rm{ice} – (A^\rm{max}_\rm{ice} – A^\rm{min}_\rm{ice} ) e^{\left ( \frac{–h_\rm{ice}}{h^0_\rm{ice}} \right )}$$
Here, $$A$$ is the albedo, $$A^\rm{max}_\rm{ice}$$ = 0.65 is the maximal albedo, $$A^\rm{min}_\rm{ice}$$ = 0.2 is the minimal albedo, $$h_\rm{ice}$$ is the ice thickness (in m) and $$h^0_\rm{ice}$$ = 0.5 m. The albedo over ice-free ocean is assumed to be 0.07. However, by definition, the Generic-PCM is meant to handle a more generalised case of planets. Since the albedo of a material (for e.g., ice or snow) is spectrally-dependent, the albedo of ice will be different around a star that is unlike the Sun. For instance, the ice/snow albedo climate feedback is significantly weaker for planets orbiting M-dwarf stars (for e.g., Proxima Centauri b) than for planets orbiting G-type stars like the Sun (see Joshi & Haberle [2012] and Shields [2014]). The Generic-PCM accounts for this by instituting a spectral dependence of albedo.
It should be noted that changes in the sea-ice model, or in parameters such as the snow albedo, can have large impacts on the simulated climate. This happens through changes in the ice extent and the global temperature through the ice-albedo feedback. This effect is, however, largely independent of the heat transport schemes presented in this page.
Heat transport by ocean circulation
The transport of heat by the ocean circulation is given by four components:
- First, the impact of subgrid-scale eddies is represented by horizontal diffusion, with a uniform diffusivity in both layers.
- Then, the mean wind-driven circulation is computed by calculating the Ekman mass fluxes in the surface layer from the surface wind stress, and taking an opposite return flow at depth. These mass fluxes are then used to advect the ocean temperature horizontally. In the case of divergent horizontal mass fluxes, the upwelling or downwelling mass flux is obtained by continuity. This simplified model reproduces the global meridional oceanic heat transport quite closely compared to an Oceanic-GCM.
Other ocean circulation features like density-driven / thermohaline circulation and horizontal gyres are not present in the model. While these features play an important role in shaping regional climates on modern Earth, their effect is weaker on the global average; gyres would be absent in the case of a global ocean.
The 2024 version of the slab ocean model
Compared to the old model, the new model has the following changes (non-exhaustive). Most of these have already been explained above:
- It can be used in parallel mode.
- It has a more realistic description of sea ice creation and evolution - simultaneous surface, side and bottom melting / freezing depending on fluxes.
- Snow has an effective heat capacity.
- Snow has "weight"; it can sink an ice block if there is too much of it.
- Snow can be blown off by wind.
- The two-layer ocean facilitates convective adjustment.
- Diffusion can follow the Gent-McWilliams scheme + Eddy diffusivity.
- Since Ekman Transport doesn't work at the equator (Coriolis term = 0), we prescribe an alternative "Sverdrup scheme" for regions close to the equator which relies on direct friction-driven wind transport.
Variables of the slab ocean model which can be written as outputs in diagfi.nc
- tslab1, tslab2: temperature of the surface and deep ocean layers respectively (K)
- pctsrf_sic: grid fraction of sea ice
- sea_ice: mass of sea ice (kg/m$$^2$$)
- tsea_ice: temperature of the sea ice surface (K)
- dt_hdiff1, dt_hdiff2: heating of the surface and deep ocean layers by horizontal diffusion (W/m$$^2$$)
- dt_ekman1, dt_ekman2: heating of the surface and deep ocean layers by Ekman transport (W/m$$^2$$)
- dt_gm1, dt_gm2: heating of the surface and deep ocean layers by Gent-McWilliams eddy transport (W/m$$^2$$)
Technical aspects
The slab ocean model is activated in callphys.def with: ok_slab_ocean = .true.
The sea ice model is activated in callphys.def with: ok_slab_sic = .true.
The horizontal heat transport is activated in callphys.def with: ok_slab_heat_transp = .true. The horizontal heat transport cannot be activated in parallel at the moment.
Values of A$$^{max}_{ice}$$ and A$$^{min}_{ice}$$ are fixed in slab_ice_h.F90
rnat is a standard diagfi output, which either represents 0 for oceanic grid points or 1 for continent grid points.
The model doesn't permit surface fluxes to enter the deep ocean. The latter is present in the model for exchanges and transport to take place.
The model also does not compute dynamical sea ice drift.
QUESTION FOR FRANCIS: It seems that the Sverdrup scheme might only be beneficial for the special case of the Earth? Need to confirm.
The model follows energy conservation. However, if an upper limit on sea ice thickness is prescribed (for e.g., 10 m, for faster convergence), energy conservation can be violated. If one wants to purely conserve energy (at the cost of computation time), no limit on sea ice thickness can be prescribed. It is OK to prescribe an upper limit for the case of modern Earth (however, if your prescription is too small, it can lead to anomalous behaviour like stronger seasonal cycles). However, it might not be accurate to do this for more "exotic" cases like Snowball Earth and tidally locked ocean worlds.
