Relationship between sea surface salinity and ocean circulation and climate change

Based on Argo sea surface salinity(SSS) and the related precipitation(P), evaporation(E), and sea surface height data sets, the climatological annual mean and low-frequency variability in SSS in the global ocean and their relationship with ocean circulation and climate change were analyzed. Meanwhile, together with previous studies, a brief retrospect and prospect of seawater salinity were given in this work. Freshwater flux(E-P) dominated the mean pattern of SSS, while the dynamics of ocean circulation modulated the spatial structure and low-frequency variability in SSS in most regions. Under global warming, the trend in SSS indicated the intensification of the global hydrological cycle, and featured a decreasing trend at low and high latitudes and an increasing trend in subtropical regions. In the most recent two decades, global warming has slowed down, which is called the"global warming hiatus". The trend in SSS during this phase, which was different to that under global warming, mainly indicated the response of the ocean surface to the decadal and multi-decadal variability in the climate system, referring to the intensification of the Walker Circulation. The significant contrast of SSS trends between the western Pacific and the southeastern Indian Ocean suggested the importance of oceanic dynamics in the cross-basin interaction in recent decades. Ocean Rossby waves and the Indonesian Throughflow contributed to the freshening trend in SSS in the southeastern Indian Ocean, while the increasing trend in the southeastern Pacific and the decreasing trend in the northern Atlantic implied a long-term linear trend under global warming. In the future, higher resolution SSS data observed by satellites, together with Argo observations, will help to extend our knowledge on the dynamics of mesoscale eddies, regional oceanography, and climate change.     

Interpreting the sea surface temperature warming trend in the Yellow Sea and East China Sea

Previous studies have demonstrated that the low-frequency sea surface temperature (SST) variability in the Yellow Sea and East China Sea (YECS) is linked to large-scale climate variability, but explanations on the mechanisms vary. This study examines the low-frequency variability and trends of some atmospheric and oceanic variables to discuss their different effects on the YECS warming. The increasing temperature trend is also observed at a hydrographic section transecting the Kuroshio. The increasing rate of ocean temperature decreases with depth, which might result in an increase in vertical stratification and a decrease in vertical mixing, and thus plays a positive role on the YECS warming. The surface net heat flux (downward positive) displays a decreasing trend, which is possibly a result of the YECS warming, and, in turn, inhibits it. Wind speeds show different trends in different datasets, such that its role in the YECS warming is uncertain. The trends in wind stress divergence and curl have large uncertainties, so their effects on SST warming are still unclear. The Kuroshio heat transport calculated in this study, displays no significantly increasing trend, so is an unlikely explanation for the SST warming. Limited by sparse ocean observations, sophisticated assimilative climate models are still needed to unravel the mechanisms behind the YECS warming.     

Multi-decadal trends in the tropical Pacific western boundary currents retrieved from historical hydrological observations

As large-scale ocean circulation is a key regulator in the redistribution of oceanic energy, evaluating the multi-decadal trends in the western Pacific Ocean circulation under global warming is essential for not only understanding the basic physical processes but also predicting future climate change in the western Pacific. Employing the hydrological observations of World Ocean Atlas 2018(WOA18) from 1955 to 2017, this study calculated the geostrophic currents, volume transport and multidecadal trends for the North Equatorial Current(NEC), the North Equatorial Countercurrent(NECC), the Mindanao Current(MC), the Kuroshio Current(KC) in the origin and the New Guinea Coastal Undercurrent(NGCUC) within tropical western Pacific Ocean over multi-decades. Furthermore, this study examined the contributions of temperature and salinity variations. The results showed significant strengthening trends in NEC, MC and NGCUC over the past six decades, which is mainly contributed by temperature variations and consistent with the tendency in the dynamic height pattern. Zonal wind stress averaged over the western Pacific Ocean in the same latitude of each current represents the decadal variation and multi-decadal trends in corresponding ocean currents, indicating that the trade wind forcing plays an important role in the decadal trend in the tropical western Pacific circulation. Uncertainties in the observed hydrological data and trends in the currents over the tropical western Pacific are also discussed. Given that the WOA18 dataset covers most of the historical hydrological sampling data for the tropical western Pacific, this paper provides important observational information on the multi-decadal trend of the large-scale ocean circulation in the western Pacific.     

Physically Consistent Responses of the Global Atmospheric Hydrological Cycle in Models and Observations

Robust and physically understandable responses of the global atmospheric water cycle to a warming climate are presented. By considering interannual responses to changes in surface temperature (T), observations and AMIP5 simulations agree on an increase in column integrated water vapor at the rate 7 %/K (in line with the Clausius–Clapeyron equation) and of precipitation at the rate 2–3 %/K (in line with energetic constraints). Using simple and complex climate models, we demonstrate that radiative forcing by greenhouse gases is currently suppressing global precipitation (P) at ～?0.15 %/decade. Along with natural variability, this can explain why observed trends in global P over the period 1988?2008 are close to zero. Regional responses in the global water cycle are strongly constrained by changes in moisture fluxes. Model simulations show an increased moisture flux into the tropical wet region at 900 hPa and an enhanced outflow (of smaller magnitude) at around 600 hPa with warming. Moisture transport explains an increase in P in the wet tropical regions and small or negative changes in the dry regions of the subtropics in CMIP5 simulations of a warming climate. For AMIP5 simulations and satellite observations, the heaviest 5-day rainfall totals increase in intensity at ～15 %/K over the ocean with reductions at all percentiles over land. The climate change response in CMIP5 simulations shows consistent increases in P over ocean and land for the highest intensities, close to the Clausius?Clapeyron scaling of 7 %/K, while P declines for the lowest percentiles, indicating that interannual variability over land may not be a good proxy for climate change. The local changes in precipitation and its extremes are highly dependent upon small shifts in the large-scale atmospheric circulation and regional feedbacks.     

Sensitivity of climate response to variations in freshwater hosing location

In a recent intercomparison of the response of general circulation models (GCMs) to high-latitude freshwater forcing (Stouffer et al., J Climate 19(8):1365–1387, 2006), a number of the GCMs investigated showed a localised warming response in the high-latitude North Atlantic, as opposed to the cooling that the other models showed. We investigated the causes for this warming by testing the sensitivity of the meridional overturning circulation (MOC) to variations in freshwater forcing location, and then analysing in detail the causes of the warming. By analysing results from experiments with HadCM3, we are able to show that the high-latitude warming is independent of the exact location of the additional freshwater in the North Atlantic or Arctic Ocean basin. Instead, the addition of freshwater changes the circulation in the sub-polar gyre, which leads to enhanced advection of warm, saline, sub-surface water into the Greenland–Iceland–Norwegian Sea despite the overall slowdown of the MOC. This sub-surface water is brought to the surface by convection, where it leads to a strong warming of the surface waters and the overlying atmosphere.     

The impact of freshwater discharges on the ocean circulation in the Skagerrak/northern North Sea area Part I: model validation

We assess the performance of an eddy-recognizing numerical ocean model in simulating the pattern and variability of the hydrography in the Skagerrak/northern North Sea area. The model we use is a version of the widely used Princeton ocean model employing a terrain-following vertical coordinate. Results from a series of five multi-year simulations of the mesoscale response are described. The simulations differ in their representation of the lateral freshwater supply to the model ocean of which the first is a reference simulation. The next four are variations in which the river discharges and/or the Baltic outflow are given more realistic representations. For validation, we have used in situ hydrographic data. A novelty is that we use the concepts of freshwater height and potential energy anomaly as objective validation tools. We find that, in general, the model faithfully reproduces many of the observed hydrographic features including their mean patterns and their variance. Not surprisingly, we find that the Baltic outflow is by far the most significant freshwater source in terms of its influence on the hydrography in the area, a result corroborating earlier findings. The best validation is obtained when all freshwater supply is made as realistic as possible, in particular the Baltic outflow. We also find that the large scale cyclonic circulation and the location of fronts are robust characteristics of the Skagerrak/northern North Sea circulation given the impact changes in the freshwater input has on the hydrography. Finally, we find that a further exploration of the impact of the lateral open boundary forcing, e.g., the input of Atlantic water, is needed.     

Köppen Climate Types in Observed and Simulated Climates

The Köppen climate classification was applied to the observed gridded climatological sets and the outputs of four general circulation models (GCMs) over the continents of the Earth. All data had been acquired via the Data Distribution Centre established by the Intergovernmental Panel on Climate Change. The ability of the GCMs to simulate the Köppen climate zones identified in the real data was explored and possible future (global warming) changes in the climate types' distribution for each GCM were assessed. Differences in the area distributions derived from the GCMs' recent climate simulations give evidence about uncertainties generally involved in climate models. As to the global warming simulations, all GCM projections of warming climate (horizon 2050) show that the zones representing tropical rain climates and dry climates become larger, and the zones identified with boreal forest and snow climates together with the polar climates are smaller.     

Effects of climate change on thermal properties of lakes and reservoirs, and possible implications

Meteorologic-driven processes exert large and diverse impacts on lakes’ internal heating, cooling, and mixing. Thus, continued global warming and climate change will affect lakes’ thermal properties, dynamics, and ecosystem. The impact of climate change on Lake Tahoe (in the states of California and Nevada in the United States) is investigated here, as a case study of climate change effects on the physical processes occurring within a lake. In the Tahoe basin, air temperature data show upward trends and streamflow trends indicate earlier snowmelt. Precipitation in the basin is shifting from snow to rain, and the frequency of intense rainfall events is increasing. In-lake water temperature records of the past 38 years (1970–2007) show that Lake Tahoe is warming at an average rate of 0.013°C/year. The future trends of weather variables, such as air temperature, precipitation, longwave radiation, downward shortwave radiation, and wind speed are estimated from predictions of three General Circulation Models (GCMs) for the period 2001–2100. Future trends of weather variables of each GCM are found to be different to those of the other GCMs. A series of simulation years into the future (2000–2040) is established using streamflows and associated loadings, and meteorologic data sets for the period 1994–2004. Future simulation years and trends of weather variables are selected so that: (1) future simulated warming trend would be consistent with the observed warming trend (0.013°C/year); and (2) future mixing pattern frequency would closely match with the historical mixing pattern frequency. Results of 40-year simulations show that the lake continues to become warmer and more stable, and mixing is reduced. Continued warming in the Tahoe has important implications for efforts towards managing biodiversity and maintaining clarity of the lake.     

Simulating precipitation in the Northeast United States using a climate-informed K-nearest neighbour algorithm

Decadal prediction using climate models faces long-standing challenges. While global climate models may reproduce long-term shifts in climate due to external forcing, in the near term, they often fail to accurately simulate interannual climate variability, as well as seasonal variability, wet and dry spells, and persistence, which are essential for water resources management. We developed a new climate-informed K-nearest neighbour (K-NN)-based stochastic modelling approach to capture the long-term trend and variability while replicating intra-annual statistics. The climate-informed K-NN stochastic model utilizes historical data along with climate state information to provide improved simulations of weather for near-term regional projections. Daily precipitation and temperature simulations are based on analogue weather days that belong to years similar to the current year's climate state. The climate-informed K-NN stochastic model is tested using 53 weather stations in the Northeast United States with an evident monotonic trend in annual precipitation. The model is also compared to the original K-NN weather generator and ISIMIP-2b GFDL general circulation model bias-corrected output in a cross-validation mode. Results indicate that the climate-informed K-NN model provides improved simulations for dry and wet regimes, and better uncertainty bounds for annual average precipitation. The model also replicates the within-year rainfall statistics. For the 1961–1970 dry regime, the model captures annual average precipitation and the intra-annual coefficient of variation. For the 2005–2014 wet regime, the model replicates the monotonic trend and daily persistence in precipitation. These improved modelled precipitation time series can be used for accurately simulating near-term streamflow, which in turn can be used for short-term water resources planning and management.     

Using the Regional Ocean Modeling System (ROMS) to improve the ocean circulation from a GCM 20<Superscript>th</Superscript> century simulation

Global coupled climate models are generally capable of reproducing the observed trends in the globally averaged atmospheric temperature. However, the global models do not perform as well on regional scales. Here, we present results from a 20-year, high-resolution ocean model experiment for the Atlantic and Arctic Oceans. The atmospheric forcing is taken from the final 20 years of a twentieth-century control run with a coupled atmosphere–ocean general circulation model. The ocean model results from the regional ocean model are validated using observations of hydrography from repeat cruises in the Barents Sea. Validation is performed for average quantities and for probability distributions in space and time. The validation results reveal that, though the regional model is forced by a coupled global model that has a noticeable sea ice bias in the Barents Sea, the hydrography and its variability are reproduced with an encouraging quality. We attribute this improvement to the realistic transport of warm, salty waters into the Barents Sea in the regional model. These lateral fluxes in the ocean are severely underestimated by the global model. The added value with the regional model that we have documented here lends hope to advance the quality of oceanic climate change impact studies.     

A low-order model for the response of the Atlantic thermohaline circulation to climate change

Concern has been expressed that anthropogenic climate change may lead to a slowdown or even collapse of the Atlantic thermohaline circulation (THC). Because of the possibly severe consequences that such an event could have on the northern North Atlantic and northwestern Europe, integrated assessment models (IAMs) are needed to explore the associated political and socioeconomic implications. State-of-the-art climate models representing the THC are, however, often too complex to be incorporated into an integrated assessment framework. In this paper we present a low-order model of the Atlantic THC which meets the main requirements of IAMs: it (1) is physically based, (2) is computationally highly efficient, (3) allows for comprehensive uncertainty analysis and (4) can be linked to globally aggregated climate models that are mostly used in IAMs. The model is an interhemispheric extension of the seminal Stommel model. Its parameters are determined by a least-squares fit to the output of a coupled climate model of intermediate complexity. Results of a number of transient global warming simulations indicate that the model is able to reproduce many features of the behaviour of coupled ocean–atmosphere circulation models such as the sensitivity of the THC to the amount, regional distribution and rate of climate change.Responsible Editor: Richard Greatbatch     

Hydrodynamics and sediment dynamics in the German Bight. A focus on observations and numerical modelling in the East Frisian Wadden Sea

This work deals with analysis of hydrographic observations and results of numerical simulations. The data base includes acoustic Doppler current profilers (ADCP) observations, continuous measurements on data stations and satellite data originating from the medium resolution imaging spectrometer (MERIS) onboard the European Space Agency (ESA) satellite ENVISAT with a spatial resolution of 300 m. Numerical simulations use nested models with horizontal resolutions ranging from 1 km in the German Bight to 200 m in the East Frisian Wadden Sea coupled with a suspended matter transport model. Modern satellite observations have now a comparable horizontal resolution with high-resolution numerical model of the entire area of the East Frisian Wadden Sea allowing to describe and validate new and so far unknown patterns of sediment distribution. The two data sets are consistent and reveal an oscillatory behaviour of sediment pools to the north of the back-barrier basins and clear propagation patterns of tidally driven suspended particulate matter outflow into the North Sea. The good agreement between observations and simulations is convincing evidence that the model simulates the basic dynamics and sediment transport processes, which motivates its further use in hindcasting, as well as in the initial steps towards forecasting circulation and sediment dynamics in the coastal zone.     

Past and projected future changes in snowpack and soil frost at the Hubbard Brook Experimental Forest,New Hampshire,USA

Long‐term data from the Hubbard Brook Experimental Forest in New Hampshire show that air temperature has increased by about 1 °C over the last half century. The warmer climate has caused significant declines in snow depth, snow water equivalent and snow cover duration. Paradoxically, it has been suggested that warmer air temperatures may result in colder soils and more soil frost, as warming leads to a reduction in snow cover insulating soils during winter. Hubbard Brook has one of the longest records of direct field measurements of soil frost in the United States. Historical records show no long‐term trends in maximum annual frost depth, which is possibly confounded by high interannual variability and infrequency of major soil frost events. As a complement to field measurements, soil frost can be modelled reliably using knowledge of the physics of energy and water transfer. We simulated soil freezing and thawing to the year 2100 using a soil energy and water balance model driven by statistically downscaled climate change projections from three atmosphere‐ocean general circulation models under two emission scenarios. Results indicated no major changes in maximum annual frost depth and only a slight increase in number of freeze–thaw events. The most important change suggested by the model is a decline in the number of days with soil frost, stemming from a concurrent decline in the number of snow‐covered days. This shortening of the frost‐covered period has important implications for forest ecosystem processes such as tree phenology and growth, hydrological flowpaths during winter, and biogeochemical processes in soil. Published in 2010 by John Wiley & Sons, Ltd.     

Impact of GODAE products on nested HYCOM simulations of the West Florida Shelf

Nested non-assimilative simulations of the West Florida Shelf for 2004–2005 are used to quantify the impact of initial and boundary conditions provided by Global Ocean Data Assimilation Experiment ocean products. Simulations are nested within an optimum interpolation hindcast of the Atlantic Ocean, the initial test of the US Navy Coupled Ocean Data Assimilation system for the Gulf of Mexico, and a global ocean hindcast that used the latter assimilation system. These simulations are compared to one that is nested in a non-assimilative Gulf of Mexico model to document the importance of assimilation in the outer model. Simulations are evaluated by comparing model results to moored Acoustic Doppler Current Profiler measurements and moored sea surface temperature time series. The choice of outer model has little influence on simulated velocity fluctuations over the inner and middle shelf where fluctuations are dominated by the deterministic wind-driven response. Improvement is documented in the representation of alongshore flow variability over the outer shelf, driven in part by the intrusion of the Loop Current and associated cyclones at the shelf edge near the Dry Tortugas. This improvement was realized in the simulation nested in the global ocean hindcast, the only outer model choice that contained a realistic representation of Loop Current transport associated with basin-scale wind-driven gyre circulation and the Atlantic Meridional Overturning Circulation. For temperature, the non-assimilative outer model had a cold bias in the upper ocean that was substantially corrected in the data-assimilative outer models, leading to improved temperature representation in the simulations nested in the assimilative outer models.     

Influences of Indian Ocean interannual variability on different stages of El Nio: A FOAM1.5 model approach

Both the tropical Indian and tropical Pacific Oceans are active atmosphere-ocean interactive regions with robust interannual variability, which also constitutes a linkage between the two basins in the mode of variability. Using a global atmosphereocean coupled model, we conducted two experiments(CTRL and PC) to explore the contributions of Indian Ocean interannual sea surface temperature(SST) modes to the occurrence of El Ni?o events. The results show that interannual variability of the SST in the Indian Ocean induces a rapid growth of El Ni?o events during the boreal autumn in an El Ni?o developing year. However, it weakens El Ni?o events or even promotes cold phase conversions in an El Ni?o decaying year. Therefore, the entire period of the El Ni?o is shortened by the interannual variations of the Indian Ocean SST. Specifically, during the El Ni?o developing years, the positive Indian Ocean Dipole(IOD) events force an anomalous Walker circulation, which then enhances the existing westerly wind anomalies over the west Pacific. This will cause a warmer El Ni?o event, with some modulations by ocean advection and oceanic Rossby and Kelvin waves. However, with the onset of the South Asian monsoon, the Indian Ocean Basin(IOB) warming SST anomalies excite low level easterly wind anomalies over the west tropical Pacific during the El Ni?o decaying years. As a result, the El Ni?o event is prompted to change from a warm phase to a cold phase. At the same time, an associated atmospheric anticyclone anomaly appears and leads to a decreasing precipitation anomaly over the northwest Pacific. In summary, with remote forcing in the atmospheric circulation, the IOD mode usually affects the El Ni?o during the developing years, whereas the IOB mode affects the El Ni?o during the decaying years.     

A simple model of seasonal open ocean convection

 Aspects of open ocean deep convection variability are explored with a two-box model. In order to place the model in a region of parameter space relevant to the real ocean, it is fitted to observational data from the Labrador Sea. A systematic fit to OWS Bravo data allows us to determine the model parameters and to locate the position of the Labrador Sea on a stability diagram. The model suggests that the Labrador Sea is in a bistable regime where winter convection can be either “on” or “off ”, with both these possibilities being stable climate states. When shifting the surface buoyancy forcing slightly to warmer or fresher conditions, the only steady solution is one without winter convection.  We then introduce short-term variability by adding a noise term to the surface temperature forcing, turning the box model into a stochastic climate model. The surface forcing anomalies generated in this way induce jumps between the two model states. These state transitions occur on the interannual to decadal time scale. Changing the average surface forcing towards more buoyant conditions lowers the frequency of convection. However, convection becomes more frequent with stronger variability in the surface forcing. As part of the natural variability, there is a non-negligible probability for decadal interruptions of convection. The results highlight the role of surface forcing variability for the persistence of convection in the ocean.     

植被覆盖状况影响中国地表气温变化的观测事实

利用NOAA/AVHRR归一化植被指数（NDVI）及观测气温与再分析地表气温的差值（Observation Minus Reanalysis, OMR）分析了植被覆盖状况对中国地表气温变化的影响.结果表明,地表气温OMR趋势值与NDVI在空间上呈现出显著的负相关关系,植被覆盖状况差（NDVI小于0.1）的区域地表升温较为显著,气温OMR趋势值超过0.2℃/10a,而植被覆盖度高（NDVI大于0.5）的区域气温OMR趋势值则变化不大,甚至出现降温.气温OMR趋势值对植被的季节变化还有着敏感的响应.不同区域植被覆盖状况的差异可能导致中国地表气温变化对全球变暖的响应不同,预测中国未来气候变化需要考虑植被覆盖状况及其动态变化的影响.     

Downscaling IPCC control run and future scenario with focus on the Barents Sea

Global atmosphere-ocean general circulation models are the tool by which projections for climate changes due to radiative forcing scenarios have been produced. Further, regional atmospheric downscaling of the global models may be applied in order to evaluate the details in, e.g., temperature and precipitation patterns. Similarly, detailed regional information is needed in order to assess the implications of future climate change for the marine ecosystems. However, regional results for climate change in the ocean are sparse. We present the results for the circulation and hydrography of the Barents Sea from the ocean component of two global models and from a corresponding pair of regional model configurations. The global models used are the GISS AOM and the NCAR CCSM3. The ROMS ocean model is used for the regional downscaling of these results (ROMS-G and ROMS-N configurations, respectively). This investigation was undertaken in order to shed light on two questions that are essential in the context of regional ocean projections: (1) How should a regional model be set up in order to take advantage of the results from global projections; (2) What limits to quality in the results of regional models are imposed by the quality of global models? We approached the first question by initializing the ocean model in the control simulation by a realistic ocean analysis and specifying air-sea fluxes according to the results from the global models. For the projection simulation, the global models’ oceanic anomalies from their control simulation results were added upon initialization. Regarding the second question, the present set of simulations includes regional downscalings of the present-day climate as well as projected climate change. Thus, we study separately how downscaling changes the results in the control climate case, and how scenario results are changed. For the present-day climate, we find that downscaling reduces the differences in the Barents Sea between the original global models. Furthermore, the downscaled results are closer to observations. On the other hand, the downscaled results from the scenario simulations are significantly different: while the heat transport into the Barents Sea and the salinity distribution change modestly from control to scenario with ROMS-G, in ROMS-N the heat transport is much larger in the scenario simulation, and the water masses become much less saline. The lack of robustness in the results from the scenario simulations leads us to conclude that the results for the regional oceanic response to changes in the radiative forcing depend on the choice of AOGCM and is not settled. Consequently, the effect of climate change on the marine ecosystem of the Barents Sea is anything but certain.     

Internal and forced variability along a section between Greenland and Portugal in the CLIPPER Atlantic model

Numerical models are used to estimate the meridional overturning and transports along the paths of two hydrographic cruises, carried out in 1997 and 2002 from Greenland to Portugal. We have examined the influence of the different paths of the two cruises and found that it could explain 0.4 to 2 Sv of difference in overturning (the precise value is model-dependent). Models show a decrease in the overturning circulation between 1997 and 2002, with different amplitudes. The CLIPPER ATL6 model reproduces well the observed weakening of the overturning in density coordinates between the cruises; in the model, the change is due to the combination of interannual and high-frequency forcing and internal variability associated with eddies and meanders. Examination of the -coordinate overturning reveals model–data discrepancies: the vertical structure in the models does not change as much as the observed one. The East Greenland current variability is mainly wind-forced in the ATL6 model, while fluctuations due to eddies and instabilities explain a large part of the North Atlantic Current variability. The time-residual transport of dense water and heat due to eddy correlations between currents and properties is small across this section, which is normal to the direction of the main current.     

Oscillating glacial northern and southern deep water formation from combined neodymium and carbon isotopes

While ocean circulation is driven by the formation of deep water in the North Atlantic and the Circum-Antarctic, the role of southern-sourced deep water formation in climate change is poorly understood. Here we address the balance of northern- and southern-sourced waters in the South Atlantic through the last glacial period using neodymium isotope ratios of authigenic ferromanganese oxides in thirteen deep sea cores from throughout the South Atlantic. The data indicate that northern-sourced water did not reach the Southern Ocean during the late glacial, and was replaced by southern-derived intermediate and deep waters. The high-resolution neodymium isotope record (~ 300 yr sample spacing) from two spliced deep Cape Basin sites indicates that over the last glacial period northern-sourced water mass export to the Southern Ocean was stronger during the major Greenland millennial warming intervals (and Southern Hemisphere cool periods), and particularly during the major interstadials 8, 12, and 14. Northern-sourced water mass export was weaker during Greenland stadials and reached minima during Heinrich Events. The benthic foraminiferal carbon isotopes in the same Cape Basin core reflect a partial control by Southern Hemisphere climate changes and indicate that deep water formation and ventilation occurred in the Southern Ocean during major Greenland cooling intervals (stadials). Together, neodymium isotopes and benthic carbon isotopes provide new information about water mass sourcing and circulation in deep Southern Ocean waters during rapid glacial climate changes. Combining carbon and neodymium isotopes can be used to monitor the relative proportion of northern- and southern-sourced waters in the Cape Basin to gain insight into the processes which control the carbon isotopic composition of deep waters. In this study we show that deep water formation and circulation was more important than biological productivity and nutrient regeneration changes for controlling the carbon isotope chemistry of Antarctic Bottom Water during millennial-scale glacial climate cycles. This observation also lends support to the hypothesis that ocean circulation is linked to interhemispheric climate changes on short timescales, and that ventilation in the glacial ocean rapidly switched between the northern and Southern Hemisphere on millennial timescales.