The role of mean ocean salinity in climate

We describe numerical simulations designed to elucidate the role of mean ocean salinity in climate. Using a coupled atmosphere-ocean general circulation model, we study a 100-year sensitivity experiment in which the global-mean salinity is approximately doubled from its present observed value, by adding 35 psu everywhere in the ocean. The salinity increase produces a rapid global-mean sea-surface warming of C within a few years, caused by reduced vertical mixing associated with changes in cabbeling. The warming is followed by a gradual global-mean sea-surface cooling of C within a few decades, caused by an increase in the vertical (downward) component of the isopycnal diffusive heat flux. We find no evidence of impacts on the variability of the thermohaline circulation (THC) or El Niño/Southern Oscillation (ENSO). The mean strength of the Atlantic meridional overturning is reduced by 20% and the North Atlantic Deep Water penetrates less deeply. Nevertheless, our results dispute claims that higher salinities for the world ocean have profound consequences for the thermohaline circulation. In additional experiments with doubled atmospheric carbon dioxide, we find that the amplitude and spatial pattern of the global warming signal are modified in the hypersaline ocean. In particular, the equilibrated global-mean sea-surface temperature increase caused by doubling carbon dioxide is reduced by 10%. We infer the existence of a non-linear interaction between the climate responses to modified carbon dioxide and modified salinity.

An assessment of upper ocean salinity content from the ocean reanalyses inter-comparison project (ORA-IP)

Many institutions worldwide have developed ocean reanalyses systems (ORAs) utilizing a variety of ocean models and assimilation techniques. However, the quality of salinity reanalyses arising from the various ORAs has not yet been comprehensively assessed. In this study, we assess the upper ocean salinity content (depth-averaged over 0–700 m) from 14 ORAs and 3 objective ocean analysis systems (OOAs) as part of the Ocean Reanalyses Intercomparison Project. Our results show that the best agreement between estimates of salinity from different ORAs is obtained in the tropical Pacific, likely due to relatively abundant atmospheric and oceanic observations in this region. The largest disagreement in salinity reanalyses is in the Southern Ocean along the Antarctic circumpolar current as a consequence of the sparseness of both atmospheric and oceanic observations in this region. The West Pacific warm pool is the largest region where the signal to noise ratio of reanalysed salinity anomalies is >1. Therefore, the current salinity reanalyses in the tropical Pacific Ocean may be more reliable than those in the Southern Ocean and regions along the western boundary currents. Moreover, we found that the assimilation of salinity in ocean regions with relatively strong ocean fronts is still a common problem as seen in most ORAs. The impact of the Argo data on the salinity reanalyses is visible, especially within the upper 500m, where the interannual variability is large. The increasing trend in global-averaged salinity anomalies can only be found within the top 0–300m layer, but with quite large diversity among different ORAs. Beneath the 300m depth, the global-averaged salinity anomalies from most ORAs switch their trends from a slightly growing trend before 2002 to a decreasing trend after 2002. The rapid switch in the trend is most likely an artefact of the dramatic change in the observing system due to the implementation of Argo.

Salinity changes in the world ocean since 1950 in relation to changing surface freshwater fluxes

Global hydrographic and air–sea freshwater flux datasets are used to investigate ocean salinity changes over 1950–2010 in relation to surface freshwater flux. On multi-decadal timescales, surface salinity increases (decreases) in evaporation (precipitation) dominated regions, the Atlantic–Pacific salinity contrast increases, and the upper thermocline salinity maximum increases while the salinity minimum of intermediate waters decreases. Potential trends in E–P are examined for 1950–2010 (using two reanalyses) and 1979–2010 (using four reanalyses and two blended products). Large differences in the 1950–2010 E–P trend patterns are evident in several regions, particularly the North Atlantic. For 1979–2010 some coherency in the spatial change patterns is evident but there is still a large spread in trend magnitude and sign between the six E–P products. However, a robust pattern of increased E–P in the southern hemisphere subtropical gyres is seen in all products. There is also some evidence in the tropical Pacific for a link between the spatial change patterns of salinity and E–P associated with ENSO. The water cycle amplification rate over specific regions is subsequently inferred from the observed 3-D salinity change field using a salt conservation equation in variable isopycnal volumes, implicitly accounting for the migration of isopycnal surfaces. Inferred global changes of E–P over 1950–2010 amount to an increase of 1 ± 0.6 % in net evaporation across the subtropics and an increase of 4.2 ± 2 % in net precipitation across subpolar latitudes. Amplification rates are approximately doubled over 1979–2010, consistent with accelerated broad-scale warming but also coincident with much improved salinity sampling over the latter period.

On the climate response of the low-latitude Pacific Ocean to changes in the global freshwater cycle

Under global warming, the predicted intensification of the global freshwater cycle will modify the net freshwater flux at the ocean surface. Since the freshwater flux maintains ocean salinity structures, changes to the density-driven ocean circulation are likely. A modified ocean circulation could further alter the climate, potentially allowing rapid changes, as seen in the past. The relevant feedback mechanisms and timescales are poorly understood in detail, however, especially at low latitudes where the effects of salinity are relatively subtle. In an attempt to resolve some of these outstanding issues, we present an investigation of the climate response of the low-latitude Pacific region to changes in freshwater forcing. Initiated from the present-day thermohaline structure, a control run of a coupled ocean-atmosphere general circulation model is compared with a perturbation run in which the net freshwater flux is prescribed to be zero over the ocean. Such an extreme experiment helps to elucidate the general adjustment mechanisms and their timescales. The atmospheric greenhouse gas concentrations are held constant, and we restrict our attention to the adjustment of the upper 1,000 m of the Pacific Ocean between 40°N and 40°S, over 100 years. In the perturbation run, changes to the surface buoyancy, near-surface vertical mixing and mixed-layer depth are established within 1 year. Subsequently, relative to the control run, the surface of the low-latitude Pacific Ocean in the perturbation run warms by an average of 0.6°C, and the interior cools by up to 1.1°C, after a few decades. This vertical re-arrangement of the ocean heat content is shown to be achieved by a gradual shutdown of the heat flux due to isopycnal (i.e. along surfaces of constant density) mixing, the vertical component of which is downwards at low latitudes. This heat transfer depends crucially upon the existence of density-compensating temperature and salinity gradients on isopycnal surfaces. The timescale of the thermal changes in the perturbation run is therefore set by the timescale for the decay of isopycnal salinity gradients in response to the eliminated freshwater forcing, which we demonstrate to be around 10-20 years. Such isopycnal heat flux changes may play a role in the response of the low-latitude climate to a future accelerated freshwater cycle. Specifically, the mechanism appears to represent a weak negative sea surface temperature feedback, which we speculate might partially shield from view the anthropogenically-forced global warming signal at low latitudes. Furthermore, since the surface freshwater flux is shown to play a role in determining the ocean's thermal structure, it follows that evaporation and/or precipitation biases in general circulation models are likely to cause sea surface temperature biases.

Toward vicarious calibration of microwave remote-sensing satellites in arid environments

The Soil Moisture and Ocean Salinity (SMOS) satellite marks the commencement of dedicated global surface soil moisture missions, and the first mission to make passive microwave observations at L-band. On-orbit calibration is an essential part of the instrument calibration strategy, but on-board beam-filling targets are not practical for such large apertures. Therefore, areas to serve as vicarious calibration targets need to be identified. Such sites can only be identified through field experiments including both in situ and airborne measurements. For this purpose, two field experiments were performed in central Australia. Three areas are studied as follows: 1) Lake Eyre, a typically dry salt lake; 2) Wirrangula Hill, with sparse vegetation and a dense cover of surface rock; and 3) Simpson Desert, characterized by dry sand dunes. Of those sites, only Wirrangula Hill and the Simpson Desert are found to be potentially suitable targets, as they have a spatial variation in brightness temperatures of <4 K under normal conditions. However, some limitations are observed for the Simpson Desert, where a bias of 15 K in vertical and 20 K in horizontal polarization exists between model predictions and observations, suggesting a lack of understanding of the underlying physics in this environment. Subsequent comparison with model predictions indicates a SMOS bias of 5 K in vertical and 11 K in horizontal polarization, and an unbiased root mean square difference of 10 K in both polarizations for Wirrangula Hill. Most importantly, the SMOS observations show that the brightness temperature evolution is dominated by regular seasonal patterns and that precipitation events have only little impact.

Impact of brine-induced stratification on the glacial carbon cycle

During the cold period of the Last Glacial Maximum (LGM, about 21 000 years ago) atmospheric CO2 was around 190 ppm, much lower than the pre-industrial concentration of 280 ppm. The causes of this substantial drop remain partially unresolved, despite intense research. Understanding the origin of reduced atmospheric CO2 during glacial times is crucial to comprehend the evolution of the different carbon reservoirs within the Earth system (atmosphere, terrestrial biosphere and ocean). In this context, the ocean is believed to play a major role as it can store large amounts of carbon, especially in the abyss, which is a carbon reservoir that is thought to have expanded during glacial times. To create this larger reservoir, one possible mechanism is to produce very dense glacial waters, thereby stratifying the deep ocean and reducing the carbon exchange between the deep and upper ocean. The existence of such very dense waters has been inferred in the LGM deep Atlantic from sediment pore water salinity and δ18O inferred temperature. Based on these observations, we study the impact of a brine mechanism on the glacial carbon cycle. This mechanism relies on the formation and rapid sinking of brines, very salty water released during sea ice formation, which brings salty dense water down to the bottom of the ocean. It provides two major features: a direct link from the surface to the deep ocean along with an efficient way of setting a strong stratification. We show with the CLIMBER-2 carbon-climate model that such a brine mechanism can account for a significant decrease in atmospheric CO2 and contribute to the glacial-interglacial change. This mechanism can be amplified by low vertical diffusion resulting from the brine-induced stratification. The modeled glacial distribution of oceanic δ13C as well as the deep ocean salinity are substantially improved and better agree with reconstructions from sediment cores, suggesting that such a mechanism could have played an important role during glacial times.

Challenges in quantifying changes in the global water cycle

Understanding observed changes to the global water cycle is key to predicting future climate changes and their impacts. While many datasets document crucial variables such as precipitation, ocean salinity, runoff, and humidity, most are uncertain for determining long-term changes. In situ networks provide long time-series over land but are sparse in many regions, particularly the tropics. Satellite and reanalysis datasets provide global coverage, but their long-term stability is lacking. However, comparisons of changes among related variables can give insights into the robustness of observed changes. For example, ocean salinity, interpreted with an understanding of ocean processes, can help cross-validate precipitation. Observational evidence for human influences on the water cycle is emerging, but uncertainties resulting from internal variability and observational errors are too large to determine whether the observed and simulated changes are consistent. Improvements to the in situ and satellite observing networks that monitor the changing water cycle are required, yet continued data coverage is threatened by funding reductions. Uncertainty both in the role of anthropogenic aerosols, and due to large climate variability presently limits confidence in attribution of observed changes.

Reconstructing paleosalinity from δ18O: coupled model simulations of the Last Glacial Maximum, Last Interglacial and Late Holocene

Reconstructions of salinity are used to diagnose changes in the hydrological cycle and ocean circulation. A widely used method of determining past salinity uses oxygen isotope (δOw) residuals after the extraction of the global ice volume and temperature components. This method relies on a constant relationship between δOw and salinity throughout time. Here we use the isotope-enabled fully coupled General Circulation Model (GCM) HadCM3 to test the application of spatially and time-independent relationships in the reconstruction of past ocean salinity. Simulations of the Late Holocene (LH), Last Glacial Maximum (LGM), and Last Interglacial (LIG) climates are performed and benchmarked against existing compilations of stable oxygen isotopes in carbonates (δOc), which primarily reflect δOw and temperature. We find that HadCM3 produces an accurate representation of the surface ocean δOc distribution for the LH and LGM. Our simulations show considerable variability in spatial and temporal δOw-salinity relationships. Spatial gradients are generally shallower but within ∼50% of the actual simulated LH to LGM and LH to LIG temporal gradients and temporal gradients calculated from multi-decadal variability are generally shallower than both spatial and actual simulated gradients. The largest sources of uncertainty in salinity reconstructions are found to be caused by changes in regional freshwater budgets, ocean circulation, and sea ice regimes. These can cause errors in salinity estimates exceeding 4 psu. Our results suggest that paleosalinity reconstructions in the South Atlantic, Indian and Tropical Pacific Oceans should be most robust, since these regions exhibit relatively constant δOw-salinity relationships across spatial and temporal scales. Largest uncertainties will affect North Atlantic and high latitude paleosalinity reconstructions. Finally, the results show that it is difficult to generate reliable salinity estimates for regions of dynamic oceanography, such as the North Atlantic, without additional constraints.

The impact of salinity perturbations on the future uptake of heat by the Atlantic Ocean

Anthropogenic ocean heat uptake is a key factor in determining climate change and sea-level rise. There is considerable uncertainty in projections of freshwater forcing of the ocean, with the potential to influence ocean heat uptake. We investigatethis by adding either -0.1 Sv or +0.1 Sv freshwater to the Atlantic in global climate model simulations, simultaneously imposing an atmospheric CO2 increase. The resulting changes in the Atlantic meridional overturning circulation are roughly equal and opposite (±2Sv). The impact of the perturbation on ocean heat content is more complex, although it is relatively small (~5%) compared to the total anthropogenic heat uptake. Several competing processes either accelerate or retard warming at different depths. Whilst positive freshwater perturbations cause an overall heating of the Atlantic, negative perturbations produce insignificant net changes in heat content. The processes active in our model appear robust, although their net result is likely model- and experiment-dependent.

Decadal predictability of the Atlantic Ocean in a coupled GCM: forecast skill and optimal perturbations using Linear Inverse Modelling

The decadal predictability of three-dimensional Atlantic Ocean anomalies is examined in a coupled global climate model (HadCM3) using a Linear Inverse Modelling (LIM) approach. It is found that the evolution of temperature and salinity in the Atlantic, and the strength of the meridional overturning circulation (MOC), can be effectively described by a linear dynamical system forced by white noise. The forecasts produced using this linear model are more skillful than other reference forecasts for several decades. Furthermore, significant non-normal amplification is found under several different norms. The regions from which this growth occurs are found to be fairly shallow and located in the far North Atlantic. Initially, anomalies in the Nordic Seas impact the MOC, and the anomalies then grow to fill the entire Atlantic basin, especially at depth, over one to three decades. It is found that the structure of the optimal initial condition for amplification is sensitive to the norm employed, but the initial growth seems to be dominated by MOC-related basin scale changes, irrespective of the choice of norm. The consistent identification of the far North Atlantic as the most sensitive region for small perturbations suggests that additional observations in this region would be optimal for constraining decadal climate predictions.

Coupled simulations of the Indian monsoon intraseasonal oscillation using a fine-resolution mixed-layer ocean model

The intraseasonal variability (ISV) of the Indian summer monsoon is dominated by a 30–50 day oscillation between “active” and “break” events of enhanced and reduced rainfall over the subcontinent, respectively. These organized convective events form in the equatorial Indian Ocean and propagate north to India. Atmosphere–ocean coupled processes are thought to play a key role the intensity and propagation of these events. A high-resolution, coupled atmosphere–mixed-layer-oceanmodel is assembled: HadKPP. HadKPP comprises the Hadley Centre Atmospheric Model (HadAM3) and the K Profile Parameterization (KPP) mixed-layer ocean model. Following studies that upper-ocean vertical resolution and sub-diurnal coupling frequencies improve the simulation of ISV in SSTs, KPP is run at 1 m vertical resolution near the surface; the atmosphere and ocean are coupled every three hours. HadKPP accurately simulates the 30–50 day ISV in rainfall and SSTs over India and the Bay of Bengal, respectively, but suffers from low ISV on the equator. This is due to the HadAM3 convection scheme producing limited ISV in surface fluxes. HadKPP demonstrates little of the observed northward propagation of intraseasonal events, producing instead a standing oscillation. The lack of equatorial ISV in convection in HadAM3 constrains the ability of KPP to produce equatorial SST anomalies, which further weakens the ISV of convection. It is concluded that while atmosphere–ocean interactions are undoubtedly essential to an accurate simulation of ISV, they are not a panacea for model deficiencies. In regions where the atmospheric forcing is adequate, such as the Bay of Bengal, KPP produces SST anomalies that are comparable to the Tropical Rainfall Measuring Mission Microwave Imager (TMI) SST analyses in both their magnitude and their timing with respect to rainfall anomalies over India. HadKPP also displays a much-improved phase relationship between rainfall and SSTs over a HadAM3 ensemble forced by observed SSTs, when both are compared to observations. Coupling to mixed-layer models such as KPP has the potential to improve operational predictions of ISV, particularly when the persistence time of SST anomalies is shorter than the forecast lead time.

The meteorology of the Western Indian Ocean, and the influence of the east African highlands

This paper reviews the meteorology of the Western Indian Ocean and uses a state–of–the–art atmospheric general circulation model to investigate the influence of the East African Highlands on the climate of the Indian Ocean and its surrounding regions. The new 44–year re–analysis produced by the European Centre for Medium range Weather Forecasts (ECMWF) has been used to construct a new climatology of the Western Indian Ocean. A brief overview of the seasonal cycle of the Western Indian Ocean is presented which emphasizes the importance of the geography of the Indian Ocean basin for controlling the meteorology of the Western Indian Ocean. The principal modes of inter–annual variability are described, associated with El Niño and the Indian Ocean Dipole or Zonal Mode, and the basic characteristics of the subseasonal weather over the Western Indian Ocean are presented, including new statistics on cyclone tracks derived from the ECMWF re–analyses. Sensitivity experiments, in which the orographic effects of East Africa are removed, have shown that the East African Highlands, although not very high, play a significant role in the climate of Africa, India and Southeast Asia, and in the heat, salinity and momentum forcing of the Western Indian Ocean. The hydrological cycle over Africa is systematically enhanced in all seasons by the presence of the East African Highlands, and during the Asian summer monsoon there is a major redistribution of the rainfall across India and Southeast Asia. The implied impact of the East African Highlands on the ocean is substantial. The East African Highlands systematically freshen the tropical Indian Ocean, and act to focus the monsoon winds along the coast, leading to greater upwelling and cooler sea–surface temperatures.

The role of lateral ocean physics in the upper ocean thermal balance of a coupled ocean-atmosphere GCM

The sensitivity of the upper ocean thermal balance of an ocean-atmosphere coupled GCM to lateral ocean physics is assessed. Three 40-year simulations are performed using horizontal mixing, isopycnal mixing, and isopycnal mixing plus eddy induced advection. The thermal adjustment of the coupled system is quite different between the simulations, confirming the major role of ocean mixing on the heat balance of climate. The initial adjustment phase of the upper ocean (SST) is used to diagnose the physical mechanisms involved in each parametrisation. When the lateral ocean physics is modified, significant changes of SST are seen, mainly in the southern ocean. A heat budget of the annual mixed layer (defined as the “bowl”) shows that these changes are due to a modified heat transfer between the bowl and the ocean interior. This modified heat intake of the ocean interior is directly due to the modified lateral ocean physics. In isopycnal diffusion, this heat exchange, especially marked at mid-latitudes, is both due to an increased effective surface of diffusion and to the sign of the isopycnal gradients of temperature at the base of the bowl. As this gradient is proportional to the isopycnal gradient of salinity, this confirms the strong role of salinity in the thermal balance of the coupled system. The eddy induced advection also leads to increased exchanges between the bowl and the ocean interior. This is both due to the shape of the bowl and again to the existence of a salinity structure. The lateral ocean physics is shown to be a significant contributor to the exchanges between the diabatic and the adiabatic parts of the ocean.

Observation and origin of an interannual salinity anomaly in the Mozambique Channel

A positive salinity anomaly of 0.2 PSU was observed between 50 and 200 m over the years 2000–2001 across the Mozambique Channel at a section at 17°S which was repeated in 2003, 2005, 2006, and 2008. Meanwhile, a moored array is continued from 2003 to 2008. This anomaly was most distinct showing an interannual but nonseasonal variation. The possible origin of the anomaly is investigated using output from three ocean general circulation models (Estimating the Circulation and Climate of the Ocean, Ocean Circulation and Climate Advanced Modeling, and Parallel Ocean Program). The most probable mechanism for the salinity anomaly is the anomalous inflow of subtropical waters caused by a weakening of the northern part of the South Equatorial Current by weaker trade winds. This mechanism was found in all three numerical models. In addition, the numerical models indicate a possible salinization of one of the source water masses to the Mozambique Channel as an additional cause of the anomaly. The anomaly propagated southward into the Agulhas Current and northward along the African coast.

Changes to Indian Ocean subantarctic mode water in a coupled climate model as CO2 forcing increases

Subantarctic mode water (SAMW) has been shown to be a good indicator of anthropogenic climate change in coupled climate models. SAMW in a coupled climate model and the response of modeled SAMW to increasing CO2 are examined in detail. How SAMW adjusts from climatological values toward a new equilibrium in the coupled model, with different climatological temperature and salinity properties, is shown. The combined formation rate of SAMW and Antarctic intermediate water is calculated as approximately 18 Sv (Sv ≡ 106 m3 s−1) in the Indian sector of the Southern Ocean, slightly lower than climatological values would suggest. When forced with increasing CO2, SAMW is produced at a similar rate but at lower densities. This result suggests that the rate of heat uptake in this part of the ocean will be unchanged by anthropogenic forcing. The important signal in the response of SAMW is the shift to colder and fresher values on isopycnals that is believed to be related to changes in thermodynamic surface forcing. It is shown that, given uniform forcing, SAMW is expected to enhance the signal relative to other water masses. Independent increases in surface heating or freshwater forcing can produce changes similar to those observed, but the two different types of forcing are distinguishable using separate forcing experiments, hodographs, and passive anomaly tracers. The changes in SAMW forced by increasing CO2 are dominated by surface heating, but changes to freshwater fluxes are also important.