In 1997, a unique hydrographic and chlorofluorocarbon (CFC: component CFC-11) dataset was obtained in the subpolar North Atlantic. To estimate the synopticity of the 1997 data, the recent temporal evolution of the CFC and Labrador Sea Water (LSW) thickness fields are examined. In the western Atlantic north of 50°N, the LSW thickness decreased considerably from 1994–97, while the mean CFC concentrations did not change much. South of 50°N and in the eastern Atlantic, the CFC concentration increased with little or no change in the LSW thickness. On shorter timescales, local anomalies due to the presence of eddies are observed, but for space scales larger than the eddies the dataset can be treated as being synoptic over the 1997 observation period.
The spreading of LSW in the subpolar North Atlantic is described in detail using gridded CFC and LSW thickness fields combined with Profiling Autonomous Lagrangian Circulation Explorer (PALACE) float trajectories. The gridded fields are also used to calculate the CFC-11 inventory in the LSW from 40° to 65°N, and from 10° to 60°W. In total, 2300 ± 250 tons of CFC-11 (equivalent to 16.6 million moles) were brought into the LSW by deep convection. In 1997, 28% of the inventory was still found in the Labrador Sea west of 45°W and 31% of the inventory was located in the eastern Atlantic.
The CFC inventory in the LSW was used to estimate the lower limits of LSW formation rates. At a constant formation rate, a value of 4.4–5.6 Sv (Sv ≡ 106 m3 s−1) is obtained. If the denser modes of LSW are ventilated only in periods with intense convection, the minimum formation rate of LSW in 1988–94 is 8.1–10.8 Sv, and 1.8–2.4 Sv in 1995–97.
Labrador Sea Water (LSW) is a main component of North Atlantic Deep Water (NADW) and contributes to the cold branch of the global thermohaline circulation. LSW is formed in the central Labrador Sea by deep convection, which occasionally reaches well below 1500 m. LSW is characterized by a minimum of salinity and potential vorticity and a maximum of chlorofluorocarbon (CFC) concentrations, and these signatures are found in the western and eastern subpolar North Atlantic. Three major pathways have been identified for LSW: spreading from the formation region to the northeast into the Irminger Sea, spreading eastward into the eastern North Atlantic, and spreading southeastward with the deep western boundary current (Talley and McCartney 1982; Sy et al. 1997).
Intensity and depth of convection in the Labrador Sea varied significantly during the past 50 years (Lazier 1988). Intensive production of LSW occured in the mid-1970s and again since the late 1980s. The modes of LSW formed in 1990–94 are the coldest and densest in this century. The intensity of convection is presumably linked to atmospheric variability, often presented as the NAO (North Atlantic Oscillation) index (e.g., Hurrell 1995). Intensive convective activity occurs during high NAO index periods, and weaker convection mostly during low NAO index periods (Dickson 1996). This simple picture might be obscured by variability of the freshwater fluxes and the variability in the stratification of the water column above the LSW. There is also some evidence that in some years deep convection might occur in the Irminger Sea, enhancing the area of LSW formation (R. Pickart 1999, personal communication).
In climate models, the variability of the LSW production rates causes variability in transport of the thermohaline circulation and thus in the meridional oceanic heat transport (e.g., Delworth et al. 1993; Delworth and Greatbatch 2000). On the other hand, in studies with high-resolution ocean-only models (Böning et al. 1996) the link between LSW formation rate and strength of the thermohaline overturning breaks down when the deep-water masses below the LSW are adequately represented. Wood et al. (1999) ran an improved oceanic climate model and obtained a more realistic representation of the main ocean currents (including the representation of overflows) than previous coupled ocean–atmosphere models. In this model, the meridional overturning in the Atlantic decreases by 20%–25% due to weakening of convection in the Labrador Sea when the greenhouse gases were quadrupled.
At present we still do not know whether in the real ocean the convection variability and thus the variability in LSW formation rate is linked to a variability in the meridional heat transport or variability in the transport of LSW out of the subpolar gyre. The knowledge of how much LSW is formed each year is crucial to address these issues. Mean LSW formation rates have been estimated from hydrographic data, from buoyancy fluxes, and from an a priori knowledge about the magnitude of the thermohaline circulation. The values range from 2 to 11 Sv (Sv ≡ 106 m3 s−1: Worthington 1976; Clarke and Gascard 1983; McCartney and Talley 1984; McCartney 1992; Speer and Tziperman 1992). Assuming that the overturning of the North Atlantic is entirely surface forced and neglecting mixing in the ocean, Marsh (2000) calculated time series of the production rate of LSW using climatologies of surface heat and freshwater fluxes, and he obtained a mean of 3.4 Sv for the 1980–97 time period, with a maximum value of 10 Sv in 1990 during the high NAO index period from 1988–94.
Smethie and Fine (2001) made the first effort to estimate rates of North Atlantic Deep Water formation from CFC data, collected between 1986 and 1992. After normalizing the observations to a common date, the CFC inventories of the various deep water masses were calculated. For LSW, a mean formation rate of 7.4 Sv was estimated for the 1970–90 time period. This estimate consisted of a formation rate that varied from 1.3 Sv to 12.5 Sv and was scaled to the observed thickness of LSW with greatest thickness and highest rates occuring during the mid 1970s and the late 1980s. Unfortunately, not enough CFC data exist to calculate the inventory for each year, making it necessary to assume a correlation with a feature that is available on a yearly timescale, such as the LSW thickness.
Here we combine hydrographic and CFC datasets obtained in the subpolar North Atlantic in 1997 (Fig. 1, Table 1), with Profiling Autonomous Lagrangian Circulation Explorer (PALACE) float trajectories. First we present the general CFC distribution in the water column and discuss the synopticity of the thickness and CFC fields. A topography-following interpolation scheme is applied to obtain a gridded dataset, which is used to discuss the spreading paths of LSW, to calculate the CFC inventory, and to estimate lower limits of LSW formation rates from 1988–97.
2. Dataset and methods
a. CFC analysis
Five different scientific groups sampled and analyzed the CFC-11 distribution in 1997 (Table 1). All water samples were analyzed on board with a gas chromatographic technique first described by Bullister and Weiss (1988). The precision of each dataset was checked by analyzing duplicates for 10%–15% of the samples and was found to be of the order of ±1%. Absolute calibration was done with standard gases calibrated against primary standards provided by R. F. Weiss, Scripps Institution of Oceanography, and D. W. R. Wallace, Institut für Meereskunde Kiel. The concentrations are reported on the SIO93 scale (Cunnold et al. 1994). The accuracies of the gas standards for CFC-11 range from 1% to 3%.
We choose the σΘ = 27.74 and 27.80 density surfaces to bound the LSW. The isopycnals correspond to σ1.5 = 34.60–34.63 and 34.695–34.705, respectively. These isopycnals are suitable to bound the LSW for the entire subpolar North Atlantic. The σΘ = 27.80 separates the LSW from the denser overflow waters, and the isopycnal is located between the salinity extrema of LSW and the overflow water. The σΘ = 27.74 separates the LSW from the lighter water masses, which are more saline than LSW.
The “classical” LSW (cLSW), from Smethie and Fine (2000), was bounded by σ1.5 = 34.62 and 34.69. These boundaries could not be taken for the 1997 dataset because the intense convection that occurred in the early 1990s created LSW with densities greater than σ1.5 = 34.69. In 1996–97, the σ1.5 = 34.69 is located in the salinity minimum of LSW, which is the characteristic core layer of LSW throughout the entire subpolar North Atlantic (Talley and McCartney 1982; Sy et al. 1997). In order to calculate the formation rates of two LSW modes separately, the LSW is split into two density classes, the sLSW (shallow LSW) bounded by the isopycnals σΘ = 27.74–27.77, and the dLSW (deep LSW), bounded by σΘ = 27.77–27.80. These should not be confused with the upper LSW (Pickart 1992), which is less dense than the LSW modes discussed here. The dLSW encloses the salinity minimum zone (about σΘ = 27.78) found everywhere in the subpolar North Atlantic, and sLSW contains the lighter modes of LSW. These lighter modes are renewed more frequently, and the ventilation of the dLSW by deep convection occured most likely in high NAO index periods (McCartney and Talley 1984; Lazier 1988; Dickson et al. 1996).
In total, 2410 individual CFC-11 measurements on 465 stations have been obtained in the density range of LSW. On average, five measurements per profile are available for the calculation of the mean CFC-11 concentration of LSW, such that the vertical resolution is sufficient to calculate the mean CFC-11 concentration on each profile for sLSW and dLSW. On 59 profiles only one or two CFC-11 samples in the LSW were available; 21 of these profiles were located at the continental shelf, where the water depth was too shallow to cover the whole density range of LSW. The remaining low-resolution profiles are found near the CGFZ (Charlie-Gibbs Fracture Zone) and above the Reykjanes Ridge. The LSW layer above the ridge is very thin (200m), so that the resolution there is sufficient, and the CGFZ profiles are surrounded by profiles with high vertical resolution. A more serious problem in calculating the mean CFC-11 concentrations for LSW might be the six low-resolution profiles clustered at 43°N, 47°W where the mean CFC values are determined by only one or two samples. For all profiles, the mean CFC-11 concentration of LSW is calculated by weighting the CFC values with the respective vertical distance they represent.
The tracer measurements in the subpolar North Atlantic in 1997 were carried out within an 8-month period. To provide a more complete dataset at the western boundary, two sections taken in August 1996 and one section in the southern Labrador Sea taken in August 1998 were added.
We checked the consistency of the CFC-11 datasets by comparing CFC concentrations west of 32.5°W corresponding to densities between σΘ = 27.77 and σΘ = 27.79 and falling in the same narrow temperature (±0.003) and salinity (±0.003) range. Convection did not reach this density layer in 1996, 1997, and 1998, and according to Sy et al. (1997) LSW from the Labrador Sea spreads into the western basin on timescales of 1 yr. The σΘ = 27.79 was chosen as a lower limit to minimize the effect of mixing with denser water masses (σΘ > 27.80), which have regionally different CFC concentrations.
For individual cruises the mean standard deviation of the CFC-11 values within the aforementioned Θ and S space is 2.5%–3%, composed of the measurement precision (about 1%) and regional differences. For the whole dataset the mean standard deviation increases to 4.5%; that is, the cruises are consistent on a 4% level. The major part of the intercruise differences are from the error in the accuracy of the CFC gas calibration standards (1%–3%).
3. Tracer distribution and temporal evolution of the tracer signal in the subpolar North Atlantic
a. CFC-11 distribution and the deep water masses
To give an overview on the vertical CFC field in the subpolar North Atlantic, the CFC-11 distribution along WOCE line A1E from 1997 is presented (Fig. 2). In the Irminger Sea, dLSW exhibits higher CFC levels than sLSW. Below the LSW, the Iceland–Scotland overflow water (ISOW) is visible in the Irminger Sea as a CFC minimum reaching from the Reykjanes Ridge to 40°W. The CFC maximum near the bottom west of 35°W represents the Denmark Strait overflow water (DSOW).
In the eastern Atlantic, the CFC levels are generally lower than in the western basins. In the mixed layer the decrease is caused by higher temperatures, allowing less CFCs to be dissolved. In middepth the CFC-11 maximum of the LSW is still a prominent signal, and the vertical CFC gradient between sLSW and dLSW is small. The LSW layer is very thin above and at the eastern flank of the Reykjanes Ridge (500 m), but thickens again in the Iceland Basin. A CFC maximum is found near 25°W at the eastern side of the Iceland Basin. The maximum is a permanent feature, which was already observed on two previous cruises in 1991 and 1994 (Sy et al. 1997; Smethie et al. 2000). At the eastern flank of the Reykjanes Ridge high CFC concentrations are found near the bottom, indicating the ISOW in the Iceland Basin. ISOW crosses into the western Atlantic through the CGFZ and, as mentioned before, forms the intermediate CFC minimum because in the western Atlantic ISOW encounters the CFC-richer deep water components LSW and DSOW. The lower CFC values of ISOW on the western flank of the Reykjanes Ridge compared to the ISOW on the eastern slope cannot be caused by mixing with LSW and DSOW, which are both richer in CFCs. The comparably weak CFC signal of ISOW at the western flank indicates that the main path of ISOW in the western Atlantic is not into the Irminger Sea as is commonly believed. Instead, ISOW spreads mainly into the Labrador Sea. A fraction of the ISOW remains in the Northeast Atlantic, flowing southward along the Mid-Atlantic Ridge (Smethie et al. 2000; Fleischmann et al. 2002). The very low CFC concentrations below the LSW in the eastern Atlantic characterize the Eastern Atlantic Bottom Water that has some Antarctic Bottom Water contributions (McCartney 1992).
b. Temporal evolution of the CFC-11 and thickness distribution of LSW
In order to get an idea whether we can treat the CFC concentrations in 1997 as being synoptic, we present here the temporal evolution of the CFC signal in the LSW since 1991. The CFC-11 concentrations in the sLSW west of the Reykjanes Ridge increased only slightly since 1991 (Fig. 3a), but in the eastern Atlantic the CFC level rose considerably between 1991 and 1994 and from 1994 to 1997. For dLSW the most dramatic change occured in the Irminger Sea between 1991 and 1994 caused by the invasion of newly ventilated CFC-rich LSW. (Fig. 3b). Since 1994 the CFC concentrations remained almost constant because the density layers of dLSW were not ventilated by deep convection. In contrast, owing to the longer travel time, formerly ventilated dLSW continued to invade the eastern Atlantic and the CFC levels increased steadily since 1991. Between 1991 and 1994 the increase was observed mainly west of 22°W (Iceland Basin), and in 1994–97 the increase was most pronounced in the eastern part of the section. This indicates the eastward propagation of LSW modes with high CFC concentrations.
The CFC-11 increase in the eastern Atlantic is about 0.1 pmol kg−1 yr−1 for the sLSW (period 1994–97). For the dLSW, the rates increase from 0.05 pmol kg−1 yr−1 near the Reykjanes Ridge to 0.3 pmol kg−1 yr−1 at 24°W and remain constant farther east. The regionally different rates of temporal CFC change from 1991 to 1997 show that it would not be appropriate to normalize the CFC data from these time periods to a common year by multiplying them with a constant factor as it was possible for the 1986–90 time period by Smethie and Fine (2001).
The LSW thickness did not change much in the eastern Atlantic from 1991 to 1997 (Fig. 3c). In the Irminger Sea the LSW was thicker in 1994 than in 1997 and in 1991. The thickness change is caused by the downward shift of the σΘ = 27.74 isopycnal. From 1994 to 1997 the LSW thickness decreased at a rate of 150–200 m yr−1. In the Labrador Sea the LSW thickness decreased from July 1996 to June 1997 by about 150–200 m, and the CFC-11 concentrations in the LSW decreased by 0.1–0.2 pmol kg−1 (differences between the observations for the cruises R/V Valdivia 161 in 1996 and M39-4 in 1997). Between 42° and 48°W along the WOCE A2 section, the CFC-11 increase from 1994 to 1997 was 0–0.1 pmol kg−1 yr−1 for the western basin, and 0.2–0.3 pmol kg−1 yr−1 for the eastern basin. In general, in 1994–97, the LSW variability in the western basin north of 50°N is dominated by a thickness decrease, and in the region south of 50°N and the eastern basins by an increase in the CFC-11 concentrations. The decrease of LSW thickness of about 150 m yr−1 in the western basins north of 50°N also means that the volume of LSW in the subpolar North Atlantic decreased since 1994. If one attributes the decrease solely to southward export and neglects mixing, the volume decrease would correspond to a southward LSW transport of 10 Sv.
4. Inferring gridded CFC-11 and thickness fields of the LSW
Having presented the CFC distribution and its temporal evolution, we now describe the methods applied to obtain gridded fields of CFC-11 and LSW thickness.
Both the CFC-11 distributions of LSW (Fig. 4) and the trajectories of the PALACE floats (Fig. 5) seem to follow the large-scale contours of the planetary potential vorticity f/H (f: Coriolis parameter, H: water depth, (Fig. 5). Thus, a topography-following mapping scheme was adapted from Davis (1998) to interpolate the horizontal CFC-11 and thickness distribution on a regular spatial grid (0.25° latitude by 0.5° longitude) from 40° to 65°N and from 60° to 10°W. The preferable method to perform such an interpolation would be objective mapping. Doing this for the CFC distribution needs an a priori knowledge of the data covariance that is not known. What is known from hydrographic and tracer observations, from float trajectories, and current measurements are the horizontal scales of these fields at key locations of the subpolar North Atlantic. Boundary current regimes are approximately 100 km wide with strong gradients along their offshore edges. There is a broad eddy field associated with the North Atlantic Current crossing the central basins (horizontal scales about 300 km). And there are some indications of the size of the convective area in the Labrador Sea (diameter about 300–500 km). Bearing these scales in mind, we constructed an interpolation scheme that takes into account the horizontal scales and the role of topography near the boundaries. “Gaussian” weights were determined as
with the “generalized” distance R defined as
In Eq. (2) a are the positions of grid points and b those of the measurements, such that the first term of the right-hand side of (2) is the squared distance between grid and data points divided by a horizontal scale factor (Λ). The second term takes account of the differences in barotropic potential vorticity:
with f the planetary vorticity and H the water depth, which is scaled by Φ. The scaling parameters were chosen (Λ = 0.5 × 10−4, Φ = 100) such that the first term in Eq. (2) has a half-width of 150 km (10% weight at 220 km) and is equivalent to a depth difference of 150 m in the second term. The interpolation method preserves the observed strong CFC gradients at the boundaries. Figure 6 gives examples of the weighting functions for various grid points. In the interior of the basins the weighting function is isotropic, and close to topographic slopes the weighting function follows the topography.
To illustrate the effect of the mapping procedure on the mean distributions and the regional variability we interpolated the data to their original positions with the same interpolation scheme. The interpolated data were virtually unchanged compared to the observations; that is, differences of the fields were negligible on average and distributed randomly in space. In the frequency distribution plot of the differences between interpolated and measured CFC-11 values and LSW thickness, the differences are symmetrically distributed around zero and the rms of the differences is 7% for CFC-11 and 12% for layer thickness relative to the mean CFC and thickness value, respectively (Fig. 7).
The calculated CFC-11 and LSW thickness fields are presented in Figs. 8 and 9. To assess the uncertainty in the two fields, we performed a variational statistical method known as jackknifing. The CFC-11 and the LSW thickness fields were calculated by randomly choosing 50% of the data, and the rms at each grid point was calculated from 100 realizations of the CFC-11 and the thickness field, respectively. The resulting rms scatter at each grid point is a measure of the combined effect of measurement errors, small-scale nonsynopticity, and low data coverage. The distributions of the rms values of the CFC-11 and the thickness field of LSW on the individual grid points are presented in Fig. 10.
For the CFC field, rms values exceeding 10% (red color) are found off Newfoundland near the continental slope and on the western side of the Mid-Atlantic Ridge between 42° and 48°N. The high rms is caused by poor sampling and/or by large CFC gradients in these regions. The rms in the thickness field is high in regions with large thickness gradients along the continental slopes north of 50°N, along the Reykjanes Ridge, and in the eastern basin at 20°W south of 45°N. Both presentations show regions with low data coverage associated with low error estimates. This somewhat counterintuitive result is caused by the small horizontal CFC and thickness gradients observed in these regions; that is, all randomly chosen observations lead to similar interpolation results in the aforementioned regions.
5. LSW pathways in the subpolar North Atlantic
Various distinct paths for the spreading of LSW can be deduced from Fig. 8. The highest CFC concentrations and the largest thickness (Fig. 9) are found in the central Labrador and Irminger Seas. High CFC values are also found along the western margin of the Newfoundland Basin. These features reflect the LSW formation in the central Labrador Sea and its export in the southeastward flowing deep western boundary current (DWBC), which has carried a CFC signal higher than 3 pmol kg−1 south to about 45°N. At 42°N, such high CFC values were restricted to a narrow 30–50 km band at the continental slope (water depth < 2500 m) and is only marginally visible in the map. The sharp CFC decrease offshore south of 48°N is caused by CFC-poorer water masses from subtropical regions flowing north with the North Atlantic Current (NAC). The limited horizontal extent of the LSW with high CFC concentrations at 42°N was also observed in August 1999 (R/V Meteor cruise M45/3). This southward export of LSW with high CFC values close to the continental slope most likely reflects the fast track of LSW into the subtropical Atlantic (Molinari et al. 1998).
CFC-11 values higher than 3 pmol kg−1 are also found in the Irminger Sea, indicative of the rapid spreading of LSW from the formation region to the Irminger Sea (timescale 6 months: Sy et al. 1997). The spreading into the Irminger Sea occurs north of 56°N in the opposite direction of the DWBC, which is located at the Greenland continental slope. The relative CFC maximum in the central Irminger Sea is found at the location of a gyre (inferred from the mean flow field derived from float trajectories (K. Lavender 1999, personal communication)], probably preserving the interior of the gyre from fast exchange with the periphery. There is evidence that deep convection in the Irminger Sea (R. Pickart 1999, personal communication) might be responsible for the CFC maximum there. On the other hand, the years 1995–97 were characterized by weak convection activities, and in the tracer fields of 1999 no CFC maximum was found in the Irminger Sea.
The relatively CFC poor water tongue spreading into the Labrador Sea to about 50°W reflects the influence of the NAC, meandering into the Labrador Sea in the “Northwest Corner” (Rossby 1996). Tongues of CFC concentrations between 2.7 and 3.0 pmol kg−1 enter the eastern Atlantic south of 53°N. Following the trajectories of the PALACE floats (Fig. 5), at least a part of this water flows first along the continental slope and at about 50°N turns north and east around topographic features (Fischer and Schott 2002) and then turns eastward and crosses into the eastern Atlantic south of 53°N. The LSW mode proceeding into the eastern Atlantic is CFC poorer and warmer than the LSW spreading into the Irminger Sea due to contributions of the warmer and CFC-poorer deep extension of the NAC coming from subtropical regions. The LSW crosses the western basin toward the CGFZ in about two years. There is a weak indication that a fraction of the LSW flows southward along the western flank of the Mid-Atlantic Ridge. In Fig. 4, one station located at the eastern flank of the ridge at about 46°N, 31°W exhibits significantly higher values than the data in the surrounding area, indicative of a southward LSW flow. Due to its small horizontal extension, this feature is not present in the interpolated field.
In the eastern Atlantic, LSW spreads mainly northeastward to the eastern flank of the Iceland Basin, following the approximate path of the North Atlantic Current. One part seems to flow farther east and south into the Porcupine Abyssal Plain, and the CFC concentrations decrease further. Smethie et al. (2000) and Fleischmann et al. (2002) reported that not all ISOW enters the western Atlantic through the CGFZ, but partly flows southward along the eastern flank of the Mid-Atlantic Ridge. It is not yet clear if the LSW also follows this path. The lowest CFC values were found farther south, presumably due to the influence of saline and CFC-poor Mediterranean Water.
The large-scale circulation features summarized by Talley and McCartney (1982) from hydrographic data collected in 1955–64 are similar to the ones reported here. Tongues of low potential vorticity values characteristic for LSW spread from the Labrador Sea northeastward into the Irminger Sea, cross into the northeast Atlantic south of 53°N, head toward the Rockall Plateau and then east and south into the Porcupine Abyssal Plain, and enter the subtropical gyre in a narrow region close to the continental slope. The southward path along the western flank of the Mid-Atlantic Ridge is not visible in the potential vorticity distribution. Owing to a very coarse resolution in this region, the CFC field there is only weakly constrained by data, and more measurements are needed to confirm this LSW pathway. In the Irminger Sea, the LSW flow scheme that Schmitz and McCartney (1993, their Fig. 11) derived is different from the one inferred from the CFC distribution. The CFC gradients show clearly that the main source of LSW in the Irminger Basin is the Labrador Sea and not the eastern Atlantic.
6. CFC-11 inventory of the subpolar North Atlantic
The CFC-11 inventory (INV) of LSW in the subpolar North Atlantic was calculated from the gridded CFC-11 and thickness fields by multiplying the average CFC-11 concentration (C) and thickness (D) at each grid point by the area (A) of the grid point and the water density and summing over all grid points. The inventories for sLSW, dLSW, and LSW were calculated separately:
The CFC-11 inventory for LSW is 16.6 ± 1.8 million moles, which is equivalent to 2300 ± 250 tons (Table 2). About 69% of the inventory is found in the western basin and 31% in the eastern basin, although the areas where LSW is present are similar for both basins; 28% of the inventory is still located in the Labrador Sea west of 45°W. The higher inventory in the western basin is caused by a greater thickness (Fig. 8) of the LSW layer (western basin: 1210 m, eastern basin: 740 m) and by the higher CFC-11 concentrations in most of the western North Atlantic (Fig. 7). The contributions of both sLSW and dLSW reflect the LSW formation and input of CFC-11 in the western basin and the fast spreading of LSW into the Irminger Sea and southward in the DWBC. The CFC-11 inventory in the dLSW is higher than in sLSW. This feature mainly reflects the different thicknesses of the layers: the sLSW has a mean thickness of 410 m, the dLSW of 600 m. Following the temporal trends, the inventory of the eastern Atlantic should increase in the next years, and as long as the deep convection in the Labrador Sea is not intensified, the CFC inventory in the western Atlantic should decrease. Smethie and Fine (2001) estimated the CFC-11 inventory in LSW for the subpolar North Atlantic in 1990 to be 12.4 million moles. This is smaller than our value and shows that CFCs continued to enter LSW between 1990 and 1997 faster than they were exported out of the subpolar gyre.
a. Uncertainties in the inventory
Compared to the CFC inventory calculations from Smethie and Fine (2001), the observations in 1997 have various advantages. The dataset is much larger, so that the CFC distribution in LSW can be studied in more detail, and the errors in the inventories can be better assessed. Since the data are synoptic, the measured CFC concentrations and the LSW thickness fields need not be normalized to a common date. We have shown that from 1991–97 the temporal evolution of the CFC and the thickness fields are regionally different and the latter procedure would be difficult to apply.
We already have assessed the uncertainty in the CFC-11 and thickness fields. The uncertainties of the CFC-11 inventories from the combined error on the CFC-11 and thickness fields was calculated to ±10% (LSW), 12% (dLSW), and 14% (sLSW).
b. Error from incompleteness of the gridded fields
CFC concentrations and LSW thickness values could not be calculated on some grid points due to lack of observations. About 12% of these grid points have water depths greater than 1000 m, so we can certainly assume that LSW is present. Here 1.4% are located in the high CFC and high thickness region of the northern Labrador Sea, but the majority are in the southernmost part of the western and eastern basins, which are both low CFC and low thickness regions (Figs. 8 and 9). Extrapolating the CFC and thickness fields from the surrounding areas to fill these regions, the northern Labrador Sea may add another 3.6% to the CFC inventory and the other areas could contribute 3% (mean thickness 700 m, mean CFC 1 pmol kg−1). Both contributions were neglected in the CFC inventory shown in Table 2. If we take them into account, the CFC inventory increases to 17.7 × 106 moles.
c. Error from the nonsynopticity of the data
Another error might arise from the large-scale nonsynopticity of the dataset. As already discussed, the LSW variability (1994–97) in the western Atlantic north of 50°N is dominated by a thickness decrease of about 100–200 m yr−1 (Irminger Sea, Labrador Sea), the region south of 50°N by a CFC increase of 0.05–0.1 pmol kg−1 yr−1, and the eastern Atlantic by a CFC-11 increase of 0.2–0.3 pmol kg−1 yr−1. However, the thickness and CFC-11 fields in the Labrador Sea are already composites of three cruises (Table 1). From the February–March cruise (1997), all profiles where the density in the upper 300 m exceeds σΘ = 27.74 were omitted from the thickness calculations and the calculation of the CFC fields. The fields in the Irminger Sea and in the eastern Atlantic north of 50°N are a composite of four cruises. For all these regions, the effect of the nonsynopticity is already covered by the jackknifing procedure.
South of 50°N, spatial resolution is coarser, and the fields are determined from only two cruises and from two sections in the western Atlantic carried out in August 1996. The contribution of the area south of 50°N to the CFC-11 inventory [σΘ = 27.74(+0.06)] is 30%, with 18% from the western basin and 12% from the eastern basin. The annual repeat sections along WOCE line A2 (nominally 48°N) do not show a temporal trend of the LSW thickness but, as mentioned above, a CFC increase that is stronger in the eastern basin. A CFC-11 increase of 0.05 pmol kg−1 per 5 months in the western Atlantic south of 50°N and of 0.15 pmol kg−1 per 5 months in the eastern Atlantic south of 50°N would increase the inventory of these areas by roughly 2.5% and 10%, respectively. The total CFC-11 inventory would then increase by about 2%. Thus, we estimate the error arising from the large-scale nonsynopticity of the data to be smaller than 3% and negligible compared to the other uncertainties.
a. Assumptions for calculating LSW formation rates
The CFC-11 inventory is used to estimate a lower limit of LSW production. For this purpose it is assumed that all CFC-11 in the LSW layer was introduced by the process of deep convection. The exact site of the formation area is irrelevant to the calculations.
If this assumption is true, the CFC-11 concentrations should increase with decreasing salinity in the subpolar North Atlantic, since the water masses surrounding LSW are in general CFC poorer and more saline than LSW. Figure 11 presents the CFC-11–salinity correlation of the observations used for the calculation of the inventory. The dominating correlation is an inverse linear relationship between CFC-11 and salinity. The two linear mixing lines between the salinity minimum of LSW and the upper water column and between the salinity minimum and the deeper waters cannot be distinguished.
However, there are three deviations from the linear correlation representing LSW (Fig. 11). The lowest branch, at CFC-11 = 1.2–1.5 pmol kg−1, is located in the eastern Atlantic south of 46°N (blue dots in Fig. 11) and indicates mixing with saline Mediterranean Water with CFC-11 concentrations around 1.5 pmol kg−1. The area where this correlation is found is small and does not significantly contribute to the inventory (<5%). Moreover, the higher CFC-11 concentrations in this area do follow the LSW correlation.
The deviation at about 2.5 pmol kg−1 (red dots in Fig. 11) indicates mixing of LSW with ISOW north of 58°N at the eastern flank of the Reykjanes Ridge. Here again, the area involved is small, the LSW thickness is only 100–200 m (Fig. 8b), and the higher CFC-11 concentrations tend to follow the LSW correlation. Farther south, the linear LSW correlation prevails.
The third deviation at about 4.2 pmol kg−1 (green dots in Fig. 11) is caused by mixing of sLSW with saline water of Atlantic origin, located between the cold and fresh East Greenland Current and the LSW. The mixing occurs at the Greenland continental slope north of 60°N in a small area and is not observed farther south.
We conclude from the dominating LSW correlation between CFC-11 and salinity that the contributions to the CFC-11 inventory by mixing with Mediterranean Water south of 46°N, with Iceland–Scotland overflow water north of 48°N (eastern Atlantic), and with Atlantic Water north of 60°N (western Atlantic) are small and can be neglected in the calculations of the LSW formation rates. The import of water from the subtropical gyre by the deep extension of the NAC is recirculated LSW plus some mixture of other water masses, but the latter are CFC poor and their impact can be neglected.
We define the subpolar North Atlantic as a closed basin where the densities that bound the LSW, the sLSW, and the dLSW (σΘ = 27.74, 27.77, and 27.80) are material surfaces outside the LSW formation region. Only one source of CFC exists for the LSW, and that is deep convection in the Labrador Sea. During convection, the water located above the LSW becomes denser than σΘ = 27.74 and is thus a new volume of LSW. Note that in this calculation the term “newly formed LSW” means the water volume that was less dense than LSW before convection starts and is then transferred to the density classes of LSW by buoyancy loss. It should not be mistaken as the conversion of a warmer, saline, and CFC poorer LSW to a colder, less saline, and CFC richer LSW mode. This is shown in Fig. 12. As an example, the density and CFC profiles from the central Labrador Sea (cruise M39-4, see Table 1) are presented in Figs. 12a and 12b. The sharp increase in the CFC concentrations is caused by the warm summer surface layer overlying the colder layers below. The LSW bounded by the isopycnals 27.74–27.80 is shaded and extends from 600 to 2300 m. Figure 12c shows the density profile after deep convection to about 1900 m. The water column above 600 m is now denser than σΘ = 27.74 and is thus a new volume of LSW that is added to the “old” LSW volume. The water below 600 m is modified—it becomes cooler, fresher and CFC richer—but it is only a conversion between LSW modes, and the higher CFC values come from the mixing with the new volume of LSW.
b. CFC saturation of new LSW
Observations in February–March in deep convection regions of the Greenland Sea have shown that in the intermediate deep mixed layer before the onset of deep convection (i.e., before deep convection reaches the density layers of LSW) the water is saturated to about 80% (Rhein 1991, 1996). The integral effect of deep convection is a vertical homogenization of the CFC profile from the surface to the convection depth. In the example of the Labrador Sea in Fig. 12d, a convection depth of 1900 m is assumed. The additional CFC introduced in the old LSW, however, originates from the upper 600 m (the new LSW volume) that was mixed with the old LSW, and the CFC saturation of the upper 600 m decreases during convection (profile 2, Fig. 12d).
Low CFC saturations of 40%–60% are observed in the convected water column—including the surface water—in the Labrador Sea (Wallace and Lazier 1988; Smethie et al. 2000) as well as other regions (Rhein 1995) because the timescale of convection (days) is too short to allow for air–sea gas exchange to significantly raise the concentration before the convection cell is capped by less dense surface water. In February–March 1997 in the central Labrador Sea, the CFC profiles with surface densities greater than σΘ = 27.74 (11 profiles) had mean CFC-11 concentrations in the upper 2000 m of 4.29 ± 0.11 pmol kg−1, and the profiles with lower surface densities had mean values of 4.07 ± 0.16 pmol kg−1, indicating that the air–sea gas exchange during convection might add about 5% to the CFC-11 content. The mean CFC-11 saturation in the upper 500 m of the profiles with surface densities less than σΘ = 27.74 (i.e., no convection below 500 m) was found to be 79.3% ± 2.1% (4.6 pmol kg−1). Combining these data led us to the assumption that the water above the LSW prior to convection is saturated to 85%. Note that the CFC saturation in the water above the LSW is used for the calculations (i.e., before LSW formation starts), not the saturation of the already convected water column (40%–60%).
The formation rates of new LSW were calculated with this boundary condition. Although the stratification within the LSW is weak and a smoothing of the vertical CFC gradient between sLSW and dLSW away from the formation region is observed, we present sLSW and dLSW formation rates separately. This is somewhat artificial since in the calculations the isopycnal σΘ = 27.77 is treated as a material surface. In doing so and assuming that the dLSW is only formed in years with a high NAO index, the difference between the production rates in years with weak LSW production (low NAO index, only sLSW is formed) and with high LSW production (high NAO index, sLSW and dLSW is formed) can be estimated.
c. Model calculations
For the sake of simplicity, we first ignore the temporal variability in the LSW formation rate. The formation rate is assumed to be constant, and the saturation of newly formed LSW is set to 85%. As already mentioned, the subpolar North Atlantic is regarded as a closed basin where the densities that bound the LSW, the sLSW, and the dLSW (σΘ = 27.74, 27.77, and 27.80) are the material surfaces outside the LSW formation region. The only CFC source of LSW is by deep convection. The atmospheric CFC concentrations are taken from R. Weiss (SIO) [for a discussion see Walker et al. (2000)]. For a constant annual LSW formation rate R the equation is
where CFCinv is the 1997 CFC inventory for LSW, sLSW, or dLSW, respectively; ρ is the density of seawater; C(t) is the CFC concentration in water (Θ = 2.8°C), which is in equilibrium with the atmospheric concentration in year t (t = 1930 … 1997); and 0.85 is the CFC-saturation degree of newly formed LSW.
At a saturation of 85%, the annual mean formation rate of new LSW lies between 4.4 Sv and 5.6 Sv. The formation rate is divided into 2.6–3.4 Sv dLSW and 1.8–2.4 Sv sLSW. These estimates are the lower limit of a mean LSW formation rate compatible with the CFC inventory in 1997. Including the estimated CFC inventory from regions where the CFC and thickness fields could not be interpolated from the observations, the mean LSW formation rate is 4.7–6.0 Sv.
In our model run with a constant LSW formation rate, 32.7% of the CFC inventory was introduced into the LSW layer during the last seven years (1991–97), and in 1990 the CFC-11 inventory was 11.2 million moles. This value is lower by 1.2 million moles than the CFC-11 inventory Smethie and Fine (2001) estimated for the subpolar North Atlantic (12.4 million moles) from the 1986–92 CFC-11 dataset.
Due to the exponential increase of the CFC concentrations from 1930 to mid-1975 and the linear increase thereafter, 56.4% of the inventory entered between 1970 and 1990, and 11% originate from the time period 1930–70. Thus, although our calculation is for 1930–97, our results are mainly representative for the period 1970–97. If a major part of the inventory had come from the years before 1970, the assumption that the subpolar North Atlantic is a closed basin would not be appropriate, since the flushing time of the subpolar North Atlantic is in the order of 20–50 years.
A formation rate of 5 Sv new LSW equals a volume of 16 × 1013 m3 yr−1. Assuming that the water (less dense than LSW before convection) above the LSW extends to depth 500 m before convection starts, the convection area in the Labrador Sea is calculated to 30 × 104 km2 (circle with radius 309 km). Lilly et al. (1999) report a region of about 500 km × 600 km where deep convection could take place in the Labrador Sea, comparable to the aforementioned area. If we assume that the newly formed LSW remains in the western subpolar North Atlantic north of 50°N for one year (area with water depth >1500 m: 2 × 1012 m2), a LSW production of 5 Sv would require a mean annual LSW thickness increase of 70 m, and thickness changes like that have been observed (Curry et al. 1998).
The monitoring of subsurface conditions in the Labrador Sea provides evidence of increased and deepening ventilation, and increasing production of LSW in 1972–76 and in 1988–94 (Dickson et al. 1996; Sy et al. 1997). We repeat the above calculations, but with the assumption that the dense LSW mode, the dLSW (σΘ = 27.77–27.80), was only formed in 1930–60, 1972–76, and 1988–94. As already mentioned, the time period before 1970 is not relevant for the CFC inventory in 1997; the contribution of this time period is marginal no matter what the convection intensity was.
The dLSW production in 1988–94 is calculated to 6.4–8.5 Sv; the sLSW formation rate is kept constant at 1.8–2.4 Sv. The total LSW production rate in the years 1988–94, are then calculated to 8.1–10.8 Sv (sLSW and dLSW), and in 1995–97 the production rate decreases to 1.8–2.4 Sv. One has to keep mind that the calculated formation rates are directly dependent on the assumed CFC saturation in newly formed LSW. A saturation degree of 75% enhances the calculated mean LSW formation rate 0.5–0.7 Sv, and a saturation degree of 95% decreases the formation rate by about the same amount (Table 3).
d. Export of LSW south of 40°N
One of the assumptions in calculating the LSW formation rate was that the subpolar North Atlantic is a closed basin for LSW. This is not entirely true. There is clear evidence of a relatively fast southward spreading of the CFC signal, and new LSW formed after 1988 reached the subtropical Atlantic (26.5°N) in 1996 (Molinari et al. 1998). In 1990, Smethie and Fine (2000) estimated that up to 15% of the CFC inventory in LSW was south of 40°N. An increase of the subtropical fraction is expected since CFC-rich LSW flowed southward with the DWBC in 1991–97.
e. CFC transfer to higher-density classes by entrainment
After overflowing Denmark Strait the transport of the DSOW increases from 2.9 Sv to 5.1 Sv by entrainment (Dickson et al. 1996). Part of the entrained water might be CFC-rich LSW, which by this process is transferred to water with a higher density than σΘ = 27.80. The major part of the entrainment, however, occurs at still depth (600 m) or slightly below (Price and O'Neil-Baringer 1994), and the depth of the upper boundary of LSW, the σΘ = 27.74 layer, is found at depth 900 m south of Denmark Strait at 64°40′N: it cannot be discounted, that part of the entrained water might be LSW.
Entrainment of LSW might also occur south of the sills between Iceland and Scotland. Smethie et al. (2000) estimated from temperature and salinity data that the composition of ISOW after entrainment contains about 35% LSW and 20% northeast Atlantic water. Harvey and Theodorou (1986) reported 25% LSW, and Swift (1984) and Smethie and Swift (1989) found no LSW contribution, but mixing of 59% ISOW and 41% northeast Atlantic water. A multivariate analysis by D.-H. Min resulted in a 16% LSW contribution to ISOW.
To summarize, the amount of LSW and the amount of CFCs transferred to higher-density classes through entrainment is open to debate. Smethie and Fine (2001) estimated that about 1.8 Sv LSW are transferred to higher-density classes by entrainment east of the Reykjanes Ridge and that no LSW is entrained into DSOW.
f. Comparison with other estimates of LSW formation
The calculated and estimated LSW formation rates are summarized in Table 3. In total, assuming constant LSW formation rates, 4.4–5.6 Sv new LSW are required to cover the CFC inventory in the LSW layer in the North Atlantic. If dLSW is only formed in high NAO index periods, the LSW formation rate in 1988–94 is 8.1–10.8 Sv, decreasing to 1.8–2.4 Sv in 1995–97.
Smethie and Fine (2001) did similar calculations for the period before 1990, with a much more limited and nonsynoptic CFC dataset collected between 1986 and 1992. Their boundary condition for CFC input was different from ours in two major ways: 1) they used a CFC-11 saturation of the deep convective layer of 60%, which was measured while deep convection took place; and 2) they excluded the CFC-11 inventory in the region of formation in the central Labrador Sea since their calculation yielded a transport out of the formation region, which was assumed to be at 60% saturation. Assuming a constant formation rate, Smethie and Fine (2000) obtained 7.6 Sv, which is higher than our estimate of 4.4–5.6 Sv. They obtained higher rates because of their high CFC-11 inventory in 1990 (12.4 million moles) compared to the model calculations here (11.2 million moles), by including the CFC-11 inventory south of 40°N (15%), and by the different choice of the boundary conditions. Smethie and Fine (2001) also estimated the formation rate assuming it varied a factor of 10 between high and low formation years and scaled to the thickness of LSW. For the 1970–90 period they obtained minimum, maximum, and average formation rates of 1.3 Sv, 12.5 Sv, and 7.4 Sv. Reducing these values for the CFC-11 inventory south of 40°N yields 1.1 Sv, 10.6 Sv, and 6.3 Sv. Our ratio between low and high formation rates turned out to be 1:5. This estimate is representative for the time period 1970–97.
Our estimates of 8.1–10.8 Sv (1988–94) and 1.8–2.4 Sv (1995–97) are comparable to the LSW overturning rates of 5–10 Sv for 1988–94 and to 1–3 Sv for 1995–97 that Marsh (2000) inferred from surface fluxes while neglecting mixing. One has to take into account that the LSW boundaries were chosen differently: Marsh (2000) defined the LSW bounded by σΘ = 27.65–27.775 compared to 27.74–27.80 (this study). As mentioned before, our calculated LSW formation rates are a lower limit. This might also be true for the Marsh (2000) values: model studies of Nurser et al. (1999) show that when mixing is neglected, the LSW formation rates tend to be lower than the true rate.
8. Summary and conclusions
The CFC-11 and CTD dataset from 1997 in the subpolar North Atlantic is discussed in conjunction with the PALACE float trajectories to study the pathways, the CFC inventory, and formation rates of LSW. Five different scientific groups were involved in the data sampling, and the CFC dataset is found to be consistent at a 4% level. The synopticity of the data was studied by examining the temporal trends for repeated sections during the 1991–98 period. Between 1994 and 1997, the LSW in the western subpolar Atlantic north of 50°N is mainly affected by thickness changes of 100–200 m yr−1, and south of 50°N and in the eastern Atlantic, the regionally different increase of CFCs of 0.05–0.3 pmol kg−1 yr−1 is more important. The LSW thickness decrease north of 50°N is roughly equivalent to a LSW volume loss of 10 Sv.
Since both the CFC contours and the float trajectories mainly follow f/H contours, a topography-following interpolation scheme (Davis 1998) is applied to calculate gridded fields of CFC concentrations and LSW thickness. The interpolation scheme preserves the sharp gradients across isobaths. We used the gridded CFC and thickness fields to calculate the CFC inventory of LSW in 1997 with unprecedented small errors (10%–14%). In 1997, about 28% of the inventory was still in the Labrador Sea west of 45°W, and 31% of the inventory had already intruded the eastern Atlantic (Fig. 13). Following the present temporal trends, the inventory of the eastern Atlantic is expected to increase in coming years, and the CFC inventory in the western Atlantic north of 50°N will decrease, as well as the volume of LSW. Since the CFCs present in the LSW layer were introduced by deep convection, the CFC inventory can be used to estimate lower limits of the LSW formation rates. For the calculations, we combined previous observations in deep convection regions with the Labrador Sea winter data of 1997 and estimated a CFC saturation in the “new” LSW volume of 85%. Furthermore, we assumed that the subpolar North Atlantic is a closed basin for LSW and that the isopycnals that bound the LSW are material surfaces. The LSW formation rate calculated for the time period 1988–97 is a lower limit. The mean LSW formation rate necessary to account for the CFC inventory in 1997 in the subpolar North Atlantic north of 40°N (2300 ± 250 tons CFC) was calculated as 4.4–5.6 Sv. For about 12% of the area with water depth greater than 1000 m the CFC and thickness fields could not be interpolated. Estimating the CFC inventory of these areas by extrapolating from neighboring regions, the CFC inventory for the subpolar North Atlantic increases by 6%–7%. If one takes these regions into account, the mean LSW formation rate is calculated as 4.7–5.9 Sv.
If the denser modes of LSW are only formed in years with high NAO index, the formation rate is calculated as 8.1–10.8 Sv in 1988–94 and 1.8–2.4 Sv in 1995–97. A shift in the saturation degree of 10% would alter the calculated formation rates by about 0.5–0.7 Sv.
Future measurements in the subpolar North Atlantic will benefit from the 1997 fields to evaluate the regionally different temporal evolution of the LSW characteristics and to study the ventilation of LSW. The interpolation method can easily be extended to other water masses in the subpolar North Atlantic, namely to study the spatial distribution of the ISOW and DSOW characteristics. One could also calculate the CO2 inventory of the subpolar North Atlantic by applying the almost linear relationship between CO2 and CFC in this region (Körtzinger et al. 1999). Tracer and inventory fields like the ones presented here can also be used to assess some key aspects in ocean models. One key aspect is the representation of water mass formation since this process links surface water variability (and thus climate variability) to the deep ocean and strongly influences the poleward heat and freshwater transports.
We thank captains, crews, and chief scientists of the cruises Knorr 147, Knorr 151/2, Discovery 230, Valdivia 161, 172, and Meteor 30, legs 2,3, and cruise 39, legs 2–5. We thank K. P. Koltermann, R. Pickart, F. Schott, and A. Sy for providing CTD data. R. Weiss thanks F. A. Van Woy, M. K. Vollmer, and P. K. Salameh for their contributions to the Knorr 151 expedition work and to the CFC data processing. W. M. Smethie would like to thank G. Mathieu for single-handedly measuring the CFC samples on the winter Labrador Sea cruise Knorr 147 and E. Gorman and H. Lee for processing the final data. D. Smythe Wright would like to thank S. M. Boswell and C. R. Harris for their substantial help with the expedition work on Discovery 230. M. Rhein wishes to thank M. Elbrächter and the late K. Bahrenfuss for analyzing the CFC samples on the Valdivia and Meteor cruises. O. Plähn processed the M39-2 dataset and K. Bulsiewicz analyzed and processed the Meteor M39-3 data. Contributions from the Scripps Institution of Oceanography were supported by the Ocean Sciences Section of the U.S. National Science Foundation. W. Smethies's work was supported by ONR Grant 0014-96-0612 and NSF Grant OCE98-11034. The German contributions from IFM Kiel and the University Bremen were funded by the Bundesministerium für Forschung und Technologie (BMBF) and by the Deutsche Forschungsgemeinschaft (DFG).
Current affiliation: Institut für Umweltphysik, Universität Bremen, Bremen, Germany
Current affiliation: School of Oceanography, University of Washington, Seattle, Washington
Corresponding author address: Dr. Monika Rhein, Universität Bremen, FB1, Institut für Umweltphysik, Abt. Ozeanographie, Kufsteiner Strasse, Geb. NW1, 28359 Bremen, Germany. Email: firstname.lastname@example.org