The roles of surface heat flux and ocean heat transport convergence in determining Atlantic Ocean temperature variability

The temperature variability of the Atlantic Ocean is investigated using an eddy-permitting (1/4°) global ocean model (ORCA-025) forced with historical surface meteorological fields from 1958 to 2001. The simulation of volume-averaged temperature and the vertical structure of the zonally averaged temperature trends are compared with those from observations. In regions with a high number of observations, in particular above a depth of 500 m and between 22° N and 65° N, the model simulation and the dataset are in good agreement. The relative contribution of variability in ocean heat transport (OHT) convergence and net surface heat flux to changes in ocean heat content is investigated with a focus on three regions: the subpolar and subtropical gyres and the tropics. The surface heat flux plays a relatively minor role in year-to-year changes in the subpolar and subtropical regions, but in the tropical North Atlantic, its role is of similar significance to the ocean heat transport convergence. The strongest signal during the study period is a cooling of the subpolar gyre between 1970 and 1990, which subsequently reversed as the mid-latitude OHT convergence transitioned from an anomalously weak to an anomalously strong state. We also explore whether model OHT anomalies can be linked to surface flux anomalies through a Hovmöller analysis of the Atlantic sector. At low latitudes, increased ocean heat gain coincides with anomalously strong northward transport, whereas at mid-high latitudes, reduced ocean heat loss is associated with anomalously weak heat transport.     

The relationships between tropical Atlantic sea level variability and major climate indices

Relationships were examined between variability in tropical Atlantic sea level and major climate indices with the use of TOPEX/POSEIDON altimeter and island tide gauge data with the aim of learning more about the external influences on the variability of the tropical Atlantic ocean. Possible important connections were found between indices related to the El Niño–Southern Oscillation (ENSO) and the sea levels in all three tropical regions (north, equatorial, and south), although the existence of only one major ENSO event within the decade of available altimetry means that a more complete investigation of the ENSO-dependence of Atlantic sea level changes has to await for the compilation of longer data sets. An additional link was found with the Indian Ocean Dipole (IOD) in the equatorial region, this perhaps surprising observation is probably an artifact of the similarity between IOD and ENSO time series in the 1990s. No evidence was obtained for significant correlations between tropical Atlantic sea level and North Atlantic Oscillation or Antarctic Oscillation Index. The most intriguing relationship observed was between the Quasi-Biennial Oscillation and sea level in a band centered approximately on 10°S. A plausible explanation for the relationship is lacking, but possibilities for further research are suggested.     

On the influence of adequate Weddell Sea characteristics in a large-scale global ocean circulation model

Global ocean circulation models usually lack an adequate consideration of high-latitude processes due to a limited model domain or insufficient resolution. Without the processes in key areas of the global thermohaline circulation, the characteristics and flow of deep and bottom waters cannot be modeled realistically. In this study, a high-resolution (~20 km) ocean model focused on the Weddell Sea sector of the Southern Ocean is combined with a low-resolution (2^° × 2^°) global ocean model applying the state estimation technique. Temperature, salinity, and velocity data on two Weddell Sea sections from the regional model are used as constraints for the large-scale model in addition to satellite altimetry and sea-surface temperatures. The differences between the model with additional constraints and without document that the Weddell Sea circulation exerts significant influence on the course of the Antarctic Circumpolar Current with consequences for Southern Ocean water mass characteristics and the spreading of deep and bottom waters in the South Atlantic. Furthermore, a warming trend in the period 1993–2001 was found in the Weddell Sea and adjacent basins in agreement with float measurements in the upper Southern Ocean. Teleconnections to the North Atlantic are suggested but need further studies to demonstrate their statistical significance.     

South Atlantic circulation in a world ocean model

The circulation in the South Atlantic Ocean has been simulated within a global ocean general circulation model. Preliminary analysis of the modelled ocean circulation in the region indicates a rather close agreement of the simulated upper ocean flows with conventional notions of the large-scale geostrophic currents in the region. The modelled South Atlantic Ocean witnesses the return flow and export of North Atlantic Deep Water (NADW) at its northern boundary, the inflow of a rather barotropic Antarctic Circumpolar Current (ACC) through the Drake Passage, and the inflow of warm saline Agulhas water around the Cape of Good Hope. The Agulhas leakage amounts to 8.7 Sv, within recent estimates of the mass transport shed westward at the Agulhas retroflection. Topographic steering of the ACC dominates the structure of flow in the circumpolar ocean. The Benguela Current is seen to be fed by a mixture of saline Indian Ocean water (originating from the Agulhas Current) and fresher Subantarctic surface water (originating in the ACC). The Benguela Current is seen to modify its flow and fate with depth; near the surface it flows north-westwards bifurcating most of its transport northward into the North Atlantic Ocean (for ultimate replacement of North Atlantic surface waters lost to the NADW conveyor). Deeper in the water column, more of the Benguela Current is destined to return with the Brazil Current, though northward flows are still generated where the Benguela Current extension encounters the coast of South America. At intermediate levels, these northward currents trace the flow of Antarctic Intermediate Water (AAIW) equatorward, though even more AAIW is seen to recirculate poleward in the subtropical gyre. In spite of the model’s rather coarse resolution, some subtle features of the Brazil-Malvinas Confluence are simulated rather well, including the latitude at which the two currents meet. Conceptual diagrams of the recirculation and interocean exchange of thermocline, intermediate and deep waters are constructed from an analysis of flows bound between isothermal and isobaric surfaces. This analysis shows how the return path of NADW is partitioned between a cold water route through the Drake Passage (6.5 Sv), a warm water route involving the Agulhas Current sheeding thermocline water westward (2.5 Sv), and a recirculation of intermediate water originating in the Indian Ocean (1.6 Sv).     

Anisotropic tomography of the Atlantic Ocean

We present the first regional three-dimensional model of the Atlantic Ocean with anisotropy. The model, derived from Rayleigh and Love wave phase velocity measurements, is defined from the Moho down to 300 km depth with a lateral resolution of about 500 km and is presented in terms of average isotropic S-wave velocity, azimuthal anisotropy and transverse isotropy.The cratons beneath North America, Brazil and Africa are clearly associated with fast S-wave velocity anomalies. The mid-Atlantic ridge (MAR) is a shallow structure in the north Atlantic corresponding to a negative velocity anomaly down to about 150 km depth. In contrast, the ridge negative signature is visible in the south Atlantic down to the deepest depth inverted, that is 300 km depth. This difference is probably related to the presence of hot-spots along or close to the ridge axis in the south Atlantic and may indicate a different mechanism for the ridge between the north and south Atlantic. Negative velocity anomalies are clearly associated with hot-spots from the surface down to at least 300 km depth, they are much broader than the supposed size of the hot-spots and seem to be connected along a north-south direction.Down to 100 km depth, a fast S-wave velocity anomaly is extenting from Africa into the Atlantic Ocean within the zone defined as the Africa superswell area. This result indicates that the hot material rising from below does not reach the surface in this area but may be pushing the lithosphere upward.In most parts of the Atlantic, the azimuthal anisotropy directions remain stable with increasing depth. Close to the ridge, the fast S-wave velocity direction is roughly parallel to the sea floor spreading direction. The hot-spot anisotropy signature is striking beneath Bermuda, Cape Verde and Fernando Noronha islands where the fast S-wave velocity direction seems to diverge radially from the hot-spots.The Atlantic average radial anisotropy is similar to that of the PREM model, that is positive down to about 220 km, but with slightly smaller amplitude and null deeper. Cratons have a lower than average radial anisotropy. As for the velocities, there is a difference between north and south Atlantic. Most hot-spots and the south-Atlantic ridge are associated with positive radial anisotropy perturbation whereas the north-Atlantic ridge corresponds to negative radial anisotropy perturbation.     

Subsidence, extension and thermal history of the West African margin in Senegal

The subsidence of the Atlantic margin in Senegal clearly shows two rapid stages related to the formation of (1) the Central Atlantic during the early Jurassic (around 200 Ma), and (2) the Equatorial Atlantic during the Cretaceous (100 Ma). A simple model of extension is used to interpret the subsidence history and to derive the thermal evolution of this basin. The present-day gravity, bathymetry, bottom hole temperatures (BHT) in oil exploration boreholes and heat flow density are used to control the validity of the model. Two cross sections from the outcropping basement to oceanic crust are used, one in Casamance and the other one at the south to latitude of Dakar. The model can fully explain the first-order subsidence history as well as the present-day observations, and therefore can provide valuable information about the thermal evolution of sediments and about the structure of the continental crust along the margin. Comparisons with the opposite margin in North America (Blake Plateau and Carolina trough) indicate a rather different evolution (the North American margin did not undergo the second stage of rifting) and a different crustal structure (crustal thinning is less important on the African margin).     

The 14C distribution in West Atlantic abyssal waters

¹⁴C specific activities in the western Atlantic show aging of about 160 years between 42°N and 30°S for southward-moving North Atlantic Deep Water. Most of the aging occurs in the North Atlantic, with a small increase in¹⁴C level for abyssal water near the equator.The northward-flowing Antarctic Bottom Water component ages about 80 years between 50°S and the equator. The rate of Atlantic bottom water formation is estimated at 18 Sverdrups; the rate of flow for northward-moving Antarctic Bottom Water at about 6 Sverdrups.     

Passive tracers in a general circulation model of the Southern Ocean

Passive tracers are used in an off-line version of the United Kingdom Fine Resolution Antarctic Model (FRAM) to highlight features of the circulation and provide information on the inter-ocean exchange of water masses. The use of passive tracers allows a picture to be built up of the deep circulation which is not readily apparent from examination of the veloCity or density fields. Comparison of observations with FRAM results gives good agreement for many features of the Southern Ocean circulation. Tracer distributions are consistent with the concept of a global “conveyor belt” with a return path via the Agulhas retroflection region for the replenishment of North Atlantic Deep Water.     

The distribution of 226Ra in the Atlantic Ocean

Based on results obtained during the GEOSECS program the primary features of the distribution of²²⁶Ra in the Atlantic Ocean can be defined. Outside the Antarctic no significant variation has been found in the²²⁶Ra content of surface waters. Eighty samples yield an average of 7.4 dpm/100 kg (normalized to a salinity of 35.00‰). Deep waters in the central Atlantic have²²⁶Ra contents several dpm/100 kg higher than expected from the mixing of Antarctic Bottom Water (21.3 dpm/100 kg) and basal North Atlantic Deep Water (10.3 dpm/100 kg). These excesses correlate well with deficiencies in O₂ and excesses in SiO₂. The intermediate water²²⁶Ra maximum in the South Atlantic is associated with the inflow of low-oxygen Circumpolar Intermediate Water beneath the Antarctic Intermediate Water.     

Heat flow and depth versus age for the Mesozoic northwest Atlantic Ocean: results from the Sohm abyssal plain and implications for the Bermuda Rise

Heat flow in the Sohm abyssal plain is measured to be 53 mW/m² at an age of 163 Ma. This is 25% higher than predicted by conductive cooling models, even though the sediment-corrected basement depth of 6.5 km at this location is normal for its age. An analysis of existing heat flow, depth and geoid anomalies in the northwest Atlantic shows that there is little correlation between heat flow and depth throughout the entire region. Depth and geoid are clearly related to the Bermuda swell while the associated heat flow anomaly, once adjusted for variations with age, is limited to 5 mW/m² and only decays to the south. This means that the Bermuda swell is probably not caused by extensive thermal reheating within the lithosphere, but instead by dynamic uplift at its lower boundary due to the convective upwelling of a mantle plume. The regionally high heat flow in the northwest Atlantic may be a thermal remanent of previous plumes which passed beneath this region early in its history. Therefore, depth and heat flow anomalies from this region cannot be used to provide constraints on steady-state parameters of the lithosphere, such as the presence or absence of a long-term boundary layer at its base.     

Regional and global effects of southern ocean constraints in a global model

Global ocean circulation models do not usually take high-latitude processes into account in an adequate form due to a limited model domain or insufficient resolution. Without the processes in key areas contributing to the lower part of the global thermohaline circulation, the characteristics and flow of deep and bottom waters often remain unrealistic in these models. In this study, various sections of the Bremerhaven Regional Ice Ocean Simulation model results are combined with a global inverse model by using temperature, salinity, and velocity constraints for the Hamburg Large Scale Geostrophic ocean general circulation model. The differences between the global model with and without additional constraints from the regional model demonstrate that the Weddell Sea circulation exerts a significant influence on the course of the Antarctic Circumpolar Current with consequences for Southern Ocean water mass characteristics and the spreading of deep and bottom waters in the South Atlantic. The influence of the Ross Sea is found to be less important in terms of global influences. However, regional changes in the Pacific sector of the Southern Ocean are found to be of Ross Sea origin. The additional constraints change the hydrographic conditions of the global model in the vicinity of the Antarctic Circumpolar Current in such a way that transport values, e.g., in Drake Passage no longer need to be prescribed to obtain observed transports. These changes not only improve the path and transport of the Antarctic Circumpolar Current but affect the meso- and large-scale circulation. With a higher (lower) mean Drake Passage transport, the mean Weddell Gyre transport is lower (higher). Furthermore, an increase (decrease) in the Antarctic Circumpolar Current leads to a decrease (increase) of the circum-Australian flow, i.e., a decrease (increase) of the Indonesian Throughflow.     

Mechanisms of the meridional heat transport in the Southern Ocean

The Southern Ocean (SO) transports heat towards Antarctica and plays an important role in determining the heat budget of the Antarctic climate system. A global ocean data synthesis product at eddy-permitting resolution from the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) project is used to estimate the meridional heat transport (MHT) in the SO and to analyze its mechanisms. Despite the intense eddy activity, we demonstrate that most of the poleward MHT in the SO is due to the time-mean fields of the meridional velocity, V, and potential temperature, θ. This is because the mean circulation in the SO is not strictly zonal. The Antarctic Circumpolar Current carries warm waters from the region south of the Agulhas Retroflection to the lower latitudes of the Drake Passage and the Malvinas Current carries cold waters northward along the Argentinian shelf. Correlations between the time-varying fields of V and θ (defined as transient processes) significantly contribute to the horizontal-gyre heat transport, but not the overturning heat transport. In the highly energetic regions of the Agulhas Retroflection and the Brazil-Malvinas Confluence the contribution of the horizontal transient processes to the total MHT exceeds the contribution of the mean horizontal flow. We show that the southward total MHT is mainly maintained by the meridional excursion of the mean geostrophic horizontal shear flow (i.e., deviation from the zonal average) associated with the Antarctic Circumpolar Current that balances the equatorward MHT due to the Ekman transport and provides a net poleward MHT in the SO. The Indian sector of the SO serves as the main pathway for the poleward MHT.     

Measurement of aerosol optical thickness over the Atlantic Ocean and in West Antarctica, 2006–2007

In this paper, optical measurements of aerosol properties made during a ship cruise from Poland to Antarctic Station in September and October 2006, and during the cruise back to Gdynia in April and May 2007 are described. A large gradient of pollution between the clear South Atlantic and the dusty North Atlantic was observed. The maximum of aerosol optical thickness at a wavelength of 500 nm reached 0.4 at 20°N in September 2006 and 0.3 at 40°N in May 2007, respectively. Strong Saharan dust transport is suggested as an explanation for the small values of Ångström exponent observed (values of 0.2 and 0.4 on these respective dates). On the Southern Hemisphere the aerosol optical thickness at 500 nm ranged from 0.05 to 0.2. Significant increases of the aerosol optical thickness were associated with strong wind and sea salt production. Good agreement was found when the in situ measurements of aerosol optical thickness were compared to satellite retrievals and modelling results.     

Particulate and soluble210Pb activities in the deep sea

Particulate and soluble,²¹⁰Pb activities have been measured by filtration of large-volume water samples at two stations in the South Atlantic. Particulate phase²¹⁰Pb (caught by a 0.4-μm filter) varies from 0.3% of total²¹⁰Pb in equatorial surface water to 15% in the bottom water. The “absolute activity” of²¹⁰Pb per unit mass of particulate matter is about 10⁷ times the activity of soluble²¹⁰Pb per unit mass of water, but because the mass ratio of particulate matter to water is about 10^?8, the particulate phase carries only about 10% of the total activity. In Antarctic surface water the particulate phase carries 40% of the total²¹⁰Pb activity; the absolute activity of this material is about the same as in other water masses and the higher fraction is due to the much larger concentration of suspended matter in surface water in this region.In the equatorial Atlantic the particulate phase²¹⁰Pb activity increases with depth, by a factor of 40 from surface to bottom, and by a factor of 4 from the Antarctic Intermediate Water core to the Antarctic Bottom Water. This increase with depth is predicted by our previously proposed particulate scavenging model which indicated a scavenging residence time of 50 years for²¹⁰Pb in the deep sea. A scavenging experiment showed that red clay sediment removes all the²¹⁰Pb from seawater in less than a week. The Antarctic particulate profile shows little or no evidence of scavenging in this region, which may be due to the siliceous nature of the particulate phase in circumpolar waters. Our previous observation that the²¹⁰Pb/²²⁶Ra activity ratio is of the order of 0.5 in the deep water is further confirmed by the two South Atlantic profiles analyzed in the present work.     

Rayleigh-wave group velocity distribution in the Antarctic region

We determined 2D group velocity distribution of Rayleigh waves at periods of 20-150 s in the Antarctic region using a tomographic inversion technique. The data are recorded by both permanent networks and temporary arrays. In East Antarctica the velocities are high at periods of 90-150 s, suggesting that the root of East Antarctica is very deep. The velocities in West Antarctica are low at all periods, which may be related to the volcanic activity and the West Antarctic Rift System. Low velocity anomalies appear at periods of 40-140 s along the Southeastern Indian Ridge and the western part of the Pacific Antarctic Ridge. The velocities are only slightly low around the Atlantic Indian Ridge, Southwestern Indian Ridge, and the eastern part of the Pacific Antarctic Ridge, where the spreading rates are small. Around two hotspots, the Mount Erebus and Balleny Islands, the velocity is low at periods of 50-150 s.     

Heat budget of the surface mixed layer south of Africa

ARGO hydrographic profiles, two hydrographic transects and satellite measurements of air–sea exchange parameters were used to characterize the properties and seasonal heat budget variations of the Surface Mixed Layer (SML) south of Africa. The analysis distinguishes the Subtropical domain (STZ) and the Subantarctic Zone (SAZ), Polar Frontal Zone (PFZ) and Antarctic Zone (AZ) of the Antarctic Circumpolar Current. While no Subantarctic Mode Water forms in that region, occurrences of deep SML (up to ∼450 m) are observed in the SAZ in anticyclones detached from the Agulhas Current retroflection or Agulhas Return Current. These are present latitudinally throughout the SAZ, but preferentially at longitudes 10–20° E where, according to previous results, the Subtropical Front is interrupted. Likely owing to this exchange window and to transfers at the Subantarctic Front also enhanced by the anticyclones, the SAZ shows a wide range of properties largely encroaching upon those of the neighbouring domains. Heat budget computations in each zone reveal significant meridional changes of regime. While air–sea heat fluxes dictate the heat budget seasonal variability everywhere, heat is mostly brought through lateral geostrophic advection by the Agulhas Current in the STZ, through lateral diffusion in the SAZ and through air–sea fluxes in the PFZ and AZ. The cooling contributions are by Ekman advection everywhere, lateral diffusion in the STZ (also favoured by the ∼10° breach in the Subtropical Front) and geostrophic advection in the SAZ. The latter likely reflects an eastward draining of water warmed through mixing of the subtropical eddies.     

A heat flow survey on anomaly M0 south of the Bermuda Rise

We have obtained a suite of 42 closely spaced, acoustically navigated, heat flow measurements on well-sedimented crust of anomaly M0 age (109 Ma) in the northwest Atlantic Ocean (25°N, 68°W; 950 km south of Bermuda). The mean and standard deviation of the values obtained are 1.13 HFU (μcal/cm² s) (47.3 mW/m²) and 0.05 HFU (2.1 mW/m²), respectively. Some of the variability is accounted for by refractive effects of the basement topography. Drill core data and our modelling suggest that the thermal conductivity contrast between sediments and basement rocks in this region is less than a factor of 1.6. The mean heat flow is close to the 1.1 HFU (46 mW/m²) predicted by both the plate and boundary layer cooling models of the oceanic lithosphere. This is the first detailed comparison with theoretical cooling models on old Atlantic Ocean crust. Since the difference in surface heat flow (0.15 HFU) predicted by the two cooling models for the oldest observed oceanic lithosphere (180 Ma) is also not much larger than the range of uncertainty in our observations, discrimination between the two models on the basis of surface heat flow data alone may prove difficult.     

Radium-barium-silica correlations and a two-dimensional radium model for the world ocean

Radium-barium-silica relationships are examined in various regions of the Pacific. Ra, Ba and Si are not linearly correlated except in the Circumpolar region. From surface water to the intermediate water, Ra usually increases linearly with Ba and Si. In the North Pacific where there is a Ra excess in the deep water, both the Ra-Ba and Ra-Si diagrams show curves concave upward with an increasing slope. The relationships are further complicated by a deep gentle Ba maximum and a broad mid-depth Si maximum.In the western boundary where the Antarctic Bottom Water and Pacific Deep Water are separated by the benthic front, the Ra-Ba and Ra-Si plots show “hooks” of various size and shape caused by a discordance in the depths of the Ra, Ba and Si maxima. The hook becomes smaller toward the south and flips from counterclockwise to clockwise in the southern basin.A two-dimensional horizontal model with a geostrophic circulation pattern is employed to calculate Ra distribution in the deep water with a constant upwelling rate but a variable mean input rate for Ra. The mean input rate or fluxF is the sum of the in-situ production rateJ by particulate dissolution in the water column and the bottom fluxQ_b. The global mean fluxF needed to balance the radiodecay in a steady-state distribution is used as a reference for discussion. It is found that anF value set at1.5 ×F fits the observed data in the deep Atlantic. TheF value needs to be about2.7 ×F in order to produce a comparable distribution in the deep Pacific. These model calculations indicate that a Ra source is required in the northeast Pacific. It appears that no uniqueF can produce a global Ra distribution in deep water comparable to both the Atlantic and Pacific observations. This is also true within the Pacific Ocean itself because of the large Ra excess observed in the northeast Pacific. Thus the Ra input rate is quite variable and increases by a factor of two from the Atlantic to the Pacific. This probably reflects the non-uniformity of the bottom fluxQ_b rather than that of the in-situ productionJ in the world oceans.     

Variability of the meridional overturning circulation of the North Atlantic: sensitivity to overflows of dense water masses

A numerical model of the Atlantic Ocean was used to study the low-frequency variability of meridional transports in the North Atlantic. The model shows a behaviour similar to those used in previous studies, and the temporal variability of certain variables compares favourably to observed time series. By changing the depth and width of the sills between the subpolar North Atlantic and the Nordic Seas, the mean horizontal and overturning circulation as well as some water mass properties are modified significantly. The reaction of meridional oceanic transports to atmospheric forcing fluctuations remains, however, unchanged. The critical role of the surface heat flux retroaction term for the meridional heat transport in stand-alone ocean models is discussed. The experiments underline the role of atmospheric variability for fluctuations of the large-scale ocean circulation on time scales from years to decades, and they support the hypothesis that the mean overturning strength is controlled by the model representation of the density of the overflow water masses.Responsible Editor: Dirk Olbers     

Western Boundary Current Separation Caused by a Deep Countercurrent

It is shown that the deep western boundary current (DWBC) in the North Atlantic may contribute to premature separation of the Gulf Stream. A simple analytical study is made based on the linear barotropic vorticity equation in which a possible effect of the DWBC is introduced as anomalous stretching of the water column in a manner consistent with certain qualitative aspects of the observed features. This additional positive vorticity may be generated by a deep downwelling associated with horizontal shears of the DWBC. It is found that separation latitudes are farther to the south as the stretching is strengthened. Vertical motion required for a southward shift of order unity (e.g., 500 km) is estimated to be of the order of 0.05 mm s m 1 , in reasonable agreement with the observed values. It is suggested that since the DWBC is intimately linked with the thermohaline circulation, a meridional migration of the Gulf Stream separation might reflect long-term variations in the North Atlantic.