A Survey of Ocean Simulation and Rendering Techniques in Computer Graphics

This paper presents a survey of ocean simulation and rendering methods in computer graphics. To model and animate the ocean's surface, these methods mainly rely on two main approaches: on the one hand, those which approximate ocean dynamics with parametric, spectral or hybrid models and use empirical laws from oceanographic research. We will see that this type of methods essentially allows the simulation of ocean scenes in the deep water domain, without breaking waves. On the other hand, physically-based methods use Navier-Stokes Equations (NSE) to represent breaking waves and more generally ocean surface near the shore. We also describe ocean rendering methods in computer graphics, with a special interest in the simulation of phenomena such as foam and spray, and light's interaction with the ocean surface

Geothermal Heating and its Influence on the Meridional Overturning Circulation

The effect of geothermal heating on the meridional overturning circulation is examined using an idealized, coarse-resolution ocean general circulation model. This heating is parameterized as a spatially uniform heat flux of 50 mW m-2 through the (flat) ocean floor, in contrast with previous studies that have considered an isolated hotspot or a series of plumes along the mid-Atlantic ridge. The equilibrated response is largely advective: a deep perturbation of the meridional overturning cell on the order of several Sv is produced, connecting with an upper-level circulation at high latitudes, allowing the additional heat to be released to the atmosphere. Risingmotion in the perturbation deep cell is concentrated near the equator. The upward penetration of this cell is limited by the thermocline, analogous to the role of the stratosphere in limiting the upward penetration of convective plumes in the atmosphere. The magnitude of the advective response is inversely proportional to the deep stratification; with a weaker background meridional overturning circulation and a less stratified abyss, the overturning maximum of the perturbation deep cell is increased. This advective response also cools the low-latitude thermocline. The qualitative behavior is similar in both a single hemisphere and double hemisphere configuration.The anomalous circulation driven by geothermal fluxes is more substantial than previously thought. We are able to understand the structure and strength of the response in the idealized geometry and further extend these ideas to explain the results of Adcroft et al. [2001], where the impact of geothermal heating was examined using a global configuration

Fast and slow Kelvin waves in the Madden-Julian Oscillation of a GCM

The structure of the Madden-Julian Oscillation (MJO) in an 1800-day integration of the Hadley Centre Unified Model was analysed, and interpreted within a Kelvin wave framework. The model was forced with constant equinoctial (March) boundary conditions so that a ``clean'' MJO signal could be separated from the effects of the seasonal cycle and forced interannual variability. The simulated MJO was fairly realistic in terms of its large-scale spatial structure and propagation characteristics, although its period of 30 days (corresponding to an average phase speed of 15 \mps) was shorter than that observed. The signal in deep convection was less coherent than in observations, and appeared to move eastward as a sequence of discrete convective anomalies, rather than by a smooth eastward propagation. Both ``fast'' and ``slow'' equatorial Kelvin waves appeared to play an important role in the eastward propagation of the simulated MJO. Enhanced convection over the Indian Ocean was associated with a ``fast'' equatorial Kelvin wave that propagated eastward at 55 m s-1 over the Pacific. On reaching the west coast of South America, a component of this Kelvin wave propagated northward and southward as a trapped wave along the mountain ranges of Central America and the Andes, in agreement with observations. The anomalous surface easterlies over the tropical eastern Pacific associated with this fast Kelvin wave enhanced the climatological mean easterlies and led to positive convective anomalies over the eastern Pacific consistent with the WISHE mechanism. However, WISHE was not able to account for the eastward development of the convective anomalies over the Indian Ocean/western Pacific region. By splitting the equatorial divergence anomalies of the simulated MJO into their du/dx and dv/dy components, the role of Kelvin wave dynamics in the ``slow'' (15 m s-1) average eastward propagation of the simulated MJO was examined. Although the two components were of comparable magnitude, the \dudx\ component exhibited a pronounced eastward propagation which tended to be disrupted by the \dvdy\ component, thus supporting the paradigm of an underlying, but strongly modified, Kelvin wave mechanism

The last glacial cycle: transient simulations with an AOGCM

A number of transient climate runs simulating the last 120kyr have been carried out using FAMOUS, a fast atmosphere-ocean general circulation model (AOGCM). This is the first time such experiments have been done with a full AOGCM, providing a three-dimensional simulation of both atmosphere and ocean over this period. Our simulation thus includes internally generated temporal variability over periods from days to millennia, and physical, detailed representations of important processes such as clouds and precipitation. Although the model is fast, computational restrictions mean that the rate of change of the forcings has been increased by a factor of 10, making each experiment 12kyr long. Atmospheric greenhouse gases (GHGs), northern hemisphere ice sheets and variations in solar radiation arising from changes in the Earth's orbit are treated as forcing factors, and are applied either separately or combined in different experiments. The long-term temperature changes on Antarctica match well with reconstructions derived from ice-core data, as does variability on timescales longer than 10 kyr. Last Glacial Maximum (LGM) cooling on Greenland is reasonably well simulated, although our simulations, which lack ice-sheet meltwater forcing, do not reproduce the abrupt, millennial scale climate shifts seen in northern hemisphere climate proxies or their slower southern hemisphere counterparts. The spatial pattern of sea surface cooling at the LGM matches proxy reconstructions reasonably well. There is significant anti-correlated variability in the strengths of the Atlantic Meridional Overturning Circulation (AMOC) and the Antarctic Circumpolar Current (ACC) on timescales greater than 10kyr in our experiments. We find that GHG forcing weakens the AMOC and strengthens the ACC, whilst the presence of northern hemisphere ice-sheets strengthens the AMOC and weakens the ACC. The structure of the AMOC at the LGM is found to be sensitive to the details of the ice-sheet reconstruction used. The precessional component of the orbital forcing induces ~20kyr oscillations in the AMOC and ACC, whose amplitude is mediated by changes in the eccentricity of the Earth's orbit. These forcing influences combine, to first order, in a linear fashion to produce the mean climate and ocean variability seen in the run with all forcings

Mechanisms affecting the overturning response in global warming simulations

Climate models used to produce global warming scenarios exhibit widely diverging responses of the thermohaline circulation (THC). To investigate the mechanisms responsible for this variability, a regional Atlantic Ocean model driven with forcing diagnosed from two coupled greenhouse gas simulations has been employed. One of the coupled models (MPI) shows an almost constant THC, the other (GFDL) shows a declining THC in the twenty-first century. The THC evolution in the regional model corresponds rather closely to that of the respective coupled simulation, that is, it remains constant when driven with the forcing from the MPI model, and declines when driven with the GFDL forcing. These findings indicate that a detailed representation of ocean processes in the region covered by the Atlantic model may not be critical for the simulation of the overall THC changes in a global warming scenario, and specifically that the coupled model’s rather coarse representation of water mass formation processes in the subpolar North Atlantic is unlikely to be the primary cause for the large differences in the THC evolution. Sensitivity experiments have confirmed that a main parameter governing the THC response to global warming is the density of the intermediate waters in the Greenland–Iceland–Norwegian Seas, which in turn influences the density of the North Atlantic Deep Water, whereas changes in the air–sea heat and freshwater fluxes over the subpolar North Atlantic are only of moderate importance, and mainly influence the interannual–decadal variability of THC. Finally, as a consequence of changing surface fluxes, the Labrador Sea convection ceases by about 2030 under both forcings (i.e., even in a situation where the overall THC is stable) indicating that the eventual breakdown of the convection is likely but need not coincide with substantial THC changes

Dynamical ocean forcing of the Madden-Julian Oscillation at lead times of up to five months

We show that a simple three-dimensional ocean model linearised about a resting basic state can accurately simulate the dynamical ocean response to wind forcing by the Madden-Julian Oscillation (MJO). This includes the propagation of equatorial waves in the Indian Ocean, from the generation of oceanic equatorial Kelvin waves to the arrival of downwelling oceanic equatorial Rossby waves in the western Indian Ocean, where they have been shown to trigger MJO convective activity. Simulations with idealised wind forcing suggest that the latitudinal width of this forcing plays a crucial role in determining the potential for such feedbacks. Forcing the model with composite MJO winds accurately captures the global ocean response, demonstrating that the observed ocean dynamical response to the MJO can be interpreted as a linear response to surface wind forcing. The model is then applied to study “primary” Madden-Julian events, which are not immediately preceded by any MJO activity nor by any apparent atmospheric triggers, but have been shown to coincide with the arrival of downwelling oceanic equatorial Rossby waves. Case study simulations show how this oceanic equatorial Rossby wave activity is partly forced by reflection of an oceanic equatorial Kelvin wave triggered by a westerly wind burst 140 days previously, and partly directly forced by easterly wind stress anomalies around 40 days prior to the event. This suggests predictability for primary Madden-Julian events on times scales of up to five months, following the re-emergence of oceanic anomalies forced by winds almost half a year earlier

Mediterranean Sea response to climate change in an ensemble of twenty first century scenarios

The Mediterranean climate is expected to become warmer and drier during the twenty-first century. Mediterranean Sea response to climate change could be modulated by the choice of the socio-economic scenario as well as the choice of the boundary conditions mainly the Atlantic hydrography, the river runoff and the atmospheric fluxes. To assess and quantify the sensitivity of the Mediterranean Sea to the twenty-first century climate change, a set of numerical experiments was carried out with the regional ocean model NEMOMED8 set up for the Mediterranean Sea. The model is forced by air–sea fluxes derived from the regional climate model ARPEGE-Climate at a 50-km horizontal resolution. Historical simulations representing the climate of the period 1961–2000 were run to obtain a reference state. From this baseline, various sensitivity experiments were performed for the period 2001–2099, following different socio-economic scenarios based on the Special Report on Emissions Scenarios. For the A2 scenario, the main three boundary forcings (river runoff, near-Atlantic water hydrography and air–sea fluxes) were changed one by one to better identify the role of each forcing in the way the ocean responds to climate change. In two additional simulations (A1B, B1), the scenario is changed, allowing to quantify the socio-economic uncertainty. Our 6-member scenario simulations display a warming and saltening of the Mediterranean. For the 2070–2099 period compared to 1961–1990, the sea surface temperature anomalies range from +1.73 to +2.97 °C and the SSS anomalies spread from +0.48 to +0.89. In most of the cases, we found that the future Mediterranean thermohaline circulation (MTHC) tends to reach a situation similar to the eastern Mediterranean Transient. However, this response is varying depending on the chosen boundary conditions and socio-economic scenarios. Our numerical experiments suggest that the choice of the near-Atlantic surface water evolution, which is very uncertain in General Circulation Models, has the largest impact on the evolution of the Mediterranean water masses, followed by the choice of the socio-economic scenario. The choice of river runoff and atmospheric forcing both have a smaller impact. The state of the MTHC during the historical period is found to have a large influence on the transfer of surface anomalies toward depth. Besides, subsurface currents are substantially modified in the Ionian Sea and the Balearic region. Finally, the response of thermosteric sea level ranges from +34 to +49 cm (2070–2099 vs. 1961–1990), mainly depending on the Atlantic forcing