The climatic response of the Arctic Ocean to Soviet river diversions

A numerical model is constructed to evaluate the effect of river diversions on the circulation of the Arctic Ocean, including the climatically important response in the extent of sea ice. The ocean model solves the primitive equations of motion in finite-difference form for the irregular geometry of the Arctic Ocean and Greenland/Norwegian Sea, using 110 km horizontal grid spacing and up to 13 unevenly spaced levels in the vertical. Annual mean atmospheric conditions and river discharges are specified from observations. The presence of sea ice is diagnosed on the basis of model ocean temperature; and the effects of sea ice on the surface fluxes of momentum, heat, and salt are included in a simplified way. Lateral exchanges at the southernmost boundary are held near observed values but respond to circulation changes in the Greenland/Norwegian Sea. Three equilibrium solutions are obtained by eighty-year integrations from simple initial conditions: the first with inflow from all rivers, the second with one-third of the inflow diverted from four major rivers (the Ob, Yenesei, Dvina, and Pechora), and the third with total diversion from those rivers. The middle case corresponds to maximal diversions which are either planned or envisioned by the Soviet Union over the next fifty years, whereas the final extreme case is run in the event that model sensitivity is low relative to that of nature.The control integration gives a good simulation of known water masses and currents. In the Central Arctic, for example, the model correctly predicts a strong shallow halocline, a relatively warm intermediate layer of Atlantic origin, and a temperature jump across the deep Lomonosov Ridge. The overall pattern of surface salinity and the margin of the pack ice are also properly simulated.When runoff into the marginal Kara and Barents Seas is diverted, either in part or in full, almost no effect on the halocline results in the Central Arctic. In particular, deep convection does not develop in the Eurasian Basin, the possibility of which was suggested by Aagaard and Coachman (1975). The vertical stability within the two marginal seas is considerably decreased by the total diversion of four rivers, but not to the point of convective overturning. The surface currents in this area change to confine the water with increased salinity to the shelf region. At deeper levels, an increased salinity tongue spreads into the deep basins of the ice-free Greenland/Norwegian Sea, where existing deep convection is slightly enhanced. As a result, there is some additional heat loss from the Atlantic layer before it enters the Central Arctic. The ice extent remains nearly the same as before within the Kara and Barents Seas. In fact, since modified bottom currents over the continental shelf bring in less heat from the Greenland Sea, an increased thickness of sea ice may result there, in spite of reduced vertical stability. These model responses are generally in agreement with those suggested by Micklin (1981) and by Soviet investigations of the effect of river diversions. These annualmean results should be regarded as tentative, pending confirmation by studies which include the seasonal cycles of runoff and atmospheric forcing.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Baltic Sea climate in the late twenty-first century: a dynamical downscaling approach using two global models and two emission scenarios

A regional ocean circulation model was used to project Baltic Sea climate at the end of the twenty-first century. A set of four scenario simulations was performed utilizing two global models and two forcing scenarios. To reduce model biases and to spin up future salinity the so-called Δ-change approach was applied. Using a regional coupled atmosphere–ocean model 30-year climatological monthly mean changes of atmospheric surface data and river discharge into the Baltic Sea were calculated from previously conducted time slice experiments. These changes were added to reconstructed atmospheric surface fields and runoff for the period 1903–1998. The total freshwater supply (runoff and net precipitation) is projected to increase between 0 and 21%. Due to increased westerlies in winter the annual mean wind speed will be between 2 and 13% larger compared to present climate. Both changes will cause a reduction of the average salinity of the Baltic Sea between 8 and 50%. Although salinity in the entire Baltic might be significantly lower at the end of the twenty-first century, deep water ventilation will very likely only slightly change. The largest change is projected for the secondary maximum of sea water age within the halocline. Further, the average temperature will increase between 1.9 and 3.2°C. The temperature response to atmospheric changes lags several months. Future annual maximum sea ice extent will decrease between 46 and 77% in accordance to earlier studies. However, in contrast to earlier results in the warmest scenario simulation one ice-free winter out of 96 seasons was found. Although wind speed changes are uniform, extreme sea levels may increase more than the mean sea level. In two out of four projections significant changes of 100-year surge heights were found.     

SIMULATION ON CLIMATIC VARIATION OF ARCTIC SEA ICE THROUGH AN ICE-OCEAN COUPLED MODEL*

An ocean model developed by the Institute of Marine Research and the University of Bergen in Norway (BOM) and a state-of-the-art sea ice model developed by NCAR (CSIM4) are coupled, Considering influences of 9 major rivers,forced by the NCEP reanalysis atmospheric fields and the Levitus surface salinity,the Arctic sea ice climatic variation from January 1949 to December.1999 was simulated through the coupled model.The comparison of simulated results and observations shows that:(1)the long-term ice concentration variation tendencies are in consistent with the observations in the divisional ocean regions;(2)simulated ice thickness horizontal distribution is reasonable.Simulated ice thickness has a decreasing tendency in the central Arctic,which agrees with the submarine observations.Simulated annually maximum ice thickness is highly related to observed fast-ice thickness off the Russian coast;and (3)sea ice area/volume fluxes through the Fram Strait are in accord with the satellite-derived data.Generally,the coupled model successfully simulated the Arctic Ocean sea ice climatic variation.     

Effects of Spectral Nudging on Oceanic States in a Coarse-Resolution Model

A 1-degree global model is used to investigate the skill of spectral nudging at coarse resolution by performing two numerical experiments, one with spectral nudging and the other without. In the spectral nudging experiment, the model temperature and salinity are nudged to an observed climatological monthly-mean field. The study compares the model mean state, as well as the interannual and decadal variability of oceanic quantities with observations, (e.g., sea surface height (SSH) and sea surface temperature (SST)). Spectral nudging is found to be effective in constraining model drift from the observed mean state of temperature and salinity in the global ocean, which has been reported in previous studies. The present study further shows that spectral nudging significantly improves the model skill of topostrophy (a measure of currents flowing along the topography) in water depth below 2000?m with no clear improvement elsewhere. Despite its known ability to damp oceanic variability at various time scales, spectral nudging can still represent the interannual and decadal variability of SSH and SST well, to a degree comparable to the other experiment.     

A possible feedback mechanism involving the Arctic freshwater, the Arctic sea ice, and the North Atlantic Drift

Model studies point to enhanced warming and to increased freshwater fluxes to high northern latitudes in response to global warming. In order to address possible feedbacks in the ice-ocean system in response to such changes, the combined effect of increased freshwater input to the Arctic Ocean and Arctic warming--the latter manifested as a gradual melting of the Arctic sea ice--is examined using a 3-D isopycnic coordinate ocean general circulation model. A suite of three idealized experiments is carried out: one control integration, one integration with a doubling of the modern Arctic river runoff, and a third more extreme case, where the river runoff is five times the modern value. In the two freshwater cases, the sea ice thickness is reduced by 1.5-2 m in the central Arctic Ocean over a 50-year period. The modelled ocean response is qualitatively the same for both perturbation experiments: freshwater propagates into the Atlantic Ocean and the Nordic Seas, leading to an initial weakening of the North Atlantic Drift.Furthermore, changes in the geostrophic currents in the central Arctic and melting of the Arctic sea ice lead to an intensified Beaufort Gyre, which in turn increases the southward volume transport through the Canadian Archipelago. To compensate for this southward transport of mass, more warm and saline Atlantic water is carried northward with the North Atlantic Drift. It is found that the increased transport of salt into the northern North Atlantic and the Nordic Seas tends to counteract the impact of the increased freshwater originating from the Arctic, leading to a stabilization of the North Atlantic Drift.     

An upper ocean general circulation model for climate studies: global simulation with seasonal cycle

The global ocean circulation with a seasonal cycle has been simulated with a two-and-a-half layer upper-ocean model. This model was developed for the purpose of coupling to an atmospheric general circulation model for climate studies on decadal time scales. The horizontal resolution is 4° latitude by 5° longitude and is thus not eddy-resolving. Effects of bottom topography are neglected. In the vertical, the model resolves the oceanic mixed layer and the thermocline. A thermodynamic sea-ice model is coupled to the mixed layer. The model is forced at the surface with seasonally varying (a) observed wind stress, (b) heat fluxes, as defined by an atmospheric equilibrium temperature, and (c) Newtonian-type surface salt fluxes. The second layer is coupled to the underlying deep ocean through Newtonian-type diffusive heat and salt fluxes, convective overturning, and mass entrainment in the upwelling regions of the subpolar gyres. The overall global distributions of mixed layer temperature, salinity and thickness are favorably reproduced. Inherent limitations due to coarse horizontal resolution result in large mixed-layer temperature errors near continental boundaries and in weak current systems. Sea ice distributions agree well with observations except in the interiors of the Ross and Weddell Seas. A realistic time rate of change of heat storage is simulated. There is also realistic heat transport from low to high latitudes.     

Sensitivity of simulated salinities in a three-dimensional ocean general circulation model to vertical mixing of destabilizing surface fluxes

 We present simulations performed with a three dimensional global ocean general circulation model which show that simulated salinities and amounts of convective mixing are very sensitive to vertical mixing of surface buoyancy fluxes. If, as usual, surface buoyancy fluxes are placed entirely in the topmost model level, our model produces excessive convective mixing in the Southern Ocean. This results in poor stimulated salinity in the Southern Ocean. In this simulation, we assume, as usual, that both surface buoyancy forcing and vertical mixing are homogeneous within each grid cell. If, on the other hand, destabilizing surface fluxes are instantaneously mixed into the subsurface ocean, the model produces much less convective mixing and much more realistic salinities. The vertical mixing of surface buoyancy fluxes performed in this simulation is equivalent to assuming that those fluxes affect only a small fraction of each grid cell, and cause vertical mixing only in that limited area. Our interpretation of these results is that the usual assumption that both surface buoyancy forcing and vertical mixing are uniform within each grid cell has a detrimental effect on model results; these results could be significantly improved by good parametrizations which treat the horizontal inhomogeneity of surface buoyancy forcing and of vertical mixing. Received: 25 February 1998 / Accepted: 9 September 1998     

The effect of Congo River freshwater discharge on Eastern Equatorial Atlantic climate variability

The surface ocean explains a considerable part of the inter-annual Tropical Atlantic variability. The present work makes use of observational datasets to investigate the effect of freshwater flow on sea surface salinity (SSS) and temperature (SST) in the Gulf of Guinea. In particular, the Congo River discharges a huge amount of freshwater into the ocean, affecting SSS in the Eastern Equatorial Atlantic (EEA) and stratifying the surface layers. The hypothesis is that an excess of river runoff emphasize stratification, influencing the ocean temperature. In fact, our findings show that SSTs in the Gulf of Guinea are warmer in summers following an anomalously high Congo spring discharge. Vice versa, when the river discharges low freshwater, a cold anomaly appears in the Gulf. The response of SST is not linear: temperature anomalies are considerable and long-lasting in the event of large freshwater flow, while in dry years they are less remarkable, although still significant. An excess of freshwater seems able to form a barrier layer, which inhibits vertical mixing and the entrainment of the cold thermocline water into the surface. Other processes may contribute to SST variability, among which the net input of atmospheric freshwater falling over EEA. Likewise the case of continental runoff from Congo River, warm anomalies occur after anomalously rainy seasons and low temperatures follow dry seasons, confirming the effect of freshwater on SST. However, the two sources of freshwater anomaly are not in phase, so that it is possible to split between atypical SST following continental freshwater anomalies and rainfall anomalies. Also, variations in air-sea fluxes can produce heating and cooling of the Gulf of Guinea. Nevertheless, atypical SSTs cannot be ascribed to fluxes, since the temperature variation induced by them is not sufficient to explain the SST anomalies appearing in the Gulf after anomalous peak discharges. The interaction processes between river runoff, sea surface salinity and temperature play an effective role in the interannual variability in the EEA region. Our results add a new source of variability in the area, which was often neglected by previous studies.     

Changes to the Canadian Arctic Archipelago Sea Ice and Freshwater Fluxes in the Twenty-First Century under the Intergovernmental Panel on Climate Change A1B Climate Scenario

A coupled ocean and sea-ice pan-Arctic model forced by the Intergovernmental Panel on Climate Change A1B climate scenario is used to study the evolution of ice and ocean surface conditions within the Canadian Arctic Archipelago (CAA) during the twenty-first century. A summer ice-free CAA is likely by the end of our simulation. Sea ice undergoes significant changes from the mid-2020s to the mid-2060s in both concentration and thickness. The simulation shows a shrinking of 65% and a thinning of 75% in summer over the 40 years, resulting in a partially open Northwest Passage by the 2050s. However, ice in central Parry Channel might increase due to a decrease in export from April to June, linked to a reduced cross-channel sea surface height (SSH) gradient, before melting thermodynamically. On a larger scale, the central CAA throughflow will experience a significant decrease in both volume and freshwater transport after 2020, which is related to the change in the SSH difference between the two ends of Parry Channel, particularly the lifting of SSH in Baffin Bay. With a lower albedo, a warmer ocean is simulated, particularly in summer. The sea surface salinity within the CAA demonstrates a strong decadal oscillation without a clear trend over the entire simulation. A north–south pattern, separated by Parry Channel, is also found in the changes of ocean temperature and salinity fields due to different ice conditions.     

ARE THERE INTERANNUAL-TO-DECADAL SCALE OSCILLATIONS ASSOCIATED WITH SEA ICE-THERMOHALINE CIRCULATION INTERACTIONS IN A SIMPLE COUPLED ATMO-SPHERE-OCEAN-SEA ICE MODEL?

1　INTRODUCTIONIn order to gain further insight into the nature of decadal- scale climate variability at highlatitudes( e.g.,Mysak et al.,1 990 ;Deser and black- mon,1 993) ,there have been a number ofrecent model studies of sea ice- thermohaline circulation interactions which exhibitoscillationson this timescale( Yang and Neelin,1 993;Zhang et al.,1 995 ;Yang and Huang,1 996 ) .Acommon feature of these studies is that the ocean models are integrated using mixedboundary conditions( MBC…     

Upward flushing of sea water through first year ice

Abstract

Observations in both the ice and slush layers suggest that sea water intrudes into the snow layer following a snow storm. Ice temperature values recorded at 1 cm below the snow‐ice interface show that the upward flux of sea water is of short duration. This is followed by a period of intense brine drainage characterized by the migration of a salty brine layer, with salinities up to 42 psu. These results suggest that a snow storm induces a complete (upward) flushing of the brine channel network and major modifications of snow and ice characteristics.

Melt rates and downward brine fluxes were calculated using salinity measured in a 40 cm deep box placed on the ice‐water interface, which isolated a 50 × 50 cm area of sea ice from ocean mixing processes. In this semi‐isolated environment, observed salinity changes allowed us to determine melt water fluxes and brine drainage or flushing even though ice thickness measurements did not show any significant change. Melt rates up to 21 cm/month and equivalent growth rates up to 32 cm/month were measured.     

Arctic climate change in 21st century CMIP5 simulations with EC-Earth

The Arctic climate change is analyzed in an ensemble of future projection simulations performed with the global coupled climate model EC-Earth2.3. EC-Earth simulates the twentieth century Arctic climate relatively well but the Arctic is about 2 K too cold and the sea ice thickness and extent are overestimated. In the twenty-first century, the results show a continuation and strengthening of the Arctic trends observed over the recent decades, which leads to a dramatically changed Arctic climate, especially in the high emission scenario RCP8.5. The annually averaged Arctic mean near-surface temperature increases by 12 K in RCP8.5, with largest warming in the Barents Sea region. The warming is most pronounced in winter and autumn and in the lower atmosphere. The Arctic winter temperature inversion is reduced in all scenarios and disappears in RCP8.5. The Arctic becomes ice free in September in all RCP8.5 simulations after a rapid reduction event without recovery around year 2060. Taking into account the overestimation of ice in the twentieth century, our model results indicate a likely ice-free Arctic in September around 2040. Sea ice reductions are most pronounced in the Barents Sea in all RCPs, which lead to the most dramatic changes in this region. Here, surface heat fluxes are strongly enhanced and the cloudiness is substantially decreased. The meridional heat flux into the Arctic is reduced in the atmosphere but increases in the ocean. This oceanic increase is dominated by an enhanced heat flux into the Barents Sea, which strongly contributes to the large sea ice reduction and surface-air warming in this region. Increased precipitation and river runoff lead to more freshwater input into the Arctic Ocean. However, most of the additional freshwater is stored in the Arctic Ocean while the total Arctic freshwater export only slightly increases.     

Dynamical contribution to sea surface salinity variations in the eastern Gulf of Guinea based on numerical modelling

In this study, we analyse the seasonal variability of the sea surface salinity (SSS) for two coastal regions of the Gulf of Guinea from 1995 to 2006 using a high resolution model (1/12°) embedded in a Tropical Atlantic (1/4°) model. Compared with observations and climatologies, our model demonstrates a good capability to reproduce the seasonal and spatial variations of the SSS and mixed layer depth. Sensitivity experiments are carried out to assess the respective impacts of precipitations and river discharge on the spatial structure and seasonal variations of the SSS in the eastern part of the Gulf of Guinea. In the Bight of Biafra, both precipitations and river runoffs are necessary to observe permanent low SSS values but the river discharge has the strongest impact on the seasonal variations of the SSS. South of the equator, the Congo river discharge alone is sufficient to explain most of the SSS structure and its seasonal variability. However, mixed layer budgets for salinity reveal the necessity to take into account the horizontal and vertical dynamics to explain the seasonal evolution of the salinity in the mixed layer. Indeed evaporation, precipitations and runoffs represent a relatively small contribution to the budgets locally at intraseasonal to seasonal time scales. Horizontal advection always contribute to spread the low salinity coastal waters offshore and thus decrease the salinity in the eastern Gulf of Guinea. For the Bight of Biafra and the Congo plume region, the strong seasonal increase of the SSS observed from May/June to August/September, when the trade winds intensify, results from a decreasing offshore spread of freshwater associated with an intensification of the salt input from the subsurface. In the Congo plume region, the subsurface salt comes mainly from advection due to a strong upwelling but for the Bight of Biafra, entrainment and vertical mixing also play a role. The seasonal evolution of horizontal advection in the Bight of Biafra is mainly driven by eddy correlations between salinity and velocities, but it is not the case in the Congo plume.     

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.     

Towards a parametrization of river discharges into ocean general circulation models: a closure through energy conservation

Diagnostic methods are defined in order to compare two numerical simulations of ocean dynamics in a region of freshwater influence. The first one is a river plume simulation based on a high resolution numerical configuration of the POM coastal ocean model in which mixing parametrizations have been previously defined. The second one is a simulation based on the NEMO Global Ocean Model used for climate simulations in its half-a-degree configuration in which a river inflow is represented as precipitation on two coastal grid cells. Both simulations are forced with the same freshwater inflows and wind stresses. The divergence of volumetric fluxes above and below the halocline are compared. Results show that when an upwelling wind blows, the two models display similar behavior although the impact of lack of precision can be observed in the NEMO configuration. When a downwelling wind blows, the NEMO Global Ocean configuration can not reproduce the coastally trapped baroclinic dynamics because its grid resolution is too coarse. To find a parametrization to help represent these dynamics in ocean general circulation models, a method based on energy conservation is investigated. This method shows that it is possible to link the energy fluxes provided by river inflows to the divergence of energy fluxes integrated over the grid cells of ocean general circulation models. A parametrization of the dynamics created by freshwater inflows is deduced from this method. This enabled creation of a box model that proved to have the same behavior as the fluxes previously computed from the high resolution configuration.     

Effects of freshwater runoff on a tropical pacific climate in the HadGEM2

This paper investigates the effects of river discharge on simulated climatology from 1979 to 1988 using the Hadley Centre Global Environmental Model version 2. Two experiments are performed with and without the inclusion of Total Runoff Integrating Pathways. The results show that the inclusion of flow routing can lead to the decrease of salinity over the coastal region due to freshwater. This reduction results in a shallower mixed layer depth, which in turn leads to the weakening of trade winds and a decrease in vertical mixing in the ocean. The enhanced sensible and latent heat fluxes over warmed SST improve the simulated precipitation and thermodynamic circulation. As a result, the experiment with flow routing is capable of improving the large-scale climate feature with an increase in precipitation over the eastern tropical equatorial Pacific region.     

Description and evaluation of the bergen climate model: ARPEGE coupled with MICOM

A new coupled atmosphere–ocean–sea ice model has been developed, named the Bergen Climate Model (BCM). It consists of the atmospheric model ARPEGE/IFS, together with a global version of the ocean model MICOM including a dynamic–thermodynamic sea ice model. The coupling between the two models uses the OASIS software package. The new model concept is described, and results from a 300-year control integration is evaluated against observational data. In BCM, both the atmosphere and the ocean components use grids which can be irregular and have non-matching coastlines. Much effort has been put into the development of optimal interpolation schemes between the models, in particular the non-trivial problem of flux conservation in the coastal areas. A flux adjustment technique has been applied to the heat and fresh-water fluxes. There is, however, a weak drift in global mean sea-surface temperature (SST) and sea-surface salinity (SSS) of respectively 0.1 °C and 0.02 psu per century. The model gives a realistic simulation of the radiation balance at the top-of-the-atmosphere, and the net surface fluxes of longwave, shortwave, and turbulent heat fluxes are within observed values. Both global and total zonal means of cloud cover and precipitation are fairly close to observations, and errors are mainly related to the strength and positioning of the Hadley cell. The mean sea-level pressure (SLP) is well simulated, and both the mean state and the interannual standard deviation show realistic features. The SST field is several degrees too cold in the equatorial upwelling area in the Pacific, and about 1 °C too warm along the eastern margins of the oceans, and in the polar regions. The deviation from Levitus salinity is typically 0.1 psu – 0.4 psu, with a tendency for positive anomalies in the Northern Hemisphere, and negative in the Southern Hemisphere. The sea-ice distribution is realistic, but with too thin ice in the Arctic Ocean and too small ice coverage in the Southern Ocean. These model deficiencies have a strong influence on the surface air temperatures in these regions. Horizontal oceanic mass transports are in the lower range of those observed. The strength of the meridional overturning in the Atlantic is 18 Sv. An analysis of the large-scale variability in the model climate reveals realistic El Niño – Southern Oscillation (ENSO) and North Atlantic–Arctic Oscillation (NAO/AO) characteristics in the SLP and surface temperatures, including spatial patterns, frequencies, and strength. While the NAO/AO spectrum is white in SLP and red in temperature, the ENSO spectrum shows an energy maximum near 3 years.     

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.     

Reducing Drift and Bias of a Global Ocean Model by Frequency-Dependent Nudging

Frequency-dependent nudging is applied to a coarse resolution (nominal 1°) global ocean model to suppress its drift and bias, and the impact of the nudging on the skill of the model is assessed. The nudging is applied to temperature and salinity in frequency bands centred on 0 and 1 cycles per year. As expected, the nudging significantly reduces the biases in the long-term mean and annual cycle of temperature, salinity, and sea level. By comparing the simulated (i) sea surface temperature with operational analyses based on observations, (ii) vertical profiles of temperature and salinity with observations made by Argo floats, and (iii) sea level with altimeter observations, it is shown that the skill of the model in simulating variability about the annual cycle is also improved. The potential benefit of applying frequency-dependent nudging to the ocean component of a coupled atmosphere–ocean model is discussed.     

Simulation and Improvements of Oceanic Circulation and Sea Ice by the Coupled Climate System Model FGOALS-f3-L

This study documents simulated oceanic circulations and sea ice by the coupled climate system model FGOALS-f3-L developed at the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, under historical forcing from phase 6 of the Coupled Model Intercomparison Project (CMIP6). FGOALS-f3-L reproduces the fundamental features of global oceanic circulations, such as sea surface temperature (SST), sea surface salinity (SSS), mixed layer depth (MLD), vertical temperature and salinity, and meridional overturning circulations. There are notable improvements compared with the previous version, FGOALS-s2, such as a reduction in warm SST biases near the western and eastern boundaries of oceans and salty SSS biases in the tropical western Atlantic and eastern boundaries, and a mitigation of deep MLD biases at high latitudes. However, several obvious biases remain. The most significant biases include cold SST biases in the northwestern Pacific (over 4°C), freshwater SSS biases and deep MLD biases in the subtropics, and temperature and salinity biases in deep ocean at high latitudes. The simulated sea ice shows a reasonable distribution but stronger seasonal cycle than observed. The spatial patterns of sea ice are more realistic in FGOALS-f3-L than its previous version because the latitude–longitude grid is replaced with a tripolar grid in the ocean and sea ice model. The most significant biases are the overestimated sea ice and underestimated SSS in the Labrador Sea and Barents Sea, which are related to the shallower MLD and weaker vertical mixing.