Coupled climate modelling of ocean circulation changes during ice age inception

Freshening of high latitude surface waters can change the large-scale oceanic transport of heat and salt. Consequently, atmospheric and sea ice perturbations over the deep water production sites excite a large-scale response establishing an oceanic "teleconnection" with time scales of years to centuries. To study these feedbacks, a coupled atmosphere-ocean-sea ice model consisting of a two dimensional atmospheric energy and moisture balance model (EMBM) coupled to a thermodynamic sea ice model and an ocean general circulation model is utilised. The coupled model reproduces many aspects of the present oceanic circulation. We also investigate the climate impact of changes in fresh water balance during an ice age initiation. In this experiment part of the precipitation over continents is stored within continental ice sheets. During the buildup of ice sheets the oceanic stratification in the North Atlantic is weakened by a reduced continental run-off leading to an enhanced thermohaline circulation. Under these conditions salinity is redistributed such that deep water is more saline than under present conditions. Once the ice sheets built up, we simulate an ice age climate without net fresh water storage on the continents. In this case the coupled model reproduces the shallow and weak overturning cell, an ice edge advance insulating the upper ocean, and many other aspects of the glacial circulation.     

Coupled ice‐ocean variability in the Greenland Sea

Abstract

A major surface feature of the Greenland Sea during winter is the frequent eastward extension of sea ice south of 75°N and an associated embayment to the north. These features are nominally connected with the East Greenland Current, and both the promontory and the embayment are readily apparent on climatic ice charts. However, there are significant changes in these features on time‐scales as short as a few days. Using a combination of satellite microwave images (SSM/I) of ice cover, meteorological data and in situ velocity, temperature and salinity records, we relate the ice distribution and its changes to the developing structure and circulation of the upper ocean during winter 1988–1989. Our measurements illustrate the preconditioning that leads to convective overturn, which in turn brings warmer water to the surface and results in the rapid disappearance of ice. In particular, the surface was cooled to the freezing point by early December and the salinity then increased through ice formation (about 0.016 m d^‐1) and brine rejection. Once the vertical density gradient was sufficiently eroded, a period of high heat flux (>300 W m^‐2) in late January provided enough buoyancy loss to convectively mix the upper water column to at least 200 m. We estimate vertical velocities at about 3 cm s^‐1 downward during the initial sinking. The deepening of the thermocline raised surface temperatures by over 1°C resulting in nearly 1.5 × 10⁵ km² of ice‐melt within two days. Average rates of ice retreat are about 11 km d^‐1 southwestward, generally consistent with a wind‐driven flow. Comparison of hydrographic surveys from before and after the overturning indicate the fresh water was advected out of the area, possibly to the south and east of our moorings.     

Application of a coupled ice‐ocean model to the Labrador Sea

Abstract

The steady, coupled ice‐ocean circulation model of Willmott and Mysak (1989) for a meridional channel is applied to the Labrador Sea for the winter season. The model consists of a thermodynamic reduced‐gravity ocean combined with a variable thickness ice cover that is in thermal equilibrium. Upon specifying the forcing fields of surface air temperature, wind stress and water temperature along the open southern boundary, the winter climatological ice‐edge position, ice thickness, ocean circulation and temperature fields are determined in the channel domain. The sensitivity of the results to the various model parameters is examined. In particular, the optimum heat exchange coefficients for the interfaces of air‐water, ice‐water and air‐ice are found.

The model ice‐edge position compares favourably with the 50% winter climatological ice concentration isoline obtained from an analysis of 32 years (1953–84) of sea‐ice concentration data. The simulations of the ocean temperature and ice thickness are also quite realistic according to the observed records available. The model is also applied to two specific winters (1981 and 1983) during which anomalous sea‐ice and weather conditions prevailed in the Labrador Sea.     

Near real-time ocean circulation assimilation and prediction in the Intra-Americas Sea with ROMS

We present the feasibility of a prototype, near real-time assimilation and ensemble prediction system for the Intra-Americas Sea run autonomously aboard a ship of opportunity based on the Regional Ocean Modeling System (ROMS). Predicting an ocean state depends upon numerical models that contain uncertainties in their modeled physics, initial conditions, and model state. An advanced model, four-dimensional variational assimilation, and ensemble forecasting techniques are used to account for each of these uncertainties. Every 3 days, data from the previous 7 days period were assimilated to generate an estimate of the circulation and to create an ensemble of 2 weeks forecasts of the ocean state. This paper presents the methods and results for a multi-resolution assimilation system and ensemble forecasts of surface fields and dominant surface circulation features. When compared to post-processed science quality observations, the state estimates suffer from our reliance on real-time, quick-look satellite observations of the ocean surface. Despite a number of issues, the ensemble forecast estimate is often superior to observational persistence. This proof-of-concept prototype performed well enough to reveal deficiencies, provide useful insights, valuable lessons, and guidance for future improvements in real-time ocean forecasting.     

Near real-time ocean circulation assimilation and prediction in the Intra-Americas Sea with ROMS

We present the feasibility of a prototype, near real-time assimilation and ensemble prediction system for the Intra-Americas Sea run autonomously aboard a ship of opportunity based on the Regional Ocean Modeling System (ROMS). Predicting an ocean state depends upon numerical models that contain uncertainties in their modeled physics, initial conditions, and model state. An advanced model, four-dimensional variational assimilation, and ensemble forecasting techniques are used to account for each of these uncertainties. Every 3 days, data from the previous 7 days period were assimilated to generate an estimate of the circulation and to create an ensemble of 2 weeks forecasts of the ocean state. This paper presents the methods and results for a multi-resolution assimilation system and ensemble forecasts of surface fields and dominant surface circulation features. When compared to post-processed science quality observations, the state estimates suffer from our reliance on real-time, quick-look satellite observations of the ocean surface. Despite a number of issues, the ensemble forecast estimate is often superior to observational persistence. This proof-of-concept prototype performed well enough to reveal deficiencies, provide useful insights, valuable lessons, and guidance for future improvements in real-time ocean forecasting.     

Sensitivity of Hudson Bay Sea ice and ocean climate to atmospheric temperature forcing

A regional sea-ice?Cocean model was used to investigate the response of sea ice and oceanic heat storage in the Hudson Bay system to a climate-warming scenario. Projections of air temperature (for the years 2041?C2070; effective CO₂ concentration of 707?C950?ppmv) obtained from the Canadian Regional Climate Model (CRCM 4.2.3), driven by the third-generation coupled global climate model (CGCM 3) for lateral atmospheric and land and ocean surface boundaries, were used to drive a single sensitivity experiment with the delta-change approach. The projected change in air temperature varies from 0.8°C (summer) to 10°C (winter), with a mean warming of 3.9°C. The hydrologic forcing in the warmer climate scenario was identical to the one used for the present climate simulation. Under this warmer climate scenario, the sea-ice season is reduced by 7?C9?weeks. The highest change in summer sea-surface temperature, up to 5°C, is found in southeastern Hudson Bay, along the Nunavik coast and in James Bay. In central Hudson Bay, sea-surface temperature increases by over 3°C. Analysis of the heat content stored in the water column revealed an accumulation of additional heat, exceeding 3?MJ?m^?3, trapped along the eastern shore of James and Hudson bays during winter. Despite the stratification due to meltwater and river runoff during summer, the shallow coastal regions demonstrate a higher capacity of heat storage. The maximum volume of dense water produced at the end of winter was halved under the climate-warming perturbation. The maximum volume of sea ice is reduced by 31% (592?km3) while the difference in the maximum cover is only 2.6% (32,350?km²). Overall, the depletion of sea-ice thickness in Hudson Bay follows a southeast?Cnorthwest gradient. Sea-ice thickness in Hudson Strait and Ungava Bay is 50% thinner than in present climate conditions during wintertime. The model indicates that the greatest changes in both sea-ice climate and heat content would occur in southeastern Hudson Bay, James Bay, and Hudson Strait.     

Modelling the Baltic Sea ocean climate on centennial time scale: temperature and sea ice

This study considers the possible use of different kinds of forcing datasets in Baltic Sea ocean climate modelling on centennial time scales, in particular for the past half millennium. We demonstrate that high-quality station data of the past century and gridded multi-proxy reconstructions for the past 500 years can be used with great success but with various levels of detail. We also demonstrate that output data from the state-of-the-art global climate model EcHo-G are not suitable for modelling the Baltic Sea ocean climate. Two climate properties were studied: the annual maximum ice extent (MIB) and the vertically and horizontally integrated annual water temperature. Centennial time scale results indicate that the seventeenth and nineteenth centuries were the coldest centuries, while the 1690s were the coldest decade and 1695 the coldest year in the last 500 years. The results also indicate that the twentieth century was the warmest century with the least MIB of the last 500 years. On a decadal time scale, the 1990s, 1930s and 1730s were the warmest decades and comparable in terms of both water temperature and MIB. The year 1989 had the minimum observed MIB of only 52,000 km², implying that the Baltic Sea has been partly ice covered in all winters of the past half millennium. Even though different climate forcing mechanisms may operate on the climate system today compared to over the last half millennium, this study cannot clearly state that the region is experiencing climate change outside the natural limits of the past 500 years.     

Effects of tropical cyclones on large-scale circulation and ocean heat transport in the South China Sea

In this study, we investigate the influence of tropical cyclones (TCs) on large-scale circulation and ocean heat transport in the South China Sea (SCS) by using an ocean general circulation model at a 1/8° resolution during 2000–2008. The model uses a data assimilation system to assimilate observations in order to improve the representation of SCS circulation. The results reveal an unexpected deep SCS circulation anomaly induced by TCs, which suggests that effects of TC can penetrate deeper into the ocean. This deep effect may result from the near inertial oscillations excited by TCs. The inertial oscillations can propagate downward to the oceanic interior. The analyses confirm that TCs have two effects on ocean heat transport of the SCS. Firstly, the wind stress curl induced by TCs affects the structure of SCS circulation, and then changes heat transport. Secondly, TCs pump surface heat downward to the thermocline, increasing the heat injection from the atmosphere to the ocean. Two effects together amplify the outflow of the surface heat southward away the SCS through the Mindoro and Karimata Straits. The TC-induced heat transports through the Mindoro, Balabac and Karimata Straits account for 20 % of the total heat transport through three straits. An implication of this study is that ocean models need to simulate the TC effect on heat transport in order to correctly evaluate the role of the SCS through flow in regulating upper ocean circulation and climate in the Indonesian maritime continent and its adjacent regions.     

Influence of coupling on atmosphere, sea ice and ocean regional models in the Ross Sea sector, Antarctica

Air–sea ice–ocean interactions in the Ross Sea sector form dense waters that feed the global thermohaline circulation. In this paper, we develop the new limited-area ocean–sea ice–atmosphere coupled model TANGO to simulate the Ross Sea sector. TANGO is built up by coupling the atmospheric limited-area model MAR to a regional configuration of the ocean–sea ice model NEMO. A method is then developed to identify the mechanisms by which local coupling affects the simulations. TANGO is shown to simulate realistic sea ice properties and atmospheric surface temperatures. These skills are mostly related to the skills of the stand alone atmospheric and oceanic models used to build TANGO. Nonetheless, air temperatures over ocean and winter sea ice thickness are found to be slightly improved in coupled simulations as compared to standard stand alone ones. Local atmosphere ocean feedbacks over the open ocean are found to significantly influence ocean temperature and salinity. In a stand alone ocean configuration, the dry and cold air produces an ocean cooling through sensible and latent heat loss. In a coupled configuration, the atmosphere is in turn moistened and warmed by the ocean; sensible and latent heat loss is therefore reduced as compared to the stand alone simulations. The atmosphere is found to be less sensitive to local feedbacks than the ocean. Effects of local feedbacks are increased in the coastal area because of the presence of sea ice. It is suggested that slow heat conduction within sea ice could amplify the feedbacks. These local feedbacks result in less sea ice production in polynyas in coupled mode, with a subsequent reduction in deep water formation.     

Wind-driven changes in Southern Ocean residual circulation, ocean carbon reservoirs and atmospheric CO2

The effect of idealized wind-driven circulation changes in the Southern Ocean on atmospheric CO₂ and the ocean carbon inventory is investigated using a suite of coarse-resolution, global coupled ocean circulation and biogeochemistry experiments with parameterized eddy activity and only modest changes in surface buoyancy forcing, each experiment integrated for 5,000 years. A positive correlation is obtained between the meridional overturning or residual circulation in the Southern Ocean and atmospheric CO₂: stronger or northward-shifted westerly winds in the Southern Hemisphere result in increased residual circulation, greater upwelling of carbon-rich deep waters and oceanic outgassing, which increases atmospheric pCO₂ by ～20 μatm; weaker or southward-shifted winds lead to the opposing result. The ocean carbon inventory in our model varies through contrasting changes in the saturated, disequilibrium and biogenic (soft-tissue and carbonate) reservoirs, each varying by O(10–100) PgC, all of which contribute to the net anomaly in atmospheric CO₂. Increased residual overturning deepens the global pycnocline, warming the upper ocean and decreasing the saturated carbon reservoir. Increased upwelling of carbon- and nutrient-rich deep waters and inefficient biological activity results in subduction of unutilized nutrients into the ocean interior, decreasing the biogenic carbon reservoir of intermediate and mode waters ventilating the Northern Hemisphere, and making the disequilibrium carbon reservoir more positive in the mode waters due to the reduced residence time at the surface. Wind-induced changes in the model carbon inventory are dominated by the response of the global pycnocline, although there is an additional abyssal response when the peak westerly winds change their latitude, altering their proximity to Drake Passage and changing the depth extent of the southward return flow of the overturning: a northward shift of the westerly winds isolates dense isopycnals, allowing biogenic carbon to accumulate in the deep ocean of the Southern Hemisphere, while a southward shift shoals dense isopycnals that outcrop in the Southern Ocean and reduces the biogenic carbon store in the deep ocean.     

Transient climate change scenario simulation of the Mediterranean Sea for the twenty-first century using a high-resolution ocean circulation model

A scenario of the Mediterranean Sea is performed for the twenty-first century based on an ocean modelling approach. A climate change IPCC-A2 scenario run with an atmosphere regional climate model is used to force a Mediterranean Sea high-resolution ocean model over the 1960–2099 period. For comparison, a control simulation as long as the scenario has also been carried out under present climate fluxes. This control run shows air–sea fluxes in agreement with observations, stable temperature and salinity characteristics and a realistic thermohaline circulation simulating the different intermediate and deep water masses described in the literature. During the scenario, warming and saltening are simulated for the surface (+3.1°C and + 0.48 psu for the Mediterranean Sea at the end of the twenty-first century) and for the deeper layers (+1.5°C and + 0.23 psu on average). These simulated trends are in agreement with observed trends for the Mediterranean Sea over the last decades. In addition, the Mediterranean thermohaline circulation (MTHC) is strongly weakened at the end of the twenty-first century. This behaviour is mainly due to the decrease in surface density and so the decrease in winter deep-water formation. At the end of the twenty-first century, the MTHC weakening can be evaluated as −40% for the intermediate waters and −80% for the deep circulation with respect to present-climate conditions. The characteristics of the Mediterranean Outflow Waters flowing into the Atlantic Ocean are also strongly influenced during the scenario.     

Climate scenarios of sea level rise for the northeast Atlantic Ocean: a study including the effects of ocean dynamics and gravity changes induced by ice melt

Here we present a set of regional climate scenarios of sea level rise for the northeast Atlantic Ocean. In this study, the latest observations and results obtained with state-of-the-art climate models are combined. In addition, regional effects due to ocean dynamics and changes in the Earth’s gravity field induced by melting of land-based ice masses have been taken into account. The climate scenarios are constructed for the target years 2050 and 2100, for both a moderate and a large rise in global mean atmospheric temperature (2 °C and 4 °C in 2100 respectively). The climate scenarios contain contributions from changes in ocean density (global thermal expansion and local steric changes related to changing ocean dynamics) and changes in ocean mass (melting of mountain glaciers and ice caps, changes in the Greenland and Antarctic ice sheets, and (minor) terrestrial water-storage contributions). All major components depend on the global temperature rise achieved in the target periods considered. The resulting set of climate scenarios represents our best estimate of twenty-first century sea level rise in the northeast Atlantic Ocean, given the current understanding of the various contributions. For 2100, they yield a local rise of 30 to 55 cm and 40 to 80 cm for the moderate and large rise in global mean atmospheric temperature, respectively.     

Breeze circulation in the Black Sea region

The modeling of atmospheric circulation in the Black Sea region using the WRF-ARW model for 30 summer days is considered to single out the contribution of the breeze component of wind speed. Assuming the quasiperiodic nature of breeze, the speed of daytime and nighttime breeze is computed for the whole region and for five separate areas on the coast. Specific features of the formation of daytime and nighttime breeze in these areas are distinguished. It is revealed that the presence of coastal mountains and the complex coastline determine the typical features of breeze development.     

A link between winter precipitation in Europe and the Arctic Sea ice,sea surface temperature,and atmospheric circulation

The spatial and temporal pattern of the link between the winter precipitation variability and variations in the North Atlantic sea surface temperature, the Arctic sea ice concentration, and 500 hPa geopotential height in the Northern Hemisphere is analyzed for the period of 1952-2012. The analysis reveals two principal modes of covariability in the analyzed characteristics. The first mode which explains the most part of covariability, is related to the impact of the North Atlantic Oscillation. The second mode indicates the significant contribution of the Atlantic Multidecadal Oscillation associated with winter precipitation anomalies of the same sign in Europe with the maxima on the East European Plain and in the Balkan region during the positive phase of AMO.     

Seasonal surface ocean circulation and dynamics in the Philippine Archipelago region during 2004–2008

The dynamics of the seasonal surface circulation in the Philippine Archipelago (117°E–128°E, 0°N–14°N) are investigated using a high-resolution configuration of the Regional Ocean Modeling System (ROMS) for the period of January 2004–March 2008. Three experiments were performed to estimate the relative importance of local, remote and tidal forcing. On the annual mean, the circulation in the Sulu Sea shows inflow from the South China Sea at the Mindoro and Balabac Straits, outflow into the Sulawesi Sea at the Sibutu Passage, and cyclonic circulation in the southern basin. A strong jet with a maximum speed exceeding 100 cm s⁻¹ forms in the northeast Sulu Sea where currents from the Mindoro and Tablas Straits converge. Within the Archipelago, strong westward currents in the Bohol Sea carry the surface water of the western Pacific (WP) from the Surigao Strait into the Sulu Sea via the Dipolog Strait. In the Sibuyan Sea, currents flow westward, which carry the surface water from the WP near the San Bernardino Strait into the Sulu Sea via the Tablas Strait.These surface currents exhibit strong variations or reversals from winter to summer. The cyclonic (anticyclonic) circulation during winter (summer) in the Sulu Sea and seasonally reversing currents within the Archipelago region during the peak of the winter (summer) monsoon result mainly from local wind forcing, while remote forcing dominates the current variations at the Mindoro Strait, western Sulu Sea and Sibutu passage before the monsoons reach their peaks. The temporal variations (with the mean removed), also referred to as anomalies, of volume transports in the upper 40 m at eight major Straits are caused predominantly by remote forcing, although local forcing can be large during sometime of a year. For example, at the Mindoro Strait, the correlation between the time series of transport anomalies due to total forcing (local, remote and tides) and that due only to the remote forcing is 0.81 above 95% significance, comparing to the correlation of 0.64 between the total and local forcing. Similarly, at the Sibutu Passage, the correlation is 0.96 for total versus remote effects, comparing to 0.53 for total versus local forcing. The standard deviations of transports from the total, remote and local effects are 0.59 Sv, 0.50 Sv, and 0.36 Sv, respectively, at the Mindoro Strait; and 1.21 Sv, 1.13 Sv, and 0.59 Sv at the Sibutu Passage. Nonlinear rectification of tides reduces the mean westward transports at the Surigao, San Bernardino and Dipolog Straits, and it also has non-negligible influence on the seasonal circulation in the Sulu Sea.